Capability of the HAWC gammaray observatory for the indirect detection of ultrahigh energy neutrinos
Abstract
The detection of ultrahigh energy neutrinos, with energies in the range or above, is a topic of great interest in modern astroparticle physics. The importance comes from the fact that these neutrinos point back to the most energetic particle accelerators in the Universe, and provide information about their underlying acceleration mechanisms. Atmospheric neutrinos are a background for these challenging measurements, but their rate is expected to be negligible above 1 . In this work we describe the feasibility to study ultrahigh energy neutrinos based on the Earthskimming technique, by detecting the charged leptons produced in neutrinonucleon interactions in a high mass target. We propose to detect the charged leptons, or their decay products, with the High Altitude Water Cherenkov (HAWC) observatory, and use as a large mass target for the neutrino interactions the Pico de Orizaba volcano, the highest mountain in Mexico. In this work we develop an estimate of the detection rate using a geometrical model to calculate the effective area of the observatory. Our results show that it may be feasible to perform measurements of the ultrahigh energy neutrino flux from cosmic origin during the expected lifetime of the HAWC observatory.
1 Introduction
The first evidence of ultrahigh energy neutrinos (in the energy range) from extraterrestrial origin was recently reported [1]. This opened a new field in astroparticle physics that will allow the identification and characterization of the most powerful particle accelerators in the Universe. Neutrinos are not affected by the electromagnetic or strong interactions, and thus point back to the source where they were produced, unlike charged cosmic rays. Gammarays are another cosmic probe that provides information about the acceleration mechanisms that occur in astrophysical sources. Due to this, there are several dedicated instruments, both ground or space based, performing a continuous survey of the Universe characterizing gammaray sources, e.g. the Fermi Gammaray Space Telescope and the imaging atmospheric Cherenkov telescopes H.E.S.S., MAGIC, VERITAS and FACT.
The HAWC observatory was designed to detect and characterize the sources of highenergy gammarays in the energy range between 100 and 100 [2] and started full operations in April 2015. It is a ground based instrument that detects atmospheric showers by measuring with high precision the arrival time of the particles that compose the air showers. This is done via the Cherenkov light produced by the air shower particles as they enter the 300 water detector tanks that constitute the observatory, with a total water volume of 54 million . At energies above 100 , gammarays suffer strong absorption from pair production with photons from the cosmic microwave background radiation that strongly reduces their mean free path [3]. For this reason, ultrahigh energy neutrinos may be a better tool to study the most energetic extragalactic particle accelerators.
The Earthskimming technique to detect ultrahigh energy neutrinos has been proposed before, see for instance [4, 5, 6, 7, 8]. The method consist on using the interaction between a neutrino and a nucleon via the exchange of a W boson to produce a charged lepton of the same flavour as the incoming neutrino. Since the neutrinonucleon cross section is very small, a largemass target is needed. The natural candidate to produce such interactions is the Earth crust, either by focusing the searches on quasihorizontal neutrinos that travel along a chord inside the Earth or by using mountains as targets. The produced charged lepton travels essentially in the same direction as the neutrino. Thus, in the considered scheme, the charged lepton will travel either upward if it was moving through a chord inside the Earth crust or quasi horizontally if the neutrino passed through a mountain.
From the three families of leptons, the electron neutrinos are unfavourable for this type of study because the produced highenergy electrons initiate electromagnetic showers that are easily absorbed by the target mass shortly after production. For neutrinos, the relatively large mean life of the produced charged lepton combined with the ultrahigh energy will produce a detectable signal only as an ultraenergetic . The neutrinos are the ones that have attracted more interest from the experimental point of view. The reason is that, because of their very short mean life, even if they are very energetic the produced charged leptons will decay into secondary particles that would make the detection of the signal easier. However, the neutrinos have proven to be one of the most elusive particles of the Standard Model, with less than 15 detections up to now [9]. Ultrahigh energy neutrinos (with searches up to 72 ) have eluded direct detection so far, even after the analysis of three years of IceCube data [10]. Even though neutrinos are disfavoured in production mechanisms at the astrophysical acceleration sites, the neutrino flavour mixing that occurs in cosmological distances is expected to produce approximately equal proportions of all neutrino flavours at the Earth.
There have been already experimental attempts to detect ultrahigh energy neutrinos using the Earthskimming technique. For instance, the Pierre Auger observatory in Argentina used their surface detectors to look for the electromagnetic signature of extensive air showers initiated by the decay of charged leptons of energies that develop close to the detector [11, 12, 13], without finding candidate signals so far. There are also studies that propose to use the fluorescence detector of the Pierre Auger observatory to detect the decay in the atmosphere of charged leptons produced by ultrahigh energy neutrinos [14]; however this idea has not been implemented yet. The Ashra1 collaboration [15] searched for neutrino emission from a GRB in the  energy range using the Earthskimming technique. The Ashra1 experiment, located on the Mauna Loa volcano and facing the Mauna Kea volcano in the Hawaii island, aimed to detect neutrinos that converted into charged leptons inside Mauna Kea. Their method consisted on measuring the Cherenkov light emitted by the particles of the atmospheric shower initiated by the decay products of the charged lepton. Their analysis of the GRB081203A did not find signals associated to neutrinos in the  energy range. There have been also studies about the feasibility to use the MAGIC telescopes for the detection of neutrinos [16, 17], by pointing their telescopes below the horizon towards the sea, or by searching for reflections of Cherenkov light by the nearby ground, the sea or clouds [18, 19]. However, no experimental results have been published yet.
In this paper we propose to adapt the EarthSkimming technique to use it with the HAWC gammaray observatory, employing the Pico de Orizaba volcano as a target for the neutrinonucleon interactions. However, we propose not to follow the method explored so far of studying the decay products of a charged lepton on the atmosphere. Instead, we propose to reconstruct directly the trajectory of the charged lepton or their decay products as they travel through the HAWC detectors. In this way we do not only restrict our studies to neutrinos, but also include the possibility for the detection of ultrahigh energy neutrinos (a first suggestion of this method was mentioned in [20]). At ultrahigh energies, the produced charged leptons will have an energy approximately equal to that of the original neutrino [7, 20]. These charged leptons, or their decay products, could travel crossing several HAWC detectors depositing large amounts of Cherenkov light, well beyond the average left from both: atmospheric muons ( 30 photoelectrons () [21]) located far from the shower core used in the gamma/hadron discrimination algorithms of HAWC and the high noise, also associated to atmospheric muons, considered by the HAWC collaboration to be in the range of 10200 [22]. Here it is important to point out that the dynamical range of the HAWC electronics goes from a fraction of a up to thousands [23].
In order to demonstrate that our proposal to search for tracks produced by charged leptons or their boosted decay products is feasible, we performed some GEANT4 [24] simulations of the decay of charged leptons with an energy of 1 (that decay approximately 50 after their creation point), and studied the shower evolution in air. We choose an energy of 1 so the charged leptons decay quickly and are in the energy regime of our studies. After decay, the opening angle of the charged products is smaller than 0.7, thus, after the decay products have travelled 2050 (the approximate distance from the edge of the volcano to the HAWC array at 4090 a.s.l is of two kilometers), they would hit at most three columns of HAWC tanks. Moreover, over 95% of the secondary charged particles are contained within an opening angle smaller than 0.2, therefore producing large Cherenkov signals only within a single row of HAWC tanks, producing a clearly identifiable track. This simplified exercise was done for 1 charged ’s, making it easy to extrapolate the results to higher energy charged leptons, since the opening angle of the decay products is inversely proportional to the energy of the primary charged lepton [25].
Based on this information, we believe that the tracking method is possible in the search for charged leptons, produced by neutrinonucleon interactions, in the energy range. The amount of Cherenkov light that could be detected by the HAWC PMTs by such energetic particles is expected to be in the range of thousands to tens of thousands . However, given the dynamic range of the current HAWC electronics one could anticipate that the collected light could be used as a proxy to at least set a lower boundary on the energy of the incoming charged lepton, in a similar manner to what IceCube does for muons that pass through the detector.
The paper is organized as follows: in Section 2 we present the calculation of the flux of charged leptons that could be produced by neutrinonucleon interactions in the Pico de Orizaba volcano. In Section 3, we describe a method to calculate the HAWC effective area based on purely geometrical considerations, and evaluate it for different trigger conditions. Then, in Section 4, we present our results for the possible detection rate and address the issue of the expected background signals. In Section 5, we discuss how the current trigger of HAWC can be useful in selecting data for these studies and discuss further a possible background rejection strategy. Finally, the conclusions are presented in Section 6.
2 Calculation of the flux of charged leptons produced by Earthskimming neutrinos
In order to obtain an estimate of the number of ultrahigh energy charged leptons that could be produced via the Earthskimming technique, we follow the formalism developed in [7, 26]. However we use further simplifications due to the detection method that we propose in this paper, and also because of the energy of the neutrinos that we plan to study. The differences between the original formalism of [7, 26] and our implementation are pointed out in the text.
The number of charged leptons () that could be detected by the observatory is given by
(1) 
where is the charged lepton flux, is the effective area of the observatory, is the live time of the experiment and is the duty cycle for observations, which basically describes which fraction of the time the experiment is actually able to take data. In this section we describe the calculation of the flux of charged leptons produced by the neutrinonucleon interactions that occur while the neutrinos traverse the volcano. We start with the differential flux of ultrahigh energy neutrinos produced by astrophysical sources, which we consider to be isotropic. This differential flux is given by
(2) 
with () being respectively the energy, polar and azimuthal angle of the neutrinos. Since we are interested on integrating the flux on only a certain region (the one covered by the Pico de Orizaba volcano), and because we are considering an isotropic flux, we can simplify the differential neutrino flux to
(3) 
where is the solid angle covered by the volcano that is being used as the target for the neutrinonucleon interactions. After the proposed isotropic neutrino flux passes through the Pico de Orizaba volcano, the produced differential charged lepton flux is given by
(4) 
with () being respectively the energy, polar and azimuthal angle of the produced charged leptons. The relation between the differential fluxes of incoming neutrinos and the produced charged leptons is thus given by
(5) 
where is a function that physically represents the convolution of the probabilities of the different processes that need to take place in order that an ultrahigh energy neutrino converts into a charged lepton inside the target material and is able to escape the mountain. Given that we are interested in studying ultrahigh energy neutrinos (), the first simplification comes from the fact that the produced charged lepton will approximately follow the same direction as the original neutrino. The angle between the original neutrino and the produced charged lepton () has been estimated to be smaller than 1 for energies above 1 for ’s [15]; so, we expect this angle to be negligible. This makes that the angular dependence of can be approximated by delta functions.
(6) 
The energy dependent part of the function can be written as the integral along the path of the neutrino, and the corresponding charged lepton, inside the volcano
(7) 
where is the survival probability for a neutrino travelling a certain distance inside the volcano. This probability can be written as
(8) 
where is the charged current interaction length. Since the average width of the volcano () is much smaller than the interaction length, almost all of the neutrinos will traverse the whole mountain. That defines the integration limit of Equation 8. The interaction length is a function of the density of the medium that the neutrino travels through. In the case of the formalism developed in [7], where the trajectories of the neutrinos were across chords inside Earth, the authors had to consider variations of the Earth density. In this particular case, since we are interested on neutrinos that travel through volcanoes, it is reasonable to consider a constant density along the path of the neutrinos. Thus, the expression for the charged current interaction length can be written as follows
(9) 
where is Avogadro’s constant and is the charged current cross section. The second probability () that enters the calculation of is that of a neutrino to convert into a charged lepton through a charged current in an infinitesimal distance . This probability is given by
(10) 
The third process () that is taken into account in the function is the probability that the produced charged lepton is able to escape the volcano, taking into account the charged lepton energy losses in the medium. This process is described by a system of two coupled differential equations:
(11) 
(12) 
Where Equation 11 describes the energy loss processes, with describing the ionization energy loss and the radiative energy loss. According to the literature, e.g. [7], the effects of are negligible at the energy regime of interest of this work; so, we consider in Equation 11 that . In Equation 12, is the mass of the charged lepton, is the speed of light and is the charged lepton lifetime. A Monte Carlo study from [20] shows that for ’s with energies of 1 , after 6 water equivalent () (6 2.3 in “standard rock”, the average path length inside the volcano), i.e. in the extreme case where the charged lepton has to travel the average width of the volcano, and this value approaches unity as the energy of the charged lepton increases. Thus, we consider that for the energy regime studied in this work, it is appropriate that for ’s we can take a value of (see Table 1).
Energy bin []  

