Measurement of Neutrino-Nucleon Neutral-Current Elastic Scattering Cross-section at SciBooNE

Size: px
Start display at page:

Download "Measurement of Neutrino-Nucleon Neutral-Current Elastic Scattering Cross-section at SciBooNE"

Transcription

1 Measurement of Neutrino-Nucleon Neutral-Current Elastic Scattering Cross-section at SciBooNE Hideyuki Takei Thesis submitted to the Department of Physics in partial fulfillment of the requirements for the degree of Doctor of Science at Tokyo Institute of Technology February, 2009

2 2

3 Abstract In this thesis, results of neutrino-nucleon neutral current (NC) elastic scattering analysis are presented. Neutrinos interact with other particles only with weak force. Measurement of cross-section for neutrino-nucleon reactions at various neutrino energy are important for the study of nucleon structure. It also provides data to be used for beam flux monitor in neutrino oscillation experiments. The cross-section for neutrino-nucleon NC elastic scattering contains the axial vector form factor G A (Q 2 ) as well as electromagnetic form factors unlike electromagnetic interaction. G A is propotional to strange part of nucleon spin ( s) in Q 2 0 limit. Measurement of NC elastic cross-section with smaller Q 2 enables us to access s. NC elastic cross-sections of neutrinonucleon and antineutrino-nucleon were measured earlier by E734 experiment at Brookheaven National Laboratory (BNL) in In this experiment, cross-sections were measured in Q 2 > 0.4 GeV 2 region. Result from this experiment was the only published data for NC elastic scattering cross-section published before our experiment. SciBooNE is an experiment for the measurement of neutrino-nucleon scattering cross-secitons using Booster Neutrino Beam (BNB) at FNAL. BNB has energy peak at 0.7 GeV. In this energy region, NC elastic scattering, charged current elastic scattering, charged current pion production, and neutral current pion production are the major reaction branches. SciBar, electromagnetic calorimeter, and Muon Range Detector are the detectors for SciBooNE. The SciBar consists of finely segmented scintillators and channels of PMTs. It has a capability to reconstruct particle track longer than 8 cm and separate proton from muons and pions using energy deposit information. Signal of NC elastic scattering is a single proton track. In νp νp process, the recoil proton is detected. On the other hand, most of νn νn is invisible because there are only neutral particles in final state, but sometimes recoil neutron is scattered by proton and recoil proton is detected. Signal of this event is also single proton track. Event selection for the single proton track events using geometrical and de/dx information of reconstructed track is performed. After the event selection, NC elastic scattering data sample is obtained. They includes νp νp and νn νn is obtained. Absolute cross-section as a function of Q 2 is evaluated from this NC elastic scattering data sample.

4

5 Contents 1 Introduction 13 2 Neutrino-Nucleon Interaction Neutrinos Charged Current Interactions Neutral Current Interactions Neutrino-Nucleon Neutral Current Elastic Scattering Cross-Section Strangeness Component of Nucleon Spin Experimental Setup of SciBooNE Experimental Setup Overview Booster Neutrino Beam Proton Beam Horn and Target Detectors SciBar Detector Electromagnetic Calorimeter Muon Range Detector Triggers Neutrino Events in the SciBar Detector Data Taking Beam Simulation Neutrino Interaction Simulation Elastic Scattering Resonant Single Meson Production Coherent Pion Production Nuclear Effects in Final State Detector Simulation Dirt Simulation

6 4 Analysis of Neutral Current Elastic Scattering Track Reconstruction Monte Carlo Normalization Event Selection Geometrical Selection Particle Identification Background Estimation Event Selection Summary Detector Unfolding Efficiency Correction Cross Section Systematic Uncertainty Neutrino Beam Neutrino Interaction Detector Response Total Systematic Uncertainty Nuclear Effect Correction Result 91 6 Conclusion 97 4

7 List of Figures 2.1 Feynman diagrams of CC elastic and CC 1-π production reaction Cross-section for charged current interactions as a function of neutrino energy. Upper plot is for neutrino and lower plot is for antineutrino Feynman diagrams of NC elastic and NC π 0 production reaction Kinematics of νp νp Feynman diagram of electron-proton deep inelastic scattering The ν µ beam flux of E734 experiment. Dots show the measured ν µ flux. The solid curve is a MC beam flux prediction The ν µ beam flux of E734 experiment. Dots show the measured ν µ flux. The solid curve is a MC beam flux prediction (a) Differential cross-section for νp νp (upper plots) and νp νp (lower plots) measured by E734 experiment. (b) Ratio of the cross-section for νp νp to the cross-section for νp νp. Solid lines and dashed lines show the result of fit. Error bars include statistic and Q 2 dependent systematic error. [1][10] Experimental setup of SciBooNE. SciBooNE detector is located 100 m downstream of the proton target. MiniBooNE detector is located 440 m downstream of the SciBooNE detector Air view of the Booster, target station, decay region and Sci- BooNE detector SciBar, an electro-magnetic calorimeter (EC), and a muon range detector (MRD) are detectors for SciBooNE. They are located 100 m downstream of the proton target

8 3.4 The SciBar detector consists of 14,336 of the combination of plastic scintillator bar and wave-length shifting (WLS) fiber. Volume of the detector is m 3, weight is 15 tons.the Scintillator bars are arranged vertically and horizontally. Each WLS fiber is attached to one anode of a multi-anode PMT The WLS fibers are attached to the MAPMT by the cookie. The light injection module is located between the scintillators and the MAPMT. Each MAPMT is mounted on the front-end board The EC consists of 1 mm diameter scintillating fibers embedded in lead foil. The calorimeter is made of modules of dimensions cm 3. Each module is read out by two 1 inch Hamamatsu PMTs per side The MRD is installed downstream of the EC. It has 12 iron plates with thickness of 5 cm and 13 plastic scintillator planes with thickness of 6 mm A neutrino event candidate on online display. Vertical (left) and horizontal (right) readout of detectors are shown Event displays of CCQE candidate. The area of dots is proportional to ADC channels. Boxes represent TDC hit information. Tracks penetrating the MRD are identified as a muon Event displays of CCQE candidate. The decay of the muon is tagged by multiple TDC hits de/dx distributions of muons (left plot) and protons (right plot). Protons are separated from muons with 90% efficiency CC 1π coherent or νn µ nπ + candidate. Muon and pion are visible CC 1π νp µ pπ + candidate. Muon, pion and proton are visible Candidate of neutral current elastic scattering. Signature of this event is single proton track. The proton is identified using de/dx information Candidate of neutral current π 0 production event. The 2 γ s from a π 0 decay produce electro-magnetic showers, which can be detected in the SciBar Accumulated number of POT (solid line) and number of POT which passed all the data quality cuts (dotted line) Number of CC event candidates normalized by POT. The flat shape for each mode and the same values for the two antineutrino mode periods shows that the SciBar took data stably throughout the beam time

9 3.18 Neutrino flux prediction at the SciBooNE detector as a function of neutrino energy E ν The Q 2 dependence of the differential cross-section for NC elastic scattering adopted in NEUT. Black line and red line shows differential cross-sections for neutrino-proton scattering and neutrino-neutron scattering respectively Probability of nucleon-nucleon interaction in 16 O after the neutrino interaction. No interaction (solid line), elastic scattering (dashed line), single π production (dot-dot line), and two π production (dash-dod line) are plotted Mean deflection angle. It has a peak around 45 degrees Proton-proton total and elastic cross sections Proton-neutron total and elastic cross sections Geometry of dirt for the dirt simulation. Gray areas indicate dirt volume. The geometry of m 3, taking detector origin (red point) as a center were considered Vertex distributions of dirt events from simulation. Only events which made one track contained in SciBar were selected. Geometry is large enough for NC elastic analysis A candidate of neutral current elastic scattering from Sci- BooNE data. Signature of this event is single proton track. The proton is identified using stopping range and de/dx information Flow chart of this analysis Total energy deposit of track-unrelated hits on top and side of SciBar. Data are compared to MC simulation. The threshold to select one track events is set to 30 MeV for both top and side view. In these plots, dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events TDC hits from muon decay in top and side view. There are multiple TDC hits on both top and side view (case 1) TDC hits from muon decay in top and side view. There are multiple TDC hits only on one view (case 2) Time distribution of muon decay sample. Horizontal axis is the time interval between primary TDC hits and secondary TDC hits that matched on top and side view

10 4.7 Track length v.s. de/dx. Red points represent protons and black points represent other particles such as µ ±,π ± and e ±. Line shows the cut for proton identification. Particles above the line are identified as protons de/dx distributions for several different track lengths. The top left is 0-10 cm, the bottom right is cm. de/dx is calculated from total energy deposit of the track divided by track length. The peak at large de/dx is from proton track while the peak at small de/dx is from other MIP particles. Data are compared to MC simulations. In these plots, dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events Angular distribution of NC elastic sample. θ is an angle of the track with respect to the z direction. Data are compared to MC simulation. Dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events. Right plot shows data/mc ratio Length distribution (left figure) and data/mc ratio (right figure) of NC elastic sample. Data are compared to MC simulation. Dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events. Right plot shows data/mc ratio True proton kinetic energy v.s. true proton track length from MC simulation. A line is a result of fit Distribution of proton kinetic energy. Proton kinetic energy is reconstructed from track length. Data are compared to MC simulation Vertex position along z axis v.s. track length of data (left), Monte-Carlo simulation in SciBar (middle) and dirt Monte- Carlo simulation (right). SciBar is divided into 8 bins along z direction to compare number of events at upstream part with that of downstream part. Fifth bin is taken as a reference and compared with first, second and third bins Track length distributions of (upstream bins)-(fifth bin). Data agree with Monte-Carlo simulation within statistic error Proton kinetic energy distribution after background subtraction

11 4.16 Proton kinetic energy reconstructed from track length using Monte-Carlo simulation v.s. proton kinetic energy reconstructed from measured track length using tracking algorithm (left figure). Both quantities are from MC simulation. From this relationship, unfolding matrix was made (right figure) Proton kinetic energy distribution after detector unfolding. Data are compared to MC simulation MC true proton kinetic energy v.s. detection efficiency Proton kinetic energy distribution after efficiency correction. Data are compared to MC simulation Differential cross section as a function of proton kinetic energy Systematic uncertainty in the NC elastic cross section due to the uncertainty in the neutrino beam The systematic uncertainty due to M A uncertainty The systematic uncertainty due to the uncertainty in the detector response Total systematic uncertainty in NC elastic scattering cross section Q 2 reconstructed from proton kinetic energy (left) v.s. Q 2 calculated from the neutrino momentum using Monte-Carlo simulation (right) Q 2 distribution after nuclear effect correction. Data are compared to MC simulation Differential cross section as a function of Q 2. Error bars include statistical error only Total systematic uncertainty in NC elastic scattering cross section Differential cross section as a function of Q 2. Error bars include statistical error only Differential cross section after corrections. Cross section from data is compared with MC simulation. MC plot is from NC elastic sample and processed in the same way as data. Error bars include statistical error only Figure 5.4 shows differential cross section compared with the experimental data from BNL E734 and MiniBooNE experiments. Error bars of E734 and MiniBooNE data include statistical and systematic errors. BNL E734 data shows the cross section for νp νp scattering while SciBooNE and Mini- BooNE data show the cross section per nucleon (νn νn)

12

13 List of Tables 2.1 The values of the fitting parameters G s A(0),F s 1(0) and F s 2(0) and M A Particles which have roots in dirt events and make NC elastic scattering-like track in SciBar Fractions of the neutrino interactions in the dirt which make a proton track in SciBar Event selection summary table of the NC elastic event selection

14

15 Chapter 1 Introduction Neutrino is a particle introduced by Pauli in Electrons emitted in nuclear beta decay were expected to have a discrete energy for two-body decay, but they were discovered to have continuous energy spectrum. To maintain energy conservation, existence of neutral, light particle was suggested by Pauli. The particle is called neutrino. Neutrino is a particle with spin 1/2. The neutrino is left-handed while anti-neutrino is right-handed. Neutrino interacts with nucleon only by weak interaction. The interaction is mediated either by a W boson or Z boson. Neutral Current (NC) interaction is an exchange of a Z boson while Charged Current (CC) interaction is an exchange of a W boson. Weak interaction only couples to left-handed particle. In other words, parity is violated. Neutrino-nucleus scattering cross sections have been measured by several experiments. However, statistics of experimental data is low in energy region around 0.7 GeV. In this energy region, several channels of neutrino interaction exist. In NC elastic scattering, there are one neutrino, one nucleon and nothing else in the final state (νn νn). This channel is one of the main components of neutrino-nucleon interaction in the neutrino energy region around 0.7 GeV. In this energy region, the fractions of NC is 28.7%, CC is 71.3%. Fractions of NC elastic scattering is 17.6%, NC pion production is 11.0%, CC elastic scattering is 40.7% and CC pion production is 30.7%. Form factor describes a difference of cross sections between point-like target and target with actual size. Cross section for neutrino-nucleon NC elastic scattering contains electro-magnetic form factor and axial form factor. Each form factor includes strangeness part. Strangeness part of axial form factor G A is proportional to strangeness component of nucleon spin ( s) in the Q 2 0 limit. Quark spin contribution to the nucleon spin was measured by EMC experiment in It was only 12±9±14 % though it was expected 13

16 to be 100 %. Cross section for NC elastic scattering was measured by E734 experiment at BNL in 1980 s [1]. From the analysis of the data from this experiment, s was evaluated to be between 0.21 and 0. This experiment proved that s can be measured by neutrino-nucleon scattering. However, the data from this experiment contains large systematic and statistic errors. One of the causes of large systematic error is that E734 was not able to measure cross section with Q 2 less than 0.4 GeV 2. New experiments are necessary to determine s with a better accuracy. SciBooNE can measure cross section with Q 2 down to 0.2 GeV 2 using fully-active, finely segmented tracking detector called SciBar. It uses a high intensity neutrino beam from Booster Neutrino Beam at FNAL. CC interactions are also measured at SciBooNE. Cross sections for CC interactions are independent of s. Taking ratio of cross sections of NC elastic events and CC events helps to cancel systematic errors in the beam profile. Furthermore, it is important for neutrino experiments such as neutrino oscillation experiments to know cross sections of signal and background channels precisely. SciBooNE can also provide neutrino flux information to neutrino oscillation experiment MiniBooNE which is on the same axis of the neutrino beam at FNAL. It improves sensitivity to the ν µ (ν µ ) disappearance experiment [12]. 14

