Pulse shape analysis in Gerda Phase II.

The GERmanium Detector Array (Gerda) collaboration searched for neutrinoless double- β decay in 76 Ge using isotopically enriched high purity germanium detectors at the Laboratori Nazionali del Gran Sasso of INFN. After Phase I (2011-2013), the experiment benefited from several upgrades, including an additional active veto based on LAr instrumentation and a significant increase of mass by point-contact germanium detectors that improved the half-life sensitivity of Phase II (2015-2019) by an order of magnitude. At the core of the background mitigation strategy, the analysis of the time profile of individual pulses provides a powerful topological discrimination of signal-like and background-like events. Data from regular 228 Th calibrations and physics data were both considered in the evaluation of the pulse shape discrimination performance. In this work, we describe the various methods applied to the data collected in Gerda Phase II corresponding to an exposure of 103.7 kg year. These methods suppress the background by a factor of about 5 in the region of interest around Q β β = 2039  keV, while preserving ( 81 ± 3 ) % of the signal. In addition, an exhaustive list of parameters is provided which were used in the final data analysis.

Introduction

Neutrinoless double- () decay is a hypothetical process in which two neutrons in a nucleus are transformed simultaneously into two protons with the emission of only two electrons. Such a process violates lepton number conservation and requires the neutrino to be its own antiparticle (Majorana particle). In combination with cosmological observations and direct neutrino mass measurement, a non-zero decay rate would highly constrain the standard light, left-handed neutrino exchange mechanism via the effective Majorana neutrino mass or shed light on alternative processes [1].

The highest half-life sensitivity to decay requires the experiments to achieve large target mass, high detection efficiency, good energy resolution and most complete elimination of background at the Q-value of the decay (). The goal of the GERmanium Detector Array (Gerda) experiment was to realize a background-free1 experiment for the first time. Gerda was located at the underground Laboratori Nazionali del Gran Sasso (LNGS) of INFN, Italy. Gerda used up to 43 kg of high purity germanium (HPGe) detectors enriched in the candidate isotope Ge up to 88%. They ensure high detection efficiency, low intrinsic background and excellent energy resolution. The bare HPGe detectors were operated in liquid argon (LAr), which served as cooling medium and as active shield against environmental backgrounds at  keV. The details of the experimental setup and its upgrade from Phase I to Phase II can be found elsewhere [2, 3].

The Phase II data taking took place between December 2015 and November 2019 with an upgrade of the detector array and the surrounding LAr instrumentation in 2018. Gerda operated three types of enriched HPGe detectors arranged in an array of 7 strings: 7 semi-coaxial detectors, referred to as coaxial detectors for brevity, from Phase I with a total mass of 15.6 kg; 30 Broad Energy Germanium (BEGe) detectors (20.0 kg) [4, 5]; and 5 inverted coaxial (IC) detectors (9.6 kg) [6, 7], which were installed in summer 2018. The accumulated exposure, product of total detector mass and respective livetime, amounts to 103.7 kg year for Phase II. In order to avoid bias in the event selection criteria, Gerda followed a strict blinding strategy, where events within a  keV energy window were processed only after the analysis had been finalized.

Signal formation, readout and processing in germanium detectors

Gerda used p-type HPGe diodes of three different geometries called coaxial, BEGe and inverted coaxial (Fig. 1). They all have a relatively thick ([0.8–2.6] mm) n electrode, formed by lithium diffusion, and a thin (300 nm) p electrode created by boron implantation. The p and n electrodes are separated by a groove with non-conducting surface. Because of their small p electrode the BEGe and IC detectors belong to the class of point-contact detectors [12] which exhibits intrinsic performance advantages with respect to energy resolution and pulse shape discrimination.

Fig. 1

Geometry and weighting potential of a typical BEGe, an inverted coaxial and a coaxial detector. The p and n contacts are indicated schematically

The Ge detectors are operated under reverse bias voltage such that almost the entire volume is depleted of free charge carriers. An interaction in the active volume creates a number of electron-hole pairs proportional to the deposited energy. The charge carriers drift according to the electric field created by both the positive bias voltage applied to the n contact and the volume charge density due to the net bulk impurity concentration. The electrons are collected at the n contact, the holes at the p contact which is used for readout. The n layer covers most of the crystal surface and features a large Li concentration. Exhibiting zero electric field beyond the p-n junction till the outer surface, a charge created in this layer will only experience thermal diffusion with two possible outcomes: recombine (loss) or reach the depleted volume (collection). As a result, two generic regions can be identified, the dead layer with high probability of no charge collection and a transition layer with partial charge collection.

