| 1 | = The 'equal perweight mode' = |
| 2 | |
| 3 | 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. |
| 4 | |
| 5 | ... |
| 6 | |
| 7 | == How to use it == |
| 8 | |
| 9 | 1) For the rejection method to work, you have to know sigma_max |
| 10 | |
| 11 | Please take your favorite jobcard, set |
| 12 | {{{ |
| 13 | &input |
| 14 | ... |
| 15 | numTimeSteps =0 |
| 16 | ... |
| 17 | / |
| 18 | }}} |
| 19 | and |
| 20 | {{{ |
| 21 | &neutrino_induced |
| 22 | ... |
| 23 | equalWeights_Mode = 1 |
| 24 | ... |
| 25 | / |
| 26 | }}} |
| 27 | Running this jobcard then should generate a file "fort.87" with a list of increasing numbers. |
| 28 | |
| 29 | 2) Take the last number of that file, add some 'security buffer', and create a second jobcard, now with final state events by setting |
| 30 | {{{ |
| 31 | &neutrino_induced |
| 32 | ... |
| 33 | equalWeights_Mode = 2 |
| 34 | equalWeights_Max = <The_Value_You_Selected> |
| 35 | ... |
| 36 | / |
| 37 | }}} |
| 38 | |