Difference between revisions of "Tamar AnalysisChapt"
Line 651: | Line 651: | ||
|- | |- | ||
|The systematic errors for the <math>A_{NH_{3}}^{\pi^+}</math> asymmetry. | |The systematic errors for the <math>A_{NH_{3}}^{\pi^+}</math> asymmetry. | ||
− | | | + | |}<br> |
Revision as of 06:20, 10 December 2012
This chapter describes the techniques used to analyze the data collected during the EG1b experiment and calculate semi-inclusive cross sections for the following reactions:
A measure of
can be used to test for independent fragmentation when it is compared to inclusive structure functions as shown in Eq. (1.54).The CLAS Data Selection
The data files from the EG1b experiment chosen for this analysis are listed in Table~\ref{tab:EG1bRuns}. During the experiment, 2.2 GeV, 4.2 GeV and 5.7 GeV longitudinally polarized electron beams were used to probe the polarized frozen ammonia NH
Run Set | Target Type | Torus Current(A) | Target Polarization | HWP |
28100 - 28102 | ND3 | +2250 | -0.18 | +1 |
28106 - 28115 | ND3 | +2250 | -0.18 | -1 |
28145 - 28158 | ND3 | +2250 | -0.20 | +1 |
28166 - 28190 | ND3 | +2250 | +0.30 | +1 |
28205 - 28217 | NH3 | +2250 | +0.75 | +1 |
28222 - 28236 | NH3 | +2250 | -0.68 | +1 |
28242 - 28256 | NH3 | +2250 | -0.70 | -1 |
28260 - 28275 | NH3 | +2250 | +0.69 | -1 |
28287 - 28302 | ND3 | -2250 | +0.28 | +1 |
28306 - 28322 | ND3 | -2250 | -0.12 | +1 |
28375 - 28399 | ND3 | -2250 | +0.25 | -1 |
28407 - 28417 | NH3 | -2250 | +0.73 | -1 |
28456 - 28479 | NH3 | -2250 | -0.69 | +1 |
Table 1.1. EG1b Runs analyzed for this work. |
Particle Identification
Additional tests were performed on the electron and a pion candidates reconstructed using the standard CLAS software package on the raw data collected during the EG1b experiment. Electrons are identified by matching the charged particle hits in the Cherenkov counter, electromagnetic calorimeter, and Time of Flight system. Geometrical and timing cuts are applied to improve electron identification. In addition, cuts are applied on the energy deposited by the particle into the calorimeter and the number of photoelectrons produced in the Cherenkov counter. Charged pions are identified by matching the hits in the drift chamber and ToF counter, along with a Cherenkov cut requiring that the number of photons for pions be less than two.
Electron Identification
The CLAS trigger system required the particle to deposit energy in the electromagnetic calorimeter and illuminate the Cherenkov counter within a
EC CUTS
The CLAS electromagnetic calorimeter was used to reduce the misidentification of electron and negative pion candidates. The electromagnetic calorimeter contains thirteen layers of lead-scintillator sandwiches composed of
Electrons interact with the calorimeter producing electromagnetic showers that release energy into the calorimeter. The deposited energy is proportional to the momentum of the electrons. Figure~\ref{fig:ECCuts} shows the correlation of the inner and outer calorimeter electron candidate's energy measured by the calorimeter and divided by the particles momentum reconstructed by the drift chamber. As shown in the Figure~\ref{fig:ECCuts}, there is an island near , which contains most of the electron candidates as well as some regions below which will be argued to be negative pions misidentified as electrons.
Pions entering the calorimeter are typically minimum ionizing particles, loosing little of their incident energy in the calorimeter at a rate of 2 . Electrons, on the other hand, deposit a larger fraction of their momentum into the calorimeter. As a result, the energy deposited into the electromagnetic calorimeter is different for electrons and pions. Pions loose about GeV of energy traversing the calorimeter independent their momentum thereby producing the constant signal in the calorimeter around GeV. In order to reduce misidentified pions from the electron sample, the following cut has been applied:
|
where
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:
|
which gives the following cut for the energy deposition into the outer layer of the calorimeter:
|
Figure 2.2 versus before and after EC cuts ( and ). After applying EC cuts about of the events have been removed from the electron sample. |
Cherenkov Counter Cut
The Cherenkov counter has been used to further reduce the negatively charged pion background in the reconstructed electron sample. 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. Cherenkov light is emitted at the critical angle
where
Taylor expanding Eq. 3.6 and keeping only the first two terms we get following
The gas used in the CLAS Cerenkov counter is perfluorobutane
Figure 2.3 Theoretical Calculation of the Number of Photoelectrons for electrons and pions. |
The distribution of the number of photoelectrons measured in the Cherenkov detector and the energy deposition dependence on number of photoelectrons are shown on Figure~\ref{fig:nphewithoutcuts} and Figure~\ref{fig:nphewithcuts}. One can see that a single photoelectron peak is caused by misidentifying pions as electrons.
Figure 2.4 The number of photoelectrons without cuts. |
Figure 2.5 The total energy deposited into the Calorimeter vs The Number of Photoelectrons. |
Geometric and Timing cuts
Negative pions may be produced when the lepton scatters at a polar angle close to zero and is not observed by the detector. In order to reduce the electron sample contamination due to those pions, geometrical cuts on the location of the particle at the entrance to the Cherenkov detector and time matching cuts have been developed by Osipenko, so called OSI cuts \cite{Osipenko}. For each CLAS Cherenkov detector segment the following cut has been applied
where
The pion contamination in an electron sample was estimated by fitting the number of photoelectron distribution using two Gaussian distributions convoluted with a Landau distribution \cite{Lanczos}:
The fits in Figure~\ref{fig:OSICuts}.(a) suggest that the pion contamination in the electron sample is
Figure 2.6 The number of photoelectrons before and after OSI Cuts. |
Pion Identification
Charged pions are identified using a coincidence hit in the drift chamber and Time-of-Flight (ToF) counter. Pions are separated from the other charged particles by looking at the particle momentum versus the
where
Using the above information (particle momentum from the drift chambers and the timing information from the ToF system), the mass squared of the charged particle was calculated and is shown on Figure~\ref{fig:PvsM2}. The pion mass band is around . To isolate charged pions from the rest of the particles, a cut on the momentum versus distribution has been applied \cite{Park2008}.
Figure 2.7 The charged particle momentum versus | distribution. The pion and proton bands are clearly separated.
Figure 2.8 The charged particle momentum versus mass squared distribution for the | electroproduction process. The bands around 0 and 1 represent pions and protons respectively <ref name="Park2008" >Park, K., Burkert, V. D., & Kim, W. (The CLAS Collaboration). (2008). Cross sections and beam asymmetries for -> in the nucleon resonance region for . Phys. Rev., C77, 015208.</ref>.
In addition to the charged particle velocity (
), the fiducial volume cuts have been applied for the charged pion identification. Since the drift chambers and scintillators are used for pion detection, the polar angle range where pions are detected is much larger than for electrons. For the EG1b experiment, pions were detected from to \cite{Park2008}. The pion identification code has been developed by Joshua Pierce \cite{JoshuaPierce}.Event Reconstruction Efficiency
The goal of this work is to measure the semi-inclusive asymmetry when an electron and a pion are detected in the final state. For this analysis, pions of opposite charge will be observed using the same scintillator paddles 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 reconstruction efficiency was evaluated in terms of the electron rate observed in two different scintillator paddles detecting the same electron kinematics.
Inclusive Electron Event Reconstruction Efficiency
The electron reconstruction efficiency for individual scintillator detectors using the 4.2 GeV EG1b data is investigated below. Only the electron is detected in the final state (inclusive case). The pion contamination in the electron sample was removed by applying the cuts described above. The electron paddle numbers 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 Figure~\ref{fig:ElectronKinematics}.
Ratios of the inclusive electron rate, normalized using the gated Faraday cup, detected in scintillator paddles # 5 and # 10 were measured. The two ratios are constructed to quantify the CLAS detector's ability to reconstruct electrons in scintillator paddle #5 using a positive Torus polarity and scintillator #10 using the negative Torus polarity.
Notice the above ratios are statistically the same. The semi-inclusive analysis to be performed in this thesis will be taking ratios using an
Figure 2.9. Electron Kinematics. |
The above ratios, which have been observed to be ammonia target independent, indicate a difference in an electron detector efficiency when the torus polarity is flipped. An electron detection efficiency "correction coefficient" is defined in terms of the above ratio and measured to be
Exclusive and Semi-Inclusive Event Reconstruction Efficiencies
After determining the electron reconstruction efficiency for the selected paddle numbers, 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 section, etc.
The MAID 2007 model has predictions of the total cross section for the following two cases that are related to our work:
+ proton + neutron |
+ neutron + proton. |
The ratio of the pions detected in the scintillator paddles, located between the Cherenkov counter and electromagnetic calorimeter, is shown in Figure~\ref{fig:PionPaddleNumbervsRatio}. The ratios were taken for four different cases. The intrinsic assumption is that, for the inbending case, positive pions and for the outbending case negative pions have the same trajectories with the same kinematics. In addition, negatively charged pions in the inbending field and positively charged pions in the outbending fields are detected by the same detector elements.
Figure 2.10. Pion paddle number versus Ratio for Semi-Inclusive case. |
Using MAID 2007, the total cross section was calculated for the following invariant mass and four momentum transferred square values: 1.7 GeV < W<1.8 GeV and
where
Pion Paddle Number versus MAID2007 - Experiment | . The Black and red data represent B>0/B<0 and B<0/B>0 cases respectively before corrections. The green and blue points represent the ratios for B>0/B<0 and B<0/B>0 after inclusive corrections.
Asymmetries
The double spin asymmetry measurements in this thesis are performed by comparing scattering events that occur when the incident probe spin and nuclear target spin are parallel to the scattering events that occur when the spins are anti-parallel.
Beam Charge Asymmetry
The helicity of the electron beam was flipped at a rate of 1 Hz. The helicity is prepared at the source such that helicity pairs are produced pseudo randomly.
Figure 2.12. 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. |
If the first electron bunch is pseudo randomly chosen to be positive (negative) then it is labeled as the original helicity state and denoted in software by a 2 (1). The next helicity state is prepared to be a complement to the first state and labeled in the software as either a 4, if the original helicity state was a 1 (negative), or 3 if the original helicity state was a 2 (positive). The helicity selection process is then repeated.
Figure~\ref{fig:HelicityState} illustrates the signals used to label the helicity states. The clock pulse (SYNC) is used to indicate that a change in the pockel cell used to define the helicity state may have occurred. The helicity bit identifies the helicity state that was set. The original/complement pulse identifies if the state is an original or complement helicity state. All three bits are recorded in the raw data file for each event and then converted to the labels 1, 2, 3, 4 during DST file production once the particles have been reconstructed.
Two scalers were used to record several ancillary detectors, such as a Faraday cup and several PMTs mounted on the beam line, according to their helicity label. One of the scalers was gated by the DAQ live time in order to record beam conditions when the DAQ was able to take data and not busy recording data. The second scaler remained ungated. Both scalers recorded the SYNC and Helicity signals from the injector along with the counts observed from ancillary detectors during the SYNC interval. The Faraday cup signal recorded by the gated helicity scaler is used to normalize the events reconstructed during the same helicity interval. The beam charge asymmetry below is measured by the gated helicity scaler. For each run number, a gaussian fit was used to fit the beam charge asymmetry distributions as in Figure~\ref{fig:BeamAsym28101}. The beam charge asymmetry is defined as
where
Figure 2.14. Beam charge asymmetry for run #28101 using the gated Faraday cup counts for two helicity pairs(1-4 and 2-3 helicity pairs). and . |
EG1b data sets, with the same half wave plane, target type, target polarization and beam torus, have been combined. The beam charge asymmetries have been calculated for the each run group and are listed in Table~\ref{tab:RunvsBeamAsym}.
Run Group | Half wave plane(HWP) | ||
28100 - 28105 | +1 | ||
28106 - 28115 | -1 | ||
28145 - 28240 | +1 | ||
28242 - 28284 | -1 | ||
28286 - 28324 | +1 | ||
28325 - 28447 | -1 | ||
28449 - 28479 | +1 |
Table 4.1. Run Group versus Beam Charge Asymmetry. |
Electron Asymmetry
A measurement of the electron cross section helicity difference needs to account for a possible helicity dependence of the incident electron flux (charge Asymmetry). Figure~\ref{fig:RunvsEAsym} shows the reconstructed electron asymmetry before it is normalized by the gated Faraday Cup as a function of the run number for the 4.2 GeV data set. The reconstructed electron asymmetry can be defined following way:
or |
where
Figure Run Number versus Electron Asymmetry before FC normalization. The black and red points represent reconstructed electron asymmetry for the helicity 1-4 pair for ND | and NH target respectively. The blue and green points represent the helicity pair 2-3 for ND and NH respectively. The green line shows the sign of the half wave plane (HWP) and the purple line is the sign of the target polarization (TPol).
Systematic effects on the asymmetry measurement may be investigated by separating the data into two groups based on which helicity state is set first. The first group (black data points) represents the electron asymmetry observed when the first (original) helicity state is negative and its complement state is positive (helicity state #1 - state #4). The second group (red data points) represents the asymmetry observed when the first state is positive and the complement state is negative (helicity state #2 - #3). Both groups were divided into two subgroups based the target type used. The diamond points on the histogram represent the data for the NH
NH3, Tpol>0 and HWP>0. | NH3, Tpol<0 and HWP>0 |
W versus | . The electron asymmetry (sign(hel42-hel13)) changes sign when the HWP or Target Polarization sign is changed.
The un-normalized reconstructed electron asymmetry has been calculated as:
and normalized by the Faraday cup
Figure Run Number versus Electron Asymmetry after applying FC normalization. The black and red points represent the reconstructed electron asymmetry for the helicity 1-4 pair for ND | and NH target respectively. The blue and green points are the helicity pair 2-3 for ND and NH respectively. The green line shows the sign of the half wave plane (HWP) and the purple line is the sign of the target polarization (TPol).
Semi-Inclusive Asymmetries
The asymmetries from semi-inclusive pion electroproduction using proton or deuteron targets can be written in terms of the difference of the yield when the electron spin is parallel and antiparallel to the spin of the nucleon. There are four combinations of semi-inclusive asymmetries:
where
The kinematic coverage for the events used in the measured asymmetries are shown on Figure~\ref{fig:InvariantMassvsQ2} and Figure~\ref{fig:MissingMassM_X}. The semi-inclusive asymmetries are listed in Table~\ref{tab:RunvsSIDISAsym}. The first group represents the asymmetry measured when the first original electron spin (hel1) is antiparallel to the target nucleon spin and its complement state (hel4) is parallel, whereas the second group represents measuring the SIDIS asymmetry when the first helicity state is parallel (hel2) and its complement state is antiparallel to the spin of the nucleon. For the final measurement, the two groups are combined into one set. Instead of looking at original and complement electron helicity states, they have been combined into the positive (hel42) and negative (he13) helicity states, which are parallel and antiparallel to the spin of the target nucleon respectively.
Invariant Mass versus | .Missing Mass. |
Target type, Beam Torus | |||
NH | , B>0,|||
NH | , B<0,|||
ND | , B>0,|||
ND | , B<0,|||
NH | , B>0,|||
NH | , B<0,|||
ND | , B>0,|||
ND | , B<0,
Figure Run Number versus SIDIS Asymmetry for Each Type Target material and Beam Torus. |
The ratio of the combined semi-inclusive deep inelastic asymmetries (
) from Table~\ref{tab:RunvsSIDISAsym}, for two different torus settings have been corrected for the electron reconstruction efficiency. The Figure~\ref{fig:TheRatioOfSIDIS} represents the asymmetry ratios before and after electron reconstruction efficiency corrections. The ratios have been calculated for each target and charged pion type. The SIDIS asymmetries before and after the correction are statistically the same. The result indicates that the electron reconstruction efficiency does not change the asymmetries.
Figure The ratio of the SIDIS asymmetries for two torus field settings | versus target and the charged pion type. The black squares represent the data before electron reconstruction efficiency and the red data - after electron reconstruction efficiency has been applied.
The data have been combined for each target type and asymmetries measured for two
SIDIS Asymmetry | ||
Figure SIDIS Asymmetries for | and .
In addition to
, the data have been subdivided according to the fraction energy of the observed final state pion( ). asymmetries for two are shown in Table~\ref{tab:SIDISz}.0.4 | ||
0.7 |
SIDIS Asymmetry. |
Dilution Factor
To exclude the contributions of polarized nucleons from the non-hydrogen nuclei in the ammonia target (Nitrogen) and the cooling material of the target (Helium) cell to the semi-inclusive rates, the measured raw double spin asymmetries have to be divided by a dilution factor. The dilution factor accounts for the fraction of events coming from the desired polarized target nucleon. The dilution factors are calculated by combining the data from runs using different target types. During the EG1b experiment, several runs were taken with Carbon C
In order to calculate the dilution factor, we need to define the number of counts with beam helicity negative and positive for ammonia, Carbon and empty target runs. All are weighted by the corresponding gated Faraday cup counts \cite{Khun2010}. Normalized rates are defined as
, |
where
The counts for all four targets (empty, Carbon and ammonia) can be expressed as the sum of counts from the entrance and exit window foils ( , , ), liquid Helium coolant ( , , ), Carbon( , , ), Nitrogen ( , , ), Hydrogen (Deuterium) ( , , ).
, |
, |
, |
where
Using above system of equations, we define two new spectra to account for the Carbon target and the difference in the amount of Helium in ammonia targets vice versa the Carbon target. For a target cell with length L and the Carbon target with length , the Carbon and LHe contributions can be written following way:
and
, |
where
We need to establish how the rate from the Carbon target is related to the rate from the N in NH and ND . The cross section for N can be written in terms of the cross sections on C target material ( ) and on a bound neutron in N ( ).
. |
It is assumed that when the scattering occurs on protons inside the target material, the ratio of
Using above quantities the background represented by the number of counts due to the non-hydrogen and non-deuterium parts of the ammonia target can be expressed as
. |
The dilution factor
. |
Item | Description | Value |
Density times target length for empty target. | Al: 167 | m; 0.045 g/cm . Kapton: 384 m; 0.055 g/cm . Total=Al + Kapton=0.0996 g/cm|
Density times target length for Carbon target. | 0.498 g/cm | |
f | The ratio of counts from foils to the C | slab in the Carbon target.|
He density. | 0.145 g/cm | |
The length of the target cell from the entrance to exit foil. | 1.90 cm | |
C | density.2.17 g/cm | |
Carbon target length. | 0.23 cm | |
NH | density0.917 g/cm | |
ND | density1.056 g/cm | |
Ammonia target length. | 0.6 cm |
Length and density values for different types of target material reproduced from the EG1b experiment \cite{Khun2010}.} |
The fractional energy of the observed final state hadron(z) | Dilution Factor( | )
NH | , &&|
NH | , &&
Calculated dilution Factor for NH | target type.
The dilution factor is compared below in Table~\ref{tab:CompareDilution}:
Reaction | Dilution Factor( | )
Resonance region | 0.11 - 0.13 |
Inclusive | 0.14 - 0.17 |
Semi-Inclusive | 0.122 - 0.182 |
{Dilution Factor compared with other results \cite{Prok2004} \cite{Fersch2008}. |
Fragmentation
A test of fragmentation can be performed by calculating the ratio of the difference of polarized to unpolarized cross sections for proton and neutron targets
= \frac{\Delta \sigma_{p}^{\pi^{+} + \pi^{-}}}{\sigma_{p}^{\pi^{+} + \pi^{-}} - \sigma_{n}^{\pi^{+} + \pi^{-}}} - \frac{\Delta \sigma_{n}^{\pi^{+} + \pi^{-}}}{\sigma_{p}^{\pi^{+} + \pi^{-}} - \sigma_{n}^{\pi^{+} + \pi^{-}}} |
. |
The left side of the fragmentation function is calculated using the measured double spin asymmetries for the ammonia targets. However, the right side can be extracted using the model.
0.4 | 14 % | 12 % |
0.7 | 21 % | 32 % |
Statistical Z - test for the data comparison with the model. The probability of not observing the data point for each | and values within the Inclusive Model.
The pion asymmetries,
Systematic errors
Systematic errors from the experimental setup and the cuts used for particle identification and background elimination were estimated. The systematic errors associated with the electromagnetic cuts and fiducial cuts for electron identification were calculated by comparing the SIDIS asymmetries before and after cuts were applied. Systematic effects related to the dilution factor and polarization were estimated following way: First, the SIDIS asymmetries were calculated for the standard values of the dilution factor and polarization. Then asymmetries were recalculated by changing the value of the each parameter by the amount of its uncertainty. The difference between these two values is the systematic effect. The systematic errors for the SIDIS asymmetries are shown in Table~\ref{tab:SystematicErrors}.
Source | Systematic Error |
Geometrical and timing cuts | |
Electromagnetic calorimeter cut | |
Dilution Factor | |
Target and beam polarization | |
Total |
The systematic errors for the | asymmetry.