TD Ddoverd 2009

From New IAC Wiki
Revision as of 18:13, 6 February 2011 by Foretony (talk | contribs) (Created page with '=5-01-2009= ==Target and Beam Polarization are positive== {| border="1" |cellpadding="20" cellspacing="0 |- | Helcode # || Negative Beam Torus || Positive Beam Torus |- | 1 …')
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

5-01-2009

Target and Beam Polarization are positive

Helcode # Negative Beam Torus Positive Beam Torus
1 h>0 BT negative & helcode 1 phi angle cm frame PbPt positive.gif BT positive & helcode 1 phi angle cm frame PbPt positive.gif
2 h<0 BT negative & helcode 2 phi angle cm frame PbPt positive.gif BT positive & helcode 2 phi angle cm frame PbPt positive.gif
3 BT negative & helcode 3 phi angle cm frame PbPt positive.gif BT positive & helcode 3 phi angle cm frame PbPt positive.gif
4 BT negative & helcode 4 phi angle cm frame PbPt positive.gif BT positive & helcode 4 phi angle cm frame PbPt positive.gif



Target Polarization Positive and Beam Polarization Negative

Helcode # Negative Beam Torus Positive Beam Torus
1 h>0 No Data BT positive & helcode 1 phi angle cm frame Pt positive Pb negative.gif
2 h<0 No Data BT positive & helcode 2 phi angle cm frame Pt positive Pb negative.gif
3 No Data BT positive & helcode 3 phi angle cm frame Pt positive Pb negative.gif
4 No Data BT positive & helcode 4 phi angle cm frame Pt positive Pb negative.gif


Inclusive Histograms For Invariant Mass

Invariant Mass Histograms for Each Sector

All helicities

InvariantMass sector 1.gifInvariantMass sector 2.gif
InvariantMass sector 3.gifInvariantMass sector 4.gif
InvariantMass sector 5.gifInvariantMass sector 6.gif


InvariantMass sector 1 gaussian.gifInvariantMass sector 2 gaussian.gif
InvariantMass sector 3 gaussian.gifInvariantMass sector 4 gaussian.gif
InvariantMass sector 5 gaussian.gifInvariantMass sector 6 gaussian.gif

Helcode=1

InvariantMass sector 1 H1.gifInvariantMass sector 2 H1.gif
InvariantMass sector 3 H1.gifInvariantMass sector 4 H1.gif
InvariantMass sector 5 H1.gifInvariantMass sector 6 H1.gif

3/20/09

InvariantMass W difference Coupleoffiles 26992.gif

InvariantMass W difference Coupleoffiles 26993.gif

InvariantMass W difference Coupleoffiles 26994.gif

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 Missing mass difference RunNumber26991 1.gif Missing mass difference RunNumber26991 1 OSICuts.gif

Missing mass difference RunNumber26991 1ECuts.gif


