Analysis of evio files
Using the evio2root converter
~/src/CLAS/evio2rootDV/bin/evio2root LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut 105
GEMC 2.4
Restricting hits for Superlayer 1, Layer 1 only
The bumps in what would otherwise be an isotropic distribution of the theta and phi values for particles is caused by two factors.
The foremost factor is that the detector has limits of near the 20degrees in phi as theta reaches 5degrees. Angles outside this limiting factor are provided in the LUND file, but if an interaction where to occur, it would not be in sector 1. These angles would correspond to neighboring sectors.
The second cause is that angles of theta outside the range 5-40 degrees would contribute to wire interactions on the left and right hand sides of the ellipse marked out on the detector plane.
Additionally, it is possible that more than one wire will register a hit. We can limit the hits by only accepting the dnhits for the first hit.
The hits on DC Sector 1 Superlayer 1 Layer 1 have an artificial bump near the middle due to the decreasing wire number as phi (or the parameter of the ellipse) moves away from its semi-major radii with respect to the ellipse of constant theta.
We can run a larger simulation or we can use the wire hit and phi angle to find theta within a certain range. This information can be used to make a plot of number of hits per angle theta bin on which the differential cross section can be overlayed. Theoretically the phi data would be obtained from additional layers in the sector as well as timing data.
Again, restricting hits to only the 1st hits within layer 1
GEMC 3.0
GEMC 3.0 has changed the particle ID for transportation from 90 to 999
GEMC 3.0 seems to not have as many interactions on layer 1 of section 1 superlayer 1
Restricting results to only the 1st hits on this layer
This lack of interactions translates to angle bins as well
Restricting results to only the 1st hits on this layer
evio2rootAnalysis
GEMC2.4
--- Header Bank --
- (10, 1) time: Fri, 09.08.2017 13:28:00
- (10, 5) beamPol: 1
- (10, 3) evn: 1
- (10, 4) evn_type: -1
- (10, 2) runNo: 1
- (10, 6) var1: 1
- (10, 7) var2: 1
- (10, 8) var3: 1
- (10, 9) var4: 0.00641143
- (10, 10) var5: -81.075
- (10, 11) var6: -20140.9
- (10, 12) var7: 1.62598
- (10, 13) var8: -10735.7
--- End of Header Bank --
--- Generated Particles Bank --
- Particle 1 pid: 11 - mom: (-9.99465,5.77041,10868.1) MeV - vert: (0,0,0) mm
- Particle 2 pid: 11 - mom: (9.99465,-5.77041,131.912) MeV - vert: (0,0,0) mm
--- End of Generated Particles Bank --
GEMC3.0
--- Header Bank --
- (10, 1) time: Fri, 09.08.2017 14:09:27
- (10, 5) beamPol: 1
- (10, 3) evn: 1
- (10, 4) evn_type: -1
- (10, 2) runNo: 1
- (10, 5) user var 1: 2
- (10, 6) user var 10: 1.62598
- (10, 7) user var 11: -10735.7
- (10, 8) user var 2: 1
- (10, 9) user var 3: 1
- (10, 10) user var 4: 1
- (10, 11) user var 5: 1
- (10, 12) user var 6: 1
- (10, 13) user var 7: 0.00641143
- (10, 14) user var 8: -81.075
- (10, 15) user var 9: -20140.9
--- End of Header Bank --
--- Generated Particles Bank --
- Particle 1 pid: 11 - mom: (-9.99465,5.77041,10868.1) MeV - vert: (0,0,0) mm
- Particle 2 pid: 11 - mom: (9.99465,-5.77041,131.912) MeV - vert: (0,0,0) mm
--- Generated User Info Particles Bank --
1 -1 1 11 0 0 -0.00999465 0.00577041 10.8681 10.8681 0.000511 0 0 0
2 -1 1 11 0 0 0.00999465 -0.00577041 0.131912 0.132417 0.000511 0 0 0
--- End of Generated Particles Bank --
GEMC2.4
EVENT # 1 GENPART = 2Detector
Gen part 0 pid = 11 px = -4.16454e-17 py = 0.190174 pz = 5.67653e-31 vx = 0 vy = 0 vz = 0
Gen part 1 pid = 11 px = -1.03524 py = -0.966768 pz = 2.58963 vx = 0 vy = 0 vz = 0
GEMC3.0
EVENT # 1 GENPART = 2Detector
Gen part 0 pid = 11 px = -4.16454e-17 py = 0.190174 pz = 5.67653e-31 vx = 0 vy = 0 vz = 0
Gen part 1 pid = 11 px = -1.03524 py = 0.966768 pz = 2.58963 vx = 0 vy = 0 vz = 0
Detector 0 dpid = -11 procID = 5 sector = 1 superlayer = 1 layer = 49 wire = 138
Detector 1 dpid = -11 procID = 5 sector = 1 superlayer = 1 layer = 50 wire = 67
Detector 2 dpid = -11 procID = 5 sector = 1 superlayer = 2 layer = 51 wire = 7
Detector 3 dpid = -11 procID = 5 sector = 1 superlayer = 3 layer = 51 wire = 55
Detector 4 dpid = -11 procID = 5 sector = 1 superlayer = 3 layer = 52 wire = 158
Detector 5 dpid = -11 procID = 5 sector = 1 superlayer = 4 layer = 52 wire = 118
Detector 6 dpid = -11 procID = 5 sector = 1 superlayer = 4 layer = 53 wire = 81
Detector 7 dpid = -11 procID = 5 sector = 1 superlayer = 5 layer = 53 wire = 9
Detector 8 dpid = -11 procID = 5 sector = 1 superlayer = 6 layer = 54 wire = 58
Detector 9 dpid = -11 procID = 5 sector = 1 superlayer = 6 layer = 55 wire = 145
Detector 10 dpid = -11 procID = 5 sector = 1 superlayer = 7 layer = 63 wire = 130
Detector 11 dpid = -11 procID = 5 sector = 1 superlayer = 7 layer = 64 wire = 54
Detector 12 dpid = -11 procID = 5 sector = 1 superlayer = 8 layer = 65 wire = 6
Detector 13 dpid = -11 procID = 5 sector = 1 superlayer = 9 layer = 65 wire = 54
Detector 14 dpid = -11 procID = 5 sector = 1 superlayer = 9 layer = 66 wire = 179
Detector 15 dpid = -11 procID = 5 sector = 1 superlayer = 10 layer = 66 wire = 114
Detector 16 dpid = -11 procID = 5 sector = 1 superlayer = 10 layer = 67 wire = 110
Detector 17 dpid = -11 procID = 5 sector = 1 superlayer = 11 layer = 66 wire = 177
Detector 18 dpid = -11 procID = 5 sector = 1 superlayer = 11 layer = 67 wire = 47
Detector 19 dpid = -11 procID = 5 sector = 1 superlayer = 12 layer = 68 wire = 13
Detector 20 dpid = -11 procID = 5 sector = 1 superlayer = 15 layer = 91 wire = 142
Detector 21 dpid = -11 procID = 5 sector = 1 superlayer = 16 layer = 92 wire = 60
Detector 22 dpid = -11 procID = 5 sector = 1 superlayer = 17 layer = 92 wire = 6
Detector 23 dpid = -11 procID = 5 sector = 1 superlayer = 18 layer = 93 wire = 79
Detector 24 dpid = -11 procID = 5 sector = 1 superlayer = 18 layer = 94 wire = 308
Detector 25 dpid = -11 procID = 5 sector = 1 superlayer = 19 layer = 106 wire = 226
Detector 26 dpid = -11 procID = 5 sector = 1 superlayer = 19 layer = 107 wire = 168
Detector 27 dpid = -11 procID = 5 sector = 1 superlayer = 20 layer = 108 wire = 15
Detector 28 dpid = -11 procID = 5 sector = 1 superlayer = 21 layer = 108 wire = 140
Detector 29 dpid = -11 procID = 5 sector = 1 superlayer = 21 layer = 109 wire = 226
Detector 30 dpid = -11 procID = 5 sector = 1 superlayer = 22 layer = 109 wire = 319
Detector 31 dpid = -11 procID = 5 sector = 1 superlayer = 22 layer = 110 wire = 43
Detector 32 dpid = -11 procID = 5 sector = 1 superlayer = 23 layer = 110 wire = 136
Detector 33 dpid = -11 procID = 5 sector = 1 superlayer = 23 layer = 111 wire = 238
Detector 34 dpid = -11 procID = 5 sector = 1 superlayer = 24 layer = 111 wire = 315
Detector 35 dpid = -11 procID = 5 sector = 1 superlayer = 24 layer = 112 wire = 41
GEMC2.4
EVENT # 38 GENPART = 2Detector
Gen part 0 pid = 11 px = -2.49137e+14 py = -1.14941e-09 pz = 5.67653e-31 vx = 0 vy = 0 vz = 0
Gen part 1 pid = 11 px = -1.04572 py = -0.945366 pz = 2.58963 vx = 0 vy = 0 vz = 0
Detector 0 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 1 wire = 1
Detector 1 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 2 wire = 2
Detector 2 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 2 wire = 1
Detector 3 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 3 wire = 1
Detector 4 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 4 wire = 1
Detector 5 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 5 wire = 1
Detector 6 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 6 wire = 1
GEMC3.0
EVENT # 38 GENPART = 2Detector
Gen part 0 pid = 11 px = -2.49137e+14 py = -1.14941e-09 pz = 5.67653e-31 vx = 0 vy = 0 vz = 0
Gen part 1 pid = 11 px = -1.04572 py = 0.945366 pz = 2.58963 vx = 0 vy = 0 vz = 0
Detector 0 dpid = -11 procID = 5 sector = 1 superlayer = 1 layer = 11 wire = 30
Detector 1 dpid = -11 procID = 5 sector = 1 superlayer = 2 layer = 12 wire = 134
Detector 2 dpid = -11 procID = 5 sector = 1 superlayer = 3 layer = 11 wire = 17
Detector 3 dpid = -11 procID = 5 sector = 1 superlayer = 4 layer = 12 wire = 149
Detector 4 dpid = -11 procID = 5 sector = 1 superlayer = 5 layer = 11 wire = 6
Detector 5 dpid = -11 procID = 5 sector = 1 superlayer = 6 layer = 11 wire = 152
Detector 6 dpid = -11 procID = 5 sector = 1 superlayer = 7 layer = 14 wire = 7
Detector 7 dpid = -11 procID = 5 sector = 1 superlayer = 8 layer = 15 wire = 170
Detector 8 dpid = -11 procID = 5 sector = 1 superlayer = 9 layer = 14 wire = 7
Detector 9 dpid = -11 procID = 5 sector = 1 superlayer = 10 layer = 15 wire = 154
Detector 10 dpid = -11 procID = 5 sector = 1 superlayer = 11 layer = 14 wire = 11
Detector 11 dpid = -11 procID = 5 sector = 1 superlayer = 12 layer = 15 wire = 154
Detector 12 dpid = -11 procID = 5 sector = 1 superlayer = 13 layer = 26 wire = 11
Detector 13 dpid = -11 procID = 5 sector = 1 superlayer = 14 layer = 27 wire = 297
Detector 14 dpid = -11 procID = 5 sector = 1 superlayer = 15 layer = 26 wire = 7
Detector 15 dpid = -11 procID = 5 sector = 1 superlayer = 16 layer = 27 wire = 277
Detector 16 dpid = -11 procID = 5 sector = 1 superlayer = 17 layer = 26 wire = 8
Detector 17 dpid = -11 procID = 5 sector = 1 superlayer = 18 layer = 27 wire = 288
Detector 18 dpid = -11 procID = 5 sector = 1 superlayer = 19 layer = 29 wire = 176
Detector 19 dpid = -11 procID = 5 sector = 1 superlayer = 20 layer = 30 wire = 82
Detector 20 dpid = -11 procID = 5 sector = 1 superlayer = 21 layer = 29 wire = 200
Detector 21 dpid = -11 procID = 5 sector = 1 superlayer = 22 layer = 30 wire = 51
Detector 22 dpid = -11 procID = 5 sector = 1 superlayer = 23 layer = 29 wire = 222
Detector 23 dpid = -11 procID = 5 sector = 1 superlayer = 24 layer = 30 wire = 63
Detector 24 dpid = -11 procID = 5 sector = 1 superlayer = 25 layer = 34 wire = 553
Detector 25 dpid = -11 procID = 5 sector = 1 superlayer = 25 layer = 35 wire = 551
Detector 26 dpid = -11 procID = 5 sector = 1 superlayer = 26 layer = 35 wire = 31
Detector 27 dpid = -11 procID = 5 sector = 1 superlayer = 27 layer = 35 wire = 445
Detector 28 dpid = -11 procID = 5 sector = 1 superlayer = 28 layer = 35 wire = 112
Detector 29 dpid = -11 procID = 5 sector = 1 superlayer = 29 layer = 35 wire = 310
Detector 30 dpid = -11 procID = 5 sector = 1 superlayer = 30 layer = 35 wire = 198
Detector 31 dpid = -11 procID = 5 sector = 1 superlayer = 31 layer = 37 wire = 153
Detector 32 dpid = -11 procID = 5 sector = 1 superlayer = 32 layer = 38 wire = 265
Detector 33 dpid = -11 procID = 5 sector = 1 superlayer = 33 layer = 37 wire = 280
Detector 34 dpid = -11 procID = 5 sector = 1 superlayer = 34 layer = 38 wire = 125
Detector 35 dpid = -11 procID = 5 sector = 1 superlayer = 35 layer = 37 wire = 449
Detector 36 dpid = -11 procID = 5 sector = 1 superlayer = 36 layer = 38 wire = 34
Detector 37 dpid = 11 procID = 1 sector = 1 superlayer = 35 layer = 38 wire = 403
GEMC 2.4 gcardFile:Eg12GEMC2 4gcard.txt
GEMC 3.0 gcardFile:Eg12GEMC2 6gcard.txt