Cosmic ray velocity and electric charge measurements with the AMS/RICH detector: prototype results

Cosmic ray velocity and electric charge measurements with the AMS/RICH detector: prototype results


The Alpha Magnetic Spectrometer (AMS) to be installed on the International Space Station (ISS) will measure charged cosmic ray spectra of elements up to iron, in the rigidity range from 1 GV to 1 TV, for at least three years. AMS is a large angular spectrometer composed of different subdetectors, including a proximity focusing Ring Imaging CHerenkov (RICH) detector. This will be equipped with a mixed radiator made of aerogel and sodium fluoride (NaF), a lateral conical mirror and a detection plane made of 680 photomultipliers coupled to light guides. The RICH detector allows measurements of particle’s electric charge up to iron, and particle’s velocity. Two possible methods for reconstructing the Čerenkov angle and the electric charge with the RICH will be discussed.

A RICH prototype consisting of a detection matrix with 96 photomultipliers, a segment of a conical mirror and samples of the radiator materials was built and its performance was evaluated using ion beam data. Results from the last test beam performed with ion fragments resulting from the collision of a 158 GeV/c/nucleon primary beam of indium ions (CERN SPS) on a lead target are reported. The large amount of collected data allowed to test and characterize different aerogel samples and the NaF radiator. In addition, the reflectivity of the mirror was evaluated. The data analysis confirms the design goals.

I The AMS-02 experiment

The Alpha Magnetic Spectrometer [1] (AMS) is a particle detector to be installed in the International Space Station (ISS) for at least three years. The spectrometer will be able to measure the rigidity (), the charge (), the velocity () and the energy () of cosmic rays from some MeV up to 1 TeV within a geometrical acceptance of 0.5 Figure 1 shows a schematic view of the AMS spectrometer. At both ends of the AMS spectrometer exist the Transition Radiation Detector (TRD) (top) and the Electromagnetic Calorimeter (ECAL) (bottom). Both will provide AMS with capability to discriminate between leptons and hadrons. Additionally the calorimeter will trigger and detect photons. The TRD will be followed by the first of the four Time-of-Flight (TOF) system scintillator planes. The TOF system [4] is composed of four roughly circular planes of 12 cm wide scintillator paddles, one pair of planes above the magnet, the upper TOF, and one pair bellow, the lower TOF. There will be a total of 34 paddles. The TOF will provide a fast trigger within 200 ns, charge and velocity measurements for charged particles, as well as information on their direction of incidence. The TOF operation at regions with very intense magnetic fields forces the use of shielded fine mesh phototubes and the optimization of the light guides geometry, with some of them twisted and bent. Moreover the system guarantees redundancy, with two photomultipliers on each end of the paddles and double redundant electronics. A time resolution of 140 ps for protons is expected [4].

The tracking system will be surrounded by veto counters and embedded in a magnetic field of about 0.9 Tesla produced by a superconducting magnet. It will consist on a Silicon Tracker [3], made of 8 layers of double sided silicon sensors with a total area of 6.7 m. There will be a total of 2500 silicon sensors arranged on 192 ladders. The position of the charged particles crossing the tracker layers is measured with a precision of 10 m along the bending plane and 30 m on the transverse direction. With a bending power (BL) of around 0.9 T.m, particles rigidity is measured with an accuracy better than 2% up to 20 GV and the maximal detectable rigidity is around 1-2 TV. Electric charge is also measured from energy deposition up to Z26.

The Ring Imaging Čerenkov Detector (RICH) [2] will be located right after the last TOF plane and before the electromagnetic calorimeter. It will be described in detail in next section.

Fig. 1: A whole view of the AMS Spectrometer.

The long stay of AMS in space and its large acceptance will allow the accumulation of a large statistic of events increasing in several orders of magnitude the sensitivity of the proposed physical measurements. With an average collection rate of 1000 events per second, a total of protons per year and around antiprotons will be accumulated.

