Mu2e Crystal Calorimeter Readout Electronics: Design and Characterisation

Abstract

The Mu2e experiment at Fermi National Accelerator Laboratory will search for the charged-lepton flavour-violating neutrinoless conversion of negative muons into electrons in the Coulomb field of an Al nucleus. The conversion electron with a monoenergetic 104.967 MeV signature will be identified by a complementary measurement carried out by a high-resolution tracker and an electromagnetic calorimeter, improving by four orders of magnitude the current single-event sensitivity. The calorimeter—composed of 1348 pure CsI crystals arranged in two annular disks—has a high granularity, 10% energy resolution and 500 ps timing resolution for 100 MeV electrons. The readout, based on large-area UV-extended SiPMs, features a fully custom readout chain, from the analogue front-end electronics to the digitisation boards. The readout electronics design was validated for operation in vacuum and under magnetic fields. An extensive radiation hardness certification campaign certified the FEE design for doses up to 100 krad and 10${}^{12}$ n${}_{1\mathrm{MeVeq}}$/cm${}^2$ and for single-event effects. A final vertical slice test on the final readout chain was carried out with cosmic rays on a large-scale calorimeter prototype.

1. Introduction

The Mu2e experiment [1] at Fermilab will search for the charged-lepton flavour-violating (CLFV) process of the muon to electron conversion in the Coulomb field of aluminium nuclei, aiming to improve by four orders of magnitude the current single-event sensitivity, by means of a high-intensity pulsed muon beam and high-performance detectors—a high-resolution (<200 keV/c) straw-tube tracker in conjunction with an electromagnetic calorimeter (EMC)—to identify the distinctive signal signature, represented by a monoenergetic conversion electron of 104.96 MeV, slightly below the muon mass. The EMC complements the tracker providing high rejection power in $\mu$/e separation, cluster-seeding for tracking algorithms, while also providing stand-alone online trigger capabilities. The calorimeter needs to maintain high levels of performance and stability while operating in a high vacuum, in the presence of a 1 T magnetic field and in a harsh radiation environment. The EMC has high granularity, good energy (<$10\%$) and timing (<500 ps) resolution for 100 MeV electrons [2]. It is composed of two disks, each housing a matrix of 674 parallelepipedonal ($34\times 34\times 200$ mm${}^3$) undoped CsI crystals. Each is read out by two large-area arrays of UV-extended silicon photomultipliers (SiPMs) and a fully custom signal chain. The design, validation and qualification steps relative to the calorimeter front-end electronics are reported in the following sections, while a detailed description of the Mu2e experiment, its calorimeter and the associated construction status can be found elsewhere in the current issue [3]. The Mu2e experiment is currently under construction at Fermilab. Detector commissioning and installation in the experimental hall will begin in late 2023.

2. Readout and Signal Chain Architecture

The Mu2e calorimeter employs a fully custom readout chain for 2716 channels. Behind the crystal matrices, 1348 readout units (ROU) are installed. A total of 136 mezzanine boards (MB) and 136 digitiser readout controllers (DiRAC) are installed in 20 crates embedded in the calorimeter’s mechanical structure (Figure 1).

2.1. Readout Units

Each ROU houses two fully independent readout channels per crystal, each comprising two custom large-area UV-extended SiPMs and the relative front-end electronics (FEE) boards. The SiPMs are glued on copper holders (Figure 2, left) which are thermally coupled to the cooling system (Figure 1, top-left) for their thermalisation down to −10 °C, to cope with the nonionising-dose-induced dark current increase.

A 2 mm gap is present between crystals and SiPMs to reduce their thermal coupling, preventing gradients along the crystal axis which could worsen their longitudinal response uniformity. A fibre-optic holder (Figure 2, bottom-left) allows for the distribution of 520 nm picosecond laser pulses for photosensor gain monitoring and timing alignment. A copper Faraday cage completes the assembly.

2.2. SiPMs

