Study on the accidental background of the JSNS 2 experiment

JSNS 2 (J-PARC Sterile Neutrino Search at J-PARC Spallation Neutron Source) is an experiment which searches for sterile neutrinos via the observation of ¯ ν µ → ¯ ν e appearance oscillations using muon decay-at-rest neutrinos. The data taking of JSNS 2 have been performed from 2021. In this manuscript, a study of the accidental background is presented. The rate of the accidental background is (9.29 ± 0 . 39) × 10 − 8 / spill with 0.75 MW beam power and comparable to the number of searching signals.


Introduction
The existence of sterile neutrinos has been a crucial issue in the neutrino physics community for over 20 years.The experimental results from [1][2][3][4][5][6] could be interpreted as indications of the existence of sterile neutrinos with a mass-square differences of around 1 eV 2 .
The JSNS 2 experiment, proposed in 2013 [7], is designed to search for neutrino oscillations due to a sterile neutrino at the Material and Life science experimental Facility (MLF) in J-PARC.MLF provides an intense and high-quality neutrino flux of 1.8×10 14 ν/year/cm 2 , from muon decay-at-rest (µDAR) produced using a 1 MW proton beam with a 25 Hz repetition rate [8].The neutrinos are produced by injecting 3 GeV protons from a arXiv:2308.02722v2[hep-ex] 22 Apr 2024 rapid cycling synchrotron onto a mercury target in the MLF.The experiment uses a Gadolinium (Gd) loaded liquid scintillator (Gd-LS) detector with 0.1 w% Gd concentration placed at 24 m from the target.
The JSNS 2 experiment aims to perform a direct test of the LSND observation [1].The same experimental principle used by the LSND experiment [1] to try and observe the νµ → νe oscillation is used: inverse beta decay (IBD).There are several improvements offered by the JSNS 2 experiment.In order to identify IBD events, a delayed coincidence between the positron signal (prompt signal: up to 53 MeV) and neutron capture signal is used for selection.Gd is used to identify neutron captures.After capturing thermal neutrons on Gd, neutron capture on Gd generates gamma-rays with higher energies and shorter capture times (∼8 MeV, ∼30 µs) than neutron capture on hydrogen (2.2 MeV, ∼200 µs).Therefore, accidental backgrounds coincident in the delayed signal region can be reduced by ∼ 6 times compared to the hydrogen capture used in the LSND experiment, due to the shorter capture time.In addition, the short-pulsed beam, two 100 ns pulses in a 600 ns interval in each spill with a repetition of 25 Hz, enables us to set a timing window for the IBD prompt signal to 2.0 to 10 µs from the proton beam collision so that the neutrinos from pion and kaon decay and fast neutrons generated at the target can be rejected efficiently.However, the efficiency for the µDAR neutrinos can be kept at 48% because of the muon lifetime (2.2 µs).The cosmogenic background is also reduced by a factor of 10 −4 .
For JSNS 2 , understanding the accidental background is essential since it is one of main backgrounds for the sterile neutrino search.Another important background, the correlated background, is described in [9].A detailed discussion about the signal detection principle and the background rejection technique can be found in [7,8].

Setup
The JSNS 2 experiment has been taking data with a single detector since 2021.Except for the beam maintenance period, which typically corresponds to a few months over the summer, the data has been accumulated continuously.The proton beam power has been increased from 600 kW in 2021 to 840 kW in 2023.There is usually a one-day facility maintenance per week and we continue to take data during that time to acquire beam-off data.The integrated number of proton-ontarget (POT) collected was 2.94 × 10 22 , corresponding to less than ∼28.0% of the required POT of the JSNS 2 experiment.
To understand the rate and properties of the accidental background, dedicated calibration data were taken.The accelerator-driven timing signal was used for the trigger [10,11], which caused a 125 µs time window for data acquisition system (FADC) to be opened, with no energy bias.The total available number of spills is 2.46 ×10 6 , therefore the expected number of IBD signal events from this calibration run is estimated to be much less than 1 [8].Thus, all observed events are likely to be background.

