Rejection of randomly coinciding events in ZnMoO$_4$ scintillating bolometers

Random coincidence of events (particularly from two neutrino double beta decay) could be one of the main sources of background in the search for neutrinoless double beta decay with cryogenic bolometers due to their poor time resolution. Pulse-shape discrimination by using front edge analysis, mean-time and $\chi^2$ methods was applied to discriminate randomly coinciding events in ZnMoO$_4$ cryogenic scintillating bolometers. These events can be effectively rejected at the level of 99% by the analysis of the heat signals with rise-time of about 14 ms and signal-to-noise ratio of 900, and at the level of 92% by the analysis of the light signals with rise-time of about 3 ms and signal-to-noise ratio of 30, under the requirement to detect 95% of single events. These rejection efficiencies are compatible with extremely low background levels in the region of interest of neutrinoless double beta decay of $^{100}$Mo for enriched ZnMoO$_4$ detectors, of the order of $10^{-4}$ counts/(y keV kg). Pulse-shape parameters have been chosen on the basis of the performance of a real massive ZnMoO$_4$ scintillating bolometer. Importance of the signal-to-noise ratio, correct finding of the signal start and choice of an appropriate sampling frequency are discussed.


Introduction
Observation of neutrinoless double beta (0ν2β) decay would imply the violation of lepton number conservation and definitely new physics beyond the Standard Model, establishing the Majorana nature of neutrino [1,2,3,4,5]. Cryogenic scintillating bolometers look the most promising detectors to search for this extremely rare process in a few theoretically favourable nuclei [6,7,8,9,10,11,12]. Zinc molybdate (ZnMoO 4 ) is one of the most promising materials to search for 0ν2β decay thanks to the absence of long-lived radioactive isotopes of constituting elements, the comparatively high percentage of molybdenum and the recently developed technique of growing large high quality radiopure ZnMoO 4 crystal scintillators [10,11,13,14,15,16].
However, a disadvantage of the low temperature bolometers is their poor time resolution, which can lead to a significant background at the energy Q 2β due to random coincidences of signals, especially of the unavoidable twoneutrino 2β decay events [17]. This issue is particularly relevant for the experiments aiming at searching for 0ν2β decay of 100 Mo, because of the short half-life of 100 Mo in comparison to the two neutrino double beta (2ν2β) decay T 1/2 = 7.1 × 10 18 yr [18]. Counting rate of two randomly coincident 2ν2β events in cryogenic Zn 100 MoO 4 detectors is expected to be on the level of 2.9 × 10 −4 counts / (keV×kg×yr) at the Q 2β energy (for 100 cm 3 crystals, under a condition that two events shifted in time in 1 ms can be resolved), meaning that randomly coincident 2ν2β decays can be even a main source of background in a future large scale high radiopurity experiment [17].
This work describes the development of pulse shape discrimination techniques to reject randomly coinciding events in ZnMoO 4 cryogenic scintillating bolometers.

Randomly coinciding events in cryogenic bolometers
The energy distribution of the randomly coinciding (rc) 2ν2β events was built using the approach described in [17], under the assumption that two events are not resolved in the time interval 45 ms. 1 The energy spectra of 10 9 Monte Carlo generated events of 100 Mo 2ν2β decay and of two randomly coinciding 2ν2β events are presented in Fig. 1. The rate for the rc events is calculated as I rc = τ · I 2 0 , where I 0 = ln 2 N/T 1/2 , N is the number of 2β decaying nuclei in Zn 100 MoO 4 crystal (with 100% enrichment by 100 Mo) of typical size of ⊘6 × 4 cm (N = 1.28 × 10 24 ), τ = 45 ms, T 1/2 = 7.1 × 10 18 yr.   [19] (upper panel), and Monte Carlo simulated energy spectra of two randomly coincident 2ν2β events, coincident 2ν2β with external gamma events, and randomly coincident external gamma events (lower panel).
In addition to the 2ν2β decay, there could be other background sources contributing to the background in the region of interest due to random coincidence. We have estimated a possible contribution of the external gamma background using the level of background already achieved in the Cuoricino detector [19]. A simplified model of the Cuoricino background (taken from [19] in the energy interval 300 -2620 keV, while exponentially extrapolated below 300 keV and equated to 0 above 2620 keV) and a Monte Carlo simulated energy spectrum of two randomly coincident background events are presented in Fig. 1.
We have also simulated coincidences of the Cuoricino background with the 2ν2β decay of 100 Mo. One can see that the main contribution to background (assuming radiopure Zn 100 MoO 4 crystal scintillators and a cryostat with a level of radioactive contamination similar to that of the Cuoricino set-up) is expected from the 2ν2β decay of 100 Mo. The total counting rate due to the random coinci-dences of 2ν2β decay events and external gamma events in the region of interest is estimated as ≈ 0.016 counts/(year keV kg) for a detector time resolution of 45 ms.