[]  
, 
The case of the survival probability for charged leptons is more complicated to evaluate, since most of the research has been done for energies above or at 100 [27, 28, 29]. We take as a base for our calculations the results presented in [28]. For , and for , P after 2.3 in “standard rock”, i.e. in the extreme case where the ’s are produced just after entering the volcano. However, the values of at lower energies cannot be obtained by a simple extrapolation since the value of below 100 is dominated by the mean life time of the . By considering this fact, we obtain approximate values of for the energy regime. Table 1 shows the average values of in the energy range of interest of this work. Finally, the factor makes sure that the charged lepton escapes the volume of the volcano with an energy . Based on Equation 11, taking , then can be written as
(13) 
Thus, we can approximate the function as
(14) 
After evaluating the integral we obtain an expression for the function.
(15) 
Equation 15 is valid in the case where we neglect variations of the density of the volcano. Going back to Equation 5, we can substitute Equation 6 and Equation 15 on it. By doing this, we obtain an equation for the differential flux with respect to the energy of the leptons
(16) 
The next step is to define a function that describes the differential flux of neutrinos as a function of energy. For this, we use the parametrization of the measured astrophysical neutrino flux made by IceCube [30, 31, 32], which has been found to be well described by an unbroken power law
(17) 
We choose to use the parametrization presented in [31] for the measurement of the astrophysical flux. For this particular parametrization we shifted the mean values of the normalization and spectral index in order to obtain the highest flux, within the allowed range given by their statistical uncertainties. These values are: and . We assume for our calculations an expected equal contribution to the astrophysical flux for , due to the neutrino flavour mixing over cosmological distances that would produce approximately equal proportions of all flavours. Then, we extrapolate the measured flux to the energy range 10 to 100 . Moreover, motivated by the most energetic event found by the diffuse flux muon neutrino search (20092015) done by IceCube (a track event that deposited in the sensible volume of the detector [33]), we also include in the flux estimations the models proposed in [34], which account for sources of multi neutrinos that are constrained by the most recent ultrahigh energy neutrino upper limits set by IceCube [35] and Pierre Auger [11]. These models, for the sum of all neutrino flavours, have the smoothlybroken power law functional forms
(18) 
The values of the parameters are: , and . take the values of: for what we refer in this work as Model A, for Model B and for Model C. The values of the normalizations for each of these models are presented in Table 2
Model  