Two adjacent custom SiPMs by Hamamatsu [4] are used for each crystal. Each photodetector (Figure 2) consists of a $2\times 3$ array of $6\times 6$ mm${}^2$ SiPMs with 50 $\mathsf{\mu }$m pixels. The photodetectors are UV-extended to better match the 315 nm spectral peak of pure CsI.

2.3. FEE Boards

FEE boards installed directly behind each photosensor (Figure 2, top-left) perform signal amplification and shaping, while also managing all slow control functions. The six cells of each photosensor are connected in two series of three elements, to decrease their equivalent capacitance, and two grounded-base BJT stages form a fast current adder with ∼30 $\mathsf{\Omega }$ input impedance for each series (Figure 3). A second amplification stage drives a pole-zero cancellation network followed by a shaping section consisting in a 3-pole pulse stretcher based on a Bessel filter. Pulse stretching allows one to shape the pulse rising edge to approximately 25 ns to allow an optimal signal reconstruction. A balanced line driver transmits SiPM pulses to the DiRAC digitiser with a 2 V dynamic range. The overall system gain (1 or 2) can be set via a dedicated jumper.

The nominal 1 MHz rate capability of the analogue section allows one to cope with the high experimental hit rate.

A custom high-stability, low-ripple high-voltage linear regulator provides individual bias voltages up to 200 V—programmable in 50 mV increments via the on-board SPI DAC—with a nominal current capability of 2 mA per SiPM. Integrated bias, current and temperature monitors are also present, with dedicated 12-bit SPI ADC channels.

2.4. Mezzanine Boards

Mezzanine boards employing ARM microprocessors perform the power distribution (low- and high-voltage) to all ROUs and manage all slow control functions to the FEE via the dedicated SPI bus and multiplexer.

2.5. DiRAC

Each DiRAC board employs a PolarFire rad-hard FPGA and allows one to digitise 20 calorimeter channels at 200 Msps via 12-bit Flash ADCs. A VTRX-based 10 Gsps optical link allows a data stream and slow control management with the Mu2e trigger and data acquisition system (TDAQ) [5]. Full signal digitisation allows one to better cope with the high expected pile-up rates. More on the design and construction status of the DiRAC boards can be found elsewhere [6,7].

2.6. Signal Reconstruction

Timing and charge reconstruction, along with pile-up disentangling, are carried out on sampled pulses by means of a template fit algorithm, using a 3-parameter optimisation. To compensate any channel-to-channel signal variation, individual waveform templates—sets of interpolated nodes with fixed proportion—are generated for each channel by aligning and averaging a large dataset of sampled hits (Figure 4). Templates can then be fitted using a three-parameter optimisation (wave baseline, wave timing offset and scale). The timing fit parameter is used as a timing reference for the characterisation studies hereby reported.

3. Qualification and QC

3.1. SiPM

An intense QC campaign was carried out on production SiPMs, including measurements of V${}_{br}$, I${}_{dark}$ and PDE, and their relative spread over the six cells of the photosensor. It was verified that the SiPM’s gain and pulse shape remained unaltered after exposure to 100 krad and 10${}^{12}$ n${}_{1\mathrm{MeVeq}}$/cm${}^2$. The maximum allowed 2.5 mA dark current following irradiation and a mean time to failure MTTF > 10${}^7$ h at 20 °C were verified on all production batches. The quality control on production SiPMs ended in late 2019 with less than 2% out-of-spec components. The main results are summarised in Figure 5. Further details can be found elsewhere [8].

3.2. Electronics Qualification Steps

General qualification steps for the FEE, MB and DiRAC comprised tests carried out to validate performance with respect to the harsh Mu2e operating environment, with a particular attention to radiation hardness performance, as discussed later. Tests with B-fields were carried out at the Argonne National Laboratory [9] on the FEE, MB and DiRAC (Figure 6, top-right), with a focus on the buck converter operation and efficiency. Tests in vacuum were carried out on the MB, FEE and DiRAC. Quality control tests on production FEEs and MBs included burn-in tests at 65 °C, along with calibration of the HV regulator and temperature/current sensing circuitry (Figure 6). Signal shape and gain linearity were also assured by pulsing the analogue section with a dedicated signal injection setup.

