Study of Very Forward Neutrons with the CMS Zero Degree Calorimeter

: orward neutrons are studied in proton-lead collisions at the CMS experiment at the CERN LHC. They provide information on the centrality and event plane of collisions and provide an opportunity to study nuclear breakup. At the CMS experiment they are detected by the Zero Degree Calorimeters (ZDCs) in the |฀|>8.5 pseudorapidity range. The ZDCs are quartz fiber Cherenkov calorimeters using tungsten as absorber. Test beam data and events with a single spectator neutron are used for the calibration of these detectors. A Fourier-based method is used correct for the effect of multiple pPb collisions. The corrected ZDC energy distribution is used to calculate centrality percentiles and unfold the neutron multiplicity distribution. Very with the CMS Zero Degree Calorimeter. Abstract: Forward neutrons are studied in proton-lead collisions at the CMS experiment at the CERN LHC. They provide information on the centrality and event plane of collisions and provide an opportunity to study nuclear breakup. At the CMS experiment they are detected by the Zero Degree Calorimeters (ZDCs) in the | η | > 8.5 pseudorapidity range. The ZDCs are quartz ﬁber Cherenkov calorimeters using tungsten as absorber. Test beam data and events with a single spectator neutron are used for the calibration of these detectors. A Fourier-based method is used correct for the effect of multiple pPb collisions. The corrected ZDC energy distribution is used to calculate centrality percentiles and unfold the neutron multiplicity distribution.


Introduction
Very forward (|η| > 8.5) neutrons are produced in hadron-nucleus and heavy ion collisions. The main physics processes involved are intranuclear cascades [1], nuclear evaporation [2], and nuclear resonances, like the giant dipole resonance [3]. The information gathered from the observation of these neutrons can be used to tag ultraperipheral collisions, calculate the event plane and estimate centrality in heavy ion and hadron-nucleus collisions. Cascade and evaporation nucleons were observed and studied by a wide range of fix target experiments [4]. Ultraperipheral collisions accompanied by nuclear resonances were studied by the STAR Collaboration at RHIC [5,6]. In the ALICE experiments at LHC cascade and evaporation neutrons are used for the estimation of centrality [7,8]. The energy spectrum of very forward neutrons produced in proton-proton collisions was measured by the LHCf experiment [9,10]. In the CMS experiment these neutrons can be observed by the Zero Degree Calorimeters (ZDC). The most important models, which include forward neutron production are DPMJet [2], URQMD [11], Geant4 [12,13], SMASH [14], JAM [15], and PHSD [16].

The Zero Degree Calorimeter of CMS Experiment
The two ZDCs are located in the neutral particle absorber (TAN), roughly ±140 m away from the CMS interaction point, between the two beampipes. They measure neutral particles at pseudorapidity values |η| > 8.5, as the charged products are removed by dipole magnets located between the central CMS detectors and the ZDCs. They are Cherenkov sampling calorimeters consisted of cladded quartz fibers and tungsten plates. They have three different sections ( Figure 1): the electromagnetic (EM) section, the hadron (HAD) section, and the reaction plane detector (RPD). The EM section is 19 radiation lengths long, which is equal to one interaction length and has five transverse segments. The RPD is a 4 × 4 array of quartz tiles and is used to determine the event plane for flow measurements in heavy ion collisions. The HAD section is 5.6 hadronic interaction length and has four longitudinal segments. The plates in the HAD section are tilted by 45 • to maximize the light yield of the fibers. This paper presents results using the EM and HAD sections of the ZDC detector on the lead going side (z < 0) Universe 2019, 5, 210 2 of 8 measured in pPb collisions at the CMS experiment in 2016. A more detailed description of the ZDC detectors can be found in references [17][18][19][20]. Every channel of the detector is read out in 10 timeslices (denoted with TS0 -TS9), each 25 ns long. The main signal arrives in TS4 as shown in the left panel of Figure 2. The bunch spacing in the 2016 pPb datataking was 100 ns, so additional collisions may happen four timeslice before or after the primary signal (out-of-time pileup). The signal value Q i for a certain i channel is extracted as: where Q i,TSX is the signal value in the Xth timeslice for channel i. The second term is used to subtract contribution from pedestal and the tail of out-of-time pileup events. In order to increase the dynamic range of the detector, the largest signals are let to slightly saturate and their signal values are calculated by scaling up the sum of TS4 and TS5 by a scale-factor R . The scale factor is calculated separately for each channel from the distribution of the values:

TS [25 ns]
which is shown in the right panel of Figure 2. The value of R is calculated from the mean of this distribution. This method ensures that both the resolution remains good for the low-energy signals and the whole energy range of the detector can be used.