During the drift of charge carriers, charge is induced on the readout contact as described by the Shockley–Ramo theorem [13]:where is the total charge carried by the holes or electrons and is the weighting potential along the drift path of holes or electrons. The weighting potential is shown in Fig. 1 for the three different detector geometries used in Gerda also indicating the geometry of the p and n surfaces. The weighting potentials have been calculated using the AGATA Detector Library [14] pulse shape simulation package.

Due to their small p contact, BEGe and IC detectors have a weighting potential distribution that is very small in most of the volume and sizable only close to the p contact. This results in similar waveforms Q(t) from interactions in a large part of the volume. Multiple energy deposits can be treated as a superposition of single interactions. In Fig. 2, normalized example pulses from a BEGe detector are shown for a SSE, a MSE, an event close to the p contact and an event near the dead layer of the n contact with incomplete charge collection. As shown, surface events produce characteristic pulse shapes being fast close to the p contact due to the strong electric field and slow near the n contact due to the weak electric field and the transition layer.

Fig. 2

Examples of different type of normalized charge pulses (red) along with the derived current pulses (blue) from a BEGe detector: single-site event, multi-site event, event near the p contact and event near the n contact with incomplete charge collection

In coaxial detectors, both the hole and electron drift play a role in the pulse formation, which result in different pulse shapes throughout the volume of the detector. In Fig. 3, simulated example pulses from different parts of the detector are shown. Similarly to BEGe detectors, coaxial detectors also show special pulse shape characteristics in case of surface events. Indeed, energy deposits near the groove or the bottom of the borehole cause faster pulses.

Fig. 3

Simulated charge pulses in indicated parts (bulk, groove, bottom/side of borehole) of a typical coaxial detector. The electronic response of the readout chain and noise are not included here

The signals induced on the p contact of the Gerda detectors are read out by charge sensitive amplifiers located in the LAr about 35 cm above the detector array. The signals are digitized at 25 MHz for s and at 100 MHz for 10 s. Both traces are centered at the rising edge of the charge pulse Q(t). The offline analysis of the digitized signals follows the procedures described in [15, 16]. The 25 MHz traces used for the energy reconstruction ensure the excellent energy resolution, while the 100 MHz traces are used in the PSD methods presented in the following sections. The energy estimator E is reconstructed with a zero-area cusp filter [17], whose parameters are optimized for each detector and calibration run.

Overview of event samples and discrimination methods

Weekly calibration runs with Th sources are performed to determine the energy scale and resolution of the detectors [18] and to calibrate and train the PSD techniques. Figure 4 shows a calibration spectrum highlighting the different event samples used in pulse shape analysis. The most prominent feature is the full energy peak (FEP) at 2615 keV from Tl decay. Its double escape peak (DEP) at 1593 keV is used as a sample of SSEs as the electron and positron from pair production deposit their energy in a small volume and both annihilation rays leave the detector. The FEP at 1621 keV from Bi is used as a sample of MSEs that is sufficiently near in energy to the DEP in order to avoid noise dependent biases. To test the performance of MSE rejection, the FEP and single escape peak (SEP) of Tl, mostly featuring MSEs are used while the Compton continuum region around  keV (CC()) serves to estimate the background rejection in the decay signal region.

Fig. 4

Calibration spectrum highlighting the different event samples used in the pulse shape analysis

In physics data, the standard neutrino accompanied double- () decay provides another sample of signal-like events that is equally distributed in the whole detector volume and used for the investigations in the following. After applying the LAr veto cut [3], about 97% of the events in the 1000–1300 keV region originate from decays. Beside MSEs from rays, the physics data have a significant amount of surface events from and decays that can be discriminated thanks to their specific pulse shape. The physics spectrum at low energies is dominated by decays of Ar up to its Q-value of 565 keV. However, these events are not used in the pulse shape analysis due to their relatively high noise. A prominent background source at is the decay of K, which is produced as a progeny of the long-lived Ar and has a Q-value of 3525 keV. Beta particles deposit their energy in germanium within a few mm resulting in events partly in the dead and transition layers of the n surface. Such n surface events can induce slow pulse shapes with incomplete charge collection. Apart from possible HPGe bulk contamination, that have been shown to be insignificant [19, 20], particles can only reach the active volume of the detector at the thin p contact or at the non-conducting groove producing pulse shapes with fast rise as shown in Figs. 2 and 3. A clean sample of surface events is found in physics data above the K Q-value. The most prominent structure at these energies is a broad peak at 5304 keV, the energy of the Po decay (U decay chain) [19, 21].

