Inter-component transformations

Thu Feb 4 11:33:36 CET 1999

Hi Kristian,

Thank you for your help.  

I tried to change the the compiler optimation,
but there was no change in the results. I will run the simulation on
another computer this weekend.

Maybe your simulation was to short to see the effect. There is no change
of x and vx for the first 
neutrons, but whenever a
neutron can pass the whole instrument without absorbtion, the factor,
which is changing x and vx between source and fermi, 
increases
 for the following neutrons.                                           

I have added the output of the first twenty neutrons with -seed1. The
fifth neutron could pass the fermi-component and after that, x
and vx values are changed by 1.07...and so on. 

If you change the variable dt in the component
TOF_monitor from 0.1 to 0.7 or bigger the simulation works fine. 

 - Philipp

-------------- next part --------------
INSTRUMENT:
COMPONENT: "a1"
POS: 0, 0, 0, 1, 0, 0, -0, 1, 0, 0, -0, 1
COMPONENT: "source"
POS: 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1
COMPONENT: "fermi"
POS: 0, 0, 0.1, 1, 0, 0, 0, 1, 0, 0, 0, 1
COMPONENT: "mon"
POS: 0, 0, 0.2, 1, 0, 0, 0, 1, 0, 0, 0, 1
INSTRUMENT END:
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.0340188, -0.0105617, 0, 4.47967, 4.72242, 1977.95, 0.000364659, 0, 0, 1
COMP: "fermi"
STATE: 0.0340188, -0.0105617, -0.1, 4.47967, 4.72242, 1977.95, 0.000364659, 0, 0, 1
STATE: 0.0341626, -0.0104101, -0.0365092, 4.47967, 4.72242, 1977.95, 0.000396758, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.0341626, -0.0104101, -0.0365092, 4.47967, 4.72242, 1977.95, 0.000396758, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.0302449, -0.0164777, 0, 4.24438, -3.51642, 1977.96, 0.000221588, 0, 0, 1
COMP: "fermi"
STATE: -0.0302449, -0.0164777, -0.1, 4.24438, -3.51642, 1977.96, 0.000221588, 0, 0, 1
STATE: -0.0301159, -0.0165845, -0.0399128, 4.24438, -3.51642, 1977.96, 0.000251966, 0, 0, 1
ABSORB:
LEAVE:
STATE: -0.0301159, -0.0165845, -0.0399128, 4.24438, -3.51642, 1977.96, 0.000251966, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.00226029, 0.0128871, 0, -2.13961, 0.212052, 1977.96, 0.000380892, 0, 0, 1
COMP: "fermi"
STATE: -0.00226029, 0.0128871, -0.1, -2.13961, 0.212052, 1977.96, 0.000380892, 0, 0, 1
STATE: -0.00231444, 0.0128925, -0.0499464, -2.13961, 0.212052, 1977.96, 0.000406198, 0, 0, 1
ABSORB:
LEAVE:
STATE: -0.00231444, 0.0128925, -0.0499464, -2.13961, 0.212052, 1977.96, 0.000406198, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.0416195, 0.0135712, 0, 3.43843, -5.67117, 1977.95, 0.000242788, 0, 0, 1
COMP: "fermi"
STATE: 0.0416195, 0.0135712, -0.1, 3.43843, -5.67117, 1977.95, 0.000242788, 0, 0, 1
STATE: 0.0417455, 0.0133634, -0.0275193, 3.43843, -5.67117, 1977.95, 0.000279432, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.0417455, 0.0133634, -0.0275193, 3.43843, -5.67117, 1977.95, 0.000279432, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.0483699, -0.0257113, 0, -5.74032, 4.81318, 1977.95, 6.26716e-05, 0, 0, 1
COMP: "fermi"
STATE: -0.0483699, -0.0257113, -0.1, -5.74032, 4.81318, 1977.95, 6.26716e-05, 0, 0, 1
STATE: -0.0486264, -0.0254963, -0.0116394, -5.74032, 4.81318, 1977.95, 0.000107344, 0, 0, 0.0775835
COMP: "mon"
STATE: -0.0486264, -0.0254963, -0.111639, -5.74032, 4.81318, 1977.95, 0.000107344, 0, 0, 0.0775835
STATE: -0.0489504, -0.0252246, 0, -5.74032, 4.81318, 1977.95, 0.000163786, 0, 0, 0.0775835
LEAVE:
STATE: -0.0489504, -0.0252246, 0, -5.74032, 4.81318, 1977.95, 0.000163786, 0, 0, 0.0775835

