Difference between revisions of "Tamar AnalysisChapt"

From New IAC Wiki
Jump to navigation Jump to search
Line 104: Line 104:
 
[[Image:e_total_vs_nphe_momentum_cut_file_dst27070.gif|200px]]<br>
 
[[Image:e_total_vs_nphe_momentum_cut_file_dst27070.gif|200px]]<br>
  
===Cherenkov Theory===
+
===Cherenkov Counter Cut===
   
+
 
 +
The Cherenkov counter has been used to separate electrons from the background negatively charged pions. These negative pions are produced when lepton goes at polar angle close to zero and is not measured by the detector. The number of photoelectron distribution measured in the cherenkov detector and the energy deposition dependence on number of photoelectrons are shown in Figure ?.?. One can see, that a single photoelectron peak is caused by the misidentified pions as electrons.
 +
 
 
When the velocity of a charged particle is greater than the local phase velocity of light or when it enters a medium  with different optical properties the charged particle will emit photons.  The Cherenkov light is emitted under a constant angle <math>\theta_c</math> - the angle of Cherenkov radiation relative to the particle's direction. It can be shown geometrically that the cosine of the Cherenkov radiation angle is anti-proportional to the velocity of the  charged particle<br>
 
When the velocity of a charged particle is greater than the local phase velocity of light or when it enters a medium  with different optical properties the charged particle will emit photons.  The Cherenkov light is emitted under a constant angle <math>\theta_c</math> - the angle of Cherenkov radiation relative to the particle's direction. It can be shown geometrically that the cosine of the Cherenkov radiation angle is anti-proportional to the velocity of the  charged particle<br>
  
 
<math>cos \theta_c=\frac{1}{n \beta}</math><br>
 
<math>cos \theta_c=\frac{1}{n \beta}</math><br>
  
where <math>\beta c</math> is the particle's velocity and n - index of refraction of the medium.The charged particle in time t travels <math>\beta c t</math> distance, while the electromagnetic waves - <math>\frac{c}{n}t</math>. For a medium with given index of refraction n there is a threshold velocity <math>\beta_{thr}=\frac{1}{n}</math>, below <math>\beta_{thr}</math> no radiation can take place.  This process may be used to observe the passage of charged particles in a detector which can measure the produced photons.<br>
+
where <math>\beta c</math> is the particle's velocity and n - index of refraction of the medium. The charged particle in time t travels <math>\beta c t</math> distance, while the electromagnetic waves - <math>\frac{c}{n}t</math>. For a medium with given index of refraction n there is a threshold velocity <math>\beta_{thr}=\frac{1}{n}</math>, below <math>\beta_{thr}</math> no radiation can take place.  This process may be used to observe the passage of charged particles in a detector which can measure the produced photons.<br>
  