The main goals of the AMS-02 experiment are:

  • A precise measurement of the charged cosmic ray spectrum between 100 MeV and 1 TeV, and the detection of photons up to a few hundred GeV;

  • A search for heavy antinuclei (), which if discovered would signal the existence of primordial antimatter;

  • Search for non-baryonic dark matter through the detection of annihilation products appearing as anomalies of the cosmic-ray spectra (e, , and );

Ii The AMS RICH detector

The RICH is a proximity focusing device with a dual radiator configuration on the top made of 92 aerogel 25 mm thick tiles with a refractive index 1.050 and sodium fluoride (NaF) tiles with a thickness of 5 mm in the center covering an area of 3434 cm. The NaF placement prevents the loss of photons in the hole existing in the center of the readout plane ( cm), in front of the ECAL calorimeter located below. The radiator tiles are supported by a 1 mm thick layer of methacrylate (n=1.5) free of UV absorbing additives.

The detection matrix is composed of 680 multiplixelized photon readout cells each consisting of a photomultiplier coupled to a light guide, HV divider plus front-end (FE) electronics, all housed and potted in a plastic shell and then enclosed in a magnetic shielding with a thickness varying from 0.8 to 1.2 mm. The photon detection is made with an array of multianode Hamamatsu tubes (R7600-00-M16) with a spectral response ranging from 300 to 650 nm and a maximum quantum efficiency at 420 nm. To increase the photon collection efficiency, a light guide consisting of 16 solid acrylic pipes glued to a thin top layer (1 mm) was produced. It is optically coupled to the active area of phototube cathode through a mm flexible optical pad. With a total height of mm and a collecting surface of mm, it presents a readout pixel size of mm and a pitch of 37 mm. The light guide is mechanically attached through nylon wires to the photomultiplier polycarbonate housing.

A high reflectivity conical mirror surrounds the whole set. The mirror was included to increase the device acceptance since around 33% of the aerogel generated photons impact on the mirror. It consists of a carbon fiber reinforced composite substrate with a multilayer coating made of aluminium and SiO deposited on the inner surface. This ensures a reflectivity higher than 85% for 420 nm wavelength photons.

The RICH has a truncated conical shape with an expansion height of 46.3 cm, a top radius of 60 cm and a bottom radius of 67 cm. The total height of the detector is 60.5 cm and it covers 80% of the AMS magnet acceptance. Figure 2 shows a schematic view of the RICH detector.

RICH was designed to measure the velocity () of singly charged particles with a resolution of 0.1, to extend the charge separation (Z) up to iron (Z=26), to contribute to separation and to albedo rejection.

In order to validate the RICH design, a prototype with an array of 911 cells filled with 96 photomultiplier readout units similar to part of the matrix of the final model was constructed. The performance of this prototype has been tested with cosmic muons and with a beam of secondary ions at the CERN SPS produced by fragmentation of a primary beam in 2002 and 2003. The light guides used were prototypes with a slightly smaller collecting area (3131 mm). Different samples of the radiator materials were tested and placed at an adjustable supporting structure. Different expansion heights were set in order to have fully contained photon rings on the detection matrix like in the flight design. A segment of a conical mirror with 1/12 of the final azimuthal coverage, which is shown in left picture of Figure 3, was also tested.

The RICH assembly has already started at CIEMAT in Spain and is foreseen to be finished in July 2007. A rectangular grid has already been assembled and has been subject to a mechanical fit test, functional tests, vibration tests and vacuum tests. The other grids will follow. The refractive index of the aerogel tiles is being measured and the radiator container was subjected to a mechanical test. The final integration of RICH in AMS will take place at CERN in 2008.

Fig. 2: Schematic view of the RICH detector.
Fig. 3: Protype with reflector (left). Top view of the test beam 2003 experimental setup using CERN SPS facility (right).

Iii Beta () and charge () reconstruction

A charged particle crossing a dielectric material of refractive index with a velocity , greater than the speed of light in that medium, emits photons. The aperture angle of the emitted photons with respect to the radiating particle is known as the Čerenkov angle, , and it is given by (see Ref. [5])