Due to their different geometries, BEGe and IC detectors are treated separately from coaxial detectors in the pulse shape analysis. In the case of the BEGe and IC detectors one parameter, A/E, is used to classify background events, where A is the maximum current amplitude as indicated in Fig. 2 and E is the energy. As MSEs and surface events at the n contact are characterized by longer, i.e. wider current pulses, they feature a lower A/E value compared to SSEs, while surface events at the p contact show a higher A/E value [22]. Therefore, rejecting events on both sides of the A/E distribution of SSEs enhances the signal to background ratio. The details of the A/E analysis are presented in Sect. 4. In the case of coaxial detectors an artificial neural network (ANN [23]) is used to discriminate SSEs from MSEs similar to the approach applied in Phase I [10]. To discard events close to the p contact, a dedicated cut on the risetime of the pulses is applied. The training and optimisation of the ANN and risetime cuts for coaxial detectors is described in Sect. 5. An additional cut is applied to all detectors to remove events with slow or incomplete charge collection. These events arise from energy depositions in a non-depleted volume (n layer or insulating groove). These events are identified through the difference between two energy estimates performed using the same digital filter but different shaping times as summarized in Sect. 6. The signal efficiency of decay is estimated using the survival fraction of DEP and decay events by taking into account the energy dependence of the different PSD techniques. The details of this extrapolation to including systematic uncertainties are found in Sects. 7 and 1.

The A/E method for BEGe and IC detectors

The amplitude of the current pulse A is computed after applying 3 times a moving window average (MWA) filter with 50 ns length and interpolating the pulse down to 1 ns sampling time. This filtering procedure optimizes the high frequency noise attenuation while preserving the pulse shape information. The energy estimator E is determined by a pseudo Gaussian filter with a shaping time of 10 s. A is then divided by E before calibration, providing the raw A/E for each pulse. The raw A/E is then corrected for time stability and energy dependence before a cut value is defined.

For each calibration run, the A/E distribution of events in the 1000–1300 keV region is fitted with a Gaussian (SSEs) and a low-side tail (MSEs) as described in [10]. The position of the Gaussian from each calibration of the four years of data taking is used to define stable time periods where the raw A/E changes by less than its resolution. Instabilities are mostly related to hardware changes and a few detectors show a small systematic drift of the raw A/E. Physics data between the stable periods are removed from the analysis causing a few percent exposure loss. After normalizing the raw A/E by the average within a given time period, the data of all calibrations are merged, only separating before and after the upgrade of 2018.

The A/E of SSEs in the Compton continuum of the merged calibration data show a small linear energy dependence of a few percent per MeV, due to the larger charge cloud size at higher energies that broadens the current pulse. The energy dependence of and is described by a linear and a type of function, respectively, as shown in Fig. 5 for a BEGe (GD61A) and an inverted coaxial (IC74A) detector as examples. In addition to the correction for the energy dependence, A/E is normalized to the mean of the A/E distribution of the DEP, which lies about 0.25% above the SSE band.

Fig. 5

Distribution of A/E as a function of energy from Th calibration data for the BEGe detector GD61A (top) and the IC detector IC74A (bottom). The dashed lines indicate the linear energy dependence and the 1- width of SSEs

The cut values for each detector are determined on the energy-independent A/E classifier defined as . Its distribution is centered around zero and has a standard deviation of one for SSEs. The low-side A/E cut against MSEs and n surface events is chosen to yield a DEP survival fraction of 90%. The resulting cut values range from to and from to for BEGe and IC detectors, respectively. The high-side A/E cut against p surface events is chosen at for each detector, in order to reject all events in physics data above 3525 keV. It has been shown that the high-side A/E cut discards events, including degraded events, from a small volume around the p contact [24] causing the survival fraction of events after the high-side A/E cut to be proportional to the detector mass.

Figure 6 shows for BEGe and IC detectors the energy distribution of calibration data before and after the A/E cuts described above as well the corresponding survival fractions. As expected the survival fraction of events in the FEPs and SEP is much lower than in the DEP. Events in the Compton continuum are discarded with about the same probability independent of their energy but depending on the overall detector size, and more generally speaking from the detector type. IC detectors discard a higher fraction of events because of the higher probability of multiple scattering of rays due to their larger size.

Fig. 6

Energy distribution of events from calibration data before and after the A/E cut (top) and their ratio (bottom) for BEGe and IC detectors

The survival fractions of events in the DEP, FEP and CC() are shown in Fig. 7 for each detector before and after the upgrade. The DEP survival fractions are slightly smaller than 90% due to the high A/E cut. The rejection of MSEs shows a small dependence on the detector position in the string because of different electronic noise conditions. This effect was reduced after the upgrade. The IC detectors (detector numbers above 35) reject MSEs more efficiently than BEGe detectors.

Fig. 7

Survival fraction of events in the Tl DEP, Bi FEP and CC() for each detector (see Table 3 for detector numbers and types). Open (filled) symbols show the calibration dataset before (after) the upgrade. The dashed lines separate the detector strings in the array. The uncertainties are only statistical and smaller than the markers

