Results from the NA61/SHINE experiment

Transcription

1 Results from the NA/SHINE experiment Alexis Häsler On Behalf of the NA/SHINE Collaboration July 3, Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

2 The NA/SHINE fixe target experiment at SPS CERN NA/SHINE SPS, CERN 3 participants from 5 institutions University of Belgrae, Belgrae, Serbia ETH, Zurich, Switzerlan Fachhochschule Frankfurt, Frankfurt, Germany Faculty of Physics, University of Sofia, Sofia, Bulgaria Karlsruhe Institute of technology, Karlsruhe, Germany Institute for Nuclear Research, Moscow, Russia Institute for Particle an Nuclear Stuies, KEK, Tsukuba, Japan Jagiellonian University, Cracow, Polan Joint Institute for Nuclear Research, Dubna, Russia University of Bern, Bern, Switzerlan LPNHE, University of Paris VI an VII, Paris, France University of Silesia, Katowice, Polan Rujer Boskovic Institute, Zagreb, Croatia National Center for Nuclear Research, Warsaw, Polan St Petersburg State University, St Petersburg, Russia University of Geneva, Geneva, Switzerlan Jan Kochanowski University in Kielce, Polan University of Athens, Athens, Greece University of Bergen, Bergen, Norway University of Frankfurt, Frankfurt, Germany University of Wroclaw, Wroclaw, Polan Faculty of Physics, University of Warsaw, Warsaw, Polan Warsaw University of Technology, Warsaw, Polan Laboratory of Astroparticle Physics, University Nova Gorica, Nova Gorica, Slovenia Wigner Research Centre for Physics of the Hungarian Acaemy of Sciences, Buapest, Hungary Alexis Ha sler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

3 The NA/SHINE experiment NA/SHINE : haron uction experiment at SPS CERN Rich physics program covering: heavy ion physics haron-uction measurements for cosmic ray experiments haron-uction measurements for accelerator neutrino experiments Large acceptance spectrometer: 5 TPCs ipole magnets (p)/p (GeV /c) (at 9Tm) (E/x)/ < E/x > 3 ToF (FTOF ) ps (TOF L/R) ps Data taken for the Neutrino Physics Program: Beam (GeV/c) graphite target year Nx p 3 cm 7 7 p 3 cm 9 5 p 3 9cm TK replica 7 p 3 9cm TK replica 9 p 3 9cm TK replica Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, 3 / 3

4 Haron Prouction measurements for the TK neutrino flux preiction Conventional ν beams uce by interactions of protons on fixe target Neutrino beam preictions rely on moeling the proton interactions an haron uction in the target Precise haron uction measurements allow to reuce uncertainties on neutrino flux preiction Examples:,%*># %*5# %*+#,%-%-/+3#!"#$%&#'(# -*+%-# /+3# )*+#,%-%-/+# 7;7%#,57# /9*:;3#*,#/<9*:;3=# *+*,-"!"!, "# µ# $ µ, $ ## /<9*:;3# /9*:;3# "!"#" $$%"#" &%!"#" &'(")#" HARP, CERN-PS 5-5 GeV Beam Mini- Sci- Micro- BooNE at Fermilab KK (KEK to Super-Kamiokane) NA & SPY/NA5, SPS - 5 GeV Beam WANF (NOMAD, CHORUS) CNGS (OPERA, ICARUS) NA9, CERN-SPS GeV Beam MIPP, FNAL-E97 GeV beam NA/SHINE, CERN-SPS 3- GeV Beam NuMi Beam Line in Fermilab MINOS, MINERvA, NOvA TK (Tokai to Super-Kamiokane) NuMI (MINOS+,MINERvA, NOvA) Future Long Base Line Experiments (LBNE, LBNO, THK) Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

5 Analysis Techniques in NA/SHINE Three analysis techniques: h analysis: analysis of π via measurements of negatively charge particles E/x [mip] 5 5 π - - K e - Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, 5 / 3

6 Analysis Techniques in NA/SHINE Three analysis techniques: h analysis: analysis of π via measurements of negatively charge particles E/x analysis at p GeV /c: π ± an protons were ientifie via energy loss in TPC E/x [mip] π + K p e + Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

