Difference between revisions of "HITs in DC"

From New IAC Wiki
Jump to navigation Jump to search
Line 33: Line 33:
  
  
 
===clas12Tags===
 
<pre>
 
<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"/>
 
</pre>
 
 
 
 
 
====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
 
<pre>
 
./Block000
 
</pre>
 
 
 
Using the Hit condition
 
<pre>
 
if(dpid[k]==11 && mpid[k]==0 && mtid[k]==0 && otid[k]==2 && procID[k]==999 && sector[k]==1)
 
</pre>
 
 
 
Using the evio2root program
 
<pre>
 
~/src/CLAS/evioRootDV/bin/evioeroot v2_6 LH2_0Sol_0Tor_11Gev_0Phi_ShieldOut 10
 
</pre>
 
 
 
 
 
 
Analyzing the hits for 4a.0.0
 
<pre>
 
root -l evio2rootAnalysis_test.C > test_v2_6_clas12Tags.out
 
</pre>
 
  
 
==Compare 2.4 and 2.6==
 
==Compare 2.4 and 2.6==

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





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]