Difference between revisions of "Tamar Electron and Pion Efficiency"

From New IAC Wiki
Jump to navigation Jump to search
 
(15 intermediate revisions by 2 users not shown)
Line 42: Line 42:
 
The EG1b data for kinematics chosen above show the same shape as E99-107 data.
 
The EG1b data for kinematics chosen above show the same shape as E99-107 data.
  
=Inclusive  and Semi-Inclusive Efficiencies=
+
=Inclusive  and Exclusive Efficiencies=
  
 
The goal of this thesis is to measure the semi-inclusive asymmetry when an electron and a pion are detected in the final state.  Pions of opposite charge will be observed using the same scintillator by flipping the CLAS Torus magnetic field direction. Although the pions will be detected by the same detector elements, the electrons will intersect different detector elements.  As a result, the electron efficiency will need to be evaluated in terms of the electron rate observed in two different scintillator paddles detecting the same electron kinematics.  Pairs of scintillator paddles  with the highest semi-inclusive rate.
 
The goal of this thesis is to measure the semi-inclusive asymmetry when an electron and a pion are detected in the final state.  Pions of opposite charge will be observed using the same scintillator by flipping the CLAS Torus magnetic field direction. Although the pions will be detected by the same detector elements, the electrons will intersect different detector elements.  As a result, the electron efficiency will need to be evaluated in terms of the electron rate observed in two different scintillator paddles detecting the same electron kinematics.  Pairs of scintillator paddles  with the highest semi-inclusive rate.
  
The goal of this paper is to compare measured single pion electroproduction efficiency to the unitary model MAID 2007 that has been developed to analyze the world data of pion photo and electroproduction. The model is well adopted for predictions of the observables for pion production, like five fold cross section, total cross secton and etc.
 
  
The MAID 2007 model has predictions of the total cross section for the following two cases, that are related to our work:<br>
+
The electron efficiency of individual scintillator detectors using the 4.2 GeV data for ND3 and NH3 targets is investigated below. Only the electron is detection in the final state (inclusive case). The pion contamination in the electron sample was removed by the applying cuts described above.  The electron paddle number 10 (B<0)  and 5 (B>0) were chosen respectively because they contained the most electron events in a first pass semi-inclusive pion analysis of the data set. The electron kinematics(Momentum, scattering angle and invariant mass) for these scintillators is shown on Fig. 3.1.
 
 
:::::::::::<math>\gamma</math> + proton<math>(NH_3) \rightarrow</math> <math>\pi^+</math> + neutron <br>
 
 
 
