Measurements of Cross Sections and Charged Pion Spectra in Proton-Carbon Interactions at 31 GeV/c

Measurements of Cross Sections and Charged Pion Spectra in Proton-Carbon Interactions at 31 GeV/c

L. S. Esposito, on behalf of the NA61/SHINE Collaboration    The NA61 Collaboration:
N. Abgrall, A. Aduszkiewicz, B. Andrieu, T. Anticic, N. Antoniou, J. Argyriades, A. G. Asryan, B. Baatar, A. Blondel, J. Blumer, M. Bogusz, L. Boldizsar, A. Bravar, W. Brooks, J. Brzychczyk, A. Bubak, S. A. Bunyatov, O. Busygina, T. Cetner, K.-U. Choi, P. Christakoglou, P. Chung, T. Czopowicz, N. Davis, F. Diakonos, S. Di Luise, W. Dominik, J. Dumarchez, R. Engel, A. Ereditato, L. S. Esposito, G. A. Feofilov, Z. Fodor, A. Ferrero, A. Fulop, X. Garrido, M. Gaździcki,, M. Golubeva, K. Grebieszkow, A. Grzeszczuk, F. Guber, H. Hakobyan, T. Hasegawa, S. Igolkin, A. S. Ivanov, Y. Ivanov, A. Ivashkin, K. Kadija, A. Kapoyannis, N. Katryńska, D. Kiełczewska, D. Kikola, J.-H. Kim, M. Kirejczyk, J. Kisiel, T. Kobayashi, O. Kochebina, V. I. Kolesnikov, D. Kolev, V. P. Kondratiev, A. Korzenev, S. Kowalski, S. Kuleshov, A. Kurepin, R. Lacey, J. Lagoda, A. Laszlo, V. V. Lyubushkin, M. Mackowiak, Z. Majka, A. I. Malakhov, A. Marchionni, A. Marcinek, I. Maris, V. Marin, T. Matulewicz, V. Matveev, G. L. Melkumov, A. Meregaglia, M. Messina, St. Mrówczyński, S. Murphy, T. Nakadaira, P. A. Naumenko, K. Nishikawa, T. Palczewski, G. Palla, A. D. Panagiotou, W. Peryt, O. Petukhov, R. Płaneta, J. Pluta, B. A. Popov,, M. Posiadała, S. Puławski, W. Rauch, M. Ravonel, R. Renfordt, A. Robert, D. Röhrich, E. Rondio, B. Rossi, M. Roth, A. Rubbia, M. Rybczyński, A. Sadovsky, K. Sakashita, T. Sekiguchi, P. Seyboth, M. Shibata, A. N. Sissakian, E. Skrzypczak, M. Słodkowski, A. S. Sorin, P. Staszel, G. Stefanek, J. Stepaniak, C. Strabel, H. Ströbele, T. Susa, P. Szaflik, M. Szuba, M. Tada, A. Taranenko, R. Tsenov, R. Ulrich, M. Unger, M. Vassiliou, V. V. Vechernin, G. Vesztergombi, A. Wilczek, Z. Włodarczyk, A. Wojtaszek, J.-G. Yi, I.-K. Yoo, W. Zipper
   ETH, Zurich, Switzerland
Fachhochschule Frankfurt, Frankfurt, Germany
Faculty of Physics, University of Sofia, Sofia, Bulgaria
Faculty of Physics, University of Warsaw, Warsaw, Poland
High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan
Institute for Nuclear Research, Moscow, Russia
Jagiellonian University, Cracow, Poland
Jan Kochanowski University in Kielce, Poland
Joint Institute for Nuclear Research, Dubna, Russia
KFKI Research Institute for Particle and Nuclear Physics, Budapest, Hungary
Karlsruhe Institute of Technology, Karlsruhe, Germany
LPNHE, University of Paris VI and VII, Paris, France
Pusan National University, Pusan, Republic of Korea
Rudjer Boskovic Institute, Zagreb, Croatia
Soltan Institute for Nuclear Studies, Warsaw, Poland
St. Petersburg State University, St. Petersburg, Russia
State University of New York, Stony Brook, USA
The Universidad Tecnica Federico Santa Maria, Valparaiso, Chile
University of Athens, Athens, Greece
University of Bergen, Bergen, Norway
University of Bern, Bern, Switzerland
University of Frankfurt, Frankfurt, Germany
University of Geneva, Geneva, Switzerland
University of Silesia, Katowice, Poland
Warsaw University of Technology, Warsaw, Poland
Present affiliation: University of Wrocław, Wrocław, Poland

