Analysis of evio files

From New IAC Wiki
Jump to navigation Jump to search

Using the evio2root converter

File:Evio2root.c

~/src/CLAS/evio2rootDV/bin/evio2root LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut 105

GEMC 2.4

Restricting hits for Superlayer 1, Layer 1 only

Layer1WireNumberTotal.pngLayer1WireNumberBinTotalWeighted.png

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.

Layer1BinHits.pngLayer1BinHitsWeighted.png


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.

1stHitsLayer1binHits.png1stHitsLayer1binHitsWeighted.png


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.

Layer1WireNumberBin0.pngLayer1WireNumberBin0Weighted.png


Layer1WireNumberBin150.pngLayer1WireNumberBin150Weighted.png


Layer1WireNumberBin299.pngLayer1WireNumberBin299Weighted.png


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.


ThetaBins.pngThetaBinsWeighted.png


Again, restricting hits to only the 1st hits within layer 1

1stHitsThetaBins.png1stHitsThetaBinsWeighted.png

GEMC 3.0

GEMC 3.0 has changed the particle ID for transportation from 90 to 999

ProcessesGEMC3.png



GEMC 3.0 seems to not have as many interactions on layer 1 of section 1 superlayer 1

Layer1BinHitsGEMC3.pngLayer1BinHitsWeightedGEMC3.png


Restricting results to only the 1st hits on this layer

1stHitsLayer1binsHitsGEMC3.png1stHitsLayer1binsHitsWeightedGEMC3.png

This lack of interactions translates to angle bins as well

ThetaBinsGEMC3.pngThetaBinsWeightedGEMC3.png

Restricting results to only the 1st hits on this layer

1stHitsThetaBinsGEMC3.png1stHitsThetaBinsWeightedGEMC3.png

evio2rootAnalysis

GEMC2.4
 EVENT # 0  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = 1.29811e+13 py  = -1.51141e+10 pz  = 5.67653e-31 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = -1.03491 py  = -0.967327 pz  = 2.58963 vx  = 0 vy  = 0 vz  = 0
          Detector 0 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 3 layer  = 1 wire  = 11
          Detector 1 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 3 layer  = 2 wire  = 11
          Detector 2 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 3 layer  = 2 wire  = 10
          Detector 3 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 3 layer  = 2 wire  = 9
          Detector 4 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 3 layer  = 1 wire  = 8
          Detector 5 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 3 layer  = 1 wire  = 7

GEMC3.0
EVENT # 0  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = 1.29811e+13 py  = -1.51141e+10 pz  = 5.67653e-31 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = -1.03491 py  = 0.967327 pz  = 2.58963 vx  = 0 vy  = 0 vz  = 0


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:Eg12GEMC1gcard.txt


GEMC 3.0 gcardFile:Eg12GEMC3gcard.txt