TD Ddoverd 2009
5-01-2009
Target and Beam Polarization are positive
Helcode # | Negative Beam Torus | Positive Beam Torus |
1 h>0 | ||
2 h<0 | ||
3 | ||
4 |
Target Polarization Positive and Beam Polarization Negative
Helcode # | Negative Beam Torus | Positive Beam Torus |
1 h>0 | No Data | |
2 h<0 | No Data | |
3 | No Data | |
4 | No Data |
Inclusive Histograms For Invariant Mass
Invariant Mass Histograms for Each Sector
All helicities
Helcode=1
3/20/09
1.) plot Wdiff
2.) plot PE using Osipenko/Josh cuts
06/6/09
Invariant Mass
Difference of Invariant Mass for two differnet cuts:
Run Number | EC Cuts+ requiring pion | OSI Cuts | EC Cuts | |
26991 |
On x-axis of FC difference between the "+" heliciy FC and the "-" helicity FC.
Run Number | W difference | FC difference | End of Run sum |
26991 | 150px | 0.00354260538 | 0.00001050400|
26994 | 0.00367 | 0.00001088263||
26993 | 0.002559223 | 0.00001120796||
26992 | 0.003781 | 0.00001090044||
26990 | 0.00331986321 | 0.00001203387||
26989 | 0.00303841557 | 0.00001267152||
26988 | 0.00380169243 | 0.00001097470||
26987 | 0.00280163240 | 0.00001515270||
26986 | 0.00352882616 | 0.00000048897||
26979 | -0.00373104126 | 0.00001531706||
26965 | 0.00454140439 | 0.00001534923||
26964 | 0.00509434283 | 0.00001557226||
26963 | 0.00405285155 | 0.00001575049||
26961 | 0.00507456102 | 0.00001682744||
26959 | 0.00523761694 | 0.00001488890||
26958 | 0.00444205478 | 0.00001577510||
26956 | 0.00504401620 | 0.00001565456||
26955 | 0.00559913829 | 0.00001425736||
26954 | 0.00494690728 | 0.00001443107||
26953 | 0.00562598448 | 0.00001525797||
26952 | 0.00485295834 | 0.00001293033||
26951 | 0.00545485906 | 0.00001558637||
26947 | 0.00465697160 | 0.00001144285||
26945 | -0.00556456389 | 0.00001389870||
26943 | -0.00656109472 | 0.00001540342||
26942 | -0.00669169106 | 0.00001526771||
26941 | -0.00676460417 | 0.00001100372||
26940 | -0.00643597022 | 0.00001397207||
26939 | -0.00650047456 | 0.00001405001||
26938 | -0.00632630313 | 0.00001411839||
26937 | -0.00533169794 | 0.00001417768||
26934 | -0.00601497156 | 0.00001422560||
26933 | -0.00586540500 | 0.00001414938||
26932 | -0.00587687192 | 0.00001383218||
26931 | -0.00630232566 | 0.00001495053||
26930 | -0.00684205010 | 0.00001541629||
26929 | -0.00521666944 | 0.00001555256||
26928 | -0.00682185316 | 0.00001536552||
26927 | -0.00599647483 | 0.00001538040||
26926 | -0.00649932378 | 0.00001302331||
26925 | -0.00575464099 | 0.00001530219||
27079 | -0.01025869470 | 0.00002322298||
27078 | -0.01128036266 | 0.00002135573||
27075 | -0.01037443017 | 0.00002113646||
27109 | -0.00876887465 | 0.00002411331||
27107 | 0.01044530752 | 0.00002188631||
27116 | 0.00173242954 | 0.00003084415||
27112 | -0.01142682761 | 0.00002505827||
27111 | -0.00872760886 | 0.00002117268||
27128 | 0.00193441756 | 0.00002177434||
27127 | 0.00215712098 | 0.00002186032||
27124 | 0.00309070808 | 0.00002235601||
27139 | -0.00288862444 | 0.00002683943||
27138 | -0.00061608634 | 0.00002396188||
27137 | -0.00275066778 | 0.00002150471||
27136 | -0.00239423237 | 0.00002151354||
27134 | -0.00355581434 | 0.00002051114||
27133 | -0.00252193724 | 0.00002181781||
27132 | 0.00096332202 | 0.00002224348||
27143 | -0.00210486278 | 0.00002838738||
27141 | -0.00301990556 | 0.00002515365||
27160 | 0.00171264744 | 0.00002170610||
27161 | 0.00224287304 | 0.00002374488||
27162 | 0.00118844445 | 0.00002480837||
27166 | 0.00077237006 | 0.00002247548||
27167 | 0.00216041281 | 0.00002188669||
27168 | 0.00083650625 | 0.00002173328||
27170 | 0.00151247511 | 0.00002430185||
27175 | -0.01217485668 | 0.00002235229||
27176 | -0.01049635167 | 0.00002127680||
27177 | -0.01147445288 | 0.00002149007||
27179 | -0.01025340153 | 0.00002206557||
27180 | -0.00868272418 | 0.00002169057||
27181 | -0.00973324188 | 0.00002374079||
27182 | -0.00983862460 | 0.00002233926||
27183 | -0.01901615252 | 0.00002135060||
27186 | 0.01079090635 | 0.00002093983||
27187 | 0.01155519076 | 0.00002157075||
27188 | 0.01086898605 | 0.00002089859||
27190 | 0.01181266490 | 0.00002159430||
27192 | -0.00988459637 | 0.00002151795||
27193 | -0.01047643036 | 0.00002159742||
27194 | -0.01185520003 | 0.00002134810
NPHE
To find out pion contamination in the electron sample i used Osipenko geometrical cuts. The number of photoelectrons before and after osipenko cuts are shown below:
No cuts | OSI Cuts |
For different fits:
06/11/09
Pion Contamination
No cuts | OSI Cuts (Gauss(0)+Landau(3)+Gauss(6)) | OSI Cuts (Gauss(0)+Gauss(3)) | OSICuts + NPHE>2.5 (Gauss) |
OSICuts (Gauss(0)+Landau(3)+Gauss(6))
Assuming that the two gaussians represent number of photoelectrons and landau number of photons produced by high energy pions, the ratio of number of pions over the sum of electrons and pion in the electron candidate sample can be calculated in the following way:
- OSIcut
- OSICut+NPHE>2.5
OSICuts (Gauss(0)+Gauss(3))
In case of only two gaussians, the number of photoelectrons produced by pions is described by Gauss(0) and the number of photoelectrons created by electrons is Gauss(3). Pion contamination is calculated below for both cases, without and with NPHE>2.5 cut.
- OSICut
- With NPHE>2.5 Cut
Number of Events after NPHE>2.5 Cut
Counts in FCup
6/12/09
1.) Improve Chi^2 in NPe fits.
2.) Calculate uncertainty in pion contamination measurement by changing mean and widths according to fit error.
3.) Pulse pair FC asymmetry, and End of Run accumulated FC asym.
pulse pair
End of Run sum
4.) Determine semi-inclusive statistic as function of X
Uncertainty in Pion Contamination
Maximum
Minimum
Pion Contamination
It appears that pion contamination in electron sample is 9.63 %
0.01 % before nphe cut and after nphe>2.5 cut contamination is about 4.029% 0.003.X_bjorken
- 1). alldataOSICuts_X.root - OSICuts applied.
- 2). alldataOSICuts_X_epx.root - OSICuts applied and electron and pion are required.
- 3). alldataOSICuts_X_epxnphe.root - OSICut and nphe>2.5 cuts applied and electron and pion are required.
- 4). alldataX_epxwithoutcuts - No cuts, electron and pion required.
- Number of Events after cuts
No Cuts | OSI Cuts | OSI+NPHE>2.5 Cuts |
68.5 % | 57.5 % |
- Error Calculation
Tthe error in the asymmetry measurement would be
X_b | X_b Asymmetry | Error | ||
0.1 | 0.00087251693 | |||
0.2 | 0.0008507530911145 | |||
0.3 | 1.0691459e-03 | |||
0.4 | 0.0014004231 | |||
0.5 | 0.0018665742 | |||
0.6 | 0.0016477095 | |||
0.7 | 0.0022190018 | |||
0.8 | 0.00291609 | |||
0.9 | 0.003592967 | |||
1 | 0.0040928449 |
I am pretty sure X_{BJ} > 0.8 is not possible with our data set
Electron theta angle and cuts
- electron theta angle for different X_b
X_b when
Number of Events for X_b>0.8
plot the vertex of the above hits with X>0.8
change below to log plots so we can see where XBj stops
10/23/09
After Months of working on detectors and writing thesis proposal it is now time to start doing some physics.
1.) Determine how pion contamination uncertainty changes when you change fit parameters by 1 S.D., 2 S.D., and 3 S.D.
2.) FC asymm plots
3.) Vertex plot for X > 0.8 events.
4.) Now that we have good electron cuts. Plot statistics for Pion cuts.
5.) After pion cuts we start looking add paddle efficiencies so we can subtract sem-inclusive rates using individual paddles but opposite magnetic fields.
Xbjorken
Vertex plot for X > 0.8 events and others
I chose 1<Q^2<4 cut because we used it to plot phi angle in cm frame vs relative rate to compare with the results in paper.
Cuts | X_bjorken | Vertex X | Vertex Y | Vertex Z | |
OSI Cuts + EC Cuts | |||||
OSI Cuts + EC Cuts + X_b>0.8 | |||||
OSI Cuts + EC Cuts + X_b<0.8 | |||||
OSI Cuts + EC Cuts + 1<Q^2<4 |
Cuts | The scattered electron energy | electron scattering angle |
FC Asymmetry
FC Asymmetry plot using the following method : End of Run sum
Pion contamination
Determine how pion contamination uncertainty changes when you change fit parameters by 1 S.D., 2 S.D., and 3 S.D.
Pion contamination(3 S.D.) in electron sample is 9.645 %
0.025%.It doesnt really change from using 1 S.D.
Pion Statistics
Before and after cuts the plot of EC_tot/P vs nphe(for pions)
10/30/09
After Months of working on detectors and writing thesis proposal it is now time to start doing some physics.
1.) Determine how pion contamination uncertainty changes when you change fit parameters by 1 S.D., 2 S.D., and 3 S.D.
2.) Do pulse pair FC asymm plot
3.) Check program's calculation of event with X > 0.8 events. and compare to similar event with X < 0.8
4.) Use statistics for Pion cuts to estimate SIDIS statistical error -vs- Xbj
5.) After pion cuts we start looking add paddle efficiencies so we can subtract sem-inclusive rates using individual paddles but opposite magnetic fields.
1.)
1.) Determine how pion contamination uncertainty changes when you change fit parameters by 1 S.D., 2 S.D., and 3 S.D.
In case of 10 S.D. :
&&3.)
I suspect the X_b >0.8 event below are pions mis-identified as electrons To figure out. Write down event number for events below as well as run number and file name. The use path length and Scintillator TDC time to determine beta under assumption that particle is a pion. Does the momentum and energy make sense? Download the cooked data file from JLab for these events so we can use CED to look at them and bosdump to look at the reconstruction.
- 1
Ebeam=5736 IBeam=4.2 ITorus=2248 ITarg=122 BeamPol=0.71 TargetPol=-0.67 BadRun=0 Target=18 PolPlate=0 Version=2 Prescalers:0:0:0:0:0:0:0 dump=13 W= 1.3504 Q= 4.78723 final electron energy= 2.6823 initial electron energy= 5.736 electron theta angle= 32.3896
- 2
Ebeam=5736 IBeam=4.2 ITorus=2248 ITarg=122 BeamPol=0.71 TargetPol=-0.67 BadRun=0 Target=18 PolPlate=0 Version=2 Prescalers:0:0:0:0:0:0:0 dump=13 W= 1.41085 Q= 5.1186 final electron energy= 2.41678 initial electron energy= 5.736 electron theta angle= 35.3749
- 3
Ebeam=5736 IBeam=4.2 ITorus=2248 ITarg=122 BeamPol=0.71 TargetPol=-0.67 BadRun=0 Target=18 PolPlate=0 Version=2 Prescalers:0:0:0:0:0:0:0 dump=13 W= 1.2575 Q= 4.83661 final electron energy= 2.7851 initial electron energy= 5.736 electron theta angle= 31.9378
- 1
W= 3.07188 Q= 0.229403 final electron energy= 1.0543 initial electron energy= 5.736 electron theta angle= 11.177
- 2
hit return for next event, q to quit: W= 2.48202 Q= 1.382 final electron energy= 2.18587 initial electron energy= 5.736 electron theta angle= 19.1106
- 3
hit return for next event, q to quit: W= 2.92788 Q= 0.274895 final electron energy= 1.49046 initial electron energy= 5.736 electron theta angle= 10.2878
4.)
&& | && | && | && |
1-12-09
Root files: alldatasector26990_4.root, alldatasector26990_5.root, alldatasector27113_4.root, alldatasector27113_5.root.
9/13/09
1.) Determine how pion contamination uncertainty changes when you change fit parameters.
Fit parameters | Pion Contamination |
1 S. D. | 9.63 % | 0.01 %
3 S. D. | 9.645 % | 0.025 %
10 S. D. | 9.446 % | 0.0233 %
2.) Do pulse pair FC asymm plot
I did it for one file(dst27113_00.B00) and it was zero.
3.) Check program's calculation of event with X > 0.8 events. and compare to similar event with X < 0.8
I suspect the X_b >0.8 event below are pions mis-identified as electrons To figure out. Write down event number for events below as well as run number and file name. The use path length and Scintillator TDC time to determine beta under assumption that particle is a pion. Does the momentum and energy make sense? Download the cooked data file from JLab for these events so we can use CED to look at them and bosdump to look at the reconstruction.
Calculation is right, need to check CED, but dont have it on daq.
4.) Use statistics for Pion cuts to estimate SIDIS statistical error -vs- Xbj
Insert table with X bj, number of reconstructed pions, statistical error.
no cut
root alldatasector27113_5_1.root
root [9] .p X_bjorken->GetBinError(2); (const Double_t)3.50713558335003626e+01 root [10] .p X_bjorken->GetBinContent(2); (const Double_t)1.23000000000000000e+03
root [11] .p X_bjorken->GetBinContent(3); (const Double_t)1.85200000000000000e+03 root [12] .p X_bjorken->GetBinError(3); (const Double_t)4.30348695827000256e+01
root [13] .p X_bjorken->GetBinError(4); (const Double_t)3.06431068920891256e+01 root [14] .p X_bjorken->GetBinContent(4); (const Double_t)9.39000000000000000e+02
root [15] .p X_bjorken->GetBinContent(5); (const Double_t)3.68000000000000000e+02 root [16] .p X_bjorken->GetBinError(5); (const Double_t)1.91833260932508765e+01
(const Double_t)1.91833260932508765e+01 root [17] .p X_bjorken->GetBinError(6); (const Double_t)1.14017542509913792e+01 root [18] .p X_bjorken->GetBinContent(6); (const Double_t)1.30000000000000000e+02
root [19] .p X_bjorken->GetBinContent(7); (const Double_t)4.10000000000000000e+01 root [20] .p X_bjorken->GetBinError(7); (const Double_t)6.40312423743284853e+00
root [21] .p X_bjorken->GetBinError(8); (const Double_t)3.74165738677394133e+00 root [22] .p X_bjorken->GetBinContent(8); (const Double_t)1.40000000000000000e+01
root [23] .p X_bjorken->GetBinContent(9); (const Double_t)5.00000000000000000e+00 root [24] .p X_bjorken->GetBinError(9); (const Double_t)2.23606797749978981e+00
root [25] .p X_bjorken->GetBinError(10); (const Double_t)1.41421356237309515e+00 root [26] .p X_bjorken->GetBinContent(10); (const Double_t)2.00000000000000000e+00
with cut
root alldatasector27113_5.root
root [3] .p X_bjorken->GetBinError(2); (const Double_t)2.89827534923788761e+01 root [4] .p X_bjorken->GetBinContent(2); (const Double_t)8.40000000000000000e+02
root [5] .p X_bjorken->GetBinError(3); (const Double_t)3.77491721763537456e+01 root [6] .p X_bjorken->GetBinContent(3); (const Double_t)1.42500000000000000e+03
root [7] .p X_bjorken->GetBinError(4); (const Double_t)2.82488937836510701e+01 root [8] .p X_bjorken->GetBinContent(4); (const Double_t)7.98000000000000000e+02 root [9] .p X_bjorken->GetBinError(5); (const Double_t)1.81107702762748346e+01 root [10] .p X_bjorken->GetBinContent(5); (const Double_t)3.28000000000000000e+02
root [11] .p X_bjorken->GetBinError(6); (const Double_t)1.04880884817015154e+01 root [12] .p X_bjorken->GetBinContent(6); (const Double_t)1.10000000000000000e+02
root [13] .p X_bjorken->GetBinError(7); (const Double_t)6.24499799839839831e+00 root [14] .p X_bjorken->GetBinContent(7); (const Double_t)3.90000000000000000e+01
root [15] .p X_bjorken->GetBinError(8); (const Double_t)3.60555127546398912e+00 root [16] .p X_bjorken->GetBinContent(8); (const Double_t)1.30000000000000000e+01
root [17] .p X_bjorken->GetBinError(9); (const Double_t)1.73205080756887719e+00 root [18] .p X_bjorken->GetBinContent(9); (const Double_t)3.00000000000000000e+00
root [19] .p X_bjorken->GetBinError(10); (const Double_t)0.00000000000000000e+00 root [20] .p X_bjorken->GetBinContent(10); (const Double_t)0.00000000000000000e+00
root [21] .p X_bjorken->GetBinError(11); (const Double_t)0.00000000000000000e+00 root [22] .p X_bjorken->GetBinContent(11); (const Double_t)0.00000000000000000e+00
root [23] .p X_bjorken->GetBinError(12); (const Double_t)0.00000000000000000e+00 root [24] .p X_bjorken->GetBinContent(12); (const Double_t)0.00000000000000000e+00
Error | |
0.1 | 0.034503278 |
0.2 | 0.026490647 |
0.3 | 0.035399616 |
0.4 | 4.255675067 |
0.5 | 0.095346259 |
0.6 | 0.160128154 |
0.7 | 0.277350098 |
0.8 | 0.577350269 |
5.) After pion cuts we start looking add paddle efficiencies so we can subtract sem-inclusive rates using individual paddles but opposite magnetic fields.
Your B-field sign change does effect paddle distribution?
The table below represents the distribution of electrons and pions on the scintillator paddles using the reaction
e(p/d,e')\pi X
File number | electron | electron | ||
26990, B<0 | ||||
27113, B>0 |
- 1.) 7< Sector_paddle <11 - (B>0, ) - 15.3% && (B<0, ) - 10.9%
- 2.) 25< Sector_paddle <29 - (B>0, ) - 7.775% && (B<0, ) - 10.97%
&&
Using Two runs: 26990(NH3, -2250) and 27124(ND3, +2250)
root files: alldatasector27124_4.root && alldatasector27124_5.root
root files: alldatasector26990_4.root && alldatasector26990_5.root
alldatasector26990_4.root
Pion Paddle Number
Negative Torus
Using runs with NH3 target
Detected particles in the final state | Pion Paddle Number | X_b vs pion paddle number | Chosen pion paddle number |
&& | 7 | ||
&& | 27 |
Detected particles in the final state and chosen pion paddle number | paddle number | vs paddle number | vs paddle number when | vs paddle number when |
&& , | ||||
&& , |
Positive Torus
Using runs with NH3 target
Detected particles in the final state | X_b vs pion paddle number | Chosen pion paddle number |
&& | 27 | |
&& | 7 |
Detected particles in the final state and chosen pion paddle number | paddle number | vs paddle number | vs paddle number when | vs paddle number when |
&& , | ||||
&& , |