---------> neutron passes the instrument -----> now x and vx are changed by 1.07...

ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.00990556, -0.037021, 0, -6.19004, 7.89478, 1977.94, 8.73028e-05, 0, 0, 1
COMP: "fermi"
STATE: -0.0106741, -0.037021, -0.1, -6.67029, 7.89478, 1977.94, 8.73028e-05, 0, 0, 1
STATE: -0.0108467, -0.0368166, -0.0488093, -6.67029, 7.89478, 1977.94, 0.000113184, 0, 0, 1
ABSORB:
LEAVE:
STATE: -0.0108467, -0.0368166, -0.0488093, -6.67029, 7.89478, 1977.94, 0.000113184, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.00129324, 0.0339112, 0, 1.78238, -3.22754, 1977.96, 0.000255021, 0, 0, 1
COMP: "fermi"
STATE: 0.00139357, 0.0339112, -0.1, 1.92066, -3.22754, 1977.96, 0.000255021, 0, 0, 1
STATE: 0.00144215, 0.0338296, -0.0499792, 1.92066, -3.22754, 1977.96, 0.00028031, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.00144215, 0.0338296, -0.0499792, 1.92066, -3.22754, 1977.96, 0.00028031, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.00242872, -0.000641701, 0, 7.48103, -3.28313, 1977.95, 0.000308543, 0, 0, 1
COMP: "fermi"
STATE: 0.00261715, -0.000641701, -0.1, 8.06144, -3.28313, 1977.95, 0.000308543, 0, 0, 1
STATE: 0.00282125, -0.000724827, -0.0499203, 8.06144, -3.28313, 1977.95, 0.000333862, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.00282125, -0.000724827, -0.0499203, 8.06144, -3.28313, 1977.95, 0.000333862, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.0026745, 0.0269914, 0, -1.57875, 6.19544, 1977.96, 0.000113326, 0, 0, 1
COMP: "fermi"
STATE: 0.00288199, 0.0269914, -0.1, -1.70123, 6.19544, 1977.96, 0.000113326, 0, 0, 1
STATE: 0.00283892, 0.0271482, -0.0499193, -1.70123, 6.19544, 1977.96, 0.000138645, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.00283892, 0.0271482, -0.0499193, -1.70123, 6.19544, 1977.96, 0.000138645, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.0147542, 0.0307725, 0, 6.63051, -6.80802, 1977.94, 0.000379731, 0, 0, 1
COMP: "fermi"
STATE: -0.0158988, 0.0307725, -0.1, 7.14493, -6.80802, 1977.94, 0.000379731, 0, 0, 1
STATE: -0.0157091, 0.0305916, -0.0474681, 7.14493, -6.80802, 1977.94, 0.00040629, 0, 0, 1
ABSORB:
LEAVE:
STATE: -0.0157091, 0.0305916, -0.0474681, 7.14493, -6.80802, 1977.94, 0.00040629, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.00259953, -0.0413944, 0, -4.87031, 2.58285, 1977.96, 0.000356093, 0, 0, 1
COMP: "fermi"
STATE: 0.00280122, -0.0413944, -0.1, -5.24817, 2.58285, 1977.96, 0.000356093, 0, 0, 1
STATE: 0.00266836, -0.041329, -0.0499287, -5.24817, 2.58285, 1977.96, 0.000381408, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.00266836, -0.041329, -0.0499287, -5.24817, 2.58285, 1977.96, 0.000381408, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.0151107, -0.0435829, 0, -7.595, -0.669311, 1977.95, 2.52383e-05, 0, 0, 1
COMP: "fermi"
STATE: -0.016283, -0.0435829, -0.1, -8.18425, -0.669311, 1977.95, 2.52383e-05, 0, 0, 1
STATE: -0.0165015, -0.0436007, -0.0471985, -8.18425, -0.669311, 1977.95, 5.19334e-05, 0, 0, 1
ABSORB:
LEAVE:
STATE: -0.0165015, -0.0436007, -0.0471985, -8.18425, -0.669311, 1977.95, 5.19334e-05, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.026172, 0.0470634, 0, 6.36439, 5.55282, 1977.95, 0.000106666, 0, 0, 1
COMP: "fermi"
STATE: -0.0282025, 0.0470634, -0.1, 6.85817, 5.55282, 1977.95, 0.000106666, 0, 0, 1
STATE: -0.0279994, 0.0472279, -0.041425, 6.85817, 5.55282, 1977.95, 0.00013628, 0, 0, 1
ABSORB:
LEAVE:
STATE: -0.0279994, 0.0472279, -0.041425, 6.85817, 5.55282, 1977.95, 0.00013628, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.00397603, -0.0124793, 0, 4.1181, 0.198356, 1977.96, 0.00026709, 0, 0, 1
COMP: "fermi"
STATE: 0.00428451, -0.0124793, -0.1, 4.4376, 0.198356, 1977.96, 0.00026709, 0, 0, 1
STATE: 0.00439712, -0.0124743, -0.0498063, 4.4376, 0.198356, 1977.96, 0.000292466, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.00439712, -0.0124743, -0.0498063, 4.4376, 0.198356, 1977.96, 0.000292466, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.00316064, -0.046072, 0, -0.9868, 6.83322, 1977.95, 0.000372324, 0, 0, 1
COMP: "fermi"
STATE: 0.00340586, -0.046072, -0.1, -1.06336, 6.83322, 1977.95, 0.000372324, 0, 0, 1
STATE: 0.00337892, -0.0458988, -0.0498857, -1.06336, 6.83322, 1977.95, 0.00039766, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.00337892, -0.0458988, -0.0498857, -1.06336, 6.83322, 1977.95, 0.00039766, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.0220952, -0.0215707, 0, 3.7745, 2.21498, 1977.96, 0.000141619, 0, 0, 1
COMP: "fermi"
STATE: 0.0238095, -0.0215707, -0.1, 4.06734, 2.21498, 1977.96, 0.000141619, 0, 0, 1
STATE: 0.0239248, -0.0215078, -0.0439045, 4.06734, 2.21498, 1977.96, 0.00016998, 0, 0, 0.290076
COMP: "mon"
STATE: 0.0239248, -0.0215078, -0.143904, 4.06734, 2.21498, 1977.96, 0.00016998, 0, 0, 0.290076
STATE: 0.0242207, -0.0213467, 0, 4.06734, 2.21498, 1977.96, 0.000242734, 0, 0, 0.290076
LEAVE:
STATE: 0.0242207, -0.0213467, 0, 4.06734, 2.21498, 1977.96, 0.000242734, 0, 0, 0.290076

