The Muon 𝑔 -2 experiment E989 [1] at Fermi National Accelerator Laboratory (Fermilab) aims to determine the anomalous magnetic moment 𝑎 𝜇 ≡ (𝑔 -2)∕2 of the muon to a relative precision of 140
where the anomalous precession frequency is defined as 𝜔 𝑎 ≡ 𝜔 𝑠 -𝜔 𝑐 .
The experiment aims to measure 𝜔 𝑎 to a statistical precision of 100 ppb, with 𝜔 𝑎 and 𝐵 systematics determined to 70 ppb each. This paper describes the final prototype instrumentation developed to measure 𝜔 𝑎 . It is organized as follows. Section 2 outlines how the measurement of the anomalous precession frequency is carried out and the consequent technical demands on the instrumentation that guided the design decisions. Section 3 describes each of the key calorimeter subsystems. Section 4 describes the experimental setup at SLAC test beam. Section 5 outlines the reconstruction and calibration methods. Section 6 provides system performance metrics such as timing resolution, linearity, energy resolution, and spatial uniformity. The measurements and analysis described herein are based on a three-week run using the End Station Test Beam (ESTB) at SLAC in June 2016. The final production of the complete calorimeter system is now installed and running in the Muon 𝑔 -2 experiment at Fermilab.
The design of the Muon 𝑔 -2 experiment at Fermilab follows largely the well-known method employed most recently in the E821 experiment [2] at Brookhaven National Laboratory (BNL). Intense bunches of polarized positive muons having a central momentum of 3.1 GeV/c are injected into the BNL 7.11 m radius, 1.45 T superconducting storage ring [3] that was relocated to Fermilab. The beam has negligible hadron contamination, but a non-negligible fraction of positrons. Once injected, the particles receive a magnetic kick to deflect them onto a stable orbit. Only a few percent of the muons remain after a few turns (cyclotron period is 149.2 ns). The unstored muons, together with the positron beam contamination, generate a prompt background in the calorimeter systems that is anticipated in the detector and electronics designs. The storage ring fill cycle is repeated at an average rate of 11.4 Hz.
Stored muons orbit the ring and their spins precess according to Eq. (1). The 𝜔 𝑎 frequency, which is the subject of the instrumentation described here, is encoded naturally in the time and energy distribution of decay positrons. Because parity violation in 𝜇 + → 𝑒 + ν𝜇 𝜈 𝑒 associates the decay positron energy in the laboratory frame to the average muon spin direction at the time of decay, the higher-energy positrons are preferentially emitted when the muon spin is aligned with its momentum. Lower-energy positrons are emitted when the spin is reversed. Therefore the number 𝑁 of higher-energy positrons striking detectors follows a functional form 𝑁(𝑡) = 𝑁 0 exp(-𝑡∕𝛾𝜏)[1 + 𝐴 cos(𝜔 𝑎 𝑡 + 𝜙)]
(
where 𝑁 0 is a normalization, 𝛾𝜏 is the time-dilated muon lifetime (≈ 64.4 μs), 𝐴 is the decay asymmetry, and 𝜙 is an arbitrary phase.
In practice, 24 calorimeter stations are positioned evenly around the inner radius of the storage ring, each being tucked into a notch of the scalloped vacuum chamber system. Fig. 1 shows the locations of two calorimeters with respect to a segment of the storage ring. The stored muons are constrained to occupy a 9-cm diameter cross-sectional area within the vacuum chamber. The decay positrons have momenta below the muon momenta and therefore curl to the inside of the ring. They exit the vacuum chamber through a thin aluminum wall directly adjacent and parallel to the front face of the calorimeter.
We define the Calorimeter System to include the following subsystems: the physical calorimeter, its readout transducers, its electronics and controls infrastructure, the waveform digitizers, a clock and control distribution network, a laser gain calibration and monitoring system, and a GPU-based frontend readout coupled to a data acquisition environment. The system must record data in a dead-time free manner throughout the ≈ 700 μs muon fills after receiving an accelerator trigger. Each calorimeter station produces ≈ 800 MB/s of raw digitized samples, which must be filtered using an online GPU farm so that only regions, dubbed ''islands'', that have at least one over-threshold sample are saved for offline analysis.
The leading 𝜔 𝑎 systematic uncertainties associated with the calorimeter are multi-particle pileup and gain stability. Both are exacerbated by the instantaneous rate that can exceed 10 MHz just after injection, but then drops by a factor of ∼ 50,000 over the course of a fill. Because of muon beam cleaning and debunching requirements, the typical start-time of a fit for 𝜔 𝑎 is 30 μs after injection. A combination of lessons learned from the BNL experiment and various Geant4 [4] simulations establishes several performance parameters that must be realized. These, along with practical considerations, motivate the following list of system requirements:
• The energy resolution at 2 GeV should be better than 5%.
• The calorimeter gain must recover by 30 μs following the intense flash at injection. • The calorimeter gain must remain stable during the 30-700 μs measurement period. • The laser calibration system must be able to correct for residual gain instabilities to better than 4 × 10 -4 during the measurement period.
• The time resolution of a reconstructed shower should be better than 100 ps for positrons with energy greater than 1.8 GeV. • The time stability of the system must be better than 7 ps throughout any fill, to ensure less than a 10 ppb shift to 𝜔 𝑎 .
• The calorimeter must be able to resolve two electromagnetic showers with impact time separations greater than 5 ns with 100% efficiency. • The calorimeter must fulfill the geometrical space limitation by fitting into the scallop of the vacuum chamber. • The calorimeter, readout, and cabling that reside adjacent to a highly uniform magnetic field must function in the field and, importantly, must not perturb the field uniformity.
An individual Muon 𝑔 -2 calorimeter station consists of fifty-four 25 × 25 × 140 mm 3 SICCASfoot_0 PbF 2 Cherenkov crystals stacked in a 9 wide by 6 high array, as shown in Fig. 2 (left). Each crystal is wrapped in a single layer of matte black, non-reflective Tedlar™ on the four long sides. A 12 × 12 mm 2 Hamamatsu MPPC 2 (SiPM) is glued to the downstream face using an index-matching optical epoxy. 3 In Ref. [5], we documented performance tests associated with the choice of PbF 2 Cherenkov crystals as a moderator for the calorimeter, and we explored two crystal wrapping alternatives. The use of Tedlar™ results in a lower light yield, but a narrower pulse shape, compared to a reflective white wrapping. In Ref. [6], we detailed the development of the SiPM electronics and the bias supply selection considerations, which had to be optimized for the high-rate and short-pulse-shape demands of this experiment.
The bias supply of the SiPMs is provided by four BK Precision 9124 commercial programmable DC power supplies. Typical values of applied bias voltages 𝑉 ap are in the range of 66-68 V, being adjusted to set the overvoltage 4 𝑉 ov in the range of approximately 1 V above the breakdown voltage 𝑉 bd . The choice of bias voltage values optimized the SiPM photo-detection efficiency (PDE), minimized the dark count rate, and allowed the gain to match the dynamic range of the electronics such that two coincident 3 GeV electrons would not saturate the digitizer dynamic range.
The low voltage for the SiPM pre-amplifier circuit boards is provided by an OTE HY3003-3 DC power supply. The SiPM pre-amplifier boards are controlled using HDMI cables that connect each to a breakout board mounted within the box service compartment. The breakout board distributes the bias voltage levels and the communication information between SiPM boards and a BeagleBone computer. The housing around the crystal SiPM ends is cooled from the bottom by air fans and a duct-work corridor internal to the box. Signals from the SiPMs are connected to the waveform digitizers through custom Samtec ECDP cables. The positioning of the calorimeter relative to the ESTB beamline is shown in Fig. 2 (right).
A SiPM is a temperature-sensitive device, and lab tests show that its breakdown voltage changes with a temperature coefficient of about 70 mV/ • C. We found that at an operational over-voltage of 2.4 V, the SiPM gain changes by 2.5% per • C; note that this overvoltage is higher than we nominally use. To initially equalize the gains of the SiPMs and to monitor their drifts, a laser calibration system has been developed. It is described in some detail in Ref. [7], and its main features are briefly recalled here. This laser system also represents an important tool for debugging functionality of the calorimeters, their electronics, and the data acquisition system prior to the experimental run. The laser calibration system deployed at this test beam is shown in Fig. 3.
In the laser box, the light from a pulsed diode laser (405 nm, PicoQuant LDH-P-C-405M, pulse length < 1 ns) is divided using an optical beam splitter with the first path directing 80% of the intensity to the calorimeter station through a 25 m long silica optical fiber. Upon The SM is designed to measure the laser power directly after the laser head, minimizing the sensitivity to any perturbation occurring in any optical elements in the light path. It utilizes a sizable fraction of the laser light in order to reduce shot noise and reach the required statistical precision rapidly. A specialized electronics module (Monitoring Board) [8] was utilized to manage photo-detectors (two pin diodes -PIN1 and PIN2 and a PMT), signal processing, and data readout for the three SM channels. It amplifies and filters the signals, digitizes them, and stores the data in local FIFO buffers. The back-end section collects the data frames from the buffer channels, performs the eventbuilding, and provides them to a PC server for data storage and further processing. The three filtered signals, made fully Gaussian with 600 ns width, are distributed to custom-designed waveform digitizers.
The LM, on the other hand, is designed to monitor the stability of the light distribution chain from the laser head to the optical fiber bundle after the diffuser. In particular, it is sensitive to any variance in the coupling and the transmission of the laser into the silica fiber and to the splitting inside the diffuser. The fast PMT of the LM sees two comparable pulses separated by 250 ns. The first pulse comes from the SM; the second corresponds to light having traveled approximately 50 m from the laser to the calorimeter and back. In this way, the LM allows for the monitoring and correction of instabilities introduced by most of the light distribution components.
The signal processing of waveforms produced by the calorimeter must provide high fidelity determination of both the time of arrival and the energy of a decay positron. The electronics must also be as robust as possible against potential rate-dependent biases in positron selection. To meet these needs, we have produced a custom-designed five-channel waveform digitizer (WFD5) [9]. The WFD5 conforms to the Advanced Mezzanine Card (AMC) MicroTCA ® standard and functions with the AMC13 module [10] designed by Boston University (BU) for the CERN CMS experiment. The AMC13 is responsible for the distribution of the synchronous clock and control signals within the MicroTCA crate [11], for the readout of the AMCs within the MicroTCA crate, and for the transmission of the data to the DAQ front-ends, which will be described in the next subsection.
Fig. 4 summarizes the architecture of the clock distribution and waveform digitization systems used for the SLAC test beam run. Apart from the master reference clock, it is reflective of that employed in the Muon 𝑔 -2 experiment. Here, an Agilent N5183A MXG Microwave Analog Signal Generator was used to provide the 40-MHz clock.
The Clock and Controls Center (CCC) takes the master 40 MHz clock and a beam-correlated signal as input and encodes both using the synchronous Timing, Trigger, and Control (TTC) protocol [12] developed for the CERN LHC experiments. The TTC system uses the master clock to realize a 160 MBaud biphase mark encoder that timedivision multiplexes two channels using a balanced DC-free code [13]. The ''A'' channel provides a trigger to signal the start of data acquisition for a muon fill or for a laser calibration run. The ''B'' channel provides a variety of synchronous commands, such as clock counter resets or the type of trigger that the ''A'' channel will represent.
The TTC Encoder within the test beam's CCC system (Fig. 4) consists of an FC7 MicroTCA module [14] outfitted with an EDA-02708 FPGA mezzanine card (FMC) to accept inputs from the clock and the accelerator controls system. The FC7 optically transmits the encoded TTC signal to an AMC13 module in the same MicroTCA crate via a second FMC, a CERN EDA-02707. The AMC13 distributes the TTC signal over the MicroTCA backplane to fanout FC7 modules. These modules further distribute the TTC signal through a pair of EDA-02707 FMCs to all MicroTCA client crates housing detector electronics. The test beam had two client MicroTCA crates, which housed the WFD5s that instrumented the PbF 2 calorimeter and the laser calibration system.
The WFD5 modules for PbF 2 detectors were custom-designed for the Muon 𝑔 -2 experiment. Each WFD5 digitization channel utilizes a Texas Instruments ADS5401 12-bit, 800-mega-samples per second (MSPS) ADC chip for the digitization. Each channel also has a dedicated Xilinx Kintex-7 FPGA that acquires the data from the ADS5401 and buffers them in a 1-Megaword DDR3 16-bit SDRAM. Each channel FPGA manages the readback of its DDR3 buffer.
A sixth Xilinx Kintex-7 FPGA (the ''master'' FPGA) controls the operation of the WFD5 module and implements the external DAQ and communication links. The master FPGA communicates with each channel FPGA over dedicated 5-Gb serial links to provide channel control and to acquire channel data. The master FPGA also transfers data at 5 Gb/s over the MicroTCA backplane to the AMC13 with a BU-designed protocol based on an 8b/10b encoding. An event builder on the AMC13 gathers the data from all 11 WFD5 modules from a calorimeter and transmits these data to a dedicated DAQ front-end via a 10-Gb Ethernet link.
The 800-MHz clock for the WFD5 module derives from the 40-MHz master clock recovered from the TTC signal by the AMC13 and transmitted over the backplane. The WFD5 module provides significant flexibility in the digitization rate of each channel, with fieldconfigurable rates ranging from 40 to 800 MHz. A Texas Instruments LMK04906 clock synthesizer provides the frequency upconversion for the five channels, with a sixth output that is routed to the WFD5 front panel. For the test beam, all channels operated at the nominal 800 MSPS.
For the SLAC electron bunches, we specified a single sample window of 560,000 ADC samples (700 μs) in length, similar to that needed for the Muon 𝑔 -2 experiment.
The data acquisition system (DAQ) used was a small-scale copy of the main DAQ designed for the Muon 𝑔 -2 experiment. The system permitted the acquisition of dead-time-free, 700 μs duration, continuously digitized waveforms from the calorimeter crystals and from the laser monitors. It was triggered by both beam and special laser-calibration events at a total trigger rate of roughly 10 Hz. An event in this language corresponds to all of the hits (typically from thousands of positrons) that are accumulated during a 700 μs muon fill in the Muon 𝑔 -2 experiment, or a sequence of laser pulses which are fired between fills. The concept was tested here, even though there was only one electron per event on average. The DAQ comprised a front-end computer and a back-end computer. The front-end computer is responsible for the readout and preprocessing of the continuously digitized waveforms from the detector systems. On the other hand, the back-end computer is responsible for the event assembly, data storage, and run control. The front-end computer reads out the raw data from the three MicroTCA crates over three point-to-point 10 GbE fiber-optic links and pre-processes them into derived datasets using a hybrid system comprising the computer's eight-core processor and two general purpose graphical processing units (GPUs). An overview of the data flow at the test beam is shown in Fig. 5. The raw data rate from the test beam electronics was approximately 2 GB/s and the processed data rate to the mass storage was approximately 200 MB/s. In the Muon 𝑔 -2 experiment, the GPU-based pre-processing is necessary to reduce the enormous rate of continuously-digitized waveforms of approximately 20 GB/s to a manageable rate of approximately 200 MB/s of stored data.
The acquisition software is based on the MIDAS [15] data acquisition framework developed at PSI and TRIUMF. The front-end readout consists of: a TCP thread that receives and re-assembles the raw data from the AMC13 controller; a GPU thread that manages the GPU-based data processing into various derived datasets; and a MIDAS thread that handles the transfer of MIDAS-formatted events to the back-end computer event builder. Mutual exclusion (Mutex) locks are used to synchronize the execution of threads and ensure the integrity of data.
The GPUs were used to capture islands of ADC samples for all of the crystals in the array, regardless of their own sample values, when a sample from any of the 54 crystals exceeds a programmable threshold. The island length includes a programmable number of pre-and postsamples around the trigger sample and is automatically extended if a second software trigger occurs during the island sample-length time period. The GPU processing was implemented on two NVIDIA Tesla K40 GPUs using custom CUDA kernels to parallelize the data processing over the 5776 available CUDA cores.
The MIDAS tools for event building, data storage, and run control were all hosted on the back-end computer. MIDAS also provided an online database (ODB) used both for saving the experimental conditions for each run and for configuring the detectors, electronics, and other subsystems.
The experimental setup during the test beam is shown in Fig. 5. The calorimeter was located about 50 cm downstream of the beamline exit window. The electron beam size on the front face of the calorimeter, after exiting the vacuum pipe and traveling through air, was ≈ 0.5×1.0 cm 2 . The mean transverse position was stable throughout the measurement period. The calorimeter, waveform digitizers, bias control, low-voltage, and front-end computer were located inside the beam tunnel. The laser and monitoring system was located just outside of the tunnel, and the DAQ backend and analysis machines were located in the control room several floors above the tunnel. This setup is representative of the distributed setup now employed at the Fermilab Muon Campus.
The End Station A secondary beamline provides a well-collimated beam of electrons at a user-defined intensity, with a typical beam-pulse frequency of 5 Hz. In each precisely timed pulse, a Poisson-distributed number of electrons arrive. During our normal mode of operation, the single-electron beam event was optimized at a probability of 37%. Energies from 2.5 to 5 GeV were used in evaluating the performance of the calorimeter system, with each energy requiring a setup period by the operators. Once the beam energy was established, the rate was stable and the energy constant to better than 1%. Choosing new energies required some hours of setup and the actual energy was typically accurate to 10% of the target value. This limited the study of linearity by sweeping beam energies. With these considerations, it was prudent to fix 3 GeV as the nominal energy for the bulk of the data taking. It corresponds to the endpoint energy of the decay positron in the Muon 𝑔 -2 experiment.
Approximately 3 μs prior to the arrival of the beam, a machine trigger was received from the SLAC accelerator control system. This trigger was then fed into the FC7 to be encoded together with the 40 MHz clock. The FC7 then sent TTC signals to the AMC13 in its MicroTCA crate and also to the MicroTCA crates of calorimeter and laser systems, which initiated digitization of the WFD5s in both crates. The FC7 also sent a NIM signal to a Laser Control Board [16] which triggered the PicoQuant laser head in a pre-configured laser firing pattern following a 1 μs delay. An example pattern implemented during the test beam is shown in Fig. 6.
The offline data analysis was performed using the art -based event processing framework [17,18] developed at Fermilab. Data products consist mainly of fits to digitizer islands to extract pulse-integrals, fit times, and pedestals. The fit results require energy calibration to MeV units and time-dependent gain corrections.
The pulse-integral is an effective measure of the number of SiPM pixels fired. It is extracted along with the definition of the hit time from a fit to the digitized trace. The pulse-fitter is based on custom pulse templates 𝑇 (𝑡) for each individual SiPM, under the assumption that the method is robust against small fluctuations in pulse shape. The procedure for template creation was described in Ref. [5]. Separate templates were built for each of the 54 crystals and for both electron showers and laser events. A comparison of electron-beam and laserbeam templates for the same crystal is shown in the upper panel of Fig. 7, while a comparison of electron-beam templates for two different crystals is shown in the lower panel.
The function used to fit the traces is:
The three free parameters include an overall scale factor (𝑠), the peak time (𝑡 0 ), and the pedestal (𝑃 ). The pulse-integral is extracted as 𝑠.
The eigen [19] linear algebra library is utilized in the fitting process owing to its computing performance. This procedure results in a pulseprocessing rate of approximately 65,000 pulses per second per CPU, 5which exceeds the expected data rate for a single calorimeter in the Muon 𝑔 -2 experiment. As described in [6], the technique permits separation of pileup events with 100% probability for separation times larger than 5 ns.
Template fitting relies on the stability of the SiPM pulse shape against various effects such as the number of photons (or energy), impact position, and impact angle of the beam. The results informed the reconstruction algorithm design. The laser calibration system can deliver light pulses with controlled relative intensity through a rotation of its neutral density (ND) filter wheel settings. A series of template pulse shapes for laser events that vary in intensity from 20-100% is shown in the upper panel of Fig. 8. There is no noticeable difference versus the number of photons over the explored range.
The template pulse shape for different electron impact angles is a critical test for Muon 𝑔 -2 experiment, because the positron events curl into the calorimeter at energy dependent angles. As shown in the middle panel of Fig. 8, there is a negligible difference between the 0 • and 15 • pulse shapes. More noticeable is the slightly changing pulse shape versus impact position of the electron beam on the crystal face. As shown in the bottom panel of Fig. 8, when the electron beam is shifted from crystal 24 to its neighbor crystal 25, the region from 5-10 samples becomes somewhat broader. The impact of this effect on the energy scale and resolution is made negligible by building templates integrating over all hit positions each crystal records.
To establish a uniform response across the calorimeter, the SiPM gains must first be matched as closely as possible. Their response depends mainly on the bias voltage and the ambient temperature. At the assembly level, the SiPMs were selected from a larger sample based on similar breakdown voltages according to the provided data sheets from Hamamatsu.
The operating temperature of the SiPMs within the calorimeter housing depends on their physical placement in the calorimeter. A relatively continuous ≈ 5 • C temperature variation exists from the leftmost column to the rightmost column owing to the airflow inlet and outlet. The SiPMs are (approximately) arranged into groups of calorimeter columns that are then connected to one of the four independent bias supplies, which provides an average 𝑉 ov appropriate to the group. Each SiPM amplifier board contains a programmable gain amplifier (PGA) with amplification adjustable over the range 2-20× (in 80 steps), which then provides for a final adjustment.
Photon calibration is the step where the fitted pulse area can be interpreted as the number of photo-electrons (npe), or pixels, fired for a given pulse. The calibration constants are obtained through a procedure using the laser and a series of runs with varying neutral density filter settings as described in [5][6][7]. An iterative sequence was used alternating between laser-calibration runs and adjustments to the PGA gain settings of the SiPMs and the SiPM group bias voltages. After several iterations, equalization of the calibration constant at the level of 8% in RMS was achieved as shown in Fig. 9. Further improvement is not necessary, as the energy scale of each calorimeter channel does not depend solely on the SiPM gain but also on the SiPM photo-detection efficiency, the PbF 2 light yield, and the PbF 2 light transmission. An energy-scale equalization and calibration technique utilizing the electron beam is described in the next subsection. Fig. 8. Upper: Peak-aligned templates created from different neutral density filter wheels corresponding to different photon yields. Middle: Fit templates created from different beam impact angles. Lower: Fit templates for various beam impact positions when the beam is moved from the center of a crystal to its direct neighbor in increments of 0.2 times a crystal width.
In order to extract the energy calibration constants that convert npe to GeV, the 3 GeV electron beam was aimed in turn at the center of each crystal. A single crystal contains approximately 85% of the shower energy, or 2.55 GeV. The number of photo-electrons fired for each of the crystals is shown in Fig. 10. The average of ≈ 2100 npe corresponds to an energy calibration constant of ≈ 0.82 npe/MeV, consistent with the value reported in [5]. The accuracy of this technique is somewhat limited by the 0.5 × 1.0 cm 2 beam spot and the precision with which the calorimeter could be aligned to the center of the beam.
During the test beam run, laser calibration filter wheel scans were performed on average every three hours to establish the stability of the SiPM gain (conversion constant of npe to pulse integral). To correct for the SiPM gain drift between two consecutive filter wheel scans, which is expected owing to environmental temperature fluctuations, a different procedure was employed. A hybrid fill structure was created as shown in Fig. 6. The light pulses that strike a SiPM are, in time order: the laser sync pulse, the Cherenkov light from an electromagnetic shower, and a series of approximately 70 in-fill laser pulses. The gain drift correction is made with respect to a reference run, where the average laser pulse area for each crystal is 𝐿 0 . The average of the SiPM response to the series of in-fill laser pulses for each fill 𝑖 is ⟨𝐿 𝑖 ⟩. A gain correction factor 𝐶 SG = 𝐿 0 ∕⟨𝐿 𝑖 ⟩ is established for each fill.
To correct for intrinsic fluctuations of the laser, information from source monitors is used similarly. The average response of PIN1 and PIN2 during the reference run is 𝑆 0 , and the average of the in-fill laser pulses measured by source PIN1 and PIN2 is ⟨𝑆 𝑖 ⟩. The laser fluctuation correction factor 𝐶 LF = ⟨𝑆 𝑖 ⟩∕𝑆 0 is then applied to the gain corrected crystal hit.
The resultant energy stability is shown in the left panel of Fig. 11, while the temperature of the corresponding period is shown in the right panel. The drop in temperature coincides with an increase in the energy scale. In the actual Muon 𝑔 -2 experiment, a laser pulse sequence is fired between normal muon storage ring fills.
As multiple WFD5 modules and multiple channels in a WFD5 were used, it was necessary to align these channels in time to enable accurate hit clustering and to optimize the rejection of multi-particle pileup. To align all the digitizer channels within a MicroTCA crate, the laser sync pulse, fired after each trigger, was used. The pulse is distributed to all 54 crystals at the beginning of each event. For example, the time difference distribution of the laser sync pulse between channel 12 and 44 is shown in Fig. 12. Two narrow distributions exist, and both of them are not centered around zero. The former is owing to the fact that as the master TTC clock is running at 40 MHz while the digitizers are running at 800 MHz, there can be small differences in the exact digitizer clock tick on which a given ADC will begin digitizing. Because the ADS5401 ADCs always send data from the ''odd internal ADC'' first, we see a clear separation of 2.5 ns or 2 clock ticks (1 clock tick = 1.25 ns) as a result. The non-zero mean value is due to the length difference in light distribution fibers. The width of the distribution will be characterized in the timing resolution section.
The hit clustering algorithm used at SLAC is based solely on time partitioning, whereby all hits on a given island within a time separation of 𝛥𝑇 are grouped. The parameter 𝛥𝑇 varied depending on the analysis type. The cluster time was given by the hit time of the highest-energy crystal within a cluster. The electron's incident position is reconstructed using a center-ofgravity method with logarithmic weights [20]. This weighting accounts for the exponential falloff of the energy deposited in each crystal across the calorimeter. It is calculated by
where 𝑥 𝑖 and 𝑦 𝑖 are the coordinates of the crystal column or row 𝑖, respectively, and 𝑤 𝑖 is its logarithmic weight given by
with 𝑤 0 being a free parameter that sets the relative importance of shower tails in the weighting. As shown in Fig. 13, the optimal value of 𝑤 0 = 3.5 was determined by minimizing the position resolution, which agrees with Geant4 simulations. The electron beam size, shown in Fig. 14, is found to be approximately twice as wide vertically as horizontally. It was verified with an ePix 20 × 20 mm silicon detector installed between the beamline and the calorimeter.
This section describes the performance of the calorimeter system as a whole in response to the electron beam. In all cases, the times or energies represent those of the shower clusters. In general, the findings using our full calorimeter with black-wrapped crystals, glued-on SiPMs, optimized SiPM pulse-summing circuits, and the 𝑔 -2 custom digitizers, are equal to or better than those published earlier using prototype devices [5].
The timing measurement of the Muon 𝑔 -2 calorimeters is driven by three factors: (a) the transit time of the SiPM device, the electronics amplification circuit, and the digitizer electronics; (b) the statistical fluctuation of photons triggering SiPM pixels; and (c) the electromagnetic shower profile fluctuation, the Cherenkov light emission, and the light propagation time. Quantitative analyses of (a) and (b) will be given, along with typical results for (c).
The timing resolution from (a) and (b) can be extracted using laser events that arrive at all SiPMs at about the same time. The time difference distribution between two WFD5 channels of the same laser event is shown in Fig. 15. The solid line is a Gaussian fit to the distribution, and the standard deviation is 27 ps for this particular set of WFD5 channels. The energy dependence of this timing resolution is studied for two different pairs: crystals read out by the same, or different, WFD5s.
The standard deviation of the time difference distribution is plotted against energy in Fig. 16. The difference in the timing resolution between ''different WFD5'' and ''same WFD5'' is attributed to the timing jitter of the 800-MHz clock between two different WFD5s. The data points from each pair are fitted using
where 𝐶 T is the constant term, 𝑆 T is the stochastic term, 𝜎 n is the noise term, and 𝐸 ef f is the effective energy term defined by
, where 𝐸 1 and 𝐸 2 are measured in the unit of npe. The use of 𝐸 ef f is motivated by the fact that there is a variation in laser intensity among calorimeter channels. The value of 𝜎 𝑛 is fixed at 5 npe based on digitizer noise level. The 𝐶 T term, which is mainly due to the timing jitter, is negligible for timing difference between two channels of the same WFD5. For two channels read out by different WFD5s, 𝐶 T is of the order of 25 ps. On the contrary, the 𝑆 T term, which is related to the photo-statistics, is consistent between the two groups having an average value of approximately 600 ps.
To extract the timing resolution of an electromagnetic shower, we studied the time difference between two calorimeter channels for electron events. The distribution of the time difference is shown in Fig. 17 for two channels from the same WFD5 and different WFD5s. The single channel timing resolution for a 3 GeV electron event is 57∕ √ 2 ≈ 40 ps (dotted line). A systematic and precise way of quantifying this resolution is difficult due to the broad beam size and uncertainty in the beam hit position.
Fig. 18 shows the energy distribution when a 3 GeV beam is centered at a particular crystal. The solid line is a Gaussian fit within a ±2𝜎 range around the mean value. The fitted energy resolution here is 2.7%.
To extract the energy resolution and linearity functions, beam energies in the range 2.5 to 5 GeV were used. The reconstructed energy distribution for each beam energy was fitted to a Gaussian function to determine the mean and width. The linearity result is shown in Fig. 19. The range extends well beyond the decay-energy endpoint for the Muon 𝑔 -2 experiment, which is at about 3.1 GeV. where 𝐶 E is a term accounting for the shower containment variance as studied in [5] and 𝑆 E is a stochastic term describing the statistical fluctuation of the electromagnetic shower. The resultant fit is shown in Fig. 20. The extracted value of 𝐶 E is (1.86 ± 0.43)% and the corresponding value for 𝑆 E is (3.56 ± 0.77)%. Based on this parameterization, the energy resolution at 2 GeV is 3.1%, surpassing the design target of 5%. While the constant term is consistent with that found in [5], the resolution is improved by about 11%.
An important aspect of the calorimeter is the spatial uniformity of the response. This aspect was studied by comparing the energy response of the calorimeter for various beam impact positions. Fig. 21 shows the deviation from the mean reconstructed cluster energy in percent for each channel when the electron beam is impinging on the center of each channel. The homogeneity is better than 2% for noncorner channels. Lower reconstructed energies for corner channels, as well as the uppermost and lowermost rows and leftmost and rightmost columns, are expected due to the leakage of electromagnetic showers. A finer scan in position is also performed by moving the beam across three channels with a step size of 2-4 mm as shown in Fig. 22. The reconstructed cluster energy as a function of beam position is stable at the 2% level.
We report on a performance study of the instrumentation for measuring the anomalous precession frequency in the Muon 𝑔 -2 experiment at Fermilab. Unique features of this performance study include fast Cherenkov crystals of PbF 2 ; exclusive use of custom 800 MSPS waveform digitizers to digitize the SiPM response; the use of hybrid CPU-GPU and MIDAS DAQ system; the use of art -based event reconstruction framework; and a laser calibration system with a high degree of pulse-to-pulse intensity stability. Through a series of test beams at SLAC and a final validation presented in this paper, we froze the design of the calorimeter system and went into production for all subsystems. The ''in-ring'' performance of the 24-times-larger calorimeter system will be reported in a future publication.
Shanghai SICCAS High Technology Corporation, 1295 Dingxi Rd., Shanghai
200050, China.2 Multi-Pixel Photon Counter Model number S12642-4040PA-50.3 Zeiss Resin OK 2030 and Hardener H 950.4 Overvoltage, 𝑉 ov = 𝑉 ap -𝑉 bd , where 𝑉 ap is the applied voltage and 𝑉 bd is the breakdown voltage.
Tested on a Macbook Pro 2015 with a quad-core processor of 2.5 GHz Intel Core i7.