The number of photons produced per unit path length of a particle with charge ze and per unit energy interval of the photons is proportional to the sine of the Cherenkov angle[http://pdg.lbl.gov/2007/reviews/passagerpp.pdf] <br>
+
The number of photons produced per unit path length of a particle with charge <math>Ze</math> and per unit energy interval of the photons is proportional to the sine of the Cherenkov angle[http://pdg.lbl.gov/2007/reviews/passagerpp.pdf] <br>
  
 
<math>\frac{d^2 N}{dEdx}=\frac{\alpha z^2}{\hbar c}sin ^2 \theta_c=\frac{\alpha z^2}{\hbar c}[1-\frac{1}{\beta^2 n^2 (E)}]</math><br>
 
<math>\frac{d^2 N}{dEdx}=\frac{\alpha z^2}{\hbar c}sin ^2 \theta_c=\frac{\alpha z^2}{\hbar c}[1-\frac{1}{\beta^2 n^2 (E)}]</math><br>
Line 124: Line 126:
 
<math>\frac{d^2 N}{dEdx}=\frac{\alpha z^2}{\hbar c}sin ^2 \theta_c=\frac{\alpha z^2}{\hbar c}[\beta^2 n^2 (E) - 1]</math>
 
<math>\frac{d^2 N}{dEdx}=\frac{\alpha z^2}{\hbar c}sin ^2 \theta_c=\frac{\alpha z^2}{\hbar c}[\beta^2 n^2 (E) - 1]</math>
  
The gas used in the CLAS Cerenkov counter is perfluorobutane <math>C_4 F_{10}</math> with index of refraction equal to 1.00153.  
+
The gas used in the CLAS Cerenkov counter is perfluorobutane <math>C_4 F_{10}</math> with index of refraction equal to 1.00153. The number of photoelectrons emitted by electrons is about 13. On the other hand, calculations show that the number of photons produced by the negatively charged pions in the Cherenkov detector is approximately 2.  
 
 
 
 
 
 
====Electrons====
 
 
 
The calculation of the number of photoelectrons emitted by electrons is shown below.
 
Electron mass <math>m_e = 0.000511GeV</math>, <math>n=1.00153</math> and <math>\beta=\frac{p}{E} = 1</math>, because mass of the electron is negligible and also <math>\frac{\alpha}{\hbar c}=370[ eV^{-1} cm^{-1}] </math><br><br>
 
The Hall B cherenkov detector is <math>\sim \; 0.7</math> m thick radiator.  We assume the PMTs used to  collect light have a constant quantum efficiency of 8% for photons with wavelength between 300 and 600 nm.
 
 
 
<math>\frac{dN}{dx} </math> <math>= 2 \pi \alpha z^2 [{\beta^2 n^2 (\lambda)} - 1]\int_{300nm}^{600nm} \frac{1}{\lambda^2} d\lambda \times (0.08)=</math>  <br>
 
<math>= 2 \pi \alpha z^2 [{\beta^2 n^2 (\lambda)} - 1] (\frac {1}{\lambda})|_{600nm}^{300nm} \times (0.08) </math>=<br>
 
<math>= 2\times 3.14 \times \frac{1}{137} \times 1^2 \times [1.00153^2 - 1] \times \frac{1}{600} \times 0.08 [nm^{-1}] = </math> <br>
 
<math>= 19 \times 10^{-9} [nm^{-1}]</math> <br>
 
 
 
;For the number of photoelectrons we have the following <br>
 
<math> N = 19 \times 10^{-9} \times 0.7 m [\frac{10^9}{m}] = 13.3 </math><br>
 
 
 
That means the number of photoelectrons should be about 13.<br>
 
 
 
Used file dst27095_05.B00 energy=5.7GeV and torus=2250(B>0). Target NH3<br>
 
 
 
[[Image:e_momentum_vs_numb_of_photoelectrons_27095_theory.gif|200px]]
 
 
 
====Pions(<math>\pi^-</math>)====
 
 
 
 
 
The threshold energy for pions is ~2.5 GeV and for electrons 9 MeV.  <br>
 
 
 
;Example for <math>\pi^-</math> <br>
 
 
 
<math>m_{\pi^-}=139.57[\frac{MeV}{c^2}]</math>, momentum <math>p=mv=3 [\frac{GeV}{c}]</math> and n=1.00153<br>
 
 
 
where <math>\frac{\alpha}{\hbar c}=370[ eV^{-1} cm^{-1}] </math><br>
 
 
 
The Hall B cherenkov detector is <math>\sim \; 0.7</math> m thick radiator.  We assume the PMTs used to  collect light have a constant quantum efficiency of 8% for photons with wavelength between 300 and 600 nm.
 
 
 
<math>\frac{dN}{dx} = \frac{2 \pi \alpha z^2}{\lambda^2}[1-\frac{1}{\beta^2 n^2 (\lambda)}]=</math> <br>
 
<math>= 2 \pi \alpha z^2 [1-\frac{1}{\beta^2 n^2 (\lambda)}]\int_{300nm}^{600nm} \frac{1}{\lambda^2} d\lambda \times (0.08)=</math>  <br>
 
<math>= 2 \pi \alpha z^2 [1-\frac{1}{\beta^2 n^2 (\lambda)}] (\frac {1}{\lambda})|_{600nm}^{300nm} \times (0.08) </math>=<br>
 
<math>= 2\times 3.14 \times \frac{1}{137} \times 1^2 \times [1-\frac{1}{0.998919^2 \times 1.00153^2}] \times \frac{1}{600} \times 0.08 [nm^{-1}] = </math> <br>
 
<math>= 5.465 \times 10^{-9} [nm^{-1}]</math> <br>
 
 
 
;For the number of photons we have the following(for pions) <br>
 
<math> N = 5.465 \times 10^{-9} \times 0.7 m [\frac{10^9}{m}] = 3.8255 </math><br>
 
 
 
Used file dst27095_05.B00 energy=5.7GeV and torus=2250(B>0). Target NH3<br>
 
  
[[Image:pi_momentum_vs_numb_of_photons_27095_theory.gif|200px]]
+
[[Image:e_momentum_vs_numb_of_photoelectrons_27095_theory.gif|200px]][[Image:pi_momentum_vs_numb_of_photons_27095_theory.gif|200px]]
  
 
===CLAS Cherenkov signal===
 
===CLAS Cherenkov signal===

Revision as of 22:25, 18 February 2011

Introduction

THIS CHAPTER IS NOT READY< PLEASE READ OTHER THREE CHAPTERS.

This chapter will describe the techniques to calculate semi inclusive double spin asymmetries for the following reactions [math]e^-p \rightarrow e^- \pi^+ n[/math] and [math]e^-n \rightarrow e^- \pi^- p[/math] for NH3 and ND3 polarized targets respectively. The goal of this work is the extraction of the ratio of polarized to unpolarized quark distribution functions using double spin asymmetries.
There are couple of main steps in the data analysis, which will be discussed bellow: 1) the data selection, 2) background elimination, 3) the target and beam polarization, 4) the detector acceptance >???????


