Difference between revisions of "HITs in DC"

From New IAC Wiki
Jump to navigation Jump to search
Line 8: Line 8:
 
   ! Version 2.6
 
   ! Version 2.6
 
|-
 
|-
   | A
+
   | Setup
   | B
+
   | Setup
 
|-
 
|-
   | C
+
   | Plots
   | D
+
   | 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

Caption
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


[math]\leftarrow[/math]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.

DeltaTheta v2 4.pngDeltaTheta v2 6.png

From 5 to 40 in Theta with Phi Isotropic

V2_6

Shield out

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


PhiThetaBinsTheoreticalOutV2 6.pngPhiThetaBinsOutV2 6.png



Shield In

<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"/>
PhiThetaBinsTheoreticalInV2 6.pngPhiThetaBinsInV2 6.png


Possible Fix

<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"/-->
PhiThetaBinsINV2 6Fix.png

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"/>
PhiThetaBinsTheoreticalOutV2 4.pngPhiThetaBinsOutV2 4.png



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"/>
PhiThetaBinsTheoreticalInV2 4.pngPhiThetaBinsInV2 4.png



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
[math]\frac{96264}{98854}=97.4\%[/math]