EPJ Web of Conferences 70, 00 033 (2014) DOI: 10.1051/epj conf/20147 000 033 C Owned by the authors, published by EDP Sciences, 2014

The Old New Frontier: Studying the CERN-SPS Energy Range with NA61/SHINE

Marek Szuba1,a for the NA61/SHINE Collaboration2

1Karlsruhe Institute of Technology, Germany 2CERN, Geneva, Switzerland

Abstract. With the Large entering its third year of granting us insight into the highest collision energies to date, one should nevertheless keep in mind the un- explored physics potential of lower energies. A prime example here is the NA61/SHINE experiment at the CERN Super . Using its large-acceptance hadronic spectrometer, SHINE aims to accomplish a number of physics goals: measuring spec- tra of identified hadrons in hadron-nucleus collisions to provide reference for accelerator experiments and cosmic-ray observatories, investigating particle properties in the large transverse-momentum range for hadron+hadron and hadron+nucleus collisions for studying the nuclear modification factor at SPS energies, and measuring hadronic ob- servables in a particularly interesting region of the phase diagram of strongly-interacting matter to study the onset of deconfinement and search for the critical point of strongly- interacting matter with nucleus-nucleus collisions. This contribution shall summarise results obtained so far by NA61/SHINE, as well as present the current status and plans of its experimental programme.

1 Introduction

There is no doubt that studies of particle collisions at the highest available energies can and do result in important scientific results; the string of discoveries from the at and the recent observation at the at CERN of what is likely the Higgs boson attest to that. Even so, one should keep in mind the unexplored physics potential of collisions at lower energies. Data from experiments such as those at the Super (SPS) at CERN, operating well below the high-energy frontier, continue to contribute to expanding our knowledge of matter, our world and the universe. Examples of projects depending on such data include:

1.1 Measurement of Neutrino Oscillation

Predicted by Pontecorvo 1957 and since then observed experimentally, neutrino oscillation is a quantum-mechanical phenomenon which can make a neutrino produced with a specific flavour be observed later as possessing different flavour [1]. It stems from differences in quantum-phase prop- agation of different neutrino mass eigenstates and as such requires to have mass. One way

ae-mail: [email protected]

This is an Open Access article distributed under the terms of the Creative Commons Attribution License 2.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Article available at http://www.epj-conferences.org or http://dx.doi.org/10.1051/epjconf/20147000033 EPJ Web of Conferences to obtain accurate measurements of neutrino oscillation is to perform it on neutrinos produced in a controlled environment of a : a high-energy beam of is collided against a graphite target to produce positive pions and kaons, which subsequently decay into muons and muon neutrinos. One experiment measuring oscillation of beam neutrinos is T2K in Japan, directing particles from the J-PARC facility in Tokai towards the Super-Kamiokande detector 295 km away in Kamioka [2]. Physics goals of T2K are to obtain one of the first measurements of the θ13 mixing angle of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino-mixing matrix, to improve precision of mea- θ ff Δ 2 surements of the 23 mixing angle and the mass di erence m23, and in the future to search for ν CP violation. Other experiments of this sort include MINOS and NOνA (NuMI beam at Fer- milab), LBNE (upcoming new ν beam at Fermilab), and OPERA and ICARUS (CNGS beam at CERN) [3][4][5][6][7]. Unfortunately even with well-defined primary beams, estimating the resulting neutrino flux is not a trivial matter. In order to achieve adequate beam intensity the graphite target must be long, leading to secondary interactions. Moreover, a non-negligible contribution to the flux is made by interactions with support structures such as the target holder and cooling systems. In short, good knowledge of properties of the resulting neutrino beam requires good knowledge of hadron production in the target — which in turn may be much easier to study using an identical target at a fixed-target facility oriented towards hadron spectroscopy than having to set up necessary detectors in situ. Given the energies of proton beams typically used for production of neutrino beams, the CERN SPS complex (indeed, in case of the CNGS beam the proton accelerator in question is the SPS) is a perfect candidate for this task.

1.2 Hadron Production in Extensive Air Showers

