Enhanced light signal for the suppression of pile-up events in Mo-based bolometers for the 0νββ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu \beta \beta $$\end{document} decay search.

Random coincidences of events could be one of the main sources of background in the search for neutrino-less double-beta decay of 100\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{100}$$\end{document}Mo with macro-bolometers, due to their modest time resolution. Scintillating bolometers as those based on Li2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}MoO4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_4$$\end{document} crystals and employed in the CROSS and CUPID experiments can eventually exploit the coincident fast signal detected in a light detector to reduce this background. However, the scintillation provides a modest signal-to-noise ratio, making difficult a pile-up pulse-shape recognition and rejection at timescales shorter than a few ms. Neganov–Trofimov–Luke assisted light detectors (NTL-LDs) offer the possibility to effectively increase the signal-to-noise ratio, preserving a fast time-response, and enhance the capability of pile-up rejection via pulse shape analysis. In this article we present: (a) an experimental work performed with a Li2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}MoO4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_4$$\end{document} scintillating bolometer, studied in the framework of the CROSS experiment, and utilizing a NTL-LD; (b) a simulation method to reproduce, synthetically, randomly coincident two-neutrino double-beta decay events; (c) a new analysis method based on a pulse-shape discrimination algorithm capable of providing high pile-up rejection efficiencies. We finally show how the NTL-LDs offer a balanced solution between performance and complexity to reach background index ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-4}$$\end{document} counts/keV/kg/year with 280 g Li2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}MoO4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_4$$\end{document} (100\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{100}$$\end{document}Mo enriched) bolometers at 3034 keV, the Qββ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{\beta \beta }$$\end{document} of the double-beta decay, and target the goal of a next generation experiment like CUPID.

Both CROSS [21] and CUPID [25] experiments develop a technology based on Li 2 MoO 4 macro-bolometers coupled to bolometric light detectors and use neutrontransmutation-doped Ge thermistors (NTD) [26] as sensors.These sensors can be produced in large numbers and they can be read out via simple, conventional, lownoise JFET-based electronics [28].Their major drawback is an intrinsic slow (1-10 ms) response time, es-sentially due to 1) their high impedance at the optimal working point, 2) the gluing interface with the crystal, and 3) the internal electron-phonon decoupling.Slow response time can lead to pile-up events mimicking single, normal events at the Q-value of the double-beta decay Q ββ [39,18].The identification of these events is extremely important when searching for 0ν2β decay of 100 Mo: in fact, having the 100 Mo a relatively "short" two-neutrino double-beta (2ν2β) half-life (T 2ν2β  1/2 = 7.1• 10 18 years [22]), this latter can bring a large background of pile-up events in the ROI.
Other experiments as AMoRE [27], also looking for 100 Mo neutrinoless double-beta decay, adopted metallicmagnetic-calorimeter (MMC) as temperature sensors [29].The major advantage of this sensor technology is to provide signals with faster response times than that of the NTDs, reducing the pile-up background in the ROI.
The counting rate of randomly coincident 2ν2β decay events in a 280 g Li 2 100 MoO 4 CUPID detector is estimated on the level of 3.3 • 10 −4 counts/keV/kg/year (ckky) at 3034 keV, assuming a time-resolving capability of 1 ms [18]: randomly-coincident 2ν2β decay events represent at the moment the main source of background for the next-generation, large-scale, high-radio-pure experiments based on the CUPID technologies [25,44].
Pile-up rejections using pulse-shape analysis of the heat channel have been presented [23]: rejection performances of 90% for pile-up times down to 2 ms were reported.However, these values are not sufficient to reach the ultimate CUPID goal.
The approach followed here is the same as that reported in Ref. [18], consisting of exploiting the light signal to improve pile-up rejection through pulse-shape discrimination.In fact, as reminded above, both CROSS and CUPID will be provided with optical bolometers to detect the scintillation light in coincidence with the heat signal measured in the Li 2 100 MoO 4 crystal.These light detectors consist of thin Ge wafers equipped with NTD's as the Li 2 100 MoO 4 crystal.The collected light is detected as a temperature pulse in the wafer.The primary function of the double heat-light readout is the rejection of the alpha background by exploiting the lower scintillation yield of alphas with respect to betas/gammas for the same deposited energy.Besides this, the light detector can play a crucial role also to mitigate the pile-up background, since its signal risetime is about ten time shorter than that of the heat signal.In order for this method to be effective, the signal-to-noise ratio needs to be enhanced, which can be obtained through the so-called Neganov-Trofimov-Luke (NTL) effect [45,46].With respect to past work, which was mainly conceptual, our results here are supported by extensive experimental results and more sophisticated simulations, providing a convincing and operational method to reject the background related to 2ν2β decay events down to the desired level.
In this study we report about an experimental and simulation work which provides a technological solution combined with simulations and analysis methods to demonstrate how to reach a background index (BI) lower than ∼10 −4 ckky in the region of interest (ROI) of 100 Mo 0ν2β decay.
It consists of three main sections: -an experiment performed with a 245 g Li 2 MoO 4 scintillating bolometer equipped with a NTL effect boosted light detector (NTL-LD) as such in [24].
The bolometer was realized in the framework of the CLYMENE project and the experiment run at the C2U (CROSS Cryogenic Underground) facility at Canfranc Underground Laboratory (Spain) as part of the CROSS/CLYMENE joint R&D program (Section 2); -a new analysis method based on an optimal-filtering pulse-shape discrimination algorithm, capable to provide enhanced pile-up rejection efficiencies (Section 3); -a simulation method to reconstruct synthetically produced random coincidences of 2ν2β decays: it is fully based on signals and noise obtained from detector measurements (Section 4).
2 The experiment