As neutrino long baseline experiments enter a new domain of precision, the careful study of systematic errors due to poor knowledge of production cross sections for pions and kaons require more dedicated measurements for precise neutrino flux predictions. The cosmic ray experiments require dedicated hadron production measurements to tune simulation models used to describe air shower profiles. Among other goals, the NA61/SHINE (SPS Heavy Ion and Neutrino Experiment) experiment at the CERN SPS aims at precision measurements (5% and below) for both neutrino and cosmic ray experiments: those will improve the prediction of the neutrino flux for the T2K experiment at J-PARC and the prediction of muon production in the propagation of air showers for the Auger and KASCADE experiments. NA61/SHINE took data during a pilot run in 2007 and in 2009 and 2010 with different carbon targets. The NA61/SHINE set-up and spectra for positive and negative pions obtained with the 2007 thin (4% interaction length) carbon target data are presented .

ETH, Institute for Particle Physics,
Zurich, Switzerland

Keywords: p+C interaction, inelastic cross section, inclusive pion spectra

1 Physics motivation

The NA61/SHINE (SPS Heavy Ion and Neutrino Experiment) experiment at the CERN SPS pursues a rich physics program in various fields . First, precise hadron production measurements are performed for improving calculations of the neutrino flux in the T2K neutrino oscillation experiment , as well as for more reliable simulations of cosmic-ray air showers in the Pierre Auger and KASCADE experiments . Second, p+p, p+Pb and nucleus+nucleus collisions will be studied extensively at SPS energies. This article presents first NA61/SHINE results on charged pion spectra in p+C interactions at 31 GeV/c which are needed for an accurate neutrino flux prediction in the T2K experiment. The results are based on the data collected during the first NA61/SHINE run in 2007.

T2K is a long baseline neutrino experiment in Japan, which uses a high intensity neutrino beam produced at J-PARCaaaJapan Proton Accelerator Research Complex organized jointly by JAEA and KEK in Tokai, Japan. It aims to precisely measure the appearance and disappearance . In order to generate the neutrino beam a high intensity 30 GeV (kinetic energy) proton beam impinging on a 90 cm long graphite target is used, where and K mesons decaying into (anti)neutrinos are produced. The neutrino fluxes and spectra are then measured both at the near detector complex, 280 m from the target, and by the Super-Kamiokande (SK) detector located 295 km away from the neutrino source and 2.5 degrees off-axis. Neutrino oscillations are probed by comparing the neutrino flux measured at SK to the predicted one. In order to predict the flux at SK one uses the near detector measurements and extrapolates them to SK with the help of Monte Carlo simulations. Up to now, these Monte Carlo predictions are based on hadron production models only. For more precise predictions, which would allow the reduction of systematic uncertainties to the level needed for the T2K physics goals, measurements of pion and kaon production off carbon targets are essential . The purpose of the NA61/SHINE measurements for T2K is to provide this information at exactly the proton extraction energy of the J-PARC Main Ring synchrotron, namely 30 GeV kinetic energy (approximately 31 GeV/c momentum). Presently, the T2K neutrino beam-line is set up to focus positively charged hadrons, in such a way that it produces a beam. Spectra of positively charged pions presented in this paper constitute directly an essential ingredient in the neutrino flux calculation.

2 The NA61/SHINE detector

The NA61/SHINE experiment is a large acceptance hadron spectrometer in the North Area H2 beam-line of the CERN SPS. The schematic layout is shown in Fig. 1 together with the overall dimensions.

Fig. 1: (Color online) The layout of the NA61/SHINE experiment at the CERN SPS (top view, not to scale). The chosen right-handed coordinate system is shown on the plot. The incoming beam direction is along the axis. The magnetic field bends charged particle trajectories in the (horizontal) plane. The drift direction in the TPCs is along the (vertical) axis.