Table 3

Detector-wise PSD cut efficiencies and the total detection efficiencies for decays used in the final analysis of [26] and corresponding exposures . For the period after the 2018 upgrade also the electron containment efficiencies are listed (see text). Quoted uncertainties account for statistics and systematics. The channels with empty entries were used in anti-coincidence only

Dec 2015–May 2018					July 2018–Nov 2019
#	Detector label	E (kg year)	ε_PSD	ε_0νββ	ε_fep	E (kg year)	ε_PSD	ε_0νββ
0	GD91A	1.167	0.890(29)	0.610(25)	0.895(2)	0.736	0.891(23)	0.612(22)
1	GD35B	1.455	0.895(27)	0.634(24)	0.899(2)	0.880	0.891(30)	0.632(26)
2	GD02B	1.181	0.878(22)	0.598(21)	0.891(2)	0.733	0.885(22)	0.603(21)
3	GD00B	1.193	0.859(27)	0.583(24)	0.893(2)	0.753	0.885(21)	0.601(21)
4	GD61A	1.340	0.876(21)	0.606(21)	0.899(2)	0.858	0.901(22)	0.624(21)
5	GD89B	0.639	0.869(21)	0.572(21)	0.885(2)	0.707	0.894(22)	0.587(22)
6	GD02D	–	–	–	–	–	–	–
7	GD91C	0.234	0.867(23)	0.592(22)	0.892(2)	0.736	0.894(29)	0.611(25)
8	ANG5	5.067	0.760(32)	0.486(37)	–	–	–	–
9	RG1	3.914	0.670(31)	0.465(39)	0.919(2)	2.392	0.661(32)	0.463(39)
10	ANG3	4.557	0.653(29)	0.449(39)	0.920(2)	2.711	0.681(32)	0.472(40)
11	GD02A	1.016	0.904(23)	0.621(21)	0.889(2)	0.640	0.886(24)	0.609(22)
12	GD32B	1.257	0.909(95)	0.620(67)	0.890(2)	0.780	0.90(20)	0.62(13)
13	GD32A	0.489	0.865(22)	0.583(23)	0.883(2)	0.529	0.886(34)	0.597(29)
14	GD32C	1.377	0.878(23)	0.609(21)	0.898(2)	0.872	0.898(57)	0.624(42)
15	GD89C	1.056	0.870(22)	0.581(23)	0.880(2)	0.698	0.894(27)	0.597(25)
16	GD61C	1.072	0.874(28)	0.594(24)	0.889(2)	0.734	0.883(24)	0.602(22)
17	GD76B	0.698	0.835(21)	0.538(20)	0.879(2)	0.426	0.858(22)	0.553(21)
18	GD00C	1.479	0.892(21)	0.618(22)	0.890(2)	0.956	0.902(22)	0.618(22)
19	GD35C	1.202	0.916(39)	0.635(30)	0.889(2)	0.739	0.895(44)	0.620(34)
20	GD76C	1.371	0.902(43)	0.614(32)	0.899(2)	0.964	0.911(47)	0.621(35)
21	GD89D	0.945	0.879(32)	0.576(26)	0.880(2)	0.617	0.891(51)	0.584(37)
22	GD00D	1.550	0.889(21)	0.613(21)	0.898(2)	0.954	0.902(31)	0.622(26)
23	GD79C	1.206	0.877(21)	0.596(20)	0.896(2)	0.953	0.905(21)	0.615(20)
24	GD35A	1.460	0.882(20)	0.618(21)	0.901(2)	0.901	0.896(22)	0.630(22)
25	GD91B	0.420	0.876(50)	0.601(37)	0.893(2)	0.734	0.879(35)	0.603(28)
26	GD61B	1.172	0.896(22)	0.615(21)	0.895(2)	0.808	0.899(21)	0.616(21)
27	ANG2	4.750	0.743(32)	0.504(41)	0.924(2)	2.751	0.740(31)	0.508(40)
28	RG2	4.017	0.630(30)	0.400(34)	0.918(2)	2.542	0.640(27)	0.412(32)
29	ANG4	4.521	0.683(32)	0.477(40)	0.920(2)	2.783	0.707(31)	0.498(39)
30	GD00A	0.945	0.885(22)	0.598(22)	0.884(2)	0.582	0.864(24)	0.584(22)
31	GD02C	1.462	0.899(22)	0.618(21)	0.898(2)	0.925	0.904(21)	0.623(21)
32	GD79B	0.861	0.896(39)	0.608(31)	0.894(2)	0.469	0.900(24)	0.612(23)
33	GD91D	1.003	0.858(44)	0.588(34)	0.895(2)	0.813	0.894(44)	0.614(34)
34	GD32D	1.232	0.874(23)	0.617(22)	0.897(2)	0.775	0.873(25)	0.618(23)
35	GD89A	0.976	0.864(21)	0.585(21)	0.888(2)	0.591	0.883(52)	0.597(39)
36	ANG1	1.812	0.649(33)	0.403(36)	–	–	–	–
37	IC50B	–	–	–	0.920(2)	2.213	0.894(18)	0.664(14)
38	IC48A	–	–	–	0.921(2)	2.159	0.910(18)	0.676(14)
39	IC50A	–	–	–	0.910(2)	1.814	0.896(17)	0.646(14)
40	IC74A				0.921(2)	2.398	0.900(15)	0.652(12)

