HITs in DC

From New IAC Wiki
Jump to navigation Jump to search

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
row two-summary
Plots Plots





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


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]