Generation of randomly coinciding signals
Sets of single and randomly coincident signals were generated by using pulse profiles and noise baselines accumulated with a real 0.3 kg ZnMoO 4 crystal scintillator operated as a cryogenic scintillating bolometer with a Ge light detector [20] in Centre de Sciences Nucléaires et de Sciences de la Matière (Orsay, France). Two measurements have been taken into account, the first one with a sampling rate of 5 kSPS (kilosamples per second) both for the light and heat channels, and the second one with the sampling rate 1.9841 kSPS for the both channels. Ten thousand of base-line samples were selected in all the cases.
The pulse profiles of heat and light signals of the detectors (sum of a few hundred pulse samples produced mainly by cosmic rays with energy of a few MeV) were obtained by fit with the following phenomenological function: where A is the amplitude, τ 1 , τ 2 , τ 3 and τ 4 are the time constants.
To generate randomly coinciding signals in the region of the Q 2β value of 100 Mo, the amplitude of the first pulse A 1 was obtained by sampling the 2ν2β distribution for 100 Mo, while the amplitude of the second pulse was chosen as A 2 = Q 2β −A 1 +∆E, where ∆E is a random component in the energy interval [−5, +5] keV (which is a typical energy resolution of a bolometer).
Ten thousand coinciding signals were randomly generated in the time interval from 0 to 3.3·τ R (∆t = [0, 3.3·τ R ], see Section 4.2), where τ R is the rise-time of the signals (defined here as the time to increase the pulse amplitude from 10% to 90% of its maximum). As it will be demonstrated in the section 4.1, the rejection efficiency of randomly coinciding signals (RE, defined as the part of the pile-up events rejected by pulse-shape discrimination) reaches almost its maximal value when the time interval of consideration exceeds (3 − 4)τ R . Ten thousand of single signals were also generated.
A signal-to-noise ratio (defined as the ratio of the maximum signal amplitude to the standard deviation of the noise baseline) was taken 30 for light signals and 900 for heat signals. These values are typical for ZnMoO 4 scintillating bolometers. Examples of the generated heat and light single pulses are presented in Fig. 2.

Methods of pulse-shape discrimination
We have applied three techniques to discriminate randomly coincident events: mean-time method, χ 2 approach, and front edge analysis. We demanded a 95% efficiency in accepting single signals. For an efficient discrimination, it is very important to develop a method for a good determination of the pulse start position.

Reconstruction of the time origin of the events
The following procedure was used to reconstruct the time origin of each signal: 1. Preliminary search for the presence of a signal by a very simple algorithm, which searches for a channel where the signal amplitude exceeds a certain level (typically about one third of the signal maximum value); 2. Summation of the data over a certain number of channels (typically over 2-6 channels for the light signals, depending on the time structure of signal and noise data; this procedure was not used for the heat signals); 3. Calculation of the standard deviation of the baseline fluctuations; 4. Search for the pulse start under the request that the signal exceeds a certain number of standard deviations of the baseline and -in a case of heat signals -the amplitude in the next several channels increases channel by channel.
The algorithm was optimized for each data set, taking into account the sampling rate, the time properties of the signals and noise, and the signal-to-noise ratio.

