The muon beam monitor for the FAMU experiment: design, simulation, test, and operation

Abstract

FAMU is an INFN-led muonic atom physics experiment based at the RIKEN-RAL muon facility at the ISIS Neutron and Muon Source (United Kingdom). The aim of FAMU is to measure the hyperfine splitting in muonic hydrogen to determine the value of the proton Zemach radius with an accuracy better than 1%. The experiment has a scintillating-fibre hodoscope for beam monitoring and data normalisation. In order to carry out muon flux estimation, low-rate measurements were performed to extract the single-muon average deposited charge. Then, detector simulation in Geant4 and FLUKA allowed a thorough understanding of the single-muon response function, which is crucial for determining the muon flux. This work presents the design features of the FAMU beam monitor, along with the simulation and absolute calibration measurements in order to enable flux determination and enable data normalisation.

1 Introduction

The aim of the FAMU experiment (Pizzolotto et al. (2020); Vacchi et al. (2023); Rossini et al. (2024a)) is to explore the magnetic structure of the proton through a measurement of the proton Zemach radius (Carlson (2015); Antognini et al. (2022)). The latter is extracted from a measurement of the hyperfine splitting energy of the muonic hydrogen (H) ground state. H atoms are produced by injecting a high-rate, low-momentum pulsed muon beam into a pressurised gaseous target. The experiment is currently in operation at the RIKEN-RAL muon facility (Matsuzaki et al. (2001); Hillier et al. (2019)) at the ISIS Neutron and Muon Source (Didcot, UK). The observable of the experiment is the number of delayed muonic oxygen (O) X-rays resulting from the transfer of the muon from H to oxygen atoms. This is clearly dependent on the number of H atoms created, which is directly related to the incoming muon flux. As a consequence, having an accurate and efficient beam monitor with minimal beam absorption is a crucial point in the data normalisation.

A beam hodoscope, composed of two planes of 32 scintillating fibres read out by silicon photomultipliers (SiPMs), has been set up for the experiment. The specific design of this detector, discussed in Section 2, is the best match among the number of available channels (64), the detector area, and its thickness. Other detector designs, such as muon cameras (Lord et al. (2011)), were avoided in order to minimise the amount of material immersed in the beam, as the beam monitor is expected to stay in the beam for the duration of the experiment. Similar detectors for higher rates and continuous beams are being developed at other muon facilities, such as PSI (Papa et al. (2015); Papa et al. (2019); Dal Maso et al. (2023)). The hodoscope serves both as a beam shape detector to optimise beam centring and focussing and as a flux meter. The latter role of the detector is made possible thanks to the analyses reported in this work.

The estimated average negative muon flux with a momentum of 55 MeV/c is on the order of muons per second (Matsuzaki et al. (2001); Hillier et al. (2019)). The beam is delivered in two 70-ns spills with an average repetition rate of 40 Hz (the synchrotron rate is 50 Hz, but one in five pulses is directed to the other target station). Therefore, during a spill, approximately 100 muons are delivered in 70 ns. Even though the system is based on a SiPM readout with fast signals ( ns), it is clearly not possible to tell single-particle signals apart. For this reason, the detector measures the total deposited charge , which is converted into muon flux using the result coming from calibration measurements (Carbone et al. (2015); Bonesini et al. (2017); Bonesini et al. (2019); Rossini et al. (2023b); Rossini et al. (2023a); Rossini et al. (2024b)).

Initially, data from cosmic muon calibrations combined with Particle Data Group (PDG) results to match the gap between energies of approximately 4 GeV and the used beam momentum ( MeV/c), were used to obtain an estimate of the muon flux vs muon beam momentum (Bonesini et al. (2017)), which compared well with previously published results (Matsuzaki et al. (2001)). In this case, two 3-mm pitch hodoscopes (Hodo-2 and Hodo-3) were used.

