Difference between revisions of "Analysis of evio files"

From New IAC Wiki
Jump to navigation Jump to search
Line 8: Line 8:
  
 
=GEMC 2.4=
 
=GEMC 2.4=
<pre>
 
~/src/CLAS/coatjava-2.4/bin/evio-dump -i LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.evio
 
</pre>
 
 
 
Looking at the evio dump, we can see that one electron is detected by several wires, several times
 
 
<pre>
 
*********************** 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 
 
</pre>
 
 
  
 
Restricting hits for Superlayer 1, Layer 1 only
 
Restricting hits for Superlayer 1, Layer 1 only
Line 80: Line 21:
  
 
[[File:Layer1BinHits.png]][[File:Layer1BinHitsWeighted.png]]
 
[[File:Layer1BinHits.png]][[File:Layer1BinHitsWeighted.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.
 +
[[File:1stHitsLayer1binHits.png]][[File:1stHitsLayer1binHitsWeighted.png]]
  
  
Line 98: Line 44:
  
 
[[File:ThetaBins.png]][[File:ThetaBinsWeighted.png]]
 
[[File:ThetaBins.png]][[File:ThetaBinsWeighted.png]]
 +
 +
 +
Again, restricting hits to only the 1st hits within layer 1
 +
[[File:1stHitsThetaBins.png]][[File:1stHitsThetaBinsWeighted.png]]
  
 
=GEMC 3.0=
 
=GEMC 3.0=

Revision as of 18:24, 31 August 2017

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 99 to 999

ProcessesGEMC3.png



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

Layer1binsHitsGEMC3.pngLayer1binsHitsWeightedGEMC3.png


This lack of interactions translates to angle bins as well

ThetaBinsGEMC3.pngThetaBinsWeightedGEMC3.png