Setup description
Figure 1 shows a picture of the scintillating macrobolometer mounting used in this work.The experiment was performed in a low radioactivity pulse tube dilution refrigerator (Hexadry-400, from Cryoconcept) at around 15 mK.The main bolometer consists of a 245 g Li 2 MoO 4 cylindrical single crystal; it is fit in a silvercoated reflective cavity and equipped with two NTD-Ge sensors (3×3×1 mm 3 and 10×3×1 mm 3 ) and one 300 kΩ Si:P heater [30].This detector was developed within the CLYMENE project [31,32] (ANR funded) and studied in synergy with the CROSS ERC granted project [21].A Neganov-Trofimov-Luke-effect boosted germanium light detector, identical to the one described in [24], is coupled to the cavity and records the scintillation emitted in coincidence with the particle events occurring in the main bolometer, allowing to reject the α background.This detector is fabricated with a 170 µm thick, = 44 mm high purity Ge wafer on top of which annular aluminum electrodes of 200 µm width and 3.8 mm pitch are evaporated.The electrodes are interconnected to produce two different electrode sets, and kept at a voltage difference V e .The NTL-LD was operated with several V e up to a maximum of 50 V: for this value it returned an enhanced sensitivity about 10 times larger than at V e = 0 V and no additional contribution in the Noise Power Spectrum (NPS).
An optical fiber is coupled with a room temperatures 880 nm wavelength LED and used to shine the NTL-LD in the cavity.The LED was driven via Keysight 33500B wave-function generator sending typically square voltage pulses of about 400 ns time-width and a voltage amplitude around 1 V.