17 Chapter 2 Neutrino-Nucleon Interaction 2.1 Neutrinos Neutrinos are neutral, light particles with spin 1/2. There are three flavors of neutrino: ν e,ν µ, and ν τ. They corresponds to the flavors of charged leptons, e,µ and τ. Neutrinos interact only via W ± and Z 0 bosons. Neutral current (NC) is a process where a Z 0 boson is exchanged. Charged current (CC) is a process where a W ± boson is exchanged. Weak interaction has much less probability to take place compared to other interactions because the masses of W ± and Z 0 are large (80 GeV for W ±, 91 GeV for Z 0 ) and needs a large amount of energy to produce. Neutrinos have two state of handedness (chirality), left-handed and righthanded state. In the weak interaction, left-handed particles and right-handed anti-particles interact in different way. Interaction via W + boson only couples to left-handed particle while interaction via W boson couples to righthanded anti-particle. In other words, parity is violated. Helicity is the projection of the spin onto its momentum vector. It is written as h = (S, p ) p (2.1) where S is the spin and p is the momentum of the particle. If the projection of the spin has the same direction as the momentum, helicity is positive. If the projection of the spin has opposite direction as the momentum, helicity is negative. For massless particles such as photon or gluon, helicity is Lorentz invariant because there is no frame that the projection of the spin is inverted. Their helicity is invariant and identical to chirality. For massive particles, helicity is not Lorentz invariant and distinguished from chirality. 15

18 Only left-handed neutrinos and right-handed antineutrinos have been observed. It is considered that only left-handed neutrinos and right-handed antineutrinos exists. This fact is included in the Standard Model. Neutrinos is predicted to be massless in Standard Model. However, the experimental data from neutrino oscillation experiments showed that neutrinos have mass. The Standard Model must be extended to include the missing chiral states. Neutrinos with the missing chiral states are called sterile neutrinos and do not interact via W ± or Z Charged Current Interactions Charged current interactions are weak interactions via W ± bosons. In neutrinonucleon scattering, neutrinos exchange charge with the nucleon and turn into charged leptons. Figure 2.1 shows Feynman diagrams of CC elastic and CC 1-π production reaction for ν µ. In an experiment, CC interactions can be distinguished from NC interactions by finding charged leptons in the final state. Cross-sections for CC interactions were measured by several experiment. However, in the neutrino energy region below 1 GeV, there are very low statistics for neutrino cross-section and no data for antineutrino cross-section (Figure 2.2). SciBooNE measures cross-sections in this energy region. It gives important information to neutrino experiments planning to use this energy region such as neutrino oscillation experiment T2K. Figure 2.1: Feynman diagrams of CC elastic and CC 1-π production reaction Neutral Current Interactions Neutral current interaction was first observed in 1973 by the Gargamelle bubble chamber at CERN [3][4]. They observed NC deep-inelastic scattering νn νx and an elastic electron scattering event νe νe using ν µ and ν µ. 16

19 Figure 2.2: Cross-section for charged current interactions as a function of neutrino energy. Upper plot is for neutrino and lower plot is for antineutrino. 17

20 Neutrino-nucleon NC elastic scattering was observed by the Columbia- Illinois-Rockfeller (CIR) [5] and Harvard-Pennsylvania-Wisconsin (HPWB) [6][7] collaborations at Brookhaven National Laboratory (BNL). They observed up to 100 events of νp νp (νp νp). These experiments were able to set bounds on electroweak theory by constraining sin 2 θ W. It was discovered by HPBW collaboration that the Lorentz structure of the neutral weak current is vector - axial-vector after taking more statistics. There are several channels of neutrino-nucleus neutral current interaction such as elastic, pion production and deep inelastic scattering. Unlike CC events, there is no charged lepton in final state. NC events can be detected by finding hadrons or their decay particles in final state. Figure 2.3 shows Feynman diagrams of NC elastic and π 0 production via resonance. Around 1 GeV region of neutrino energy, the major fractions of interactions are elastic and pion production. Figure 2.3: Feynman diagrams of NC elastic and NC π 0 production reaction. 2.2 Neutrino-Nucleon Neutral Current Elastic Scattering Cross-Section The neutral weak current of the nucleon is written as J µ = < N(p ) F 1 (Q 2 )γ µ + F 2 (Q 2 )σ µν q ν + G A (Q 2 )γ µ γ 5 N(p) > (2.2) where F 1 (Q 2 ),F 2 (Q 2 ), and G A (Q 2 ) are the nucleon form factors. The first and second terms are the vector part and the last term is the axial vector part. The differential cross-section as a function of Q 2 is written as dσ dq 2 = 1 J 2 (2.3) 64πEνM 2 N 2 18

21 where the elements of the matrix J is the neutral weak current of the nucleon from Equation 2.2. The differential cross-section using Llewellyn-Smith formalism [36] is written dσ dq 2 = M2 NG 2 F 8πE 2 ν [ A(Q 2 ) ± B(Q 2 ) s u ] + C(Q 2 (s u)2 ) MN 2 MN 4 (2.4) where the G F is the Fermi constant, s u = 4M N E ν Q 2. ± in this equation is + for neutrinos and for antineutrinos. A(Q 2 ),B(Q 2 ) and C(Q 2 ) contain the form factors. A(Q 2 ) = B(Q 2 ) = C(Q 2 ) = 1 4 Q2 M 2 N +F 2 2 Q2 M 2 N [ G 2 A (1 + Q2 (1 Q2 4M 2 N 4MN 2 ) Q 2 ) F1 (1 2 Q2 4M 2 N + F 1 F 2 Q 2 M 2 N 4MN 2 ] ) (2.5) G A (F 1 + F 2 ) (2.6) [ G 2 A + F F 2 2 Q 2 4M 2 N At low Q 2, C(Q 2 ) term is dominant. Q 2 can be written using recoil proton kinetic energy as ] (2.7) Q 2 = 2M p T p (2.8) where T p is the kinetic energy of recoil proton. In SciBooNE experiment, proton track is distinguishable from other charged particles and kinetic energy of proton is reconstructed using track length. The incident neutrino energy is determined as E ν = M N cosθ p (1 + 2M p /T p ) 1/2 1 (2.9) where θ p is track angle of the recoil proton with respect to the beam direction. Q 2 and E ν can be determined without measuring incident energy nor scattering angle of neutrino Strangeness Component of Nucleon Spin The Q 2 dependence of the form factors is parameterized using a dipole form. The Q 2 dependence of the axial form factor is given as G A (Q 2 ) = 1 G A (0) 2 (1 + Q 2 /MA) γ G s A(Q 2 ) (2.10) 19

22 Figure 2.4: Kinematics of νp νp. where M A is the dipole cutoff mass, γ 3 is +1 for proton and 1 for neutron. G A (0) is determined from beta decay. When we consider strange quark contributions to the neutral weak nucleon current, additional form factors need to be introduced. Parameterized form factors are F1(Q s 2 ) = 1 γsq 2 2 (2.11) 6 (1 + Q 2 /M1) 2 2 µ s F2(Q s 2 ) = (2.12) (1 + Q 2 /M2) 2 2 G s A(Q 2 ) = 1 s (2.13) 2 (1 + Q 2 /MA) 2 2 where M 1 and M 2 are the relevant masses of the strange vector form factors, γ s is the strange radius of the nucleon (analogous to the nucleon charge radius), µ s is the strange anomalous magnetic moment of the nucleon, and s is the strange component of nucleon spin. Measurement of s using electron beams s can be measured using electron-nucleon deep inelastic scattering (DIS). Figure 2.5 shows a Feynman diagram of electron-proton deep inelastic scattering. Q 2 of this interaction is written as Q 2 = q 2 (2.14) = (k k ) 2 (2.15) 20

23 = 4EE sin 2θ 2 (2.16) where k = (E, k ) and k = (E, k ) are the momentum of electron before and after scattering, P = (M, 0) is the momentum of target proton, and θ is the scattering angle. Invariant mass W is W 2 = (P + q) 2 (2.17) = M 2 + 2Pq + q 2 (2.18) = M 2 + 2Mν Q 2 (2.19) where ν is the energy transfered from electron to proton by photon. ν is written ν = Pq M Bjorken scaling parameter x is defined as (2.20) x = Q2 2Pq. (2.21) For DIS, the value of x is 0 < x < 1. Spatial resolution increases as Q 2 become larger and electrons are scattered by point-like particles in the nucleon as was confirmed by past experiments [9]. This point-like particles were called partons that are known as quarks today. x represents a fraction of the momentum of a quark inside nucleon at the limit of Q 2. In DIS interaction, a single quark directly interacts with a virtual photon. It does not interacts with any other quarks. Probability of electron-quark interaction with x is represented by quark distribution function. Sum of quark distribution functions of strange quark s(x) and its antiparticle s(x) is written as S(x) = s(x) + s(x). Strange quark contribution to the nucleon spin is obtained by performing integration of S(x) from x = 0 to x = 1. Measurement of s using neutrino beams s can be measured using νn νn interaction. This interaction is neutral current elastic scattering. Neutrino-nucleon and antineutrino-nucleon elastic scattering cross-section were measured by BNL E734 experiment [1]. The E734 experiment used the neutrino beam with the peak energy of 0.6 GeV, the mean energy of 1.2 GeV. 79% of protons in the detector is bound in carbon and aluminum. 21% of protons are free. Figure2.6 and Figure2.7 shows the ν µ and ν µ beam flux of E734 experiment. 21

24 Figure 2.5: Feynman diagram of electron-proton deep inelastic scattering. 22

25 Figure 2.6: The ν µ beam flux of E734 experiment. Dots show the measured ν µ flux. The solid curve is a MC beam flux prediction. Figure 2.7: The ν µ beam flux of E734 experiment. Dots show the measured ν µ flux. The solid curve is a MC beam flux prediction. 23

26 Fit G s A(0) F1(0) s F2(0) s M A χ/n DOF ± / ± ± / ± ± ± ± / ± ± ± ± /11 Table 2.1: The values of the fitting parameters G s A(0), F s 1(0) and F s 2(0) and M A. Re-analysis of the E734 data was performed [10]. Figure 2.8 (a) shows the differential cross-section for νp νp (upper plots) and νp νp (lower plots) measured by E734 experiment. Figure 2.8 (b) shows the ratio of the cross-section for νp νp to the cross-section for νp νp. Solid lines and dashed lines show the result of fit. Error bars include statistic and Q 2 dependent systematic error. From this data, form factors G s A(0), F s 1(0) and F s 2(0) are determined by fitting. The χ 2 is defined as χ 2 = ν (σ f i σ d i /s ν ) 2 δσ 2 i + ν (σ f i σ d i /s ν ) δσ 2 i + (1 s νs ν ) (1 s ν/s ν ) (M A 1.061) (2.22) where ν and ν are summations over all measurement points. σ f i is the differential cross section calculated from the form factors considered as parameters, σ d i is the measured differential cross section. s ν and s ν are scale factors of neutrino and antineutrino data. δσ i is statistical and Q 2 dependent systematic error, and are Q 2 independent systematic error ±0.026 GeV is axial dipole mass. Fitting was performed changing G s A(0),F s 1(0) and F s 2(0). The result of the fitting is shown in Table 2.1. Solid line in Figure 2.8 is the result of fit 1 and dashed line is the result of fit 4. From this fitting, s turned out to have a value from 0.21 to 0. This experiment showed that s can be measured using neutrino-nucleon NC elastic interaction. However, more experiments with new techniques are needed for more precise measurements since the E734 expeirment had large systematic error. 24

27 Figure 2.8: (a) Differential cross-section for νp νp (upper plots) and νp νp (lower plots) measured by E734 experiment. (b) Ratio of the crosssection for νp νp to the cross-section for νp νp. Solid lines and dashed lines show the result of fit. Error bars include statistic and Q 2 dependent systematic error. [1][10] 25

28

29 Chapter 3 Experimental Setup of SciBooNE 3.1 Experimental Setup Overview SciBooNE uses the FNAL Booster Neutrino Beam (BNB) with a peak energy of 0.7 GeV. Protons with 8 GeV from the Booster hit a beryllium target and produces mesons, i.e. pions and kaons. The target is surrounded by a magnetic focusing horn. The produced mesons are focussed, and decay inside a 50 m long decay volume: π ± µ ± + ν(ν), K ± µ ± + ν(ν). The polarity of the magnetic focusing horn can be changed to select the neutrino or antineutrino mode which enables us to study both neutrino and antineutrino cross-sections. SciBar, an electro-magnetic calorimeter (EC) [16], and a muon range detector (MRD) [17] are detectors for SciBooNE. They are located 100 m downstream of the proton target. MiniBooNE detector is located 440 m downstream of the SciBooNE detector. Experimental setup of SciBooNE is shown in Figure 3.1. Air view of the Booster, target station, decay region and SciBooNE detector are shown in Figure Booster Neutrino Beam Proton Beam 8 GeV protons from the Booster as the primary beam are provided to the target station. The proton beam has approximately protons per 1.6 micro second beam batch at a maximum rate of 5 Hz. The timing structure within each beam batch has 81 bunches, each bunch is approximately 6 ns wide, separated by 19 ns. 27

