Difference between revisions of "Radphi MonteCarlo"

From UConn PAN
Jump to navigation Jump to search
Line 10: Line 10:
 
* The calculation below is based on 300,000 omega events that were  generated using mcwrap and simulated using gradphi on gluey.
 
* 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.
 
*To simulate 300,000 events it took 6840 seconds.
*:<table border ="5">
+
<ul>
 +
<table border ="5">
 
   <tr>
 
   <tr>
 
     <td>Number generated(millions)</td>  
 
     <td>Number generated(millions)</td>  
Line 23: Line 24:
 
   </tr>
 
   </tr>
 
  </table>
 
  </table>
 
+
</ul>
  
 
===Results of the simulation===
 
===Results of the simulation===
 
* &gamma; p -> p (&omega; -> &pi;<sup>0</sup>&gamma;)
 
* &gamma; p -> p (&omega; -> &pi;<sup>0</sup>&gamma;)
 
* &gamma; p -> (&Delta;<sup>+</sup> -> p &pi;<sup>0</sup>) (&omega; -> &pi;<sup>0</sup>&gamma;)
 
* &gamma; p -> (&Delta;<sup>+</sup> -> p &pi;<sup>0</sup>) (&omega; -> &pi;<sup>0</sup>&gamma;)

Revision as of 11:00, 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γ)