Difference between revisions of "Radphi MonteCarlo"
Jump to navigation
Jump to search
Line 24: | Line 24: | ||
</tr> | </tr> | ||
</table> | </table> | ||
+ | |||
===Results of the simulation=== | ===Results of the simulation=== | ||
* γ p -> p (ω -> π<sup>0</sup>γ) | * γ p -> p (ω -> π<sup>0</sup>γ) | ||
* γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ) | * γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ) |
Revision as of 10:59, 11 June 2014
To generate Radphi Monte Carlo events steps below are followed:
- Edit the decayt.dat[1] file and specify the reaction of interest and number of events to generate.
- Run mcwrap[2], mcwrap requires cebaf.geom [3] to run, and specify output name.
- ./mcwrap -l4.390 -u5.390 -comega
- Output file generated with output_name.itape (omega.itape).
- Move this file to the folder Gradphi.
- Edit events.in[4] and control.in[5] files and the run gradphi.x[6].
Resources needed to generate 100000000 events.
- The calculation below is based on 300,000 omega events that were generated using mcwrap and simulated using gradphi on gluey.
- To simulate 300,000 events it took 6840 seconds.
Number generated(millions) | CPU time per event (s) | Output file size (MB) | Total time in days |
100 | 0.0216 | 59000 | 25 |
Results of the simulation
- γ p -> p (ω -> π0γ)
- γ p -> (Δ+ -> p π0) (ω -> π0γ)