The main components of the current detector were constructed and used by the NA49 collaboration . A set of scintillation and Cherenkov counters as well as beam position detectors (BPDs) upstream of the spectrometer provide timing reference, identification and position measurements of the incoming beam particles. The main tracking devices of the spectrometer are large volume Time Projection Chambers (TPCs). Two of them, the vertex TPCs (VTPC-1 and VTPC-2 in Fig. 1), are located in a free gap of 100 cm between the upper and lower coils of the two superconducting dipole magnets. Their maximum combined bending power is 9 Tm. In order to optimize the acceptance of the detector at 31 GeV/c beam momentum, the magnetic field used during the 2007 data taking period was set to a bending power of 1.14 Tm. Two large TPCs (MTPC-L and MTPC-R) are positioned downstream of the magnets symmetrically to the beam line. The TPCs are filled with Ar:CO gas mixtures in proportions 90:10 for VTPCs and 95:5 for MTPCs. The particle identification capability of the TPCs based on measurements of the specific energy loss, , is augmented by time-of-flight measurements using Time-of-Flight (ToF) detectors. The ToF-L and ToF-R arrays of scintillator pixels have a time resolution of better than 90 ps . Before the 2007 run the experiment was upgraded with a new forward time-of-flight detector (ToF-F) in order to extend the acceptance. The ToF-F consists of 64 scintillator bars with photomultiplier (PMT) readout at both ends resulting in a time resolution of about 115 ps. The target under study is installed 80 cm in front of the VTPC-1. The results presented here were obtained with an isotropic graphite target of dimensions 2.5(W)2.5(H)2(L) cm and with a density of  g/cm. The target thickness along the beam is equivalent to about 4% of a nuclear interaction length ().

3 Analysis techniques

This section presents the procedures used for data analysis to extract pion cross sections. Crucial for this analysis is the identification of the produced pions. Depending on the momentum interval, different approaches have been adopted, which lead also to different track selection criteria. The calibrated distributions as a function of particle momentum for positively and negatively charged particles are presented in Fig. 2.

Fig. 2: (Color online) Specific energy loss in the TPCs for positively () and negatively () charged particles as a function of momentum. Curves show parameterizations of the mean calculated for different particle species.

The task is facilitated for the negatively charged pions, by the observation that more than 90% of primary negatively charged particles produced in p+C interactions at this energy are , and thus the analysis of spectra can also be carried out without additional particle identification.

In the low momentum region (less than about 1 GeV/c), it is sufficient to distinguish pions from electrons/positrons, kaons and protons by means of particle identification via measurements of specific energy loss () in the TPCs. A reliable identification of mesons was not possible at momenta above 1 GeV/c where the Bethe-Bloch (BB) curves for pions, kaons and protons cross each other. On the other hand, for mesons, where the contribution of and antiprotons is almost negligible, the analysis could be extended in momentum up to 3 GeV/c allowing consistency checks with the other analysis methods in the region of overlap.

Fig. 3: (Color online) Right: Mass squared, derived from the ToF-F measurement and the fitted path length and momentum, versus momentum . The lines show the expected mass squared values for different particles. Left: Example of two-dimensional plots for positively charged particles in the momentum range 2-3 GeV/c. contour around fitted pion peaks are shown.

High purity particle identification can be performed by combining the and information. Moreover, in the momentum range 1-4 GeV/c, where bands for different particle species overlap, particle identification is in general only possible using the method, see Figs. 3. The ToF-F detector was designed to cover the necessary acceptance in momentum and polar angle required by the T2K experiment, although limited to particle momenta above about 0.8 GeV/c.