3.3. Readout Units QC

An automated testing station was developed at LNF-INFN to perform QC steps on production ROUs [10]. As shown in Figure 7, a 420 nm nanosecond UV LED source was pulsed at kHz repetition rate over two assembled ROUs at a time. A motorised neutral density filter selector was used to vary the source intensity.

The SiPM gain was evaluated as $G=\sqrt{p_0}/\left(q_e\left[pC\right]\right)$ from the fit $\frac{{\sigma }_Q}{Q}=\frac{p_0}{\sqrt{Q}}\oplus \frac{p_1}{Q}\oplus p_2$, also shown in Figure 7. The SiPM response was calibrated—and stored in a dedicated ROU database—for different overvoltages and operating temperatures. An overall gain spread <3% RMS was observed in production. A reproducibility better than 1.5% was observed during the station operation.

4. Irradiation Tests of the Readout Electronics

4.1. Radiation Environment and Requirements

A full simulation of the Mu2e radiation environment was carried out via MARS15 and FLUKA. The electronic boards were mocked by equivalent boxes containing a compound of Si, Au, Al, Cu and polyethylene. Worst-case requirements for the EMC electronics (upstream disk, inner region) were set to 100 krad (inclusive of a safety factor of 12) and 3 $\times {10}^{12}$ n${}_{1\mathrm{MeVeq}}$/cm${}^2$ for (inclusive of a safety factor of 6) for 5 years of operation. The reference high-energy (>20 MeV) hadron flux for single-event effects (SEE) validation was set to 1 $\times {10}^{10}$ p/cm${}^2$.

4.2. Ionising Dose Campaign

All total ionising dose (TID) qualification steps for the FEE and MB were carried out at ENEA-Calliope, a pool-type, high-intensity ${}^{60}$Co source. After a 2-year-long campaign, the final choice of FEE rad-hard SPI ADCs and DACs (ADC1280S022CIMT and 121S101CIMKX, both from Texas Instruments) was confirmed with a final test in 2020, in which a total of six FEE boards from controlled production lots were irradiated up to 105 krad (0.85 krad/h dose-rate) in different operating configurations (Figure 8). All slow control monitoring channels exhibited excellent stability throughout the test with a maximum absolute drift of 2 LSB, comparable with the case of nonirradiated boards. The stability of the regulated high voltage was better than 0.1 % at a nominal 160 V. Similar tests with ionising dose on the MB were also performed to validate the hardness of the microprocessor and the SPI communication, along with the stability and efficiency of the onboard buck converters. No detrimental effects were observed for the 25 krad TID MB requirement.

4.3. Nonionising Dose Campaign

Tests were carried out at ENEA-FNG [11], utilising a 14 MeV neutron gun based on T(d,n)$\alpha$ fusion. Separate runs were carried out on the FEE and MB boards. No failure or performance degradation was observed in both cases due to cumulative nonionising dose damage after the delivery of an equivalent $3\times {10}^{12}$ n${}_{1\mathrm{MeVeq}}$/cm${}^2$ fluence over the course of 5 h.

4.4. Single-Event Effects

Dedicated studies for SEE have been carried out and are in progress on the MB processor and the DiRAC FPGA, to study firmware functionality, memory and communication corruption in the DiRAC–MB–FEE chain. The measurement of single-event upset (SEU) and latch-up (SEL) cross sections is crucial for the estimation of the expected functional interrupt rate under the high-energy hadrons exposures during experimental runs, and to optimise the implementation of hardware (power supply latch-up protection, watchdogs and dedicated reset lines) and firmware (scrubbing, memory and data integrity checks) countermeasures.