Experimental setup
The JSNS 2 detector is a cylindrical liquid scintillator detector with 4.6 m diameter and 3.5 m height located at a distance of 24 m from the mercury target of the MLF.It consists of 17 tonnes of Gd-LS contained in an acrylic vessel, and 33 tonnes unloaded liquid scintillator (LS) in a layer between the acrylic vessel and a stainless steel tank.The LS and the Gd-LS consist of LAB (linear alkyl benzene) as the base solvent, 3 g/L PPO (2,5-diphenyloxazole) as the fluor, and 15 mg/L bis-MSB (1,4-bis(2-methylstyryl) benzene) as the wavelength shifter.The LS volume is separated into two independent volumes by an optical separator.The region inside the optical separator, called the "inner detector", consists of the entire volume of the Gd-LS and ∼25 cm thick LS layer.Scintillation light from the inner detector is observed by 96 Hamamatsu R7081 photomultiplier tubes (PMTs) each with a 10-inch diameter.The outer layer, called the "veto layer", is used to detect cosmic-ray induced particles coming into the detector.A total of 24 of 10-inch PMTs are set in the veto layer.On the whole inner surfaces of the veto layer, reflection sheets are attached in order to improve the collection efficiency of the scintillation light.

Data acquisition and trigger system
PMT signal waveforms from both the inner detector and the veto layer are digitized and recorded at a 500 MHz sampling rate by 8-bit flash analog-to-digital converters (FADCs).As a trigger, we utilize a 25 Hz periodic signal from the accelerator scheduled timing which directs the proton beam towards the MLF target, called the "kicker trigger".The width of the acquired waveform in this trigger scheme is set to 125 µs , which fully covers the prompt and delayed signal timing window of the IBD events.The rate limitation of the JSNS 2 data acquisition system under this condition is 5 Hz, thus, a pre-scaling factor of 5 is required.We mainly used this trigger and data acquisition to obtain the beam-on data for an accidental background estimation within the sterile neutrino search.
Detailed description of the detector and the triggers are given in [10,11], respectively.

Event selection
The data set used for the background estimation was obtained using the kicker trigger.The total number of beam spills was 2460399, equivalent to ∼6 days of data taking with a pre-scale factor of five.Given the readout 125 µs time, the obtained waveforms contain multiple events.We therefore used an event definition based on the number of hit PMTs in order to extract each event.We constructed a hit time and charge series at each trigger by accumulating hit information along the FADC window with 60 ns coincidence width over all the PMTs.The coincidence width was determined by considering a typical PMT pulse shape and a safety factor from timing calibration.The event discrimination threshold is set to 10 hits and 50 p.e., which corresponds to an energy well below 1 MeV.
The event vertex position and energy reconstruction is performed simultaneously based on a maximumlikelihood algorithm using the charge response of each PMT.Both vertex position and energy were calibrated by deploying a 252 Cf source and using the 8 MeV peak in the energy spectrum resulting from thermal neutron capture on Gd (nGd).The reconstruction performance can be found in [12] for nGd events and in [13] for events with up to 60 MeV using Michel electrons.
The selection criteria and estimated efficiency are given in Table 1 and a more detailed description is found in [8].Each efficiency in Table 1 is for each corresponding criterium only.
The single rates of the prompt and the delayed candidates are separately estimated using the energy and the timing selections shown in Table 1.Note that a time difference from the beam collision timing, 5 ≤ ∆t beam−d ≤ 105 µs is used for the delayed single rate calculation instead of the selection in Table 1.These single rates are multiplied to predict the number of the accidental backgrounds.The fiducial volume is defined with R < 140 cm and |z| < 100 cm region to avoid external backgrounds.Note that origin of the coordinate system is the center of the detector, and R is defined as R = x 2 + y 2 .
The prompt IBD candidates are selected using a time difference from beam collision timing (∆t beam−p ) and its energy (E p ).We applied the following requirements; 2 ≤ ∆t beam−p ≤ 10 µs and 20 ≤ E p ≤ 60 MeV, in order to fully cover µDAR neutrinos from the mercury target.The timing selection rejects beam-induced Table 1 The IBD selection criteria and their efficiencies in the JSNS 2 experiment [8].The single rates of the prompt and the delayed candidates are estimated separately using the energy and the timing selection in this table.Note that the timing window: 5 ≤ ∆t beam−d ≤ 105 µs is used for the delayed single rate calculation.Each efficiency is for to each corresponding criterium only.

