Difference between revisions of "Radphi MonteCarlo"

From UConn PAN
Jump to navigation Jump to search
Line 29: Line 29:
 
===Results of the simulation===
 
===Results of the simulation===
 
* 10<sup>8</sup> events of &gamma; p -> p (&omega; -> &pi;<sup>0</sup>&gamma;)
 
* 10<sup>8</sup> events of &gamma; p -> p (&omega; -> &pi;<sup>0</sup>&gamma;)
* 10<sup>8</sup> events of &gamma; p -> (&Delta;<sup>+</sup> -> p &pi;<sup>0</sup>) (&omega; -> &pi;<sup>0</sup>&gamma; &pi;<sup>+</sup> &pi;<sup>-</sup>)
+
* 10<sup>8</sup> events of &gamma; p -> (&Delta;<sup>+</sup> -> p &pi;<sup>0</sup>) (&omega; -> &pi;<sup>0</sup> &pi;<sup>+</sup> &pi;<sup>-</sup>)
 
* 10<sup>8</sup> events of &gamma; p -> (&Delta;<sup>+</sup> -> p &pi;<sup>0</sup>) (&omega; -> &pi;<sup>0</sup>&gamma;)
 
* 10<sup>8</sup> events of &gamma; p -> (&Delta;<sup>+</sup> -> p &pi;<sup>0</sup>) (&omega; -> &pi;<sup>0</sup>&gamma;)

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γ)