Measurements
We artificially produced pile-up events in the light detector via LED photo-pulses sent through the optical fiber.The voltage amplitudes, A1 and A2, of the wavefunction generator were tuned, together with the pulse duration (typically a few µs, much shorter than the risetime τ r of the light detector of the order of 1 ms), to provide single pulses of energy E = 450 eV.Hence, a perfect synchronous two LED pulses (∆t = 0 ms) delivers 900 eV, i.e. the amount of scintillating energy recovered on a single light detector for an event at the 100 Mo Q ββ (3 MeV), for a light yield of 0.3 keV/MeV as in CROSS or CUPID setup [19,33].For simplicity, we worked with A1 = A2 as the simulation presented hereafter shows that this pulse amplitude combination gives the main pile-up contributions, and sending two subsequent pulses spaced by given ∆t time interval.We adopted the following ∆t time pattern: 0.1; 0.3; 0.5; 0.7; 1.0; 3.0; 5.0; 0.05 ms.We delivered 100 (pile-up) signals for each ∆t, at a repetition rate of 0.2 Hz.We have performed several set of measurements, by varying the NTL-LD electrode voltage bias V e (10, 30 and 50 V).As a precaution, between each set of voltage we used the LED to charge-reset the NTL-LD, and clean the residual, spurious electric field built-up in the semiconductor due to charge trapping as described in [34] 1 .
We measured an average performance on the light channel of σ ∼ 90 eV (∼ 9 eV) baseline noise when operated at V e = 0 V (50 V) and sensitivity of 1 µV/keV (10 µV/keV).We explored different NTL-LD detector working points by scanning with respect to different sensor bias currents (Table 1).We finally chose to work at a bias of 1 nA, which provides a signal rise-time of 1.2 ms.
The NTL-LD signals were filtered with an analog low-pass Bessel filter: the cut-off frequency was set at 500 Hz [35,36,37], above the intrinsic detector cut-off frequency F intr = 0.35/τ r .All the measurements where 24 bit sampled at a frequency of 5 kHz and continuously recorded (streaming) during the course of the experiment.The streaming files were analyzed offline via an analysis software using Gatti-Manfredi matched opti-Table 1 Characterization of different NTL-LD working points with respect to several I N T D thermal sensor bias current, for a mixing chamber temperature regulated at 15 mK.We finally performed all the measurements at 1 nA bias current, which is chosen to obtain a good compromise between fast response and decent signal amplitude.mum filter [38]: it returns the energy estimate and pileup rejection parameters, as defined hereafter in Section 3.

Data analysis: method and results
To optimally reconstruct the energy of an event over the noise by exploiting the maximum of a signal pulse we employed the optimal filtering algorithm reported in [38].Few elements are required: 1) the typical average baseline power spectrum of the measurements, obtained by using data streams where no pulses are present (more than 200 baseline windows); 2) a mean pulse m(t), obtained by averaging 100 pulses with the lowest used ∆t of 50 µs.This value being much shorter than the rise-time we can consider these events as singles (Figure 2).We then built a transfer function H(f ) accordingly and filtered each pulse in the frequency domain; the maximum A of the filtered pulse s f (t) in the time domain is then, by construction, the pulse amplitude (energy) estimator which provides the best signal-to-noise ratio (S/N).
In order to evaluate the difference in shape of each pulse (pile-up) with respect to the average pulse (single), every filtered signal is compared in the time-domain with the filtered average pulse m f (t) by fitting with s f (t) = A f it • m f (t+δt); the χ2 is minimized with respect to the A f it and δt.The minimization returns a second pulse amplitude estimator A f it .
For each signal we can determine a pulse shape parameter discriminator as PSD = A A f it : it will have value of 1 for pulses identical to that of the average pulse, lower if the shape is different as for pile-up events, and measures how much a signal differs from the reference pulse.
Particular attention was devoted to compare the pulse shape of particle-induced scintillation events with LED-generated ones; to this end we analyzed 232 Thsource calibration run (48 hours) and selected the 208 Tl (2615 keV γs) scintillation events.We observed very minor differences in the shape, whose effect is negligible with respect to the scope of this work.
To better distinguish signals differing in the risetime part it is convenient to consider the derivative of the signal instead of the signal itself.In the following we will define as PSD sig the shape parameter based on the signal itself, and PSD der on the derivative of the signal.Figure 2 shows the average signal and its derivative, whose maximum corresponds to the maximum slope of the signal in its rising part.With simple, geometrical consideration we can determine that a single event has a derivative maximum larger than that of a pile-up event at the same energy.Figure 3 shows an example of two pile-up signals separated by 0.7 ms and 3 ms, together with their derivatives.
From now on we will work exploiting PSD der instead of PSD sig .Figure 4 reports the PSD der as a function of the running time for a full measurement acquired at V e = 50 V and LED pulses injected according to the ∆t time pattern as defined in Section 2.2.A PSD cut value, shown as a pink horizontal line in Fig. 4, is determined and adjusted such as to keep 90% of singlepulse events (∆t = 0.05 ms).For each group of pile-up pulses (100 consecutive) the number of events N rej with PSD der < PSD cut is counted and compared to the total number N inj of events injected.The rejection factor is defined as r ∆t = N rej /N inj and the pile-up rejection power curves are traced for each set of measurements.Figure 5 shows all the experimental results together 2 ; for completeness, for each set the S/N ratio is Fig. 4 Scatter-plot of the PSD der as a function of the acquisition time (the detector is biased at V e = 50 V).Every 500 s the ∆t between the two pile-up pulses is changed, according to the pattern reported in the text.The distribution allows to trace the pile-up rejection curve with respect to ∆t. also given, which shows that the rejection capabilities increases with increasing S/N ratio, as also reported in [39,18], and expected from signal processing in general.