On x-axis of FC difference between the "+" heliciy FC  and the "-" helicity FC.
Run Number W difference FC difference End of Run sum [math]\equiv \frac {\sum_i(FC(i)^+) - \sum(FC(i)^-)}{\sum (FC(i)^+) + \sum(FC(i)^-)})[/math]
26991 Missing mass difference RunNumber26991 1 OSICuts.gif 150px 0.00354260538 [math]\pm[/math] 0.00001050400
26994 Missing mass difference RunNumber26994 1 OSICuts.gif FC difference RunNumber26994 OSICuts.gif 0.00367 [math]\pm[/math] 0.00001088263
26993 Missing mass difference RunNumber26993 1 OSICuts.gif FC difference RunNumber26993 OSICuts.gif 0.002559223 [math]\pm[/math] 0.00001120796
26992 Missing mass difference RunNumber26992 1 OSICuts.gif FC difference RunNumber26992 OSICuts.gif 0.003781[math]\pm[/math]0.00001090044
26990 Missing mass difference RunNumber26990 1 OSICuts.gif FC difference RunNumber26990 OSICuts.gif 0.00331986321[math]\pm[/math]0.00001203387
26989 Missing mass difference RunNumber26989 1 OSICuts.gif FC difference RunNumber26989 OSICuts.gif 0.00303841557[math]\pm[/math]0.00001267152
26988 Missing mass difference RunNumber26988 1 OSICuts.gif FC difference RunNumber26988 OSICuts.gif 0.00380169243[math]\pm[/math]0.00001097470
26987 Missing mass difference RunNumber26987 1 OSICuts.gif FC difference RunNumber26987 OSICuts.gif 0.00280163240[math]\pm[/math]0.00001515270
26986 Missing mass difference RunNumber26986 1 OSICuts.gif FC difference RunNumber26986 OSICuts.gif 0.00352882616[math]\pm[/math]0.00000048897
26979 Missing mass difference RunNumber26979 1 OSICuts.gif FC difference RunNumber26979 OSICuts.gif -0.00373104126[math]\pm[/math]0.00001531706
26965 Missing mass difference RunNumber26965 1 OSICuts.gif FC difference RunNumber26965 OSICuts.gif 0.00454140439[math]\pm[/math]0.00001534923
26964 Missing mass difference RunNumber26964 1 OSICuts.gif FC difference RunNumber26964 OSICuts.gif 0.00509434283[math]\pm[/math]0.00001557226
26963 Missing mass difference RunNumber26963 1 OSICuts.gif FC difference RunNumber26963 OSICuts.gif 0.00405285155[math]\pm[/math]0.00001575049
26961 Missing mass difference RunNumber26961 1 OSICuts.gif FC difference RunNumber26961 OSICuts.gif 0.00507456102 [math]\pm[/math]0.00001682744
26959 Missing mass difference RunNumber26959 1 OSICuts.gif FC difference RunNumber26959 OSICuts.gif 0.00523761694[math]\pm[/math]0.00001488890
26958 Missing mass difference RunNumber26958 1 OSICuts.gif FC difference RunNumber26958 OSICuts.gif 0.00444205478[math]\pm[/math] 0.00001577510
26956 Missing mass difference RunNumber26956 1 OSICuts.gif FC difference RunNumber26956 OSICuts.gif 0.00504401620[math]\pm[/math]0.00001565456
26955 Missing mass difference RunNumber26955 1 OSICuts.gif FC difference RunNumber26955 OSICuts.gif 0.00559913829[math]\pm[/math]0.00001425736
26954 Missing mass difference RunNumber26954 1 OSICuts.gif FC difference RunNumber26954 OSICuts.gif 0.00494690728[math]\pm[/math]0.00001443107
26953 Missing mass difference RunNumber26953 1 OSICuts.gif FC difference RunNumber26953 OSICuts.gif 0.00562598448[math]\pm[/math]0.00001525797
26952 Missing mass difference RunNumber26952 1 OSICuts.gif FC difference RunNumber26952 OSICuts.gif 0.00485295834[math]\pm[/math]0.00001293033
26951 Missing mass difference RunNumber26951 1 OSICuts.gif FC difference RunNumber26951 OSICuts.gif 0.00545485906[math]\pm[/math]0.00001558637
26947 Missing mass difference RunNumber26947 1 OSICuts.gif FC difference RunNumber26947 OSICuts.gif 0.00465697160[math]\pm[/math]0.00001144285
26945 Missing mass difference RunNumber26945 1 OSICuts.gif FC difference RunNumber26945 OSICuts.gif -0.00556456389[math]\pm[/math] 0.00001389870
26943 Missing mass difference RunNumber26943 1 OSICuts.gif FC difference RunNumber26943 OSICuts.gif -0.00656109472[math]\pm[/math]0.00001540342
26942 Missing mass difference RunNumber26942 1 OSICuts.gif FC difference RunNumber26942 OSICuts.gif -0.00669169106[math]\pm[/math]0.00001526771
26941 Missing mass difference RunNumber26941 1 OSICuts.gif FC difference RunNumber26941 OSICuts.gif -0.00676460417[math]\pm[/math]0.00001100372
26940 Missing mass difference RunNumber26940 1 OSICuts.gif FC difference RunNumber26940 OSICuts.gif -0.00643597022 [math]\pm[/math]0.00001397207
26939 Missing mass difference RunNumber26939 1 OSICuts.gif FC difference RunNumber26939 OSICuts.gif -0.00650047456[math]\pm[/math]0.00001405001
26938 Missing mass difference RunNumber26938 1 OSICuts.gif FC difference RunNumber26938 OSICuts.gif -0.00632630313[math]\pm[/math]0.00001411839
26937 Missing mass difference RunNumber26937 1 OSICuts.gif FC difference RunNumber26937 OSICuts.gif -0.00533169794[math]\pm[/math]0.00001417768
26934 Missing mass difference RunNumber26934 1 OSICuts.gif FC difference RunNumber26934 OSICuts.gif -0.00601497156[math]\pm[/math]0.00001422560
26933 Missing mass difference RunNumber26933 1 OSICuts.gif FC difference RunNumber26933 OSICuts.gif -0.00586540500[math]\pm[/math]0.00001414938
26932 Missing mass difference RunNumber26932 1 OSICuts.gif FC difference RunNumber26932 OSICuts.gif -0.00587687192[math]\pm[/math]0.00001383218
26931 Missing mass difference RunNumber26931 1 OSICuts.gif FC difference RunNumber26931 OSICuts.gif -0.00630232566[math]\pm[/math]0.00001495053
26930 Missing mass difference RunNumber26930 1 OSICuts.gif FC difference RunNumber26930 OSICuts.gif -0.00684205010[math]\pm[/math]0.00001541629
26929 Missing mass difference RunNumber26929 1 OSICuts.gif FC difference RunNumber26929 OSICuts.gif -0.00521666944[math]\pm[/math]0.00001555256
26928 Missing mass difference RunNumber26928 1 OSICuts.gif FC difference RunNumber26928 OSICuts.gif -0.00682185316[math]\pm[/math]0.00001536552
26927 Missing mass difference RunNumber26927 1 OSICuts.gif FC difference RunNumber26927 OSICuts.gif -0.00599647483[math]\pm[/math]0.00001538040
26926 Missing mass difference RunNumber26926 1 OSICuts.gif FC difference RunNumber26926 OSICuts.gif -0.00649932378[math]\pm[/math]0.00001302331
26925 Missing mass difference RunNumber26925 1 OSICuts.gif FC difference RunNumber26925 OSICuts.gif -0.00575464099[math]\pm[/math]0.00001530219
27079 Missing mass difference RunNumber27079 1 OSICuts.gif FC difference RunNumber27079 OSICuts.gif -0.01025869470[math]\pm[/math]0.00002322298
27078 Missing mass difference RunNumber27078 1 OSICuts.gif FC difference RunNumber27078 OSICuts.gif -0.01128036266[math]\pm[/math]0.00002135573
27075 Missing mass difference RunNumber27075 1 OSICuts.gif FC difference RunNumber27075 OSICuts.gif -0.01037443017[math]\pm[/math]0.00002113646
27109 Missing mass difference RunNumber27109 1 OSICuts.gif FC difference RunNumber27109 OSICuts.gif -0.00876887465[math]\pm[/math]0.00002411331
27107 Missing mass difference RunNumber27107 1 OSICuts.gif FC difference RunNumber27107 OSICuts.gif 0.01044530752[math]\pm[/math]0.00002188631
27116 Missing mass difference RunNumber27116 1 OSICuts.gif FC difference RunNumber27116 OSICuts.gif 0.00173242954[math]\pm[/math]0.00003084415
27112 Missing mass difference RunNumber27112 1 OSICuts.gif FC difference RunNumber27112 OSICuts.gif -0.01142682761[math]\pm[/math]0.00002505827
27111 Missing mass difference RunNumber27111 1 OSICuts.gif FC difference RunNumber27111 OSICuts.gif -0.00872760886[math]\pm[/math]0.00002117268
27128 Missing mass difference RunNumber27128 1 OSICuts.gif FC difference RunNumber27128 OSICuts.gif 0.00193441756[math]\pm[/math]0.00002177434
27127 Missing mass difference RunNumber27127 1 OSICuts.gif FC difference RunNumber27127 OSICuts.gif 0.00215712098[math]\pm[/math]0.00002186032
27124 Missing mass difference RunNumber27124 1 OSICuts.gif FC difference RunNumber27124 OSICuts.gif 0.00309070808[math]\pm[/math]0.00002235601
27139 Missing mass difference RunNumber27139 1 OSICuts.gif FC difference RunNumber27139 OSICuts.gif -0.00288862444[math]\pm[/math]0.00002683943
27138 Missing mass difference RunNumber27138 1 OSICuts.gif FC difference RunNumber27138 OSICuts.gif -0.00061608634[math]\pm[/math]0.00002396188
27137 Missing mass difference RunNumber27137 1 OSICuts.gif FC difference RunNumber27137 OSICuts.gif -0.00275066778[math]\pm[/math]0.00002150471
27136 Missing mass difference RunNumber27136 1 OSICuts.gif FC difference RunNumber27136 OSICuts.gif -0.00239423237[math]\pm[/math]0.00002151354
27134 Missing mass difference RunNumber27134 1 OSICuts.gif FC difference RunNumber27134 OSICuts.gif -0.00355581434[math]\pm[/math]0.00002051114
27133 Missing mass difference RunNumber27133 1 OSICuts.gif FC difference RunNumber27133 OSICuts.gif -0.00252193724[math]\pm[/math]0.00002181781
27132 Missing mass difference RunNumber27132 1 OSICuts.gif FC difference RunNumber27132 OSICuts.gif 0.00096332202[math]\pm[/math]0.00002224348
27143 Missing mass difference RunNumber27143 1 OSICuts.gif FC difference RunNumber27143 OSICuts.gif -0.00210486278[math]\pm[/math]0.00002838738
27141 Missing mass difference RunNumber27141 1 OSICuts.gif FC difference RunNumber27141 OSICuts.gif -0.00301990556[math]\pm[/math]0.00002515365
27160 Missing mass difference RunNumber27160 1 OSICuts.gif FC difference RunNumber27160 OSICuts.gif 0.00171264744[math]\pm[/math]0.00002170610
27161 Missing mass difference RunNumber27161 1 OSICuts.gif FC difference RunNumber27161 OSICuts.gif 0.00224287304[math]\pm[/math]0.00002374488
27162 Missing mass difference RunNumber27162 1 OSICuts.gif FC difference RunNumber27162 OSICuts.gif 0.00118844445[math]\pm[/math]0.00002480837
27166 Missing mass difference RunNumber27166 1 OSICuts.gif FC difference RunNumber27166 OSICuts.gif 0.00077237006[math]\pm[/math]0.00002247548
27167 Missing mass difference RunNumber27167 1 OSICuts.gif FC difference RunNumber27167 OSICuts.gif 0.00216041281[math]\pm[/math]0.00002188669
27168 Missing mass difference RunNumber27168 1 OSICuts.gif FC difference RunNumber27168 OSICuts.gif 0.00083650625[math]\pm[/math]0.00002173328
27170 Missing mass difference RunNumber27170 1 OSICuts.gif FC difference RunNumber27170 OSICuts.gif 0.00151247511[math]\pm[/math]0.00002430185
27175 Missing mass difference RunNumber27175 1 OSICuts.gif FC difference RunNumber27175 OSICuts.gif -0.01217485668[math]\pm[/math]0.00002235229
27176 Missing mass difference RunNumber27176 1 OSICuts.gif FC difference RunNumber27176 OSICuts.gif -0.01049635167[math]\pm[/math]0.00002127680
27177 Missing mass difference RunNumber27177 1 OSICuts.gif FC difference RunNumber27177 OSICuts.gif -0.01147445288[math]\pm[/math]0.00002149007
27179 Missing mass difference RunNumber27179 1 OSICuts.gif FC difference RunNumber27179 OSICuts.gif -0.01025340153[math]\pm[/math]0.00002206557
27180 Missing mass difference RunNumber27180 1 OSICuts.gif FC difference RunNumber27180 OSICuts.gif -0.00868272418[math]\pm[/math]0.00002169057
27181 Missing mass difference RunNumber27181 1 OSICuts.gif FC difference RunNumber27181 OSICuts.gif -0.00973324188[math]\pm[/math]0.00002374079
27182 Missing mass difference RunNumber27182 1 OSICuts.gif FC difference RunNumber27182 OSICuts.gif -0.00983862460[math]\pm[/math]0.00002233926
27183 Missing mass difference RunNumber27183 1 OSICuts.gif FC difference RunNumber27183 OSICuts.gif -0.01901615252[math]\pm[/math]0.00002135060
27186 Missing mass difference RunNumber27186 1 OSICuts.gif FC difference RunNumber27186 OSICuts.gif 0.01079090635[math]\pm[/math]0.00002093983
27187 Missing mass difference RunNumber27187 1 OSICuts.gif FC difference RunNumber27187 OSICuts.gif 0.01155519076[math]\pm[/math]0.00002157075
27188 Missing mass difference RunNumber27188 1 OSICuts.gif FC difference RunNumber27188 OSICuts.gif 0.01086898605[math]\pm[/math]0.00002089859
27190 Missing mass difference RunNumber27190 1 OSICuts.gif FC difference RunNumber27190 OSICuts.gif 0.01181266490[math]\pm[/math]0.00002159430
27192 Missing mass difference RunNumber27192 1 OSICuts.gif FC difference RunNumber27192 OSICuts.gif -0.00988459637[math]\pm[/math]0.00002151795
27193 Missing mass difference RunNumber27193 1 OSICuts.gif FC difference RunNumber27193 OSICuts.gif -0.01047643036[math]\pm[/math]0.00002159742
27194 Missing mass difference RunNumber27194 1 OSICuts.gif FC difference RunNumber27194 OSICuts.gif -0.01185520003[math]\pm[/math]0.00002134810


