Difference between revisions of "Benchmark GEMC Occupancy Prediction for 11GeV Electron Beam with 0T Solenoid for LH2"
Line 28: | Line 28: | ||
[[DV_RunGroupC_Moller#Benchmark_GEMC_Occupancy_Prediction_for_11GeV_Electron_Beam_with_0T_Solenoid_for_LH2|Back]] | [[DV_RunGroupC_Moller#Benchmark_GEMC_Occupancy_Prediction_for_11GeV_Electron_Beam_with_0T_Solenoid_for_LH2|Back]] | ||
+ | |||
+ | Using the evio2rootAnalysis.C | ||
+ | |||
+ | Checking the output versus the LUND_Spread_LH2_00.LUND | ||
+ | |||
+ | <pre> | ||
+ | Processing evio2rootAnalysis.C... | ||
+ | EVENT # 0 GENPART = 2Detector | ||
+ | Gen part 0 pid = 11 px = -0.0113655 py = -0.00200404 pz = 10.8681 vx = 0 vy = 0 vz = 0 | ||
+ | Gen part 1 pid = 11 px = 0.0113655 py = 0.00200404 pz = 0.131912 vx = 0 vy = 0 vz = 0 | ||
+ | Detector 0 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 2 wire = 1 | ||
+ | Detector 1 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 1 wire = 1 | ||
+ | Detector 2 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 2 wire = 1 | ||
+ | EVENT # 1 GENPART = 2Detector | ||
+ | Gen part 0 pid = 11 px = -0.0113654 py = -0.00200403 pz = 10.8681 vx = 0 vy = 0 vz = 0 | ||
+ | Gen part 1 pid = 11 px = 0.0113654 py = 0.00200403 pz = 0.13191 vx = 0 vy = 0 vz = 0 | ||
+ | Detector 0 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 1 wire = 1 | ||
+ | Detector 1 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 2 wire = 1 | ||
+ | Detector 2 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 4 wire = 1 | ||
+ | Detector 3 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 1 wire = 1 | ||
+ | Detector 4 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 2 wire = 1 | ||
+ | Detector 5 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 3 wire = 1 | ||
+ | Detector 6 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 4 wire = 1 | ||
+ | </pre> | ||
+ | |||
+ | <pre> | ||
+ | 2 1 1 1 1 0.006411440670 -81.074954045618540 -20140.877319449689821 1.625975760324849 -10735.677605378126827 | ||
+ | 1 -1 1 11 0 0 -0.011365497628610 -0.002004043881097 10.868088162529032 10.868094302146620 0.000511 0 0 0 | ||
+ | 2 -1 1 11 0 0 0.011365497628610 0.002004043881097 0.131911825606992 0.132416696768492 0.000511 0 0 0 | ||
+ | 2 1 1 1 1 0.006411635317 -81.076265062091281 -20140.885108455677255 1.625924392583043 -10735.681784681821227 | ||
+ | 1 -1 1 11 0 0 -0.011365408694758 -0.002004028199659 10.868090252277955 10.868096391798467 0.000511 0 0 0 | ||
+ | 2 -1 1 11 0 0 0.011365408694758 0.002004028199659 0.131909735858069 0.132414607116645 0.000511 0 0 0 | ||
+ | </pre> |
Revision as of 18:35, 3 August 2016
Splitting an Isotropic distribution for an 11GeV Electron bean at 0T for the Solenoid field with an LH2 target, we start with 5 files LUND_Spread_00 - LUND_Spread_04. LUND_Spread_00 - LUND_Spread_03 are equally divided with 75000 lines or 25000 Moller events in each file. LUND_Spread_04 has 37731 lines or 12577 Moller events. These LUND files are run in GEMC with:
~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,LUND_Spread_LH2_00.LUND" -N=75000 eg12.gcard > tracking.dat
Using the evio file created from the GEMC run:
~/src/CLAS/evioRoot/bin/evio2root eg12.evio eg12.root 25000
LUND_Spread_02 has too many events with respect to the 75000 events GEMC can handle. Initially, LUND_Spread_02 had 25000 Moller events or 75000 lines(Header, Scattered, and Moller) in the file. Do to the scattering and secondary particles, this file is divided into two parts LUND_Spread_02_part_00 and LUND_Spread_02_part_01 of 37500 lines or 12500 Moller events.
The second file is also two large and is similarly split into two parts LUND_Spread_02_part_01_00 and LUND_Spread_02_part_01_01 of 18750 lines or 6250 Moller events.
The root files can be combined with:
hadd 0T_LH2.root ./LUND_Spread_00/eg12.root LUND_Spread_01/eg12.root LUND_Spread_02/eg12_00.evio LUND_Spread_02/eg12_01_00.evio LUND_Spread_02/eg12_01_01.evio LUND_Spread_03/eg12.root LUND_Spread_04/eg12.root
Links
Using the evio2rootAnalysis.C
Checking the output versus the LUND_Spread_LH2_00.LUND
Processing evio2rootAnalysis.C... EVENT # 0 GENPART = 2Detector Gen part 0 pid = 11 px = -0.0113655 py = -0.00200404 pz = 10.8681 vx = 0 vy = 0 vz = 0 Gen part 1 pid = 11 px = 0.0113655 py = 0.00200404 pz = 0.131912 vx = 0 vy = 0 vz = 0 Detector 0 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 2 wire = 1 Detector 1 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 1 wire = 1 Detector 2 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 2 wire = 1 EVENT # 1 GENPART = 2Detector Gen part 0 pid = 11 px = -0.0113654 py = -0.00200403 pz = 10.8681 vx = 0 vy = 0 vz = 0 Gen part 1 pid = 11 px = 0.0113654 py = 0.00200403 pz = 0.13191 vx = 0 vy = 0 vz = 0 Detector 0 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 1 wire = 1 Detector 1 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 2 wire = 1 Detector 2 dpid = 11 procID = 90 sector = 1 superlayer = 1 layer = 4 wire = 1 Detector 3 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 1 wire = 1 Detector 4 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 2 wire = 1 Detector 5 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 3 wire = 1 Detector 6 dpid = 11 procID = 90 sector = 1 superlayer = 2 layer = 4 wire = 1
2 1 1 1 1 0.006411440670 -81.074954045618540 -20140.877319449689821 1.625975760324849 -10735.677605378126827 1 -1 1 11 0 0 -0.011365497628610 -0.002004043881097 10.868088162529032 10.868094302146620 0.000511 0 0 0 2 -1 1 11 0 0 0.011365497628610 0.002004043881097 0.131911825606992 0.132416696768492 0.000511 0 0 0 2 1 1 1 1 0.006411635317 -81.076265062091281 -20140.885108455677255 1.625924392583043 -10735.681784681821227 1 -1 1 11 0 0 -0.011365408694758 -0.002004028199659 10.868090252277955 10.868096391798467 0.000511 0 0 0 2 -1 1 11 0 0 0.011365408694758 0.002004028199659 0.131909735858069 0.132414607116645 0.000511 0 0 0