Analysis of evio files

From New IAC Wiki
Jump to navigation Jump to search

File:Evio2root.c

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

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


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

GEMC 3.0