Difference between revisions of "HITs in DC"
Jump to navigation
Jump to search
| Line 8: | Line 8: | ||
! Version 2.6 | ! Version 2.6 | ||
|- | |- | ||
| − | | | + | | Setup |
| − | | | + | | Setup |
|- | |- | ||
| − | | | + | | Plots |
| − | | | + | | Plots |
|} | |} | ||
Revision as of 17:23, 14 December 2017
Before wall padding is taken into effect. Only the limits on phi as a function of theta is considered.
5 to 15 degrees in Theta at phi 0
| Version 2.4 | Version 2.6 |
|---|---|
| Setup | Setup |
| Plots | Plots |
Version 2.4
Setup
setenv JLAB_ROOT /home/vanwdani/src/CLAS/GEMC
setenv JLAB_VERSION 1.3
source $JLAB_ROOT/$JLAB_VERSION/ce/jlab.csh
setenv CLASROOT /home/vanwdani/src/CLAS/clas-root
#setenv LD_LIBRARY_PATH /home/vanwdani/src/CLAS/clas-root/lib:${LD_LIBRARY_PATH}
setenv CCDB_CONNECTION sqlite:///home/vanwdani/src/CLAS/DB/clas12.sqlite
setenv CCDB_DATABASE sqlite:///home/vanwdani/src/CLAS/DB/clas12.sqlite
> 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
clas12Tags
<detector name="home/vanwdani/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="home/vanwdani/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/dc/dc" factory="TEXT" variation="java"/>
Older banks
<detector name="../../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/dc/dc\" factory=\"TEXT\" variation=\"java\"/>
Comparing Banks
Running the first ten sets of LUND files, for a total of 10000 events, from 5 degrees to 15 degrees in </math>\theta</math> by increments of 0.001
./Block000
Using the Hit condition
if(dpid[k]==11 && mpid[k]==0 && mtid[k]==0 && otid[k]==2 && procID[k]==90 && sector[k]==1)
Using the evio2root program
~/src/CLAS/evioRootDV/bin/evioeroot v2_4 LH2_0Sol_0Tor_11Gev_0Phi_ShieldOut 10
Analyzing the hits for the older Banks
root -l evio2rootAnalysis_test.C > test_v2_4.out
Analyzing the hits for 4a.0.0
root -l evio2rootAnalysis_test.C > test_v2_4_clas12Tags.out
Comparing the differences
diff -y --width=190 test_v2_4_clas12Tags.out test_v2_4.out
This shows that the two methods look to have similar results
Processing evio2rootAnalysis_test.C... Processing evio2rootAnalysis_test.C...
{θ(Th)=5.0000, φ(Th)=0.00, θ(Hit)=5.0000 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=0} {θ(Th)=5.0000, φ(Th)=0.00, θ(Hit)=5.0000 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=0}
{θ(Th)=5.0010, φ(Th)=0.00, θ(Hit)=5.0009 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=1} {θ(Th)=5.0010, φ(Th)=0.00, θ(Hit)=5.0009 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=1}
{θ(Th)=5.0020, φ(Th)=0.00, θ(Hit)=5.0021 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=2} {θ(Th)=5.0020, φ(Th)=0.00, θ(Hit)=5.0021 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=2}
.
.
.
{θ(Th)=14.9970, φ(Th)=0.00, θ(Hit)=14.9961 , φ(Hit)=0.01 ,layer= 1 , wire=33 hit=9997} {θ(Th)=14.9970, φ(Th)=0.00, θ(Hit)=14.9961 , φ(Hit)=0.01 ,layer= 1 , wire=33 hit=9997}
{θ(Th)=14.9980, φ(Th)=0.00, θ(Hit)=14.9986 , φ(Hit)=0.00 ,layer= 1 , wire=33 hit=9998} {θ(Th)=14.9980, φ(Th)=0.00, θ(Hit)=14.9986 , φ(Hit)=0.00 ,layer= 1 , wire=33 hit=9998}
{θ(Th)=14.9990, φ(Th)=0.00, θ(Hit)=14.9973 , φ(Hit)=0.00 ,layer= 1 , wire=33 hit=9999} {θ(Th)=14.9990, φ(Th)=0.00, θ(Hit)=14.9973 , φ(Hit)=0.00 ,layer= 1 , wire=33 hit=9999}
Total hits=10000 Total hits=10000
Testing to see if there is any differences at all
diff -q test_v2_4_clas12Tags.out test_v2_4.out
NO DIFFERENCES
Analyzing Angle change
Using evio2RootAnalysis_AngleChange.C
root -l evio2rootAnalysis_AngleChange.C
<pre>
The total number of hits that deviate from the LUND value of Theta by 0.5 degrees is collected
<pre>
{Δθ=-0.00, Δφ(Th)=-0.00, hit=0}
{Δθ=0.00, Δφ(Th)=-0.00, hit=1}
{Δθ=-0.00, Δφ(Th)=-0.00, hit=2}
.
.
.
{Δθ=0.00, Δφ(Th)=-0.01, hit=9997}
{Δθ=-0.00, Δφ(Th)=-0.00, hit=9998}
{Δθ=0.00, Δφ(Th)=-0.00, hit=9999}
Total hits=10000
Total Change in abs(Theta) >= .5 degrees=0
Version 2.6
Setup
setenv JLAB_ROOT /home/lds/src/CLAS/GEMC #setenv JLAB_VERSION 1.3 #setenv JLAB_VERSION 2.0 setenv JLAB_VERSION 2.1 source $JLAB_ROOT/$JLAB_VERSION/ce/jlab.csh setenv CLARA_HOME /home/lds/src/CLAS/CLARA/myCLARA #setenv CCDB_CONNECTION mysql://clas12reader@clasdb.jlab.org/clas12 setenv CCDB_CONNECTION sqlite:///home/lds/src/CLAS/DB/clas12.sqlite setenv CCDB_DATABASE sqlite:///home/lds/src/CLAS/DB/clas12.sqlite
> 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
clas12Tags
<detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
<detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/dc/dc" factory="TEXT" variation="java"/>
Analyzing Hits
Running the first ten sets of LUND files, for a total of 10000 events, from 5 degrees to 15 degrees in </math>\theta</math> by increments of 0.001
./Block000
Using the Hit condition
if(dpid[k]==11 && mpid[k]==0 && mtid[k]==0 && otid[k]==2 && procID[k]==999 && sector[k]==1)
Using the evio2root program
~/src/CLAS/evioRootDV/bin/evioeroot v2_6 LH2_0Sol_0Tor_11Gev_0Phi_ShieldOut 10
Analyzing the hits for 4a.0.0
root -l evio2rootAnalysis_test.C > test_v2_6_clas12Tags.out
Analyzing Angle change
Using evio2RootAnalysis_AngleChange.C
root -l evio2rootAnalysis_AngleChange.C
<pre>
The total number of hits that deviate from the LUND value of Theta by 0.2 degrees is collected
<pre>
{Δθ=0.26, Δφ(Th)=-0.64, hit=0}
{Δθ=-0.09, Δφ(Th)=-2.68, hit=1}
{Δθ=0.15, Δφ(Th)=-0.88, hit=2}
.
.
.
{Δθ=-0.09, Δφ(Th)=-0.37, hit=9997}
{Δθ=-0.86, Δφ(Th)=-5.89, hit=9998}
{Δθ=0.00, Δφ(Th)=-0.33, hit=9999}
Total hits=9978
Total Change in abs(Theta) >= .5 degrees=1983
Compare 2.4 and 2.6
Comparing the differences
diff -y --width=190 LH2_0Sol_0Tor_11GeV_0Phi_ShieldOut_v2_6/test_v2_6_clas12Tags.out LH2_0Sol_0Tor_11GeV_0Phi_ShieldOut_v2_4_sqlite/test_v2_4_clas12Tags.out
This shows that the two methods look to have similar results, but 2.6 has less hits
Processing evio2rootAnalysis_test.C... Processing evio2rootAnalysis_test.C...
{θ(Th)=5.0000, φ(Th)=0.00, θ(Hit)=4.7438 , φ(Hit)=0.64 ,layer= 7 , wire=1 hit=0} | {θ(Th)=5.0000, φ(Th)=0.00, θ(Hit)=5.0000 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=0}
{θ(Th)=5.0010, φ(Th)=0.00, θ(Hit)=5.0928 , φ(Hit)=2.68 ,layer= 1 , wire=1 hit=1} | {θ(Th)=5.0010, φ(Th)=0.00, θ(Hit)=5.0009 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=1}
{θ(Th)=5.0020, φ(Th)=0.00, θ(Hit)=4.8534 , φ(Hit)=0.88 ,layer= 7 , wire=1 hit=2} | {θ(Th)=5.0020, φ(Th)=0.00, θ(Hit)=5.0021 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=2}
{θ(Th)=5.0030, φ(Th)=0.00, θ(Hit)=5.3135 , φ(Hit)=1.89 ,layer= 1 , wire=1 hit=3} | {θ(Th)=5.0030, φ(Th)=0.00, θ(Hit)=5.0029 , φ(Hit)=0.00 ,layer= 1 , wire=1 hit=3}
.
.
.
{θ(Th)=14.9970, φ(Th)=0.00, θ(Hit)=15.0912 , φ(Hit)=0.37 ,layer= 1 , wire=32 hit=9997} | {θ(Th)=14.9970, φ(Th)=0.00, θ(Hit)=14.9961 , φ(Hit)=0.01 ,layer= 1 , wire=33 hit=9997}
{θ(Th)=14.9980, φ(Th)=0.00, θ(Hit)=15.8566 , φ(Hit)=5.89 ,layer= 1 , wire=34 hit=9998} | {θ(Th)=14.9980, φ(Th)=0.00, θ(Hit)=14.9986 , φ(Hit)=0.00 ,layer= 1 , wire=33 hit=9998}
{θ(Th)=14.9990, φ(Th)=0.00, θ(Hit)=14.9983 , φ(Hit)=0.33 ,layer= 1 , wire=32 hit=9999} | {θ(Th)=14.9990, φ(Th)=0.00, θ(Hit)=14.9973 , φ(Hit)=0.00 ,layer= 1 , wire=33 hit=9999}
Total hits=9978 | Total hits=10000
Testing to see if there is any differences at all
diff -q LH2_0Sol_0Tor_11GeV_0Phi_ShieldOut_v2_6/test_v2_6_clas12Tags.out LH2_0Sol_0Tor_11GeV_0Phi_ShieldOut_v2_4_sqlite/test_v2_4_clas12Tags.out
Files LH2_0Sol_0Tor_11GeV_0Phi_ShieldOut_v2_6/test_v2_6_clas12Tags.out and LH2_0Sol_0Tor_11GeV_0Phi_ShieldOut_v2_4_sqlite/test_v2_4_clas12Tags.out differ
Some of the differences come from the difference in the LUND values of theta and the angle calculated from the average x, y, and z positions of the DC hit.
From 5 to 40 in Theta with Phi Isotropic
V2_6
gcard
<detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ec/ec" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/pcal/pcal" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/dc/dc" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ftof/ftof" factory="TEXT" variation="java"/>
<detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ec/ec" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/pcal/pcal" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/dc/dc" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ftof/ftof" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ft/ft" factory="TEXT" variation="FTOn"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/beamline/beamline" factory="TEXT" variation="FTOn"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/beamline/vacuumLine" factory="CAD"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/beamline/cadBeamline" factory="CAD"/>
<detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ec/ec" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/pcal/pcal" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/dc/dc" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ftof/ftof" factory="TEXT" variation="java"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/ft/ft" factory="TEXT" variation="FTOn"/> <detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/beamline/beamline" factory="TEXT" variation="FTOn"/> <!--detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/beamline/vacuumLine" factory="CAD"/--> <!--detector name="/home/lds/src/CLAS/GEMC/clas12Tags/4a.2.2/experiments/clas12/beamline/cadBeamline" factory="CAD"/-->
V2_4
Shield out
<detector name="../../../../../../clas12/ec/ec" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/pcal/pcal" factory="TEXT" variation="javageom"/> <detector name="../../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/dc/dc" factory="TEXT" variation="ccdb"/> <detector name="../../../../../../clas12/ftof/ftof" factory="TEXT" variation="java"/>
Shield In
<detector name="../../../../../../clas12/ec/ec" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/pcal/pcal" factory="TEXT" variation="javageom"/> <detector name="../../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/dc/dc" factory="TEXT" variation="ccdb"/> <detector name="../../../../../../clas12/ftof/ftof" factory="TEXT" variation="java"/> <detector name="../../../../../../clas12/beamline/beamline" factory="TEXT" variation="FTOn"/> <detector name="../../../../../../clas12/ft/ft" factory="TEXT" variation="FTOn"/>
insert molller weighted graphs + binned efficiency graph (Recons/Theory), a slice at phi=0 and phi=30 (2-D plot)
insert beamline shielding into g-card and graph rate difference (With shield - no shield) in units of Hz assuming 10 nA and 11 GeV