Synthetic data generation, analysis and results
In order to reproduce, control and understand our pileup experimental data, we adopted a method which permits to generate synthetic data, indistinguishable from the real, measured data [40].The method is hereafter described.
Starting from a record of data (streaming) we proceed by: 1) constructing an average signal m(t) in the exact same condition as previously performed using the experimental data (see above); 2) summing the average signal re-scaled to an amplitude A1 with an another average signal of re-scaled amplitude A2, but time-shifted of ∆t; 3) inserting the previously constructed pile-up pulses to the experimental data streaming (the full measurement was taken into account and we made sure that in the selected region no pulses were present).
The full procedure is sketched in Fig. 6.We built the synthetic data streaming conforming with our experimental measurements, namely equal amplitudes of the two signals, identical ∆t pattern and the same number of pile-up events and repetition rate.The synthetic data are then elaborated using the same analysis pipeline of the experimental data.The pile-up rejection values r ∆t are hence extracted as described in the previous section.
To benchmark the quality of the synthetic data we compared in Fig. 7 the synthetic pile-up rejection performances with the experimental one.An excellent agreement is observed in the whole range, within the uncertainties 3 , between the measurement and the synthetic data.
In the next section we will show how we performed wider synthetic data simulations by modifying the pileup time step, number of occurrences, and A1/A2 ratio to best comply with a real double-beta decay experiment.

Fig. 7
Rejection power curves versus ∆t for an electrode bias V e = 30 V. Blue dots are from the measurements (see Fig. 5), green crosses are from the synthetic data, red plain curve is from the full simulation, i.e.A1&A2 ∈ 2ν2β-spectrum, and orange dashed curve is from the approximate simulation with A1 = A2.Error bars are within the symbols. 3Bayesian statistics calculation was used to infer the error bars on the experimental data and are within the symbols.The variance for a having k out of n occurrences reads: We finally have 68% within k ± σ where σ = √ V .For a rejection of 0.1 ( k n = 0.1) and a number of occurrences of 100, σ = 0.03.They are equal as for the synthetic simulation as the number of occurrences is the same, namely 100 events for each ∆t [41].

Refined synthetic pile-up simulation
First of all we refined our synthetic data by producing pile-up events with a ∆t running from 0 to 3 ms with a time step of 0.1 ms, and injecting 1800 events for each, assuming equal amplitudes for the two signals (A1 = A2).This produces the orange curve shown in Fig. 7.
However, a realistic synthetic simulation (called hereafter full simulation) must take into account that the amplitudes of the two signals, A1 and A2, should comply with the 2ν2β energy distribution [42].Therefore, pile-up events are constructed by picking up twice randomly inside this energy distribution and summing them up with a ∆t ranging from 0 to 3 ms.We finally select those events having an energy in the ROI, i.e. 3034 keV ± 50 keV (Fig. 8).The full simulation pile-up rejection curve is shown in Fig. 7 as red solid line.It should be noticed that we expect a lower pile-up rejection efficiency than for the case where pulses have identical amplitudes (orange curve).This is due to the fact that for A1 << A2 (A1 >> A2) pile-up events may be barely distinguishable from single-pulses and may not be rejected.This is precisely what is observed in Fig. 7 as the red plain curve (A1 = A2) is below the orange dashed one (A1 = A2).The full simulation and the approximate simulation returning almost the same rejection power curves motivates the simplified approach employed in the experimental work.