Mean-time method
The following formula was applied to calculate the parameter t (mean-time) for each pulse f (t k ): where the sum is over time channels k, starting from the origin of a pulse and up to a certain time.
As a first step we have chosen the time interval ∆t to analyze efficiency of the pulse-shape discrimination. Six sets of single and randomly coinciding light (with τ R = 3 ms) and heat signals (with τ R = 13.6 ms) were generated in the time intervals ranging from 0 to a maximum value, varying from 1 to about 6 pulse rise-times. The results of this analysis are presented in Fig. 3. The uncertainties of the rejection efficiency were estimated by analysis of three sets of data generated using three sets of different noise baseline profiles (about 3300 profiles in the each set). One can see that the rejection efficiency of randomly coinciding signals reaches its maximal value when the time interval ∆t is larger than (3 − 4)τ R . All the further analysis was done by using data generated in the time interval ∆t = [0, 3.3τ R ]. A typical distributions of the mean time parameters for single and pile-up events are presented in Fig. 4. The rejection efficiency of randomly coinciding pulses, under the requirement to detect 95% of single events, is 92.2%.
One could expect that the rejection efficiency of pulseshape discrimination depends on the choice of the time interval used to calculate a discrimination parameter. For instance, in Fig. 5 the results of the mean-time method optimization are presented. The rejection efficiency has a maximum when the mean-time parameter is calculated from the signal origin to the 30th channel which approximately corresponds to ∼ τ D . All the discrimination methods were optimized in a similar way.
The dependences presented in Fig. 5 demonstrate also the importance to optimize an algorithm to find the start of the signals (see Section 4.1), particularly for the light pulses with comparatively low signal-to-noise ratio. Rejection efficiency is substantially higher when the start position of each pulse is known from the generation algorithm. The distributions of start positions for the sets of single and randomly coinciding events shown in Inset of Fig. 5 demonstrate that it is more problematic to correctly find  the start position of randomly coinciding signals than that of single events.

χ 2 method
The approach is based on the calculation of the χ 2 parameter defined as where the sum is over time channels k, starting from the origin of pulse and up to a certain time, and f S (t) is defined by Eq. (1). The number of channels to calculate the χ 2 has been optimized to reach a maximal rejection efficiency.

Front edge analysis
The front edge parameter can be defined as the time between two points on the pulse front edge with amplitudes Y 1 % and Y 2 % of the pulse amplitude. The parameters Y 1 and Y 2 should be optimized to provide maximal rejection efficiency. For instance, the highest rejection efficiency for heat pulses with τ R = 13.6 ms was reached with the front edge parameter determined as time between the signal origin and the time where the signal amplitude is Y 2 = 90% of its maximum (RE = 98.4%). However, the rejection efficiency of the front edge method is limited due to the fraction of randomly coinciding events with a small first (with the amplitude A 1 below Y 1 ) or second pulse (with a low amplitude, and appearing well after the first signal maximum).

Results and discussion
The methods of pulse-shape discrimination are compared in Table 1. The data were obtained with start positions of the signals found by our algorithms, and using a priori known signal start positions from the generation procedure (to estimate a maximum achievable efficiency). All the methods give a 86%−92% rejection efficiency by using the light signals with a rise-time of 3 ms and 98%−99% for the much slower heat signals with a rise-time of 13.6 ms. One can conclude that the signal-to-noise ratio (set to 30 for the light and to 900 for the heat signals at the energy Q 2β , as observed in real bolometers) plays a crucial role in the pulse-shape discrimination of randomly coinciding events in cryogenic bolometers. Analysis of signals with lower level of noise allows to reach much higher rejection efficiency even with slower heat signals. Dependence of the rejection efficiency (by using the mean-time method) on the signal-to-noise ratio for heat signals confirms the assumption (see Fig. 6).
We tried also to analyse a dependence of the rejection efficiency on the rise-time of light pulses (see Fig. 7). We have assumed that the pulse amplitude does not change by shortening the rise-time. One could expect that for Table 1. Rejection efficiency of randomly coinciding 2ν2β events by pulse-shape discrimination of light and heat signals for the two conditions of the signal start determination, i.e. (i) start of the signals known from the generation procedure, (ii) and start position found by the pulse profile analysis.

