Difference between revisions of "Analysis of evio files"
(35 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: | ||
− | These can be expressed in terms of Theta bins | + | |
+ | 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]] | [[File:ThetaBins.png]][[File:ThetaBinsWeighted.png]] | ||
+ | |||
+ | With the edited gcard we only have 68114 hits per Theta bins | ||
Line 40: | Line 68: | ||
+ | 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: | + | [[File:Layer1WireNumberBinTotalAdditional.png]] [[File:Layer1WireNumberBinTotalWeightedAdditional.png]] |
− | [[File:Layer1binsHitsAdditional.png]][[File:Layer1binsHitsWeightedAdditional.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 73: | 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). The databank for the DC has also eliminated the superlayer variable from the Digitialized output, moving up the wire number as a result. | ||
− | [[File: | + | [[File:Layer1BinHitsv2-6local.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: | + | [[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 98: | Line 190: | ||
[[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=== | ===Event 1=== | ||
Line 368: | Line 494: | ||
</pre> | </pre> | ||
− | + | =evio2rootAnalysis= | |
The DC.xml files have changed | The DC.xml files have changed | ||
Line 396: | Line 522: | ||
<pre> | <pre> | ||
EVENT # 3 GENPART = 2Detector | EVENT # 3 GENPART = 2Detector | ||
− | Gen part 0 pid = 11 px = | + | 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 = | + | 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 | ||
− | EVENT # 23 GENPART = 2Detector | + | EVENT # 23 GENPART = 2Detector |
− | Gen part 0 pid = 11 px = - | + | 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 = | + | 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 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 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 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 | Detector 3 dpid = 11 procID = 90 sector = 1 superlayer = 5 layer = 4 wire = 1 | ||
− | + | EVENT # 26 GENPART = 2Detector | |
− | Gen part 0 pid = 11 px = - | + | 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 = | + | Gen part 1 pid = 11 px = 10.4765 py = -4.84082 pz = 131.912 vx = 0 vy = 0 vz = 0 |
</pre> | </pre> | ||
GEMC2.6 | GEMC2.6 | ||
<pre> | <pre> | ||
− | + | EVENT # 3 GENPART = 2Detector | |
− | Gen part 0 pid = 11 px = | + | 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 = | + | 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 0 dpid = 11 procID = 1 sector = 1 layer = 25 wire = 14 | ||
Detector 1 dpid = 11 procID = 1 sector = 1 layer = 26 wire = 15 | Detector 1 dpid = 11 procID = 1 sector = 1 layer = 26 wire = 15 | ||
Line 435: | Line 561: | ||
EVENT # 23 GENPART = 2Detector | EVENT # 23 GENPART = 2Detector | ||
− | Gen part 0 pid = 11 px = - | + | 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 = | + | 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 = - | + | 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 = | + | 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 0 dpid = 11 procID = 999 sector = 1 layer = 25 wire = 1 | ||
Detector 1 dpid = 11 procID = 999 sector = 1 layer = 26 wire = 1 | Detector 1 dpid = 11 procID = 999 sector = 1 layer = 26 wire = 1 |
Latest revision as of 02:50, 14 October 2017
Using the evio2root converter
~/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
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.
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.
Similarly, editing the gcard only results in 92761 1st hits via wires
These wire hits can be expressed in terms of Theta bins
With the edited gcard we only have 68114 hits per Theta bins
Similarly, these can account for only the 1st interactions
With the edited gcard, we find 61102 hits
Overlaying this with the expected differential cross-section expected from the LUND file
![EfficiencyThetaBins.png](/./images/b/be/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](/./images/b/b5/Ratio.png)
This can be partially explained in that Superlayer 1, Layer 1 is not always hit first
Lower Fix
Determining what values of phi are not within the sector following the Mathematica equation developed earlier
Upper Fix
Fixing the upper tail of hits on wires from larger theta values
Combined Results
Combining these plots
hadd combined.root NonIsotropic_BinAdditions.root LH2_0Sol_0Tor_11GeV_IsotropicPhi_ShieldOut.root
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
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.
Restricting results to only the 1st hits on this layer
This lack of interactions translates to angle bins as well
Restricting results to only the 1st hits on this layer
Upper fix
Combined Results
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
The GEANT simulation results look to be identical