Requirement
Efficiency /% fast neutrons in the beam on-bunch timing (0 ≤ ∆t beam−p ≤ 1.5 µs) as well as neutrino backgrounds from kaon and pion decay whose lifetimes are 12 ns and 26 ns, respectively.
In order to identify neutron captures on Gd, the 8 MeV peak of the delayed signal energy (E d ) is selected using a requirement of 7 ≤ E d ≤ 12 MeV [14].There are nGd events associated with fast neutrons induced by the beam contributing to the IBD delayed candidates as an accidental background [15].Pulse Shape Discrimination (PSD) is used within the IBD prompt timing region to reduce neutrons induced by cosmic rays as shown in [9].However, it is optimized for the energy range of 20 ≤ E p ≤ 60 MeV.Thus, we need another strategy to reject those nGd events.Since these nGd events correlate spatially with an activity made by beam-induced fast neutrons, we can reject them with the spatial correlation.In particular, we applied a requirement on the spatial difference between on-bunch events and delayed candidates, ∆VTX OB−d ≥ 110 cm.The distribution of ∆VTX OB−d and the efficiency estimation can be found in [14].The on-bunch event tagging condition is set to 0 ≤ ∆t beam−OB < 1.5 µs and 1 ≤ E OB ≤ 200 MeV.
Figure 1 demonstrates the energy and timing selection windows in a two-dimensional distribution of energy and timing.The red and green boxes represent the prompt and delayed signal regions.The beam on-bunch event are also defined as an orange dashed box.
The events induced from muons or incoming particles from outside are rejected from both prompt and delayed candidates by the veto region of the detector.The events which have more than 30 photo-electrons Fig. 1 A two-dimensional distribution of energy and timing before event selection, used to demonstrate the IBD selection region.The selected regions for the IBD prompt candidate from positron (red solid box), the IBD delayed candidate from gamma-rays resulted in thermal neutron capture on Gd (green dashed box) and the beam on-bunch event (orange dashed box) are overlaid.Note that the events are shown around the prompt signal timing region.There are two event clusters within 0 to 1.5 µs , which reflect the proton beam structure of the MLF.They are caused by neutrons produced at the mercury target.One can see that the IBD prompt signal region is well separated from the on-bunch region.The events in the IBD delayed signal region must also satisfy ∆t p−d > 0 in the actual delayed coincidence.
(p.e.) of a total charge for the top-side 12 PMTs or 40 p.e. for the bottom-side 12 PMTs are rejected.
The decayed Michel electrons (ME) from muon decay are also crucial background.The IBD candidates are rejected if parent muons are found during 10 µs before the candidate.The events that have more than 100 p.e. for the top 12 or bottom 12 PMTs are categorized as parent muons except for the beam timing.At the beam timing, incoming neutrons induced by the beam (0.7464 ± 0.0006 / spill) could pass the criteria for selecting parent muons; thus we applied a different condition only during the beam timing.The equation E + Q veto /9 > 200 is used to define the parent muons of Michel electrons during the beam timing, where E is the deposited energy in the target and catcher regions and Q veto is the total charge of veto region with p.e. unit."Q veto /9" converts p.e. to MeV units in the equation.For example, Q veto = 9 p.e. corresponds to 1 MeV in the veto region.The number of selected events of the IBD prompt and the delayed candidates are 542 and 44336, respectively, thus the single rates for those are (2.20±0.09)×10 −4 / spill (prompt) and (1.80±0.01)×10 −2 / spill (delayed), respectively.Note that the uncertainties men-  tioned in this paper contain both the statistical and systematic uncertainties.For the systematic uncertainties, the energy scale, the FADC timing and the reproducibility are considered.Figure 2 (d) shows that the rejection using spatial correlation between beam neutrons and delayed candidates (∆VTX OB−d ) works well.

Beam on-off comparison of single rates
The comparison of single rates between beam-on and beam-off are also performed.The beam-off data is also taken using the the kicker trigger with 125 µs, and the IBD event selection is identical to that in beam-on data.The total number of the beam-off triggers is 466,348, which is equivalent to about 1 day of data taking.Table 2 shows the comparison of the IBD prompt and delayed candidates.This result shows the rate of the IBD prompt background are independent from the beam power, which indicates that the cosmic ray induced particles are dominated in this region.On the other hand, that of the delayed background has large dependence of the beam.