The CLAS Data Selection

The data files from the EG1b run were chosen for this analysis. During the run 2.2 GeV, 4.2 GeV and 5.7 GeV longitudinally polarized electron beam were used on the polarized frozen ammonia. This work will discuss the analysis of 4.2 GeV electron beam on hydrogen and deuteron targets.

Table ?? EG1b Runs used for Analysis.(Beam Energy, Target Type, Torus Current, Target and Beam Polarization)

The collected data during the experiment has been filtered by applying restrictions, which will be discussed in the Data Analysis chapter.

Beam Helicity and Charge Asymmetry

In double spin asymmetry measurements it is important to eliminate any source of false asymmetry, like charge asymmetry.

The helicity of the electron beam was flipped with a rate of 1 HZ and its original state was chosen at the injector in a pseudo-random way. The next state of the electron helicity, so called complement state, is always with opposite helicity. The diagram on Fig.?.? represents the helicity state.

Figure ??. The Helicity State

A one bit signal from the beam injector gives the helicity information, whereas a sync bit with a 2 HZ frequency is generated at the same time and is equal to the helicity flip time. However, the electron beam helicity sequence can be interrupted by the Data Acquisition System(DAQ), like leaving unpaired helicity states and dead time problems.

In order to eliminate the asymmetry caused by the DAQ, the Faraday cup was installed to measure the charge simultaneously with the helicity state flip. The stored data files include accumulated charge value for each helicity state.

Data Analysis

Radiative Corrections

Analysis chapter is not ready, Please start reading Experimental setup, theory and CEBAF. Those are the chapters that i finished, and would like them corrected. Ill rewrite the chapter below, it is not argent and include the diagrams.

In double spin asymmetry analysis the electron nucleon scattering process is given as an one photon exchange event, so called the Born approximation. In reality, there are multiple photon effects during the experiments. These high order processes, also called radiative effects, can be calculated and used to correct the cross section data.

The Born Approximation



There are two types of radiative corrections, internal and external. Internal radiative corrections describe the contributions which took place during the lepton-hadron interaction. In first order approximation they include vertex photon exchange, self energy and vacuum polarization.<ref name="RadiativeCorrections"> Nucleon Form Factors. In Scholarpedia, from http://www.scholarpedia.org/article/Nucleon_Form_factors#History </ref>