30 Figure 3.1: Experimental setup of SciBooNE. SciBooNE detector is located 100 m downstream of the proton target. MiniBooNE detector is located 440 m downstream of the SciBooNE detector. Trajectories and positions of the proton beam are measured on a pulseby-pulse basis. The beam position monitors located near the target measure the beam alignment and divergence. The beam position and angle are automatically adjusted to keep the proton beam centered on the target. The beam profile monitors measure beam focusing on the target. The number of protons delivered to the target is measured for each proton batch with two toroids located near the target along the beam line Horn and Target The proton target is located in a target station. The target is 71 cm long, 1 cm diameter beryllium. Primary protons hit the target and secondary mesons (pions and kaons) are produced. The target is enclosed in the magnetic focusing horn. The horn makes coaxial magnetic field in the direction of the beam to focus the secondary mesons produced on the target. The horn current pulse is approximately 174 ka, 143µs long, synchronized to each beam spill. The horn can change the polarity of magnetic field to focus positive charged particles and defocus negative charged particles, and vice versa. The secondary mesons from the target decay into neutrinos in the decay region. The decay region is 50 m long and its radius is 90 cm. The cylindrical volume is filled with air at atmospheric pressure. A beam absorber is located at the end of the decay region. It stops the hadronic and muonic component of the beam. 28

31 Figure 3.2: Air view of the Booster, target station, decay region and Sci- BooNE detector. 29

32 3.3 Detectors SciBar, an electro-magnetic calorimeter (EC) [16], and a muon range detector (MRD) [17] are detectors for SciBooNE (Figure 3.3). They are located 100 m downstream of the proton target. A right-handed Cartesian coordinate system is used. z axis is the beam direction and y axis is the vertical upward direction. The origin of the coordinate system is defined as the center of SciBar in xy plane and the most upstream edge of SciBar in z direction. Figure 3.3: SciBar, an electro-magnetic calorimeter (EC), and a muon range detector (MRD) are detectors for SciBooNE. They are located 100 m downstream of the proton target SciBar Detector The SciBar detector [15] is a fully-active neutrino detector with finely segmented structure. Carbon and hydrogen in the plastic scintillators are the target nuclei. Its volume is m 3 and weight is 15 tons. The SciBar detector consists of 14,336 plastic scintillator bars (Figure 3.4). Dimensions 30

33 of each plastic scintillator bar are cm 3. The SciBar detector consists of 64 layers arranged perpendicular to the beam direction. Each layer has a vertical plane and a horizontal plane. One plane consists of 112 bars of plastic scintillator arranged vertically or horizontally. The scintillator bars are made of polystyrene, infused with PPO (1%) and POPOP (0.03 %) with 0.25 mm thickness TiO 2 coating. The emission wave length is 420 nm at peak position. One wave length shifting (WLS) fiber is inserted in each plastic scintillator bar. The diameter of the WLS fiber is 1.5 mm. The attenuation length of each WLS fiber is measured, and is about 350 cm on average. The scintillation light from the plastic scintillators is sent to a 64-channel multi-anode photo-multiplier tube (MAPMT) through the WLS fibers. In total, 224 of MAPMT (H8804 made by Hamamatsu Photonics K.K.) are used in the SciBar detector. Each has 64 anodes with 2 2 mm 2 area. The typical gain is with a linear response up to 200 photoelectrons (p.e.); the gain uniformity over all channels is 20% in RMS. The fibers and the MAPMT are connected using a dedicated jig, called a cookie (Figure3.5) to align the fibers to the anodes. There is cross-talk mainly caused by light entering adjacent channels on the surface of the MAPMT. The average level of the cross-talk from an adjacent channel is measured to be 3.5%. The gains of all the SciBar channels are monitored using a gain monitoring system which consists of LEDs, clear fibers, and light injection modules. The light from four LEDs is sent via clear fibers to each light injection module located between the plastic scintillators and every MAPMT. The light injection modules are designed to uniformly distribute light from the clear fibers to each WLS fiber. PIN-photo diodes are located near each LED to monitor the LED luminosity. The gain and pedestal data for each MAPMT channel are taken between every beam spill to monitor the detector stability. Light yield and timing resolution of the SciBar detector are evaluated using cosmic-ray muons. Light yield is 20 p.e./1.3 cm for minimum ionizing particles at the detector edge close to the MAPMT. Timing resolution is 1.6 ns. The number of dead channels of the SciBar detector is one out of 14,336. Each MAPMT is mounted on a front-end board as shown in Figure 3.5. Each front-end board has two sets of Application Specific Integrated Circuit (ASIC), called VA/TA. The VA is an ASIC with a preamplifier, a shaper and a multiplexer which serialize charge from 32 channels of the MAPMT to a single output. The TA is an ASIC which provides inclusive timing information for the 32 channels. The front-end boards process signals from the MAPMTs, and send them to VME-9U backend electronics boards. They process the charge and timing information from the MAPMT with 12-bit flash ADCs and 64-ch multi-hit TDCs. Charged particles whose track length exceeds 8 cm inside the detector are 31

34 Figure 3.4: The SciBar detector consists of 14,336 of the combination of plastic scintillator bar and wave-length shifting (WLS) fiber. Volume of the detector is m 3, weight is 15 tons.the Scintillator bars are arranged vertically and horizontally. Each WLS fiber is attached to one anode of a multi-anode PMT. 32

35 reconstructed with 99% efficiency. Figure 3.5: The WLS fibers are attached to the MAPMT by the cookie. The light injection module is located between the scintillators and the MAPMT. Each MAPMT is mounted on the front-end board Electromagnetic Calorimeter The Electromagnetic Calorimeter called Electron Catcher (EC) is located downstream of SciBar. The EC is designed to measure the electron neutrino (ν e ) contamination in the beam and to tag γ s from π 0 decay. The EC consists of 1 mm diameter scintillating fibers embedded in lead foil (Figure 3.6). The calorimeter is made of modules of dimensions cm 3. Each module is read out by two 1 inch Hamamatsu PMTs per side. Total number of PMTs is 256. The EC modules were built for the CHORUS neutrino experiment at CERN [18] and later used in HARP and then K2K. The modules construct one vertical and one horizontal plane, and each plane has 32 modules. The EC has a thickness of 11 radiation length along the beam direction. An active area of the plane is m 2. A minimum ionizing particle with a minimal path length deposits approximately 85 MeV. The energy resolution for electrons was measured to be 14% / E(GeV) using a test beam. 33

36 Figure 3.6: The EC consists of 1 mm diameter scintillating fibers embedded in lead foil. The calorimeter is made of modules of dimensions cm 3. Each module is read out by two 1 inch Hamamatsu PMTs per side Muon Range Detector The MRD (figure 3.7) is installed downstream of the EC. The MRD is designed to measure the momentum of muons. It has 12 iron plates with thickness of 5 cm and 13 plastic scintillator planes with thickness of 6 mm. Iron plates are sandwiched between scintillator planes. Scintillator counters are aligned altering horizontal and vertical planes with 2 inch PMTs. Total number of PMTs are 362. Iron plates cover areas of cm 2. The MRD measures the momentum of muons up to 1.2 GeV/c using the observed muon range. Hit finding efficiency was continuously monitored using cosmic ray data taken between beam spills. The average hit finding efficiency is 99% Triggers There are two triggers, for neutrino data taking (beam trigger) and for detector calibration data (off-beam trigger). They are taken in every beam cycle. One cycle is about 2 seconds. Beam trigger is a fast timing signal from the extraction magnet on BNB. Off-beam trigger is set after the beam trigger. Each detector takes pedestal and cosmic ray data. SciBar and the EC use a common cosmic ray trigger generated using fast signals from the TA. The MRD has independent self-generated cosmic ray trigger. SciBar also takes LED data for monitoring relative gain of all channels. Events taken by beam trigger and off-beam trigger is checked by shifter once per hour using online event display. Neutrino events can be seen occa- 34

37 Figure 3.7: The MRD is installed downstream of the EC. It has 12 iron plates with thickness of 5 cm and 13 plastic scintillator planes with thickness of 6 mm. 35

38 sionally on online monitor in the beam trigger, but most of the event display of the beam trigger have only pedestal (no hits). Figure 3.8 is a neutrino event candidate on online display. Vertical (left) and horizontal (right) readout of detectors are shown. Figure 3.8: A neutrino event candidate on online display. Vertical (left) and horizontal (right) readout of detectors are shown. 3.4 Neutrino Events in the SciBar Detector Figures 3.9 and 3.10 show charged current quasi elastic (CCQE) event candidates from actual data. The area of dots is propotional to the ADC channel (i.e. energy deposit). Boxes represent TDC hit information by their colors. The track penetrating the MRD in Figure 3.9 is identified as muon. One of the tracks in Figure 3.10 is identified as muon which decayed inside the SciBar using the multi-hit TDC information. Existence of muon is the sign of charged current events. Difference in de/dx between muons and protons can be seen. de/dx distributions of muons and protons are shown in Figure Protons are separated from muons with 90% efficiency. CC 1π production can be categorized into three different groups: nuclear coherent interaction νa µaπ, neutron scattering νn µ nπ + and proton scattering νp µ pπ +. As shown in Figure 3.12, scattering on neutrons or coherent interactions appears as two charged tracks with a common vertex in the SciBar, while scattering on proton appears as three charged tracks as in Figure If one of the tracks is not long enough to be tracked in the proton interaction event, i.e. as is the case for the proton track in Figure 36

39 Figure 3.9: Event displays of CCQE candidate. The area of dots is proportional to ADC channels. Boxes represent TDC hit information. Tracks penetrating the MRD are identified as a muon. Figure 3.10: Event displays of CCQE candidate. The decay of the muon is tagged by multiple TDC hits. Figure 3.11: de/dx distributions of muons (left plot) and protons (right plot). Protons are separated from muons with 90% efficiency. 37

40 3.13, the event can be distinguished from the neutron or coherent scattering events by looking at large energy deposit near the vertex. Figure 3.12: CC 1π coherent or νn µ nπ + candidate. Muon and pion are visible. Figure 3.13: CC 1π νp µ pπ + candidate. Muon, pion and proton are visible. There is no muon in the final state of neutral current events. Figure 3.14 shows a candidate of neutral current elastic scattering. The signature of this event is a single proton track. The proton is identified using de/dx information. Figure 3.15 shows a candidate of neutral current π 0 production event. The 2 γ s from a π 0 decay produce electro-magnetic showers which can be detected in SciBar. These features of SciBar enable us to make precise measurements of neutrino-nucleus interactions in the region of neutrino energy below 1 GeV. 3.5 Data Taking We started data taking in June 2007 with antineutrino mode, switched to neutrino mode in October We completed taking data of projected protons on target (POT) for neutrino mode in April 2008, then switched back to antineutrino mode. We completed taking data of antineutrino mode in August So, the total beam time was about 13 months. 38

41 Figure 3.14: Candidate of neutral current elastic scattering. Signature of this event is single proton track. The proton is identified using de/dx information. Figure 3.15: Candidate of neutral current π 0 production event. The 2 γ s from a π 0 decay produce electro-magnetic showers, which can be detected in the SciBar. Accumulated number of POT and number of POT which passed all data quality cuts are shown in Figure POT were delivered all together. Beam quality cuts are applied to all data. Beam intensity is required to be at least protons per spill. The agreement between the readout of two troids need to be within 10%. The absolute peak horn current is required to be greater than 170 ka. The targeting efficiency is required to be greater than 95%. The efficiency of all the beam quality cuts is greater than 99%. After applying all (beam and detector) data quality cuts, data of POT are selected for physics analysis. As for neutrino data, are selected. The rest is antineutrino data. The number of CC event candidates is evaluated from tracks starting inside the SciBar every week during the data taking period. With the obtained number, one can monitor the SciBar status. Figure 3.17 shows number of CC event candidate normalized by POT. The flat shape for each mode and same values for the two antineutrino mode periods shows that the SciBar took data stably since we started data taking. Neutrino event rate is 4.5 times larger than antineutrino event rate because of the difference between cross-sections of neutrino-nucleus and antineutrino-nucleus scattering, and production cross-sections of π + and π that are parent particles of neutrino and antineutrino beam. 39

42 Figure 3.16: Accumulated number of POT (solid line) and number of POT which passed all the data quality cuts (dotted line). 40

43 Figure 3.17: Number of CC event candidates normalized by POT. The flat shape for each mode and the same values for the two antineutrino mode periods shows that the SciBar took data stably throughout the beam time. 41

44 3.6 Beam Simulation Neutrino flux at the SciBooNE is predicted using GEANT4-based [19] beam Monte Carlo simulation. SciBooNE uses the same simulation code developed by the MiniBooNE Collaboration [20] at FNAL. The primary protons are generated according to the beamline optics upstream of the target. The interactions of primary protons with the beryllium target are simulated according to state-of-the-art hadron interaction data. The geometry and materials of all objects in the target hall and decay region are modeled in the simulation code. Predictions of the secondary particles production on beryllium target are important for neutrino flux predictions. Experimental data of hadron interaction from HARP [21] and BNL E910 [22] are used. Production of secondary protons, neutrons, charged pions, and charged and neutral kaons is taken into account. Elastic and quasi-elastic scatterings of protons in the target are also simulated. Particles from the primary proton-beryllium target interactions are simulated by the GEANT4 framework taking into account all the physics process. Hadronic re-interactions of pions and nucleons with beryllium target and aluminum materials are particularly important. They are described by recently updated models. A Sanford-Wang parameterization [23] describes the double-differential inclusive production cross-section for pions and neutral kaons. The Sanford-Wang parameterization of the inclusive production cross section is given by d 2 σ dpdω = c p p c4 1(1 )exp[ c 3 p B c 9 p c c 5 6 θ(p c 7 p B cos c 8 θ)] (3.1) B where p and θ are the meson momentum and angle, p B is the beam momentum, and the constants c 1,c 2,..., c 9 are parameters determined from fits to meson production data from HARP and BNL E910. The Feynman scaling parameterization describes the K + production cross section [24]. The Feynman scaling parameterization is given by d 2 σ dpdω = P 2 K E K c 1 (1 x F ) c 3 exp[ c 3 x F c 4 (3.2) c 7 p t x F c 6 c 2 p t c 5 p 2 t] where x F = p CM /p CM,max is the Feynman scaling parameter (p is longitudinal momentum), p t is the transverse momentum of the K + and the constants c 1,c 2,..., c 7 are parameters determined from fits to experimental data [25][26][27][28][29][30][31][32][33]. The other hadronic processes and all electromagnetic processes are described by default GEANT4 physics lists. 42

