Version 1 (modified by 14 months ago) (diff) | ,
---|
The 'equal perweight mode'
As described elsewhere, elementary reactions introduced on nucleons are connected with a statistical weight, the so called 'perWeight', which is directly connected with the total interaction cross section.
...
How to use it
1) For the rejection method to work, you have to know sigma_max
Please take your favorite jobcard, set
&input ... numTimeSteps =0 ... /
and
&neutrino_induced ... equalWeights_Mode = 1 ... /
Running this jobcard then should generate a file "fort.87" with a list of increasing numbers.
2) Take the last number of that file, add some 'security buffer', and create a second jobcard, now with final state events by setting
&neutrino_induced ... equalWeights_Mode = 2 equalWeights_Max = <The_Value_You_Selected> ... /