RunNumber vs FCAsymmetry EndofRunSum.jpg

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
The number of photoelectrons without cuts
The number of photoelectrons with OSI cuts


For different fits:

The number of photoelectrons after OSICuts with two Gaussian fits
The number of photoelectrons after OSICuts with Landau+Gaussian 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)
The number of photoelectrons without cuts
The number of photoelectrons with OSI cuts(gauss+landau+gauss)
The number of photoelectrons with OSI cuts two gaussian fits
The number of photoelectrons with OSI+NPHE cuts


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

Electrons nphe with OSIcuts all data Gauss0.gifElectrons nphe with OSIcuts all data Landau3.gifElectrons nphe with OSIcuts all data Gauss6.gif

[math] Pion Contamination = [/math]
[math] = \frac {Integral(landau(3))}{Integral(gauss(0) + landau(3) + gauss(6))} =[/math]
[math] = \frac{1.643 \times 10^9}{ 7.682 \times 10^9 + 1.643 \times 10^9 + 7.732 \times 10^9} = [/math]
[math] = 9.6324 % [/math]

OSICut+NPHE>2.5

Electrons nphe with OSIcuts all data Gauss0nphe2-5.gifElectrons nphe with OSIcuts all data Landau3nphe2-5.gifElectrons nphe with OSIcuts all data Gauss6nphe2-5.gif