In order to check the validity of the A/E corrections and cuts in the whole dataset, survival fractions of the usual event samples are studied for each calibration run. The average survival fraction from all BEGe detectors is shown for each calibration in Fig. 8 as a function of time. The A/E cut shows, for both BEGe and IC detector types, a stable behaviour at the 3% relative level during the whole data collection period when applied to Compton continuum events at . Residual instabilities stem on one side from changes in the detector gain or leakage current and on the other side from statistical fluctuations.

Fig. 8

Average survival fractions of events in the Tl DEP, Bi FEP and CC() for BEGe detectors as a function of time. Each data point represents a calibration run with its statistical uncertainty

The raw A/E is corrected in the same way for physics data and the same cut values are applied as in calibration data. In order to check the validity of the corrections and the cut values, Fig. 9 shows a comparison of the A/E classifier between DEP events from calibration data and events from physics data for BEGe and IC detectors. By construction, the A/E classifier peaks at 0 and has a width of about 1 for these SSEs. The agreement between physics and calibration data is satisfactory and confirms the applied correction procedure. For each detector, the residual difference of the average A/E, between calibration and physics data, is included in the systematic uncertainties by shifting the cut value accordingly.

Fig. 9

A/E classifier distributions after LAr veto comparing events in the DEP from calibration data and decay events from physics data from BEGe (left) and IC (right) detectors. The histograms are normalized to their integrals

The energy distributions of events before and after the A/E cut are shown in Fig. 10 for the whole physics data corresponding to 61.8 kg year exposure from BEGe (53.3 kg year) and IC (8.6 kg year) detectors. Events in the decay region survive the cut with a high probability, while the K and K peaks at 1461 keV and 1525 keV, respectively, mostly featuring MSEs, are reduced significantly. High energy events above 3525 keV coming from p surface events are all discarded by the high A/E cut by definition. In a 240 keV wide window around , only 7 events in BEGe detectors and 1 event in IC detectors survive the LAr veto and A/E cuts. This results in the corresponding unique background indices2 of Gerda in the region of interest of  counts/(keV kg year) and  counts/(keV kg year) for BEGe and IC detectors, respectively.

Fig. 10

Energy distributions of all Gerda Phase II physics data from the BEGe and IC detectors before (grey) and after (blue) the A/E cut. The position of is indicated. The prominent K and K lines are located at 1461 keV and 1525 keV, respectively

The ANN and risetime methods for coaxial detectors

The 100 MHz waveform traces are used to compute the artificial neural network input variables (IVs) and the risetime. They are first filtered with a MWA of 30 ns width three times. They are subsequently baseline subtracted and normalized between 0 and 1 by the amplitude of the 25 MHz pole-zero corrected traces. This amplitude is provided by a trapezoidal filter with a typical precision of 0.2% at 2 MeV.

The ANN IVs are a list of the 50 times at which the resulting waveform reaches [1%, 3%, , 99%] of its amplitude (see Fig. 11). An interpolation of the 10 ns wide gaps between data points allows for a more precise estimation. These IVs are computed for all physics and calibration events found in Gerda. However, given the degraded signal-to-noise ratio at low energy, only events above 1000 keV are considered in this pulse shape analysis to avoid loosely constrained energy dependence correction. Calibration runs are used to optimize the discrimination of SSEs from MSEs. The network is built on two hidden layers with 50 and 51 neurons, using the TMVA-MLPBFGS algorithm [23].

Fig. 11

Zoom on a typical normalized 100 MHz trace of a coaxial detector. The 50 ANN input variables (red circles) and rise time estimates (dashed green) are indicated

Figure 12 shows as an example the classifier distributions from the ANN training of the ANG5 detector with the events from the indicated DEP and FEP peaks. The distributions from the Compton events under the peaks are statistically subtracted using the distributions of the events in the energy side-bands of the peaks. The lower and upper side-bands are defined by selecting events falling in the [,] and [,] energy regions where is the full width at half maximum used to quantify the energy resolution of the Ge detectors. The indicated ANN cut keeps 90% of the events in the DEP peak.

