Difference between revisions of "Analysis of evio files"
Jump to navigation
Jump to search
| Line 103: | Line 103: | ||
[[File:ProcessesGEMC3.png]] | [[File:ProcessesGEMC3.png]] | ||
| + | |||
| + | |||
GEMC 3.0 seems to not have as many interactions on layer 1 of section 1 superlayer 1 | GEMC 3.0 seems to not have as many interactions on layer 1 of section 1 superlayer 1 | ||
| + | |||
[[File:Layer1binsHitsGEMC3.png]][[File:Layer1binsHitsWeightedGEMC3.png]] | [[File:Layer1binsHitsGEMC3.png]][[File:Layer1binsHitsWeightedGEMC3.png]] | ||
| Line 111: | Line 114: | ||
This lack of interactions translates to angle bins as well | This lack of interactions translates to angle bins as well | ||
| + | |||
[[File:ThetaBinsGEMC3.png]][[File:ThetaBinsWeightedGEMC3.png]] | [[File:ThetaBinsGEMC3.png]][[File:ThetaBinsWeightedGEMC3.png]] | ||
Revision as of 21:28, 30 August 2017
Using the evio2root converter
~/src/CLAS/evio2rootDV/bin/evio2root LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut 105
GEMC 2.4
~/src/CLAS/coatjava-2.4/bin/evio-dump -i LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.evio
Looking at the evio dump, we can see that one electron is detected by several wires, several times
*********************** EVENT # 7 ***********************
+------------------------------------------------+------------+------------+
| bank| nrows| ncols|
+------------------------------------------------+------------+------------+
| DC::dgtz| 27| 11|
| DC::true| 27| 24|
| GenPart::true| 2| 7|
+------------------------------------------------+------------+------------+
Press Enter for Next Event or Bank Name: DC::dgtz
*****>>>>> BANK DC::dgtz >>>> SIZE = 11
sector (int) : 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
superlayer (int) : 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1
layer (int) : 1 2 3 4 5 6 1 2 3 4 4 5 6 6 6 5 5 4 3 2 1 6 5 4 3 2 1
wire (int) : 47 48 47 48 47 48 50 51 51 51 52 51 52 60 61 60 61 61 61 62 61 61 60 61 60 61 60
LR (int) : 1 -1 1 -1 1 -1 1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 -1 1 -1 1 -1 1 -1 1
tdc (int) : Exception in thread "main" java.lang.NullPointerException
at org.jlab.data.io.BasicDataBank.show(BasicDataBank.java:356)
at org.jlab.clasrec.rec.EvioDump.main(EvioDump.java:75)
Press Enter for Next Event or Bank Name: DC::true
*****>>>>> BANK DC::true >>>> SIZE = 24
pid (int) : 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11
mpid (int) : 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
tid (int) : 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
mtid (int) : 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
otid (int) : 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
hitn (int) : 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27
trackE (double) : 5.35903 5.35723 5.35324 5.35107 5.34911 5.34711 5.33412 5.33094 5.32836 5.32684 5.32177 5.32119 5.31732 5.07671 5.07648 5.07334 5.07157 5.07076 5.06764 5.06206 5.05933 5.05164 5.04940 5.04743 5.04519 5.04133 5.03954
totEdep (double) : 0.00178 0.00396 0.00210 0.00195 0.00195 0.00264 0.00315 0.00248 0.00152 0.00506 0.00059 0.00385 0.00931 0.00024 0.00314 0.00177 0.00073 0.00251 0.00558 0.00272 0.00410 0.00218 0.00197 0.00220 0.00382 0.00172 0.00284
avgX (double) : 794.10983 797.17733 802.45218 807.35462 812.87136 818.87821 842.85191 848.07336 857.10227 865.11391 868.14070 873.60264 882.39062 995.97798 995.51549 994.79675 994.16597 993.37823 991.14933 988.25963 985.34151 969.95437 964.78943 960.23769 955.72371 948.55600 942.78265
avgY (double) : 162.76457 162.20898 161.15242 159.88527 158.50215 157.10496 151.84826 151.08776 149.78433 148.70277 148.44552 148.02561 147.61420 113.90546 112.62839 109.72058 107.67506 105.68644 101.27035 97.24497 93.37926 79.47363 75.47866 72.19536 69.12653 64.60327 61.31620
avgZ (double) : 2161.06397 2168.39078 2180.10583 2190.29767 2201.15050 2212.25981 2254.21993 2262.31716 2275.10088 2285.04437 2288.59301 2294.53303 2303.46285 2255.06346 2249.45389 2238.58371 2231.30695 2224.69130 2210.31875 2197.39957 2184.94882 2139.92904 2129.18428 2119.75783 2110.99867 2098.76854 2089.41744
avgLx (double) : -125.16538 -125.31973 -123.21014 -122.70399 -120.33064 -119.75937 -188.10069 -186.43297 -185.93573 -183.76592 -183.38010 -183.79703 -182.17539 -135.75253 -134.27846 -131.74892 -129.45291 -126.16799 -122.12677 -116.73682 -113.39315 -60.06076 -55.33306 -53.09397 -49.26141 -45.66517 -41.52230
avgLy (double) : -126.28963 -119.85697 -126.60508 -119.64840 -125.78195 -118.20176 -84.81801 -76.59357 -80.96566 -71.03576 -69.82598 -74.44436 -63.31871 55.88135 57.68871 54.30481 56.58089 65.44454 62.01375 71.41930 66.61718 56.48520 50.07427 56.96588 50.21041 56.05215 48.43307
avgLz (double) : 1.81252 -1.83675 -0.57606 -0.85325 -0.27178 0.74928 0.19779 -2.38295 0.89282 1.16455 5.65988 1.22568 0.90681 5.04616 -0.23330 1.73723 -5.12434 0.67293 -1.16900 -1.97299 -2.36446 -0.95711 -1.29197 -0.17289 1.56693 -0.96054 -0.28944
px (double) : 2.00485 2.03220 2.18854 2.29291 2.36567 2.50810 2.71171 2.87839 3.05683 3.28934 3.47055 3.48460 3.71766 -0.36708 -0.38977 -0.33289 -0.43757 -0.50413 -0.71220 -1.08034 -1.05085 -1.81485 -2.12195 -2.01850 -2.18842 -2.43490 -2.54570
py (double) : -0.34518 -0.39116 -0.43566 -0.61155 -0.62272 -0.57946 -0.40833 -0.42691 -0.45366 -0.42929 -0.35057 -0.28296 -0.28643 -1.01900 -1.01203 -1.32990 -1.38436 -1.36427 -1.45018 -1.49218 -1.43464 -1.59224 -1.63474 -1.50370 -1.47441 -1.57627 -1.44044
pz (double) : 4.93148 4.91487 4.83907 4.76880 4.72945 4.65877 4.54661 4.43739 4.31049 4.13642 3.98655 3.97887 3.75629 -4.93343 -4.93289 -4.85780 -4.83237 -4.83075 -4.77595 -4.68717 -4.70885 -4.40784 -4.24974 -4.34511 -4.26963 -4.09151 -4.07191
vx (double) : 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
vy (double) : 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
vz (double) : 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
mvx (double) : 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
mvy (double) : 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
mvz (double) : 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
avgT (double) : 7.74108 7.76844 7.81097 7.84911 7.89017 7.93275 8.09568 8.12807 8.18071 8.22554 8.23932 8.26718 8.31502 10.97539 10.99475 11.03256 11.05799 11.08131 11.13277 11.17925 11.22388 11.38987 11.43204 11.46883 11.50368 11.55334 11.59211
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.
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.
GEMC 3.0
GEMC 3.0 has changed the particle ID for transportation from 99 to 999
GEMC 3.0 seems to not have as many interactions on layer 1 of section 1 superlayer 1
This lack of interactions translates to angle bins as well
