INFO: Using directory: "1"
INFO: Using existing c-file: ./SNS_analytic_test.c
INFO: Using existing binary: ./SNS_analytic_test.out
INFO: ===
Simulation 'SNS_analytic_test' (SNS_analytic_test.instr): running on 20 nodes (master is 'hypatia.fysik.dtu.dk', MPI version 3.1).
Using all default parameter values



















Using the input file: /u/data/pkwi/WORK/micromamba/envs/mcstas-3.x/share/mcstas/resources/data/a1Gw2-5-f5_fit_fit.dat
File opened...
## begin version 0.98 fit file header             
## generated at 39
## from binned performance data file a1Gw2-5-f5.dat
## output reported is for 10x12 viewed area.       
## normalized to  2.0807e+14 p/pulse and  1.000 GeV
## (2MW, 1GeV, 60 Hz)
## headerend
##
## section 1 key spectral fit
##
## Fitting function:    
##      I(E) = I*1e12 * exp(-c/sqrt(E))
##                * ( R1*E/(k*T1)**2*exp(-E/(k*T1)) + R2*E/(k*T2)**2*exp(-E/(k*T2)
##                    + R3*E/(k*T3)**2*exp(-(E/(k*T3))**b) + D(E)*rho(E)/E**(1-a) )
##      with 
##      D(E) = 1/(1+(Ecut/E)**s)    
##      rho(E) = 1 + delta*exp(-x)(1 + x +0.5*x**2)
##      x(E) = g*(E-2B); for E>2B
##           = 0;        for E<=2B
##
##      constants:
##      k = 1.3805e-23 J/K = 8.617e-5 eV/K
##      B = 7.36e-3 eV
##     
##      parameters:
##        c  ,   R1   ,   T1,   R2,   T2,   R3,   T3,   a,   b,   Ecut ,   s,   delta,   g,   I
## data start
##
## end section 1
## 
## section 2: pulse shape fit
## using the fit function
##      f(E,t) = (1-R)*A/2*(A*(t-to*10))**2 *exp(-A*(t-to*10)) 
##                    +R*B*(A/(A-B))**3 *[ exp(-B*(t-to*10)) 
##                      - exp(-A*(t-to*10))*(1+(A-B)*(t-to*10)+0.5*(A-B)**2*(t-to*10)**2) ]
##   with 
##      A, B, R, to
##
## headerend
##
## subsection 1 fit of parameter A=A(E)
## Fitting function: Pade' type   
##      f(x) = a*x**b*(1+c*x+d*x**2+(x/f)**g)/(1+h*x+i*x**2+(x/j)**k)
##   with 
##      a,b,c,d,f,g,h,i,j,k
## data start
## 
## end subsection 1
## 
## subsection 2 fit of parameter B=B(E)
## Fitting function: Pade' type   
##      f(x) = a*x**b*(1+c*x+d*x**2+(x/f)**g)/(1+h*x+i*x**2+(x/j)**k)
##   with 
##      a,b,c,d,f,g,h,i,j,k
## data start
## 
## end subsection 2
## 
## 
## subsection 3 fit of parameter R=R(E)
## fitting function Padee type   
##      f(x) = a*x**b*(1+c*x+d*c**2+(x/f)**g)/(1+h*x+i*x**2+(x/j)**k)
##   with 
##      a,b,c,d,f,g,h,i,j,k
## data start
## 
## end subsection 3
## 
## 
## subsection 4 fit of parameter to-to(E)
## fitting function Padee type   
##      f(x) = a*x**b*(1+c*x+d*c**2+(x/f)**g)/(1+h*x+i*x**2+(x/j)**k)
##   with 
##      a,b,c,d,f,g,h,i,j,k
## data start

   para_sp(0)=1.225883e-02
   para_sp(1)=5.013027e-01
   para_sp(2)=7.393808e+00
   para_sp(3)=4.918982e+01
   para_sp(4)=7.906112e+01
   para_sp(5)=-4.362656e+01
   para_sp(6)=1.149904e+02
   para_sp(7)=1.806666e-02
   para_sp(8)=1.173401e+00
   para_sp(9)=2.500000e-02
   para_sp(10)=8.095194e-01
   para_sp(11)=6.141977e-01
   para_sp(12)=2.226451e+01
   para_sp(13)=4.267161e-01

   para_a(0)=2.490382e+00
   para_a(1)=4.816776e-01
   para_a(2)=-9.941653e+03
   para_a(3)=6.624781e+04
   para_a(4)=7.977572e-05
   para_a(5)=1.079246e+00
   para_a(6)=2.619925e+04
   para_a(7)=1.034427e+05
   para_a(8)=3.935614e-03
   para_a(9)=1.974287e+00

   para_b(0)=2.235995e-06
   para_b(1)=-6.894272e-01
   para_b(2)=1.159062e+04
   para_b(3)=7.411655e+04
   para_b(4)=2.609572e-02
   para_b(5)=5.868703e+00
   para_b(6)=1.447994e+01
   para_b(7)=2.976349e+01
   para_b(8)=1.145569e-01
   para_b(9)=4.699988e+00

   para_R(0)=4.009638e-01
   para_R(1)=0.000000e+00
   para_R(2)=-3.162201e+03
   para_R(3)=1.050311e+03
   para_R(4)=2.964059e-04
   para_R(5)=9.949269e-01
   para_R(6)=-4.571167e+00
   para_R(7)=3.076276e+03

   para_to(0)=1.300000e+00
   para_to(1)=-5.000002e-02
   para_to(2)=-1.950000e+02
   para_to(3)=5.640000e+02
   para_to(4)=3.588075e-07
   para_to(5)=1.000000e+00
   para_to(6)=-1.000000e+02
   para_to(7)=7.200000e+01
   para_to(8)=5.283668e-06
   para_to(9)=1.400000e+00

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12

   CItot =   3.5504e+12
*** TRACE end *** 
Detector: dett_I=2.13987e+10 dett_ERR=1.35025e+07 dett_N=2.5116e+06 "tout.det"
Detector: det_I=8.51301e+10 det_ERR=2.69315e+07 det_N=9.99183e+06 "det_1775709523.x_y"
INFO: Placing instr file copy SNS_analytic_test.instr in dataset 1
INFO: Placing generated c-code copy SNS_analytic_test.c in dataset 1