:::::::::::<math>\gamma</math> + neutron<math>(ND_3) \rightarrow</math> <math>\pi^-</math> + proton <br>
 
  
The electron efficiency of individual scintillator detectors using the 4.2 GeV data for ND3 and NH3 targets is described below. In the final state, only electron detection was required(inclusive case). The contamination in the electron sample was removed by applying cuts described in the first chapter.  The electron paddle number 10 (B<0)  and 5 (B>0) were chosen  respectively because they contained the most electron events in a first pass semi-inclusive pion analysis of the data set. The electron kinematics(Momentum, scattering angle and invariant mass) for these scintillators is shown on Fig. 3.1.
+
FC is not a unit, you need to use the calibration to convert FC counts to nC
  
 
{| border="1"  |cellpadding="20" cellspacing="0  
 
{| border="1"  |cellpadding="20" cellspacing="0  
 
|-
 
|-
|[[File:EmomInclusiveoverlay4-2GeV.gif|300px]] ||[[File:EthetaInclusiveoverlay4-2GeV.gif|300px]] || [[File:WInclusiveoverlay4-2GeV.gif|300px]]
+
|[[File:EmomInclusiveoverlay4-2GeVfcup.gif|300px]] ||[[File:EthetaInclusiveoverlay4-2GeVfcup.gif|300px]] || [[File:WInclusiveoverlay4-2GeVfcup.gif|300px]]
 
|-
 
|-
 
| Electron Momentum((NH3,B>0),  (NH3,B<0), (ND3,B>0) && (ND3,B<0))|| Electron Scattering Angle <math>\theta</math>((NH3,B>0),  (NH3,B<0), (ND3,B>0) && (ND3,B<0))|| W Invariant mass((NH3,B>0),  (NH3,B<0), (ND3,B>0) && (ND3,B<0))
 
| Electron Momentum((NH3,B>0),  (NH3,B<0), (ND3,B>0) && (ND3,B<0))|| Electron Scattering Angle <math>\theta</math>((NH3,B>0),  (NH3,B<0), (ND3,B>0) && (ND3,B<0))|| W Invariant mass((NH3,B>0),  (NH3,B<0), (ND3,B>0) && (ND3,B<0))
Line 66: Line 61:
  
  
To be able to compare our results to the MAID 2007, the ratio of the number of electrons were taken for different types of target for inbending and outbending cases. The ratios are following and the are normalized by faraday cup count:
+
Ratios of the inclusive electron rate, normalized using the FC, of scintillator paddles 5 and 10 were measured.   The two ratios are designed to measure the CLAS detectors ability to observed the same electron kinematics using different detector elements positioned for opposite Torus polarities. 
 +
 
 +
<math>\frac{ND3,B>0, PaddleNumber^{e^-}=5}{ND3,B<0, PaddleNumber^{e^-}=10}=1.57 \pm 0.16</math>
 +
 
 +
<math>\frac{NH3,B>0, PaddleNumber^{e^-}=5}{NH3,B<0, PaddleNumber^{e^-}=10}=1.76 \pm 0.17</math>
 +
 
 +
 
 +
Notice the ratio is statistically the same if the NH3 and ND3 targets are used for this ratio in a manner which will be similar to the semi-inclusive analysis.
  
 
<math>\frac{ND3,B>0, PaddleNumber^{e^-}=5}{NH3,B<0, PaddleNumber^{e^-}=10}=1.55 \pm 0.15</math>
 
<math>\frac{ND3,B>0, PaddleNumber^{e^-}=5}{NH3,B<0, PaddleNumber^{e^-}=10}=1.55 \pm 0.15</math>
  
<math>\frac{ND3,B<0, PaddleNumber^{e^-}=10}{NH3,B>0, PaddleNumber^{e^-}=5}=0.55 \pm 0.06</math>
 
  
In order to make detector efficiency the same for electrons, the ratios need to be adjusted so that it equals to one. Each of the two cases are multiplied by so called "correction coefficient".
+
 
 +
The above ratios, which have been observed to be ammonia target independent, indicate a difference in electron detector efficiency when the Torus polarity is flipped.  In order to make detector efficiency the same for electrons, the ratios were used as a "correction coefficient".   The "correction coefficient" for the case <math>\frac{ND3,B>0,E_{PaddleNumber}=5}{NH3,B<0,E_{PaddleNumber}=10}</math> is <math>0.645</math> and for the <math>\frac{ND3,B<0,E_{PaddleNumber}=10}{NH3,B>0,E_{PaddleNumber}=5}</math> it is <math>1.82</math>.
  
  
In this chapter we discuss the pion efficiency in semi-inclusive case. The ratio of the pions detected in the scintillator paddles, located between the Cherenkov counter and electromagnetic calorimeter, is shown on Fig. 3.2. The ratios were taken for four different cases. Assuming that, for the inbending case positive pions and for the outbending case negative pions have the same trajectories(the same kinematics) and vice versa((the inbending,negative pion) and (the outbending, positive pions)).<br>
+
 
 +
 
 +
The measured single pion electroproduction rate was compared to the MAID 2007 unitary model that has been developed using the world data of pion photo and electro-production to determine the impact of using the above "correction coefficient". The model is well adopted for predictions of the observables for pion production, like five fold cross section, total cross secton and etc.
 +
 
 +
The MAID 2007 model has predictions of the total cross section for the following two cases, that are related to our work:<br>
 +
 
 +
:::::::::::<math>\gamma</math> + proton<math>(NH_3) \rightarrow</math> <math>\pi^+</math> + neutron <br>
 +
 
 +
:::::::::::<math>\gamma</math> + neutron<math>(ND_3) \rightarrow</math> <math>\pi^-</math> + proton <br>
 +
 
 +
 
 +
 
 +
The ratio of the pions detected in the scintillator paddles, located between the Cherenkov counter and electromagnetic calorimeter, is shown on Fig. 3.2. The ratios were taken for four different cases. Assuming that, for the inbending case positive pions and for the outbending case negative pions have the same trajectories(the same kinematics) and vice versa((the inbending,negative pion) and (the outbending, positive pions)).<br>
  
 
[[File:paddenumbvsratio.jpg|500px]]<br>
 
[[File:paddenumbvsratio.jpg|500px]]<br>
Line 86: Line 100:
 
After applying correction coefficients from inclusive cases, the ratios have been compared to the results from MAID2007.
 
After applying correction coefficients from inclusive cases, the ratios have been compared to the results from MAID2007.
  
[[File:black_red_maid_inclusiveerror.png|350px]] [[File:green_blue_inclusiveerror.png|350px]]
+
[[File:negativepionND3beforecorrection.png|350px]] [[File:negativepionND3aftercorrection.png|350px]]
 +
 
 +
[[File:positivepionND3beforecorrection.png|350px]] [[File:positivepionND3aftercorrection.png|350px]]
 +
 
  
'''Figure 3.3.  Pion paddle number vs Ratio after correction.'''
+
'''Figure 3.3.  Pion paddle number vs Ratio before and after correction.'''
  
 
Applied corrections are following:
 
Applied corrections are following:

Latest revision as of 15:01, 10 August 2011

Electron-Pion separation

The requirements of the measured particles attributes that are used to differentiate electrons from negative pions is described below. Separating a negatively charged pion from an electron becomes more difficult when a negative pion has sufficient energy to generate cerenkov radiation and as a result may be misidentified as an electron. This contamination of the electron sample is reduced using cuts on the energy deposited in the electromagnetic calorimeter and the number of photoelectrons in the Cherenkov counter.

The energy deposited in the calorimeter for electrons and pions is different. Pions are minimum ionizing charged particles when their momentum is above ~ 0.8 GeV and therefore the energy they deposit in the calorimeter is momentum independent. On the other hand, electrons produce photoelectrons and generate an electromagnetic shower releasing energy into the calorimeter that is proportional to their momentum. The following cut is introduced to take advantage of this feature: [math]EC_{total}\gt 0.2*p[/math]. The cut was also applied to the energy collected in the inner part of the calorimeter: [math]EC_{inner}\gt 0.06*p[/math] , because the ratio of the total energy deposited to the energy deposited in the inner calorimeter depends on the thickness of the detector and is a constant.
E total vs e inner1 before cuts file dst27070.gifE total vs e inner1 after cuts file dst27070.gif

Figure 1.1. The energy deposition into the total and inner calorimeter before and after cuts.

In addition to the cut on the energy deposited into the electromagnetic calorimeter, misidentified electrons are reduced by requiring a signal in the threshold CLAS Cherenkov detector. Pion's misidentified as electrons have been shown to produce around ~1.5 photoelectons (PEs) in the Cherenkov detector, as shown below. Geometrical cuts on the location of the particle at the entrance to the Cerenkov detector were applied to reduce the pion contamination. The second histogram below shows that, after these cuts, the peak around 1.5 PEs is substantially reduced.

No cuts OSI Cuts
The number of photoelectrons without cuts
The number of photoelectrons with OSI cuts

Figure 1.2. The number of photoelectrons before and after cuts.

Exclusive Pi+ production using NH3

A comparison was made with the results from experiment E99-107 <ref name="Park2008" >Park, K., Burkert, V. D., & Kim, W. (The CLAS Collaboration) Phys. Rev., C77, 015208.</ref> as one benchmark to evaluate the integrity of this semi-inclusive analysis. Experiment E99-107 measured exclusive pion production ([math]\vec{e} p[/math]->[math]en\pi^+[/math]) using a liquid hydrogen target and the CLAS. The EG1 experiment performed similar measurements using a polarized NH3 target allowing for an asymmetry measurement in addition to the differential cross-section in the nucleon resonance region for [math]1.7 \lt Q^2 \lt 4.5[/math] [math](GeV)^2[/math]. In the EG1 experiment polarized target and a 5.7 GeV continuous electron beam resulted in a luminosity of [math] 1.09209 \times 10^{33} \frac{[ \# of electrons]}{second \times cm^2}[/math].

The kinematics of single pion electroproduction can be described by five variables: the virtual photon negative four-momentum transfered squared [math]Q^2[/math], W invariant mass of the photon-nucleon system, [math]{\theta_{\pi}}^*[/math] the polar and [math]{\varphi_{\pi}}^*[/math] the azimuthal angle of the outgoing pion in center of mass frame and [math]{\varphi}_e[/math] the scattered electron azimuthal angle.

Figure 2.1. Kinematics of single [math]\pi^+[/math] electroproduction


The five-fold differential cross section can be written in the following way for a single pion electroproduction:

[math]\frac{\partial^5 \sigma}{\partial E_f \partial \Omega_e \partial {\Omega_{\pi}}^*} = \frac{1}{2 \pi} \Sigma \frac{1}{L_{int} A_{cc} \epsilon_{CC} \Delta W \Delta Q^2 \Delta cos {\theta_{\pi}}^{\!*} \Delta {\phi_{\pi}}^{\!*}} \frac{d(W, Q^2)}{d(E_f, cos \theta_e)}[/math]

where [math]L_{int}[/math] - the integrated luminosity, [math]A_{cc}[/math] is the acceptance factor, [math]\epsilon_{CC} [/math] represents the efficiency of the cherenkov detector and the Jacobian term can be expressed in terms of the initial and final energy of lepton:

[math]\frac{d(W, Q^2)}{d(E_f, cos \theta_e)} = \frac{2 M_p E_i E_f}{W}[/math]

Because of wide range of kinematics, [math]{\varphi_{\pi}}^*[/math] was measured only for certain lepton scattering [math]\theta_e[/math] angles and invariant mass (W). Applying above described cuts: EC_inner>0.06, EC_tot/p>0.2, nphe>2.5 and [math]0.9\lt M_x\lt 1.1[/math], for the following invariant mass [math]1.44\lt W\lt 1.46[/math] and [math]0.4\lt cos\theta_{pion}^{CM}\lt 0.6[/math] the [math]{\varphi_{\pi}}^*[/math] vs relative rate distribution is shown below on the graph and compared with E99-107 data, which by itself is in agreement with the models.

Figure 2.2. [math]{\varphi_{\pi}}^*[/math] vs Relative rate for fixed [math]cos\theta_{pion}^{CM}=0.5[/math] and [math]W=1.45 GeV[/math] <ref name="Park2008" > Bogus text</ref> <ref> http://clasweb.jlab.org/cgi-bin/clasdb/msm.cgi?eid=14&mid=16&data=on </ref>


The EG1b data for kinematics chosen above show the same shape as E99-107 data.

Inclusive and Exclusive Efficiencies

The goal of this thesis is to measure the semi-inclusive asymmetry when an electron and a pion are detected in the final state. Pions of opposite charge will be observed using the same scintillator by flipping the CLAS Torus magnetic field direction. Although the pions will be detected by the same detector elements, the electrons will intersect different detector elements. As a result, the electron efficiency will need to be evaluated in terms of the electron rate observed in two different scintillator paddles detecting the same electron kinematics. Pairs of scintillator paddles with the highest semi-inclusive rate.


The electron efficiency of individual scintillator detectors using the 4.2 GeV data for ND3 and NH3 targets is investigated below. Only the electron is detection in the final state (inclusive case). The pion contamination in the electron sample was removed by the applying cuts described above. The electron paddle number 10 (B<0) and 5 (B>0) were chosen respectively because they contained the most electron events in a first pass semi-inclusive pion analysis of the data set. The electron kinematics(Momentum, scattering angle and invariant mass) for these scintillators is shown on Fig. 3.1.

FC is not a unit, you need to use the calibration to convert FC counts to nC
EmomInclusiveoverlay4-2GeVfcup.gif EthetaInclusiveoverlay4-2GeVfcup.gif WInclusiveoverlay4-2GeVfcup.gif
Electron Momentum((NH3,B>0), (NH3,B<0), (ND3,B>0) && (ND3,B<0)) Electron Scattering Angle [math]\theta[/math]((NH3,B>0), (NH3,B<0), (ND3,B>0) && (ND3,B<0)) W Invariant mass((NH3,B>0), (NH3,B<0), (ND3,B>0) && (ND3,B<0))

Figure 3.1. Electron Kinematics(The histograms presented above are not normalized by faraday cup. We just want to show the kinematic range for chosen electron events).


Ratios of the inclusive electron rate, normalized using the FC, of scintillator paddles 5 and 10 were measured. The two ratios are designed to measure the CLAS detectors ability to observed the same electron kinematics using different detector elements positioned for opposite Torus polarities.

[math]\frac{ND3,B\gt 0, PaddleNumber^{e^-}=5}{ND3,B\lt 0, PaddleNumber^{e^-}=10}=1.57 \pm 0.16[/math]

[math]\frac{NH3,B\gt 0, PaddleNumber^{e^-}=5}{NH3,B\lt 0, PaddleNumber^{e^-}=10}=1.76 \pm 0.17[/math]


Notice the ratio is statistically the same if the NH3 and ND3 targets are used for this ratio in a manner which will be similar to the semi-inclusive analysis.

[math]\frac{ND3,B\gt 0, PaddleNumber^{e^-}=5}{NH3,B\lt 0, PaddleNumber^{e^-}=10}=1.55 \pm 0.15[/math]


The above ratios, which have been observed to be ammonia target independent, indicate a difference in electron detector efficiency when the Torus polarity is flipped. In order to make detector efficiency the same for electrons, the ratios were used as a "correction coefficient". The "correction coefficient" for the case [math]\frac{ND3,B\gt 0,E_{PaddleNumber}=5}{NH3,B\lt 0,E_{PaddleNumber}=10}[/math] is [math]0.645[/math] and for the [math]\frac{ND3,B\lt 0,E_{PaddleNumber}=10}{NH3,B\gt 0,E_{PaddleNumber}=5}[/math] it is [math]1.82[/math].



The measured single pion electroproduction rate was compared to the MAID 2007 unitary model that has been developed using the world data of pion photo and electro-production to determine the impact of using the above "correction coefficient". The model is well adopted for predictions of the observables for pion production, like five fold cross section, total cross secton and etc.

The MAID 2007 model has predictions of the total cross section for the following two cases, that are related to our work:

[math]\gamma[/math] + proton[math](NH_3) \rightarrow[/math] [math]\pi^+[/math] + neutron
[math]\gamma[/math] + neutron[math](ND_3) \rightarrow[/math] [math]\pi^-[/math] + proton


The ratio of the pions detected in the scintillator paddles, located between the Cherenkov counter and electromagnetic calorimeter, is shown on Fig. 3.2. The ratios were taken for four different cases. Assuming that, for the inbending case positive pions and for the outbending case negative pions have the same trajectories(the same kinematics) and vice versa((the inbending,negative pion) and (the outbending, positive pions)).

Paddenumbvsratio.jpg
Figure 3.2. Pion paddle number vs Ratio.

Using MAID 2007 total cross section was calculated for the following invariant mass and four momentum transferred square: [math]1.7\lt W\lt 1.8[/math] and [math]Q^2=1.1 GeV^2[/math]. <ref name="MAID2007" > http://wwwkph.kph.uni-mainz.de/MAID//maid2007/maid2007.html</ref>. After applying correction coefficients from inclusive cases, the ratios have been compared to the results from MAID2007.

[math]\sigma = \sigma_{T} + \epsilon \sigma_{L} + \sqrt{2\epsilon(1 + \epsilon)}\sigma_{LT} cos{\phi_{\pi}}^{CM} + \epsilon \sigma{TT} cos2{\phi_{\pi}}^{CM} + h \sqrt{2\epsilon (1-\epsilon)}\sigma_{LT^{\prime}}sin{\phi_{\pi}}^{CM} [/math]

After applying correction coefficients from inclusive cases, the ratios have been compared to the results from MAID2007.

NegativepionND3beforecorrection.png NegativepionND3aftercorrection.png

PositivepionND3beforecorrection.png PositivepionND3aftercorrection.png


Figure 3.3. Pion paddle number vs Ratio before and after correction.

Applied corrections are following:

[math]\frac{N(\pi^+,ND_3,B\lt 0)}{N(\pi^-,NH_3,B\gt 0)} \times 1.82[/math]

[math]\frac{N(\pi^+,ND_3,B\gt 0)}{N(\pi^-,NH_3,B\lt 0)} \times 0.645[/math]

[math]\frac{N(\pi^-,ND_3,B\gt 0)}{N(\pi^+,NH_3,B\lt 0)} \times 0.645[/math]

[math]\frac{N(\pi^-,ND_3,B\lt 0)}{N(\pi^+,NH_3,B\gt 0)} \times 1.82[/math]

Notes

<references/>


Media:NH3positivepionMAID.pdfMedia:ND3negativepionMAID.pdf

Nomenclature

[math]N(\pi^+,ND_3,B\lt 0[/math]) - Number of events in which [math]\pi^+[/math] and electron are detected in the final state using [math]ND_3[/math] target, when [math]B\lt 0[/math].
[math]N(\pi^+,ND_3,B\gt 0[/math]) - Number of events in which [math]\pi^+[/math] and electron are detected in the final state using [math]ND_3[/math] target, when [math]B\gt 0[/math].
[math]N(\pi^-,ND_3,B\lt 0[/math]) - Number of events in which [math]\pi^-[/math] and electron are detected in the final state using [math]ND_3[/math] target, when [math]B\lt 0[/math].
[math]N(\pi^-,ND_3,B\gt 0[/math]) - Number of events in which [math]\pi^-[/math] and electron are detected in the final state using [math]ND_3[/math] target, when [math]B\gt 0[/math].
[math]N(\pi^+,NH_3,B\lt 0[/math]) - Number of events in which [math]\pi^+[/math] and electron are detected in the final state using [math]NH_3[/math] target, when [math]B\lt 0[/math].
[math]N(\pi^+,NH_3,B\gt 0[/math]) - Number of events in which [math]\pi^+[/math] and electron are detected in the final state using [math]NH_3[/math] target, when [math]B\gt 0[/math].
[math]N(\pi^-,NH_3,B\lt 0[/math]) - Number of events in which [math]\pi^-[/math] and electron are detected in the final state using [math]NH_3[/math] target, when [math]B\lt 0[/math].
[math]N(\pi^-,NH_3,B\gt 0[/math]) - Number of events in which [math]\pi^-[/math] and electron are detected in the final state using [math]NH_3[/math] target, when [math]B\gt 0[/math].

[Go Back https://wiki.iac.isu.edu/index.php/Delta_D_over_D]