INFO: Using directory: "1" INFO: Using existing c-file: ./ILL_IN6.c INFO: Using existing binary: ./ILL_IN6.out INFO: === Instrument Simulation ILL_IN6 (/home/nexmap/pkwi/TESTS/2022-03-01/McStas_CPU_GCC_KISS_5e6/ILL_IN6/ILL_IN6.instr) using computed monochromator take-off angles: 75.4011 76.1933 76.9908 [deg] Wavelength [AA] 4.14 Neutron velocity [m/s] 955.564 Monochr. Bragg angle [deg] 38.0966 Incident Energy [meV] 4.77282 Focusing Energy transfer[meV] 0 Travel time: Supp./Fermi [us] 209.301 Travel time: Supp./Det. [us] 3294.39 TOF Delay [us] 282.75 TOF Dead Time [us] 2.37202 TOF Period (1 cycle) [us] 6034.37 TOF Channel width [us] 47.125 CHOP Fermi Phase [deg] -6.24325 CHOP Suppressor Phase [deg] 0 CHOP Fermi Speed [rpm] 4971.52 CHOP Suppressor Speed [rpm] 4971.52 Number of time channels 128 Current Elastic Peak Ch. 64 Elast. peak ch. for 0-delay 70 Time window: min=0.56114 max=6.59314 delay=0.28275 tof-width=6.032 [ms] [ILL_IN6] Initialize Source_gen: component VCS (square 0.03 x 0.2) focusing distance dist=2.55 area=0.03 x 0.2 spectra 4.065 to 4.215 AA (4.605 to 4.951 meV) T1=216.8 K (2.093 AA), T2=33.9 K (5.292 AA), T3=16.7 K (7.539 AA) Flux is dPhi/dlambda in [n/s/cm2]. Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/HOPG.rfl' (Table_Read_Offset) Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/HOPG.trm' (Table_Read_Offset) Reusing input file 'HOPG.rfl' (Table_Read_Offset) Reusing input file 'HOPG.trm' (Table_Read_Offset) Reusing input file 'HOPG.rfl' (Table_Read_Offset) Reusing input file 'HOPG.trm' (Table_Read_Offset) FermiChopper: Suppressor: Frequency nu=82.8587 [Hz] 4971.52 [rpm], time frame=0.0241375 [s] phase=0 [deg] FermiChopper: FermiM: Frequency nu=82.8587 [Hz] 4971.52 [rpm], time frame=0.0241375 [s] phase=-6.24325 [deg] Isotropic_Sqw: Sample: is a hollow cylinder: radius=0.010000 thickness=0.005000 height=0.055000 [J Comp Phys 228 (2009) 5251] Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_inc.sqw' (Table_Read_Offset) Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_inc.sqw' (Table_Read_Offset) Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_inc.sqw' (Table_Read_Offset) Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_inc.sqw' (Table_Read_Offset) Isotropic_Sqw: Sample: Temperature set to T=350 [K] Isotropic_Sqw: Sample: S(q,w) data from Rb_liq_inc.sqw (inc) assumed to be classical (symmetrised in energy). Isotropic_Sqw: Sample: Applying exp(hw/2kT) factor with T=350 [K] on Rb_liq_inc.sqw file (classical/symmetric) using 'Frommhold' quantum correction Isotropic_Sqw: Sample: Generated Rb_liq_inc.sqw incoherent Sqw q=[0.1:10 Angs-1] w=[0:101.849 meV] |S|=0.985192 size=[101x10185] sigma=0.5 [barns] normalization factor S(q,w)*0.492586 (auto) \int q^2 S(q) dq=676.701 Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_coh.sqw' (Table_Read_Offset) Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_coh.sqw' (Table_Read_Offset) Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_coh.sqw' (Table_Read_Offset) Opening input file '/u/data/pkwi/McStas/mcstas/3.x-dev/tools/Python/mcrun/../mccodelib/../../../data/Rb_liq_coh.sqw' (Table_Read_Offset) Isotropic_Sqw: Sample: Temperature set to T=350 [K] Isotropic_Sqw: Sample: S(q,w) data from Rb_liq_coh.sqw (coh) assumed to be classical (symmetrised in energy). Isotropic_Sqw: Sample: Applying exp(hw/2kT) factor with T=350 [K] on Rb_liq_coh.sqw file (classical/symmetric) using 'Frommhold' quantum correction Isotropic_Sqw: Sample: Generated Rb_liq_coh.sqw coherent Sqw q=[0.1:10 Angs-1] w=[0:101.849 meV] |S|=0.887933 size=[101x10185] sigma=6.32 [barns] normalization factor S(q,w)*0.477012 (auto) \int q^2 S(q) dq=698.369 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.633013 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.961699 in slit does not match its length=0.012 (slit exit expected). Trace ETA 9.95 [min] % 1 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.971742 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.694046 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.648148 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.682192 in slit does not match its length=0.012 (slit exit expected). 10 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.657376 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.965614 in slit does not match its length=0.012 (slit exit expected). 20 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.628732 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.961564 in slit does not match its length=0.012 (slit exit expected). 30 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.639501 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.957871 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.662009 in slit does not match its length=0.012 (slit exit expected). 40 50 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.947826 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.979484 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.67558 in slit does not match its length=0.012 (slit exit expected). 60 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.964208 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.958054 in slit does not match its length=0.012 (slit exit expected). 70 80 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.650225 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.972371 in slit does not match its length=0.012 (slit exit expected). FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.632807 in slit does not match its length=0.012 (slit exit expected). 90 FermiChopper: Suppressor: ABSORB Neutron propagation time v*dt/length=0.946394 in slit does not match its length=0.012 (slit exit expected). *** TRACE end *** Save [ILL_IN6] Detector: Mon_SuppInL_I=7.29436e+09 Mon_SuppInL_ERR=1.99899e+06 Mon_SuppInL_N=1.33579e+07 "Mon_SuppInL_1646092546.L" Detector: Mon_SuppInT_I=7.68347e+09 Mon_SuppInT_ERR=2.05531e+06 Mon_SuppInT_N=1.4021e+07 "Mon_SuppInT_1646092546.t" FermiChopper: Suppressor Mean velocity v = 954.888 [m/s] Mean wavelength lambda= 4.14293 [Angs] Mean energy omega = 4.76607 [meV] Mean arrival time t = 0.00228245 [s] Mean phase = 70.2703 [deg] (use phase=-(this value) to optimize) Slit pack divergence = 154.011 [deg] (full width) Opening time dt = 0.295824 [s] Intensity reaching FC = 3.63336e+09 [n/s] Intensity passing FC = 9.71129e+08 [n/s] FermiChopper: Suppressor: Lost events anaylsis Already inside: 0 By Top/Bottom of cylinder: 764170 At cylinder entrance: 4767031 Hit cyl. entrance sides: 0 Can't prop. to slit pack: 144808 At slit pack entrance: 1201597 On absorbing slit coating: 941316 Exiting slit pack: 0 Too many iterations: 0 Prop. in wrong direction : 0 Mad neutron (no control): 41 At cylinder exit: 2315677 Warning: Can not reach slit wall: 0 Warning: Can not reach slit output: 86 Detector: Mon_SuppOutT_I=9.63706e+08 Mon_SuppOutT_ERR=346700 Mon_SuppOutT_N=7.75271e+06 "Mon_SuppOutT_1646092546.t" FermiChopper: FermiM Mean velocity v = 957.074 [m/s] Mean wavelength lambda= 4.13347 [Angs] Mean energy omega = 4.78792 [meV] Mean arrival time t = 0.000161106 [s] Mean phase = 6.05234 [deg] (use phase=-(this value) to optimize) Slit pack divergence = 2.10061 [deg] (full width) Opening time dt = 0.00403485 [s] Intensity reaching FC = 6.36867e+08 [n/s] Intensity passing FC = 2.57397e+07 [n/s] FermiChopper: FermiM: Lost events anaylsis Already inside: 0 By Top/Bottom of cylinder: 129361 At cylinder entrance: 117742 Hit cyl. entrance sides: 0 Can't prop. to slit pack: 7 At slit pack entrance: 5530 On absorbing slit coating: 4791005 Exiting slit pack: 0 Too many iterations: 0 Prop. in wrong direction : 0 Mad neutron (no control): 0 At cylinder exit: 388 Warning: Can not reach slit wall: 0 Warning: Can not reach slit output: 3 Detector: Mon_FermiOutdT_I=2.52893e+07 Mon_FermiOutdT_ERR=56142.9 Mon_FermiOutdT_N=203589 "Mon_FermiOutdT_1646092546.t" Detector: MonokMonitor_I=66731.2 MonokMonitor_ERR=144.179 MonokMonitor_N=214946 "MonokMonitor_1646092546.U1_L" Detector: Mon_SampleInT_I=830707 Mon_SampleInT_ERR=2034.8 Mon_SampleInT_N=167236 "Mon_SampleInT_1646092546.t" Detector: Mon_SampleInXY_I=444875 Mon_SampleInXY_ERR=961.194 Mon_SampleInXY_N=214946 "Mon_SampleInXY_1646092546.x_y" Detector: Out_I=705334 Out_ERR=939.161 Out_N=684904 "out.dat" Detector: M_theta_t_all_I=121390 M_theta_t_all_ERR=391.094 M_theta_t_all_N=117869 "M_theta_t_all_1646092563.A_E" Detector: Detector_nM_I=228511 Detector_nM_ERR=536.182 Detector_nM_N=222090 "Detector_nM.dat" Finally [ILL_IN6: 1]. Time: 9.76667 [min] Isotropic_Sqw: Sample: 2605477 neutron events (out of 3427148) that should have scattered were transmitted because scattering conditions WARNING could not be satisfied after 100 tries. Isotropic_Sqw: Sample: 2445 neutron events (out of 821671) reached WARNING unrealistic weight. The S(q,w) norm might be too high. Isotropic_Sqw: Sample: Scattering fraction=0.0751957 of incoming intensity Absorption fraction =0.12305 Single scattering intensity =669507 (coh=615360 inc=54147.1) Multiple scattering intensity =35827 INFO: Placing instr file copy /home/nexmap/pkwi/TESTS/2022-03-01/McStas_CPU_GCC_KISS_5e6/ILL_IN6/ILL_IN6.instr in dataset 1 loading system configuration loading override configuration from /home/nexmap/pkwi/McCode/tools/Python/mctest/../mccodelib/mcstas-test/McStas_CPU_GCC_KISS/mccode_config.json