[math] Pion Contamination = [/math]
[math] = \frac {Integral(landau(3))}{Integral(gauss(0) + landau(3) + gauss(6))} =[/math]
[math] = \frac{5.905 \times 10^8}{ 6.656 \times 10^9 + 5.905 \times 10^8 + 7.395 \times 10^9} = [/math]
[math] = 4.03305% [/math]

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

Electrons nphe with OSIcuts all data twogaussians Gauss0.gifElectrons nphe with OSIcuts all data twogaussians Gauss3.gif

[math] Pion Contamination = [/math]
[math] = \frac {Integral(gauss(0))}{Integral(gauss(0) + gauss(3))} =[/math]
[math] = \frac{2.151 \times 10^8}{ 2.151 \times 10^8 + 1.594 \times 10^{10} } [/math]
[math] = 1.3 % [/math]

With NPHE>2.5 Cut

[math] Pion Contamination = 0 [/math]

Number of Events after NPHE>2.5 Cut

[math] Number of Events after NPHE\gt 2.5 Cut [/math] =
[math] = \frac{2.978 \times 10^8}{3.496 \times 10^8} = [/math]
[math] = 85.18 % [/math]

Counts in FCup

FCupCountsForHelicity+ 29679 13 OSICuts Alldata.gifFCupCountsForHelicity- 29679 24 OSICuts Alldata.gif