Background index extraction
To extract the background index resulting from pile-up events from these simulations we need to compute the probability, Q ββ (∆t), to obtain an event at 3034 keV resulting from a pile-up of two events coming from the 2ν2β decay for different values of ∆t between the two single events.
These probabilities are shown in Fig. 9 where the distributions of the amplitude of the pile-up events are presented for ∆t = 0 and 1 ms.The rightmost distribution corresponds to ∆t = 0 ms and does not depend on the signal shape since it is identical to the initial average signal.Hence the amplitude of the signal does not vary after the optimal filtering whereas for ∆t = 0 it does.The presented distributions were obtained for an analytical signal of 1.2 ms rise-time, the description of which is presented hereafter.The leftmost curve (blue) corresponds to the amplitude of the pile-up events for ∆t = 1 ms and the middle (red) one corresponds to its distribution, slightly distorted by the optimal filtering procedure.
Fig. 9 Distributions of the sum of two 2ν2β decays separated by ∆t equal 0 and 1 ms after the optimal filtering procedure.The probability is given at 3034 keV.
We then calculated the background indexes associated to pile-ups of 100 Mo 2ν2β events for every simulated configuration (Table 2) as follows: where τ is the time step (here 0.1 ms as in the simulations), N100 M o the number of 100 Mo nuclei and is the half-life of the 100 Mo 2ν2β decay.
In the next section we study how the signal rise-time changes the pile-up rejection efficiency.

Signal rise-time and pile-up rejection efficiency
It is of major interest to investigate and predict how the light detector rise-time τ r will affect the pile-up rejection capability, combined with the S/N ratio.To this end, we generated an analytical average pulse, which is profiled by adding a Gaussian function, to account for the rising part of the signal, with a sum of exponential functions, to account for the decay part of the signal.
Figure 10 shows three analytic average pulses with different rise-time values; even though the blue curve (τ r = 1.2 ms) does reproduce fairly well the experimental average pulse m(t) (dashed gray line) extracted from the measurements, all the simulations with these rise-time were performed with the experimental average pulse.Analytical pulses enable to make faster/slower average pulses in a simple manner, to produce different synthetic simulation data and extract, accordingly, the power rejection curves r ∆t = f (S/N, τ r ).The variation of the pile-up rejection power curves with respect to the pulse decay time is less important and therefore is not taken into account in this study.
We have performed full simulations with the pulses of Fig. 10, varying for each the S/N ratio accordingly; results are displayed in Fig. 11.It should be noticed that the S/N achieved correspond to the 1.2 ms rise-time signal (as in the experiment) and for NTL-LD electrode biases V e = 10, 30, 50 V, the three left round blue dots of Fig. 11.The S/N corresponding to 70 and 90 V were extrapolated.
The method here presented is a powerful tool to predict the pile-up rejection power and hence the pileup background index achievable in the ROI, by simply combining few, generic detector specification as S/N ratio and rise-time.

Approximate evaluation of the background index
As we performed the full simulations for every bias, we were able to determine the ∆t value corresponding to 50% pile-up rejection (∆t 50% ).This latter was Fig. 11 Background index associated to 100 Mo 2ν2β decay events, after a PSD der selection, for three values of the signal rise-time and several S/N ratios (Table 2).The simulations of the blue circles corresponding to τ r = 1.2 ms were done with the average experimental signal whereas the two others, τ r = 2 ms and τ r = 0.8 ms, used the analytical signals described in Fig. 10.calculated using a linear interpolation between the two consecutive ∆t leading to rejection below and above 50%.Corresponding values are reported in Table 2 as well as the background index.The background index associated to pile-ups of 100 Mo 2ν2β events in a 280 g bolometer, 95% 100 Mo-enriched and using the PSD der parameter can be approximated as: where ∆t 50% is expressed in ms.
Because of the high slope of the rejection power curve around ∆t 50% , the above equation gives a better feeling of the BI that can be achieved than by simply extrapolating from the ∆t value where 90% rejection is obtained as often expressed [23].With this rule of thumb hereby expressed, in order to obtain 5•10 −5 ckky [25] we should reach 50% rejection for signal separated by ∆t = 0.18 ms.