A  
B  
C 
According to [34], Model A & B would correspond to the spectra produced by BL Lac AGNs, and combinations of Models A & C would follow the expected shape of GZK neutrinos produced from EBL and CMB interactions.
In order to calculate the flux of charged leptons produced by neutrinonucleon interactions we need to define the input parameters that enter the calculation of the number of produced charged leptons (Equation 16). For the value of the parameter that describes the radiative energy losses of the charged leptons as they travel trough the rock (). One can find in the literature parametrizations that may differ by up to a factor of two. For instance, it ranges from 0.26 to 0.59 at 100 in [36]. For our calculations we decided to use the values obtained with the ASW structure functions calculated in [36]; and for the results obtained in [20]. For this latter case the results are available up to an energy of 1 , but these values are enough for our calculations. For the numerical integrations, we used intermediate values of and in different bins of energy as shown in Table 3
evaluated at []  Energy Bin []  []  [] 

50  [, ]  2.496  4.960 
500  [, ]  2.987  5.143 
5  [, ]  3.554  N/A 
5  [, ]  4.184  N/A 
For the density of the volcano, we take that of “standard rock”, [20]. For the calculation of the charged current interaction length we use the results from [37], with the cross section for neutrinonucleon interactions via charged currents in the energy range given by
(19) 
This parametrization is a well known result. However, as a cross check, we compared the parametrization given by [37] with more recent calculations [38, 39]. The result of this is presented in Figure 1. One can notice that the result from [37] follows the general trend predicted by the work of [38] and Sarkar et. al. Due to this, and the fact that there is an analytical expression for the cross section of [37], we decided to use Equation 19 in our numerical calculations. The anti neutrinonucleon cross section is taken to be the same [37]. Using these results we get a value of , much smaller than the lightyear interaction length for average neutrinos in lead [40].
For the average width of the Pico de Orizaba, we calculated the typical length of a chord that goes through a cone that follows the geometry of the volcano, as it will be shown in the following section of the paper. At 4100 a.s.l., the volcano has a width of 6 and at the summit at 5 500 a.s.l., a width of 0.35 (see Figure 5 and Figure 5). The average path length for the neutrinos is of 2.33 inside the volcano. Finally, we assume a duty cycle =95%, which agrees with the reported value observed during actual HAWC operations [41]. The energy range for the charged leptons that exit the volcano is taken in the range (, ). In this way we consider an energy range that contains a good fraction of the charged leptons that exit the volcano. This is because the charged leptons escape with an energy distribution, after loosing some of their energy in the medium. For a study of the propagation of a monoenergetic beam of charged leptons through rock see for instance the work from [28].
We integrate numerically Equation 16 using the differential neutrinos fluxes from Equation 17 and 18 to obtain the number of charged leptons () that escape the mountain. The results are presented in Tables 4 and 5.
Neutrino Energy []  Astrophysical  Model A  Model B  Model C 

10  100  10.7471  0.4732  0.8269  0.1008 
100  10  11.0342  0.0233  0.4251  0.3690 
Neutrino Energy []  Astrophysical  Model A  Model B  Model C 

10  100  1.0816  0.0476  0.0832  0.0101 
100  1  0.4685  0.0015  0.0267  0.0185 
One can notice that the flux of charged leptons is more than an order of magnitude lower than that of ’s. The reason for this can be seen in Equation 16, where the charged lepton flux depends inversely on the radiative energy losses of the charged leptons, which are more than an order of magnitude higher for the ’s with respect to ’s (see Table 3). This does not contradict the fact that, at the highest energies, for both ’s and ’s the survival probability . quantifies the probability that the charged leptons are able to escape the volcano, while the factor in Equation 16 appears because of the form of and is independent of . This result is consistent with the arguments developed in [42], where it is pointed out that charged leptons have a higher probability to escape the Earth crust compared to ’s.
In the following section we present a simple method that allows to approximate the effective area of the HAWC observatory to ultrahigh energy charged leptons.
3 Effective area determination
In general, the effective area or acceptance of an observatory is calculated using a detailed Monte Carlo simulation, see for instance in [43, 44] studies of the sensitivity of air Cherenkov and fluorescence observatories to neutrinos using the Earthskimming technique and in [18] an early proposal about using fluorescence measurements on observatories as Pierre Auger in order to search for charged lepton showers. The full Monte Carlo method of calculating the effective area represents a complex and time consuming task. In the early results obtained by the Pierre Auger collaboration, their effective area calculation neglected the effect of the topography that surrounds the observatory, and took into account its effect in the systematic error of their observation limits [13]. However, their most recent results have incorporated the topography into their effective area calculations in two of their three analysis channels [11]. In our case we are specifically interested on the effect of the largest volcano that surrounds the HAWC observatory as a target for the neutrinonucleon interactions. In order to reproduce the topography that surrounds HAWC, we use data from the Instituto Nacional de Estadística y Geografía (INEGI, México) [45]. Figure 2 shows the topography that surrounds the HAWC site (indicated by a small red rectangle). The observatory is located between two volcanoes: in the direction NorthingEasting by the Pico de Orizaba and in the opposite direction by the Sierra Negra volcano. Due to its much larger volume, we focus our attention on the Pico de Orizaba volcano as the target for the neutrinonucleon interactions.
Figure 3 shows the topography of the Pico de Orizaba volcano as seen in a coordinate system centred at the location of the HAWC array at 4090 a.s.l. The approximate symmetry of the Pico de Orizaba volcano makes it easy to motivate a simplification in the calculation of the effective area. We can assume that the geometry of the mountain can be modelled to be conical. Figure 5 and Figure 5 show the profile of the volcano along the Easting and Northing directions, and how we can approximate this profile using a cone with a height of 1500 m ( taken at the altitude above sea level of the HAWC observatory) and a diameter of 6000 m.
Our method to calculate the effective area to detect ultrahigh energy charged leptons produced by neutrinonucleon interaction is as follows:

We consider that the ultrahigh energy charged lepton, or their highly boosted decay products, will travel following a straight trajectory along the direction of the initial neutrino.

We approximate the region where the charged leptons are produced as a triangular surface, a 2D projection of the volcano, that faces the HAWC observatory. We can use this simplification since here we are only interested on studying the trajectories of the produced charged leptons. The width of the mountain is considered in the charged lepton flux calculation, as presented in Section 2. On the triangular surface we draw a grid made of lines parallel to the & axis in the coordinate system shown in Figure 7. The lines of the grid have a separation of 0.5 in each direction. The surface is located at a distance of 5.7 to the center of the HAWC observatory, approximately the distance between the center of the detector array and that of the volcano.

The detection volume of HAWC is modelled as a rectangular prism, with dimensions () 140 140 4.5 . This approximation is motivated by the actual configuration of the HAWC array. The observatory is made of a compact group of 300 Water Cherenkov Detectors (WCDs). Each WCD is a cylinder having 7.3 in diameter and 4.5 in height. The array covers an area of approximately 22 000 [41].

From each cell of the grid ( million), we generate vectors that point towards our model of HAWC, approximated as a rectangular prism. Figure 7 shows a diagram of the definition of the angles used in this work. The azimuth angle covers the range from to (with steps of 1 degree), pointing towards the  direction in the coordinate system shown in Figure 7. This is done for the different orientations of the polar angle . The polar angle is defined such that the vector orientation that lies in the plane corresponds to 90, and increases as the orientation gets closer to the plane (see Figure 7).