[math]Ratio of Counts 24/13 = \frac{2.131 \times 10^6}{2.116 \times 10^6} = 1.007 [/math]

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 [math]\equiv \sum_i(\frac{FC(i)^+ - FC(i)^-}{FC(i)^+ + FC(i)^-})[/math]

End of Run sum [math]\equiv \frac {\sum_i(FC(i)^+) - \sum(FC(i)^-)}{\sum (FC(i)^+) + \sum(FC(i)^-)})[/math]

4.) Determine semi-inclusive statistic as function of X

Uncertainty in Pion Contamination

Maximum

Electrons nphe with OSIcuts all data Gauss0MAX.gifElectrons nphe with OSIcuts all data Landau3MAX.gifElectrons nphe with OSIcuts all data Gauss6MAX.gif

[math] Pion Contamination = [/math]
[math] = \frac {Integral(landau(3))}{Integral(gauss(0) + landau(3) + gauss(6))} =[/math]
[math] = \frac{1.645 \times 10^9}{ 7.695 \times 10^9 + 1.645 \times 10^9 + 7.745 \times 10^9} = [/math]
[math] = 9.6283 % [/math]

Electrons nphe with OSIcuts all data Gauss0MAX nphe2-5cut.gifElectrons nphe with OSIcuts all data Landau3MAX nphe2-5cut.gifElectrons nphe with OSIcuts all data Gauss6MAX nphe2-5cut.gif