Then, a 1-mm pitch hodoscope with adjacent fibres (Hodo-1) was calibrated at the CNAO synchrotron in Pavia (Italy) with a low-rate proton beam with energy loss comparable to FAMU muons (Rossini et al. (2024b)). A proton beam of kinetic energy 150 MeV was tuned to allow single-particle events and directed against the hodoscope for testing.

The latest FAMU hodoscope (Hodo-4), that is, the position-sensitive muon beam monitor detector, is here thoroughly discussed. The design of the detector, discussed in Section 2, fits best with the number of available channels (64) and the required thickness, active area, and space resolution. It is composed of 1 mm fibres, spaced by 1 mm. The simulation and tests of the detector are presented in this work. Calibration measurements were carried out in the FAMU setup, exploiting the RIKEN-RAL Port1 muon beam with a modified configuration to obtain a low-rate muon beam, as later explained in Section 3. This has been crucial in order to disclose single-particle signals. In addition, the detector has been simulated in Geant4 (Agostinelli et al. (2003)) and Flair-FLUKA (Battistoni et al. (2015); Ahdida et al. (2022); Vlachoudis (2009)) to understand its response and extract crucial parameters and information for its calibration.

The equation to extract the muon rate from the hodoscope reading is the following:where Hz is the beam repetition rate, is the total deposited charge during a full-rate beam spill, is the average charge deposited by a 55 MeV/c muon interacting with both planes of fibres, for muons interacting with one fibre only, and are the fractions of muons interacting with one or two fibres, respectively. The charge deposited in one fibre is written as because this value is not directly measured, and is calculated from the simulation. In particular, is extracted from the simulation, while is determined using low-rate data. The main aim of this work, that is, the calibration of the FAMU beam monitor, is to compute the value of the calibration factor . This work provides a general method to calibrate other fibre-based hodoscopes that will be used as charged particle beam monitors.

Simulation, measurements, and analysis techniques are presented in Sections 3, 4 and 5, whereas the results are extensively presented and discussed in Section 6. Eventually, the value of is computed, and a test estimation of flux during full-rate beam is shown.

2 Design of the 1-mm hodoscope (Hodo-4)

A channel (XY configuration) beam monitor has been set up for the FAMU experiment by INFN Milano-Bicocca and INFN Pavia. The hodoscope consists of two crossed planes of 32 single-clad Saint-Gobain/Luxium BCF-12 polystyrene scintillating fibres. Each fibre is squared, with a pitch of 1 mm, and each fibre is covered with a nominal m-thick layer of TiO₂-based extramural absorber (EMA, or coating) to avoid inter-fibre optical cross-talk.

Previous hodoscope versions (Carbone et al. (2015); Bonesini et al. (2019); Rossini et al. (2023b); Rossini et al. (2023a); Rossini et al. (2024b)) had either too much material immersed in the beam ( fibres with a 3 mm pitch, that is, active area of cm² and a total thickness of 6 mm) or too small active area ( fibres with 1-mm pitch, that is, active area of cm² and a total thickness of 2 mm). The key point of the detector described in this work is having a spacing of 1 mm between adjacent fibres, as shown in Figure 1. This allows a cm² detector area despite keeping the thickness below 2 mm. However, this hodoscope has inhomogeneous volume, which slightly complicates its response function to the muon beam, as further discussed in Sections 5 and 4. In fact, each muon can interact with zero, one, or two fibres depending on its interaction position in the XY plane. The geometric features of this model compared to the previous ones are presented in Table 1.

FIGURE 1

This model of the hodoscope has been considered the best compromise among all features, and it has, therefore, been installed in the final FAMU setup for 2023 and 2024 runs. Consequently, carrying out single-particle calibration was crucial for its operation in the FAMU experiment. The position of the hodoscope in the FAMU setup is shown in Figure 2.

FIGURE 2

Rossini et al. (2024b).