For each orientation , we calculate a differential element of effective area
(20) where, is the total number of vectors whose directions points towards the detection volume of HAWC with a minimum length of the trajectory across the plane shown in Figure 7 (inside the volume of the rectangular prism). We will refer to this minimum length as the trigger condition, that approximately represents how many WCDs would measure Cherenkov light produced by the incoming lepton or their collimated decay products. is the total number of vectors that point towards the rectangular prism, regardless of the trajectory that they have inside the prism. is the section of the area of the triangular surface over which the vectors that point towards the rectangular prism () were generated. Finally, is the azimuthal angle covered by the area , as seen from the center of the rectangular prism.
By following the procedure described above, we can obtain the differential elements of effective area as a function of the polar angle orientation. Figure 8 shows the results for different orientations of the polar angle.
One can see that the differential effective area increases with increasing values of the polar angle, up to a maximum at 93, and then decreases because the area over which the vectors are generated decreases as well, as it corresponds to the upper sections of the triangular grid where the vectors are generated. The different line styles used in the plot show the results obtained for different trigger conditions. For instance, the curve with 14 as trigger condition indicates the results obtained when the requirement of the path of the charged lepton through the detection volume corresponds to at least a trajectory in the plane inside the HAWC rectangular prism of 14 . This would correspond roughly to have the ultrahigh energy charged lepton (or its boosted decay products) to traverse through at least two WCDs. This may not seem enough to even provide a rough idea of the direction of the incoming lepton. However, one should keep in mind that each WCD is equipped with four photomultiplier tubes (PMTs), recording the amplitude and arrival time of the Cherenkov light. By combining the information from the eight PMTs from two WCDs it may be possible to have enough information to reconstruct a trajectory, although details of the angular reconstruction of the tracks is beyond the scope of this work. In Figure 8, we present the results for trigger conditions that correspond roughly to have the charged lepton (or its collimated decay products) traversing from two up to six WCDs. As expected, the effective area that corresponds to a given trigger condition decreases as one increases the number of detectors required to have signals. We are interested in the total effective area that goes into the calculation of the number of possible detections of Equation 1. This total effective area is given by
(21) 
where the integration range in goes from 90 to 105 (the angular height of the Pico de Orizaba volcano as seen from the center of the HAWC observatory is of 15). Our results of the effective area calculation as a function of the different trigger conditions are shown in Table 6.
Trigger []  [] 

14  0.013867 
21  0.012544 
28  0.010578 
35  0.008418 
42  0.006782 
With the results from Sections 2 and 3 we have: the flux of charged leptons given a certain isotropic flux of neutrinos and the effective area of the HAWC observatory obtained with a simple trigger condition. Using this, in the next section we present our results for the flux of ultrahigh energy charged leptons that could be detected by the HAWC observatory.
4 Results, possible background signals and discussion
Table 7 shows the number of detectable charged ’s produced by neutrinonucleon interactions in the Pico de Orizaba volcano, obtained using the effective area results presented in Table 6 and the ultrahigh energy lepton fluxes from Table 4. We restrict the results to those from charged leptons since, as it was shown in Section 2, the flux from charged leptons is approximately an order of magnitude lower, making their detection not feasible with the method proposed in this work. The detection estimates are shown as a function of the different trigger conditions and of the different ultrahigh energy neutrino fluxes discussed in Section 2. The eight column presents an estimate of the background signals that are expected using this method. We estimated the background considering that the main contribution is that of ultrahigh energy atmospheric muons, coming from the direction of the Pico de Orizaba. In order to calculate this, we used the characterization of the atmospheric muon flux above 15 measured by IceCube [46], that can be modelled as an unbroken power law:
(22) 
In this case, we shifted the mean values (within the statistical and systematic uncertainties of the IceCube measurement) of the normalization and spectral index in order to obtain the worst case scenario of the background, and then extrapolated the muon flux to the energy range of our interest. The values of the parameters that maximize the background muon flux are: and a spectral index of 3.73. Thus, the background flux can be calculated using Equation 1, where in this case the charged lepton flux corresponds to the integral of Equation 22 within the energy bins used in this work. Of course this is an overestimation, since these atmospheric muons would have to be able to survive their travel through the volcano, however we consider this estimation of the background appropriate because the study of the propagation of atmospheric muons through the volcano is beyond the scope of this work. Moreover, it has already been noticed that the atmospheric muon background is relevant up to energies of 100 [47]. The last two columns of Table 7 show the estimated time to get a detection for each trigger condition, first for the case of only having as a source the neutrino flux given by the extrapolation of the measurement done by IceCube [30], and also for the case of including GZK neutrinos from [34].
From Table 7, we find that it would take around eight years of data taking in order to be able to find a tracklike signal traversing three WCDs based on the extrapolation of the measured astrophysical flux, for neutrinos in the energy range from 10 to 100 . In order to find a tracklike signal that propagates through four WCDs it could take more than nine years of data taking, but such a long track may allow to reconstruct with better accuracy the direction of the original neutrino. We have verified the predicted detection rate by doing a back of the envelope comparison with the neutrino measurements reported by IceCube [1], and found that, taking into account the different experimental conditions, our results fall in the expected order of magnitude.
The HAWC observatory is currently planned to operate for at least five years [48], but given its importance as a trigger for future instruments such as CTA, it is not unlikely that it may operate for at least 10 years. Table 7 shows also for completeness the results expected for neutrinos in the energy range from 100 to 100 . Charged leptons that do not decay before arriving or inside the HAWC array (energies of 100 or larger) will leave almost the same signal because the values of (see Table 3) change very slowly with energy [8]. This condition will make very difficult to estimate the energy of the incoming . The results also show that the background from ultrahigh energy muons at the energy regime of interest of this work is very low and would allow a clean identification of a signal coming from the direction of the Pico de Orizaba volcano.
Neutrino energy []  Trigger []  :IC  :  :  :  :IC  Bckg  T IC []  T IC [] 