[math] Pion Contamination = [/math]
[math] = \frac {Integral(landau(3))}{Integral(gauss(0) + landau(3) + gauss(6))} =[/math]
[math] = \frac{5.914 \times 10^8}{ 6.666 \times 10^9 + 5.914 \times 10^8 + 7.427 \times 10^9} = [/math]
[math] = 4.02740 % [/math]

Minimum

Electrons nphe with OSIcuts all data Gauss0MIN.gifElectrons nphe with OSIcuts all data Landau3MIN.gifElectrons nphe with OSIcuts all data Gauss6MIN.gif

[math] Pion Contamination = [/math]
[math] = \frac {Integral(landau(3))}{Integral(gauss(0) + landau(3) + gauss(6))} =[/math]
[math] = \frac{1.642 \times 10^9}{ 7.67 \times 10^9 + 1.642 \times 10^9 + 7.719 \times 10^9} = [/math]
[math] = 9.64124 % [/math]

Electrons nphe with OSIcuts all data Gauss0MIN nphe2-5cut.gifElectrons nphe with OSIcuts all data Landau3MIN nphe2-5cut.gifElectrons nphe with OSIcuts all data Gauss6MIN nphe2-5cut.gif

[math] Pion Contamination = [/math]
[math] = \frac {Integral(landau(3))}{Integral(gauss(0) + landau(3) + gauss(6))} =[/math]
[math] = \frac{5.895 \times 10^8}{ 6.645 \times 10^9 + 5.895 \times 10^8 + 7.401 \times 10^9} = [/math]
[math] = 4.02788 % [/math]

Pion Contamination

It appears that pion contamination in electron sample is 9.63 % [math]\pm[/math] 0.01 % before nphe cut and after nphe>2.5 cut contamination is about 4.029% [math]\pm[/math] 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.

X bjorken withoutcuts electronpionrequired.gifX bjorken OSICuts electronpionrequired.gifX bjorken OSINPHECuts electronpionrequired.gif
Qsqrd withoutcuts electronpionrequired.gifQsqrd OSICuts electronpionrequired.gifQsqrd OSINPHECuts electronpionrequired.gif


Number of Events after cuts
No Cuts OSI Cuts OSI+NPHE>2.5 Cuts
[math]6.724 \times 10^7[/math] [math]4.606 \times 10^7[/math] [math]3.868 \times 10^7[/math]
68.5 % 57.5 %


Error Calculation

Tthe error in the asymmetry measurement would be [math]\frac{\Delta A}{A} = \frac{2}{\sqrt{N}}[/math]

X_b [math] Number(X_b)^+[/math] [math] Number(X_b)^-[/math] X_b Asymmetry Error
0.1 [math]5.25 \times 10^6[/math] [math]5.25 \times 10^6[/math] [math]-4.001 \times 10^{-4}[/math] 0.00087251693
0.2 [math]5.52 \times 10^6[/math] [math]5.3 \times 10^6[/math] [math]-7.86 \times 10^{-4}[/math] 0.0008507530911145
0.3 [math]3.496 \times 10^6[/math] [math]3.50 \times 10^6[/math] [math]-9.025 \times 10^{-4}[/math] 1.0691459e-03
0.4 [math]2.0379 \times 10^6[/math] [math]2.04 \times 10^6[/math] [math]-8.37 \times 10^{-4}[/math] 0.0014004231
0.5 [math]1.14 \times 10^6[/math] [math]1.15 \times 10^6[/math] [math]-2.978 \times 10^{-4}[/math] 0.0018665742
0.6 [math]6.37\times 10^5[/math] [math]6.38 \times 10^5[/math] [math]-1.115 \times 10^{-3}[/math] 0.0016477095
0.7 [math]3.514 \times 10^5[/math] [math]3.519 \times 10^5[/math] [math]-8.26 \times 10^{-3}[/math] 0.0022190018
0.8 [math]2.022 \times 10^5[/math] [math] 2.042 \times 10^5[/math] [math]-5.039 \times 10^{-3}[/math] 0.00291609
0.9 [math]1.348 \times 10^5[/math] [math]1.34 \times 10^5[/math] [math]2.935 \times 10^{-3}[/math] 0.003592967
1 [math]1.038 \times 10^5[/math] [math]1.033 \times 10^5[/math] [math]2.47 \times 10^{-3}[/math] 0.0040928449