Each fibre is read out by a mm² Hamamatsu S12571-050P SiPM (cell size 50 m) on one side. SiPMs are supplied with a positive bias of V (finely tuned for each group of 16 SiPMs to optimise the uniformity of the detector response)¹. SiPM signals are fanned out through 4 m cables with MCX connectors and digitised using two CAEN V1742 (32 channels each). The trigger is supplied to all digitisers (including these) through the FAMU system, which provides a beam trigger coming from the synchrotron with a rate of 50 Hz.

3 Hodoscope measurements at RIKEN-RAL Port1

The FAMU experiment is installed at the RIKEN-RAL Port1 muon beamline at the ISIS Neutron and Muon Source in Didcot, United Kingdom. The experiment consists of a pressurised cryostat holding bar of a hydrogen-oxygen mixture at a temperature of approximately 90 K. The gas chamber is the target of the muon beam, with the aim of forming muonic hydrogen atoms.

The ISIS synchrotron accelerates protons with an energy of 800 MeV with a pulse rate of 50 Hz. Four consecutive pulses are sent to Target Station 1 (TS1), and one is sent to Target Station 2 (TS2). The graphite target for muon beamlines is located in the beampipe connecting the synchrotron to TS1. Hence, it receives protons and produces muon pulses at the same rate as TS1 ( pulses per second). Negative pions are directed in the RIKEN beampipe, where they decay to negative muons and are delivered to the four RIKEN-RAL Ports alternatively. RIKEN-RAL Port1 is currently dedicated to the FAMU experiment. The beam time structure, which is the same for high and low-rate measurements, is shown in Figure 3 as measured by the FAMU hodoscope in a low-rate test beam. Each beam pulse consists of two 70 ns spills separated by 320 ns.

FIGURE 3

3.1 High-rate measurements at RIKEN-RAL

The beam is generally set to work at the highest available rate. At the momentum value used in FAMU, the average rate is the order of muons/s. The quadrupole and bending magnet configuration have been optimised to deliver the best beam rate and geometry that suit the experiment during the FAMU beam commissioning in July 2023 (4 beam days, dataset RAL202303).

After that, two FAMU data-taking runs were carried out in October 2023 (6 beam days, RAL202305) and December 2023 (12 beam days, RAL202306).

3.2 Low-rate measurements at RIKEN-RAL

After carrying out some tests with protons at the CNAO synchrotron in Pavia (Rossini et al. (2024b)), it was decided to characterise the current hodoscope directly on the FAMU setup as its particular design was expected to be more sensitive to beam geometry changes.

In order to obtain a single-particle beam, the currents of some quadrupole and bending magnets were de-tuned to minimise the amount of pions directed in the beampipe. In particular, the first two quadrupoles (RQ1 and RQ2) were turned off to widen the pion bunch, and the first bending magnet (RB1) was de-tuned to direct the beam halo, rather than its central part, into the beampipe. See Figure 4 for a detailed map of the path followed by the beam from the target to Port1, including the magnets encountered. The choice of which magnets had to be tuned was made to not compromise the beam optics, which would result in not delivering the beam to Port1. The shut-down of the two quadrupoles resulted in a 90% beam intensity drop. The optimisation of the bending magnet current was carried out progressively to make sure that the rate would be as low as required. Figure 5 shows the effect of the progressive variation of the RB1 current out of its optimal value. The muon current (proportional to ) decreases, whereas the number of events marked as single coincidences increases, reaches a maximum, and starts to decrease. This latter behaviour means that the muon flux is so low it allows events with only one coincidence, that is, single-muon spills.

FIGURE 4

FIGURE 5

The low-rate data acquisition RAL202306 consisted of several hours of beam optimisation and final -event measurements with and without the beam.

4 Data analysis

The data analysis technique, described in Figure 6, is based on imposing single coincidences between each detector plane.

FIGURE 6

