Radphi MonteCarlo

From UConn PAN
Jump to navigation Jump to search

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

All results from this simulation are stored on the UConn srm. They are available on the UConn cluster in the nfs mounted directory /pnfs/phys.uconn.edu/data/Gluex/radphi/omega_sim-6-2014 or remotely at the srm://grinch.phys.uconn.edu/Gluex/radphi/omega_sim-6-2014 url.

  • 108 events of γ p -> p (ω -> π0γ)
  • 108 events of γ p -> (Δ+ -> p π0) (ω -> π0 π+ π-)
  • 108 events of γ p -> (Δ+ -> p π0) (ω -> π0γ)