First observation of the cosmic ray shadow of the Moon and the Sun with KM3NeT/ORCA

This article reports the first observation of the Moon and the Sun shadows in the sky distribution of cosmic-ray induced muons measured by the KM3NeT/ORCA detector. The analysed data-taking period spans from February 2020 to November 2021, when the detector had 6 Detection Units deployed at the bottom of the Mediterranean Sea, each composed of 18 Digital Optical Modules. The shadows induced by the Moon and the Sun were detected at their nominal position with a statistical significance of 4.2σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} and 6.2σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}, and an angular resolution of σres=0.49∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{res}=0.49^\circ $$\end{document} and σres=0.66∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{res}=0.66^\circ $$\end{document}, respectively, consistent with the prediction of 0.53∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.53^\circ $$\end{document} from simulations. This early result confirms the effectiveness of the detector calibration, in time, position and orientation and the accuracy of the event direction reconstruction. This also demonstrates the performance and the competitiveness of the detector in terms of pointing accuracy and angular resolution.


Abstract
This article reports the first observation of the Moon and the Sun shadows in the sky distribution of cosmic-ray induced muons measured by the KM3NeT/ORCA detector.The analysed data-taking period spans from February 2020 to November 2021, when the detector had 6 Detection Units deployed at the bottom of the Mediterranean Sea, each composed of 18 Digital Optical Modules.The shadows induced by the Moon and the Sun were detected at their nominal position with a statistical significance of 4.2σ and 6.2σ, and an angular resolution of σ res = 0.49 • and σ res = 0.66 • , respectively, consistent with the prediction of 0.53 • from simulations.This early result confirms the effectiveness of the detector calibration, in time, position and orientation and the accuracy of the event direction reconstruction.This also demonstrates the performance and the competitiveness of the detector in terms of pointing accuracy and angular resolution.Cosmic rays (CR) are charged particles mainly composed of protons and light nuclei.Since they can be deflected by irregular Galactic magnetic fields, their arrival directions at the Earth are almost isotropic.Once a primary CR particle reaches the Earth's upper atmosphere and interacts with an air nucleus, it produces secondary particles.The most penetrating component of these are muons, which can be detected at the surface of the Earth but also at underground or underwater detectors.In order to reach the KM3NeT/ORCA detector [1] at a depth of 2500 meters below sea level [2], vertically down-going muons need a minimal energy of around 500 GeV at sea level and thus they must have originated from primary CRs of energies exceeding several TeV/nucleon as described in [3].90% of the CR particles yielding muons used in this analysis have energies between 3 and 330 TeV as illustrated in Figure 1.The angle between the primary CR and the secondary muon is on average within 0.1°at these energies [4] which is around 5 times lower than the expected angular resolution of the detector for muons.At the KM3NeT/ORCA detector the muons are dominantly minimal-ionizing with energies ranging from few tens of GeV to few hundred GeV.CRs are blocked by nearby celestial bodies such as the Moon and the Sun.This induces a deficit in the atmospheric muon flux and in other secondary CR particles coming from the direction of these objects.This effect had been predicted by Clark in 1957 [5].Its observation can be used to verify the pointing accuracy and angular resolution of detectors which are able to measure secondary particles from CR interactions.The Moon and the Sun shadows in CRs have been observed by multiple experiments (IceCube [6], ANTARES [7,8], MACRO [9], L3 [10], MINOS [11], ARGO-YBJ [12], HAWC [13] and others).As a permanent bright high-energy neutrino source has yet to be found, the observation of the Moon and Sun shadows in CRs is an important calibration tool for neutrino telescopes such as KM3NeT/ORCA and helps to demonstrate their pointing accuracy and to measure their angular resolution.

KM3NeT/ORCA detector
KM3NeT is a research infrastructure consisting of undersea Cherenkov neutrino telescopes currently under construction at the bottom of the Mediterranean Sea off-shore the Italian Sicily coast (KM3NeT/ARCA) and 40 km offshore Toulon, France (KM3NeT/ORCA) [1].The two detectors are optimised for different neutrino energy ranges.They are composed of vertical Detection Units (DUs).Six of them had been operational in KM3NeT/ORCA when the data used in this analysis were acquired.Each DU consists of 18 spherical Digital Optical Modules (DOMs), with 31 photomultiplier tubes (PMTs) distributed almost isotropically within each DOM [14].These PMTs detect the Cherenkov light emitted along the path of relativistic charged particles propagating through water.In the data acquisition a hit is produced when a photon impinging on a PMT induces an electrical signal above a defined threshold.A hit consists of a time stamp and a time over threshold.An event is created when the trigger algorithm identifies a series of causally-connected hits.For the analysis presented here, these events are processed by a track reconstruction algorithm.The hits are fitted with a model of a Cherenkov light emitting muon.This particle is assumed to follow a long, straight trajectory and to propagate practically at the speed of light in vacuum through water.The position, time and direction of such a track is determined by using a maximum-likelihood method based on a set of causally-connected hit times and positions [15].A time-dependent calibration of the detector, with the monitored positions and orientations of every DOM computed and interpolated every 10 minutes, is used in this work to account for movements of the strings with the sea current.

