Difference between revisions of "Radphi MonteCarlo"
| (2 intermediate revisions by one other user not shown) | |||
| Line 50: | Line 50: | ||
To make this change edit the file decay.dat for each reaction of interest. Then the generation and simulation proceed as before. | To make this change edit the file decay.dat for each reaction of interest. Then the generation and simulation proceed as before. | ||
| − | + | These data files are stored in the group dcache area under /pnfs/phys.uconn.edu/data/Gluex/radphi/omega_sim-6-2014. | |
| − | |||
# 10<sup>8</sup> events of γ p -> p (ω -> π<sup>0</sup>γ) | # 10<sup>8</sup> events of γ p -> p (ω -> π<sup>0</sup>γ) | ||
#* 100 files of 1000000 events each | #* 100 files of 1000000 events each | ||
| − | + | #* files are named omega2p_sim_<i>nnn</i> where <i>nnn</i> = 1..100 | |
# 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup> π<sup>+</sup> π<sup>-</sup>) | # 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup> π<sup>+</sup> π<sup>-</sup>) | ||
#* 100 files of 1000000 events each | #* 100 files of 1000000 events each | ||
| − | + | #* files are named omega3D+p_sim_<i>nnn</i> where <i>nnn</i> = 101..200 | |
# 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ) | # 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> p π<sup>0</sup>) (ω -> π<sup>0</sup>γ) | ||
#* 100 files of 1000000 events each | #* 100 files of 1000000 events each | ||
| − | + | #* files are named omega2D+p_sim_<i>nnn</i> where <i>nnn</i> = 201..300 | |
# 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> n π<sup>+</sup>) (ω -> π<sup>0</sup>γ) | # 10<sup>8</sup> events of γ p -> (Δ<sup>+</sup> -> n π<sup>+</sup>) (ω -> π<sup>0</sup>γ) | ||
#* 100 files of 1000000 events each | #* 100 files of 1000000 events each | ||
| − | + | #* files are named omega2D+n_sim_<i>nnn</i> where <i>nnn</i> = 301..400 | |
# 10<sup>8</sup> events of γ n -> (Δ<sup>0</sup> -> p π<sup>-</sup>) (ω -> π<sup>0</sup>γ) | # 10<sup>8</sup> events of γ n -> (Δ<sup>0</sup> -> p π<sup>-</sup>) (ω -> π<sup>0</sup>γ) | ||
#* 100 files of 1000000 events each | #* 100 files of 1000000 events each | ||
| + | #* files are named omega2D0p_sim_<i>nnn</i> where <i>nnn</i> = 401..500 | ||
Latest revision as of 13:36, 22 September 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
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. Files in this directory are either raw simulation output or reconstructed simulation results in the form of root trees. Files of the former type have the .itape extension and those of the latter type have the .root extension. The file names indicate the simulated channel, as given in the list below. The events below were generated using a slope parameter of 8.0/GeV2. These data files are stored in the group dcache area under /pnfs/phys.uconn.edu/data/Gluex/radphi/omega_sim-6-2014.
- 108 events of γ p -> p (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2p_sim_nnn where nnn = 1..100
- 108 events of γ p -> (Δ+ -> p π0) (ω -> π0 π+ π-)
- 100 files of 1000000 events each
- files are named omega3D+p_sim_nnn where nnn = 101..200
- 108 events of γ p -> (Δ+ -> p π0) (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2D+p_sim_nnn where nnn = 201..300
- 108 events of γ p -> (Δ+ -> n π+) (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2D+n_sim_nnn where nnn = 301..400
- 108 events of γ n -> (Δ0 -> p π-) (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2D0p_sim_nnn where nnn = 401..500
Monte Carlo Events with Slope Parameter b=5.5/GeV2
Monte Carlo events are generated with a new slope parameter of 5.5/GeV2. The previous results (shown above) were generated using a slope parameter of 8.0. 108 events are generated and then simulated in the detector. To make this change edit the file decay.dat for each reaction of interest. Then the generation and simulation proceed as before.
These data files are stored in the group dcache area under /pnfs/phys.uconn.edu/data/Gluex/radphi/omega_sim-6-2014.
- 108 events of γ p -> p (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2p_sim_nnn where nnn = 1..100
- 108 events of γ p -> (Δ+ -> p π0) (ω -> π0 π+ π-)
- 100 files of 1000000 events each
- files are named omega3D+p_sim_nnn where nnn = 101..200
- 108 events of γ p -> (Δ+ -> p π0) (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2D+p_sim_nnn where nnn = 201..300
- 108 events of γ p -> (Δ+ -> n π+) (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2D+n_sim_nnn where nnn = 301..400
- 108 events of γ n -> (Δ0 -> p π-) (ω -> π0γ)
- 100 files of 1000000 events each
- files are named omega2D0p_sim_nnn where nnn = 401..500