Difference between revisions of "Analysis of evio files"

From New IAC Wiki
Jump to navigation Jump to search
 
(51 intermediate revisions by 2 users not shown)
Line 8: Line 8:
  
 
=GEMC 2.4=
 
=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
 
Restricting hits for Superlayer 1, Layer 1 only
  
Line 23: Line 44:
  
  
 +
 +
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.
 
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.
Line 29: Line 52:
  
  
 +
 +
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
 
Fixing the upper tail of hits on wires from larger theta values
  
  
−
[[File:Layer1WireNumberTotalAdditionals.png]] [[File:Layer1WireNumberBinTotalWeightedAdditional.png]]
+
[[File:Layer1WireNumberBinTotalAdditional.png]] [[File:Layer1WireNumberBinTotalWeightedAdditional.png]]
  
  
  
  
−
[[File:Layer1binsHitsAdditional.png]][[File:1stHitsLayer1binsHitsWeightedAdditional.png]]
+
[[File:Layer1binsHitsAdditional.png]] [[File:Layer1binsHitsWeightedAdditional.png]]
  
  
  
  
−
[[File:1stHitsLayer1binsHitsAdditional.png]][[File:1stHitsLayer1binsHitsWeightedAdditional.png]]
+
[[File:1stHitsLayer1binsHitsAdditional.png]] [[File:1stHitsLayer1binsHitsWeightedAdditional.png]]
 +
 
 +
==Combined Results==
  
 
Combining these plots
 
Combining these plots
Line 60: Line 142:
 
[[File:1stHitsLayer1binsHitsCombined.png]][[File:1stHitsLayer1binsHitsWeightedCombined.png]]
 
[[File:1stHitsLayer1binsHitsCombined.png]][[File:1stHitsLayer1binsHitsWeightedCombined.png]]
  
−
=GEMC 3.0=
+
=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
 
GEMC 3.0 has changed the particle ID for transportation from 90 to 999
  
 
[[File:ProcessesGEMC3.png]]
 