Data and Monte Carlo samples
The data used in this analysis were collected between February 11, 2020, and November 18, 2021 for a total of 499.3 days.Quality cuts on the number of used hits and the likelihood of the track reconstruction were applied to remove poorly reconstructed events keeping 83% of the initial event sample.An average event density of 3000 events per square degree is measured in the vicinity of Moon and Sun.It is expected that about 640 CR events are blocked by each of the two sky objects.The position of the Moon/Sun in the sky is obtained using the astropy package [16] that relies on the International Celestial Reference System (ICRS) coordinates described in [17].The latitude of 43 • North of the detector and the data taking period of more than one year lead to broad zenith angle distributions of the selected events in the vicinity of Sun and Moon as shown in Figure 2. Monte Carlo (MC) simulations are used to optimise the track selection, and predict the angular resolution and amplitude of the CR deficit induced by the Moon and the Sun.Secondary muons at the KM3NeT/ORCA detector are simulated with the MUPAGE package [18] [19].These muons are then propagated through sea water and Cherenkov photons are created within a cylindrical volume surrounding the simulated detector.Finally the detector response is simulated by producing digitized hits from photons detected by PMTs and by adding noise hits from environmental background mainly due to 40 K decay and bioluminescence.These are derived from real data runs in a time-dependent way following the run-by-run approach previously introduced by the ANTARES Collaboration [3].The time-dependent PMT efficiencies are monitored and used in the simulation.The last two steps are done with KM3NeT custom software [15].The resulting hit patterns are passed through the trigger and reconstruction steps using the same software that is used to process real data.The simulated track sample is about 4 times larger compared to the real data sample.
The optimisation of the track selection is performed by varying the cut on the angular error estimate from the track reconstruction to maximize the significance of the shadow signal.By requiring a more and more stringent cut a smaller but higher quality event sample is selected.The chosen value keeps 45% of the original MC sample.

Analysis method
The search for the shadow of the Moon/Sun is done in the phase space of the angular differences between the reconstructed track coordinates and the celestial object.Both a 1-dimensional and a 2-dimensional analysis have been performed.The 1-dimensional analysis uses the space angle between the direction of the Moon/Sun and the track.The 2-dimensional analysis uses two Cartesian coordinates (x, y), starting from the zenith and azimuth angles of the sky object θ sk , α sk and the tracks θ tr , α tr respectively, 2D maps that extend to ±6 • in (x, y) have been constructed.Data are only used for times when the full ±6 • field around the Moon/Sun is above the horizon.An angular range of ±6 • contains the detectable shadow signal entirely and it allows for a simultaneous fit of the shadow and background parameters while analysing just a moderately-sized atmospheric muon sample.It has been verified that the results do not depend on the precise choice of the chosen angular range.A constant 0.
is used, where n i stands for the event count in the i-th space angle bin to be compared with the expectations N i,H under the H 0 and H 1 hypotheses.The difference in ∆χ 2 H1/H0 = χ 2 (H 1 ) − χ 2 (H 0 ) values is used to determine the probability to reject the null hypothesis and to extract the significance of the observation from it.The background event distribution in azimuth is found to be uniform, while the zenith angle dependency can be conveniently parametrized with a 2 nd order polynomial function.This yields for the 2D maps with ρ a constant track density per space angle and a 1 , a 2 parameters which are determined during minimisation.For the 1D maps, each bin contains events from an almost symmetric zenith angle range above and below the Moon/Sun position resulting in a uniform exposure, i.e.N i,H0 = ρ.
The event expectation in bin i for the signal plus background hypothesis H 1 is defined as where G i describes the deficit of secondary CR events due to the shadowing effects of Moon/Sun in bin i with coordinates x i , y i as a bi-dimensional Gaussian where A is the relative shadow amplitude.For A = 1 the number of blocked CR events correspond to ρπR 2 s .The angular resolution of the detector for the selected sample of CR events is measured by σ res , the angular width of the Gaussian shadow.The R s , x s , y s parameters are the apparent angular radius of the celestial object and the relative angular position of the Moon/Sun shadow with respect to their nominal positions.For the 1D maps, the term 5is replaced by δ 2 , the square of the angular distance between the track and the sky object from Equation 6with V tr and V sk the direction vector of the track and the sky object respectively.
The significance of the shadow is found by fitting ρ, a 1 , a 2 , A and σ res at (x s , y s ) = (0, 0) on the 2D or 1D event distribution respectively.Results from these fits are summarized in section 5.2.The position of the shadow is obtained by simultaneously fitting ρ, a 1 , a 2 , (x s , y s ) and A, with σ res fixed to its expectation value from MC. Results from these fits are found in Section 5.3.
The assumption of Gaussianity of the shadowing effect is the result of a few approximations.Firstly, the influence of the size of the Moon/Sun is neglected.This is acceptable as long as the angular resolution of the detector is larger than the angular radius of the Moon/Sun, a condition which is amply satisfied in the present case.In addition, the real point spread function (PSF) of the detector has a slightly different radial shape compared to a Gaussian function (see Figure 3, right) which will be accounted for by fitting the shadow amplitude A when comparing data to the H 1 hypothesis.It is observed, that the PSF is identical for the data samples selected for Moon and Sun.
Further, the PSF is perfectly symmetric in x and y as shown in Figure 3 (left), allowing for a reliable fit of x s and y s .The Moon radius varies between 0.245 • and 0.279 • resulting in an amplitude variation of ±14% around the mean value.However our data sample covers several Moon cycles, so the average value is used.A similar statement can be made for the Sun whose apparent radius varies between 0.262 • and 0.271 • during the year.