45 To simulate the neutrinos from the secondary particles (mesons), the output of the GEANT4 codes is used by the second MC code as input. The second MC generates the neutrino kinematics distributions from meson and muon decays and obtains the final neutrino fluxes extrapolated to the SciBooNE detector. Current best knowledge of neutrino-producing meson and muon decay branching fractions, and decay form factors in three-body semi-leptonic decays are used. Polarization effects in muon decays are also accounted for. All neutrinos that enter the detectors are considered for SciBooNE flux predictions. Weights for each neutrino event is calculated using information of neutrino flavor, energy, parent type and kinematics. Information from the interaction and detector simulation such as neutrino interaction probability, detailed detector geometry and specifications are also used. Neutrino flux prediction at the SciBooNE detector as a function of neutrino energy E ν is shown in figure A total neutrino flux per proton on target of cm 2 is expected at the SciBooNE detector location in neutrino mode. Mean neutrino energy is predicted to be 0.7 GeV. 93% of total flux is muon neutrinos. Fraction of muon antineutrinos is 6.4% and electron neutrinos and antineutrinos is 0.6%. 43

46 Figure 3.18: Neutrino flux prediction at the SciBooNE detector as a function of neutrino energy E ν. 44

47 3.7 Neutrino Interaction Simulation NEUT [34][35] is a program for neutrino interaction simulation. Elastic scattering, single meson production, single gamma production, coherent pion production and deep inelastic scattering in both neutral and charged currents are simulated by NEUT. Nuclear effects are also considered. NEUT handles protons, oxygen, carbon and iron as nuclear targets. The energy range for the simulation is from 100 MeV to 100 TeV. Output from the beam simulation is used as input for the NEUT Elastic Scattering Elastic scattering is implemented in NEUT using the model of Llewellyn- Smith [36]. Fermi gas model of Smith and Moniz [37] is used for scattering off nucleons in the nucleus. The Fermi motion of nucleons along with the Pauli exclusion principle is taken into account. The momentum distribution of the target nucleon is assumed to be flat up to a fixed Fermi surface momentum of 217 MeV/c for carbon and 250 MeV/c for iron. The same Fermi momentum distribution is also used for all other nuclear interactions. The nuclear potential is set to be 27 MeV for carbon and 32 MeV for iron. Both vector and axial-vector form factor are assumed to be dipole. The vector mass in elastic scattering is set to be 0.84 GeV/c 2. The axial vector mass is set to be 1.11 GeV/c 2. The Q 2 dependence of the differential cross-section for NC elastic scattering in NEUT is shown in Figure Black line and red line shows differential cross-sections for neutrino-proton scattering and neutrino-neutron scattering respectively Resonant Single Meson Production Resonant production of pions, kaons and etas via baryon resonances are described by the model of Rein and Sehgal [38]. The model assumes an intermediate baryon resonance, N, in the reaction of νn ln, N N m. The differential cross-section for single meson production depends on the amplitude for the production of a given resonance. The probability of the baryon resonances with mass less than 2 GeV/c 2 is included. Those baryon resonances with mass greater than 2 GeV/c 2 are simulated as deep inelastic scattering. Lepton mass effects from the non-conservation of lepton current and the pion-pole term in the hadronic axial vector current are included in the simulation [39][40]. To determine the angular distribution of a pion in the final state, Rein s method [41] is used for the P 33 (1232) resonance. For other resonances, the directional distribution of the generated pion is set to 45

48 Figure 3.19: The Q 2 dependence of the differential cross-section for NC elastic scattering adopted in NEUT. Black line and red line shows differential cross-sections for neutrino-proton scattering and neutrino-neutron scattering respectively. 46

49 be isotropic in the resonance rest frame. The angular distribution of π + has been measured for ν µ p µ pπ + [42] and the results agree well with NEUT s prediction. Pauli blocking is accounted for in the decay of the baryon resonance by requiring the momentum of the nucleon to be larger than the Fermi surface momentum. Pion-less decay is also taken into account, where 20% of the events do not have a pion [43] Coherent Pion Production Coherent pion production is a neutrino interaction with a nucleus which remains intact, creating one pion with the same charge as the incoming weak current. Rein and sehgal model [44] with lepton mass correction [46] and Kartavtsev et al. model [47] are used to simulate the coherent pion production Nuclear Effects in Final State NEUT simulates the interactions of the generated mesons and nucleons in nucleus. Nucleons which have obtained recoil energy in nucleus interacts with other nucleons bound in the same nucleus. The nucleon-nucleon elastic scattering cross section is based on the measurements by Bertini et al. [48]. The pion production from the decay of produced is also taken into account according to the isobar production model [49]. Probability of nucleon-nucleon interaction in 16 O after the neutrino interaction is shown in figure3.20. No interaction, elastic scattering, single π production, and two π production are plotted. The production of nucleons with nucleon momentum below 225 MeV/c is suppressed by the Pauli blocking effect. Momentum above 300 MeV/c, about a half of the nucleons interact mainly by elastic scattering. Mean deflection angle is shown in figure3.21. It has a peak around 45 degrees. The concept of formation zone is also considered for all hadrons generated in nucleus. This is a distance or time from the neutrino interaction point to the hadron production point. The intermediate states are assumed to be non-bound quark states. The formation length of each hadron is expressed as L = p/µ 2 (3.3) where p is the momentum of the hadron and µ 2 =0.08 GeV 2 is a fitted constant from the SKAT experiment [50]. We can see the effect of nucleon-nucleon 47

50 Figure 3.20: Probability of nucleon-nucleon interaction in 16 O after the neutrino interaction. No interaction (solid line), elastic scattering (dashed line), single π production (dot-dot line), and two π production (dash-dod line) are plotted. 48

51 Figure 3.21: Mean deflection angle. It has a peak around 45 degrees. 49

52 interactions in the nucleus in the significant distortion of the geometry of particles observed in detector. 3.8 Detector Simulation The GEANT4 framework is used for the detector simulation. Output from neutrino interaction simulation is used as input for the detector simulation. Geometries of the detectors, the detector frame and experimental hall are based on survey measurements which were carried out during detector construction period. Particle Simulation in the Detectors For the hadronic interactions in the detectors, the Bertini cascade model in GEANT4 [51] is used. The details of the nuclear model are discribed in reference [52][53]. The motion of the bound nucleons, the Pauli exclusion principle, the diffuseness of the nuclear edge, and a local potential for nucleons and pions are included in the model. Data from several free-particle experiments are used for the crosssections for nucleon-nucleon and nucleon-pion reactions in the Bertini cascade model [54][55][56][57][58][59][60][61][62][63]. Figure 3.22 shows proton-proton total and elastic cross sections. Figure 3.23 shows the same plot for protonneutron. Where these data were lacking, they were estimated by use of isospin symmetry or phase-shift analysis. Detector Response For the SciBar detector simulation, information of particles generated by NEUT is used as input to the simulation. The energy loss of a charged particle in a single strip of SciBar is simulated by GEANT. This energy scale is tuned using cosmic ray data. Scintillator quenching is simulated using Birk s law [64] with a value of Birk s constant of cm/mev measured at K2K [65]. The energy deposited by a charged particle is converted to photoelectrons using conversion factors. Conversion factors are measured for each channel of SciBar with cosmic muons. The measured light attenuation length of each fiber is used in the simulation. Average light attenuation length is approximately 350 cm. Cross-talk between nearby MAPMT channels is measured in laboratory and used in the simulation. The number of photoelectrons is smeared by Poisson statistics. The single photoelectron resolution of the MAPMT is simulated. The number of photoelectrons is 50

53 Figure 3.22: Proton-proton total and elastic cross sections. Figure 3.23: Proton-neutron total and elastic cross sections. 51

54 converted to ADC channels, and then electronics noise and threshold effects of the TA are simulated. TDC hit simulation includes light propagation delays in the WLS fibers. A logical OR of 32 MAPMT channels is made for each TDC channel, and the time of each hit is converted to TDC channels. Multiple TDC hits in each channel are simulated. In the EC detector simulation, true energy deposition in scintillating fibers in the detector is converted to the number of photoelectrons using the conversion factor. The conversion factor is measured for each channel with cosmic muons. The attenuation of light in the fiber is simulated using the measured attenuation light value. The number of photoelectrons is smeared by Poisson statistics and by the PMT resolution, and then converted to ADC counts. The time-dependent ADC gain due to the overshoot of the PMT signal is simulated based on a measurement with cosmic muons. Electronics noise is also simulated. For the MRD detector simulation, true energy deposition in each scintillator is converted to ADC channels using the conversion factor measured with cosmic muons. The attenuation of light in the scintillator and electronics noise are simulated. Gaps between scintillator counters in each plane are included in the simulation. The time of energy deposition is digitized and converted into TDC channels. 3.9 Dirt Simulation Recoiled particles recoiled by neutrinos inside dirt or wall of experimental hall around the detector sometimes enter the detector and make hits or interact with particles inside the detector. These events are called dirt events. Dirt events mimic neutrino interactions in SciBar and can become backgrounds. For NC elastic analysis, the major background is protons recoiled by neutrons from the dirt (dirt neutron). Since both the signals and dirt neutron events are single proton track, it is difficult to distinguish the signals from dirt neutron events. To estimate backgrounds from dirt events, they are simulated using the beam simulation, NEUT and the detector simulation. The beam simulation is extended to cover a volume of m 3, taking the detector origin as a center (figure 3.24). Gray areas indicate the dirt volume. For NEUT, carbon is used for the material of the dirt with density 2.15 g/cm 3. For detector simulation, default composition of concrete which is available in GEANT4 is used for the material of the dirt. Density is 2.15 g/cm 3 based on a survey measurement. In this simulation, same material and density are used for dirt and walls of the experimental hall. 52

55 Figure 3.24: Geometry of dirt for the dirt simulation. Gray areas indicate dirt volume. The geometry of m 3, taking detector origin (red point) as a center were considered. 53

56 Figure 3.25: Vertex distributions of dirt events from simulation. Only events which made one track contained in SciBar were selected. Geometry is large enough for NC elastic analysis. 54

57 Chapter 4 Analysis of Neutral Current Elastic Scattering Figure 4.1 shows a candidate of neutral current elastic scattering from Sci- BooNE data. Signature of this event is single proton track. The proton is identified using the stopping range and de/dx informaiton. Analysis for evaluating the cross section for neutrino-nucleon neutral current elastic scattering is performed. Flow chart of this analysis is shown in Figure Track Reconstruction Tracks in SciBar are reconstructed as two-dimentional (2D) tracks in each view of SciBar using a cellular automaton algorithm [66]. Three dimentional (3D) tracks are reconstructed by matching the timing of 2D tracks and z position of track edges. The conditions for matching requires the timing difference between two 2D tracks to be less than 50 ns and the z-edge difference to be less than 6.6 cm for both upstream and downstream edges. The hit threshold for tracking is set to two photoelectrons. Reconstructed tracks are required to penetrate at least three layers. Due to this requirement, the minimum length of a reconstructed track is 8 cm. 4.2 Monte Carlo Normalization All MC events used in this analysis is normalized with a data over MC ratio of number of CC events N cc /NMC. cc CC events are identified by matching a track in SciBar and a track or hits in the MRD. The track which satisfied this condition is called SciBar-MRD matched track. The most energetic 55

58 Figure 4.1: A candidate of neutral current elastic scattering from SciBooNE data. Signature of this event is single proton track. The proton is identified using stopping range and de/dx information. 56

59 Figure 4.2: Flow chart of this analysis. 57

60 SciBar-MRD matched track is considered as a muon candidate. A MRD track is required to have the upstream track edge on the first two layers of the MRD. The transverse distance between the two tracks at the first layer of the MRD is required to be less than 30 cm. The requirement on the difference between track angles with respect to the beam direction is given by θ MRD θ SciBar < θ MAX where θ MAX is a function of the length of the MRD track. θ MAX varies between 0.4 to 1.1 radians. At least two layers with hits in each view are required for track reconstruction in the MRD. If no MRD track is found, the SciBar track is extrapolated to the MRD and search for nearby contiguous hits in the MRD identifying a short muon track. For matching MRD hits to a SciBar track, the MRD hit is required to be within a cone with an aperture of ±0.5 radian and a transverse offset within 10 cm of the extrapolated SciBar track at the upstream edge of the MRD. The timing difference between the SciBar track and the track or hits in the MRD is required to be within 100 ns. Lower limit of the momentum of muon with matched track is 350 MeV/c due to these matching criteria. 4.3 Event Selection Event selections using geometrical information and particle identification are performed to pick up NC elastic events in fiducial volume (FV). FV is defined from 3rd layer to 62th layer, ±130 cm from detector center along the axis perpendicular to the beam direction Geometrical Selection Veto Outside the FV of the SciBar is used for veto to avoid charged particles entering the SciBar from outside and escaping from inside the SciBar. Veto uses a clustering method. Definition of cluster is more than two hits with photoelectrons larger than 10 in a distance less than two cells. If there is any cluster in veto region and time difference between the track inside the FV and the cluster is less than 100 ns, event is rejected. This method automatically selects tracks contained in FV. Tracks need to be contained so that kinetic energy of recoil proton is reconstructed correctly. EC is also used for veto. The purpose of the EC veto is to reject NC π 0 events in which 2γ from π 0 decay escape from SciBar without making electromagnetic shower. If there is a coincidence of a pair of PMTs of at least one EC module, event is rejected. It reduces 11.9% of backgrounds from NC pion production. 58