For each beam trigger, the FAMU DAQ system opens an acquisition window and digitises the signals coming from each hodoscope channel at a rate of 1 G/s. The resulting data packet is called an event. In the previous analysis procedure, extensively described and tested by Rossini et al. (2024b), an event was considered a single-particle hit if only one fibre per plane had integrated charge over a given threshold. This method has been considered valid as the hodoscopes had no interspacing, and the measurements were particle-triggered. This means that most muons arrived at the same time (given the fixed pre-trigger window) and passed by one fibre per plane. However, in this case, the trigger comes from the synchrotron, and the full-rate beam shape is complex (every synchrotron trigger corresponds to two ns spills separated by a ns gap). The low-rate time structure of the beam is the same: even though the beam is tuned to allow single-particle events, they might come from either the first or the second spill, and some events might have more than one muon. In addition, the spacing between adjacent fibres makes it less probable to have muons hitting one fibre per plane, as better discussed and quantified in Section 5. For this reason, it was decided to use a time coincidence-based approach.

During data processing, for each event and for each hodoscope fibre , the hodoscope low-rate data processing system retrieves the 64 waveforms and looks for peaks, returning the total integrated charge and, for every peak , the time-of-arrival and the pulse height . At this point, the coincidence is imposed, with a tolerance of 50 ns (small enough to exclude particles coming from two different spills), for hodoscope peaks having over a certain threshold to be determined. Events with only one coincidence are selected as single-particle events and used for the hodoscope characterisation.

The value of total deposited energy for every event is the sum on all fibres of the integrated charge . This holds for both low- and high-rate measurements. As one can see in Section 6, the shape of the histogram for low-rate measurements is asymmetric. After exploring some possibilities (combinations of Landau and Gauss profiles), it has been decided to fit this histogram with the convolution of a Gaussian and a decreasing exponential profile; that is:having five free parameters: additive constant , amplitude , Gaussian mean , Gaussian sigma and exponential decay constant . The fit boundaries are chosen by looking for optimal and stable reduced . The maximum, which corresponds to the estimate for , has no known analytical expression. Consequently, it was determined on the fit function through the numerical Brent method². The uncertainty on is obtained by variations of the fit boundaries around the optimum. This is done recursively to select a region in the two fitting boundaries in which . The variation of in this region is then used as an estimation for its uncertainty.

The value obtained by analysing the data taken at RIKEN-RAL with low-rate muons is ADC channels. The uncertainty is taken by varying the fitting boundaries and imposing . The histogram and fit are shown in Figure 7. This is the most probable value of deposited charge by 55 MeV/c muons interacting with two scintillating fibres, one per plane. As discussed, the muon flux with this geometry cannot be simply obtained as , as most muons do not interact with two fibres. However, the fraction of muons interacting with one fibre per plane and with one fibre only is mostly geometric and must be extracted from the simulation. Therefore, the flux can be estimated with Eq. 1, that is, weighting by a factor , which is obtained from the simulation in Section 5.

FIGURE 7

5 Hodoscope simulation

In order to understand the energy loss of muons in the detector and, therefore, its theoretical response function, the hodoscope has been simulated using the Geant4 (Agostinelli et al. (2003)) toolkit. The geometry consists of the fibres, coatings and entrance windows as described in Section 2; that is, each fibre (polystyrene, 1 mm pitch, 6.4 cm length) is coated with a 15 m layer of TiO₂ and positioned in a 32-fibre plane with 1 mm interspacing between adjacent fibres (measured coating-to-coating); two planes are juxtaposed with crossing fibre direction, separated from the world volume with a 0.1-mm-thick PVC window.

The muon beam simulated for this work is a 55 MeV/c negative muon beam with 2-dimensional Gaussian shape, with mm and mm. The reproduces the beam configuration optimised for the experiment, as measured with the hodoscope during 12 h of full-rate data acquisition.

In order to obtain an uncertainty budget, all simulations were repeated with Gaussian dispersion of momentum (%), variable beam size within the - uncertainties, and variable coating thickness, considering a 5 m coating thickness tolerance. This resulted in a geometric systematic uncertainty, which was added to the uncertainty budget as an independent contribution.