It follows that the velocity of the particle, , is straightforward derived from the Čerenkov angle reconstruction, which is based on a fit to the pattern of the detected photons. Complex photon patterns can occur at the detector plane due to mirror reflected photons, as can be seen on the left display of Figure 4. The event shown is generated by a simulated beryllium nucleus crossing the NaF radiator.

Fig. 4: Beryllium event display generated in a NaF radiator. The reconstructed photon pattern (full line) includes both reflected and non-reflected branches. The outer circular line corresponds to the lower boundary of the conical mirror. The square is the limit of the non-active region (left). Hit’s residuals distribution belonging to photon rings generated in an aerogel radiator, n=1.05, 2.5 cm thick.

The Čerenkov angle reconstruction procedure relies on the higly accurate information (see Section I) of the particle direction (,) provided by the tracker. The tagging of the hits signaling the passage of the particle through the solid light guides in the detection plane provides an additional track element, however, those hits are excluded from the reconstruction. The best value of will result from the maximization of a Likelihood function, built as the product of the probabilities, , that the detected hits belong to a given (hypothetical) Čerenkov photon pattern ring,


This probability takes into account , the closest distance of the hit to the Čerenkov pattern and the signal strength. The probability of a hit belonging to the pattern is obtained by taking into account that it can either be part of the noise (essentially flat) or of the Čerenkov pattern (Gaussian distributed). Expressing as the photon background fraction, as the detector’s dimensions and as the width of the residuals distribution (as can be seen in right plot of Figure 4), we can write:


where and respectively for aerogel and NaF in the prototype setup. The detector’s dimension is defined as D=100 cm for flight setup and 20 cm for prototype setup. For a more complete description of the method see Ref. [5].

The Čerenkov photons produced in the radiator are uniformly emitted along the particle path inside the dielectric medium, , and their number per unit of energy depends on the particle’s charge, , and velocity, , and on the refractive index, , according to the expression:


Therefore to reconstruct the charge the following procedure is required:

  • Čerenkov angle reconstruction.

  • Estimation of the particle path, , which relies on the information of the particle direction provided by the tracker.

  • Counting the number of photoelectrons. The number of photoelectrons related to the Čerenkov ring has to be counted within a fiducial area, in order to exclude the uncorrelated background. In particular it ensures the exclusion of photons which are scattered in the radiator. A distance of 13 mm to the ring was defined as the limit for photoelectron counting, corresponding to a ring width of 5 pixels.

  • Evaluation of the photon detection efficiency. The number of radiated photons () which will be detected () is reduced due to the interactions with the radiator (), the photon ring acceptance (), light guide efficiency () and photomultiplier efficiency ().


The charge is then calculated according to expression 4, where the normalization constant can be evaluated from a calibrated beam of charged particles. For a more complete description of the charge reconstruction method see Ref. [5].

Iv Results with the RICH prototype

The RICH prototype was subject to cosmic muons and to in-beam tests using secondary nuclei from fragmentation of 20 GeV/c/nucleon lead (Pb) ions in a beryllium target and 158 GeV/c/nucleon indium nuclei in a Pb target from the CERN SPS in 2002 and 2003, respectively [6].

In 2003 a monocromatic particle beam with a momentum resolution 0.15%1.5% was obtained. The optics of the line was tuned to provide a beam as parallel as possible, with a divergence less than 1 mrad. The beam section was 1 mm for the narrow beam runs and 1 cm for the spread beam runs.

The beam nuclear composition could be selected according to the desired A/Z value of the fragmentation products by setting the beam line rigidity at the appropriate value. Three main selection values were established: A/Z=2 (He, Li, B, C,…); A/Z=2.25 to enhance the Be peak and A/Z=2.35 to enhance the indium peak.

Figure 3 (right) shows a general view of the 2003 test beam setup in the experimental area H8-SPS at CERN. The prototype was placed inside a light-tight container. The setup was completed with AMS silicon tracker layers placed upstream in the beam, a TOF prototype placed downstream, two multi-wire proportional chambers, two organic scintillator counters and during a certain period a plastic Čerenkov counter. The two scintillators placed 1 m apart in front of the prototype cointainer, provided the DAQ trigger as well as an independent charge measurement. The silicon tracker prototype provided a very precise measurement of the particle track parameters for the event reconstruction as well as an external selection of charge.

