Difference between revisions of "HITs in DC"

From New IAC Wiki
Jump to navigation Jump to search
Line 68: Line 68:
 
<pre>
 
<pre>
 
root -l evio2rootAnalysis_test.C > test_v2_6_clas12Tags.out
 
root -l evio2rootAnalysis_test.C > test_v2_6_clas12Tags.out
−
</pre>
 
−
 
−
===Analyzing Angle change===
 
−
Using evio2RootAnalysis_AngleChange.C
 
−
 
−
 
−
<pre>
 
−
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
 
 
</pre>
 
</pre>
  

Revision as of 18:05, 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

5 to 15 degrees in Theta at phi 0 GEMC Version Tests
Version 2.4 Version 2.6
GEMC V2.4 Trial Setup Conditions GEMC V2.6 Trial Setup Conditions
GEMC V2.4 JLAB 1.3 Banks GEMC V2.6 clas12Tags Banks
GEMC V2.4 clas12Tags Banks
GEMC V2.4 Bank Comparison
GEMC V2.4 Angle Change GEMC V2.6 Angle Change
Plots Plots





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

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]
Retrieved from "https://wiki.iac.isu.edu/index.php?title=HITs_in_DC&oldid=120209"