Indeed, three analysis methods were applied to obtain pion spectra:

  1. [(i)]

  2. analysis of mesons via measurements of negatively charged particles ( analysis

  3. analysis of and mesons identified via measurements in the TPCs ( analysis at low momentum and

  4. analysis of and mesons identified via time-of-flight and measurements in the ToF-F and TPCs, respectively ( analysis.

Each analysis yields fully corrected pion spectra with independently calculated statistical and systematic errors. The spectra were compared in overlapping phase-space domains to check their consistency. Complementary domains were combined to reach maximum acceptance.

The agreement between the spectra obtained by different methods is, in general, better than 10%. Note, that data points in the same () bin from different analysis methods are statistically correlated as they result from the analysis of the same data set. In order to obtain the final spectra consisting of statistically uncorrelated points the measurement with the smallest total error was selected.

The analysis was based on a sample of 521 k events selected from the total sample of 667 k of registered and reconstructed events. This criterion essentially removes a contamination by interactions upstream of the target. For the event sample with the target removed, the selection reduces the number of events from 46 k to 17 k.

4 Results

This section presents results on inelastic and production cross sections as well as on differential spectra of and mesons in p+C interactions at 31 GeV/c.

4.1 Inelastic and production cross sections

The total inelastic cross section is measured to be 

The production cross section was calculated from the inelastic cross section by subtracting the quasi-elastic contribution . The result is :

The production cross section is compared to previous measurements in Fig. 4.

Fig. 4: (Color online) Beam momentum dependence of the production cross section for p+C interactions. The NA61/SHINE result (filled square) is compared with previous measurements: Bellettini et al. (circle) , Carroll et al. (triangle)  and Denisov et al. (open squares) . For the NA61/SHINE point, the error bar indicates statistical and systematic uncertainties added in quadrature. The result from Ref.  was recalculated by subtracting from the measured inelastic cross section a quasi-elastic contribution at 20 GeV/c of () mb.

4.2 Spectra of and mesons

Fig. 5: (Color online) Laboratory momentum distributions of mesons produced in production p+C interactions at 31 GeV/c in different intervals of polar angle (). The spectra are normalized to the mean multiplicity in all production p+C interactions. Error bars indicate statistical and systematic uncertainties added in quadrature. The overall uncertainty () due to the normalization procedure is not shown. Predictions of hadron production models, Fluka2008  (solid line), Urqmd1.3.1  (dashed line) and Venus4.12  (dotted line) are also indicated.
Fig. 6: (Color online) Laboratory momentum distributions of mesons produced in production p+C interactions at 31 GeV/c in different intervals of polar angle (). The spectra are normalized to the mean multiplicity in all production p+C interactions. Error bars indicate statistical and systematic uncertainties added in quadrature. The overall uncertainty () due to the normalization procedure is not shown. Predictions of hadron production models, Fluka2008  (solid line), Urqmd1.3.1  (dashed line) and Venus4.12  (dotted line) are also indicated.

The and spectra presented in this section refer to pions produced in strong and electromagnetic processes in p+C interactions at 31 GeV/c.

The spectra are presented as a function of particle momentum in ten intervals of the polar angle. Both quantities are calculated in the laboratory system. The chosen binning takes into account the available statistics of the 2007 data sample, detector acceptance and particle production kinematics.

The final spectra are plotted in Figs. 5 and  6. For the purpose of a comparison of the data with model predictions the spectra were normalized to the mean multiplicity in all production interactions by dividing by . This avoids uncertainties due to the different treatment of quasi-elastic interactions in models as well as problems due to the absence of predictions for inclusive cross sections.

5 Conclusion

This work presents inelastic and production cross sections as well as positively and negatively charged pion spectra in p+C interactions at 31 GeV/c on the 2007 thin (4% interaction length) carbon target. These data are essential for precise predictions of the neutrino flux for the T2K long baseline neutrino oscillation experiment in Japan. Furthermore, they provide important input to improve hadron production models needed for the interpretation of air showers initiated by ultra high energy cosmic particles.

A much larger data set with both the thin and the T2K replica carbon targets was recorded in 2009 and 2010 and is presently being analysed. This will lead to results of higher precision for pions and extend the measurements to other hadron species such as charged kaons, protons, and . Analysis of the data collected with the T2K replica target during the 2007 run is in progress . The new data will allow a further significant reduction of the uncertainties in the prediction of the neutrino flux in the T2K experiment.

6 Acknowledgments

This work was supported by the following funding agencies: the Hungarian Scientific Research Fund (grants OTKA 68506 and 79840), the Polish Ministry of Science and Higher Education (grants 667/N-CERN/2010/0, N N202 1267 36, N N202 287838 (PBP 2878/B/H03/2010/38), DWM/57/T2K/2007), the Federal Agency of Education of the Ministry of Education and Science of the Russian Federation (grant RNP, the Russian Academy of Sciences and the Russian Foundation for Basic Research (grants 08-02-00018 and 09-02-00664), the Ministry of Education, Culture, Sports, Science and Technology, Japan, Grant-in-Aid for Scientific Research (grants 18071005, 19034011, 19740162, 20740160 and 20039012), the Toshiko Yuasa Lab. (France-Japan Particle Physics Laboratory), the Institut National de Physique Nucléaire et Physique des Particules (IN2P3, France), the German Research Foundation (grant GA 1480/2-1), the Swiss National Science Foundation (Investigator-Driven projects and SINERGIA) and the Swiss State Secretariat for Education and Research (FORCE grants).
The authors also wish to acknowledge the support provided by the collaborating institutions, in particular, the ETH Zurich (Research Grant TH-01 07-3), the University of Bern and the University of Geneva.

Finally, it is a pleasure to thank the European Organization for Nuclear Research for a strong support and hospitality and, in particular, the operating crews of the CERN SPS accelerator and beam lines who made the measurements possible.


  • 1. NAbgrall et al. [ NA61/SHINE Collaboration ], “Measurements of Cross Sections and Charged Pion Spectra in Proton-Carbon Interactions at 31 GeV/c,” Submitted to: Phys.Rev.C. [arXiv:1102.0983 [hep-ex]].
  • 2. N. Antoniou et al. [NA49-future Collaboration], CERN-SPSC-2006-034, (2006).
  • 3. N. Antoniou et al. [NA61/SHINE Collaboration], CERN-SPSC-2007-004, (2007).
  • 4. N. Antoniou et al. [NA61/SHINE Collaboration], CERN-SPSC-2007-019, (2007).
  • 5. N. Abgrall et al. [NA61/SHINE Collaboration], CERN-SPSC-2008-018, (2008).
  • 6. Y. Itow et al. The JHF-Kamiokande neutrino project, arXiv:hep-ex/0106019, (2001).
  • 7. J. Abraham et al. [Pierre Auger Collaboration], Nucl. Instrum. Meth. A 523, 50 (2004).
  • 8. T. Antoni et al. [KASCADE Collaboration], Nucl. Instrum. Meth. A 513, 490 (2003).
  • 9. S. Afanasiev et al. [NA49 Collaboration], Nucl. Instrum. Meth. A 430, 257 (1999).
  • 10. T. Palczewski, Ph.D. Thesis in preparation, Soltan Institute for Nuclear Studies, 2011, Warsaw Poland.
  • 11. M. Posiadala, Ph.D. Thesis in preparation, University of Warsaw, 2011, Warsaw Poland.
  • 12. S. Murphy, Ph.D. Thesis in preparation, University of Geneva, 2011, Geneva Switzerland.
  • 13. C. Strabel, Ph.D. Thesis, ETH, January 2011, Zurich, Switzerland.
  • 14. R. J. Glauber and G. Matthiae, Nucl. Phys. B 21, 135 (1970).
  • 15. G. Bellettini et al., Nucl. Phys. 79, 609 (1966).
  • 16. A. S. Carroll et al., Phys. Lett. B 80, 319 (1979).
  • 17. S. P. Denisov et al., Nucl. Phys. B 61, 62 (1973).
  • 18. A. Fasso et al., CERN-2005-10 (2005); G. Battistoni et al., AIP Conf. Proc. 896, 31, (2007).
  • 19. S.M. Bass et al., Prog. Part. Nucl. Phys. 41 225 (1998); M. Bleicher et al., J. Phys. G: Nucl. Part. Phys. 25, 1859 (1999).
  • 20. K. Werner, Nucl. Phys. A525, 501c (1991);
    K. Werner, Phys. Rep. 232, 87 (1993).
  • 21. N. Abgrall, Ph.D. Thesis in preparation, University of Geneva, 2011, Geneva Switzerland.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description