10  100  14  0.1416  0.0063  0.0109  0.0014  0.1492  0.0145  7.1  6.7 
21  0.1281  0.0057  0.0099  0.0012  0.1349  0.0131  7.8  7.4  
28  0.1080  0.0048  0.0083  0.0010  0.1138  0.0111  9.3  8.8  
35  0.0860  0.0038  0.0066  0.0008  0.0906  0.0088  11.6  11.0  
42  0.0693  0.0031  0.0054  0.0007  0.0730  0.0071  14.4  13.7  
100  100  14  0.1454  0.0003  0.0056  0.0049  0.1505  6.9  6.6  
21  0.1315  0.0003  0.0050  0.0044  0.1321  7.6  7.6  
28  0.1109  0.0003  0.0042  0.0037  0.1148  9.0  8.7  
35  0.0882  0.0002  0.0034  0.0030  0.0914  11.3  10.9  
42  0.0711  0.0002  0.0027  0.0023  0.0736  14.1  13.6 
5 Triggering of signals and background rejection
The current software trigger used by the HAWC observatory consists on requiring a certain number of PMTs to be above a charge threshold (it actually consist on two thresholds used to estimate the charge with the TimeoverThreshold method; the thresholds are of 0.25 and 4 ) within a time sliding window of the order of 150 [21]. For the HAWC results presented in [41], the number of PMTs required to trigger an event was of 15. Thus, the multiplicity trigger used by the HAWC Collaboration is already useful for the collection of signals that are needed for the neutrino searches proposed in this work. To show this, we can estimate the average number of PMTs that have signals during any 150 time window. From [49] we know that the sum of the signals from 112 HAWC’s PMTs in a 60 time window was of 2.285. Then, we can use this number to estimate the number of PMTs fired in any 150 triggering window, which is of 6. Thus, only 9 additional PMTs would be required to have signals within the specified time window in order to produce a trigger useful for the neutrino search proposed in this work. This condition can be easily fulfilled if an ultrahigh energy lepton or its decay products pass through 3 WCDs, firing at least nine of its 12 PMTs (with a propagation time of the corresponding signals along 3 WCDs of 70 , well within the current trigger window). Of course, a specialized trigger can be developed in order to obtain a better trigger efficiency. This could be done, for instance, by selecting events with topologies consistent with tracks propagating through the WCDs. This would also have the consequence of reducing the bandwidth needed for this additional trigger.
A second point of interest to discuss is the ability of the instrument to separate the signals produced by ultrahigh energy charged leptons from those from the background of lower energy atmospheric muons. The amount of emitted Cherenkov photons per unit length inside a HAWC WCD is given by the FrankTamm equation; and the total amount of light will be proportional to the the sum of the path lengths of the primary and all the produced secondary charged particles (from bremsstrahlung and pair production processes) that travel inside a WCD and have a velocity larger than the threshold for Cherenkov light production. According to the study presented in [50] this sum of path lengths increases linearly with the energy of the primary particle. Thus, the way to discriminate background signals could be based on the total light yield detected in the WCDs that belong to a given track (e.g. a 1 primary electron would produce an order of magnitude more photons than those from a 100 primary electron). Based on the results presented in [51], the mean momentum of nearly horizontal tracks ( = 75 ) is of 200 /c, and in this work we aim to detect the tracks produced by leptons of energies, so we could expect a factor larger than in the deposited light in the WCDs for the ultrahigh energy neutrino initiated signals relative to the average muon background tracks. Moreover, the results presented in [52] show that the integral intensity of muons of 1 /c is larger than that of 200 /c muons by a factor of 52, making reasonable to easily handle the data rate if a specialized trigger is set to keep the data of tracks related to large deposits of Cherenkov light.
An additional experimental proof that this analysis could be performed, is the ability to observe the cosmic ray shadow produced by the volcano, as for instance pointed out in [18]. However, this requires the use of data and algorithms property of the HAWC Collaboration that are no publicly available. Nonetheless we are confident, because of the interest on this result, that the Pico de Orizaba cosmic ray shadow will be publicly available soon.
Due to this facts, we believe that it is feasible to separate the huge atmospheric background from the signals we are looking for. We are aware that the precision that could be achieve in determining the energy of the primary tau lepton will be low, nonetheless the primary interest of this work is the detection technique and we leave the development of a possible energy estimator for future work.
6 Conclusions
We presented an estimate of the detection capabilities of the HAWC observatory to study ultrahigh energy neutrinos interacting in the Pico de Orizaba volcano. We based our study in the analytic method developed by [7, 26] and modified it to a simpler case where the neutrino conversion takes place inside a medium of constant density. The effective area of the HAWC observatory for the detection of ultrahigh energy charged leptons was calculated geometrically for different trigger conditions. We used the astrophysical neutrino flux measured by IceCube [30] and extrapolated it to energy range from 10 up to 100 , and also considered models for multi neutrinos [34] that are constrained by the most recent data from both IceCube [35] and the Pierre Auger observatory [11]. With this, we found that in order to find a signal consistent with the propagation of an ultrahigh energy charged lepton coming from the direction of the Pico de Orizaba volcano, will require approximately nine years of data taking if the signal is required to propagate in four WCDs of HAWC. We estimated the expected background for this analysis using the measured atmospheric muon flux above 15 [46], and found that it should be feasible to perform this study with a reasonable signal to background ratio during the lifetime of the HAWC observatory. We also showed that the current software trigger used by HAWC should be sufficient to acquire data for this analysis, and that given the dynamic range of the HAWC electronics it is feasible to be able to discriminate the background signals produced by lower energy muons using the light yield detected by the HAWC PMTs.
As an anonymous referee pointed out to us, the detection rate should be taken with caution since there are several quantities that are uncertain to at least some degree, such as the flux of neutrinos at the ultrahigh energy regime, the neutrinonucleon cross section, and the parameters that describes the radiative energy loss processes for the leptons traversing the Earth crust. By selecting a different set of the input parameters for the calculations one could find a detection rate that could decrease by a factor of three or even more. However we also did not consider scenarios where the neutrinonucleon cross section at ultrahigh energies could have enhanced values relative to the standard model predictions, due to the presence of new physics, as described for instance by [53]. Further work should be done using a complete Monte Carlo simulation that incorporates the detailed topography of the volcanoes that surround the HAWC observatory and the complete detector simulation to study the detector response. However this first step shows encouraging results to pursue more detailed studies. The detection rate is certainly low, not comparable to the one that dedicated neutrino experiment can achieve, but as pointed out in [34], a single detection of a neutrino with an energy higher than 10 , would give evidence of a flux beyond what is firmly established. Our results indicate that, although extremely challenging, it is worth trying to detect ultrahigh energy neutrinos, interacting in the Pico de Orizaba with the HAWC observatory.
Competing Interests
The authors declare that there is no conflict of interest regarding the publication of this article.
Acknowledgments
The authors thank an anonymous referee for very useful discussions that allowed to greatly improve the quality of this work. This work was supported by the programme UNAMDGAPAPAPIIT IA102715, Consejo Nacional de Ciencia y Tecnología (CONACyT), Mexico grant 254964 and Coordinación de la Investigación Científica UNAM.
References
 IceCube Collaboration [2013] IceCube Collaboration. Evidence for HighEnergy Extraterrestrial Neutrinos at the IceCube Detector. Science, 342:1242856, November 2013.
 Abeysekara et al. [2013] A. U. Abeysekara, R. Alfaro, C. Alvarez, J. D. Álvarez, R. Arceo, J. C. ArteagaVelázquez, H. A. Ayala Solares, A. S. Barber, B. M. Baughman, N. BautistaElivar, E. Belmont, S. Y. BenZvi, D. Berley, M. Bonilla Rosales, J. Braun, R. A. CaballeroLopez, A. Carramiñana, M. Castillo, U. Cotti, J. Cotzomi, E. de la Fuente, C. De León, T. DeYoung, R. Diaz Hernandez, J. C. DiazVelez, B. L. Dingus, M. A. DuVernois, R. W. Ellsworth, A. Fernandez, D. W. Fiorino, N. Fraija, A. Galindo, J. L. GarciaLuna, G. GarciaTorales, F. Garfias, L. X. González, M. M. González, J. A. Goodman, V. Grabski, M. Gussert, Z. HampelArias, C. M. Hui, P. Hüntemeyer, A. Imran, A. Iriarte, P. Karn, D. Kieda, G. J. Kunde, A. Lara, R. J. Lauer, W. H. Lee, D. Lennarz, H. León Vargas, E. C. Linares, J. T. Linnemann, M. Longo, R. LunaGarcía, A. Marinelli, O. Martinez, J. MartínezCastro, J. A. J. Matthews, P. MirandaRomagnoli, E. Moreno, M. Mostafá, J. Nava, L. Nellen, M. Newbold, R. NoriegaPapaqui, T. OcegueraBecerra, B. Patricelli, R. Pelayo, E. G. PérezPérez, J. Pretz, C. Rivière, J. Ryan, D. RosaGonzález, H. Salazar, F. Salesa, A. Sandoval, E. Santos, M. Schneider, S. Silich, G. Sinnis, A. J. Smith, K. Sparks, R. W. Springer, I. Taboada, P. A. Toale, K. Tollefson, I. Torres, T. N. Ukwatta, L. Villaseñor, T. Weisgarber, S. Westerhoff, I. G. Wisher, J. Wood, G. B. Yodh, P. W. Younk, D. Zaborov, A. Zepeda, and H. Zhou. Sensitivity of the high altitude water Cherenkov detector to sources of multiTeV gamma rays. Astroparticle Physics, 50:26–32, December 2013.
 Gould and Schréder [1966] Robert J. Gould and Gerald Schréder. Opacity of the universe to highenergy photons. Phys. Rev. Lett., 16:252–254, Feb 1966.
 Fargion et al. [1999] D. Fargion, Aiello A., and Conversano R. Horizontal Tau air showers from mountains in deep vally :Traces of Ultrahigh neutrino tau. International Cosmic Ray Conference, 2:396, August 1999.
 Fargion [2002] D. Fargion. Discovering ultrahighenergy neutrinos through horizontal and upward Ï air showers: Evidence in terrestrial gamma flashes? The Astrophysical Journal, 570(2):909, 2002.
 LetessierSelvon [2001] A. LetessierSelvon. Establishing the gzk cutoff with ultra high energy tau neutrinos. In AIP Conference Proceedings, pages 157–171, 2001.
 Feng et al. [2002] J. L. Feng, P. Fisher, F. Wilczek, and T. M. Yu. Observability of EarthSkimming Ultrahigh Energy Neutrinos. Physical Review Letters, 88(16):161102, April 2002.
 Fargion et al. [2004] D. Fargion, P. G. De Sanctis Lucentini, M. De Santis, and M. Grossi. Tau Air Showers from Earth. Astrophysical Journal, 613:1285–1301, October 2004.
 Conrad et al. [2010] J. Conrad, A. de Gouvêa, S. Shalgar, and J. Spitz. Atmospheric tau neutrinos in a multikiloton liquid argon detector. Physical Review D, 82(9):093012, November 2010.
 Aartsen et al. [2016a] M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, I. Ansseau, and et al. Search for astrophysical tau neutrinos in three years of IceCube data. Physical Review D, 93(2):022001, January 2016a.
 Aab et al. [2015] A. Aab, P. Abreu, M. Aglietta, E. J. Ahn, I. Al Samarai, I. F. M. Albuquerque, I. Allekotte, P. Allison, A. Almela, J. Alvarez Castillo, and et al. Improved limit to the diffuse flux of ultrahigh energy neutrinos from the Pierre Auger Observatory. Physical Review D, 91(9):092008, May 2015.
 Abraham et al. [2008] J. Abraham, P. Abreu, M. Aglietta, C. Aguirre, D. Allard, I. Allekotte, J. Allen, P. Allison, J. AlvarezMuñiz, M. Ambrosio, and et al. Upper Limit on the Diffuse Flux of Ultrahigh Energy Tau Neutrinos from the Pierre Auger Observatory. Physical Review Letters, 100(21):211101, May 2008.
 Abraham et al. [2009] J. Abraham, P. Abreu, M. Aglietta, C. Aguirre, E. J. Ahn, D. Allard, I. Allekotte, J. Allen, P. Allison, J. AlvarezMuñiz, and et al. Limit on the diffuse flux of ultrahigh energy tau neutrinos with the surface detector of the Pierre Auger Observatory. Physical Review D, 79(10):102001, May 2009.
 Miele et al. [2006] G. Miele, S. Pastor, and O. Pisanti. The aperture for UHE tau neutrinos of the Auger fluorescence detector using a Digital Elevation Map. Physics Letters B, 634:137–142, March 2006.
 Aita et al. [2011] Y. Aita, T. Aoki, Y. Asaoka, T. Chonan, M. Jobashi, M. Masuda, Y. Morimoto, K. Noda, M. Sasaki, J. Asoh, N. Ishikawa, S. Ogawa, J. G. Learned, S. Matsuno, S. Olsen, P.M. Binder, J. Hamilton, N. Sugiyama, Y. Watanabe, and Ashra1 Collaboration. Observational Search for PeVEeV Tau Neutrino from GRB081203A. Astrophysical Journal, Letters, 736:L12, July 2011.
 Gaug et al. [2008] M. Gaug, C. Hsu, J. K. Becker, A. Biland, M. Mariotti, W. Rhode, and M. Teshima. Tau neutrino search with the MAGIC telescope. volume 3, pages 1273–1276, 2008.
 Becker et al. [2008] Julia K. Becker, Markus Gaug, ChingâCheng Hsu, Wolfgang Rhode, and MAGIC Collaboration. Grb neutrino search with magic. volume 1000, pages 245–248, 2008.
 Fargion et al. [2008a] D. Fargion, P. Oliva, F. Massa, and G. Moreno. Cherenkov flashes and fluorescence flares on telescopes: New lights on UHECR spectroscopy while unveiling neutrinos astronomy. Nuclear Instruments and Methods in Physics Research A, 588:146–150, April 2008a.
 Fargion et al. [2008b] D. Fargion, M. Gaug, and P. Oliva. Reflecting on Čerenkov reflections. In Journal of Physics Conference Series, volume 110 of Journal of Physics Conference Series, page 062008, May 2008b.
 Lipari and Stanev [1991] P. Lipari and T. Stanev. Propagation of multiTeV muons. Physical Review D, 44:3543–3554, December 1991.
 Smith for the HAWC Collaboration [2015] Smith for the HAWC Collaboration. HAWC: Design, Operation, Reconstruction and Analysis. August 2015.
 Abeysekara et al. [2017] A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, J. D. Álvarez, R. Arceo, J. C. ArteagaVelázquez, H. A. Ayala Solares, A. S. Barber, N. BautistaElivar, A. Becerril, E. BelmontMoreno, S. Y. BenZvi, D. Berley, J. Braun, C. Brisbois, K. S. CaballeroMora, T. Capistrán, A. Carramiñana, S. Casanova, M. Castillo, U. Cotti, J. Cotzomi, S. Coutiño de León, E. de la Fuente, C. De León, T. DeYoung, B. L. Dingus, M. A. DuVernois, J. C. DíazVélez, R. W. Ellsworth, D. W. Fiorino, N. Fraija, J. A. GarcíaGonzález, M. Gerhardt, A. González Munoz, M. M. González, J. A. Goodman, Z. HampelArias, J. P. Harding, S. Hernandez, A. HernandezAlmada, J. Hinton, C. M. Hui, P. Hüntemeyer, A. Iriarte, A. JardinBlicq, V. Joshi, S. Kaufmann, D. Kieda, A. Lara, R. J. Lauer, W. H. Lee, D. Lennarz, H. León Vargas, J. T. Linnemann, A. L. Longinotti, G. L. Raya, R. LunaGarcía, R. LópezCoto, K. Malone, S. S. Marinelli, O. Martinez, I. MartinezCastellanos, J. MartínezCastro, H. MartínezHuerta, J. A. Matthews, P. MirandaRomagnoli, E. Moreno, M. Mostafá, L. Nellen, M. Newbold, M. U. Nisa, R. NoriegaPapaqui, R. Pelayo, J. Pretz, E. G. PérezPérez, Z. Ren, C. D. Rho, C. Rivière, D. RosaGonzález, M. Rosenberg, E. RuizVelasco, H. Salazar, F. Salesa Greus, A. Sandoval, M. Schneider, H. Schoorlemmer, G. Sinnis, A. J. Smith, R. W. Springer, P. Surajbali, I. Taboada, O. Tibolla, K. Tollefson, I. Torres, T. N. Ukwatta, L. Villaseñor, T. Weisgarber, S. Westerhoff, I. G. Wisher, J. Wood, T. Yapici, G. B. Yodh, P. W. Younk, A. Zepeda, and H. Zhou. Observation of the Crab Nebula with the HAWC GammaRay Observatory. ArXiv eprints, January 2017.
 Ayala Solares et al. [2015] H. A. Ayala Solares, M. Gerhardt, C. M. Hui, R. J. Lauer, Z. Ren, F. Salesa Greus, H. Zhou, and for the HAWC Collaboration. The Calibration System of the HAWC GammaRay Observatory. August 2015.
 Agostinelli et al. [2003] S. Agostinelli et al. GEANT4: A Simulation toolkit. Nucl. Instrum. Meth., A506:250–303, 2003. doi: 10.1016/S01689002(03)013688.
 Atağ and Gürkanlı [2016] S. Atağ and E. Gürkanlı. Prediction for CP violation via electric dipole moment of lepton in process at CLIC. Journal of High Energy Physics, 6:118, June 2016.
 Feng et al. [2001] J. L. Feng, P. Fisher, F. Wilczek, and T. M. Yu. Observability of EarthSkimming Ultrahigh Energy Neutrinos. preprint, hepph/0105067v1, 2001.
 Dutta et al. [2005] S. I. Dutta, Y. Huang, and M. H. Reno. Tau neutrino propagation and tau energy loss. Physical Review D, 72(1):013005, July 2005.
 Blanch Bigas et al. [2008] O. Blanch Bigas, O. Deligny, K. Payet, and V. van Elewyck. Tau energy losses at ultrahigh energy: Continuous versus stochastic treatment. Physical Review D, 77(10):103004, May 2008.
 Iyer Dutta et al. [2001] S. Iyer Dutta, M. H. Reno, I. Sarcevic, and D. Seckel. Propagation of muons and taus at high energies. Physical Review D, 63(9):094020, May 2001.
 The IceCube Collaboration et al. [2015a] The IceCube Collaboration, M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, and et al. The IceCube Neutrino Observatory  Contributions to ICRC 2015 Part II: Atmospheric and Astrophysical Diffuse Neutrino Searches of All Flavors: Combined Analysis of the HighEnergy Cosmic Neutrino Flux at the IceCube Detector. October 2015a.
 The IceCube Collaboration et al. [2015b] The IceCube Collaboration, M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, and et al. The IceCube Neutrino Observatory  Contributions to ICRC 2015 Part II: Atmospheric and Astrophysical Diffuse Neutrino Searches of All Flavors: A measurement of the diffuse astrophysical muon neutrino flux using multiple years of IceCube data. October 2015b.
 Aartsen et al. [2015] M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, M. Archinger, and et al. A Combined Maximumlikelihood Analysis of the Highenergy Astrophysical Neutrino Flux Measured with IceCube. Astrophysical Journal, 809:98, August 2015.
 Schoenen and Raedel [2015] S. Schoenen and L. Raedel. Detection of a multiPeV neutrinoinduced muon event from the Northern sky with IceCube. The Astronomer’s Telegram, 7856, July 2015.
 Kistler and Laha [2016] M. D. Kistler and R. Laha. MultiPeV Signals from a New Astrophysical Neutrino Flux Beyond the Glashow Resonance. May 2016.
 The IceCube Collaboration et al. [2015c] The IceCube Collaboration, M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, and et al. The IceCube Neutrino Observatory  Contributions to ICRC 2015 Part II: Atmospheric and Astrophysical Diffuse Neutrino Searches of All Flavors: Observation of Astrophysical Neutrinos in Four Years of IceCube Data. October 2015c.
 Armesto et al. [2008] N. Armesto, C. Merino, G. Parente, and E. Zas. Charged current neutrino cross section and tau energy loss at ultrahigh energies. Physical Review D, 77(1):013001, January 2008.
 Gandhi et al. [1998] R. Gandhi, C. Quigg, M. H. Reno, and I. Sarcevic. Neutrino interactions at ultrahigh energies. Physical Review D, 58(9):093009, November 1998.
 Connolly et al. [2011] Amy Connolly, Robert S. Thorne, and David Waters. Calculation of high energy neutrinonucleon cross sections and uncertainties using the martinstirlingthornewatt parton distribution functions and implications for future experiments. Physical Review D, 83:113009, Jun 2011.
 Argüelles et al. [2015] C. A. Argüelles, F. Halzen, L. Wille, M. Kroll, and M. H. Reno. Highenergy behavior of photon, neutrino, and proton cross sections. Physical Review D, 92(7):074040, October 2015.
 Cho [2013] Adrian Cho. Physicists snare a precious few neutrinos from the cosmos. Science, 342(6161):920–920, 2013. ISSN 00368075.
 Abeysekara et al. [2014] A. U. Abeysekara, R. Alfaro, C. Alvarez, J. D. Álvarez, R. Arceo, J. C. ArteagaVelázquez, H. A. Ayala Solares, A. S. Barber, B. M. Baughman, N. BautistaElivar, E. Belmont, S. Y. BenZvi, D. Berley, M. Bonilla Rosales, J. Braun, K. S. CaballeroMora, A. Carramiñana, M. Castillo, U. Cotti, J. Cotzomi, E. de la Fuente, C. De León, T. DeYoung, R. Diaz Hernandez, J. C. DíazVélez, B. L. Dingus, M. A. DuVernois, R. W. Ellsworth, D. W. Fiorino, N. Fraija, A. Galindo, F. Garfias, M. M. González, J. A. Goodman, M. Gussert, Z. HampelArias, J. P. Harding, P. Hüntemeyer, C. M. Hui, A. Imran, A. Iriarte, P. Karn, D. Kieda, G. J. Kunde, A. Lara, R. J. Lauer, W. H. Lee, D. Lennarz, H. León Vargas, J. T. Linnemann, M. Longo, R. LunaGarcía, K. Malone, A. Marinelli, S. S. Marinelli, H. Martinez, O. Martinez, J. MartínezCastro, J. A. J. Matthews, J. McEnery, E. Mendoza Torres, P. MirandaRomagnoli, E. Moreno, M. Mostafá, L. Nellen, M. Newbold, R. NoriegaPapaqui, T. OcegueraBecerra, B. Patricelli, R. Pelayo, E. G. PérezPérez, J. Pretz, C. Rivière, D. RosaGonzález, E. RuizVelasco, J. Ryan, H. Salazar, F. Salesa Greus, A. Sandoval, M. Schneider, G. Sinnis, A. J. Smith, K. Sparks Woodle, R. W. Springer, I. Taboada, P. A. Toale, K. Tollefson, I. Torres, T. N. Ukwatta, L. Villaseñor, T. Weisgarber, S. Westerhoff, I. G. Wisher, J. Wood, G. B. Yodh, P. W. Younk, D. Zaborov, A. Zepeda, H. Zhou, and HAWC Collaboration. Observation of Smallscale Anisotropy in the Arrival Direction Distribution of TeV Cosmic Rays with HAWC. Astrophysical Journal, 796:108, December 2014.
 Kusenko and Weiler [2002] A. Kusenko and T. J. Weiler. Neutrino Cross Sections and Future Observations of UltrahighEnergy Cosmic Rays. Physical Review Letters, 88(16):161101, April 2002.
 Góra et al. [2015] D. Góra, E. Bernardini, and A. Kappes. Searching for tau neutrinos with Cherenkov telescopes. Astroparticle Physics, 61:12–16, February 2015.
 Cao et al. [2005] Z Cao, M A Huang, P Sokolsky, and Y Hu. Ultra high energy Î½ Ï detection with a cosmic ray tau neutrino telescope using fluorescence/cerenkov light technique. Journal of Physics G: Nuclear and Particle Physics, 31(7):571, 2005.
 INEGI [2012] INEGI. Continuo de elevaciones mexicano 3.0 (cem 3.0), 2012. URL http://www.inegi.org.mx/geo/contenidos/datosrelieve/continental/continuoelevaciones.aspx. Cited 25 Dec 1999.
 Aartsen et al. [2016b] M. G. Aartsen, K. Abraham, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, M. Archinger, and et al. Characterization of the atmospheric muon flux in IceCube. Astroparticle Physics, 78:1–27, May 2016b.
 Anchordoqui et al. [2014] L. A. Anchordoqui, V. Barger, I. Cholis, H. Goldberg, D. Hooper, A. Kusenko, J. G. Learned, D. Marfatia, S. Pakvasa, T. C. Paul, and T. J. Weiler. Cosmic neutrino pevatrons: A brand new pathway to astronomy, astrophysics, and particle physics. Journal of High Energy Astrophysics, 1:1–30, May 2014.
 Abdo et al. [2015] A. A. Abdo, A. U. Abeysekara, R. Alfaro, B. T. Allen, C. Alvarez, J. D. Álvarez, R. Arceo, J. C. ArteagaVelázquez, T. Aune, H. A. Ayala Solares, A. S. Barber, B. M. Baughman, N. BautistaElivar, J. Becerra Gonzalez, E. Belmont, S. Y. BenZvi, D. Berley, M. Bonilla Rosales, J. Braun, R. A. CaballeroLopez, K. S. CaballeroMora, A. Carramiñana, M. Castillo, G. E. Christopher, U. Cotti, J. Cotzomi, E. de la Fuente, C. De León, T. DeYoung, R. Diaz Hernandez, L. DiazCruz, J. C. DíazVélez, B. L. Dingus, M. A. DuVernois, R. W. Ellsworth, D. W. Fiorino, N. Fraija, A. Galindo, F. Garfias, M. M. González, J. A. Goodman, V. Grabski, M. Gussert, Z. HampelArias, J. P. Harding, E. Hays, C. M. Hoffman, C. M. Hui, P. Hüntemeyer, A. Imran, A. Iriarte, P. Karn, D. Kieda, B. E. Kolterman, G. J. Kunde, A. Lara, R. J. Lauer, W. H. Lee, D. Lennarz, H. León Vargas, E. C. Linares, J. T. Linnemann, M. Longo, R. LunaGarcIa, J. H. MacGibbon, A. Marinelli, S. S. Marinelli, H. Martinez, O. Martinez, J. MartínezCastro, J. A. J. Matthews, J. McEnery, E. Mendoza Torres, A. I. Mincer, P. MirandaRomagnoli, E. Moreno, T. Morgan, M. Mostafá, L. Nellen, P. Nemethy, M. Newbold, R. NoriegaPapaqui, T. OcegueraBecerra, B. Patricelli, R. Pelayo, E. G. PérezPérez, J. Pretz, C. Rivière, D. RosaGonzález, E. RuizVelasco, J. Ryan, H. Salazar, F. Salesa, A. Sandoval, P. M. Saz Parkinson, M. Schneider, S. Silich, G. Sinnis, A. J. Smith, D. Stump, K. Sparks Woodle, R. W. Springer, I. Taboada, P. A. Toale, K. Tollefson, I. Torres, T. N. Ukwatta, V. Vasileiou, L. Villaseñor, T. Weisgarber, S. Westerhoff, D. A. Williams, I. G. Wisher, J. Wood, G. B. Yodh, P. W. Younk, D. Zaborov, A. Zepeda, and H. Zhou. Milagro limits and HAWC sensitivity for the ratedensity of evaporating Primordial Black Holes. Astroparticle Physics, 64:4–12, April 2015.
 Abeysekara et al. [2015] A. U. Abeysekara, R. Alfaro, C. Alvarez, J. D. Álvarez, R. Arceo, J. C. ArteagaVelázquez, H. A. Ayala Solares, A. S. Barber, B. M. Baughman, N. BautistaElivar, S. Y. BenZvi, M. Bonilla Rosales, J. Braun, K. S. CaballeroMora, A. Carramiñana, M. Castillo, U. Cotti, J. Cotzomi, E. de la Fuente, C. De León, T. DeYoung, R. Diaz Hernandez, B. L. Dingus, M. A. DuVernois, R. W. Ellsworth, D. W. Fiorino, N. Fraija, A. Galindo, F. Garfias, M. M. González, J. A. Goodman, M. Gussert, Z. HampelArias, J. P. Harding, P. Hüntemeyer, C. M. Hui, A. Imran, A. Iriarte, P. Karn, D. Kieda, G. J. Kunde, A. Lara, R. J. Lauer, W. H. Lee, D. Lennarz, H. León Vargas, J. T. Linnemann, M. Longo, R. LunaGarcía, K. Malone, A. Marinelli, S. S. Marinelli, H. Martinez, O. Martinez, J. MartínezCastro, J. A. J. Matthews, E. Mendoza Torres, P. MirandaRomagnoli, E. Moreno, M. Mostafá, L. Nellen, M. Newbold, R. NoriegaPapaqui, T. O. OcegueraBecerra, B. Patricelli, R. Pelayo, E. G. PérezPérez, J. Pretz, C. Rivière, D. RosaGonzález, H. Salazar, F. Salesa Greus, A. Sandoval, M. Schneider, G. Sinnis, A. J. Smith, K. Sparks Woodle, R. W. Springer, I. Taboada, K. Tollefson, I. Torres, T. N. Ukwatta, L. Villaseñor, T. Weisgarber, S. Westerhoff, I. G. Wisher, J. Wood, G. B. Yodh, P. W. Younk, D. Zaborov, A. Zepeda, H. Zhou, and HAWC Collaboration. Search for GammaRays from the Unusually Bright GRB 130427A with the HAWC GammaRay Observatory. Astrophysical Journal, 800:78, February 2015.
 Rädel and Wiebusch [2013] L. Rädel and C. Wiebusch. Calculation of the Cherenkov light yield from electromagnetic cascades in ice with Geant4. Astroparticle Physics, 44:102–113, April 2013.
 Patrignani et al. [2016] C. Patrignani et al. Review of Particle Physics. Chin. Phys., C40(10):100001, 2016.
 Jokisch et al. [1979] H. Jokisch, K. Carstensen, W. D. Dau, H. J. Meyer, and O. C. Allkofer. Cosmicray muon spectrum up to 1 tev at 75 zenith angle. Phys. Rev. D, 19:1368–1372, Mar 1979.
 Sarcevic [2007] Ina Sarcevic. Ultrahigh energy cosmic neutrinos and the physics beyond the standard model. volume 60, page 175, 2007.