Difference between revisions of "Analysis of evio files"

From New IAC Wiki
Jump to navigation Jump to search
 
(98 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 +
Using the evio2root converter
 +
 
[[File:Evio2root.c]]
 
[[File:Evio2root.c]]
  
 +
<pre>
 +
~/src/CLAS/evio2rootDV/bin/evio2root LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut 105
 +
</pre>
 +
 +
=GEMC 2.4=
 +
<pre>
 +
> Common Environment Version: <1.3>  (Thu March 31 2016)
 +
> Running as vanwdani on minerve2
 +
> OS Release:    Linux__Ubuntu14.04-x86_64-gcc4.8
 +
> JLAB_ROOT set to:    /home/vanwdani/src/CLAS/GEMC
 +
> JLAB_SOFTWARE set to: /home/vanwdani/src/CLAS/GEMC/1.3/Linux__Ubuntu14.04-x86_64-gcc4.8
 +
 +
> CLHEP  version: 2.2.0.8
 +
> XERCES  version: 3.1.1
 +
> QTDIR  version: 5.3.2
 +
> Geant4  version: 4.10.01.p03
 +
> ROOT    version: 5.34.34
 +
> GEMC    version: 2.3
 +
> JANA    version: 0.7.4p2
 +
> Build  version: 1.3
 +
> EVIO    version: 4.4.4
 +
> CCDB    version: 1.06
 +
> Banks  version: 1.2
 +
</pre>
 +
 +
 +
==Base Results==
 +
Restricting hits for Superlayer 1, Layer 1 only
 +
 +
[[File:Layer1WireNumberTotal.png]][[File:Layer1WireNumberBinTotalWeighted.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.
 +
 +
[[File:Layer1BinHits.png]][[File:Layer1BinHitsWeighted.png]]
 +
 +
 +
 +
Editing the eg12.gcard so that there is no energy cut and that energy hits of 0 are recorded only results in 103883 hits
 +
 +
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]]
 +
 +
 +
 +
Similarly, editing the gcard only results in 92761 1st hits via wires
 +
 +
These wire hits can be expressed in terms of Theta bins
 +
 +
[[File:ThetaBins.png]][[File:ThetaBinsWeighted.png]]
 +
 +
 +
With the edited gcard we only have 68114 hits per Theta bins
 +
 +
 +
Similarly, these can account for only the 1st interactions
 +
 +
[[File:1stHitsThetaBins.png]][[File:1stHitsThetaBinsWeighted.png]]
 +
 +
 +
With the edited gcard, we find 61102 hits
 +
 +
 +
Overlaying this with the expected differential cross-section expected from the LUND file
 +
 +
<center>[[File:EfficiencyThetaBins.png]]</center>
 +
 +
We can find the efficiency of the detector with respect to the particles input via the LUND file
 +
 +
<pre>
 +
TH1F *Ratio=new TH1F("TheoryExperiment","Theoretical and Experimental Differential Cross-Section Lab Frame",90,0,45);
 +
Ratio->Add(L1stHitsThetaBinsWeighted);
 +
Ratio->Divide(MolThetaLabweighted);
 +
Ratio->Draw();
 +
</pre>
 +
 +
<center>[[File:Ratio.png]]</center>
 +
 +
 +
This can be partially explained in that Superlayer 1, Layer 1 is not always hit first
 +
 +
[[File:LayersHitFirst.png]]
 +
 +
==Lower Fix==
 +
 +
Determining what values of phi are not within the sector following the Mathematica equation developed earlier
 +
 +
[[File:RejectedTheta.png]][[File:RejectedThetaWeighted.png]]
 +
 +
 +
 +
[[File:CorrectedMolThetaLab.png]][[File:CorrectedMolThetaLabWeighted.png]]
 +
 +
 +
[[File:EfficiencyThetaBinsCorrected.png]]
 +
 +
 +
 +
[[File:RatioCorrected.png]]
 +
 +
==Upper Fix==
 +
 +
Fixing the upper tail of hits on wires from larger theta values
 +
 +
 +
[[File:Layer1WireNumberBinTotalAdditional.png]] [[File:Layer1WireNumberBinTotalWeightedAdditional.png]]
 +
 +
 +
 +
 +
[[File:Layer1binsHitsAdditional.png]] [[File:Layer1binsHitsWeightedAdditional.png]]
 +
 +
 +
 +
 +
[[File:1stHitsLayer1binsHitsAdditional.png]] [[File:1stHitsLayer1binsHitsWeightedAdditional.png]]
 +
 +
==Combined Results==
 +
 +
Combining these plots
 +
<pre>
 +
hadd combined.root NonIsotropic_BinAdditions.root LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.root
 +
</pre>
 +
 +
[[File:Layer1WireNumberTotalCombined.png]][[File:Layer1WireNumberBinTotalWeightedCombined.png]]
 +
 +
 +
 +
[[File:Layer1binsHitsCombined.png]][[File:Layer1binsHitsWeightedCombined.png]]
 +
 +
 +
 +
[[File:1stHitsLayer1binsHitsCombined.png]][[File:1stHitsLayer1binsHitsWeightedCombined.png]]
 +
 +
=GEMC 2.6=
 +
 +
<pre>
 +
> Common Environment Version: <2.1>  (Wed January 25 2017)
 +
> Running as lds on lds3
 +
> OS Release:    Linux__Ubuntu16.04-x86_64-gcc5.4.0
 +
> JLAB_ROOT set to:    /home/lds/src/CLAS/GEMC
 +
> JLAB_SOFTWARE set to: /home/lds/src/CLAS/GEMC/2.1/Linux__Ubuntu16.04-x86_64-gcc5.4.0
 +
 +
> BANKS version:  1.3
 +
> CCDB version:  1.06.02
 +
> CLHEP version:  2.3.4.3
 +
> EVIO version:  5.1
 +
> GEANT4 version:  4.10.02.p03
 +
> GEMC version:  devel
 +
> JANA version:  0.7.7p1
 +
> MLIBRARY version:  1.1
 +
> MYSQL installed in  /home/lds/src/CLAS/GEMC/2.1/Linux__Ubuntu16.04-x86_64-gcc5.4.0/mysql/lib
 +
> QT  version:  5.8.0
 +
> ROOT version:  6.10.02
 +
> SCONS version:  1.5
 +
> XERCESC version:  3.1.4
 +
</pre>
 +
 +
GEMC 3.0 has changed the particle ID for transportation from 90 to 999
 +
 +
[[File:ProcessesGEMC3.png]]
 +
 +
==Base Results==
 +
 +
 +
GEMC 3.0 seems to use a new LUND format that has an additional variable on the header (11 variables instead of 10).  The databank for the DC has also eliminated the superlayer variable from the Digitialized output, moving up the wire number as a result.
 +
 +
[[File:Layer1BinHitsv2-6local.png]][[File:Layer1BinHitsWeightedGEMC3.png]]
 +
 +
 +
Restricting results to only the 1st hits on this layer
 +
 +
[[File:1stHitsLayer1binHitsGEMC3.png]][[File:1stHitsLayer1binsHitsWeightedGEMC3.png]]
 +
 +
This lack of interactions translates to angle bins as well
 +
 +
[[File:ThetaBinsGEMC3.png]][[File:ThetaBinsWeightedGEMC3.png]]
 +
 +
Restricting results to only the 1st hits on this layer
 +
 +
[[File:1stHitsThetaBinsGEMC3.png]][[File:1stHitsThetaBinsWeightedGEMC3.png]]
 +
 +
==Upper fix==
 +
 +
[[File:Layer1WireNumberBinTotaAdditionalv2-6.png]][[File:Layer1WireNumberBinTotalWeightedAdditionalv2-6.png]]
 +
 +
 +
 +
[[File:Layer1binsHitsAdditionalv2-6.png]][[File:Layer1binsHitsWeightedAdditionalv2-6.png]]
  
−
Looking at the evio dump, we can see that one electron is detected by several wires, several times
 
  
 +
 +
 +
[[File:1stHitsLayer1binsHitsAdditionalv2-6.png]][[File:1stHitsLayer1binsHitsWeightedAdditionalv2-6.png]]
 +
 +
 +
==Combined Results==
 +
 +
[[File:Layer1WireNumberTotalCombinedv2-6.png]][[File:Layer1WireNumberTotalWeightedCombinedv2-6.png]]
 +
 +
 +
 +
[[File:Layer1BinHitsCombinedv2-6.png]][[File:Layer1BinHitsWeightedCombinedv2-6.png]]
 +
 +
 +
 +
 +
[[File:1stHitsLayer1BinsCombinedv2-6.png]][[File:1stHitsLayer1BinsWeightedCombinedv2-6.png]]
 +
 +
 +
 +
[[File:ThetaBinsCombinedv2-6.png]][[File:ThetaBinsWeightedCombinedv2-6.png]]
 +
 +
 +
 +
 +
[[File:1stHitsThetaBinsCombinedv2-6.png]][[File:1stHitsThetaBinsWeightedCombinedv2-6.png]]
 +
 +
=Comparing GEMC txt outputs with LUND events=
 +
 +
===Event 1===
 +
LUND File Event #1
 +
<pre>
 +
2      1      1      1      1      0.006411431256  -81.074954045618540    -20140.877319449689821  1.625975760324849      -10735.677605378126827
 +
1      -1      1      11      0      0      -0.009994650877718      0.005770414374707      10.868088162529032      10.868094302146620      0.000511        0      0      0
 +
2      -1      1      11      0      0      0.009994650877718      -0.005770414374707      0.131911825606992      0.132416696768492      0.000511        0      0      0
 +
</pre>
 +
 +
GEMC 2.4 txt Output
 
<pre>
 
<pre>
−
*********************** EVENT # 7  ***********************
 
  
−
        +------------------------------------------------+------------+------------+
+
--- Header Bank --
−
        |                                            bank|      nrows|      ncols|
+
    - (10, 1) time:    Fri, 09.08.2017 13:28:00
−
        +------------------------------------------------+------------+------------+
+
    - (10, 5) beamPol: 1
−
        |                                        DC::dgtz|          27|          11|
+
    - (10, 3) evn:     1
−
        |                                        DC::true|          27|          24|
+
    - (10, 4) evn_type: -1
−
        |                                  GenPart::true|          2|          7|
+
    - (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 --
 +
</pre>
  
  
−
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)
 
  
 +
GEMC 2.6 txt output
 +
<pre>
 +
--- 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 --
 +
</pre>
  
−
Press Enter for Next Event or Bank Name: DC::true
+
===Event 4===
−
*****>>>>> BANK DC::true  >>>> SIZE = 24
+
<pre>
−
          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 
+
2       1      1       1       1       0.006411431256 -81.074954045618540     -20140.877319449689821 1.625975760324849      -10735.677605378126827
−
          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 
+
1      -1     1       11      0      0      -0.010054529398262      0.005665436152917       10.868088162529032      10.868094302146620      0.000511       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 
+
2      -1     1      11     0      0      0.010054529398262       -0.005665436152917      0.131911825606992       0.132416696768492       0.000511       0      0      0
−
          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>
 
</pre>
 +
 +
GEMC 2.4 txt output
 +
<pre>
 +
--- Header Bank --
 +
    - (10, 1) time:    Fri, 09.08.2017 13:28:02
 +
    - (10, 5) beamPol:  1
 +
    - (10, 3) evn:      4
 +
    - (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 --
 +
--- dc  (1300, 0) ----
 +
  -- integrated true infos bank  (1301, 0) --
 +
    - (1301, 1) pid:    11      11      11
 +
    - (1301, 2) mpid:  11      11      11
 +
    - (1301, 3) tid:    4      4      4
 +
    - (1301, 4) mtid:  2      2      2
 +
    - (1301, 5) otid:  2      2      2
 +
    - (1301, 6) trackE: 0.634108034368  0.633793371728  0.619215626462
 +
    - (1301, 7) totEdep:        0.000314662639331      0.0144607410525 0.0243402284788
 +
    - (1301, 8) avg_x:  131.791999321  123.566944653  94.9968565589
 +
    - (1301, 9) avg_y:  -470.15169746  -483.948184006  -507.869160412
 +
    - (1301, 10) avg_z: 3853.02700968  3860.9752654    3878.33257645
 +
    - (1301, 11) avg_lx:        -41.3093673896  -26.9736441028  10.6351028328
 +
    - (1301, 12) avg_ly:        -1140.84164507  -1138.58540668  -1132.5132615
 +
    - (1301, 13) avg_lz:        -9.30733652917  1.20768893133  -0.10133225752
 +
    - (1301, 14) px:    -0.106403185734 -0.115813395757 -0.248484313361
 +
    - (1301, 15) py:    -0.356534471221 -0.355047989359 -0.168066679021
 +
    - (1301, 16) pz:    0.0503452691474 0.0332014252875 0.179769973716
 +
    - (1301, 17) vx:    287.397369224  287.397369224  287.397369224
 +
    - (1301, 18) vy:    -166.271496667  -166.271496667  -166.271496667
 +
    - (1301, 19) vz:    3780.19056065  3780.19056065  3780.19056065
 +
    - (1301, 20) mvx:  0      0      0
 +
    - (1301, 21) mvy:  0      0      0
 +
    - (1301, 22) mvz:  0      0      0
 +
    - (1301, 23) avg_t: 14.5511687073  14.6575148415  14.92162824
 +
    - (1301, 24) nsteps:        1      38      52
 +
    - (1301, 25) procID:        1      1      1
 +
    - (1301, 99) hitn:  1      2      3
 +
  -- end of integrated raw bank.
 +
  -- step by step true infos bank  (1303, 0) --
 +
    - (1303, 1) pid:    11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11    $
 +
    - (1303, 2) mpid:  11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11    $
 +
    - (1303, 3) tid:    4      4      4      4      4      4      4      4      4      4      4      4      4      4      4      4      4      4      4      4      4    $
 +
    - (1303, 4) mtid:  2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2    $
 +
    - (1303, 5) otid:  2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2    $
 +
    - (1303, 6) trackE: 0.634108034368  0.633793371728  0.633477569588  0.633182010133  0.632899209368  0.632534038225  0.632363194029  0.631604098661  0.631089621082  0.630523245566  0.6304$
 +
    - (1303, 7) edep:  0.000314662639331      0.000315802140496      0.000295559454684      0.000282800765013      0.000365171142964      0.000170844195944      0.000759095368027    $
 +
    - (1303, 8) x:      131.791999321  131.525946939  131.062023807  130.574215386  130.136897089  129.581705142  129.125400977  128.705282768  128.425166483  128.216463957  127.95$
 +
    - (1303, 9) y:      -470.15169746  -471.096052992  -471.9383109    -472.739810395  -473.625688688  -474.42588539  -475.156216969  -475.884647232  -476.559136567  -477.246613989  -477.8$
 +
    - (1303, 10) z:    3853.02700968  3853.20644572  3853.47611913  3853.81856913  3853.95122507  3854.14990847  3854.65623004  3855.18651845  3855.85861592  3856.55186646  3857.3$
 +
    - (1303, 11) lx:    -41.3093673896  -40.5536816238  -39.7000779726  -38.853860213  -37.9905883128  -37.083929907  -36.3118208728  -35.5713242748  -34.9822898585  -34.444790485  -33.93$
 +
    - (1303, 12) ly:    -1140.84164507  -1140.37325444  -1140.12423758  -1139.94842081  -1139.59564698  -1139.39682119  -1139.32466331  -1139.24437984  -1139.18721707  -1139.09083666  -1139.$
 +
    - (1303, 13) lz:    -9.30733652917  -8.85529917147  -8.40065841982  -7.90002435215  -7.54797677835  -7.19235507837  -6.56259312245  -5.90415929155  -5.10736140828  -4.27154832868  -3.421$
 +
    - (1303, 14) px:    -0.106403185734 -0.115813395757 -0.136916080085 -0.150814762509 -0.145010754298 -0.199136002024 -0.13384355983  -0.154250037512 -0.115259477193 -0.110002800753 -0.116$
 +
    - (1303, 15) py:    -0.356534471221 -0.355047989359 -0.335251367992 -0.323395619811 -0.334703747574 -0.310022376893 -0.29547682199  -0.258428989113 -0.24757977418  -0.230527788047 -0.226$
 +
    - (1303, 16) pz:    0.0503452691474 0.0332014252875 0.0950497365759 0.111668045703  0.0799174930503 0.0566604095974 0.183146586919  0.217313678672  0.250147420835  0.266825785753  0.2672$
 +
    - (1303, 17) vx:    287.397369224  287.397369224  287.397369224  287.397369224  287.397369224  287.397369224  287.397369224  287.397369224  287.397369224  287.397369224  287.39$
 +
    - (1303, 18) vy:    -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.2$
 +
    - (1303, 19) vz:    3780.19056065  3780.19056065  3780.19056065  3780.19056065  3780.19056065  3780.19056065  3780.19056065  3780.19056065  3780.19056065  3780.19056065  3780.1$
 +
    - (1303, 20) mvx:  0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0    $
 +
  - (1303, 20) mvx:  0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0    $
 +
    - (1303, 21) mvy:  0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0    $
 +
    - (1303, 22) mvz:  0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0    $
 +
    - (1303, 3) t:      14.5511687073  14.5567644189  14.5623651239  14.5679705093  14.5735803796  14.5791960506  14.5848144389  14.5904450786  14.5960842096  14.6017327148  14.607$
 +
    - (1303, 98) stepn: 1      1      2      3      4      5      6      7      8      9      10      11      12      13      14      15      16      17      18      19      20    $
 +
    - (1303, 99) hitn:  1      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2    $
 +
  -- end of integrated raw bank.
 +
  -- integrated digitized bank  (1302, 0) --
 +
    - (1302, 1) sector: 6      6      6
 +
    - (1302, 2) superlayer:    4      4      4
 +
    - (1302, 3) layer:  1      1      2
 +
    - (1302, 4) wire:  6      7      7
 +
    - (1302, 5) LR:    1      -1      1
 +
    - (1302, 6) doca:  14.632819648    8.92145482021  2.82512271753
 +
    - (1302, 7) sdoca:  15.8954043687  9.16801612686  2.66655073379
 +
    - (1302, 8) time:  -562.800755694  343.1328777    108.658566059
 +
    - (1302, 9) stime:  -611.361706488  352.616004879  102.559643607
 +
    - (1302, 99) hitn:  1      2      3
 +
  -- End of integrated digitized bank.
 +
--- Generated Particles Bank --
 +
    - Particle 1 pid: 11  -  mom: (-10.0545293983,5.66543615292,10868.0881625) MeV  -  vert: (0,0,0) mm
 +
    - Hit >dc< Has 0 hit with nphe 52 nphe and time -1 ns
 +
    - Particle 2 pid: 11  -  mom: (10.0545293983,-5.66543615292,131.911825607) MeV  -  vert: (0,0,0) mm
 +
    - Hit >dc< Has 3 hits with etot 0.0391156321706 MeV and time 14.5511687073 ns
 +
--- End of Generated Particles Bank --
 +
---- End of Event  ----
 +
</pre>
 +
 +
 +
GEMC 2.6 txt output
 +
<pre>
 +
--- Header Bank --
 +
    - (10, 1) time:    Fri, 09.08.2017 14:09:27
 +
    - (10, 5) beamPol:  1
 +
    - (10, 3) evn:      4
 +
    - (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 --
 +
--- dc  (1300, 0) ----
 +
  -- integrated true infos bank  (1301, 0) --
 +
    - (1301, 1) pid:    11      11      11      11      11      11      11      11      11      11      11      11      11
 +
    - (1301, 2) mpid:  11      11      11      11      11      11      11      11      11      11      11      11      11
 +
    - (1301, 3) tid:    5      5      5      5      5      5      5      5      5      5      5      5      5
 +
    - (1301, 4) mtid:  2      2      2      2      2      2      2      2      2      2      2      2      2
 +
    - (1301, 5) otid:  2      2      2      2      2      2      2      2      2      2      2      2      2
 +
    - (1301, 6) trackE: 0.949388582004  0.940332214077  0.934494025085  0.927573842663  0.920814497501  0.914354487067  0.891081092992  0.84639023391  0.844990872798  0.838522744792  0.834$
 +
    - (1301, 7) totEdep:        0.00904606715855        0.00579501812076        0.00689659522427        0.00675934516173        0.00635791522896        0.0232733940759 0.00365813275209    $
 +
    - (1301, 8) avg_x:  759.379619692  765.484437268  774.275314803  776.456064884  784.315115793  790.550305531  804.181137618  883.421094035  890.879392559  897.152007278  899.8$
 +
    - (1301, 9) avg_y:  -298.328628898  -312.149860636  -330.381887207  -349.193864221  -368.751809558  -386.968027511  -397.80904012  -563.001316729  -573.311839997  -592.318661765  -605.$
 +
    - (1301, 10) avg_z: 4642.90906326  4665.45564462  4691.00331355  4717.34285842  4743.93926011  4757.88244598  4767.56722985  4802.90145032  4811.74475553  4832.5905734    4856.$
 +
    - (1301, 11) avg_lx:        489.743049701  502.25852125    523.004731632  541.022911591  563.222282761  579.716352502  589.634492075  -672.892638993  -671.656432199  -669.09811563$
 +
    - (1301, 12) avg_ly:        -1260.62082054  -1248.25460744  -1264.05828281  -1256.30151967  -1273.24729023  -1256.68812158  -1247.33941777  -1194.94271276  -1171.52565482  -1178.6382185$
 +
    - (1301, 13) avg_lz:        3.02815593863  0.242305589297  1.31154226572  0.304901701047  1.93080608734  -8.59727092904  5.94076266839  13.035201988    -0.600402606294 -0.4257531565$
 +
    - (1301, 14) px:    0.572402965856  0.168370376263  0.262418954714  -0.0347256988354        0.13373947209  0.246632288352  0.442502772887  0.253740192919  0.313486114187  0.12458741774$
 +
    - (1301, 15) py:    -0.500354741341 -0.362756831882 -0.429013599434 -0.50446479439  -0.405846971334 -0.557445528752 -0.422858975539 -0.495849460111 -0.458541616585 -0.483597315406 -0.31$
 +
    - (1301, 16) pz:    0.249436908086  0.680561300788  0.599368685039  0.586159316767  0.635752658076  0.450944829016  0.397853253172  0.380796230867  0.379941747065  0.438875972667  0.577$
 +
    - (1301, 17) vx:    333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.5$
 +
    - (1301, 18) vy:    -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.$
 +
    - (1301, 19) vz:    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.$
 +
    - (1301, 20) mvx:  0      0      0      0      0      0      0      0      0      0      0      0      0
 +
    - (1301, 21) mvy:  0      0      0      0      0      0      0      0      0      0      0      0      0
 +
    - (1301, 22) mvz:  0      0      0      0      0      0      0      0      0      0      0      0      0
 +
    - (1301, 23) avg_t: 16.588993761    16.6975009768  16.8281981999  16.9592766      17.099372919    17.2379652052  17.2731483319  18.0611350419  18.1264992527  18.2519112872  18.36$
 +
    - (1301, 24) nsteps:        32      32      33      32      35      26      21      2      29      28      28      28      30
 +
    - (1301, 25) procID:        1      1      1      1      1      1      1      1      1      1      1      1      1
 +
    - (1301, 99) hitn:  1      2      3      4      5      6      7      8      9      10      11      12      13
 +
  -- end of integrated raw bank.
 +
  -- step by step true infos bank  (1303, 0) --
 +
    - (1303, 1) pid:    11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11  $
 +
    - (1303, 2) mpid:  11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11  $
 +
    - (1303, 3) tid:    5      5      5      5      5      5      5      5      5      5      5      5      5      5      5      5      5      5      5      5      5    $
 +
    - (1303, 4) mtid:  2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2    $
 +
    - (1303, 5) otid:  2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2      2    $
 +
    - (1303, 6) trackE: 0.949388582004  0.949210228712  0.949059210299  0.948963656871  0.948819742953  0.948384748492  0.948141850466  0.947981442509  0.947847140282  0.947812874904  0.947$
 +
    - (1303, 7) edep:  0.000178353291944      0.000151018413463      9.55534272704e-05      0.000143913918743      0.000434994460428      0.000242898026525      0.000160407956567    $
 +
    - (1303, 8) x:      751.161065374  751.860717041  752.559370943  753.275632335  753.999182044  754.759518445  755.134418785  755.55021289    755.915476237  756.277413233  756.6$
 +
    - (1303, 9) y:      -289.835644655  -290.440224957  -291.019505524  -291.579579811  -292.138373284  -292.676598067  -293.090925158  -293.509671237  -293.927687466  -294.329965222  -294.$
 +
    - (1303, 10) z:    4629.94072447  4630.31963845  4630.73866614  4631.15306938  4631.55670308  4631.91908604  4632.74829221  4633.55510899  4634.3866782    4635.22749687  4636.$
 +
    - (1303, 11) lx:    481.502290699  482.05401634    482.582447413  483.089905087  483.594922755  484.074177145  484.487349052  484.900052603  485.317910787  485.720840437  486.1$
 +
    - (1303, 12) ly:    -1263.46568185  -1262.93111863  -1262.4169592  -1261.88699333  -1261.34606558  -1260.75679261  -1260.7240878  -1260.64465144  -1260.62124044  -1260.60636032  -1260$
 +
    - (1303, 13) lz:    -12.1984616258  -11.5593633653  -10.8843314137  -10.2060493848  -9.53444770052  -8.88468515353  -7.97472941215  -7.06778290217  -6.15975828958  -5.24475659136  -4.33$
 +
    - (1303, 14) px:    0.572402965856  0.571141727226  0.566961862957  0.562680310117  0.564586470154  0.594940037941  0.314602977748  0.308381710342  0.275065800282  0.291970211906  0.311$
 +
    - (1303, 15) py:    -0.500354741341 -0.48025925577  -0.471530546853 -0.4528699342  -0.456998840478 -0.452200886559 -0.354012683085 -0.346677580854 -0.360611472749 -0.350521754092 -0.33$
 +
    - (1303, 16) pz:    0.249436908086  0.288146398449  0.309526554546  0.343091572828  0.333965860704  0.282620690327  0.643080919628  0.649818670846  0.656956989883  0.655070419941  0.656$
 +
    - (1303, 17) vx:    333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.518384964  333.5$
 +
    - (1303, 18) vy:    -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.$
 +
    - (1303, 19) vz:    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.$
 +
    - (1303, 20) mvx:  0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0    $
 +
    - (1303, 21) mvy:  0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0    $
 +
    - (1303, 22) mvz:  0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0      0    $
 +
    - (1303, 3) t:      16.5164322252  16.5203865135  16.5243410622  16.5282957759  16.5322507379  16.5362064507  16.5401625829  16.544118992    16.5480756331  16.5520323334  16.55$
 +
    - (1303, 98) stepn: 1      2      3      4      5      6      7      8      9      10      11      12      13      14      15      16      17      18      19      20      21  $
 +
    - (1303, 99) hitn:  1      1      1      1      1      1      1      1      1      1      1      1      1      1      1      1      1      1      1      1      1    $
 +
  -- end of integrated raw bank.
 +
  -- integrated digitized bank  (1302, 0) --
 +
    - (1302, 1) sector: 1      1      1      1      1      1      1      6      6      6      6      6      6
 +
    - (1302, 2) layer:  25      26      27      28      29      30      30      31      32      33      34      35      36
 +
    - (1302, 3) wire:  14      15      14      14      14      14      15      18      19      19      19      19      19
 +
    - (1302, 4) tdc:    150.170275628  422.332559783  64.7021667673  489.984305976  227.451422915  -nan    378.909951474  449.615130971  31.4386354597  276.078058361  208.135598613$
 +
    - (1302, 5) LR:    1      -1      1      1      -1      1      -1      1      -1      -1      1      -1      1
 +
    - (1302, 6) doca:  5.76363904831  12.4786140173  2.76879069637  9.33072695957  7.88209072488  12.1510073193  12.2891326513  14.0770483845  1.59297893001  9.50112017251  7.393$
 +
    - (1302, 7) sdoca:  5.99412035875  12.4746059601  2.63254007937  9.17250481486  7.78364663978  12.1366797612  12.3360086408  14.2065554099  1.56833379855  9.4943429161    7.418$
 +
    - (1302, 8) time:  150.128974537  425.630835655  64.0505236738  499.930699204  221.66413493    -nan    376.624594825  435.115380612  35.1365854396  278.640463385  206.227280385$
 +
    - (1302, 9) stime:  150.170275628  422.332559783  64.7021667673  489.984305976  227.451422915  -nan    378.909951474  449.615130971  31.4386354597  276.078058361  208.135598613$
 +
    - (1302, 99) hitn:  1      2      3      4      5      6      7      8      9      10      11      12      13
 +
  -- End of integrated digitized bank.
 +
--- Generated Particles Bank --
 +
    - Particle 1 pid: 11  -  mom: (-10.0545293983,5.66543615292,10868.0881625) MeV  -  vert: (0,0,0) mm
 +
    - Hit >dc< Has 0 hit with nphe 30 nphe and time -1 ns
 +
    - Particle 2 pid: 11  -  mom: (10.0545293983,-5.66543615292,131.911825607) MeV  -  vert: (0,0,0) mm
 +
    - Hit >dc< Has 13 hits with etot 0.0944646746061 MeV and time 16.5164322252 ns
 +
--- Generated User Info Particles Bank --
 +
      1 -1 1 11 0 0 -0.0100545293983 0.00566543615292 10.8680881625 10.8680943021 0.000511 0 0 0
 +
      2 -1 1 11 0 0 0.0100545293983 -0.00566543615292 0.131911825607 0.132416696768 0.000511 0 0 0
 +
--- End of Generated Particles Bank --
 +
---- End of Event  ----
 +
</pre>
 +
 +
=evio2rootAnalysis=
 +
 +
The DC.xml files have changed
 +
 +
coat-java version/etc/bankdefs/DC.xml
 +
 +
coatjava 2.4
 +
<pre>
 +
  <section name="dgtz" tag="1302" num="0" info="uploading the hit definition">
 +
            <column name="sector"    num="1" type="int32" info="sector number"/>
 +
            <column name="superlayer" num="2" type="int32" info="layer number"/>
 +
            <column name="layer"      num="3" type="int32" info="view"/>
 +
            <column name="wire"      num="4" type="int32" info="strip number"/>
 +
</pre>
 +
 +
 +
coatjava 4a.0.0
 +
<pre>
 +
<section name="dgtz" tag="1302" num="0" info="uploading the hit definition">
 +
            <column name="sector"    num="1" type="int32" info="sector number"/>
 +
            <column name="layer"      num="2" type="int32" info="layer number (1-36)"/>
 +
            <column name="wire"      num="3" type="int32" info="wire number"/>
 +
</pre>
 +
The arrangement of how the wire number and layers are read has changed.  Eliminating the superlayer from GEMC2.6
 +
 +
GEMC2.4
 +
<pre>
 +
EVENT # 3  GENPART = 2Detector
 +
          Gen part 0 pid  = 11 px  = -10.0545 py  = 5.66544 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
 +
          Gen part 1 pid  = 11 px  = 10.0545 py  = -5.66544 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
 +
          Detector 0 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 4 layer  = 1 wire  = 6
 +
          Detector 1 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 4 layer  = 1 wire  = 7
 +
          Detector 2 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 4 layer  = 2 wire  = 7
 +
 +
EVENT # 23  GENPART = 2Detector
 +
          Gen part 0 pid  = 11 px  = -10.4252 py  = 4.95027 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
 +
          Gen part 1 pid  = 11 px  = 10.4252 py  = -4.95027 pz  = 131.912 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  = 1 
 +
          Detector 2 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 3 wire  = 1 
 +
          Detector 3 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 4 wire  = 1
 +
 +
EVENT # 26  GENPART = 2Detector
 +
          Gen part 0 pid  = 11 px  = -10.4765 py  = 4.84082 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
 +
          Gen part 1 pid  = 11 px  = 10.4765 py  = -4.84082 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
 +
</pre>
 +
 +
GEMC2.6
 +
<pre>
 +
EVENT # 3  GENPART = 2Detector
 +
          Gen part 0 pid  = 11 px  = -10.0545 py  = 5.66544 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
 +
          Gen part 1 pid  = 11 px  = 10.0545 py  = -5.66544 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
 +
          Detector 0 dpid  = 11 procID  = 1 sector  = 1 layer  = 25 wire  = 14
 +
          Detector 1 dpid  = 11 procID  = 1 sector  = 1 layer  = 26 wire  = 15
 +
          Detector 2 dpid  = 11 procID  = 1 sector  = 1 layer  = 27 wire  = 14
 +
          Detector 3 dpid  = 11 procID  = 1 sector  = 1 layer  = 28 wire  = 14
 +
          Detector 4 dpid  = 11 procID  = 1 sector  = 1 layer  = 29 wire  = 14
 +
          Detector 5 dpid  = 11 procID  = 1 sector  = 1 layer  = 30 wire  = 14
 +
          Detector 6 dpid  = 11 procID  = 1 sector  = 1 layer  = 30 wire  = 15
 +
          Detector 7 dpid  = 11 procID  = 1 sector  = 6 layer  = 31 wire  = 18
 +
          Detector 8 dpid  = 11 procID  = 1 sector  = 6 layer  = 32 wire  = 19
 +
          Detector 9 dpid  = 11 procID  = 1 sector  = 6 layer  = 33 wire  = 19
 +
          Detector 10 dpid  = 11 procID  = 1 sector  = 6 layer  = 34 wire  = 19
 +
          Detector 11 dpid  = 11 procID  = 1 sector  = 6 layer  = 35 wire  = 19
 +
          Detector 12 dpid  = 11 procID  = 1 sector  = 6 layer  = 36 wire  = 19
 +
 +
EVENT # 23  GENPART = 2Detector
 +
          Gen part 0 pid  = 11 px  = -10.4252 py  = 4.95027 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
 +
          Gen part 1 pid  = 11 px  = 10.4252 py  = -4.95027 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
 +
 +
EVENT # 26  GENPART = 2Detector
 +
          Gen part 0 pid  = 11 px  = -10.4765 py  = 4.84082 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
 +
          Gen part 1 pid  = 11 px  = 10.4765 py  = -4.84082 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
 +
          Detector 0 dpid  = 11 procID  = 999 sector  = 1 layer  = 25 wire  = 1
 +
          Detector 1 dpid  = 11 procID  = 999 sector  = 1 layer  = 26 wire  = 1
 +
          Detector 2 dpid  = 11 procID  = 999 sector  = 1 layer  = 28 wire  = 1
 +
</pre>
 +
 +
 +
GEMC 2.4 gcard[[File:eg12GEMC2_4gcard.txt]]
 +
 +
 +
GEMC 2.6 gcard[[File:eg12GEMC2_6gcard.txt]]
 +
 +
 +
CED shows
 +
 +
 +
GEMC 2.4
 +
[[File:GEMC2_4CED.png]]
 +
 +
 +
GEMC 2.6
 +
[[File:GEMC2_6CED.png]]
 +
 +
 +
The GEANT simulation results look to be identical
 +
 +
GEMC 2.6
 +
[[File:GEANT_GEMC2_6.png]]
 +
 +
 +
GEMC 2.4
 +
[[File:GEANT_GEMC2_4.png]]

Latest revision as of 02:50, 14 October 2017

Using the evio2root converter

File:Evio2root.c 

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

GEMC 2.4

 > Common Environment Version: <1.3>  (Thu March 31 2016)
 > Running as vanwdani on minerve2
 > OS Release:    Linux__Ubuntu14.04-x86_64-gcc4.8
 > JLAB_ROOT set to:     /home/vanwdani/src/CLAS/GEMC
 > JLAB_SOFTWARE set to: /home/vanwdani/src/CLAS/GEMC/1.3/Linux__Ubuntu14.04-x86_64-gcc4.8

 > CLHEP   version: 2.2.0.8
 > XERCES  version: 3.1.1
 > QTDIR   version: 5.3.2
 > Geant4  version: 4.10.01.p03
 > ROOT    version: 5.34.34
 > GEMC    version: 2.3
 > JANA    version: 0.7.4p2
 > Build   version: 1.3
 > EVIO    version: 4.4.4
 > CCDB    version: 1.06
 > Banks   version: 1.2


Base Results

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


Editing the eg12.gcard so that there is no energy cut and that energy hits of 0 are recorded only results in 103883 hits

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


Similarly, editing the gcard only results in 92761 1st hits via wires

These wire hits can be expressed in terms of Theta bins

ThetaBins.pngThetaBinsWeighted.png


With the edited gcard we only have 68114 hits per Theta bins


Similarly, these can account for only the 1st interactions

1stHitsThetaBins.png1stHitsThetaBinsWeighted.png


With the edited gcard, we find 61102 hits


Overlaying this with the expected differential cross-section expected from the LUND file

EfficiencyThetaBins.png

We can find the efficiency of the detector with respect to the particles input via the LUND file 

TH1F *Ratio=new TH1F("TheoryExperiment","Theoretical and Experimental Differential Cross-Section Lab Frame",90,0,45);
Ratio->Add(L1stHitsThetaBinsWeighted);
Ratio->Divide(MolThetaLabweighted);
Ratio->Draw();
Ratio.png


This can be partially explained in that Superlayer 1, Layer 1 is not always hit first

LayersHitFirst.png

Lower Fix

Determining what values of phi are not within the sector following the Mathematica equation developed earlier

RejectedTheta.pngRejectedThetaWeighted.png


CorrectedMolThetaLab.pngCorrectedMolThetaLabWeighted.png


EfficiencyThetaBinsCorrected.png


RatioCorrected.png

Upper Fix

Fixing the upper tail of hits on wires from larger theta values


Layer1WireNumberBinTotalAdditional.png Layer1WireNumberBinTotalWeightedAdditional.png



Layer1binsHitsAdditional.png Layer1binsHitsWeightedAdditional.png



1stHitsLayer1binsHitsAdditional.png 1stHitsLayer1binsHitsWeightedAdditional.png

Combined Results

Combining these plots 

hadd combined.root NonIsotropic_BinAdditions.root LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.root

Layer1WireNumberTotalCombined.pngLayer1WireNumberBinTotalWeightedCombined.png


Layer1binsHitsCombined.pngLayer1binsHitsWeightedCombined.png


1stHitsLayer1binsHitsCombined.png1stHitsLayer1binsHitsWeightedCombined.png

GEMC 2.6

> Common Environment Version: <2.1>  (Wed January 25 2017)
 > Running as lds on lds3
 > OS Release:    Linux__Ubuntu16.04-x86_64-gcc5.4.0
 > JLAB_ROOT set to:     /home/lds/src/CLAS/GEMC
 > JLAB_SOFTWARE set to: /home/lds/src/CLAS/GEMC/2.1/Linux__Ubuntu16.04-x86_64-gcc5.4.0

 > BANKS 	 version:  1.3
 > CCDB 	 version:  1.06.02
 > CLHEP 	 version:  2.3.4.3
 > EVIO 	 version:  5.1
 > GEANT4 	 version:  4.10.02.p03
 > GEMC 	 version:  devel
 > JANA 	 version:  0.7.7p1
 > MLIBRARY 	 version:  1.1
 > MYSQL 	 installed in  /home/lds/src/CLAS/GEMC/2.1/Linux__Ubuntu16.04-x86_64-gcc5.4.0/mysql/lib
 > QT   	 version:  5.8.0
 > ROOT 	 version:  6.10.02
 > SCONS 	 version:  1.5
 > XERCESC 	 version:  3.1.4

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

ProcessesGEMC3.png

Base Results

GEMC 3.0 seems to use a new LUND format that has an additional variable on the header (11 variables instead of 10). The databank for the DC has also eliminated the superlayer variable from the Digitialized output, moving up the wire number as a result.

Layer1BinHitsv2-6local.pngLayer1BinHitsWeightedGEMC3.png


Restricting results to only the 1st hits on this layer

1stHitsLayer1binHitsGEMC3.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

Upper fix

Layer1WireNumberBinTotaAdditionalv2-6.pngLayer1WireNumberBinTotalWeightedAdditionalv2-6.png


Layer1binsHitsAdditionalv2-6.pngLayer1binsHitsWeightedAdditionalv2-6.png



1stHitsLayer1binsHitsAdditionalv2-6.png1stHitsLayer1binsHitsWeightedAdditionalv2-6.png


Combined Results

Layer1WireNumberTotalCombinedv2-6.pngLayer1WireNumberTotalWeightedCombinedv2-6.png


Layer1BinHitsCombinedv2-6.pngLayer1BinHitsWeightedCombinedv2-6.png



1stHitsLayer1BinsCombinedv2-6.png1stHitsLayer1BinsWeightedCombinedv2-6.png


ThetaBinsCombinedv2-6.pngThetaBinsWeightedCombinedv2-6.png



1stHitsThetaBinsCombinedv2-6.png1stHitsThetaBinsWeightedCombinedv2-6.png

Comparing GEMC txt outputs with LUND events

Event 1

LUND File Event #1 

2       1       1       1       1       0.006411431256  -81.074954045618540     -20140.877319449689821  1.625975760324849       -10735.677605378126827
1       -1      1       11      0       0       -0.009994650877718      0.005770414374707       10.868088162529032      10.868094302146620      0.000511        0       0       0
2       -1      1       11      0       0       0.009994650877718       -0.005770414374707      0.131911825606992       0.132416696768492       0.000511        0       0       0

GEMC 2.4 txt Output 


 --- 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 --


GEMC 2.6 txt output 

 --- 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 --

Event 4

2       1       1       1       1       0.006411431256  -81.074954045618540     -20140.877319449689821  1.625975760324849       -10735.677605378126827
1       -1      1       11      0       0       -0.010054529398262      0.005665436152917       10.868088162529032      10.868094302146620      0.000511        0       0       0
2       -1      1       11      0       0       0.010054529398262       -0.005665436152917      0.131911825606992       0.132416696768492       0.000511        0       0       0

GEMC 2.4 txt output 

 --- Header Bank --
    - (10, 1) time:     Fri, 09.08.2017 13:28:02
    - (10, 5) beamPol:  1
    - (10, 3) evn:      4
    - (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 --
 --- dc  (1300, 0) ----
   -- integrated true infos bank  (1301, 0) --
    - (1301, 1) pid:    11      11      11
    - (1301, 2) mpid:   11      11      11
    - (1301, 3) tid:    4       4       4
    - (1301, 4) mtid:   2       2       2
    - (1301, 5) otid:   2       2       2
    - (1301, 6) trackE: 0.634108034368  0.633793371728  0.619215626462
    - (1301, 7) totEdep:        0.000314662639331       0.0144607410525 0.0243402284788
    - (1301, 8) avg_x:  131.791999321   123.566944653   94.9968565589
    - (1301, 9) avg_y:  -470.15169746   -483.948184006  -507.869160412
    - (1301, 10) avg_z: 3853.02700968   3860.9752654    3878.33257645
    - (1301, 11) avg_lx:        -41.3093673896  -26.9736441028  10.6351028328
    - (1301, 12) avg_ly:        -1140.84164507  -1138.58540668  -1132.5132615
    - (1301, 13) avg_lz:        -9.30733652917  1.20768893133   -0.10133225752
    - (1301, 14) px:    -0.106403185734 -0.115813395757 -0.248484313361
    - (1301, 15) py:    -0.356534471221 -0.355047989359 -0.168066679021
    - (1301, 16) pz:    0.0503452691474 0.0332014252875 0.179769973716
    - (1301, 17) vx:    287.397369224   287.397369224   287.397369224
    - (1301, 18) vy:    -166.271496667  -166.271496667  -166.271496667
    - (1301, 19) vz:    3780.19056065   3780.19056065   3780.19056065
    - (1301, 20) mvx:   0       0       0
    - (1301, 21) mvy:   0       0       0
    - (1301, 22) mvz:   0       0       0
    - (1301, 23) avg_t: 14.5511687073   14.6575148415   14.92162824
    - (1301, 24) nsteps:        1       38      52
    - (1301, 25) procID:        1       1       1
    - (1301, 99) hitn:  1       2       3
   -- end of integrated raw bank.
   -- step by step true infos bank  (1303, 0) --
    - (1303, 1) pid:    11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11    $
    - (1303, 2) mpid:   11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11    $
    - (1303, 3) tid:    4       4       4       4       4       4       4       4       4       4       4       4       4       4       4       4       4       4       4       4       4     $
    - (1303, 4) mtid:   2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2     $
    - (1303, 5) otid:   2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2     $
    - (1303, 6) trackE: 0.634108034368  0.633793371728  0.633477569588  0.633182010133  0.632899209368  0.632534038225  0.632363194029  0.631604098661  0.631089621082  0.630523245566  0.6304$
    - (1303, 7) edep:   0.000314662639331       0.000315802140496       0.000295559454684       0.000282800765013       0.000365171142964       0.000170844195944       0.000759095368027     $
    - (1303, 8) x:      131.791999321   131.525946939   131.062023807   130.574215386   130.136897089   129.581705142   129.125400977   128.705282768   128.425166483   128.216463957   127.95$
    - (1303, 9) y:      -470.15169746   -471.096052992  -471.9383109    -472.739810395  -473.625688688  -474.42588539   -475.156216969  -475.884647232  -476.559136567  -477.246613989  -477.8$
    - (1303, 10) z:     3853.02700968   3853.20644572   3853.47611913   3853.81856913   3853.95122507   3854.14990847   3854.65623004   3855.18651845   3855.85861592   3856.55186646   3857.3$
    - (1303, 11) lx:    -41.3093673896  -40.5536816238  -39.7000779726  -38.853860213   -37.9905883128  -37.083929907   -36.3118208728  -35.5713242748  -34.9822898585  -34.444790485   -33.93$
    - (1303, 12) ly:    -1140.84164507  -1140.37325444  -1140.12423758  -1139.94842081  -1139.59564698  -1139.39682119  -1139.32466331  -1139.24437984  -1139.18721707  -1139.09083666  -1139.$
    - (1303, 13) lz:    -9.30733652917  -8.85529917147  -8.40065841982  -7.90002435215  -7.54797677835  -7.19235507837  -6.56259312245  -5.90415929155  -5.10736140828  -4.27154832868  -3.421$
    - (1303, 14) px:    -0.106403185734 -0.115813395757 -0.136916080085 -0.150814762509 -0.145010754298 -0.199136002024 -0.13384355983  -0.154250037512 -0.115259477193 -0.110002800753 -0.116$
    - (1303, 15) py:    -0.356534471221 -0.355047989359 -0.335251367992 -0.323395619811 -0.334703747574 -0.310022376893 -0.29547682199  -0.258428989113 -0.24757977418  -0.230527788047 -0.226$
    - (1303, 16) pz:    0.0503452691474 0.0332014252875 0.0950497365759 0.111668045703  0.0799174930503 0.0566604095974 0.183146586919  0.217313678672  0.250147420835  0.266825785753  0.2672$
    - (1303, 17) vx:    287.397369224   287.397369224   287.397369224   287.397369224   287.397369224   287.397369224   287.397369224   287.397369224   287.397369224   287.397369224   287.39$
    - (1303, 18) vy:    -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.271496667  -166.2$
    - (1303, 19) vz:    3780.19056065   3780.19056065   3780.19056065   3780.19056065   3780.19056065   3780.19056065   3780.19056065   3780.19056065   3780.19056065   3780.19056065   3780.1$
    - (1303, 20) mvx:   0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0     $
   - (1303, 20) mvx:   0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0     $
    - (1303, 21) mvy:   0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0     $
    - (1303, 22) mvz:   0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0     $
    - (1303, 3) t:      14.5511687073   14.5567644189   14.5623651239   14.5679705093   14.5735803796   14.5791960506   14.5848144389   14.5904450786   14.5960842096   14.6017327148   14.607$
    - (1303, 98) stepn: 1       1       2       3       4       5       6       7       8       9       10      11      12      13      14      15      16      17      18      19      20    $
    - (1303, 99) hitn:  1       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2     $
   -- end of integrated raw bank.
   -- integrated digitized bank  (1302, 0) --
    - (1302, 1) sector: 6       6       6
    - (1302, 2) superlayer:     4       4       4
    - (1302, 3) layer:  1       1       2
    - (1302, 4) wire:   6       7       7
    - (1302, 5) LR:     1       -1      1
    - (1302, 6) doca:   14.632819648    8.92145482021   2.82512271753
    - (1302, 7) sdoca:  15.8954043687   9.16801612686   2.66655073379
    - (1302, 8) time:   -562.800755694  343.1328777     108.658566059
    - (1302, 9) stime:  -611.361706488  352.616004879   102.559643607
    - (1302, 99) hitn:  1       2       3
   -- End of integrated digitized bank.
 --- Generated Particles Bank --
    - Particle 1 pid: 11   -  mom: (-10.0545293983,5.66543615292,10868.0881625) MeV   -  vert: (0,0,0) mm
    - Hit >dc< Has 0 hit with nphe 52 nphe and time -1 ns
    - Particle 2 pid: 11   -  mom: (10.0545293983,-5.66543615292,131.911825607) MeV   -  vert: (0,0,0) mm
    - Hit >dc< Has 3 hits with etot 0.0391156321706 MeV and time 14.5511687073 ns
 --- End of Generated Particles Bank --
 ---- End of Event  ----


GEMC 2.6 txt output 

 --- Header Bank --
    - (10, 1) time:     Fri, 09.08.2017 14:09:27
    - (10, 5) beamPol:  1
    - (10, 3) evn:      4
    - (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 --
 --- dc  (1300, 0) ----
   -- integrated true infos bank  (1301, 0) --
    - (1301, 1) pid:    11      11      11      11      11      11      11      11      11      11      11      11      11
    - (1301, 2) mpid:   11      11      11      11      11      11      11      11      11      11      11      11      11
    - (1301, 3) tid:    5       5       5       5       5       5       5       5       5       5       5       5       5
    - (1301, 4) mtid:   2       2       2       2       2       2       2       2       2       2       2       2       2
    - (1301, 5) otid:   2       2       2       2       2       2       2       2       2       2       2       2       2
    - (1301, 6) trackE: 0.949388582004  0.940332214077  0.934494025085  0.927573842663  0.920814497501  0.914354487067  0.891081092992  0.84639023391   0.844990872798  0.838522744792  0.834$
    - (1301, 7) totEdep:        0.00904606715855        0.00579501812076        0.00689659522427        0.00675934516173        0.00635791522896        0.0232733940759 0.00365813275209     $
    - (1301, 8) avg_x:  759.379619692   765.484437268   774.275314803   776.456064884   784.315115793   790.550305531   804.181137618   883.421094035   890.879392559   897.152007278   899.8$
    - (1301, 9) avg_y:  -298.328628898  -312.149860636  -330.381887207  -349.193864221  -368.751809558  -386.968027511  -397.80904012   -563.001316729  -573.311839997  -592.318661765  -605.$
    - (1301, 10) avg_z: 4642.90906326   4665.45564462   4691.00331355   4717.34285842   4743.93926011   4757.88244598   4767.56722985   4802.90145032   4811.74475553   4832.5905734    4856.$
    - (1301, 11) avg_lx:        489.743049701   502.25852125    523.004731632   541.022911591   563.222282761   579.716352502   589.634492075   -672.892638993  -671.656432199  -669.09811563$
    - (1301, 12) avg_ly:        -1260.62082054  -1248.25460744  -1264.05828281  -1256.30151967  -1273.24729023  -1256.68812158  -1247.33941777  -1194.94271276  -1171.52565482  -1178.6382185$
    - (1301, 13) avg_lz:        3.02815593863   0.242305589297  1.31154226572   0.304901701047  1.93080608734   -8.59727092904  5.94076266839   13.035201988    -0.600402606294 -0.4257531565$
    - (1301, 14) px:    0.572402965856  0.168370376263  0.262418954714  -0.0347256988354        0.13373947209   0.246632288352  0.442502772887  0.253740192919  0.313486114187  0.12458741774$
    - (1301, 15) py:    -0.500354741341 -0.362756831882 -0.429013599434 -0.50446479439  -0.405846971334 -0.557445528752 -0.422858975539 -0.495849460111 -0.458541616585 -0.483597315406 -0.31$
    - (1301, 16) pz:    0.249436908086  0.680561300788  0.599368685039  0.586159316767  0.635752658076  0.450944829016  0.397853253172  0.380796230867  0.379941747065  0.438875972667  0.577$
    - (1301, 17) vx:    333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.5$
    - (1301, 18) vy:    -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.$
    - (1301, 19) vz:    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.$
    - (1301, 20) mvx:   0       0       0       0       0       0       0       0       0       0       0       0       0
    - (1301, 21) mvy:   0       0       0       0       0       0       0       0       0       0       0       0       0
    - (1301, 22) mvz:   0       0       0       0       0       0       0       0       0       0       0       0       0
    - (1301, 23) avg_t: 16.588993761    16.6975009768   16.8281981999   16.9592766      17.099372919    17.2379652052   17.2731483319   18.0611350419   18.1264992527   18.2519112872   18.36$
    - (1301, 24) nsteps:        32      32      33      32      35      26      21      2       29      28      28      28      30
    - (1301, 25) procID:        1       1       1       1       1       1       1       1       1       1       1       1       1
    - (1301, 99) hitn:  1       2       3       4       5       6       7       8       9       10      11      12      13
   -- end of integrated raw bank.
   -- step by step true infos bank  (1303, 0) --
    - (1303, 1) pid:    11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11   $
    - (1303, 2) mpid:   11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11      11   $
    - (1303, 3) tid:    5       5       5       5       5       5       5       5       5       5       5       5       5       5       5       5       5       5       5       5       5    $
    - (1303, 4) mtid:   2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2    $
    - (1303, 5) otid:   2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2       2    $
    - (1303, 6) trackE: 0.949388582004  0.949210228712  0.949059210299  0.948963656871  0.948819742953  0.948384748492  0.948141850466  0.947981442509  0.947847140282  0.947812874904  0.947$
    - (1303, 7) edep:   0.000178353291944       0.000151018413463       9.55534272704e-05       0.000143913918743       0.000434994460428       0.000242898026525       0.000160407956567    $
    - (1303, 8) x:      751.161065374   751.860717041   752.559370943   753.275632335   753.999182044   754.759518445   755.134418785   755.55021289    755.915476237   756.277413233   756.6$
    - (1303, 9) y:      -289.835644655  -290.440224957  -291.019505524  -291.579579811  -292.138373284  -292.676598067  -293.090925158  -293.509671237  -293.927687466  -294.329965222  -294.$
    - (1303, 10) z:     4629.94072447   4630.31963845   4630.73866614   4631.15306938   4631.55670308   4631.91908604   4632.74829221   4633.55510899   4634.3866782    4635.22749687   4636.$
    - (1303, 11) lx:    481.502290699   482.05401634    482.582447413   483.089905087   483.594922755   484.074177145   484.487349052   484.900052603   485.317910787   485.720840437   486.1$
    - (1303, 12) ly:    -1263.46568185  -1262.93111863  -1262.4169592   -1261.88699333  -1261.34606558  -1260.75679261  -1260.7240878   -1260.64465144  -1260.62124044  -1260.60636032  -1260$
    - (1303, 13) lz:    -12.1984616258  -11.5593633653  -10.8843314137  -10.2060493848  -9.53444770052  -8.88468515353  -7.97472941215  -7.06778290217  -6.15975828958  -5.24475659136  -4.33$
    - (1303, 14) px:    0.572402965856  0.571141727226  0.566961862957  0.562680310117  0.564586470154  0.594940037941  0.314602977748  0.308381710342  0.275065800282  0.291970211906  0.311$
    - (1303, 15) py:    -0.500354741341 -0.48025925577  -0.471530546853 -0.4528699342   -0.456998840478 -0.452200886559 -0.354012683085 -0.346677580854 -0.360611472749 -0.350521754092 -0.33$
    - (1303, 16) pz:    0.249436908086  0.288146398449  0.309526554546  0.343091572828  0.333965860704  0.282620690327  0.643080919628  0.649818670846  0.656956989883  0.655070419941  0.656$
    - (1303, 17) vx:    333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.518384964   333.5$
    - (1303, 18) vy:    -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.080756828  -181.$
    - (1303, 19) vz:    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.9164769    4330.$
    - (1303, 20) mvx:   0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0    $
    - (1303, 21) mvy:   0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0    $
    - (1303, 22) mvz:   0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0       0    $
    - (1303, 3) t:      16.5164322252   16.5203865135   16.5243410622   16.5282957759   16.5322507379   16.5362064507   16.5401625829   16.544118992    16.5480756331   16.5520323334   16.55$
    - (1303, 98) stepn: 1       2       3       4       5       6       7       8       9       10      11      12      13      14      15      16      17      18      19      20      21   $
    - (1303, 99) hitn:  1       1       1       1       1       1       1       1       1       1       1       1       1       1       1       1       1       1       1       1       1    $
   -- end of integrated raw bank.
   -- integrated digitized bank  (1302, 0) --
    - (1302, 1) sector: 1       1       1       1       1       1       1       6       6       6       6       6       6
    - (1302, 2) layer:  25      26      27      28      29      30      30      31      32      33      34      35      36
    - (1302, 3) wire:   14      15      14      14      14      14      15      18      19      19      19      19      19
    - (1302, 4) tdc:    150.170275628   422.332559783   64.7021667673   489.984305976   227.451422915   -nan    378.909951474   449.615130971   31.4386354597   276.078058361   208.135598613$
    - (1302, 5) LR:     1       -1      1       1       -1      1       -1      1       -1      -1      1       -1      1
    - (1302, 6) doca:   5.76363904831   12.4786140173   2.76879069637   9.33072695957   7.88209072488   12.1510073193   12.2891326513   14.0770483845   1.59297893001   9.50112017251   7.393$
    - (1302, 7) sdoca:  5.99412035875   12.4746059601   2.63254007937   9.17250481486   7.78364663978   12.1366797612   12.3360086408   14.2065554099   1.56833379855   9.4943429161    7.418$
    - (1302, 8) time:   150.128974537   425.630835655   64.0505236738   499.930699204   221.66413493    -nan    376.624594825   435.115380612   35.1365854396   278.640463385   206.227280385$
    - (1302, 9) stime:  150.170275628   422.332559783   64.7021667673   489.984305976   227.451422915   -nan    378.909951474   449.615130971   31.4386354597   276.078058361   208.135598613$
    - (1302, 99) hitn:  1       2       3       4       5       6       7       8       9       10      11      12      13
   -- End of integrated digitized bank.
 --- Generated Particles Bank --
    - Particle 1 pid: 11   -  mom: (-10.0545293983,5.66543615292,10868.0881625) MeV   -  vert: (0,0,0) mm
    - Hit >dc< Has 0 hit with nphe 30 nphe and time -1 ns
    - Particle 2 pid: 11   -  mom: (10.0545293983,-5.66543615292,131.911825607) MeV   -  vert: (0,0,0) mm
    - Hit >dc< Has 13 hits with etot 0.0944646746061 MeV and time 16.5164322252 ns
 --- Generated User Info Particles Bank --
      1 -1 1 11 0 0 -0.0100545293983 0.00566543615292 10.8680881625 10.8680943021 0.000511 0 0 0
      2 -1 1 11 0 0 0.0100545293983 -0.00566543615292 0.131911825607 0.132416696768 0.000511 0 0 0
 --- End of Generated Particles Bank --
 ---- End of Event  ----

evio2rootAnalysis

The DC.xml files have changed

coat-java version/etc/bankdefs/DC.xml

coatjava 2.4 

  <section name="dgtz" tag="1302" num="0" info="uploading the hit definition">
            <column name="sector"     num="1" type="int32" info="sector number"/>
            <column name="superlayer" num="2" type="int32" info="layer number"/>
            <column name="layer"      num="3" type="int32" info="view"/>
            <column name="wire"       num="4" type="int32" info="strip number"/>


coatjava 4a.0.0 

 <section name="dgtz" tag="1302" num="0" info="uploading the hit definition">
            <column name="sector"     num="1" type="int32" info="sector number"/>
            <column name="layer"      num="2" type="int32" info="layer number (1-36)"/>
            <column name="wire"       num="3" type="int32" info="wire number"/>

The arrangement of how the wire number and layers are read has changed. Eliminating the superlayer from GEMC2.6

GEMC2.4 

 EVENT # 3  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = -10.0545 py  = 5.66544 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = 10.0545 py  = -5.66544 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
          Detector 0 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 4 layer  = 1 wire  = 6
          Detector 1 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 4 layer  = 1 wire  = 7
          Detector 2 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 4 layer  = 2 wire  = 7

 EVENT # 23  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = -10.4252 py  = 4.95027 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = 10.4252 py  = -4.95027 pz  = 131.912 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  = 1   
          Detector 2 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 3 wire  = 1   
          Detector 3 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 4 wire  = 1

EVENT # 26  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = -10.4765 py  = 4.84082 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = 10.4765 py  = -4.84082 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0

GEMC2.6 

EVENT # 3  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = -10.0545 py  = 5.66544 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = 10.0545 py  = -5.66544 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
          Detector 0 dpid  = 11 procID  = 1 sector  = 1 layer  = 25 wire  = 14
          Detector 1 dpid  = 11 procID  = 1 sector  = 1 layer  = 26 wire  = 15
          Detector 2 dpid  = 11 procID  = 1 sector  = 1 layer  = 27 wire  = 14
          Detector 3 dpid  = 11 procID  = 1 sector  = 1 layer  = 28 wire  = 14
          Detector 4 dpid  = 11 procID  = 1 sector  = 1 layer  = 29 wire  = 14
          Detector 5 dpid  = 11 procID  = 1 sector  = 1 layer  = 30 wire  = 14
          Detector 6 dpid  = 11 procID  = 1 sector  = 1 layer  = 30 wire  = 15
          Detector 7 dpid  = 11 procID  = 1 sector  = 6 layer  = 31 wire  = 18
          Detector 8 dpid  = 11 procID  = 1 sector  = 6 layer  = 32 wire  = 19
          Detector 9 dpid  = 11 procID  = 1 sector  = 6 layer  = 33 wire  = 19
          Detector 10 dpid  = 11 procID  = 1 sector  = 6 layer  = 34 wire  = 19
          Detector 11 dpid  = 11 procID  = 1 sector  = 6 layer  = 35 wire  = 19
          Detector 12 dpid  = 11 procID  = 1 sector  = 6 layer  = 36 wire  = 19

EVENT # 23  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = -10.4252 py  = 4.95027 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = 10.4252 py  = -4.95027 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0

EVENT # 26  GENPART = 2Detector
          Gen part 0 pid  = 11 px  = -10.4765 py  = 4.84082 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
          Gen part 1 pid  = 11 px  = 10.4765 py  = -4.84082 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
          Detector 0 dpid  = 11 procID  = 999 sector  = 1 layer  = 25 wire  = 1
          Detector 1 dpid  = 11 procID  = 999 sector  = 1 layer  = 26 wire  = 1
          Detector 2 dpid  = 11 procID  = 999 sector  = 1 layer  = 28 wire  = 1


GEMC 2.4 gcardFile:Eg12GEMC2 4gcard.txt


GEMC 2.6 gcardFile:Eg12GEMC2 6gcard.txt


CED shows


GEMC 2.4 GEMC2 4CED.png


GEMC 2.6 GEMC2 6CED.png


The GEANT simulation results look to be identical

GEMC 2.6 GEANT GEMC2 6.png


GEMC 2.4 GEANT GEMC2 4.png

Retrieved from "https://wiki.iac.isu.edu/index.php?title=Analysis_of_evio_files&oldid=118696"