X b vs Asymmetry OSICUTs+NPHE2-5 1.jpg


I am pretty sure X_{BJ} > 0.8 is not possible with our data set

Electron theta angle and [math]Q^2[/math] cuts

electron theta angle for different X_b

ElectronThetaAgle less0-8X b.gifElectronThetaAgle above0-8X b.gif

[math]Q^2[/math]

Q^2 less0-8X b.gifQ^2 above0-8X b.gif

X_b when [math]Q^2 \lt 1[/math]

X b forQ2less1 onefile.gif


X b NumberofEventsAbove0-8X b.gif

Number of Events for X_b>0.8 [math] = \frac{3531}{313537} = 1.1 %[/math]

plot the vertex of the above hits with X>0.8

X b OSI+EC+NPHE allcuts.gif

change below to log plots so we can see where XBj stops

X b OSI+EC+NPHE allcuts and withoutcuts.gifX b OSI+EC+NPHE allcuts and withoutcuts LogScale.gif

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 [math]Q^2[/math] Vertex X Vertex Y Vertex Z
OSI Cuts + EC Cuts X b Run26994 OSI+EC allx.gif Q 2 Run26994 OSI+EC allx.gif Vertex x Run26994 OSI+EC allx.gif Vertex y Run26994 OSI+EC allx.gif Vertex z Run26994 OSI+EC allx.gif
OSI Cuts + EC Cuts + X_b>0.8 X bgreater0-8 Run26994 OSI+EC allx.gif Q 2greater0-8 Run26994 OSI+EC allx.gif Vertex xgreater0-8 Run26994 OSI+EC allx.gif Vertex ygreater0-8 Run26994 OSI+EC allx.gif Vertex zgreater0-8 Run26994 OSI+EC allx.gif
OSI Cuts + EC Cuts + X_b<0.8 X bless0-8 Run26994 OSI+EC allx.gif Q 2less0-8 Run26994 OSI+EC allx.gif Vertex xless0-8 Run26994 OSI+EC allx.gif Vertex yless0-8 Run26994 OSI+EC allx.gif Vertex zless0-8 Run26994 OSI+EC allx.gif
OSI Cuts + EC Cuts + 1<Q^2<4 X bQcut Run26994 OSI+EC allx.gif Q 2Qcut Run26994 OSI+EC allx.gif Vertex xQcut Run26994 OSI+EC allx.gif Vertex yQcut Run26994 OSI+EC allx.gif Vertex zQcut Run26994 OSI+EC allx.gif


Cuts The scattered electron energy [math]\theta[/math] electron scattering angle
[math]X_b \gt 0.8[/math] Scattered electron energy 1.gif Theta electron scattering angle 1.gif
[math]X_b \lt 0.8[/math] Scattered electron energy 2.gif Theta electron scattering angle 2.gif

FC Asymmetry

FC Asymmetry plot using the following method : End of Run sum [math]\equiv \frac {\sum_i(FC(i)^+) - \sum(FC(i)^-)}{\sum (FC(i)^+) + \sum(FC(i)^-)})[/math]

RunNumber vs FCAsymmetry EndofRunSum.jpg


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 % [math]\pm[/math] 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)

Ectotpvsnphebefore.gifEctotpvsnpheafter.gif


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. : [math]9.446 % \pm 0.233 %[/math] && [math]3.76 % \pm 0.08 %[/math]

3.)

[math]X_b\gt 0.8[/math]

 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

[math]X_b\lt 0.8[/math]

  • 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.)

[math]\pi^+[/math] && [math]B^+[/math] [math]\pi^-[/math] && [math]B^+[/math] [math]\pi^-[/math] && [math]B^-[/math] [math]\pi^+[/math] && [math]B^-[/math]
Positivetorusmagnet positivepion.gif Positivetorusmagnet negativepion.gif Negativetorusmagnet negativepion.gif Negativetorusmagnet positivepion.gif

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 % [math]\pm[/math] 0.01 %
3 S. D. 9.645 % [math]\pm[/math] 0.025 %
10 S. D. 9.446 % [math]\pm[/math] 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