Calibration
The relative gains of the different channels are matched (intercalibrated) by using distributions of signal ratios between different channels. These are compared to the same distributions from the 2010 PbPb data collection period, following the test beam calibration of ZDCs. An example of two such distributions is shown in Figure 3. The low energy part of the total ZDC signal distribution is shown in Figure 4. Since the slow neutrons are close to being monoenergetic, the peaks correspond to events with single, double, and triple neutrons. The widths of the peaks are due to the slight variation in the neutron energy and the finite resolution of the detector. The spectrum is fitted with the sum of Gaussian shapes, with the n-neutron peak position µ n and width σ n constrained as: where µ 1 and σ 1 are the position and width of the single neutron peak respectively. The results from the fit are used for the absolute energy calibration of the detector, since the position of the first peak corresponds to 2.56 TeV, the energy of nucleons in the Pb ion in pPb collisions with 8.16 TeV collision energy [22]. The detector resolution is found to be 24%, which is slightly larger than the 15% resolution extrapolated from the test beam measurement [20].

Pileup Correction
The ZDC energy spectrum is shown in Figure 5 for two different data collection periods with different average number of collisions in a single event (pileup). An important feature of these spectra is a shoulder at high energies, which is larger in the data set with a larger pileup. It is caused by pileup collisions. This effect is corrected by the following mathematical procedure, based on the compound Poisson distribution.
Assuming that n is the number of simultaneous pPb collisions in a single event and it follows a Poisson distribution: where only the n > 0 case is considered, therefore 1 − e −µ is included in the denominator to ensure the normalization of the above expression. µ is the average number of collisions in an event, that provide at least one neutron in the ZDC. The ZDC energy deposit is described by a random variable X: where Y i is the random variable describing ZDC energy deposit for an event with single collision. The probability density functions of Y i and X are denoted by g(x) and f (x) respectively. The f (x) function can be expressed as: where the * operation stands for convolution. Then using Fourier transform and the convolution theorem one may write: After substituting the definition of p n into the equation above, the infinite series becomes a geometrical series and can be summed up to: By rearranging this equation, it is possible to express g(x) as: providing a formula to calculate the pileup corrected distribution. The result of this calculation is shown at three different µ values in Figure 6. As it is expected the shoulder at high energies disappeared. The method gives a similar result when the value of µ is varied; there is only a moderate variation in the tail of the distributions.

Using ZDC as a Centrality Estimator in pPb Collisions
The typical quantities used to estimate centrality in heavy ion collisions, such as the multiplicity of charged particles, cannot be used in hadron-nucleus collisions, since they are only loosely correlated with the quantities N coll and N part . A good alternative is to use the very forward energy to estimate centrality, since it is correlated with the number of spectator nucleons. The centrality percentiles, calculated from the corrected energy distribution, are shown in Figure 7. For a proper usage of this zero degree energy as a centrality estimator in a physics measurement, a model of spectator neutron production is needed, which connects the zero degree energy with N coll and N part . The current models are not valid in the LHC collision energies [1,2,4,12], therefore results from ZDC detectors give a useful input for the development of these models.

Unfolding Neutron Multiplicity Distribution
From the fits in Section 3, neutron multiplicity distribution can be calculated for events where only a few neutrons are produced. But the fit cannot be extended to the larger ZDC energy range, as the Gaussian distributions would overlap causing rapid oscillation of the amplitudes. This problem can be circumvented by constraining the amplitude values to change smoothly with the neutron number. This property is achieved by considering the calculation as an unfolding problem. The ZDC response is calculated for N-neutron events and a response matrix R is constructed (left panel of Figure 8), assuming the linearity of the detector. Using the response matrix, the original fit can be performed by minimalizing a χ 2 term. The smoothness of the neutron number distribution is ensured by a linear regularization term (Tikhonov regularization [23]), which requires the first derivative of the distribution to be small. The full χ 2 term to minimalize is: where n and e are vectors, whose elements represent the unknown neutron distribution and the measured ZDC spectrum respectively, V is the covariance matrix of e, D is the first derivative matrix and λ is the regularization strength.
The optimal solution is calculated by taking the derivative of the above equation, with respect to n. The derivative can be rearranged into the form of linear system of equations: Universe 2019, 5, 210 The optimal n is calculated from this expression using the LU decomposition method. The result of this calculation is shown in the right panel of Figure 8. This result serves as a strong constraint and may challenge the models of hadron-nucleus collisions and spectator neutron production [2,[11][12][13][14][15][16].

Conclusions
Very forward neutrons are produced in hadron-nucleus and heavy ion collisions via intranuclear cascades, evaporation and nuclear resonances. In the CMS experiment they are observed by the ZDCs. The individual channels of ZDCs in 2016 were gain matched using data collected in 2010. Peaks are observed in the ZDC energy spectrum, which correspond to events with one, two and three nearly monoenergetic spectator neutrons. This provides an opportunity to calibrate the detector and measure the resolution, which is 24%. The pileup effects are corrected using a Fourier deconvolution method. As the very forward energy is suitable to be a centrality estimator, the corrected ZDC energy distribution is used to calculate centrality percentiles. The neutron multiplicity distribution can be unfolded using Tikhonov regularization and the results serve as a constraint for the theoretical nuclear break-up models.