Vertex Photon Exchange
Vacuum Polarization
Self Energy
Self Energy



On the other hand, the external radiative corrections account for the Bremsstrahlung by the incoming and scattered electron and by the recoiling target nucleon.

Externall Radiation



One of the major advantages of the double polarization experiments is the minimum contribution from the radiative corrections. For the ~5.6 GeV incident electron beam data the radiative corrections are less than 5%<ref>http://www.jlab.org/Hall-B/secure/eg1/EG2000/fersch/QUALITY_CHECKS/file_quality/runinfo.txt</ref>. Due to negligible contributions from the radiative corrections they are not included in the double spin asymmetry analysis.

EC CUTS

During the EG1b experiment, the electromagnetic calorimeter was used to separate electrons from the pions. The energy deposition into the electromagnetic calorimeter is different for electrons and pions.

Electrons interact with the matter , producing electromagnetic showers and releasing the energy into the matter. The deposited energy into matter is linearly proportional to the momentum of the electrons. On the other hand , pions are minimum ionizing particles(MIP), losing energy with the rate of approximately [math]2 MeV g^{-1} cm^2[/math].

Electromagnetic calorimeter(EC) consists of 39 of 10 thick scintillator and ~2 mm thick lead. After passing through the calorimeter detector pions lose about 0.08 GeV energy independent their momentum. The energy loss by the electrons is linearly proportional to the momentum of the electron. Pions produce the constant signal in the calorimeter around 0.08 GeV. In order to eliminate misidentified pions from the electron sample, following cut has been applied:

[math]EC_{inner}\gt 0.08 \times p[/math]

where p represents particle momentum and [math]EC_{innner}[/math]- inner part of the calorimeter.

Electromagnetic calorimeter contains 13 layers of lead-scintillator sandwich. Each set of 13 layers are subdivided into 5 inner and 8 outer layers. Since the energy loss of pions is related to the detector thickness the correlation can be established between the energy deposition into the inner and outer layers of the detector:

[math]\frac{EC_{tot}}{EC_{inner}} = \frac{13}{5}[/math]

which gives the following cut for the energy deposition into the outer layer of the calorimeter:

[math]EC_{tot} \gt 0.2 \times p[/math]

E total vs e inner1 before cuts file dst27070.gif E total vs e inner1 after cuts file dst27070.gif

Figure ?.?. [math]EC_{inner}/p[/math] vs [math]EC_{tot}/p[/math] before and after EC cuts ([math]EC_{tot}\gt 0.2p[/math], for EC inner - [math]EC_{inner}\gt 0.08p[/math]).

NPHE CUTS

the misidentified electrons were excluded requiring a signal in the threshold CLAS Cherenkov detector. Pion's misidentified as electrons have been shown to produce around ~1.5 photoelectons 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 cuts the peak around 1.5 is substantially reduced.

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


E total vs nphe momentum cut file dst27070.gif

Cherenkov Counter Cut

The Cherenkov counter has been used to separate electrons from the background negatively charged pions. These negative pions are produced when lepton goes at polar angle close to zero and is not measured by the detector. The number of photoelectron distribution measured in the cherenkov detector and the energy deposition dependence on number of photoelectrons are shown in Figure ?.?. One can see, that a single photoelectron peak is caused by the misidentified pions as electrons.

When the velocity of a charged particle is greater than the local phase velocity of light or when it enters a medium with different optical properties the charged particle will emit photons. The Cherenkov light is emitted under a constant angle [math]\theta_c[/math] - the angle of Cherenkov radiation relative to the particle's direction. It can be shown geometrically that the cosine of the Cherenkov radiation angle is anti-proportional to the velocity of the charged particle

[math]cos \theta_c=\frac{1}{n \beta}[/math]