61 One Track Event Selection For NC elastic sample, only events with one track are selected. 2D track in both view are required to be single track. In case tracks do not penetrate three layers because of large angle or short length, they are not reconstructed. To detect such tracks, two additional methods are used. For large angle tracks, more than five continuous hits with 5 photoelectrons in the same layer are considered as one track. Time difference between reconstructed track and all hits are required to be less than 100 ns. Searching area of large angle track is from two layers upstream to two layers downstream of the vertex. For short tracks and other tracks that are not reconstructed, total energy deposit from track-unrelated hits is used. Total energy deposit is calculated for each view. Hits with more than 10 photoelectrons are used for energy deposit calculation to pick up minimum ionizing particles, considering light attenuation by WLS fibers. Figure shows total energy deposit of track-unrelated hits on top and side of SciBar. The threshold to select one track events is set to 30 MeV for both top and side view. Figure 4.3: Total energy deposit of track-unrelated hits on top and side of SciBar. Data are compared to MC simulation. The threshold to select one track events is set to 30 MeV for both top and side view. In these plots, dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events. 59

62 4.3.2 Particle Identification To identify protons, two types of process were performed: decay particle rejection and de/dx cut. Decay particle rejection uses multi-hit TDC information. de/dx cut uses energy deposit along the track and track length of reconstructed tracks. Decay Particle Rejection Muons and charged pions stopped inside SciBar FV is identified using multiple TDC hit. If a muon decays inside SciBar, there are two TDC hits: one from the muon, the other from the electron from decay muon. TDC information in areas of channels around upstream and downstream track edges is checked for both top and side view. TDC hits within 100 ns from the first hit are ignored due to TDC dead time. There are two criterias for tagging decay particles corresponding to two cases. Case 1: There are multiple TDC hits on both top and side view. To avoid mis-identification by multiple TDC hits from noise, the same time interval between first hit and ith hit are required for both top and side view. If the difference of the time interval on top and side is less than 20 ns, they are identified as hits from decay particles (figure 4.4). Case 2: There are multiple TDC hits only on one view. The time of the ith hit in one view and the time of the first hit in the other view are compared. If the time difference of those hits is less than 20 ns and the first hit which is not followed by second hit have more than 10 photoelectrons, it is also identified as a signal from decay particles (figure 4.5). Muon lifetime is measured using data tagged as muon decay. Figure 4.6 shows the time distribution of muon decay sample. Horizontal axis is the time interval between primary TDC hits and secondary TDC hits that matched on top and side view. From the result of fitting, muon lifetime is 2.09 µs. de/dx Cut Protons have larger energy deposit compared to muons, charged pions and electrons for the same track length. Thus, protons can be identified using de/dx information. de/dx is calculated by total energy deposit along the track divided by the track length. Figure 4.7 shows track length v.s. de/dx. Red points represent protons and black points represent other particles such as muons, charged pions and electrons. Line shows the cut for proton identification. Particles above the line are identified as protons. de/dx distributions for several track length are shown in figure de/dx is calculated from 60

63 Figure 4.4: TDC hits from muon decay in top and side view. There are multiple TDC hits on both top and side view (case 1). Figure 4.5: TDC hits from muon decay in top and side view. There are multiple TDC hits only on one view (case 2). 61

64 Figure 4.6: Time distribution of muon decay sample. Horizontal axis is the time interval between primary TDC hits and secondary TDC hits that matched on top and side view. 62

65 total energy deposit of the track divided by the track length. The peak with large de/dx is from proton track while the peak with small de/dx is from muons, charged pions and electrons. Using this cut, protons are distinguished from other particles with 95% purity and 95% efficiency where efficiency is defined by number of protons above cut divided by total number of protons. Figure 4.7: Track length v.s. de/dx. Red points represent protons and black points represent other particles such as µ ±,π ± and e ±. Line shows the cut for proton identification. Particles above the line are identified as protons. Angular distribution of NC elastic sample is shown in figure 4.9. θ is an angle of the track with respect to the z direction. Length distribution (left figure) and data/mc ratio (right figure) of NC elastic sample are shown in figure There are data excesses in regions of the angle greater than 60 and the track length shorter than 30 cm. 63

66 Figure 4.8: de/dx distributions for several different track lengths. The top left is 0-10 cm, the bottom right is cm. de/dx is calculated from total energy deposit of the track divided by track length. The peak at large de/dx is from proton track while the peak at small de/dx is from other MIP particles. Data are compared to MC simulations. In these plots, dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events. 64

67 Figure 4.9: Angular distribution of NC elastic sample. θ is an angle of the track with respect to the z direction. Data are compared to MC simulation. Dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events. Right plot shows data/mc ratio. 65

68 Figure 4.10: Length distribution (left figure) and data/mc ratio (right figure) of NC elastic sample. Data are compared to MC simulation. Dots with error bars is data and others are MC simulation. Red is NC elastic events, yellow is other interaction, and green is dirt events. Right plot shows data/mc ratio. Proton kinetic energy is reconstructed from track length. To obtain kinetic energy from proton track length, a relationship between them are studied using true proton kinetic energy and true proton track length inside SciBar in MC simulation. Figure 4.11 shows true proton kinetic energy v.s. true proton track length. A line is a result of fitting. Distribution of proton kinetic energy is shown in figure There is a data excess with proton kinetic energy less than 0.2 GeV which corresponds to proton track length of 30 cm and Q 2 of about 0.4 (GeV/c) 2. There is a possibility that tracking algorithm works differently for data and MC when a track is short. Currently the discrepancy between data and MC in the region of proton kinetic energy less than 0.2 GeV is being studied. In this region, there is also a large contribution of dirt events. Data with proton kinetic energy greater than 0.2 GeV are used for further analysis. 66

69 Figure 4.11: True proton kinetic energy v.s. true proton track length from MC simulation. A line is a result of fit. 67

70 Figure 4.12: Distribution of proton kinetic energy. Proton kinetic energy is reconstructed from track length. Data are compared to MC simulation. 68

71 4.4 Background Estimation There are two major backgrounds in the NC elastic sample. One is from NC pion production without electromagnetic shower caused by gammas from decay neutral pion. The other background is from protons recoiled by dirt neutrons. Both of backgrounds are estimated using MC simulation. NC pion production candidates are, in usual case, events with two tracks of electromagnetic shower caused by two gammas from decay neutral pion. NC pion production at SciBooNE is now being studied [67]. In this analysis, NC pion production candidates shows good agreement between data and MC in number of event, geometrical information of the signal and reconstructed mass of neutral pions. Current MC simulation for NC pion production can be used for background estimation. Dirt MC simulation was checked whether tuning is needed or not. There are more dirt events in upstream part of SciBar than downstream part. Figure 4.13 shows distributions of vertex position along z axis v.s. track length of data (left), SciBar MC (middle) and dirt MC simulation (right). SciBar is divided into 8 bins along z direction to compare number of events at upstream part with that of downstream part. Fifth bin is taken as a reference and is compared with first, second and third bins. Figure 4.13: Vertex position along z axis v.s. track length of data (left), Monte-Carlo simulation in SciBar (middle) and dirt Monte-Carlo simulation (right). SciBar is divided into 8 bins along z direction to compare number of events at upstream part with that of downstream part. Fifth bin is taken as a reference and compared with first, second and third bins. Track length distributions of (upstream bins)-(fifth bin) are shown in 69

72 figure Data and sum of SciBar MC and dirt MC are compared. Data agrees with MC simulation within statistic error. Figure 4.14: Track length distributions of (upstream bins)-(fifth bin). Data agree with Monte-Carlo simulation within statistic error. Particles which have roots in dirt events and make NC elastic scatteringlike track in SciBar are shown in Table % of the dirt background is proton track. Fractions of the neutrino interactions in the dirt which make proton track in SciBar are shown in Table4.2. Proton tracks from dirt events are mainly recoiled by neutrons from the dirt. Neutrino-neutron NC elastic scattering (νn νn) in the dirt have the largest contribution to the dirt background. 70

73 particle number of events proton 1814 electron 64 muon 4 pion 19 total 1901 Table 4.1: Particles which have roots in dirt events and make NC elastic scattering-like track in SciBar. interaction fraction (%) NC elastic (νp νp) 3 NC elastic (νn νn) 30 NC pion production 13 CC elastic (νn µ p) 16 CC pion production 27 Table 4.2: Fractions of the neutrino interactions in the dirt which make a proton track in SciBar. 4.5 Event Selection Summary Summary of the NC elastic event selection is shown in Table 4.1. MC (signal) is νp νp and νn νn in MC simulation. MC (B. G.) is a total of non-nc elastic events in SciBar and all dirt events. After event selection, there are 8441 events in data and 7675 events in MC. event selection data MC (signal) MC (B. G.) purity (%) total veto track decay particle rejection de/dx cut Table 4.3: Event selection summary table of the NC elastic event selection. 71

74 4.6 Detector Unfolding Proton kinetic energy distribution after background subtraction is shown in figure Bin width is changed considering the bin width of Q 2 distribution which is reconstructed from proton kinetic energy. Figure 4.15: Proton kinetic energy distribution after background subtraction. Proton kinetic energy is smeared by track length resolution of the detector. Detector unfolding is applied to the NC elastic sample after background subtraction. Figure 4.16 shows distribution of true proton kinetic energy v.s. proton kinetic energy reconstructed from measured track length using tracking algorithm (left figure). Both quantities are from MC simulation. From this relationship, detector unfolding matrix was made (right figure). 72

75 Figure 4.16: Proton kinetic energy reconstructed from track length using Monte-Carlo simulation v.s. proton kinetic energy reconstructed from measured track length using tracking algorithm (left figure). Both quantities are from MC simulation. From this relationship, unfolding matrix was made (right figure). 73

76 Detector unfolding matrix U det is defined as U det = Number of events in the ith bin of true proton kinetic energy is obtained using the unfolding matrix as follows: n true i = j u ji n rec j (4.1) where j is the bin of reconstructed proton kinetic energy, n rec j is the number of events in jth bin, and u ji is the (j, i) element of the unfolding matrix. Figure 4.17 shows proton kinetic energy distribution after detector unfolding. 4.7 Efficiency Correction After detector unfolding, proton kinetic energy distribution is still biased by detector acceptance and tracking efficiency. To obtain unbiased distribution, correction for detection efficiency is needed. Detection efficiency is defined using MC for each true proton kinetic energy bin as number of events in NC elastic sample divided by total number of events generated in SciBar. Figure 4.18 shows the detection efficiencies for each proton kinetic energy bin. Yield of ith bin Y i is written as Y i = n true i ǫ i (4.2) where ǫ i is the detection efficiency for ith true proton kinetic energy bin. Proton kinetic energy distribution after efficiency correction is shown in figure

77 Figure 4.17: Proton kinetic energy distribution after detector unfolding. Data are compared to MC simulation. 75

78 Figure 4.18: MC true proton kinetic energy v.s. detection efficiency. 76

79 Figure 4.19: Proton kinetic energy distribution after efficiency correction. Data are compared to MC simulation. 77

80 4.8 Cross Section Cross section is evaluated from the yield, the number of neutrinos that entered the detector, and the number of targets: dσ i de p = Y i /( E p ) N POT N target φν (E ν )de ν (4.3) where dσ i /de p is ith E p bin of the differential cross section. Y i is a yield of ith bin after several corrections. E p is the bin width of the E p distribution. N POT and N target are the number of proton on beryllium target and the number of nucleons in SciBar FV. φ ν (E ν ) is the neutrino flux. For this analysis, data sample of POT are used. Neutrino flux cannot be measured directly. It is evaluated using neutrino beam MC simulation. Integrated flux φ ν (E ν )de ν is determined to be (1/(cm 2 POT)) using beam MC simulation. N target is calculated by N target = N A M SB where N A is Avogadro s number, M SB = g is the mass of SciBar FV. Differential cross section as a function of proton kinetic energy is shown in figure Error bars include statistical error only. Systematic error is estimated in Section

81 Figure 4.20: Differential cross section as a function of proton kinetic energy. 79

82 4.9 Systematic Uncertainty Neutrino Beam The sources of the systematic uncertainty in the neutrino beam mainly come from secondary particle production, hadronic interactions, horn magnetic field and number of protons on target. Details are written in reference [20]. Secondary Particle Production Secondary particles such as pions and kaons from the proton-beryllium target interaction is simulated by the beam simulation using two models. The Sanford-Wang model is used for π +,π and K 0 productions. The Feynman scaling model is used for K + production. The Sanford-Wang parameterization and the Feynman scaling parameterization are discussed in Section 3.6. The uncertainty in the fitted parameters is evaluated by varying the parameters of the models. Hadronic Interactions: The uncertainty in the hadronic interactions such as proton-beryllium interactions and the probability for mesons to survive possible hadronic interactions in the target or horn and decay to produce neutrinos. Uncertainties in the cross sections of these interactions affect both the rate and shape of the flux. Horn Magnetic Field: The uncertainty in the modeling of the current within the inner cylinder in the magnetic focusing horn is taken into account. The horn current has uncertainty of ±1 ka. Protons On Target: The uncertainty in the measurement of the number of protons delivered to the target is evaluated to be 2%. This uncertainty is considered as a normalization error of the cross section. The change in the neutrino flux due to these uncertainties is calculated by using random parameter vectors and weighting each event by a factor corresponding to the variation of the yield of the parent meson with the given 80

83 momentum and angle. 1,000 parameter vector sets are used for the calculation. Systematic uncertainty in the NC elastic cross section due to the uncertainty in the neutrino beam is shown in Figure Figure 4.21: Systematic uncertainty in the NC elastic cross section due to the uncertainty in the neutrino beam Neutrino Interaction Axial vector mass M A is set to 1.1 GeV/c 2. The uncertainty in M A is estimated to be approximately ±0.1 GeV/c 2 based on results from K2K [68] and MiniBooNE [69] experiments. M A is varied to 1.2 GeV/c 2 considering 81

