Difference between revisions of "Radphi MonteCarlo"
Jump to navigation
Jump to search
| Line 1: | Line 1: | ||
| − | + | ===To generate Radphi Monte Carlo events steps below are followed:=== | |
*Edit the decayt.dat[http://zeus.phys.uconn.edu/~mokaya/decayt.dat] file and specify the reaction of interest and number of events to generate. | *Edit the decayt.dat[http://zeus.phys.uconn.edu/~mokaya/decayt.dat] file and specify the reaction of interest and number of events to generate. | ||
*Run mcwrap[http://zeus.phys.uconn.edu/~mokaya/mcwrap], mcwrap requires cebaf.geom [http://zeus.phys.uconn.edu/~mokaya/cebaf.geom] to run, and specify output name. | *Run mcwrap[http://zeus.phys.uconn.edu/~mokaya/mcwrap], mcwrap requires cebaf.geom [http://zeus.phys.uconn.edu/~mokaya/cebaf.geom] to run, and specify output name. | ||
| Line 7: | Line 7: | ||
*Edit events.in[http://zeus.phys.uconn.edu/~mokaya/events.in] and control.in[http://zeus.phys.uconn.edu/~mokaya/control.in] files and the run gradphi.x[http://zeus.phys.uconn.edu/~mokaya/gradphi.x]. | *Edit events.in[http://zeus.phys.uconn.edu/~mokaya/events.in] and control.in[http://zeus.phys.uconn.edu/~mokaya/control.in] files and the run gradphi.x[http://zeus.phys.uconn.edu/~mokaya/gradphi.x]. | ||
| − | + | ===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. | * 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"> | <table border ="5"> | ||
<tr> | <tr> | ||
| Line 18: | Line 17: | ||
<td>Output file size (MB)</td> | <td>Output file size (MB)</td> | ||
<td> Total time in days</td> | <td> Total time in days</td> | ||
| − | + | </tr> <tr> | |
| − | |||
<td>100</td> | <td>100</td> | ||
<td>0.0216</td> | <td>0.0216</td> | ||
<td>59000 </td> | <td>59000 </td> | ||
<td> 25 </td> | <td> 25 </td> | ||
| − | + | </tr> | |
| − | |||
</table> | </table> | ||
| + | |||
| + | ===Results of the simulation=== | ||
| + | * γ p -> p (ω -> π<sup>0</sup>γ) | ||
| + | * γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ) | ||
Revision as of 10:58, 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γ)