where [math]\beta c[/math] is the particle's velocity and n - index of refraction of the medium. The charged particle in time t travels [math]\beta c t[/math] distance, while the electromagnetic waves - [math]\frac{c}{n}t[/math]. For a medium with given index of refraction n there is a threshold velocity [math]\beta_{thr}=\frac{1}{n}[/math], below [math]\beta_{thr}[/math] no radiation can take place. This process may be used to observe the passage of charged particles in a detector which can measure the produced photons.

The number of photons produced per unit path length of a particle with charge [math]Ze[/math] and per unit energy interval of the photons is proportional to the sine of the Cherenkov angle[1]

[math]\frac{d^2 N}{dEdx}=\frac{\alpha z^2}{\hbar c}sin ^2 \theta_c=\frac{\alpha z^2}{\hbar c}[1-\frac{1}{\beta^2 n^2 (E)}][/math]


[math]\frac{d^2 N}{d\lambda dx} = \frac{2 \pi \alpha z^2}{\lambda^2}[1-\frac{1}{\beta^2 n^2 (\lambda)}][/math]

[math]\beta=\frac{v}{c}=\frac{pc}{\sqrt{(pc)^2 + (mc^2)^2}}[/math]

after deriving the Taylor expansion of our function and considering only the first two terms, we get
[math]\frac{d^2 N}{dEdx}=\frac{\alpha z^2}{\hbar c}sin ^2 \theta_c=\frac{\alpha z^2}{\hbar c}[\beta^2 n^2 (E) - 1][/math]

The gas used in the CLAS Cerenkov counter is perfluorobutane [math]C_4 F_{10}[/math] with index of refraction equal to 1.00153. The number of photoelectrons emitted by electrons is about 13. On the other hand, calculations show that the number of photons produced by the negatively charged pions in the Cherenkov detector is approximately 2.

E momentum vs numb of photoelectrons 27095 theory.gifPi momentum vs numb of photons 27095 theory.gif

CLAS Cherenkov signal

Electrons

The cherenkov signal measured in CLAS for particles identified as electrons by the tracking algorithm is shown below. There are two distributions present. One distribution is centered around 1.5 PEs and the second distribution is at 8 PEs when two gaussians and a Landau distribution are combined and fit to the spectrum. As we will show below, the first peak is due to the misidentification of a negative pion as an electron.

PE Fit equation
[math]N_{pe}= p_0 e^{-0.5 \left (\frac{x-p_1}{p_2} \right )^2} + p4\frac{1}{1-\left(\frac{x-p5}{p6}\right )} + p_6 e^{-0.5 \left (\frac{x-p_7}{p_8} \right )^2}[/math]
File:Gaussian fitting function.pdf
C.Lanczos, SIAM Journal of Numerical Analysis B1 (1964), 86.
File:C.Lanczos SIAM Journal of Numerical Analysis B1 1964 86.pdf


normalized fitting functions Landau and Landaun difference for me