84 the results from these experiments and the other results from the past experiments [70] to estimate the systematic uncertainty due to M A value. Figure 4.22 shows the systematic uncertainty due to M A uncertainty. Figure 4.22: The systematic uncertainty due to M A uncertainty Detector Response The uncertainty in detector response comes from the measurement error of the cross talk, the scintillator quenching, single photoelectron resolution and the hit threshold for track reconstruction. Since only tracks with large energy deposit are used for NC elastic analysis, the uncertainty in the crosstalk and 82

85 the scintillator quenching are dominant. Crosstalk Crosstalk of the MAPMT is measured for an adjacent channel, a diagonal channel and a next-to-next channel. The amount of the crosstalk for an adjacent channel and an diagonal channel are 3.15% and 0.7%, respectively. The amount of the crosstalk for a next-to-next channel is from 0.1% to 0.3%. The absolute error of the crosstalk to the adjacent channel is estimated to be 0.4%. Scintillator Quenching Scintillator quenching is simulated using Birk s law [64]. Birk s constant of the SciBar scintillator is measured to be ± (stat.) ± (sys.) cm/mev [65]. The uncertainty in NC elastic cross section due to the crosstalk is estimated by varying crosstalk value to 2.75% and 3.55% in the simulation. The uncertainty in NC elastic cross section due to the scintillator quenching is estimated by varying the Birk s constant by ± cm/mev in the simulation. Figure 4.23 shows the systematic uncertainty due to the uncertainty in the detector response Total Systematic Uncertainty All the systematic uncertainties discussed in this section are combined. Each systematic uncertainty is independent to others. Figure 5.1 shows total systematic uncertainty in NC elastic scattering cross section as a function of proton kinetic energy. 83

86 Figure 4.23: The systematic uncertainty due to the uncertainty in the detector response. 84

87 Figure 4.24: Total systematic uncertainty in NC elastic scattering cross section. 85

88 4.10 Nuclear Effect Correction Q 2 is calculated correctly from recoiled proton kinetic energy if the target proton is at rest in the initial state and does not interact with any nucleons. However, in the real case, the target protons have Fermi momentum. They are also bound with nuclear binding energy. Furthermore, they are occasionally rescattered by other nucleons in the same nucleus. Recoiled proton kinetic energy or scattering angle change due to these effects. These effects are called nuclear effects. It mainly comes from Fermi momentum in the nuclei. Since Q 2 is reconstructed from recoiled proton kinetic energy, Q 2 is smeared. Figure 4.25 shows a distribution of Q 2 reconstructed from proton kinetic energy v.s. Q 2 calculated from the neutrino momentum using Monte-Carlo simulation. Figure 4.25: Q 2 reconstructed from proton kinetic energy (left) v.s. Q 2 calculated from the neutrino momentum using Monte-Carlo simulation (right). Unsmearing matrix for nuclear correction U nucl is made from this distribution. U nucl is defined as 86

SciBar and future K2K physics. F.Sánchez Universitat Aútonoma de Barcelona Institut de Física d'altes Energies

SciBar and future K2K physics. F.Sánchez Universitat Aútonoma de Barcelona Institut de Física d'altes Energies SciBar and future K2K physics F.Sánchez Universitat Aútonoma de Barcelona Institut de Física d'altes Energies ICRR, 29 th October 2003 Outline Introduction: K2K SciBar detector: Physics goals Design Electron

More information

Single π production in MiniBooNE and K2K

Single π production in MiniBooNE and K2K Single π production in MiniBooNE and K2K Yoshinari Hayato (ICRR, Univ. Tokyo) 1. Introduction 2. Study of neutral current single π 0 production at K2K 3. Study of charged current coherent π + production

More information

Neutrino-Nucleus Scattering at MINERvA

Neutrino-Nucleus Scattering at MINERvA 1 Neutrino-Nucleus Scattering at MINERvA Elba XIII Workshop: Neutrino Physics IV Tammy Walton Fermilab June 26, 2014 2 MINERvA Motivation Big Picture Enter an era of precision neutrino oscillation measurements.

More information

Neutrino interaction at K2K

Neutrino interaction at K2K Neutrino interaction at K2K KEK Yoshinari Hayato for the K2K collaboration Contents Introduction Neutrino flux measurements and neutrino interaction studies in K2K Charged current quasi-elastic scattering

More information

Long Baseline Neutrinos

Long Baseline Neutrinos Long Baseline Neutrinos GINA RAMEIKA FERMILAB SLAC SUMMER INSTITUTE AUGUST 5-6, 2010 Lecture 1 Outline Defining Long Baseline Experiment Ingredients Neutrino Beams Neutrino Interactions Neutrino Cross

More information

Charged current single pion to quasi-elastic cross section ratio in MiniBooNE. Steven Linden PAVI09 25 June 2009

Charged current single pion to quasi-elastic cross section ratio in MiniBooNE. Steven Linden PAVI09 25 June 2009 Charged current single pion to quasi-elastic cross section ratio in MiniBooNE Steven Linden PAVI09 25 June 2009 Motivation Oscillation searches needed for leptonic CP violation. One approach: search for

More information

HARP (Hadron Production) Experiment at CERN

HARP (Hadron Production) Experiment at CERN HARP (Hadron Production) Experiment at CERN 2nd Summer School On Particle Accelerators And Detectors 18-24 Sep 2006, Bodrum, Turkey Aysel Kayιş Topaksu Çukurova Üniversitesi, ADANA Outline The Physics

More information

MiniBooNE Progress and Little Muon Counter Overview

MiniBooNE Progress and Little Muon Counter Overview MiniBooNE Progress and Little Muon Counter Overview Neutrino Introduction and MiniBooNE Motivation MiniBooNE Detector and LMC Calibration and Performance Progress Toward Physics Results Terry Hart, University

More information

Analysis of muon and electron neutrino charged current interactions in the T2K near detectors