All primaries have been tracked and assigned flags depending on whether they passed by front and back plane fibres. In fact, given the geometry of the hodoscope, muons can pass by zero, one, or two fibres. The contributions of these muons are plotted separately and also jointly in Figure 8. In the zero-hit case, the energy deposit different from zero in some events is caused by secondary particles (e.g., delta rays generated in the coating and decay electrons) interacting with the fibres. The probability of hitting zero, one, or two fibres is about 25%, 50%, and 25%, respectively, as one can derive from geometrical considerations from Figure 9. However, the exact values of and depend on beam geometry, scattering processes, and coating thickness; as a consequence, they must be extracted from the simulation. is equal to the fraction of particles interacting with only one fibre, while is the fraction of muons passing by two fibres, one per plane.

FIGURE 8

FIGURE 9

In the Geant4 simulation, negative muons with momentum 55 MeV/c were launched. The beam geometry is the one extracted from 12 h of full-rate hodoscope measurements. As expected, about of the muons were marked as passing by two fibres, with a statistical counting uncertainty of around 0.2%. The uncertainty budget was completed by repeating the simulation with variations in the fibre pitch within its tolerance (30 m), the coating thickness within 5 m, the beam momentum within 10%, and the beam shape within the measured uncertainty. Other effects, such as small fibre misplacements, are expected to be averaged and cancelled due to the beam spot size. The total contribution, which is dominated by the uncertainty on the coating thickness, is about 1.8% on the number of muons passing by two fibres. The final estimate for the double-hit fraction is . This value is consistent with the heuristic value of 25% estimated from the geometry of the detector elementary cells. Similarly, the single-hit fraction estimate is .