Another domain which can benefit from hadroproduction measurements in the energy range of the SPS are studies of extensive air showers (EAS) produced in Earth’s atmosphere by cosmic rays. De- tecting and measuring EAS is the standard technique for studying ultra high-energy cosmic rays, as such particles reach our planet so infrequently and from so many different directions that it is vir- tually impossible to gather a statistically significant sample through direct measurements (i.e. using satellite- or balloon-based detectors to observe cosmic-ray particles themselves); using EAS detectors makes it both technically and financially easier to cover a much wider area. Examples of modern EAS experiments include KASCADE-Grande in Germany, Pierre Auger Observatory in Argentina and the Telescope Array in the United States [8][9][1]. Improved observation rates of EAS-based cosmic-ray experiments come at a cost of having to accurately reconstruct properties of the original particle from what has undergone multiple levels of interactions with the atmosphere before reaching the detectors. In particular, evolution of the hadronic component of showers remains relatively poorly understood. With primary energies of ultra high-energy cosmic rays remaining largely out of range of man-made accelerators, one is required to employ simulations to provide reference for such evolution — which in turn requires careful tuning of models used for this purpose. This is where measuring hadron spectra at the SPS comes in. Although its energy range falls 15 several orders of magnitude short of the energy of cosmic rays observed by EAS detectors (E0 = 10 − 1020 eV, it provides an excellent match to the energy of the last generation of hadronic interactions in the shower. As many experimental observables (e.g. the number of muons detected by surface detectors) are directly tied to this stage of shower evolution, studying it at the SPS is expected to yield significantly improved predictions.

00033-p.2 ICFP 2012

1.3 High-pT Hadrons in p+p and p+A Colliisions

Proton-proton and proton-nucleus collisions constitute important reference systems for a wide range of different studies (for instance spectra, scaling with the number of wounded nucleons or binary collisions, nuclear modification factor and Cronin effect, among others) in nucleus-nucleus reactions. However, in the past data sets of this sort at SPS energies and below were relatively small. This prob- lem becomes particularly pressing in case of high-pT hadrons, originating from hard scatterings and thus useful for probing the perturbative-quantum chromodynamics (QCD) regime of strong interac- tions — they are produced so rarely that most such elementary collisions contain none at all, further reducing effective sample size. Fortunately, development of sophisticated detectors and readout com- ponents capable of coping with high event rates of high-luminosity machines such as the LHC has made it possible to also improve event rates of lower-energy detectors, allowing at last collection of large-statistics p+p and p+A data sets. Analysis of this data is expected to greatly improve precision of many existing results as well as allow new, hitherto-infeasible studies.

1.4 Critical Point and the Onset of Deconfinement of QCD Matter

Theoretical predictions based on QCD tell us that above certain energy density, quarks and gluons previously confined to hadrons can undergo a phase transition into the quark-gluon plasma — a state in which they can be considered free. The energy range of the SPS is very important for studying this transition because at these energies one can produce collisions both right below and right above the energy threshold for deconfinement [1]. Moreover, it is also in this energy range that we now expect the critical point of the QCD phase diagram [1]. A number of experimental observables were proposed to be sensitive to the onset of deconfine- ment; these include changes in energy dependence of the pion yield per wounded nucleon (“the kink”), K+ of the π+ ratio (“the horn”) and of the mean transverse mass (“the step”) [1]. On the other hand, the emergence of the critical point is predicted to be observable in event-by-event fluctuations of e.g. multiplicity and average transverse momentum [1]. Indications in favour of this hypothesis has been observed both at the SPS and elsewhere [1][1][1].

2 The NA61/SHINE Experiment

2.1 Overview and Physics Goals

NA61/SHINE (SPS Heavy Ion and Nutrino Experiment) is a fixed-target experiment located in the North Area of the CERN SPS, using a large-acceptance hadronic spectrometer to study a wide range of phenomena in a number of different hadron-hadron, hadron-nucleus and nucleus-nucleus reac- tions [1]. It is the successor of the NA49 experiment, which took data in the years 1994–2002, and reuses most of its predecessor’s hardware and software [1]. Its large acceptance (around 50 % for 2 −4 −1 pT ≤ 2.5 GeV/c), high momentum resolution (σ(p)/p ≈ 10 (GeV/c) ) and tracking efficiency σ dE / dE ≈ ,σ ≈ (over 95 %), and excellent particle-identification capabilities ( ( dx ) dx 4% (tToF) 100 ps) make it an excellent tool for investigating hadron spectra. NA61/SHINE has the following physics goals: • search for the critical point on the phase diagram of strongly-interacting (QCD) matter and study in detail the onset of deconfinement by performing a two-dimensional scan of the phase diagram (see Figure 1),

• investigate production of high-pT hadrons in p+p and p+Pb collisions,

00033-p.3 EPJ Web of Conferences

200 158A GeV 13A GeV

180 NA61 future data T (MeV) p+p NA61 collected data C+C 160 Si+Si Pb+Pb p+p 140 NA49 Be+Be Ar+Ca top low Xe+La 120 SPS SPS

100 200 300 400 500 600 μ (MeV) B

Figure 1. Estimated [1] (squares) and extrapolated (circles) chemical freeze-out points corresponding to NA49 and NA61/SHINE for both collected and future data.

• provide reference measurements for the T2K neutrino experiment by measuring hadron production in collisions of a proton beam at 31 GeV with a carbon platelet and a replica of the T2K target, • provide reference measurements for current and future Fermilab neutrino experiments by measuring hadron production in collisions of a proton beam at 60, 90 and 120 GeV with a carbon platelet and a replica of the NuMI target, • perform reference measurements of hadron production in p+p, p+C and π+C interactions, for cosmic-ray extended air-shower experiments: KASCADE, KASCADE-Grande and Pierre Auger Observatory.

2.2 Detector Set-up

Figure 2 shows a diagram of the NA61/SHINE apparatus. The primary, tracking detectors are four large-volume Time Projection Chambers (TPCs) inherited from NA49. Two of them (Vertex TPCs) are located inside superconducting magnets along the beam line, the other two (Main TPCs) are placed symmetrically to the beam line further downstream. Another, smaller chamber (Gap TPC) is placed between the VTPCs. Together, these chambers measure trajectories and momentum of particles as well as allow their identification through ionisation energy loss. Just downstream of the MTPCs there exist three Time-of-Flight (ToF) walls, the side two inherited from NA49 and the central one added by NA61 in 2007. Measurements from these detectors comple- dE ment particle-identification information in the momentum range where dx alone is ambiguous.

00033-p.4 ICFP 2012

 

   

 

  

                

          

     !    # "

Figure 2. Layout of the NA61/SHINE experimental set-up (top view, not to scale). For details see [1].

The most downstream detector of the apparatus is the Projectile Spectator Detector (PSD), in- stalled in 2011. The purpose of this calorimeter is to provide precise, high-granularity event-by-event measurements of the energy of non-interacting fragments of the projectile nuclei, making it possible to determine centrality of the collision and the orientation of the reaction plane. Finally, a number of additional detectors installed in the beam line both upstream and downstream of the target monitor beam properties and provide trigger information. It is also worth mentioning at this point that a sophisticated beam line set-up was designed and put in place by the CERN accelerators-and-beams department in collaboration with NA61 for the purpose of providing beryllium beams for SHINE. At the moment the only primary beams the SPS can provide are protons and lead ions, with argon and xenon to be made possible during the long shutdown of the accelerator complex in 2013-2014. As adding new species of primary beams is an immensely complicated task, beryllium beams for SHINE were instead produced in a fragmentation- ion beam line. Primary lead ions from the SPS were collided with a beryllium block far upstream in the beam line and the resulting fragments of the projectile were guided through a sequence of magnets, collimators, degraders and beam detectors in order to filter out beryllium ions. This configuration was proven to produce highly pure, stable beams of 7Be over a wide range of Pb-beam energies.

2.3 Data Taking

NA61/SHINE has been taking data since the year 2007 and so far has acquired over 150 million events in a variety of data sets. We focused primarily on data for cosmic-ray and neutrino physics in the first three years of running, and progressed to high-pT data and the exploration of the phase diagram in the following years. The presently-approved physics programme is foreseen to be concluded in the year 2015.

00033-p.5 EPJ Web of Conferences

    θ θ

σ

!"# $$ %&' $ σ   θ θ

  θ θ

σ

!"# $$ %&' $ σ   θ θ

(( ((

Figure 3. Inclusive spectra of positive (top) and negative (bottom) pions, in several polar-angle bins, measured by NA61/SHINE in p+C collisions at 31 GeV (points), compared to predictions (lines) from FLUKA 2008, UrQMD 1.3.1 and VENUS 4.12 [2].

3 Results So Far

The first π± (see Figure 3) and K+ (Figure 4) spectra from p+C interactions at 31 GeV, based on data acquired in 2007, are published [2][2]. These results have already been used by the T2K experiment as reference for calculations of neutrino flux in the J-PARC beam, leading to their first measure- ment of the θ13 ν mixing angle [2]. Moreover, comparison of pion spectra with models have made it possible to introduce improvements to UrQMD as well as to the implementation of FRITIOF in the toolkit [2][2]. Results from data taken in 2009, offering a factor-of-ten increase in statistics as well improved particle identification through extension of the forward Time-of-Flight wall, are being prepared for release.

00033-p.6 ICFP 2012



θ θ             π θ     θ         σ     

 

σ    !"    #$% 

 

       

Figure 4. Left: Inclusive spectra of positively-charged kaons measured by NA61/SHINE in p+C collisions at 31 GeV (points), compared to predictions by various models (lines). Right: Comparison of the K+/π+ ratio, using the above data along with results shown in Fig. 3, again compared to models [2].

 *  π   #%&'  (%) π 

π   ! "#$ π

             

        

  

   

  #$!#$

          !" #$# !      %&% #     %&'( )*) '%   "#  !" #$#  '"  ()% *  !" ($+!   +*, '%   "#                        π  π   

Figure 5. Left: Rapidity spectra of π− in p+C collisions at 31 GeV and (scaled) central Pb+Pb collisions at 30A GeV, for pT < 1 GeV/c. Corresponding gaussian fits are also shown as lines. Middle: Transverse-mass spectra of π− in the same reactions, for 0

Additional, preliminary results from p+C interactions at 31 GeV have been compared to measure- ments of NA49 in Pb+Pb collisions [2]. Transverse-mass and forward rapidity spectra of negative pions were extracted and compared to NA49 results from central Pb+Pb collisions at 30A GeV (see Λ 0 Δ++ Figure 5). We have also begun investigating yields of , KS and particles in these interactions. Last but not least, we have recently presented the first preliminary π±, K− and p spectra from p+p interactions at 40, 80 and 158 GeV — reference data sets for the phase-diagram scan [2]; an overviews of these results can be found in Figure 6. Analogous measurements from p+p events at 13, 20 and 31 GeV are to follow suit.

00033-p.7 EPJ Web of Conferences

1.6 p+p at 40 GeV/cPreliminary+ 1.6 p+p at 80 GeV/cPreliminary+ 1.6 p+p at 158 GeV/cPreliminary+ 1.4 1.4 1.4 y=0.5 y=0.3 y=0.1 1.2 y=0.7 1.2 y=0.5 1.2 y=0.3 ......

dy [1/(GeV/c)] 1 y=2.7 dy [1/(GeV/c)] 1 y=2.7 dy [1/(GeV/c)] 1 y=2.7 T y=2.9 T y=2.9 T y=2.9 0.8 0.8 0.8 n/dp n/dp n/dp 2 2 2

d 0.6 d 0.6 d 0.6 0.4 0.4 0.4 0.2 0.2 0.2

00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 p [GeV/c] p [GeV/c] p [GeV/c] T T T 1.8 1.8 1.8 p+p at 40 GeV/cPreliminary- p+p at 80 GeV/cPreliminary- p+p at 158 GeV/cPreliminary- 1.6 1.6 1.6

1.4 y=0.1 1.4 y=0.1 1.4 y=0.1 1.2 y=0.3 1.2 y=0.3 1.2 y=0.3 ...... dy [1/(GeV/c)] y=2.1 dy [1/(GeV/c)] y=2.1 dy [1/(GeV/c)] y=2.1 T 1 y=2.3 T 1 y=2.3 T 1 y=2.3 y=2.6 y=2.6 y=2.6 0.8 0.8 0.8 n/dp n/dp n/dp 2 2 2 d 0.6 d 0.6 d 0.6 0.4 0.4 0.4 0.2 0.2 0.2

00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 p [GeV/c] p [GeV/c] p [GeV/c] T T T 0.8 p+p at 40 GeV/cPreliminaryp 0.8 p+p at 80 GeV/cPreliminaryp 0.8 p+p at 158 GeV/cPreliminaryp

0.6 y=1.7 0.6 y=1.9 0.6 y=1.9 y=1.5 y=1.7 y=1.7 ......

dy [1/(GeV/c)] y=0.3 dy [1/(GeV/c)] y=0.3 dy [1/(GeV/c)] y=0.3 T y=0.1 T y=0.1 T y=0.1 0.4 0.4 0.4 n/dp n/dp n/dp 2 2 2 d d d 0.2 0.2 0.2

00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 p [GeV/c] p [GeV/c] p [GeV/c] T T T 0.1 0.1 0.1 - - - p+p at 40 GeV/cPreliminaryK p+p at 80 GeV/cPreliminaryK p+p at 158 GeV/cPreliminaryK 0.08 0.08 0.08 y=0.5 y=0.3 y=0.1 y=0.7 y=0.5 y=0.3

0.06 ... 0.06 ... 0.06 ... dy [1/(GeV/c)] y=1.5 dy [1/(GeV/c)] y=2.3 dy [1/(GeV/c)] y=2.3 T y=1.7 T y=2.5 T y=2.5

n/dp 0.04 n/dp 0.04 n/dp 0.04 2 2 2 d d d

0.02 0.02 0.02

00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 00 0.2 0.4 0.6 0.8 1 1.2 p [GeV/c] p [GeV/c] p [GeV/c] T T T

Figure 6. Double-differential spectra of π+, π−, p and K+ (from top to bottom, respectively) in inelastic p+p inter- actions at 40, 80 and 158 GeV. Only statistical uncertainties are shown; average systematic error is approximately 8 % [2].

4 Summary

NA61/SHINE results have already provided important input for the calculation of the neutrino flux in the T2K beam which was used for the measurement of θ13. Moreover, models employed for simu- lating cosmic-ray showers were significantly constrained. The scan of the phase diagram of strongly interacting matter has begun and first results from p+p collisions were presented.

00033-p.8 ICFP 2012

References

[1] U. Dore, D. Orestano, Rept.Prog.Phys. 71, 106201 (2008), arXiv:0811.1194 [hep-ex] [2] Y. Itow et al. (T2K Collaboration), KEK-REPORT-2001-4, ICRR-REPORT-477-2001-7, TRI- PP-01-05 (2001), arXiv:hep-ex/0106019 [3] E. Ables et al. (MINOS Collaboration), FERMILAB-PROPOSAL-0875, NUMI-L-79 (1995) [4] D. Ayres et al. (NOvA Collaboration), FERMILAB-PROPOSAL-0929 (2004), arXiv:hep-ex/0503053 [5] M.C. Sanchez (LBNE DUSEL Collaboration), AIP Conf.Proc. 1222, 479 (2010) [6] M. Guler et al. (OPERA Collaboration), CERN-SPSC-2000-028, CERN-SPSC-P-318, LNGS- P25-00 (2000) [7] F. Arneodo et al. (ICARUS collaboration), LNGS-P28-2001, ICARUS-TM-2001-03 (2001), arXiv:hep-ex/0103008 [8] A. Badea et al. (KASCADE-Grande Collaboration and the LOPES Collaboration), Nucl.Phys.Proc.Suppl. 136, 384 (2004), arXiv:astro-ph/0409319 [9] J. Abraham et al. (Pierre Auger Collaboration), Nucl.Instrum.Meth. A523, 50 (2004) [10] K. Martens (Telescope Array Collaboration), Nucl.Phys.Proc.Suppl. 165, 33 (2007) [11] C. Alt et al. (NA49 Collaboration), Phys.Rev. C77, 024903 (2008), arXiv:0710.0118 [nucl-ex] [12] Z. Fodor, S. Katz, JHEP 0404, 050 (2004), arXiv:hep-lat/0402006 [13] M. Gazdzicki, M.I. Gorenstein, Acta Phys.Polon. B30, 2705 (1999), arXiv:hep-ph/9803462 [14] M.A. Stephanov, K. Rajagopal, E.V. Shuryak, Phys.Rev. D60, 114028 (1999), arXiv:hep-ph/9903292 [15] K. Grebieszkow (NA49 Collaboration), Nucl.Phys. A830, 547C (2009), arXiv:0907.4101 [nucl-ex] [16] L. Kumar (STAR Collaboration), J.Phys.G G38, 124145 (2011), arXiv:1106.6071 [nucl-ex] [17] N. Antoniou et al. (NA61/SHINE Collaboration), CERN-SPSC-2006-034, CERN-SPSC-P-330 (2006) [18] S. Afanasev et al. (NA49 Collaboration), Nucl.Instrum.Meth. A430, 210 (1999) [19] F. Becattini, J. Manninen, M. Gazdzicki, Phys.Rev. C73, 044905 (2006), arXiv:hep-ph/0511092 [20] N. Abgrall et al. (NA61/SHINE Collaboration), Phys.Rev. C84, 034604 (2011), arXiv:1102.0983 [hep-ex] [21] N. Abgrall et al. (NA61/SHINE Collaboration), Phys.Rev. C85, 035210 (2012), arXiv:1112.0150 [hep-ex] [22] K. Abe et al. (T2K Collaboration), Phys.Rev.Lett. 107, 041801 (2011), arXiv:1106.2822 [hep-ex] [23] V. Uzhinsky (2011), arXiv:1107.0374 [hep-ph] [24] V. Uzhinsky (2011), arXiv:1109.6768 [hep-ph] [25] K. Grebieszkow (NA61/SHINE Collaboration) (2012), arXiv:1201.0742 [nucl-ex] [26] S. Puławski (NA61/SHINE Collaboration), to be published in the proceedings of Quark Matter 2012

00033-p.9