With respect to the MB processor, a preliminary simulation was carried out on-bench by fault-injection, accessing the processor memory via JTAG and performing random single or multiple bit flips while running the usual slow control functions, while first SEE data were taken during the aforementioned runs at ENEA-FNG (Figure 9, left) to validate the mitigation techniques effectiveness: the occurrence of approximately one SEU-related functional interrupt every $2\times {10}^{10}$ n${}_{1\mathrm{MeVeq}}$/cm${}^2$ delivered on the MB processor was observed, and it was verified that the board could correctly detect and recover operations via reset or power cycling. Similar tests were carried out for the FEE and DiRAC boards. The final validation tests for SEE on the FEE–MB–DiRAC chain are being carried out, at the moment of writing, using a 200 MeV proton beam at the Northwestern Medicine Warrenville Proton Center, US (Figure 9, right). These data will provide the correct estimation of the expected mean time between functional interrupts for each board, along with the validation of the SEL protection and other countermeasures.

5. Readout Chain Vertical Slice Test

In late 2021 a vertical slice test (VST) of the full calorimeter read out chain was performed with cosmic muons at LNF-INFN on the calorimeter Module-0 [3] (Figure 10, right), a large-scale prototype counting 51 crystals and designed with the final configuration of all detector services and subsystems (mechanics, cooling system, laser calibration system and vessel for vacuum operation).

For the VST, an electronics crate with a DiRAC and an MB board—and the respective cooling systems—were installed inside the Module-0 vacuum vessel to read out a total of 20 FEE channels (Figure 10, left). The 17 innermost crystals were instrumented with production ROUs, and a dual readout was used for the central and the central uppermost/lowermost crystals.

Triggering was performed either via the coincidence of two scintillator plates placed above and below the Module-0, or via the DiRAC onboard autotrigger, for which a minimum hit number of three crystals was set. The energy response was equalised online acting on SiPM biases to a level <5% and further offline equalisation was performed on the minimum ionizing particles (MIP) most probable 21 MeV deposit (Figure 11, left). The pedestals were characterised for all channels (Figure 12, right), yielding a ∼200 keV worst-case equivalent noise. The SiPM gain stability was checked against the MIP peak, and its fluctuations were compatible with the ±1 °C stability of the temperature regulation system (Figure 12, left).

The MIP tracks inside the calorimeter were reconstructed with a 2D fit in the plane orthogonal to the crystal axes (Figure 11, right) and an iterative algorithm was used to calculate an energy-weighted time centroid for each event and adjust the relative timing offset between channels, accounting for the relative 2D time-of-flight (ToF) corrections, resulting in a relative alignment better than 5 ps RMS after five iterations (Figure 13, left). A mean time resolution of order 210 ps was evaluated with MIPs for individual crystals, from the time difference between the left and right channels of the corresponding ROU, as shown for the central crystal in Figure 13, top-left.

The timing difference distributions between channels, after ToF correction, were well centred around zero (Figure 13, bottom-left), and their sigma scales well within the quadrature sum of the individual timing resolution terms previously obtained, apart from a minor spread, associated to the ToF contribution due the unknown track slope along the calorimeter axis. An algorithm to reconstruct this slope from the optical transport timing differences on the crystal matrix is currently being developed for increased performance. The timing results were similar for the other instrumented channels. Finally, the mean timing resolution dependence on energy over the MIP energy range is shown in Figure 14 for the central channel.

6. Summary

The Mu2e calorimeter has fully custom readout electronics, from SiPM to digitiser. The design of the FEE, mezzanine boards, and DiRAC digitiser boards were validated after an intense and prolonged qualification campaign, to assess their performance and compliance to the requirements under the Mu2e operating conditions, featuring a high vacuum, 1 T B-Field and a harsh radiation environment. The production of all FEE was completed in 2021, and the ROU assembly and QC are in progress. The final SEE tests are in progress for the MB and DiRAC boards, prior to their production launch.