The purposes of the tests were testing flight front-end electronics, characterize the performance of the aerogel and the NaF radiators, estimating the mirror reflectivity and evaluating the global functionality of the prototype. A total number of 11 million events were recorded during eleven days. Different particle incidences were obtained by rotating the prototype setup with respect to the beam line (0, 5, 10,15,20).

The event selection was mainly intended to remove wrongly reconstructed tracks, events with clusterized hits and events arising from later fragmentation. First, consistency between the external determination of the track transverse coordinates and the estimation from the reconstructed ring is required. Then events with more than one particle cluster in the detection matrix are rejected. Futhermore, the Kolmogorov probability of the event is calculated requiring an uniform azimuthal distribution of the ring hits.

The evaluation of the aerogel samples in order to make a final radiator choice was one of the key issues of these tests. Different production batches from two manufacturers, Matsushita Electric Co. (MEC) and Catalysis Institute of Novosibirsk (CIN) with different refractive indexes, 1.03 and 1.05, were analyzed. The required criteria for a good candidate were a high photon yield, in order to ensure a good ring reconstruction efficiency and accurate and charge measurements.

The aerogel light yield depends on the tile thickness and on its optical properties (refractive index and clarity). The light yield has been evaluated from the analysis of helium samples collected in 2003 and from the analysis of proton data samples gathered in 2002 with different beam momentum between 5 and 13 GeV/c [7].

Fig. 5: Light yield as function of proton beam momentum for the different aerogel samples tested in 2002 (left). Light yield comparison based on test beam data 2002 and 2003. All values were extrapolated for fully cointained rings generated by a particle with 1 and in an aerogel radiator with a common thickness of 3 cm (right).

Figure 5 (left) shows the evolution of the light yield of the different aerogel samples tested in 2002 with the proton beam momentum. A fit to each set of data was applied and the light yield for a proton with 1, generating fully contained rings in a radiator with a common thickness of 3 cm was extrapolated. Right plot of the same figure shows the normalized to 3 cm thickness light yield for the different aerogel samples tested in 2002 and 2003. Two interesting features are enhanced. In one hand, the same sample of CINy02.103 was used in both years and its light yield analysis shows the same value which proves the setup stability and the aerogel good performance after one year period; on the other hand it is notorious the highest signal comes from a CIN sample produced in 2003 with refractive index reflecting the very good clarity (m/cm) of the aerogel batch.

The resolution of the measurement, obtained as explained in Section III, was estimated using a Gaussian fit to the reconstructed spectrum, shown in left plot of Figure 6 for helium nuclei. Data were collected with the aerogel radiator CINy03.105, 2.5 cm thick together with an expansion height of 35.31 cm. The events shown correspond to particles inciding vertically and generating fully contained rings. The beta reconstructed from a simulated helium data is also shown (shaded histogram) superimposed with good agreement between data and Monte Carlo.

{overpic} DATAMC

Fig. 6: Comparison of the distribution for helium data (black dots) and simulation (shaded) (left). Evolution of the resolution with the charge obtained for the same aerogel radiator. Simulated points for Z=2, 6, 16 are marked with full squares (right).

The charge dependence of the velocity relative resolution for the same radiator is shown in the right plot of Figure 6. The different charges were selected using external and independent measurements performed by the silicon tracker prototype and by the two scintillators. The observed resolution varies according to a law 1/Z, as it is expected from the charge dependence of the photon yield in the Čerenkov emission, up to a saturation limit set by the pixel size of the detection unit cell. The function used to perform the fit is the following:


where, means the resolution for a singly charged particle while means the resolution for a very high charge generating a large number of hits, here the resolution is dominated by the pixel size as mentioned before. The fitted values are and for the run conditions stated above. Simulated data points for Z=2, 6, 16 are marked upon the same plot with full squares. Once more the agreement between data and Monte Carlo measurements for charges different of Z=2 is good.

