Difference between revisions of "Radphi MonteCarlo"
Jump to navigation
Jump to search
| Line 28: | Line 28: | ||
===Results of the simulation=== | ===Results of the simulation=== | ||
| − | * γ p -> p (ω -> π<sup>0</sup>γ) | + | * 10<sup>8</sup> events of γ p -> p (ω -> π<sup>0</sup>γ) |
| − | * γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ) | + | * 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ π<sup>+</sup> π<sup>-</sup>) |
| + | * 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ) | ||
Revision as of 11:24, 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
- 108 events of γ p -> p (ω -> π0γ)
- 108 events of γ p -> (Δ+ -> p π0) (ω -> π0γ π+ π-)
- 108 events of γ p -> (Δ+ -> p π0) (ω -> π0γ)