Fig. 12

Trained ANN classifier values for events in the Tl DEP (SSEs) and Bi FEP (MSEs) from the ANG5 detector. The histograms are normalized to their integrals

Due to the significant change in hardware, the data taken before and after the 2018 upgrade periods have been trained separately. Finer splittings of the data have yielded signal efficiencies and background rejection values which agreed on the one-percent level. As a result of the limited statistics, a minimal number of two datasets was preferred. Typically, about 10,000 and 15,000 events enter the signal (DEP) and background (Bi FEP) samples, respectively, for the ANN training of each detector. Similarly to the A/E method, the ANN cut for each detector and each training period is set on the Tl DEP classifier distribution such that 90% of the events survive.

The risetime is estimated after interpolating the waveform with a 1 ns time step. From studies on the rejection of particles [25], the [10–90%] amplitude signal risetime was selected (see Fig. 11). This parameter is used to reject events on the p contact that develop faster signals (see Fig. 3). The cut definition relies on the maximization of the following figure of merit:where is the survival fraction at risetime cut x and is the corresponding survival fraction of events. Only physics data after ANN-MSE and LAr veto, to increase purity of samples, are used for this figure of merit that allows to reject most of the particles while preserving a high decay signal survival fraction. Figure 13 depicts an example of such a cut definition. On average, about 90% of the high energy particles are rejected.

Fig. 13

Survival fractions after the rise time cut for decay and particles events from the Phase II data before the upgrade of the ANG4 detector. Also shown are the figure of merit and the chosen cut value

The calibration energy spectrum of all coaxial detectors before and after applying the PSD cuts is shown in Fig. 14. As targeted, the ANN cut removes preferentially the regions highly populated by MSEs (FEPs and SEP in particular) and preserves 90% of the Tl DEP. On the contrary, the risetime cut deployed to reject events with fast risetime is insensitive to these types of events and hence has a high survival fraction for both SSEs and MSEs.

Fig. 14

Energy distributions of calibration data events from the coaxial detectors before and after the ANN and risetime cuts (top) and the corresponding survival fractions (bottom)

Various survival fractions for all coaxial detectors, after applying both ANN and risetime cuts for the two data taking periods, are depicted in Fig. 15. In addition, the stabilities of these cuts in three calibration energy regions are plotted in Figs. 16 and 17. Apart from a slight improvement of the ANN rejection of the Compton continuum events at after the 2018 upgrade due to an improved signal cable management, an overall 3% level stability in PSD performance is observed over the course of Gerda Phase II. Compared to Gerda Phase I [10], on average a 7% relative worsening is observed, mostly attributed to the different electronics scheme, hence different noise. The risetime cut also shows a very stable behavior during this period.

Fig. 15

Survival fractions of events in the Tl DEP, Bi FEP and CC() for each coaxial detector after ANN and risetime cuts. Open (filled) symbols show the calibration dataset before (after) the 2018 upgrade. The uncertainties are only statistical and smaller than the markers

Fig. 16

Average survival fractions of events after the ANN cut in the Tl DEP, Bi FEP and CC() for coaxial detectors as a function of time. Each data point represents a calibration run with its statistical uncertainty

Fig. 17

Average survival fractions of events after the risetime cut in the Tl DEP, Bi FEP and CC() for coaxial detectors as a function of time. Each data point represents a calibration run with its statistical uncertainty

The result of the pulse shape analysis applied on the full 41.8 kgyr physics data exposure, after the LAr veto cut, is summarized in Fig. 18. The ANN is preserving 80.5(3)% of the decay event sample while removing 62(1)% of the K and K lines. The reduction of the Bi FEPs at 1806 keV and 2204 keV is also visible. However, it is in general inefficient at rejecting fast and degraded events from about 3500 keV to 5500 keV, as 58.4(8)% of them survive. After applying also the risetime cut, 67.8(4)% of decay events remain. The background is suppressed by a factor of 23 above 3500 keV, only 93 events out of 2169 survive the PSD cuts. The complementarity between the two cuts yields for the coaxial detectors a dataset that remains background-free in the region of interest with a background index of  counts/(keV kg year).

Fig. 18

Top: ANN classifier of physics events from coaxial detectors as a function of energy. Middle: Risetime of events surviving the ANN cut as a function of energy. Bottom: Energy distribution of events from coaxial detectors before and after the ANN and risetime cuts

Events with incomplete charge collection