The beta resolution for a 1, helium nuclei impacting in each of the aerogel samples tested in 2003 are summarized in Table I. The results are for a common expansion height extrapolated from the values measured at the adjusted heights.

All the tested radiators fulfill the RICH requirement for measurement.

resolution for Z=2, H=33.5 cm
radiator CIN103 MEC103 CIN105
0.4210.003 0.4350.002 0.4590.004
TABLE I: Beta resolution for a helium particle with 1 obtained for all the aerogel samples tested in 2003 and extrapolated for a common expansion height of 33.5 cm.

The distribution of the reconstructed charges in an aerogel radiator of n=1.05, 2.5 cm thick is shown in left plot of Figure 7. The reconstruction method used was the one described in Section III. The spectrum enhances a structure of well separated individual charge peaks over the whole range up to iron (Z=26). This spectrum was measured for a beam selection of A/Z=9/4.

The charge resolution for each nuclei, shown in right panel of Figure 7, was evaluated through individual Gaussian fits to the reconstructed charge peaks selected by the independent measurements performed by the scintillators and silicon tracker detectors. A charge resolution for proton events slightly better than charge units is achieved and as expected the best charge resolution is provided by this radiator due to its higher photon yield.

The charge resolution as function of the charge Z of the particle follows a curve that corresponds to the error propagation on Z which can be expressed as:


This expression describes the two distinct types of uncertainties that affect Z measurement: the statistical and the systematic. The statistical term is independent of the nuclei charge and depends essentially on the amount of Čerenkov signal detected for singly charged particles () and on the resolution of the single photoelectron peak (). The systematic uncertainty scales with Z, dominates for higher charges and is around 1%. It appears due to non-uniformities at the radiator level coming from variations in the refractive index, tile thickness or clarity or due to non-uniformities at the photon detection efficiency like PMT temperature effects or light guide non-uniformities.

The RICH goal of a good charge separation in a wide range of nuclei charges implies a good mapping and monitoring of the potential non-uniformities present on the detector. In order to keep the systematic uncertainties below 1%, the aerogel tile thickness, the refractive index and the clarity should not have a spread greater than 0.25 mm, and , respectively; at the detection level a precise knowledge ( level) of the single unit cell photo-detection efficiency and gains is required.

Fig. 7: Charge peaks distribution measured with the RICH prototype using a n=1.05 aerogel radiator, 2.5 cm thick. Individual peaks are identified up to Z26 (left). Charge resolution versus particle Z for the same aerogel radiator. The curve gives the expected value estimated as explained in the text (right).

Runs with a mirror prototype were also performed and its reflectivity was derived from data analysis. The obtained value is in good agreement with the design value.

V Conclusions

AMS-02 will be equipped with a proximity focusing RICH detector based on a mixed radiator of aerogel and sodium fluoride, enabling velocity measurements with a resolution of about 0.1% and extending the charge measurements up to the iron element. Velocity reconstruction is made with a likelihood method. Charge reconstruction is made in an event-by-event basis. Evaluation of both algorithms on real data taken with in-beam tests at CERN, in October 2003 was done. The detector design was validated and a refractive index 1.05 aerogel was chosen for the radiator, fulfilling both the demand for a large light yield and a good velocity resolution. The RICH detetector is being constructed and its assembling to the AMS complete setup is foreseen for 2008.


  1. C. Lechanoine-Leluc, Proc. 29th ICRC (Pune) 3, 381-384 (2005).
  2. F. Barão, L. Arruda, et al., Proc. 29th ICRC (Pune) 9, 299-302 (2005).
  3. C. Lechanoine-Leluc, Proc. 29th ICRC (Pune) 9, 299-302 (2005).
  4. D. Casadei et al., Nuclear Physics B (Proc. Suppl.) 113,133 (2002).
  5. F. Barao, L. Arruda et al., NIM A502,310 (2003).
  6. P. Aguayo, L.Arruda et al., NIM A560, 291-302 (2006).
  7. R. Pereira, L.Arruda et al., In-beam aerogel light yield characterization. for the AMS/RICH detector , AMS internal note in preparation
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description