Background
The background distributions in x and y are shown in Figure 4 for the Sun as an example.Fits of the H 0 hypothesis using a constant, and 2 nd order polynomial function describe the data well.The MC predictions are also compatible with these functions.Figure 5 shows the 1D distributions of the event density as a function of the angular distance from the Moon/Sun.The fit results from the 1D and 2D fits at (x s , y s ) = (0, 0) are summarized in Table 1.The significances are derived from the ∆χ 2 H1/H0 with two degrees of freedom (A, σ res ).The values obtained for A and σ res can be compared to MC expectations of A = 0.90 ± 0.09 and σ res = (0.53 ± 0.04) • for the combined Moon/Sun sample.The fitted values of σ res are found compatible with the prediction from simulations.The deeper amplitude and higher significance of the Sun shadow is consistent with the effects of the particular structure of the Sun's magnetic field during the periods of low solar activity, whose dipole shape is expected to enhance the Sun shadowing effect [21].

Positional fits
The results from the fits of the 2D map in (x s , y s ) with σ res = 0.53 • are shown in Figure 6.The plot illustrates χ 2 (H 1 ) − χ 2 (H 0 ) in colour coding.The image of the shadow is clearly visible.Figure 7 shows the 1σ (68.3%), 2σ (95.4%), 3σ (99.7%) confidence contours for the two parameters (x s , y s ) around the best fit point x s = (0.11 ± 0.21) • , y s = (0.04 ± 0.13) • for the Moon and x s = (−0.01± 0.11) • , y s = (0.10 ± 0.12) • for the Sun.The true position of the Moon and the Sun in Figure 7 are contained within the 68% contours, yielding a 84% and 67% compatibility between the nominal and the best fit positions, calculated from the corresponding ∆χ 2 with two degrees of freedom (x s , y s ).The slightly different shapes of the contours for Moon and Sun can be entirely attributed to statistical fluctuations.

Figure 1 :
Figure 1: Distribution of CR events which yield a reconstructed event in the KM3NeT/ORCA detector.The white region indicates the energy range which provides 90% of the event sample.

Figure 2 :
Figure 2: Zenith angle distribution of selected events for the Moon (left) and the Sun (right).Data (black crosses) and simulations rescaled to data (blue histograms) are shown.

Figure 3 :
Figure 3: PSF from MC, i.e. the difference between the true muon direction and the reconstructed one in 2-dimensional angular space (left) and as function of the space angle difference (right, blue and magenta for Sun and Moon, respectively).A Gaussian PSF with σ = 0.53 • is shown in green for comparison.

Figure 4 :
Figure 4: Event distribution in x and y for the Sun data sample (black crosses) compared to the polynomial fits from the H 0 hypothesis (red lines) and MC predictions rescaled to data (blue histograms)

Figure 5 :
Figure 5: Event density as a function of the distance to the Moon on the left and the Sun on the right.Data (black crosses) are compared to the H 0 fit (red) and the H 1 fit (magenta).

Figure 6 :
Figure 6: colour-coded ∆χ 2 H1/H0 as a function of xs and ys for the Moon (left) and the Sun (right).

Figure 7 :
Figure 7: Confidence intervals derived from Figure 6 for the Moon (left) and the Sun (right).The black cross indicates the best fit point.
be University of Münster, Institut für Kernphysik, Wilhelm-Klemm-Str.9, Münster, 48149 Germany bf Utrecht University, Department of Physics and Astronomy, PO Box 80000, Utrecht, 3508 TA Netherlands bg UCLouvain, Centre for Cosmology, Particle Physics and Phenomenology, Chemin du Cyclotron, 2, Louvain-la-Neuve, 1349 Belgium bh INFN, CNAF, v.le C. Berti-Pichat, 6/2, Bologna, 40127 Italy [20]binning in x and y is chosen, however results are reproduced when choosing a smaller binning.The 1D histogram contains tracks that are closer than 4 • to the Moon/Sun position.The significance of the shadowing effect of Moon/Sun is determined with a likelihood ratio test, by comparing the likelihood of a background hypothesis model H 0 , with the likelihood of a signal plus background model H 1 which includes a shadowing effect.The Poisson likelihood with the definition in Ref[20]