In addition, the ratio between the double-hit and the single-hit mean deposited energy had to be computed. In fact, the data selection described in Section 4 enables only the inclusion of double-hit events. As a consequence, the value of calculated from the measurements is only the mean charge deposited by particles hitting two fibres. In principle, with a local linear approximation of the energy-loss curve, one could assume that the energy released by single-hitting muons is (i.e., that , but this must be verified, as the linear approximation might not hold. To do so, the single- and double-hit simulated spectra (see Figure 8) were fitted with Eq. 2. The fit stability was tested and used to determine the uncertainties as in the case of data coming from low-rate measurements (see Section 4), along with parameter variation. The estimated values of deposited energy are MeV for the 2-hit and MeV for the 1-hit. The resulting beam momentum straggling, which comes from the sparse detector geometry, is comparable to the momentum bite of the incoming beam (% (Matsuzaki et al. (2001))) and small compared to the stopping range in the FAMU apparatus. As a consequence, the presence of the hodoscope does not spoil the FAMU data.

The estimate for the ratio between the double-hit and the single-hit mean deposited energy is .

In parallel, an independent simulation based on the FLUKA-CERN (Battistoni et al. (2015); Ahdida et al. (2022)) toolkit was also developed using the Flair interface (Vlachoudis (2009)) for comparison (reported just as FLUKA in this work, for simplicity). The FLUKA simulation has been modelled to match the exact geometry and beam characteristics with the one in Geant4. The geometric factors determined with the FLUKA simulation with the same number of events are: and . The uncertainties on the FLUKA predictions are underestimated as they only comprise the statistical component. Both values differ by less than three standard deviations from the values estimated in Geant4. The qualitative comparison between the FLUKA and Geant4 histograms of the energy deposited in the hodoscope active volumes is shown in Figure 10. The deposited energies in the 1-hit and 2-hit cases, obtained from fitting with Eq. 2, are MeV and MeV. These energy deposits differ by about 7% and 6%, respectively, from the Geant4 values. Such a difference is generally considered a sign of accordance between the results retrieved from two independent codes. Their ratio is , which is consistent with the value extracted from the Geant4 simulation.

FIGURE 10

The FLUKA simulation was also used to estimate effective backscattering, that is, the fraction of particles interacting with one fibre in the first plane, one in the second plane, and then one back in the first plane, as a result of backscattering. Ideally, this value should be minimal for better hodoscope accuracy. The effective backscattering rate is ‰, which is negligible.

All the useful parameters extracted from Geant4 and FLUKA simulations are reported in Table 2. The coefficients for Eq. 1 used in this work for flux estimation are taken from the Geant4 simulation because the applied physics lists and transport thresholds had already been tuned and validated for the FAMU experiment.

Figure 11 shows some 2D distributions obtained using the Flair interface to FLUKA-CERN. In particular, the first row shows the probability density for the particles present in the simulation (primary muons, decay electrons, and photons/electrons resulting from elastic and inelastic processes). The second row shows the space distribution of the energy deposition, which is higher in the detector fibres than in air, reflecting the Gaussian shape of the beam on the plane. Muon inelastic processes are negligible for the sake of the experiment. The main mechanism of muon energy loss is by delta ray emission and elastic scattering processes.

FIGURE 11

6 Results

Following the procedures explained in the previous sections of this work, it has been possible to estimate the values of ADC channels, , , and . The values of , and are extracted from the Geant4 simulation as it has already been optimised for muons in this energy range. The physics in the FLUKA simulation is currently being tuned, but the first results presented in Section 5 are qualitatively promising. Hence, the value of the calibration factor in Eq. 1 is

Fitting the full-rate beam histogram of the charge deposited in the hodoscope in each muon spill with a Gaussian profile (see Figure 12) made it possible to extract the average value of for the analysed run. However, this flux estimation can also be carried out event-by-event simply taking the punctual value of and converting it into punctual muon flux. Taking the mean value and converting it into mean muon rate by applying Eq. 1 with the value of in Eq. 3, one gets muons/s. This value has been obtained with synchrotron current % the maximum value. It is consistent with the expected order of magnitude for the 55 MeV/c negative muon flux at full synchrotron current ( muons/s) (Matsuzaki et al. (2001); Hillier et al. (2019)).

FIGURE 12

By extracting weighing factors for at other momenta from the simulation, it has also been possible to estimate the muon flux during high-rate measurements at momentum different from 55 MeV/c. The result is presented in Figure 13, and the trend is increasing with momentum, as expected.

FIGURE 13

7 Conclusion

A full calibration protocol for a beam hodoscope that will be used as a flux monitor has been explained, applied, and tested.

When the beam repetition rate is known, the values of , and are extracted from simulation, and the value of is extracted from low-rate measurements allowing single-particle events. In particular, by modelling the detector in a simulation toolkit (Geant4/FLUKA), and are obtained by counting the primary muons interacting with fibres in both planes of the detector or with one fibre only, respectively, whereas is the ratio of the mean deposited energies in the two cases. On the other hand, is obtained by tuning the beampipe (bending and quadrupole) magnets to deliver a small fraction of the beam. By tuning the magnets to allow single-particle spills to reach the detector, it was possible to determine the amount of charge deposited by muons hitting two fibres by imposing coincidence between the two planes during data analysis. This allowed to calculate the calibration constant to convert the high-rate deposited charge into the muon rate . The muon rate plays a crucial role in the data normalisation for the FAMU experiment.

This protocol can be applied to similar detectors to use them for the same particle counting task in cases where single-particle discrimination is not possible.

This procedure can be carried out with any scintillating fibre-based hodoscope that will be used as a muon beamline monitor with various applications. For example, muonic atom X-ray spectroscopy (-XES), used for elemental and isotopic analysis (Clemenza et al. (2019); Cataldo et al. (2022); Rossini et al. (2023c)), is a non-destructive technique for the depth-dependent characterisation of materials of interest such as Cultural Heritage samples. It consists of a spectroscopic analysis of the muonic atom X-rays emitted by a sample. Hence, knowing the injected beam rate would give important information about the number of atoms created, helping quantify the elements and isotopes present in the sample. The presence of a calibrated beam hodoscope in such applications would, therefore, help improve the technique.

References

• 1

  AgostinelliS.AllisonJ.AmakoK.ApostolakisJ.AraujoH.ArceP.et al (2003). Geant4—a simulation toolkit. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.506, 250–303. 10.1016/s0168-9002(03)01368-81016/S0168-9002(03)01368-8

  • CrossRef
  • Google Scholar

• 2

  AhdidaC.BozzatoD.CalzolariD.CeruttiF.CharitonidisN.CimminoA.et al (2022). New capabilities of the fluka multi-purpose code. Front. Phys.9, 788253. 10.3389/fphy.2021.788253

  • CrossRef
  • Google Scholar

• 3

  AntogniniA.HagelsteinF.PascalutsaV. (2022). The proton structure in and out of muonic hydrogen. Annu. Rev. Nucl. Part. Sci.72, 389–418. 10.1146/annurev-nucl-101920-024709annurev-nucl-101920-024709

  • CrossRef
  • Google Scholar

• 4

  BattistoniG.BoehlenT.CeruttiF.ChinP. W.EspositoL. S.FassòA.et al (2015). Overview of the fluka code. Ann. Nucl. Energy82, 10–18. 10.1016/j.anucene.2014.11.007

  • CrossRef
  • Google Scholar

• 5

  BonesiniM.BenocciR.BertoniR.MazzaR.deBariA.MenegolliA.et al (2019). The upgraded beam monitor system of the famu experiment at riken–ral. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.936, 592–594. 10.1016/j.nima.2018.08.092092

  • CrossRef
  • Google Scholar

• 6

  BonesiniM.BertoniR.ChignoliF.MazzaR.CerviT.de BariA.et al (2017). The construction of the fiber-sipm beam monitor system of the r484 and r582 experiments at the riken-ral muon facility. J. Instrum. (JINST)12, C03035. 10.1088/1748-0221/12/03/c030351748-0221/12/03/C03035

  • CrossRef
  • Google Scholar

• 7

  BrunR.RademakersF. (1997). Root — an object oriented data analysis framework. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.389, 81–86. 10.1016/s0168-9002(97)00048-x1016/S0168-9002(97)00048-X

  • CrossRef
  • Google Scholar

• 8

  CarboneR.BonesiniM.BertoniR.MazzaR.RossellaM.TortoraL.et al (2015). The fiber-sipmt beam monitor of the r484 experiment at the riken-ral muon facility. J. Instrum. (JINST)10, C03007. 10.1088/1748-0221/10/03/C03007

  • CrossRef
  • Google Scholar

• 9

  CarlsonC. (2015). The proton radius puzzle. Prog. Part. Nucl. Phys.82, 59–77. 10.1016/j.ppnp.2015.01.002

  • CrossRef
  • Google Scholar

• 10

  CataldoM.ClemenzaM.IshidaK.HillierA. (2022). A novel non-destructive technique for cultural heritage: depth profiling and elemental analysis underneath the surface with negative muons. Appl. Sci.12, 4237. 10.3390/app12094237app12094237

  • CrossRef
  • Google Scholar

• 11

  ClemenzaM.BaldazziG.BalleriniG.BonesiniM.CarpinelliM.CremonesiO.et al (2019). Chnet-tandem experiment: use of negative muons at riken-ral port4 for elemental characterization of “nuragic votive ship” samples. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.936, 27–28. 10.1016/j.nima.2018.11.076

  • CrossRef
  • Google Scholar

• 12

  Dal MasoG.BarchettiF.FrancesconiM.GalliL.GreuterU.HildebrandtM.et al (2023). Beam monitoring detectors for high intensity muon beams. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.1047, 167739. 10.1016/j.nima.2022.167739

  • CrossRef
  • Google Scholar

• 13

  HillierA. D.LordJ. S.IshidaK.RogersC. (2019). Muons at isis. Philosophical Trans. R. Soc. A Math. Phys. Eng. Sci.377, 20180064. 10.1098/rsta.2018.0064

  • CrossRef
  • Google Scholar

• 14

  LordJ.McKenzieI.BakerP.BlundellS. J.CottrellS. P.GiblinS. R.et al (2011). Design and commissioning of a high magnetic field muon spin relaxation spectrometer at the isis pulsed neutron and muon source. Rev. Sci. Instrum.82, 073904. 10.1063/1.3608114

  • CrossRef
  • Google Scholar

• 15

  MatsuzakiT.IshidaK.NagamineK.WatanabeI.EatonG.WilliamsW. (2001). The riken-ral pulsed muon facility. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.465, 365–383. 10.1016/S0168-9002(01)00694-5

  • CrossRef
  • Google Scholar

• 16

  PapaA.BarchettiF.GrayF.RipicciniE.RutarG. (2015). A multi-purposed detector with silicon photomultiplier readout of scintillating fibers. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.787, 130–133. 10.1016/j.nima.2014.11.074

  • CrossRef
  • Google Scholar

• 17

  PapaA.RutarG.BarchettiF.HildebrandtM.KettleP. (2019). A fast and quasi non-invasive muon beam monitor working at the intensity frontier. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.936, 634–635. 10.1016/j.nima.2018.10.145

  • CrossRef
  • Google Scholar

• 18

  PizzolottoC.AdamczakA.BakalovF. D.et al (2020). The famu experiment: muonic hydrogen high precision spectroscopy studies. Eur. Phys. J.A 56, 185. 10.1140/epja/s10050-020-00195-9

  • CrossRef
  • Google Scholar

• 19

  RossiniR.AdamczakA.BakalovD.BaldazziG.BanfiS.BaruzzoM.et al (2024a). Status of the detector setup for the famu experiment at riken-ral for a precision measurement of the zemach radius of the proton in muonic hydrogen. J. Instrum. (JINST)19, C02034. 10.1088/1748-0221/19/02/C02034

  • CrossRef
  • Google Scholar

• 20

  RossiniR.BenocciR.BertoniR.BonesiniM.CarsiS.ClemenzaM.et al (2024b). Characterisation of a low-momentum high-rate muon beam monitor for the famu experiment at the cnao-xpr beam facility. J. Instrum. (JINST)19, C01024. 10.1088/1748-0221/19/01/C01024

  • CrossRef
  • Google Scholar

• 21

  RossiniR.BenocciR.BertoniR.BonesiniM.ClemenzaM.De VecchiC.et al (2023a). Characterisation of a scintillating fibre-based hodoscope exposed to the cnao low-energy proton beam. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.1046, 167746. 10.1016/j.nima.2022.167746

  • CrossRef
  • Google Scholar

• 22

  RossiniR.BenocciR.BertoniR.BonesiniM.ClemenzaM.MenegolliA.et al (2023b). Characterisation of muon and proton beam monitors based on scintillating fibres with a sipm read-out. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip.1046, 167684. 10.1016/j.nima.2022.167684

  • CrossRef
  • Google Scholar

• 23

  RossiniR.Di MartinoD.AgoroT.CataldoM.GoriniG.HillierA. D.et al (2023c). A new multidisciplinary non-destructive protocol for the analysis of stony meteorites: gamma spectroscopy, neutron and muon techniques supported by Raman microscopy and sem-eds. J. Anal. At. Spectrom.38, 293–302. 10.1039/D2JA00263A

  • CrossRef
  • Google Scholar

• 24

  VacchiA.MocchiuttiE.AdamczakA.BakalovD.BaldazziG.BaruzzoM.et al (2023). Investigating the proton structure: the famu experiment. Nucl. Phys. News33, 9–16. 10.1080/10619127.2023.21989132198913

  • CrossRef
  • Google Scholar

• 25

  VlachoudisV. (2009). “Flair: a powerful but user friendly graphical interface for FLUKA,” in International Conference on Mathematics, Computational Methods & Reactor Physics 2009, Saratoga Springs, New York, 3-7 May 2009, 790–800.

  • Google Scholar