Conclusions
Some future bolometric experiments, like CROSS and CUPID, rely on a detector technology based on NTDequipped scintillating bolometers to search for the 0ν2β of 100 Mo.Due to the poor time resolution of the main heat bolometer, randomly-coincident (pile-up) 100 Mo 2ν2β decay events will eventually represent the main source of background in the next-generation, large-scale, 0ν2β decay searches based on this technology.To circumvent this problem we have operated a Li 2 MoO 4 scintillating bolometer in combination with a Neganov-Trofimov-Luke light detector.Thanks to the fast time response and the enhanced S/N ratio performances of this latter, we investigated the rejection of pile-ups from single-pulse events utilizing the mere scintillation signal.We studied the pile-up rejection performance of the setup in a validation run hosted in the CROSS cryogenic underground facility.We demonstrated for first time with an experimental pile-up measurement coupled with a simulation, that it is possible to reject randomly-coincident events in massive Li 2 100 MoO 4 cryogenic scintillating bolometers, as those of CROSS and CUPID experiments, down to background indexes of ∼ 10 −4 ckky (at the Q ββ of 100 Mo), via the scintillation signal.
Together with the experimental work, we presented (i) an approximate and full simulation method to reconstruct, synthetically, randomly coincident events from the 2ν2β decay mode, (ii) a new pulse-shape discrimination algorithm based on a derivative, optimal-filtering signal processing technique capable to provide superior pile-up rejection efficiencies.The experimental and simulation pile-up rejection results have been compared: we observed an excellent agreement between them.As we have shown that there is a strong dependence with the signal rise-time, it is strongly recommended that we use the real signal such as the one coming from calibrations to extract the average signal to be used in the simulation to predict the achievable background index.
With respect to CUPID, the CROSS detectors foresee to have the additional capability to reject surface events [43].The reduction of the random-coincidence background enables to fully exploit the CROSS surface sensitive technology.In fact, the background induced by beta particles emitted by the radioactive contamination of the passive materials facing the detectors is a sub-dominant contribution of the background in the ROI that will emerge only after the mitigation of the random-coincidence component.
New investigations are ongoing to further enhance the light detector S/N ratio via higher NTL electrode bias, shorten the light detector rise-time and hopefully reach background indexes in the ROI of 5 • 10 −5 ckky, which would fully comply with the CUPID experiment goal.

Fig. 1 (
Fig.1(LEFT) Main heat bolometer, accommodated in an silver-coated copper case and operated at 15 mK.It is equipped with 2 different NTD-Ge thermistors and a 300 kΩ silicon doped heater; (RIGHT) a NTL-LD germanium detector, equipped with interconnected, concentric electrodes: it faces the heat bolometer and closes the cavity on top.It has an effective signalto-noise gain of the order of 10 when the differential of potential of the two set of interconnected electrodes is kept at 50 V.An anti-reflecting SiO coating (bluish) is deposited on the germanium wafer, between the electrodes.

Fig. 2
Fig.2Mean, normalized LED pulse m(t) and its derivative d dt m(t).A 500 ms time-window was used for the Gatti-Manfredi optimal filter analysis.

Fig. 3
Fig. 3 Example of two 0.45 keV pile-up pulses (top) and derivatives (bottom, time-shifted to have the maximum centered at t=0.03 s) recorded by the NTL-LD at V e = 50 V, for a pile-up time ∆t of 3 ms (red) and ∆t = 0.7 ms (blue), respectively.Our pulse shape parameter algorithm utilizes the derivative of the signals and allows to reject 95% of pile-up pulses with ∆t = 0.7 ms.

Fig. 5
Fig.5Experimental rejection curves traced for several V e (continuous lines are used to guide the eye).The S/N on the filtered signal s f (t) spans from 9.6 at V e = 0 V to 101.8 at V e = 50 V, corresponding to an effective NTL gain of about 10.

Fig. 6
Fig. 6 Two average signals of different amplitudes are summed together with a time difference of ∆t and with a noise sample to give the synthetic data shown in the bottom of the figure.

Fig. 8
Fig. 8 Energy spectrum of 100 Mo 2ν2β decay.The two red bars correspond to the energies of randomly picked up values which sum up to an energy in the ROI of 100 Mo 0ν2β decay.

Fig. 10
Fig. 10 Average signals coming from the fit as described in the text for rise times equal to 0.8, 1.2 and 2.0 ms.The gray dashed line corresponds to the average signal for 1.2 ms rise time.

Table 2
Background index values for different S/N ratio values and signal rise-time as shown in Fig. 11.∆t 50% are the ∆t corresponding to 50% of pile-up rejection for signals of τ r = 1.2 ms.Bias S/N BI [in 10 −4 ckky] ∆t 50% [ms] [V]for τ r = 2 ms for τ r = 1.2 ms for τ r = 0.8 ms for τ r = 1.2 ms