Events from the n layer or the groove featuring slow or incomplete charge collection (see Fig. 2) have an uncertainty on energy reconstruction because the ZAC filter [17] is optimized for the FEP resolution with a relatively short integration time. These particular events, especially for coaxial detectors, can survive the ANN and risetime cuts. In order to discard these events with uncertain ZAC energy reconstruction, an additional classification based on energies reconstructed with pseudo Gaussian filters with different integration times is performed. The energy is reconstructed with an integration time of 4 s () and 20 s (). The ratio is then normalized to its average observed in events from the 2615 keV line in calibration data. The classifier is defined as:where E is the default ZAC-reconstructed energy for each event. With this definition, the classifier has an average of 0 keV.

The / distribution normalization is performed time dependently in order to account for possible instabilities of the readout electronics: for each calibration run the mean of the fitted Gaussian to the FEP distribution is used as a normalization factor for the following physics events.

Figure 19 shows the resulting distributions for calibration events in the Tl FEP and in the CC() as well as for physics events in the region for the ANG2 detector as an example. Events in the Compton continuum have a higher fraction with large negative values. This is due to the higher fraction of pulses with incomplete charge collection that is not present in a peak where the whole energy has been collected. Physics data in the region shows a similar Gaussian distribution as calibration data with no significant energy dependence.

Fig. 19

Distribution of the classifier for calibration events in the Tl FEP and in the Compton region around as well as for physics events in the region for the coaxial ANG2 detector. The histograms are normalized to their integrals

A cut value is applied to the lower side of the distribution and it is set for each detector separately away from zero.3 The cut value is loose enough that more than 99% of signal events are kept whereas those with uncertain energy (significant difference between energies reconstructed by the ZAC and Gauss algorithms) are rejected.

For the BEGe and IC detectors, most of the decay events rejected by the cut are also rejected by the A/E cut. For BEGe detectors out of 3477 events cut by either method only 7 are rejected by the cut (for energies above 1000 keV 12 events in total are cut by ). For the IC detectors none of the 521 decay events is removed by the cut alone.

For the coaxial detectors the correlation between ANN, risetime and cuts is weaker. Analyzing the decay region for the full Phase II dataset, one gets 4970 events (out of 15,433) removed by either method, 1 event cut by all three methods and 57 by only. 83 events from coaxial detectors with  keV are cut by only. Many of these events show slow pulses that could originate from the detector surface but are not cut by ANN or risetime cuts.

The survival fractions for different energy regions of physics and calibration data have been studied for each detector separately. Survival fractions of decay events for the different detector types are presented in Table 1 before and after the other PSD cuts. Values for DEP, FEP, CC() and decay events for the two datasets before and after the 2018 upgrade are very close to 100%. The Compton region as well as the region show lower acceptance because of the contribution of slow pulses.

Table 1

Survival fractions of decay events for the cut without and in combination with other PSD methods (A/E or ANN+risetime). The uncertainties are statistical only

Detector type	Survival fraction [%]
	Before PSD	After PSD
Coaxial	99.57 ± 0.05	99.46 ± 0.07
BEGe	98.47 ± 0.09	99.96 ± 0.02
IC	98.58 ± 0.20	>99.90 (95% C.I.)

For BEGe and IC detectors the impact of the cut is negligible, while for coaxial detectors a small correction of the efficiency has to be taken into account.

PSD detection efficiencies at Q

In the absence of signal proxies at in calibration and physics data and due to a not sufficiently accurate modeling of the pulse shape response in simulation, is estimated for point contact and coaxial detectors from the extrapolation of the survival fractions of the Tl DEP and of the decay events at 1593 keV and 1150 keV to (see circles and squares in Fig. 20).

Fig. 20

Extrapolation of reference survival fractions (circles) to (squares) using the energy dependence deduced from indicated two samples of down-scaled waveforms (see text). Examples are given for BEGe (top), coaxial (middle) and IC detectors (bottom) for the A/E, ANN and methods, respectively. As to the BEGe-A/E example on top: the reference survival fraction is corrected for Compton events below the Tl DEP while for the rescaled distribution (red) this correction is missing

For that purpose, sets of down-scaled waveforms are used to evaluate the energy dependence of the various PSD methods. This is based on the observation that the energy dependence of PSD classifiers is dominated by the electronic noise in Gerda. Two samples from calibration data have been considered, namely Tl DEP and CC () events (see Fig. 4). The waveform at energy E is produced by down-scaling the waveform , of true energy , from one of these two samples and superimposing it to a baseline sample b(k) obtained from random triggers:where k is the index running from 0 s to 160 s and from 0 to 10 s for the low and high frequency traces introduced in Sect. 2, respectively. By applying the complete set of digital processing steps and analysis methods described earlier, such a procedure allows to estimate the contribution of the energy dependence, induced by the signal-to-noise ratio variation, on the PSD methods results (for more details see [25]).