------------> next neutron passed -------------> x and vy are now changed by 1.36.....

ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.0187861, -0.0334026, 0, -0.947765, 6.0142, 1977.96, 0.00033168, 0, 0, 1
COMP: "fermi"
STATE: 0.025693, -0.0334026, -0.1, -1.29622, 6.0142, 1977.96, 0.00033168, 0, 0, 1
STATE: 0.0256556, -0.033229, -0.0429161, -1.29622, 6.0142, 1977.96, 0.00036054, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.0256556, -0.033229, -0.0429161, -1.29622, 6.0142, 1977.96, 0.00036054, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: -0.0169663, -0.0271032, 0, 6.2246, -2.36785, 1977.95, 0.000274668, 0, 0, 1
COMP: "fermi"
STATE: -0.0232041, -0.0271032, -0.1, 8.51314, -2.36785, 1977.95, 0.000274668, 0, 0, 1
STATE: -0.0229649, -0.0271697, -0.0444141, 8.51314, -2.36785, 1977.95, 0.000302771, 0, 0, 1
ABSORB:
LEAVE:
STATE: -0.0229649, -0.0271697, -0.0444141, 8.51314, -2.36785, 1977.95, 0.000302771, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.0456468, 0.00886401, 0, 2.48912, 5.67558, 1977.96, 0.000175824, 0, 0, 1
COMP: "fermi"
STATE: 0.0624293, 0.00886401, -0.1, 3.40428, 5.67558, 1977.96, 0.000175824, 0, 0, 1
STATE: 0.0624293, 0.00886401, -0.1, 3.40428, 5.67558, 1977.96, 0.000175824, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.0624293, 0.00886401, -0.1, 3.40428, 5.67558, 1977.96, 0.000175824, 0, 0, 1
ENTER:
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "a1"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
COMP: "source"
STATE: 0, 0, 0, 0, 0, 1, 0, 0, 0, 1
STATE: 0.042397, -0.0101563, 0, 4.98077, 2.91501, 1977.96, 0.000364389, 0, 0, 1
COMP: "fermi"
STATE: 0.0579846, -0.0101563, -0.1, 6.812, 2.91501, 1977.96, 0.000364389, 0, 0, 1
STATE: 0.0579846, -0.0101563, -0.1, 6.812, 2.91501, 1977.96, 0.000364389, 0, 0, 1
ABSORB:
LEAVE:
STATE: 0.0579846, -0.0101563, -0.1, 6.812, 2.91501, 1977.96, 0.000364389, 0, 0, 1
Neutrons: 2, Deviation: 81.6172, Deviation: 81.617%
Detector: mon_I=0 mon_ERR=0