//__[2] [3] ____________________________________________________________________________
Double_t TMath::Landau(Double_t x, Double_t mpv, Double_t sigma, Bool_t norm)
{
  // The LANDAU function with mpv(most probable value) and sigma.
  // This function has been adapted from the CERNLIB routine G110 denlan.
  // If norm=kTRUE (default is kFALSE) the result is divided by sigma
  static Double_t p1[5] = {0.4259894875,-0.1249762550, 0.03984243700, -0.006298287635,   0.001511162253};
  static Double_t q1[5] = {1.0         ,-0.3388260629, 0.09594393323, -0.01608042283,    0.003778942063};
  static Double_t p2[5] = {0.1788541609, 0.1173957403, 0.01488850518, -0.001394989411,   0.0001283617211};
  static Double_t q2[5] = {1.0         , 0.7428795082, 0.3153932961,   0.06694219548,    0.008790609714};
  static Double_t p3[5] = {0.1788544503, 0.09359161662,0.006325387654, 0.00006611667319,-0.000002031049101};
  static Double_t q3[5] = {1.0         , 0.6097809921, 0.2560616665,   0.04746722384,    0.006957301675};
  static Double_t p4[5] = {0.9874054407, 118.6723273,  849.2794360,   -743.7792444,      427.0262186};
  static Double_t q4[5] = {1.0         , 106.8615961,  337.6496214,    2016.712389,      1597.063511};
  static Double_t p5[5] = {1.003675074,  167.5702434,  4789.711289,    21217.86767,     -22324.94910};
  static Double_t q5[5] = {1.0         , 156.9424537,  3745.310488,    9834.698876,      66924.28357};
  static Double_t p6[5] = {1.000827619,  664.9143136,  62972.92665,    475554.6998,     -5743609.109};
  static Double_t q6[5] = {1.0         , 651.4101098,  56974.73333,    165917.4725,     -2815759.939};
  static Double_t a1[3] = {0.04166666667,-0.01996527778, 0.02709538966};
  static Double_t a2[2] = {-1.845568670,-4.284640743};
  if (sigma <= 0) return 0;
  Double_t v = (x-mpv)/sigma;
  Double_t u, ue, us, den;
  if (v < -5.5) {
     u   = TMath::Exp(v+1.0);
     if (u < 1e-10) return 0.0;
     ue  = TMath::Exp(-1/u);
     us  = TMath::Sqrt(u);
     den = 0.3989422803*(ue/us)*(1+(a1[0]+(a1[1]+a1[2]*u)*u)*u);
  } else if(v < -1) {
     u   = TMath::Exp(-v-1);
     den = TMath::Exp(-u)*TMath::Sqrt(u)*
            (p1[0]+(p1[1]+(p1[2]+(p1[3]+p1[4]*v)*v)*v)*v)/
            (q1[0]+(q1[1]+(q1[2]+(q1[3]+q1[4]*v)*v)*v)*v);
  } else if(v < 1) {
     den = (p2[0]+(p2[1]+(p2[2]+(p2[3]+p2[4]*v)*v)*v)*v)/
           (q2[0]+(q2[1]+(q2[2]+(q2[3]+q2[4]*v)*v)*v)*v);
  } else if(v < 5) {
     den = (p3[0]+(p3[1]+(p3[2]+(p3[3]+p3[4]*v)*v)*v)*v)/
           (q3[0]+(q3[1]+(q3[2]+(q3[3]+q3[4]*v)*v)*v)*v);
  } else if(v < 12) {
     u   = 1/v;
     den = u*u*(p4[0]+(p4[1]+(p4[2]+(p4[3]+p4[4]*u)*u)*u)*u)/
               (q4[0]+(q4[1]+(q4[2]+(q4[3]+q4[4]*u)*u)*u)*u);
  } else if(v < 50) {
     u   = 1/v;
     den = u*u*(p5[0]+(p5[1]+(p5[2]+(p5[3]+p5[4]*u)*u)*u)*u)/
               (q5[0]+(q5[1]+(q5[2]+(q5[3]+q5[4]*u)*u)*u)*u);
  } else if(v < 300) {
     u   = 1/v;
     den = u*u*(p6[0]+(p6[1]+(p6[2]+(p6[3]+p6[4]*u)*u)*u)*u)/
               (q6[0]+(q6[1]+(q6[2]+(q6[3]+q6[4]*u)*u)*u)*u);
  } else {
     u   = 1/(v-v*TMath::Log(v)/(v+1));
     den = u*u*(1+(a2[0]+a2[1]*u)*u);
  }
  if (!norm) return den;
  return den/sigma;
}

Fitting the Histograms

root [13] e_numb_of_photoelectrons->Draw();   
root [2] g3= new TF1("g3","gaus(0)+landau(3)+gaus(6)",0,20);

To get fit parameters for g3 we should fit individually each of them(gaus(0),landau(3),gaus(6))

root [3] g3->SetParameters(2.3e3,8,8.1,8.6e+3,8.6e-1,1.6,8.6e+3,8.6e-1,1.6,8.6e3)
root [11] e_numb_of_photoelectrons->Fit("g3","R+");

Notes

<references/>

[Go Back]