Alldatasector27113 5 1 root.gif

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


Alldatasector27113 5 root.gif

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



[math]x_b[/math] 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 [math]\pi^+[/math] [math]\pi^-[/math]
26990, B<0 Electron26990 4.gif Electron26990 5.gif 4pion26990.gif 5pion26990.gif
27113, B>0 Electron27113 4.gif Electron27113 5.gif 4pion27113.gif 5pion27113.gif
  • 1.) 7< Sector_paddle <11 - (B>0, [math]\pi^-[/math]) - 15.3% && (B<0, [math]\pi^+[/math]) - 10.9%
  • 2.) 25< Sector_paddle <29 - (B>0, [math]\pi^+[/math]) - 7.775% && (B<0, [math]\pi^-[/math]) - 10.97%

[math]\frac{\Delta d}{d}[/math] && [math]\frac{\Delta d}{d}[/math]

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
[math]e^-[/math] && [math]\pi^+[/math] PositivepionpaddlenumberNH3.gif X b vs positivepionpaddlenumberNH3.gif 7
[math]e^-[/math] && [math]\pi^-[/math] NegativepionpaddlenumberNH3.gif X b vs negativepionpaddlenumberNH3.gif 27


Detected particles in the final state and chosen pion paddle number [math]e^-[/math] paddle number [math]X_b[/math] vs [math]e^-[/math] paddle number [math]X_b[/math] vs [math]e^-[/math] paddle number when [math]x_b\lt 0.3[/math] [math]X_b[/math] vs [math]e^-[/math] paddle number when [math]x_b\gt 0.3[/math]
[math]e^-[/math] && [math]\pi^+[/math], [math]PaddleNumber_{\pi^+}=7[/math] Electronpaddlenumberforpositivepionpaddlenumber7.gif X b vs electronpaddlenumberforpositivepionpaddlenumber7.gif X b vs electronpaddlenumberforpositivepionpaddlenumber7lowX b.gif X b vs electronpaddlenumberforpositivepionpaddlenumber7highX b.gif
[math]e^-[/math] && [math]\pi^-[/math], [math]PaddleNumber_{\pi^-}=27[/math] Electronpaddlenumberfornegativepionpaddlenumber27.gif X b vs electronpaddlenumberfornegativepionpaddlenumber27.gif X b vs electronpaddlenumberfornegativepionpaddlenumber27lowX b.gif X b vs electronpaddlenumberfornegativepionpaddlenumber27highX b.gif

Positive Torus

Using runs with NH3 target

Detected particles in the final state X_b vs pion paddle number Chosen pion paddle number
[math]e^-[/math] && [math]\pi^+[/math] X b vs positivepionpaddlenumberNH3positivetorus.gif 27
[math]e^-[/math] && [math]\pi^-[/math] X b vs negativepionpaddlenumberNH3positivetorus.gif 7


Detected particles in the final state and chosen pion paddle number [math]e^-[/math] paddle number [math]X_b[/math] vs [math]e^-[/math] paddle number [math]X_b[/math] vs [math]e^-[/math] paddle number when [math]x_b\lt 0.3[/math] [math]X_b[/math] vs [math]e^-[/math] paddle number when [math]x_b\gt 0.3[/math]
[math]e^-[/math] && [math]\pi^+[/math], [math]PaddleNumber_{\pi^+}=27[/math] Electronpaddlenumberforpositivepionpaddlenumber7positivetorus.gif X b vs electronpaddlenumberforpositivepionpaddlenumber27positivetorus.gif X b vs electronpaddlenumberforpositivepionpaddlenumber27lowX bpositivetorus.gif X b vs electronpaddlenumberforpositivepionpaddlenumber27highX bpositivetorus.gif
[math]e^-[/math] && [math]\pi^-[/math], [math]PaddleNumber_{\pi^-}=7[/math] Electronpaddlenumberfornegativepionpaddlenumber27positivetorus.gif X b vs electronpaddlenumberfornegativepionpaddlenumber7positivetorus.gif X b vs electronpaddlenumberfornegativepionpaddlenumber7lowX bpositivetorus.gif X b vs electronpaddlenumberfornegativepionpaddlenumber7highX bpositivetorus.gif