Figure 20 shows the survival fractions of the Tl DEP and CC() datasets as a function of the down-scaled energy for the A/E, ANN and methods.4 Each plot depicts the energy dependence of a given detector type, as a result of a weighted average of the individual detector responses, to illustrate the general trends. The reference efficiency is extrapolated to using the average energy dependence of Tl DEP and CC() datasets.5 Table 2 shows how the different PSD procedures contribute to the resulting average detection efficiencies of the various detector types for each data taking period. The efficiencies of the point contact detectors are almost 30% larger than those of the coaxial detectors.

Table 2

Average signal detection efficiencies at of individual PSD methods for the different detector types and data taking periods. The corresponding total PSD detection efficiencies and their systematic uncertainties are shown in bold; they are estimated via a MC sampling of individual values (see Table 3) and are thus different from the product of the individual PSD methods efficiencies reported here

	Dec 2015–May 2018		July 2018–Nov 2019
	Coaxial	BEGe	Coaxial	BEGe	Inverted coaxial
Exposure E	28.6 kg year	31.5 kg year	13.2 kg year	21.9 kg year	8.5 kg year
ε_ANN	82.5%	–	81.8%	–	–
ε_risetime	85.7%	–	85.0%	–	–
ε_A/E	–	88.4%	–	89.3%	90.0%
ε_δE	99.5%	100.0%	99.7%	100.0%	99.7%
ε_PSD	(69.1 ± 5.6)%	(88.2 ± 3.4)%	(68.8 ± 4.1)%	(89.0 ± 4.1)%	(90.0 ± 1.8)%

The analysis of the final Gerda results [26] is based on time-dependent and detector-wise datasets. Table 3 (Appendix A) reports the overall PSD efficiency for each detector separately. Since the central values of the individual detector efficiencies exhibit significant shifts, the average overall PSD efficiencies of Table 2 have been obtained, in fact, from an exposure weighted Monte-Carlo sampling, and not from a simple average.

The systematic uncertainty of the extrapolation is estimated from the difference of the slopes of the two down-scaled waveform samples. It is on average about 0.5% and 1.3% for point contact and coaxial detectors, respectively, the latter detector type being more sensitive to the noise due to its larger p contact capacitance.

In addition, two other systematic effects have been taken into account: (1) the difference between the calibration and physics data and (2) the difference between signal proxies and decay events. The former applies only to point contact detectors as the Tl DEP from calibration data is used as signal proxy. Shifting the A/E cut, for each detector, by the A/E distribution bias observed between the Tl DEP and decay (Fig. 9) leads to an average relative systematic uncertainty of 1.9%. The latter makes use of pulse shape simulation (see Appendix B) to quantify the PSD performance bias between decay and Tl DEP events and those coming from decay. Indeed, the signal proxies feature much lower energies hence different Bremsstrahlung contribution. In addition, Tl DEP events have a higher probability to happen close to the detector surface while decays would occur homogeneously throughout the detector bulk. A 2.3% and 1.5% absolute uncertainty has been estimated for BEGe and IC detectors, respectively, while it amounts to 4% for coaxial detectors. This last estimate is larger due to the difficulty to match the ANN training performance in simulation with the data.

Summary

Nowadays, running a background-free decay experiment is essential to boost the sensitivity on a reasonable time scale. Over the past years, the Gerda collaboration demonstrated the feasibility of such a program by upgrading its initial setup with additional point contact detectors (BEGe and IC), a LAr veto instrumentation and lower mass holders. As a consequence, the sensitivity linearly increased with the exposure. The interplay between passive and active shielding techniques has proven to be highly effective. In this paper, we focused on the ability offered by germanium detectors to analyze the topological structure of the recorded events. This topology, distinct for signal-like and background-like events at a 100 ns time scale, is best scrutinized with the high frequency based data acquisition system of Gerda. Using the Phase II dataset, we confirmed the superior discriminating power of point contact detectors (BEGe and IC) over the historical coaxial detectors thanks to the simple A/E parameter. Their decay PSD efficiency is 89% and 69%, respectively for a similar background index of about  counts/(keV kg year) after all cuts. We also demonstrated the high and stable performance in LAr of the newly produced enriched inverted-coaxial point contact germanium detectors. The LEGEND collaboration will deploy these IC detectors for its future Ge decay search program.

The total detection efficiencies for decay are discussed and listed for each detector in Appendix A.

Table 4

Active volume fractions of IC detectors used in the final analysis of [26]. The values are pulled from a preliminary analysis found in [27]. Quoted uncertainties account for statistics and systematics

IC detectors
#	Detector label	f_av
36	IC48B	0.935(5)
37	IC50B	0.938(4)
38	IC48A	0.936(7)
39	IC50A	0.920(4)
40	IC74A	0.913(6)