Channel, Start
Mean-time Front edge χ 2 method, rise-time position method, % analysis, % % Signal-to-noise ratio Fig. 6. Dependence of the rejection efficiency (using the meantime method) on the signal-to-noise ratio for the heat channel.
faster signal any methods should give a higher efficiency of pulse-shape discrimination. However, the trend of the rejection efficiency improvement for faster signals is rather weak. Furthermore, the rejection efficiency for pulses with a rise-time of 2 ms is even worse in comparison to slower signals with rise-times of 3 ms and 4.5 ms. This feature can be explained by a rather low sampling rate (1.9841 kSPS) used for the data acquisition. Indeed, the pulse profiles acquired with this sampling rate are too discrete: for instance, the front edge of the signals with the rise-time 2 ms is represented by only 4 points. Such a low discretization even provides difficulties to set the acceptance factor of single events at a certain level (95% in our case), particularly in the front edge analysis. The importance of the data acquisition sampling rate for the pulse-shape analysis was proved by taking the noise baselines data with 2 times lower sampling rate (we have transformed the data by simply averaging two adjacent channels to one). The rejection efficiency decreases in this case, as one can see in Fig. 7.
Following the previous discussion, we will calculate now how our pile-up rejection procedure would improve the background figure estimated at the end of Section 2, and corresponding to ≈ 0.016 in an enriched Zn 100 MoO 4 detector based on a crystal size of ⊘6 × 4 cm. We remind that this value was obtained assuming that pulses separated by an interval longer than 45 ms would be far enough to be analysed independently and on the contrary all the concomitant pulses within this time interval would give rise to a single amplitude equal to the sum of the individual ones. We are authorised to abate the resulting background level by 99.3% in the heat channel and 92.3% in the light channel, as these are our best pile-up rejection efficiencies with unknown pulse start position (see Table  1). 2 The final estimated background level after pile-up rejection is ≈ 1.1 × 10 −4 counts/(year keV kg) using pulseshape discrimination in the heat channel. A higher value is obtained using the light channel because of the worse rejection efficiency. An improvement of the speed and of the signal-to-noise ratio in cryogenic bolometers -both in the heat and light channel -is an important experimental goal to further enhance the rejection of RC-generated background.

Conclusions
Random coincidence of events (especially but not only from 2ν2β decay) could be one of the main sources of background in cryogenic bolometers to search for 0ν2β decay because of their poor time resolution, particularly for 100 Mo due to the short half-life in comparison to that of the two-neutrino mode. However, this background can be effectively suppressed with the help of pulse-shape discrimination. The randomly coinciding 2ν2β decay signals were discriminated with an efficiency at the level of 99% by applying the mean-time approach to the heat signals from ZnMoO 4 cryogenic bolometer with a rise-time of about 14 ms and a signal-to-noise ratio of 900, and at the level of 97% for the light signals with 3 ms rise-time and signalto-noise ratio of 30 (however, the last estimation was obtained for the signals with a priory known start position). χ 2 approach provides comparable rejection efficiencies, while the results of the front edge analysis are slightly worse due to problems in discriminating pile-up events when one of the randomly coinciding signals is too small to be detected by this method.
The signal-to-noise ratio looks the most important feature to reject randomly coinciding events, particularly in ZnMoO 4 due to the comparatively low light yield, which leads to a rather low signal-to-noise ratio in the light channel.
Development of algorithms to find the origin of a signal with as high as possible accuracy is requested to improve the rejection capability of any pulse-shape discrimination technique. The sampling rate of the data acquisition should be high enough to provide effective pulseshape discrimination of randomly coinciding events. Finally, any pulse-shape discrimination methods should be optimized taking into account certain detector performance to reduce the background effectively.
The analysis performed in this work proves further that the counting rate due to the random coincidence of events can be reduced to a level ≈ 10 −4 counts/(year keV kg), which makes ZnMoO 4 cryogenic scintillating bolometers very promising to search for neutrinoless double beta decay at a level of sensitivity high enough to probe the inverted hierarchy region of the neutrino mass pattern [10].

Acknowledgements
This work was supported in part by the project "Cryogenic detector to search for neutrinoless double beta decay of molybdenum" in the framework of the Programme "Dnipro" based on the Ukraine-France Agreement on Cultural, Scientific and Technological Cooperation. The study of ZnMoO 4 scintillating bolometers is part of the program of ISOTTA, a project receiving funds from the ASPERA 2nd Common Call dedicated to R&D activities.