Accidental background rates
The accidental background rate is evaluated by multiplying the single rates as follows: where R acc is the accidental background rate (/spill), R p is the single rate of the IBD prompt candidates (/spill), and R d is the rate of of delayed candidates (/spill).If the additional selections are applied, the selection efficiency to the accidental background (ϵ cut ) should be also multiplied.As shown in Table 1, the spatial correlation selection between the prompt and delayed candidates, ∆VTX p−d , and timing likelihood should also be considered.In this manuscript, only the spatial correlation cut is described.The likelihood will be considered in a different future paper.In order to evaluate the efficiency of ∆VTX p−d of the pure accidental uncorrelated background, a novel technique is invented, which we call "spill shift method".Within the same beam spill, the correlated background is dominating, as discussed in [9].However, if we use subsequent spills, beam correlated events disappear after more than 40 ms and thus a pure sample of uncorrelated accidental events is obtained.Once the IBD prompt candidate in one spill is found, the different beam spills are analyzed to find the paired delayed candidates.At first, the next beam spill of the spill that have the IBD prompt candidate is used.Secondly, the beam spill that is different by two spills from the spill with IBD prompt candidate is used.This process is repeated from the 1 st to the 10,000 th next spills to get higher statistics.Figure 3 shows the cartoon illustrating the principle of the spill shift in the case of 1 spill shift.
Figure 4 shows the estimated ∆VTX p−d distribution.The estimated efficiency of the ∆VTX p−d ≤ 60 cm cut is 5.1±0.1%.

Accidental background rates
The expected accidental background rate is (9.29±0.42)× 10 −8 / spill / 0.75MW using the Equation 1.The numbers used in this calculation of the accidental back- ground rate are summarized in Table 3.The single rates, the spatial correlation cut efficiency mentioned above and the assumed efficiency of timing likelihood (∼46%) described in [8] are used.
Fig. 5 The expected number of accidental background events from the spill shift method during this calibration run with 2,460,399 spills after the ∆VTX p−d ≤ 60 cm cut.This is for the complementary method to estimate the accidental background rate.
The spill shift method can also be utilized to estimate the rate of the accidental backgrounds independently from the single rates method because each nth spill shift case can find number of accidental background events in this control sample, individually.Figure 5 shows the number of accidentally paired background events in each spill shift cases (1,2, ... 10,000 spills shift) during 2,460,399 spills after the ∆VTX p−d cut, i.e.; one certain n-th spill shift case provides one entry in this histogram.The mean of the Poisson dis-tribution is 0.494 ± 0.008 events over 2,460,399 spills, thus the estimated accidental background rate is (9.24± 0.15) × 10 −8 / spill / 0.75MW.This is consistent with the single rates method.
The expected oscillation IBD signal rate is 4.59×10 −8 / spill / MW [8] with the LSND best fit oscillation parameters.Thus, JSNS 2 has a comparable accidental background rate as the expected neutrino oscillation signal with the LSND best bit oscillation parameters.

Summary
JSNS 2 aims to perform a direct test of the positive result of the LSND experiment using a decay-at-rest neutrino source at the MLF and a Gd-LS detector.We started data taking in 2021.
The calibration runs with the accelerator scheduled timing to study the accidental background have been performed.As a result, the single rates of the IBD prompt and delayed candidates at 0.75 MW of averaged beam power are (2.20±0.09)×10 −4 / spill and (1.80±0.01)×10 −2 / spill, respectively.The expected rate of the accidental background is (9.29±0.39)× 10 −8 / spill, which is a similar level with that of the IBD neutrino oscillation signal with the LSND best fit oscillation parameters.We anticipate the future improvements to reduce the accidental background with sophisticated likelihood and/or machine-learning techniques.

Fig. 2
shows the distributions of E p (a), E d (b) and T p from beam (c), T d from beam (d) after applying the IBD selection.

Fig. 2
Fig. 2 The distributions of the variables used in the IBD selection: (a) E p , (b) E d , (c) T p from beam timing and (d) T d from beam.Each distribution has all of the selection criteria applied, except for the selection on the variable displayed.The red boxes show the selection criteria for E p and E d .

Fig. 3
Fig.3The principle of the spill shift.

Fig. 4
Fig. 4 The estimated ∆VTX p−d distribution of the accidental background with the spill shift method.The red line shows a cut value, and the remained fraction of the accidental background is 5.1±0.1%.

Table 3
The numbers used in Equation 1 and the calculated accidental background.