Difference between revisions of "Analysis of evio files"

From New IAC Wiki
Jump to navigation Jump to search
 
(76 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:
  
  
−
The hits on DC Sector 1 Superlayer 1 Layer 1 have an artificial bump near the middle due to the decreasing wire number as phi (or the parameter of the ellipse) moves away from its semi-major radii with respect to the ellipse of constant theta.
 
  
−
[[File:Layer1WireNumberBin0.png]][[File:Layer1WireNumberBin0Weighted.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
  
−
[[File:Layer1WireNumberBin150.png]][[File:Layer1WireNumberBin150Weighted.png]]
 
  
 +
Similarly, these can account for only the 1st interactions
  
−
[[File:Layer1WireNumberBin299.png]][[File:Layer1WireNumberBin299Weighted.png]]
+
[[File:1stHitsThetaBins.png]][[File:1stHitsThetaBinsWeighted.png]]
  
  
−
We can run a larger simulation or we can use the wire hit and phi angle to find theta within a certain range.  This information can be used to make a plot of number of hits per angle theta bin on which the differential cross section can be overlayed.  Theoretically the phi data would be obtained from additional layers in the sector as well as timing data.
+
With the edited gcard, we find 61102 hits
  
  
 +
Overlaying this with the expected differential cross-section expected from the LUND file
  
−
[[File:ThetaBins.png]][[File:ThetaBinsWeighted.png]]
+
<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=
  
−
Again, restricting hits to only the 1st hits within layer 1
+
<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
  
−
[[File:1stHitsThetaBins.png]][[File:1stHitsThetaBinsWeighted.png]]
+
> 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=
 
 
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 not have as many interactions on layer 1 of section 1 superlayer 1
+
[[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 75: Line 189:
  
 
[[File:1stHitsThetaBinsGEMC3.png]][[File:1stHitsThetaBinsWeightedGEMC3.png]]
 
[[File:1stHitsThetaBinsGEMC3.png]][[File:1stHitsThetaBinsWeightedGEMC3.png]]
 +
 +
==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>
 +
 +
--- 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 --
 +
</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:      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>
 +
 +
===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 --
 +
    - (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"