[[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.
  
−
GEMC 3.0 seems to use a new LUND format that has an additional variable on the header (11 variables instead of 10).
+
[[File:Layer1BinHitsv2-6local.png]][[File:Layer1BinHitsWeightedGEMC3.png]]
−
 
 
−
[[File:Layer1BinHitsGEMC3.png]][[File:Layer1BinHitsWeightedGEMC3.png]]
 
  
  
 
Restricting results to only the 1st hits on this layer
 
Restricting results to only the 1st hits on this layer
  
−
[[File:1stHitsLayer1binsHitsGEMC3.png]][[File:1stHitsLayer1binsHitsWeightedGEMC3.png]]
+
[[File:1stHitsLayer1binHitsGEMC3.png]][[File:1stHitsLayer1binsHitsWeightedGEMC3.png]]
  
 
This lack of interactions translates to angle bins as well
 
This lack of interactions translates to angle bins as well
Line 85: Line 190:
 
[[File:1stHitsThetaBinsGEMC3.png]][[File:1stHitsThetaBinsWeightedGEMC3.png]]
 
[[File:1stHitsThetaBinsGEMC3.png]][[File:1stHitsThetaBinsWeightedGEMC3.png]]
  
−
=evio2rootAnalysis=
+
==Upper fix==
 +
 
 +
[[File:Layer1WireNumberBinTotaAdditionalv2-6.png]][[File:Layer1WireNumberBinTotalWeightedAdditionalv2-6.png]]
 +
 
 +
 
 +
 
 +
[[File:Layer1binsHitsAdditionalv2-6.png]][[File:Layer1binsHitsWeightedAdditionalv2-6.png]]
 +
 
 +
 
 +
 
 +
 
 +
[[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>
−
GEMC2.4
+
 
 
  --- Header Bank --
 
  --- Header Bank --
 
     - (10, 1) time:    Fri, 09.08.2017 13:28:00
 
     - (10, 1) time:    Fri, 09.08.2017 13:28:00
Line 107: Line 258:
 
     - Particle 2 pid: 11  -  mom: (9.99465,-5.77041,131.912) 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 --
 
  --- End of Generated Particles Bank --
 +
</pre>
 +
  
−
GEMC2.6
+
 
 +
GEMC 2.6 txt output
 +
<pre>
 
  --- Header Bank --
 
  --- Header Bank --
 
     - (10, 1) time:    Fri, 09.08.2017 14:09:27
 
     - (10, 1) time:    Fri, 09.08.2017 14:09:27
Line 134: Line 289:
 
       2 -1 1 11 0 0 0.00999465 -0.00577041 0.131912 0.132417 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 --
 
  --- End of Generated Particles Bank --
 +
</pre>
  
−
GEMC2.4
+
===Event 4===
 +
<pre>
 +
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
 +
</pre>
 +
 
 +
GEMC 2.4 txt output
 +
<pre>
 
  --- Header Bank --
 
  --- Header Bank --
 
     - (10, 1) time:    Fri, 09.08.2017 13:28:02
 
     - (10, 1) time:    Fri, 09.08.2017 13:28:02
Line 227: Line 391:
 
  --- End of Generated Particles Bank --
 
  --- End of Generated Particles Bank --
 
  ---- End of Event  ----
 
  ---- End of Event  ----
 +
</pre>
  
  
−
GEMC2.6
+
GEMC 2.6 txt output
 +
<pre>
 
  --- Header Bank --
 
  --- Header Bank --
 
     - (10, 1) time:    Fri, 09.08.2017 14:09:27
 
     - (10, 1) time:    Fri, 09.08.2017 14:09:27
Line 326: Line 492:
 
  --- End of Generated Particles Bank --
 
  --- End of Generated Particles Bank --
 
  ---- End of Event  ----
 
  ---- End of Event  ----
 +
</pre>
 +
 +
=evio2rootAnalysis=
  
 +
The DC.xml files have changed
  
−
The arrangement of how the wire number and layers are read has 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
 
GEMC2.4
 +
<pre>
 
  EVENT # 3  GENPART = 2Detector
 
  EVENT # 3  GENPART = 2Detector
−
           Gen part 0 pid  = 11 px  = 8.89462e-25 py  = -1.44143e+29 pz  = 5.67653e-31 vx  = 0 vy  = 0 vz  = 0
+
           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  = -1.03587 py  = 0.965647 pz  = 2.58963 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 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 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
 
           Detector 2 dpid  = 11 procID  = 1 sector  = 6 superlayer  = 4 layer  = 2 wire  = 7
  
−
GEMC3.0
+
  EVENT # 23 GENPART = 2Detector
−
  EVENT # 3 GENPART = 2Detector
+
           Gen part 0 pid  = 11 px  = -10.4252 py  = 4.95027 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
−
           Gen part 0 pid  = 11 px  = 8.89462e-25 py  = -1.44143e+29 pz  = 5.67653e-31 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
−
           Gen part 1 pid  = 11 px  = -1.03587 py  = 0.965647 pz  = 2.58963 vx  = 0 vy  = 0 vz  = 0
+
           Detector 0 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 1 wire  = 1
−
           Detector 0 dpid  = 11 procID  = 1 sector  = 1 superlayer  = 25 layer  = 14 wire  = 150
+
           Detector 1 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 2 wire  = 1 
−
           Detector 1 dpid  = 11 procID  = 1 sector  = 1 superlayer  = 26 layer  = 15 wire  = 422
+
           Detector 2 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 3 wire  = 1 
−
           Detector 2 dpid  = 11 procID  = 1 sector  = 1 superlayer  = 27 layer  = 14 wire  = 64
+
           Detector 3 dpid  = 11 procID  = 90 sector  = 1 superlayer  = 5 layer  = 4 wire  = 1
−
           Detector 3 dpid  = 11 procID  = 1 sector  = 1 superlayer  = 28 layer  = 14 wire  = 489
+
 
−
           Detector 4 dpid  = 11 procID  = 1 sector  = 1 superlayer = 29 layer  = 14 wire  = 227
+
EVENT # 26  GENPART = 2Detector
−
           Detector 5 dpid  = 11 procID  = 1 sector  = 1 superlayer = 30 layer  = 14 wire  = -2147483648
+
          Gen part 0 pid  = 11 px  = -10.4765 py  = 4.84082 pz  = 10868.1 vx  = 0 vy  = 0 vz  = 0
−
           Detector 6 dpid  = 11 procID  = 1 sector  = 1 superlayer = 30 layer  = 15 wire  = 378
+
          Gen part 1 pid  = 11 px  = 10.4765 py  = -4.84082 pz  = 131.912 vx  = 0 vy  = 0 vz  = 0
−
           Detector 7 dpid  = 11 procID  = 1 sector  = 6 superlayer = 31 layer = 18 wire  = 449
+
</pre>
−
           Detector 8 dpid  = 11 procID  = 1 sector  = 6 superlayer = 32 layer = 19 wire  = 31
+
 
−
           Detector 9 dpid  = 11 procID  = 1 sector  = 6 superlayer = 33 layer = 19 wire  = 276
+
GEMC2.6
−
           Detector 10 dpid  = 11 procID  = 1 sector  = 6 superlayer = 34 layer  = 19 wire  = 208
+
<pre>
−
           Detector 11 dpid  = 11 procID  = 1 sector  = 6 superlayer = 35 layer  = 19 wire  = 425
+
EVENT # 3  GENPART = 2Detector
−
           Detector 12 dpid  = 11 procID  = 1 sector = 6 superlayer = 36 layer  = 19 wire  = 57
+
          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>
 
</pre>
  

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"