7 hmp Entries 59 Mean x 9 Mean y 3 RMS x 3 RMS y 3 Analysis Techniques in NA/SHINE Three analysis techniques: h analysis: analysis of π via measurements of negatively charge particles E/x analysis at p GeV /c: π ± an protons were ientifie via energy loss in TPC ToF E/x analysis at p GeV /c: information from E/x an ToF is combine to ientify π ±, K ± an protons E/x [mip] 5 5 ] /c mmass versus versus momentum m [GeV 5 p K π 3 ] + π + K p e + /c [GeV m < < 3 ] /c [GeV m momentum [GeV/c] - 5 E/x [mip] Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, 7 / 3

8 Seven Particle Species Release for p+c at 3 GeV/c p + C π + + X p + C π + X < θ < mra < θ < mra 3 < θ < mra 3 < θ < mra < θ < mra < θ < mra 5 5 < θ < mra 5 5 < θ mra < 5 < θ < mra 5 < θ < mra 5 5 < mra < θ 5 5 θ mra < < < θ < mra 5 < θ < mra Data < θ < mra 5 < θ < mra Data < θ < 3 mra 5 3 < θ < 3 mra FTF_BIC - G95 FTF_BIC - G9 FTF_BIC - G QGSP_BERT - G 5 < θ < 3 mra 5 3 < θ < 3 mra FTF_BIC - G95 FTF_BIC - G9 FTF_BIC - G QGSP_BERT - G Typical uncertainty in region of interest for TK: % Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

9 Seven Particle Species Release for p+c at 3 GeV/c 5 p + C π + + X < θ < mra < θ < mra 5 5 < θ < mra < θ < mra < θ < mra < θ < mra 5 5 < θ < mra < θ < mra 5 K < θ < mra < θ < θ < mra < θ < mra 5 mra < ra GeV/c p + C π + X 3 < θ < mra 5 5 < θ < mra 5 5 < θ < mra 3 < θ < mra 3 5 K5 < θ < mra 5 5 < θ < mra < θ < mra 5 5 < θ < mra 5 5 < θ mra < < θ < mra 5 5 θ mra < < 5 5 θ mra < < < θ < mra < θ < 3 mra < θ < mra 5 < θ < mra Data < θ < 3 mra Data - 9 FTF_BIC - G95 < FTF_BIC θ < 3 mra - G9 < θ < 3 mra 3 FTF_BIC FTF_BIC - G95- G 5 FTF_BIC QGSP_BERT - G9 - G FTF_BIC - G QGSP_BERT - G < θ < mra < θ < mra 3 < mra < θ 3 < mra < θ < θ < 3 mra < θ < mra Data - 9 Data - 9 FTF_BIC - G95 FTF_BIC - G9 3 < θ < 3 mra 5 FTF_BIC - G FTF_BIC - G95 QGSP_BERT FTF_BIC - G - G9 FTF_BIC - G QGSP_BERT - G Typical uncertainty in region of interest for TK: 5% Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

10 Seven Particle Species Release for p+c at 3 GeV/c p + C π + + X < θ < mra < θ < mra < θ < mra < θ < mra 3 < θ < mra 5 5 < θ < mra < θ < mra p + C π < θ < mra + X < θ < mra 5 5 < θ < mra θ < mra θ < mra 5 3 < θ < 3 mra Data Data FTF_BIC - G95 FTF_BIC - G9 FTF_BIC - G < θ < 3 mra 3 < QGSP_BERT θ < 3 mra- G [GeV/c] FTF_BIC - G95 5 p FTF_BIC - G9 FTF_BIC - G QGSP_BERT - G 3 < θ < mra 5 5 < θ < mra 5 < K s 9 - stat FTF_BIC - G95 FTF_BIC - G9 FTF_BIC - G 5 QGSP_BERT - G 3 < θ < mra 3 < mra < θ 3 < mra < θ 3 < θ < mra 3 < θ < mra 5 5 < θ < mra < θ < mra 5 Data - 9 Data - 9 FTF_BIC - G95 FTF_BIC - G9 FTF_BIC - G 3 < θ < 3 mra 5 QGSP_BERT - G FTF_BIC - G95 FTF_BIC - G9 FTF_BIC - G QGSP_BERT - G Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

11 Seven Particle Species Release for p+c at 3 GeV/c p + C π + + X < θ < mra 5 5 < θ < mra < θ < mra < θ < mra p + C π + X < θ < mra < θ < mra < θ < mra < θ < mra 3< θ < mra 5 5 < θ < mra 3 < θ < mra < θ < mra 5 5 < θ < mra < < θ < mra 5 < θ < mra θ < mra < θ < mra < θ < mra < θ < mra K < θ < mra < θ < mra < θ < mra K K < θ < mra < θ < mra protons < θ < mra < θ < mra < θ < mra 5 < θ < mra < θ < mra mra 5 < θ < mra 5 5 θ < < [GeV/c] p mra < θ < mra θ < 3 < θ < mra < θ < 3 mra < mra < θ < mra < θ < mra < θ < mra < < θ < θ < mra 5 5 < θ < Data mra - 9 K s 9 - stat [GeV/c] Data - 9 p 3 < θ < 3 mra 3 < θ < 3 mra Data - 9 FTF_BIC - G95 < θ < 3 mra - G9 < FTF_BIC - G95 < 3 mra < θ < 3 mra 3 θ 3 < < θ < θ 3 < 3 mra mra 5 FTF_BIC FTF_BIC - G95- G FTF_BIC - G9 FTF_BIC - G95 FTF_BIC QGSP_BERT - G9 - G FTF_BIC - G9 [GeV/c] FTF_BIC - G FTF_BIC - G FTF_BIC - G 5 p QGSP_BERT - G QGSP_BERT - G ra GeV/c QGSP_BERT - G 3 3 < θ < mra 5 5 < θ < mra 5 < θ < mra < θ < mra Data - 9 Data - 9 FTF_BIC - G95 FTF_BIC - G9 3 < θ < 3 mra FTF_BIC - G FTF_BIC - G95 QGSP_BERT FTF_BIC - G - G9 5 FTF_BIC - G QGSP_BERT - G Use to tune seconary nucleon uction Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

12 Seven Particle Species Release for p+c at 3 GeV/c p + C π + + X < θ < mra 5 5 < θ < mra < θ < mra < θ < mra ra GeV/c p + C π + X < θ < mra < θ < mra < θ < mra < θ < mra 3< θ < mra 3 < θθ< < mra mra 5 < θ < 5 mra 5 < θ < mra < θ < mra < θ < mra < θ < mra < θ < mra < θ < mra < θ < mra K 5 + < θ < 5 mra 5 < θ < 5 mra < θ < mra K K < θ < mra5 5 < θ < 5 mra 5 < θ < mra protons 5 5 < θ < 55 mra 5 5 < θ < mra < θ < mra < θ < mra < θ < mra < θ < mra < θ < mra < θ < mra p 5[GeV/c] 55 5 p 5[GeV/c] < θ < mra < θ < 3 mra < θ < mra < θ < 3 mra Use to constrain fee own corrections 5 5 < θ < mra < θ < mra < θ < mra 5 < mra 5 5 < θ < θ < mra < θ < mra < θ < mra < θ < mra < θ < mra < θ < 5 mra < θ < mra 5 < θ < 3 mra < θ < mra < θ K s 9 - stat Data - 9 < mra Data - 9 [GeV/c] Data - 9 p 3 < θ < 3 mra 3 < θ < 3 mra Data - 9 Data - 9 FTF_BIC - G95 < FTF_BIC - G9 FTF_BIC < 3 mra - G95 FTF_BIC - G95 FTF_BIC - G9 θ 3 < θ < 3 mra < θ < 3 mra FTF_BIC QGSP_BERT - G95 [GeV/c] FTF_BIC - G9 G FTF_BIC - 3 G < θ < FTF_BIC - G9 FTF_BIC mra - G Λ 9 QGSP_BERT FTF_BIC - G - G95 FTF_BIC - G9 - stat 3 < θ < 3 mra < θ < 3 mra 3 < θ < 3 mra FTF_BIC - G FTF_BIC 5 - G9 FTF_BIC - G FTF_BIC - G - 5 p FTF_BIC - G QGSP_BERT - G FTF_BIC - G95 QGSP_BERT - G QGSP_BERT - G FTF_BIC - G QGSP_BERT - G - G9 FTF_BIC p FTF_BIC - G QGSP_BERT - G Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

13 The TK Beam Tuning with NA/SHINE Measurements p+c at 3 GeV/c ata in TK Major part of the TK phase space is covere by NA/SHINE ata Interaction chain is store at the simulation level to be tune later with measurements Tuning of tertiary particles requires extrapolation from NA/SHINE ata: Extrapolation to ifferent incient nucleon momenta Extrapolation to ifferent material (carbon to aluminum) Larger uncertainty to the flux ue to tertiary particle uction p Horn 9 Coverage 7 Coverage Arb Units only moel p+c π + + X [ra] θ π + 5 p [GeV/c] π + Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, 9 / 3

14 The TK Beam Tuning with NA/SHINE Measurements p+c at 3 GeV/c ata in TK [ra] θ π + Major part of the TK phase space is covere by NA/SHINE ata Interaction chain is store at the simulation level to be tune later with measurements Tuning of tertiary particles requires extrapolation from NA/SHINE ata: Extrapolation to ifferent incient nucleon momenta Extrapolation to ifferent material (carbon to aluminum) Larger uncertainty to the flux ue to tertiary particle uction 9 Coverage 7 Coverage 5 Arb Units Tuning weight Fractional Tuning weight Error E ν (GeV) - E ν E (GeV) ν (GeV) SK ν π + FIG : Ratio of e Flux the haron (TK) interaction KAbe et al, PhysRevD7 3 re-weighte (3) flux ove Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE VietnamTotal July e 3, (lower left), 9 / 3 e rror Pion Tuning Kaon Tuning Int Rate Tuning Total Tuning SK ν µ Flux 3 Pion Total Tuning Pion Prouction Kaon Kaon Tuning Prouction Seconary Nucleon Prouction Int Prouction Rate Tuning Cross Section Total Tuning Fractional Error rror

15 categories: the seconary component relates to neutrinos coming from ecays of pions uce in primary interactions of the protons in the target p + C π + + X; represents % of pions uce in the 9 cm long TK The TK Beam Tuning with NA/SHINE Measurements target; this component can be containe by measurements of pion an kaon cross-sections with a thin graphite target [, 5] Alternative approach to the neutrino flux preiction: the tertiary component relates to neutrinos coming from ecays of harons uce in interactions of x-z backwar positions of outgoing particles target surface extrap position onon tgt surface Measurement of haron multiplicities seconary particles either in the target at orthe out of the target (elements in the beam line) This component surface of a TK replica target place in the 3 represent % of the neutrino flux an can only be constraine through measurements of haron uction NA/SHINE experimental set-up with a replica of the TK target Analysis in (p, θ, z) bins ( longituinal z bins) Three set of measurements with a replica of the TK target have been taken in NA/SHINE M triggers for Re-weighting multiplicities of harons exiting the the 7 perio, M trigger for the 9 perio an M triggers for the perio target in the TK beam simulation - Re-weighting up to 9% as compare to % in - the stanar approach hxz_tgt_surface 59-9 x [cm] Entries Mean x Mean y RMS x RMS y 3 z -3 Relative contribution - - z z 3 z - - z - z z [cm] Seconary Haron Seconary+Tertiary Haron - In Target Alexis Ha sler (University of Geneva) Flavour Physics Conference / ICISE Vietnam 3 5 Eν (GeV) July 3, / 3

16 The TK Beam Tuning with NA/SHINE Measurements Pilot Analysis conucte on 7 ataset emonstrate a proof of principle to tune neutrino flux preiction with replica target measurements Due to very low statistics, final uncertainty on the neutrino flux preiction coul not be reuce with respect to the thin target tuning Weights were compute as ω NA π (p, θ, z) = n NA πata /nna πmc within the etector acceptance without any further corrections pot] / MeV/ flux [cm ν µ 5 replica target, pion multiplicity error thin target: total error pion multiplicity error E ν (GeV) thin/replica re weighting 5 χ /nf 399/ p ± E ν (GeV) Pion emission from the TK replica target: metho, results an application Nuclear Inst an Methos in Physics Research, A (3) pp 99- Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

17 The TK Beam Tuning with NA/SHINE Measurements Large statistics 9 ataset of the TK replica target measurements allow to extract precise pion spectra Monte Carlo base Corrections are applie to account for geometrical acceptance of the etector, particle loss, fee-own, Example of π + spectra for the six longituinal bins an two theta intervals as a function of momentum 3 <θ< <θ< n/p x /pot n/p x /pot Z Z <θ< <θ< <θ< <θ< <θ< Z Z 3 Z Z 5 Z NA stat&sys NA stat&sys FLUKA Spline FLUKA Spline 7 9 <θ< <θ< < θ< < θ< < θ< Z Z 3 Z Z 5 Z NA stat&sys NA stat&sys FLUKA Spline FLUKA Spline Statistical errors are of the orer of 5- % Systematic uncertainties of the orer of 5% for the center of the target, up to % for the most upstream part an ownstream face Ongoing work on the implementation of these ata in the TK flux simulation will allow to reuce the uncertainties on the preictions Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, / 3

18 Summary The NA/SHINE experiment has prove it s ability to eliver high quality ata use for Neutrino Physics programs All ata requeste by the TK long baseline experiment have been recore 7 particle species are extracte for p+c at 3 GeV/c: π ±, K ±, K s, protons an Λ release 3 articles have been publishe from p+c at 3 GeV/c Measurement of Prouction Properties of Positively Charge Kaons in Proton-Carbon Interactions at 3 GeV/c; PhysRev C5 () Measurements of Cross Sections an Charge Pion Spectra in Proton-Carbon Interactions at 3 GeV/c; PhysRev C () Pion emission from the TK replica target: metho, results an application; Nuclear Inst an Methos in Physics Research, A (3) NA/SHINE publishe results have been extensively use an have contribute to the recent achievements of the TK experiment Detaile explanations have been publishe in an article by the TK collaboration (TK neutrino flux preiction; Phys Rev D 7,(3)) Future results with the thin (cm) an TK replica graphite target analyses will allow to further improve the TK neutrino flux preiction NA/SHINE plans to continue taking haron uction measurements for future neutrino program at Fermilab an for long base line neutrino experiments Alexis Häsler (University of Geneva) Flavour Physics Conference / ICISE Vietnam July 3, 3 / 3