Radphi MonteCarlo

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

• 10⁸ events of γ p -> p (ω -> π⁰γ)
• 10⁸ events of γ p -> (Δ⁺ -> p π⁰) (ω -> π⁰ π⁺ π^-)
• 10⁸ events of γ p -> (Δ⁺ -> p π⁰) (ω -> π⁰γ)