Analysis of muon and electron neutrino charged current interactions in the T2K near detectors Nuclear Physics B Proceedings Supplement (14) 1 6 Nuclear Physics B Proceedings Supplement Analysis of muon and electron neutrino charged current interactions in the TK near detectors A. Hillairet (On

More information

The Development of Particle Physics. Dr. Vitaly Kudryavtsev E45, Tel.:

The Development of Particle Physics. Dr. Vitaly Kudryavtsev E45, Tel.: The Development of Particle Physics Dr. Vitaly Kudryavtsev E45, Tel.: 0114 4531 v.kudryavtsev@sheffield.ac.uk The structure of the nucleon Electron - nucleon elastic scattering Rutherford, Mott cross-sections

More information

Neutrino Cross Sections and Scattering Physics

Neutrino Cross Sections and Scattering Physics Neutrino Cross Sections and Scattering Physics Bonnie Fleming Yale University, New Haven, CT. Abstract. Large flux uncertainties and small cross sections have made neutrino scattering physics a challenge.

More information

The achievements of the CERN proton antiproton collider

The achievements of the CERN proton antiproton collider The achievements of the CERN proton antiproton collider Luigi DiLella Scuola Normale Superiore, Pisa, Italy Motivation of the project The proton antiproton collider UA1 and UA2 detectors Discovery of the

More information

Introduction to Neutrino Interaction Physics NUFACT08 Summer School June 2008 Benasque, Spain Paul Soler

Introduction to Neutrino Interaction Physics NUFACT08 Summer School June 2008 Benasque, Spain Paul Soler Introduction to NUCT08 Summer School -3 June 008 Benasque, Spain Paul Soler 4. uasi-elastic, resonant, coherent and diffractive scattering 4. Motivation 4. Charged current quasi-elastic scattering 4.3

More information

Design, construction, and initial performance of SciBar detector in K2K experiment

Design, construction, and initial performance of SciBar detector in K2K experiment Design, construction, and initial performance of SciBar detector in K2K experiment Shimpei Yamamoto for the K2K-SciBar group Abstract-- A new near detector in K2K long baseline neutrino experiment, the

More information

Neutrino Cross Section Measurements for Long-Baseline Acceleratorbased Neutrino Oscillation Experiments

Neutrino Cross Section Measurements for Long-Baseline Acceleratorbased Neutrino Oscillation Experiments Neutrino Cross Section Measurements for Long-Baseline Acceleratorbased Neutrino Oscillation Experiments Katori, T arxiv record: http://arxiv.org/abs/0805.476 For additional information about this publication

More information

Charged Current Quasielastic Analysis from MINERνA

Charged Current Quasielastic Analysis from MINERνA Charged Current Quasielastic Analysis from MINERνA Anushree Ghosh CBPF - Centro Brasileiro de Pesquisas F sicas, Rua Dr. Xavier Sigaud 150, Urca, Rio de Janeiro, Rio de Janeiro, 22290-180, Brazil on behalf

More information

arxiv: v1 [nucl-ex] 7 Sep 2009

arxiv: v1 [nucl-ex] 7 Sep 2009 CCπ Event Reconstruction at MiniBooNE R.H. Nelson University of Colorado, Dept. of Physics, 39 UCB, Boulder, CO 839, USA arxiv:99.238v [nucl-ex] 7 Sep 29 Abstract. We describe the development of a fitter

More information

MINOS. Luke A. Corwin, for MINOS Collaboration Indiana University XIV International Workshop On Neutrino Telescopes 2011 March 15

MINOS. Luke A. Corwin, for MINOS Collaboration Indiana University XIV International Workshop On Neutrino Telescopes 2011 March 15 MINOS Luke A. Corwin, for MINOS Collaboration Indiana University XIV International Workshop On Neutrino Telescopes 2011 March 15 2 Overview and Current Status Beam Detectors Analyses Neutrino Charged Current

More information

Fall Quarter 2010 UCSB Physics 225A & UCSD Physics 214 Homework 1

Fall Quarter 2010 UCSB Physics 225A & UCSD Physics 214 Homework 1 Fall Quarter 2010 UCSB Physics 225A & UCSD Physics 214 Homework 1 Problem 2 has nothing to do with what we have done in class. It introduces somewhat strange coordinates called rapidity and pseudorapidity

More information

PMT Signal Attenuation and Baryon Number Violation Background Studies. By: Nadine Ayoub Nevis Laboratories, Columbia University August 5, 2011

PMT Signal Attenuation and Baryon Number Violation Background Studies. By: Nadine Ayoub Nevis Laboratories, Columbia University August 5, 2011 PMT Signal Attenuation and Baryon Number Violation Background Studies By: Nadine Ayoub Nevis Laboratories, Columbia University August 5, 2011 1 The Standard Model The Standard Model is comprised of Fermions

More information

Neutrino Physics at Short Baseline

Neutrino Physics at Short Baseline Neutrino Physics at Short Baseline E. D. Zimmerman University of Colorado Lepton-Photon 2005 Uppsala, Sweden 2 July 2005 Neutrino Physics at Short Baseline Parameter regions accessible with short baselines

More information

K2K and T2K experiments

K2K and T2K experiments K2K and T2K experiments Issei Kato TRIUMF, Canada for the K2K and T2K collaborations Neutrino Oscillation Workshop 2006 at Conca Specchiulla,, Otranto, Italy Outline K2K experiment, shortly What s updated

More information

Neutrino Scattering in Liquid Argon TPC Detectors

Neutrino Scattering in Liquid Argon TPC Detectors Bonnie T. Fleming NuFACT 05 INFN Neutrino Scattering in Liquid Argon TPC Detectors What can we learn? Scattering at 1 GeV Neutrino Scattering using Liquid Argon TPCS with conventional neutrino beams in

More information

Outline. Charged Leptonic Weak Interaction. Charged Weak Interactions of Quarks. Neutral Weak Interaction. Electroweak Unification

Outline. Charged Leptonic Weak Interaction. Charged Weak Interactions of Quarks. Neutral Weak Interaction. Electroweak Unification Weak Interactions Outline Charged Leptonic Weak Interaction Decay of the Muon Decay of the Neutron Decay of the Pion Charged Weak Interactions of Quarks Cabibbo-GIM Mechanism Cabibbo-Kobayashi-Maskawa

More information

The First Results of K2K long-baseline Neutrino Oscillation Experiment

The First Results of K2K long-baseline Neutrino Oscillation Experiment The First Results of K2K long-baseline Neutrino Oscillation Experiment Taku Ishida, representing K2K collaboration arxiv:hep-ex/0008047v1 23 Aug 2000 Institute for Particle and Nuclear Studies(IPNS) High

More information

NA62: Ultra-Rare Kaon Decays

NA62: Ultra-Rare Kaon Decays NA62: Ultra-Rare Kaon Decays Phil Rubin George Mason University For the NA62 Collaboration November 10, 2011 The primary goal of the experiment is to reconstruct more than 100 K + π + ν ν events, over

More information

Model independent extraction of the axial mass parameter in CCQE anti neutrino-nucleon scattering

Model independent extraction of the axial mass parameter in CCQE anti neutrino-nucleon scattering Model independent extraction of the axial mass parameter in CCQE anti neutrino-nucleon scattering Jerold Eugene Young III Department of Physics, Western Illinois University dvisor Dr. Gil Paz Wayne State

More information

Measuring Form Factors and Structure Functions With CLAS

Measuring Form Factors and Structure Functions With CLAS Measuring Form Factors and Structure Functions With CLAS Jerry Gilfoyle for the CLAS Collaboration Physics Department, University of Richmond, Virginia Outline: 1. Jefferson Lab and the CLAS Detector..

More information

The reaction p(e,e'p)π 0 to calibrate the Forward and the Large Angle Electromagnetic Shower Calorimeters

The reaction p(e,e'p)π 0 to calibrate the Forward and the Large Angle Electromagnetic Shower Calorimeters The reaction p(e,e'p)π 0 to calibrate the Forward and the Large Angle Electromagnetic Shower Calorimeters M.Battaglieri, M.Anghinolfi, P.Corvisiero, A.Longhi, M.Ripani, M.Taiuti Istituto Nazionale di Fisica

More information

Lecture 8. CPT theorem and CP violation

Lecture 8. CPT theorem and CP violation Lecture 8 CPT theorem and CP violation We have seen that although both charge conjugation and parity are violated in weak interactions, the combination of the two CP turns left-handed antimuon onto right-handed

More information

Recent Results from Alysia Marino, University of Colorado at Boulder Neutrino Flux Workshop, University of Pittsburgh, Dec 6 8,2012

Recent Results from Alysia Marino, University of Colorado at Boulder Neutrino Flux Workshop, University of Pittsburgh, Dec 6 8,2012 Recent Results from Alysia Marino, University of Colorado at Boulder Neutrino Flux Workshop, University of Pittsburgh, Dec 6 8,2012 Outline T2K Beamline and Detectors Uncertainties from prior measurements

More information

CCQE and SciBooNE, T2K. 12/3/2013 K Mahn, INT

CCQE and SciBooNE, T2K. 12/3/2013 K Mahn, INT CCQE and SciBooNE, T2K 12/3/2013 K Mahn, INT Neutrino flux Neutrino flux neutrino energy (E ) in GeV Predominantly ~6%, ~0.5% ε + ε In MiniBooNE: ~1 per 1e15 POT In SciBooNE: ~0.5 per 1e15 POT ~5x closer,

More information

Results from HARP. Malcolm Ellis On behalf of the HARP collaboration DPF Meeting Riverside, August 2004

Results from HARP. Malcolm Ellis On behalf of the HARP collaboration DPF Meeting Riverside, August 2004 Results from HARP Malcolm Ellis On behalf of the HARP collaboration DPF Meeting Riverside, August 2004 The HAdRon Production Experiment 124 physicists 24 institutes2 Physics Goals Input for precise calculation

More information

An Introduction to Modern Particle Physics. Mark Thomson University of Cambridge

An Introduction to Modern Particle Physics. Mark Thomson University of Cambridge An Introduction to Modern Particle Physics Mark Thomson University of Cambridge Science Summer School: 30 th July - 1 st August 2007 1 Course Synopsis Introduction : Particles and Forces - what are the

More information

Charged Current Inclusive Scattering in MINERnA

Charged Current Inclusive Scattering in MINERnA Charged Current Inclusive Scattering in MINERnA What is Minerna? Why Minerna? n beam and n flux n / n inclusive x-sections x-section ratios (A-depndence) Outlook NUFACT 2013 August 21 st 13 Alessandro

More information

ω γ Neutral Current Single Photon Production (NCγ) Outline 1. Oscillation physics 2. NOMAD 3. T2K/MINERvA 4. MicroBooNE 5. MiniBooNE+ 6.

ω γ Neutral Current Single Photon Production (NCγ) Outline 1. Oscillation physics 2. NOMAD 3. T2K/MINERvA 4. MicroBooNE 5. MiniBooNE+ 6. Neutral Current Single Photon Production (NCγ) Outline physics ν Z ν N ω γ Teppei Katori Queen Mary University of London INT workshop, Seattle, USA, Dec. 12, 2013 N 1 2 NuSTEC protocol - way to avoid Donkey

More information

Neutrino Physics at Short Baseline

Neutrino Physics at Short Baseline Neutrino Physics at Short Baseline E. D. Zimmerman University of Colorado Boulder, Colorado XXVI Physics in Collision Búzios, 6 July 2006 The Amazon from 33,000 feet Neutrino Physics at Short Baseline

More information

SHiP: a new facility with a dedicated detector for neutrino physics

SHiP: a new facility with a dedicated detector for neutrino physics SHiP: a new facility with a dedicated detector for neutrino physics Università "Federico II" and INFN, Naples, Italy E-mail: giovanni.de.lellis@cern.ch The SHiP facility recently proposed at CERN copiously

More information

First neutrino beam and cosmic tracks. GDR neutrino

First neutrino beam and cosmic tracks. GDR neutrino Magali Besnier T2K status tt of the experiment First neutrino beam and cosmic tracks GDR neutrino 28 04 09 T2K Tokai 2 Kamiokande Main goals : 1) Measurement of 13 mixing angle studying the e oscillation

More information

Weak interactions, parity, helicity

Weak interactions, parity, helicity Lecture 10 Weak interactions, parity, helicity SS2011: Introduction to Nuclear and Particle Physics, Part 2 2 1 Weak decay of particles The weak interaction is also responsible for the β + -decay of atomic

More information

Particle Production Measurements at Fermilab

Particle Production Measurements at Fermilab Particle Production Measurements at Fermilab Dr. Nickolas Solomey, IIT and Fermilab co Spokesman of E907 TEV II Astroparticle Physics Conference Univ. of Wisconsin, Madison 28 31 Aug., 2006 Particle Production

More information

Near detector tracker concepts. D. Karlen / U. Vic. & TRIUMF T2K ND280m meeting August 22, 2004

Near detector tracker concepts. D. Karlen / U. Vic. & TRIUMF T2K ND280m meeting August 22, 2004 Near detector tracker concepts D. Karlen / U. Vic. & TRIUMF T2K ND280m meeting August 22, 2004 Longitudinal extent of tracker modules Consensus has developed that the near detector should consist of a

More information

PHYS 5326 Lecture #2. Wednesday, Jan. 24, 2007 Dr. Jae Yu. Wednesday, Jan. 24, 2007 PHYS 5326, Spring 2007 Jae Yu

PHYS 5326 Lecture #2. Wednesday, Jan. 24, 2007 Dr. Jae Yu. Wednesday, Jan. 24, 2007 PHYS 5326, Spring 2007 Jae Yu PHYS 5326 Lecture #2 Wednesday, Jan. 24, 2007 Dr. 1. Sources of Neutrinos 2. How is neutrino beam produced? 3. Physics with neutrino experiments 4. Characteristics of accelerator based neutrino experiments

More information

Experimental Aspects of Deep-Inelastic Scattering. Kinematics, Techniques and Detectors

Experimental Aspects of Deep-Inelastic Scattering. Kinematics, Techniques and Detectors 1 Experimental Aspects of Deep-Inelastic Scattering Kinematics, Techniques and Detectors 2 Outline DIS Structure Function Measurements DIS Kinematics DIS Collider Detectors DIS process description Dirac

More information

Results From The HARP Experiment

Results From The HARP Experiment Results From The HARP Experiment Lara Howlett University of Sheffield Overview of HARP and forward analysis Al analysis for K2K experiment Be analysis for MiniBooNE experiment Motivations Systematic study

More information

Physics at Hadron Colliders

Physics at Hadron Colliders Physics at Hadron Colliders Part 2 Standard Model Physics Test of Quantum Chromodynamics - Jet production - W/Z production - Production of Top quarks Precision measurements -W mass - Top-quark mass QCD

More information

Neutrinos Induced Pion Production in MINERvA

Neutrinos Induced Pion Production in MINERvA Neutrinos Induced Pion Production in MINERvA Aaron Higuera Universidad de Guanajuato, Mexico On behalf of the MINERvA collaboration Outline Neutrinos Induced Pion Production MINERvA Experiment E-938 CC

More information

Kaon Identification at NA62. Institute of Physics Particle, Astroparticle, and Nuclear Physics groups Conference 2015

Kaon Identification at NA62. Institute of Physics Particle, Astroparticle, and Nuclear Physics groups Conference 2015 Kaon Identification at NA62 Institute of Physics Particle, Astroparticle, and Nuclear Physics groups Conference 2015 Francis Newson April 2015 Kaon Identification at NA62 K πνν NA48 and NA62 K + π + νν

More information

The US Liquid Argon Time Projection Chamber (LArTPC) experiments

The US Liquid Argon Time Projection Chamber (LArTPC) experiments The US Liquid Argon Time Projection Chamber (LArTPC) experiments outline 1. Introduction 2. ArgoNeuT experiment 3. MicroBooNE experiment 4. Future large LArTPC experiments 5. Conclusion Teppei Katori for

More information

Study of Strange Quark in the Nucleon with Neutrino Scattering

Study of Strange Quark in the Nucleon with Neutrino Scattering July 28, 2004 NuFact 04, Osaka Study of Strange Quark in the Nucleon with Neutrino Scattering T.-A. Shibata Tokyo Institute of Technology Contents: 3. Physics Motivation --- Quark Structure of the Nucleon

More information

Neutrino Scattering Results from MiniBooNE and SciBooNE

Neutrino Scattering Results from MiniBooNE and SciBooNE Neutrino Scattering Results from MiniBooNE and SciBooNE Outline: introduction, motivation MiniBooNE/SciBooNE experiments measurements, results interpretations future conclusions R. Tayloe, Indiana U. e-n

More information

Neutrino Helicity Measurement

Neutrino Helicity Measurement PHYS 851 Introductory Nuclear Physics Instructor: Chary Rangacharyulu University of Saskatchewan Neutrino Helicity Measurement Stefan A. Gärtner stefan.gaertner@gmx.de December 9 th, 2005 2 1 Introduction

More information

Neutrino Physics at MiniBooNE

Neutrino Physics at MiniBooNE Neutrino Physics at MiniBooNE G. P. Zeller Columbia University Department of Physics 538 W. 120th St., New York, NY 10027, USA for the MiniBooNE Collaboration As its primary design goal, the MiniBooNE

More information

CLAS12 at Jefferson Lab

CLAS12 at Jefferson Lab CLAS12 at Jefferson Lab Daria Sokhan University of Glasgow, UK IPPP/NuSTEC Topical Meeting on Neutrino-Nucleus Scattering IPPP, Durham, UK 19 April 2017 Jefferson Lab 6 GeV era Jefferson Lab CEBAF: Continuous

More information

Standard Model of Particle Physics SS 2013

Standard Model of Particle Physics SS 2013 Lecture: Standard Model of Particle Physics Heidelberg SS 013 Weak Interactions II 1 Important Experiments Wu-Experiment (1957): radioactive decay of Co60 Goldhaber-Experiment (1958): radioactive decay

More information

Cross Section of Exclusive π Electro-production from Neutron. Jixie Zhang (CLAS Collaboration) Old Dominion University Sep. 2009

Cross Section of Exclusive π Electro-production from Neutron. Jixie Zhang (CLAS Collaboration) Old Dominion University Sep. 2009 Cross Section of Exclusive π Electro-production from Neutron Jixie Zhang (CLAS Collaboration) Old Dominion University Sep. 2009 Exclusive π electro-production Detect e`, π and at least ONE of the two final

More information

RECENT RESULTS FROM THE OPERA EXPERIMENT

RECENT RESULTS FROM THE OPERA EXPERIMENT RECENT RESULTS FROM THE OPERA EXPERIMENT P. DEL AMO SANCHEZ on behalf of the OPERA Collaboration LAPP, 9, Chemin du Bellevue, BP 110 74941 Annecy-le-Vieux, FRANCE The OPERA neutrino experiment recently

More information

The CNGS neutrino beam

The CNGS neutrino beam 10th Topical Seminar on Innovative Particle and Radiation Detectors (IPRD06) 1-5 October 2006 Siena, Italy ν The CNGS neutrino beam G. Sirri INFN Bologna CNGS (CERN Neutrinos to Gran Sasso) The project

More information

Latest Results from the OPERA Experiment (and new Charge Reconstruction)

Latest Results from the OPERA Experiment (and new Charge Reconstruction) Latest Results from the OPERA Experiment (and new Charge Reconstruction) on behalf of the OPERA Collaboration University of Hamburg Institute for Experimental Physics Overview The OPERA Experiment Oscillation

More information

COMET muon conversion experiment in J-PARC

COMET muon conversion experiment in J-PARC Institute for Basic Science, Daejeon, Korea E-mail: myeongjaelee@ibs.re.kr COMET is an experiment at J-PARC, Japan, which will search for neutrinoless conversion of muons into electrons in the field of

More information

energy loss Ionization + excitation of atomic energy levels Mean energy loss rate de /dx proportional to (electric charge) 2 of incident particle

energy loss Ionization + excitation of atomic energy levels Mean energy loss rate de /dx proportional to (electric charge) 2 of incident particle Lecture 4 Particle physics processes - particles are small, light, energetic à processes described by quantum mechanics and relativity à processes are probabilistic, i.e., we cannot know the outcome of

More information

A Measurement of the Induced polarization of electro-produced Λ(1116) with CLAS

A Measurement of the Induced polarization of electro-produced Λ(1116) with CLAS A Measurement of the Induced polarization of electro-produced Λ(1116) with CLAS Marianna Gabrielyan Florida International University HUGS 2008 Why study electromagnetic production of kaons? Formalism.

More information

MINOS. Physics Program and Construction Status. Karol Lang The University of Texas at Austin. YITP: Neutrinos and Implications for Physics Beyond

MINOS. Physics Program and Construction Status. Karol Lang The University of Texas at Austin. YITP: Neutrinos and Implications for Physics Beyond MINOS Physics Program and Construction Status Karol Lang The University of Texas at Austin YITP: Neutrinos and Implications for Physics Beyond YITP Conference: Neutrinos and Implications The Standard for

More information

Recent results from MiniBooNE on neutrino oscillations

Recent results from MiniBooNE on neutrino oscillations Recent results from MiniBooNE on neutrino oscillations Alexis [for the MiniBooNE collaboration] IX International Symposium on High Energy Physics 1 SILAFAE 2012,, December 10 14, 2013 Outline LSND and

More information

Physics Quarknet/Service Learning

Physics Quarknet/Service Learning Physics 29000 Quarknet/Service Learning Lecture 3: Ionizing Radiation Purdue University Department of Physics February 1, 2013 1 Resources Particle Data Group: http://pdg.lbl.gov Summary tables of particle

More information

STATUS OF ATLAS TILE CALORIMETER AND STUDY OF MUON INTERACTIONS. 1 Brief Description of the ATLAS Tile Calorimeter

STATUS OF ATLAS TILE CALORIMETER AND STUDY OF MUON INTERACTIONS. 1 Brief Description of the ATLAS Tile Calorimeter STATUS OF ATLAS TILE CALORIMETER AND STUDY OF MUON INTERACTIONS L. E. PRICE Bldg 362, Argonne National Laboratory, Argonne, IL 60439, USA E-mail: lprice@anl.gov (For the ATLAS Tile Calorimeter Collaboration)

More information

PoS(KAON13)012. R K Measurement with NA62 at CERN SPS. Giuseppe Ruggiero. CERN

PoS(KAON13)012. R K Measurement with NA62 at CERN SPS. Giuseppe Ruggiero. CERN CERN E-mail: giuseppe.ruggiero@cern.ch A precision measurement of the ratio R K of the rates of kaon leptonic decays K e ν and K µ ν with the full minimum bias data sample collected with low intensity

More information

Thin Calorimetry for Cosmic-Ray Studies Outside the Earth s Atmosphere. 1 Introduction

Thin Calorimetry for Cosmic-Ray Studies Outside the Earth s Atmosphere. 1 Introduction Thin Calorimetry for Cosmic-Ray Studies Outside the Earth s Atmosphere Richard WIGMANS Department of Physics, Texas Tech University, Lubbock TX 79409-1051, USA (wigmans@ttu.edu) Abstract Cosmic ray experiments

More information

PoS(KAON09)023. Beam Hole Photon Veto For J-PARC K O TO experiment. Yosuke Maeda Kyoto University

PoS(KAON09)023. Beam Hole Photon Veto For J-PARC K O TO experiment. Yosuke Maeda Kyoto University Beam Hole Photon Veto For J-PARC K O TO experiment Kyoto University E-mail: maeda_y@scphys.kyoto-u.ac.jp The Beam Hole Photon Veto counter (BHPV) for the J-PARC K O TO experiment was designed by MC simulation.

More information

1 Introduction. KOPIO charged-particle vetos. K - RARE Meeting (Frascati) May Purpose of CPV: veto Kl

1 Introduction. KOPIO charged-particle vetos. K - RARE Meeting (Frascati) May Purpose of CPV: veto Kl Introduction - Purpose of CPV: veto Kl decay modes with a real or apparent π and a pair of charged particles - Examples of background modes: (i) K l π π + π (ii) K l π π ± eν there are always (iii) K l

More information

INTRODUCTION TO THE STANDARD MODEL OF PARTICLE PHYSICS

INTRODUCTION TO THE STANDARD MODEL OF PARTICLE PHYSICS INTRODUCTION TO THE STANDARD MODEL OF PARTICLE PHYSICS Class Mechanics My office (for now): Dantziger B Room 121 My Phone: x85200 Office hours: Call ahead, or better yet, email... Even better than office

More information

PoS(NOW2016)003. T2K oscillation results. Lorenzo Magaletti. INFN Sezione di Bari

PoS(NOW2016)003. T2K oscillation results. Lorenzo Magaletti. INFN Sezione di Bari INFN Sezione di Bari E-mail: lorenzo.magaletti@ba.infn.it The TK (Tokai-to-Kamioka) experiment is a second generation long baseline neutrino oscillation experiment that probes physics beyond the Standard

More information

The HARP Experiment. G. Vidal-Sitjes (INFN-Ferrara) on behalf of the HARP Collaboration

The HARP Experiment. G. Vidal-Sitjes (INFN-Ferrara) on behalf of the HARP Collaboration The HARP Experiment (INFN-Ferrara) on behalf of the HARP Collaboration New Views in Particle Physics Outline Goals for a HAdRon Production experiment Example: KEK PS Neutrino beam-line Detector layout

More information

Quasi-Elastic Scattering in MINERvA

Quasi-Elastic Scattering in MINERvA Quasi-Elastic Scattering in MINERvA Kevin S. McFarland Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 USA and Fermi National Accelerator Laboratory, Batavia, IL 60510

More information

arxiv: v1 [hep-ex] 1 Oct 2015

arxiv: v1 [hep-ex] 1 Oct 2015 CIPANP2015-Galati October 2, 2015 OPERA neutrino oscillation search: status and perspectives arxiv:1510.00343v1 [hep-ex] 1 Oct 2015 Giuliana Galati 1 Università degli Studi di Napoli Federico II and INFN

More information

Outline. Charged Leptonic Weak Interaction. Charged Weak Interactions of Quarks. Neutral Weak Interaction. Electroweak Unification

Outline. Charged Leptonic Weak Interaction. Charged Weak Interactions of Quarks. Neutral Weak Interaction. Electroweak Unification Weak Interactions Outline Charged Leptonic Weak Interaction Decay of the Muon Decay of the Neutron Decay of the Pion Charged Weak Interactions of Quarks Cabibbo-GIM Mechanism Cabibbo-Kobayashi-Maskawa

More information

Neutrino Interaction Physics for Oscillation Analyses

Neutrino Interaction Physics for Oscillation Analyses Neutrino Interaction Physics for Oscillation Analyses Patrick Stowell on behalf of the T2K Neutrino Interactions Working Group NuIntUK 20th July 2015 20th July 2015 Neutrino Interaction Physics for Oscillation

More information

DESY Summer Students Program 2008: Exclusive π + Production in Deep Inelastic Scattering

DESY Summer Students Program 2008: Exclusive π + Production in Deep Inelastic Scattering DESY Summer Students Program 8: Exclusive π + Production in Deep Inelastic Scattering Falk Töppel date: September 6, 8 Supervisors: Rebecca Lamb, Andreas Mussgiller II CONTENTS Contents Abstract Introduction.

More information

Neutrino Nucleus Deep Inelastic Scattering at MINERvA

Neutrino Nucleus Deep Inelastic Scattering at MINERvA Neutrino Nucleus Deep Inelastic Scattering at MINERvA University of Florida E-mail: wospakrk@fnal.gov Neutrino-nucleus charged-current deep inelastic scattering (DIS) provides a complementary probe to

More information

Electroweak Physics. Krishna S. Kumar. University of Massachusetts, Amherst

Electroweak Physics. Krishna S. Kumar. University of Massachusetts, Amherst Electroweak Physics Krishna S. Kumar University of Massachusetts, Amherst Acknowledgements: M. Grunewald, C. Horowitz, W. Marciano, C. Quigg, M. Ramsey-Musolf, www.particleadventure.org Electroweak Physics

More information

HERA e-p scattering events observed in the H1Detector. H1 Events Joachim Meyer DESY

HERA e-p scattering events observed in the H1Detector. H1 Events Joachim Meyer DESY HERA e-p scattering events observed in the H1Detector H1 Events Joachim Meyer DESY 2005 1 The idea The realisation The Physics The events H1 Events Joachim Meyer DESY 2005 2 What we think what happens,

More information

Study of ν τ production by measuring D s τ events in 400 GeV proton interactions: Test of lepton universality in neutrino charged-current interactions

Study of ν τ production by measuring D s τ events in 400 GeV proton interactions: Test of lepton universality in neutrino charged-current interactions Letter of intent Study of ν τ production by measuring D s τ events in 400 GeV proton interactions: Test of lepton universality in neutrino charged-current interactions S. Aoki 1, A. Ariga 2, T. Ariga 2,

More information

Results and Status from HARP and MIPP

Results and Status from HARP and MIPP Results and Status from HARP and MIPP M. Sorel (IFIC, CSIC Valencia U.) Neutrino 08, May 25 31, Christchurch (New Zealand) Outline The experiments The data Hadron production for neutrino physics: Results

More information

On the measurements of neutrino energy spectra and nuclear effects in neutrino-nucleus interactions

On the measurements of neutrino energy spectra and nuclear effects in neutrino-nucleus interactions On the measurements of neutrino energy spectra and nuclear effects in neutrino-nucleus interactions Xianguo LU ( 盧 顕国 ) University of Oxford Kyoto HE Seminar 1 Outline 1. Introduction 2. Measuring nuclear

More information

Concepts of Event Reconstruction

Concepts of Event Reconstruction August 3, 2007 Directly Detectable Particles electrons, positrons: e ±, lightest charged lepton photons: γ, gauge boson for electromagnetic force pions: π ±, lightest mesons kaons: K ±, K L, lightest strange

More information

Muon reconstruction performance in ATLAS at Run-2

Muon reconstruction performance in ATLAS at Run-2 2 Muon reconstruction performance in ATLAS at Run-2 Hannah Herde on behalf of the ATLAS Collaboration Brandeis University (US) E-mail: hannah.herde@cern.ch ATL-PHYS-PROC-205-2 5 October 205 The ATLAS muon

More information

PoS(HEP2005)177. Status of the OPERA experiment. Francesco Di Capua. Universita Federico II and INFN.

PoS(HEP2005)177. Status of the OPERA experiment. Francesco Di Capua. Universita Federico II and INFN. Francesco Di Capua Universita Federico II and INFN E-mail: francesco.di.capua@cern.ch The OPERA experiment has been designed to observe an unambiguous signal of ν µ ν τ oscillation in the parameter region

More information

The MINERnA Experiment

The MINERnA Experiment The MINERnA Experiment What is Minerna? Why Mienrna? n / n CCQE sections Inclusive n sections n beam and n flux Outlook INPC 2013 June 4 th 13 Alessandro Bravar for the Mienrna Collaboration The MINERnA

More information

Recent results from MINERvA

Recent results from MINERvA Recent results from MINERvA Jonathan Miller for the MINERvA Collaboration Universidad Tecnica Federico Santa Maria 1 PANIC NOW Outline Introduction to MINERvA Recent important results: charged pion production

More information

New Limits on Heavy Neutrino from NA62

New Limits on Heavy Neutrino from NA62 CERN E-mail: michal.koval@cern.ch The NA6 experiment at CERN collected large samples of charged kaon decays in flight with a minimum bias trigger configuration in 7 and in 15 using a completely new detector

More information

Lecture 11. Weak interactions

Lecture 11. Weak interactions Lecture 11 Weak interactions 1962-66: Formula/on of a Unified Electroweak Theory (Glashow, Salam, Weinberg) 4 intermediate spin 1 interaction carriers ( bosons ): the photon (γ) responsible for all electromagnetic

More information

OUTLINE. CHARGED LEPTONIC WEAK INTERACTION - Decay of the Muon - Decay of the Neutron - Decay of the Pion

OUTLINE. CHARGED LEPTONIC WEAK INTERACTION - Decay of the Muon - Decay of the Neutron - Decay of the Pion Weak Interactions OUTLINE CHARGED LEPTONIC WEAK INTERACTION - Decay of the Muon - Decay of the Neutron - Decay of the Pion CHARGED WEAK INTERACTIONS OF QUARKS - Cabibbo-GIM Mechanism - Cabibbo-Kobayashi-Maskawa

More information

High Energy Physics. Lecture 9. Deep Inelastic Scattering Scaling Violation. HEP Lecture 9 1

High Energy Physics. Lecture 9. Deep Inelastic Scattering Scaling Violation. HEP Lecture 9 1 High Energy Physics Lecture 9 Deep Inelastic Scattering Scaling Violation HEP Lecture 9 1 Deep Inelastic Scattering: The reaction equation of DIS is written e+ p e+ X where X is a system of outgoing hadrons

More information

THE FLORIDA STATE UNIVERSITY COLLEGE OF ARTS & SCIENCES W BOSON PRODUCTION CHARGE ASYMMETRY IN THE ELECTRON CHANNEL ASHLEY S HUFF

THE FLORIDA STATE UNIVERSITY COLLEGE OF ARTS & SCIENCES W BOSON PRODUCTION CHARGE ASYMMETRY IN THE ELECTRON CHANNEL ASHLEY S HUFF THE FLORIDA STATE UNIVERSITY COLLEGE OF ARTS & SCIENCES W BOSON PRODUCTION CHARGE ASYMMETRY IN THE ELECTRON CHANNEL By ASHLEY S HUFF A Thesis submitted to the Department of Physics In partial fulfillment

More information

Time-like Compton Scattering with transversely polarized target

Time-like Compton Scattering with transversely polarized target Time-like Compton Scattering with transversely polarized target Vardan Tadevosyan AANSL (YerPhI) Foundation Arthur Mkrtchyan CUA Outline Physics case and motivation Experimental setup Simulation results

More information

CLEO c. Anders Ryd Cornell University June 7, e e cc D D. D K,D K e

CLEO c. Anders Ryd Cornell University June 7, e e cc D D. D K,D K e CLEOc Cornell University June 7, 25 e e cc D D D K,D K e K K e 1 Outline CLEOc experiment and the physics program Some early results D> Absolute hadronic branching fractions Semileptonic decays 2 Testing

More information

Neutrino Oscillations and the Matter Effect

Neutrino Oscillations and the Matter Effect Master of Science Examination Neutrino Oscillations and the Matter Effect RAJARSHI DAS Committee Walter Toki, Robert Wilson, Carmen Menoni Overview Introduction to Neutrinos Two Generation Mixing and Oscillation

More information

Lecture 3. lecture slides are at:

Lecture 3. lecture slides are at: Lecture 3 lecture slides are at: http://www.physics.smu.edu/ryszard/5380fa16/ Proton mass m p = 938.28 MeV/c 2 Electron mass m e = 0.511 MeV/c 2 Neutron mass m n = 939.56 MeV/c 2 Helium nucleus α: 2 protons+2

More information

Variation in MC prediction of MB nu flux

Variation in MC prediction of MB nu flux Hadron Production Results from the HARP Experiment Linda Coney For the HARP Collaboration Beach 2008 June 27, 2008 Outline Introduction Experiment & Detector Data Hadron Production Results Results for

More information

Discovery of muon neutrino

Discovery of muon neutrino Discovery of muon neutrino Maria Lorenzon July 22, 2010 The article I will mainly refer to is Observation of high-energy neutrino reactions and the existence of two kinds of neutrinos, written by G. Danby,

More information