Benchmark GEMC Occupancy Prediction for 11GeV Electron Beam with 0T Solenoid for LH2

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

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

Links

Back