\documentclass{nuws}
\usepackage{amsmath,amssymb,url}
%-------------------------------------------------------------------%
 \newif\ifpdf % Check for the existence of the variable \pdfoutput:
 \ifx\pdfoutput\undefined
 \pdffalse    % we are not running PDFLaTeX
 \else
 \pdfoutput=1 % we are running PDFLaTeX
 \pdftrue
 \fi \ifpdf   % Then use your new variable \ifpdf
     \usepackage[pdftex]{graphicx} % when using pdfLatex (preferred)
     \pdfcompresslevel=9
 \else
     \usepackage{graphicx}         % when using Latex and dvips (the
                                   % latter best with option -Pcmz,
                                   % if available, to invoke Type
                                   % 1 cm fonts)
 \fi
%-------------------------------------------------------------------%

\newcommand{\Red}{\special{color cmyk 0.50 0.85 1.00 0.0}}
\newcommand{\Black}{\special{color cmyk 0 0 0 1.}}
\newcommand{\Blue}{\special{color cmyk 1.0 0.9 0.1 0.0}}
\newcommand{\Green}{\special{color cmyk 0.92 0 0.59 0.7}}
\def\w{4.38cm} \def\gw{8.82cm} \def\W{8.51cm} \def\gW{11.40cm}
\urlstyle{rm}
\urldef{\FLUKA-3D}\url{http://www.mi.infn.it/~battist/neutrino.html}
\urldef{\Workshop}\url{http://www-zeuthen.desy.de/nuastro/
                   publications/sources/workshop2001/icpost_tr.html}


\begin{document}

\title{Atmospheric muons and neutrinos}
\author[1]{Vadim A. Naumov}%\thanksref{Supp}
\affil[1]{Dipartimento di Fisica and Sezione INFN di Ferrara,
          Via del Paradiso 12, I-44100 Ferrara, Italy and \\
          Laboratory for Theoretical Physics, Irkutsk State
          University, Gagarin boulevard 20, RU-664003 Irkutsk,
          Russia}
%\thanks[Supp]{Partially supported by grant from Ministry of
%              Education of Russian Federation, within the
%              framework of the program ``Universities of Russia
%              -- Basic Researches'', grant No.015.02.01.004.}
\correspondence{V.\,A.\,Naumov (naumov@fe.infn.it)}

\firstpage{1}
\pubyear{2001}

%\titleheight{11cm} % uncomment and adjust in case your title block
                    % does not fit into the default and minimum 7.5cm

\maketitle

\begin{abstract}
This paper is a mini-review of the atmospheric muon and neutrino
flux calculations based upon a recent version of {CORT} code and
up-to-date data on primary cosmic rays and hadronic interactions.
A comparison of calculations with a representative set of
atmospheric muon data for momenta below $\sim1$ TeV/c is presented.
The overall agreement between the calculated muon fluxes and the
data provides an evidence in favor of the validity of adopted
description of hadronic interactions and shower development.
In particular, this supports the low-energy atmospheric neutrino
fluxes predicted with CORT which are essentially lower than those
used in current analyses of the sub-GeV and multi-GeV neutrino
induced events in underground neutrino detectors.
\end{abstract}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\protect\section{Introduction}
\label{sec:Introduction}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

In recent paper by \citet{Fiorentini01a} an updated version of
FORTRAN code {CORT} ({\bf C}osmic-{\bf O}rigin {\bf R}adiati\-on
{\bf T}ransport) has been described and some results of low and
intermediate energy muon and neutrino flux calculations performed
with {CORT} have been discussed. A detailed comparison with the muon
fluxes and charge ratios measured in several modern balloon-borne
experiments suggests that the atmospheric neutrino flux is
essentially lower than one used in the current analyses of the
sub-GeV and multi-GeV neutrino induced events in underground
detectors. Some additional results which confirm this conclusion were
discussed in \citep{Fiorentini01b}.

Present work can be considered as an {\em addendum} to the
above-mentioned papers. Its main goal is to provide a comparison of
the predicted muon fluxes with a more representative set of data from
balloon-borne and ground-based experiments, including rather old (but
by no means outdated), which cover extensive ranges of muon momenta,%
\footnote{In this short write-up we limit ourselves with a selected
          set of muon data at momenta $p\lesssim1$ TeV/c (some of
          these data were taken from the recent compilations by
          \citet{Vulpescu98,Vulpescu01} and \citet{Hebbeker01}) and
          do not consider the underground and underwater muon data.
          Additional data can be found in the comprehensive
          reference manual by \citet{Grieder01}, and also in Refs.
          \citep{Bugaev98,Naumov98a,Naumov00a,%
          Sinegovskaya01,Naumov01b}.}
zenith angles and atmospheric depths. Such a comparison may be useful
to obtain further insight into the atmospheric neutrino problem.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% The array of cosmic-ray muon data accumulated by many balloon-borne
% and ground-based experiments are virtually unlimited. In this paper
% we consider only a small part of the world data.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\protect\section{The CORT code}
\label{sec:CORT}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Like the earlier versions of {CORT}
\citep{Naumov84,Bugaev84,Bugaev85,Bugaev86,Bugaev87a,Bugaev87b,%
Bugaev88,Bugaev89a,Bugaev90,Naumov93,Bugaev98},
the new Fortran 90 code implements a numerical integration of a
system of one-dimensional kinetic equations describing the production
and transport of nuclei, nucleons, light mesons, muons, neutrinos,
and antineutrinos of low and intermediate energies in the atmosphere.
It takes into account solar modulation, geomagnetic and
meteorological effects, energy loss of charged particles, muon
polarization and depolarization effects. The exact relativistic
kinematics is applied for description of particle interactions and
decays.  The new code has a number of options and run modes. It is
rather flexible and permits fast modification of many input data with
intrinsic switch keys.

In order to evaluate geomagnetic effects and to take into account the
anisotropy of the primary cosmic-ray flux in the vicinity of the
Earth, CORT uses the method of Ref. \citep{Naumov84} and detailed
maps of the effective vertical cutoff rigidities by \citet{Dorman71}.
The maps are corrected for the geomagnetic pole drift and compared
with the later results reviewed by \citet{Smart94} and with the
recent AMS % (Alpha Magnetic Spectrometer)
data the on the proton flux in near Earth orbit \citep{Alcaraz00a}.
The interpolation between the reference points of the maps is
performed by means of two-dimensional local B-spline. The Quenby-Wenk
relation \citep{Dorman71}, re-normalized to the vertical cutoffs, is
applied for evaluating the effective cutoffs for oblique directions.
More sophisticated effects, like the short-period variations of the
geomagnetic field, Forbush decrease, re-entrant cosmic-ray albedo
contribution, etc., are neglected. We also neglect the geomagnetic
bending of the trajectories of charged secondaries and multiple
scattering effects. Validity of our treatment of propagation of
secondary nucleons and nuclei has been verified using all available
data on secondary proton and neutron spectra in the atmosphere
\citep{Naumov84,Bugaev84,Bugaev85,Naumov00b,Naumov01a}.

The meteorological effects are included using the Dorman model of the
atmosphere \citep{Dorman72} which assumes an isothermal stratosphere
and constant gradient of temperature (as a function of depths) below
the tropopause. Ionization, radiative and photonuclear muon energy
losses are treated as continuous processes. This approximation is
quite tolerable for atmospheric depths $h\lesssim2\times10^3$
g/cm${}^2$ at all energies of interest \citep{Naumov94}. Propagation
of $\mu^+$ and $\mu^-$ originating from every source (pion or kaon
decay) is described by separate kinetic equations for muons with
definite polarization at production. These equations automatically
account for muon depolarization through the energy loss (but not
through the Coulomb scattering).


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\protect\section{Primary cosmic ray spectrum and composition}
\label{sec:Primaries}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

In the present calculations, the nuclear component of primary cosmic
rays is broken up into 5 principal groups: H, He, CNO, Ne-S and Fe
with average atomic masses $A$ of 1, 4, 15, 27 and 56, respectively.
We do not take into account the isotopic composition of the primary
nuclei and assume $Z=A/2$ for $A>1$, since the expected effect on the
secondary lepton fluxes is estimated to be small with respect to
present-day experimental uncertainties in the absolute cosmic-ray
flux and chemical composition.

We parametrize the spectra of the H and He groups at $E<120$
GeV/nucleon by fitting the data of the balloon-borne experiment BESS
obtained by a flight in 1998 \citep{Sanuki00,Sanuki01a}.
For higher energies (but below the knee) we use the data by a series
of twelve balloon flights of JACEE \citep{Asakimori98} and the result
of an analysis by \citet{Wiebel-Sooth98} based upon a representative
compilation of world data on primaries. Since the fits obtained by
\citet{Wiebel-Sooth98} turned out to be very close to the combined
result of the JACEE experiments, the model is called ``BESS+JACEE''
fit.

Figure \ref{f:PS} shows a comparison between the BESS+JACEE fit,
the data from \citet{Sanuki00,Sanuki01a} (BESS) and
\citet{Asakimori98} (JACEE), the fit from \citet{Wiebel-Sooth98}
(shaded areas) and the data from several other experiments
\citep{Mason72,Ryan72,Seo91,Beatty93,Ichimura93,Buckley94,Diehl97,%
      Menn00,Aharonian99,Bellotti99,Boezio99a,Alcaraz00b,Alcaraz00c},
performed in different periods of solar activity. The filled
areas in Fig. \ref{f:PS} represent power-type parametrizations
of the spectra derived by
\citet{Asakimori98,Ryan72,Ichimura93,Aharonian99,Alcaraz00b,%
      Alcaraz00c}
from the original data. The legends indicate the publication dates
and (when relevant) the years of measurements.
%-------------------------------------------------------------------%
\begin{figure}[h]
\centering
\ifpdf
\includegraphics[width=\W]{ps_pr.pdf} % for PDFLaTeX
\includegraphics[width=\W]{ps_he.pdf} % for PDFLaTeX
\else
\includegraphics[width=\W]{ps_pr.eps} % for LaTeX
\includegraphics[width=\W]{ps_he.eps} % for LaTeX
\fi
\protect\caption{Differential spectra of protons (top panel) and
                 helium nuclei (bottom panel) as a function of
                 kinetic energy per nucleon.
                 The data points and filled/shaded areas are for
                 the experimental data and power-law fits of the
                 data as is shown in the legends.
                 The solid curves represent the BESS+JACEE best fit
                 from \citet{Fiorentini01a} (see text for details).
\label{f:PS}}
\end{figure}
%-------------------------------------------------------------------%

We assume that the spectra of the remaining three nuclear groups are
similar to the helium spectrum. This assumption does not contradict
the world data for the CNO and Ne-S nuclear groups but works a bit
worse for the iron group.
Nevertheless, a more sophisticated model would be unpractical since
the corresponding correction would affect the secondary lepton fluxes
by a negligible margin.

In this paper we do not consider the effects of solar modulation.
Therefore the predicted muon and neutrino fluxes are to some extent
the maximum ones possible within our approach.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\protect\section{Nucleon-nucleus and nucleus-nucleus interactions}
\label{sec:interactions}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

All calculations with the earlier version of CORT
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%\citep{Naumov84,Bugaev84,Bugaev85,Bugaev86,Bugaev87a,Bugaev87b,%
%%Bugaev88,Bugaev89a,Bugaev90,Naumov93,Bugaev98},
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
were based on semiempirical models for inclusive nucleon and meson
production in collisions of nucleons with nuclei by
\citet{Kimel'74,Kimel'75,Serov73}.%
\footnote{See also Refs. \citet{Kalinovsky85,Sychev99} for the most
          recent versions.}
The Kimel'-Mokhov (KM) model is valid for projectile nucleon momenta
above $\sim4$ GeV/c and for the secondary nucleon, pion and kaon
momenta above 450, 150 and 300 MeV/c, respectively. Outside these
ranges (that is mainly within the region of resonance production of
pions) the Serov-Sychev (SS) model was used.

Both models are in essence comprehensive parametrizations of the
relevant accelerator data. It is believed that the combined
``KM+SS'' model provides a safe and model-inde\-pendent basis to the
low-energy atmospheric muon and neutrino calculations.
However it is not free of uncertainties.
For the present calculation, the fitting
parameters of the KM model for meson and nucleon production off
different nuclear targets were updated using accelerator data not
available for the original analysis
\citep{Kimel'74,Kimel'75,Kalinovsky85}. The values of the parameters
were extrapolated to the air nuclei (N, O, Ar, C).  The overall
correction is less than 10-15\% within the kinematic regions
significant to atmospheric cascade calculations. Besides that
energy-dependent correction factors were introduced into the model to
tune up the output $\pi^+/\pi^-$ ratio taking into account the
relevant new data.

The processes of meson regeneration and charge exchange
($\pi^\pm+\mathrm{Air}\to\pi^{\pm(\mp)}+X$ etc.) are not of critical
importance for production of leptons with energies of our interest
and can be considered in a simplified way. Here we use a proper
renormalization of the meson interaction lengths, which was deduced
from the results by \citet{Vall86} obtained for high-energy cascades.

The next important ingredient of any cascade calculations is a model
for nucleus-nucleus collisions. Here we consider a modest
generalization of a simple ``Glauber-like'' model used in
\citep{Naumov84,Bugaev85}. Namely, we write the inclusive spectrum of
secondary particles $c$ ($c=p$, $n$, $\pi^\pm$, $K^\pm$,
$K^0,\ldots$) produced in collisions of nuclei as
\begin{align*}
\frac{dN_{\mathrm{AB} \to cX}}{d x} = & \xi_{\mathrm{AB}}^c
\left[Z \frac{dN_{p\mathrm{B} \to cX}}{dx}+
   (A-Z)\frac{d N_{n\mathrm{B} \to cX}}{dx}\right]\nonumber\\
&+\left(1-\xi_{\mathrm{AB}}^c\right)
\left[Z\delta_{cp}+(A-Z)\delta_{cn}\right]\delta(1-x).
\end{align*}
Here $dN_{N\mathrm{B} \to cX}/dx$ is the spectrum of particles $c$
produced in a collision of a nucleon $N$ ($N=p,n$) with nucleus
B, $x$ is the Feynman variable, and $\xi_{\mathrm{AB}}^c$ is the
average fraction of inelastically interacting nucleons of
the projectile nucleus A. The term proportional to delta function
describes the contribution of ``spectator'' nucleons from the
projectile nucleus.

In the standard Glauber-Gribov multiple scattering theory the
quantity $\xi_{\mathrm{AB}}^c$ is certainly independent of the type
of inclusive particle $c$. On the other hand, it depends of the type
of nucleus collision. Indeed, essentially all nucleons participate in
the central AB collisions ($\xi_{\mathrm{AB}}^c\simeq1$)%
\footnote{Here we assume for simplicity that the atomic weight of
          the projectile nucleus is not much larger than that of
          the target nucleus.}
while, according to the well-known Bialas--Bleszy\'nski--Czy\.z (BBC)
relation \citep{Bialas76},
$\xi_{\mathrm{AB}}^c=\sigma_{N\mathrm{B}}^{\mathrm{inel}}/
                     \sigma_{\mathrm{AB}}^{\mathrm{inel}}$
for the minimum bias collisions.

To use the above model in a cascade calculation one should take into
account that nucleons and mesons are effectively produced in nuclear
collisions of different kind. Namely, the contribution from central
collisions is almost inessential for the nucleon component of the
cascade but quite important for light meson production. Thus one can
expect that effectively
$\xi_{\mathrm{AB}}^{\pi,K}>\xi_{\mathrm{AB}}^{p,n}$. We use the BBC
relation for nucleon production by any nucleus while for meson
production we put $\xi_{\mathrm{He-Air}}^{\pi,K}=\xi$, where $\xi$ is
a free parameter. Variations of this parameter within the
experimental limits yield a comparatively small effect to the muon
fluxes (except for very high altitudes) and inessential
($\lesssim6$\%) effect to the neutrino fluxes at sea level
\citep{Fiorentini01a}. Effect of similar variations of the parameters
$\xi_{\mathrm{A-Air}}^{\pi,K}$ for other nuclear groups is completely
negligible.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\protect\section{Numerical results for muons}
\label{sec:NumRes}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Figure \ref{f:CAPRICE94-98} displays a comparison of the
calculated momentum spectra of $\mu^+$ and $\mu^-$ for 10 atmospheric
depth ranges $\Delta h_i= 15-35, 35-65, 65-90, 90-120, 120-150,
150-190, 190-250, 250-380, 380-580$, and $580-890$ g/cm${}^2$ with
the data of two balloon-borne experiments CAPRICE\,94
\citep{Boezio99b,Boezio00} and CAPRICE\,98 \citep{Hansen01}
(the latter experimental data are {\em preliminary}).
The nominal geomagnetic cutoff rigidity $R_c$ is about 0.5 GV
(4.5 GV) and the detection cone is about $20^\circ$ ($14^\circ$)
around the vertical direction with the average incident angle of
about $9^\circ$ ($8^\circ$) for CAPRICE\,94 (CAPRICE\,98).
The data points in Fig. \ref{f:CAPRICE94-98} are the muon
intensities at the Flux-weighted Average Depths (FAD) defined in
\citep{Boezio99b}, while the filled areas display the calculated
variations of the muon spectra inside the ranges $\Delta h_i$
($i=1,\ldots,10$). Namely, they are obtained by considering the
minimal and maximal muon fluxes within each range $\Delta h_i$.
The calculations are performed for the conditions of the
experiment. The parameter $\xi$ is taken to be 0.685
(the best value, according to \citet{Fiorentini01a}).%
\footnote{It has been shown \citep{Fiorentini01a} that
          the indetermination of $\xi$ is only
          significant for $h<(15-20)$ g/cm${}^2$.}

%-------------------------------------------------------------------%
\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
\begin{figure*}[t]
\centering
\ifpdf
\includegraphics[width=\gw]{c94f.pdf}~% for PDFLaTeX
\includegraphics[width=\gw]{c98f.pdf} % for PDFLaTeX
\else
\includegraphics[width=\gw]{c94f.eps}~% for LaTeX
\includegraphics[width=\gw]{c98f.eps} % for LaTeX
\fi
\protect\caption{Differential momentum spectra of $\mu^+$ and $\mu^-$
                 for 10 atmospheric depth ranges $\Delta h_i$
                 (are indicated on the left of each panel).
                 The data points are from the CAPRICE\,94 experiment
                 \citep{Boezio99b,Boezio00} and from the
                 {\em preliminary} analysis of the CAPRICE\,98
                 experiment \citep{Hansen01}.
                 The filled areas display the expected variations
                 of the muon fluxes within the ranges $\Delta h_i$.
                 All the data are scaled with the factors
                 indicated at the first and third panels on the
                 right.
\label{f:CAPRICE94-98}}
\end{figure*}
%-------------------------------------------------------------------%
\begin{figure*}[b]
\centering
\ifpdf
\includegraphics[width=\linewidth]{cgl.pdf} % for PDFLaTeX
\else
\includegraphics[width=\linewidth]{cgl.eps} % for LaTeX
\fi
\protect\caption{Differential momentum spectra of positive and
                 negative muons at the ground level.
                 The data points are from the experiments
                 CAPRICE\,94 \citep{Boezio99b,Boezio00,Kremer99} and
                 CAPRICE\,97 \citep{Kremer99}. The $\mu^+$ data are
                 scaled with a factor of 3. The solid (narrow) bands
                 and dotted curves are calculated with CORT as is
                 explained in legends.
\label{f:CAPRICE94-97}}
\end{figure*}
\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
%-------------------------------------------------------------------%

It is important to note that the thickness of the bands is relatively
small just for the region of effective muon and neutrino production
that is in the neighborhood of the broad maxima of the muon flux
($100-300$~g/cm${}^2$). By this is meant that, in this region,
an error in evaluation of the FAD cannot introduce essential
uncertainty. Outside the region of effective production of leptons,
the amplitude of the muon flux variations increases with decreasing
muon momenta on account for the strong dependence of the meson
production rate upon the depth and the growing role of the muon
energy loss and decay at $h\gtrsim300$~g/cm${}^2$.

Figure \ref{f:CAPRICE94-97} shows a comparison of the calculated
differential momentum spectra of $\mu^+$ and $\mu^-$ with the
most accurate current data obtained at ground level in the
experiments CAPRICE\,94 \citep{Boezio99b,Boezio00,Kremer99}
($h=1000$ g/cm$^2$, $R_c=0.5$ GV) and CAPRICE\,97 \citep{Kremer99}
($h=886$ g/cm$^2$, $R_c=4.2$ GV).
The detection cone in both experiments was the same as described
above for CAPRICE\,94. The calculation are done by using the KM+SS
interaction model. Variation of the parameter $\xi$ from 0.517
(the BBC value) to 0.710 (an experimental upper limit derived from
the data on interactions of $\alpha$ particles with light nuclei).%
\footnote{See, e.g., Ref. \citep{Kowalski81} and references
         therein.}
leads to the almost negligible ($\lesssim3$\%) effect in the
ground-level muon flux.
For reference, the spectra calculated by using the superposition
model for nucleus-nucleus interactions are also added. The data
below $\sim10$ GeV/c seem to be precise enough in order to
conclude that they disfavor the superposition model.

As is seen from Figs. \ref{f:CAPRICE94-98} and \ref{f:CAPRICE94-97},
there is a substantial agreement between the calculations with CORT
and the data of the CAPRICE experiments within a wide range
atmospheric depths. This agreement is {\em none the worse} than it
is for the recent 3D Monte Carlo calculations \citep{Battistoni01a,%
Engel01,Honda01b,Liu01,Poirier01,Wentz01a,Wentz01b}.%
\footnote{See also Ref. \citep{Fiorentini01b} for the direct
          comparison with the FLUKA 3D result.}

Figure \ref{f:Geo} shows expected geomagnetic effect for vertical
differential momentum spectra of $\mu^+$ and $\mu^-$ and for muon
charge ratio $\mu^+/\mu^-$ at sea level. A comparison of these
predictions with experiment will be discussed elsewhere.

Figure \ref{f:FslSurvey} collects the data on the near-vertical
differential momentum spectrum of $\mu^++\mu^-$ at ground level
from many experiments performed over the past five decades
% over the course of the past semicentennial period
\citep{Caro50,Owen55,Pine59,Pak61,Holmes61,Hayman62a,Hayman62b,%
       Aurela67,Baber68a,Allkofer70a,Allkofer70b,Allkofer71,%
       Bateman71,Nandi72a,Allkofer75,Ayre75,Thompson77,Baschiera79,%
       Green79,Rastin84a,Tsuji98,Kremer99,LeCoultre01}.
Only the data for muon momenta below $\sim1$ TeV/c are shown.

The compilation, of course, is not exhaustive. In particular, it
does not involve the low-energy data for large geomagnetic cutoffs,
the ground-level data obtained at the altitudes far different from
sea level, and also indirect muon data from underground
measurements.%
\footnote{These were discussed in detail by \citet{Bugaev98}.}
The spectra calculated for $h=10^3$ g/cm$^2$, $\theta=0^\circ$,
and $R_c=0$ with the new and old (from \citet{Bugaev98}) versions
of CORT are plotted. The best-fit spectrum to Nottingham data
\citep{Rastin84a} and the spectrum calculated by \citet{Agrawal96}
are also shown in bottom panel.
%-------------------------------------------------------------------%
\begin{figure}[htb]
\centering
\ifpdf
\includegraphics[width=8.4cm]{geof.pdf} % for PDFLaTeX
\includegraphics[width=8.4cm]{geor.pdf} % for PDFLaTeX
\else
\includegraphics[width=8.4cm]{geof.eps} % for LaTeX
\includegraphics[width=8.4cm]{geor.eps} % for LaTeX
\fi
\protect\caption{Vertical differential momentum spectra of positive
                 and negative muons and muon charge ratios at
                 sea level calculated with CORT for several
                 geomagnetic cutoffs $R_c$. The $\mu^+$ spectra on
                 the top panel are scaled with a factor of 2.
\label{f:Geo}}
\end{figure}
%-------------------------------------------------------------------%

%-------------------------------------------------------------------%
\begin{figure*}
\centering
\ifpdf
\includegraphics[width=0.90\linewidth]{fsltot.pdf} % for PDFLaTeX
\includegraphics[width=0.90\linewidth]{fsl-he.pdf} % for PDFLaTeX
\else
\includegraphics[width=0.90\linewidth]{fsltot.eps} % for LaTeX
\includegraphics[width=0.90\linewidth]{fsl-he.eps} % for LaTeX
\fi
\protect\caption{Near-vertical differential momentum spectrum of
                 $\mu^++\mu^-$ at ground level (top panel) and its
                 rescaled fragment for $p>10$ GeV/c (bottom panel).
                 The thickness of the band calculated with CORT
                 reflects the  indetermination in the parameter
                 $\xi$. %The data are described in the text.
\label{f:FslSurvey}}
\end{figure*}
%-------------------------------------------------------------------%

In fact, not all experiments listed above measure at sea level.
Moreover, the experiments were performed in different locations
and periods. Therefore the solar modulation effects, meteorological
and geomagnetic conditions are generally different in these
experiments.
%-------------------------------------------------------------------%
\begin{figure*}
\centering
\ifpdf
\includegraphics[width=\linewidth]{crsl.pdf} % for PDFLaTeX
\else
\includegraphics[width=\linewidth]{crsl.eps} % for LaTeX
\fi
\protect\caption{Muon charge ratio at sea level for near-vertical
                 direction.
%                The data of many experiments were obtained at
%                different geomagnetic locations.
                 Geomagnetic cutoffs $R_c$ for different experiments
                 are shown in the legend.
\label{f:CRSL}}
\end{figure*}
%-------------------------------------------------------------------%
The same is true for the detection cones, shielding factors, adopted
procedures for the momentum spectrum unfolding, and so on. All these
factors might cause a bias when all the data are compared together
\citep{Hebbeker01}. Here we do not pursue any correction procedure
(like one suggested by \citet{Hebbeker01}) and plot all the
experimental data as they are. Instead we show the world average fit
derived by \citet{Hebbeker01} from an analysis of (almost) the same
data set at $p\geq10$~GeV/c.
It can be seen that the ground-level experiments are in rather poor
agreement to one another. A large share of the data are arranged
outside (or even far from) the borderline of the world average fit.
However the data of the most recent experiments, specifically
CAPRICE\,94 \citep{Kremer99} and L3+C \citep{LeCoultre01} are
placed close to or within these bounds.%
\footnote{Note that the L3+C data on muon flux and charge ratio
          reported in ICRC'27 \citep{LeCoultre01} (see also Refs.
          \citep{Petersen01,Ladron01}), are preliminary and may be
          subject to refinement (within the systematic error of
          7.7\%) when analyses are complete.}

Fig. \ref{f:FslSurvey} suggests that the muon flux calculated
with the new CORT is systematically lower than that obtained by
\citet{Bugaev98} using the earlier version. The deviation is about
18\% for $p=10$~GeV/c and about a few per cent for
$p\gtrsim80$~GeV/c.
The main reason of this difference is that the BESS+JACEE primary
spectrum below $\sim200$~GeV/nucleon is essentially lower than
that used by \citet{Bugaev98} and in the earlier calculations
with CORT. Another reasons pertain to the changes in the inclusive
cross sections and improved description of muon propagation.
Below $\sim400$ GeV/c the new spectrum is close to the world
average fit and to the data from CAPRICE\,94 and L3+C experiments
as well as to the data from several experiments carried out by
the BESS Collaboration \citep{Sanuki01a,Motoki01,Sanuki01b}.%
\footnote{Unfortunately, the tabulated muon data are still
          unavailable from the BESS Collaboration.}

Let us now dwell to the muon charge ratio. It has been shown
\citep{Fiorentini01a} that calculations with CORT describe well
the modern data on the charge ratio obtained at different atmospheric
depths. Here we consider only the ground-level data for near-vertical
directions. The world survey is presented in Fig. \ref{f:CRSL}.
The data are taken from many experiments
\citep{Owen49,Conversi50,Owen51,Moroney54,Filosofo54,Fukui55,%
       Holmes61,Pak61,Coates62,Hayman62a,Hayman62b,Kawaguchi65,%
       Rastin65,Baber68b,Allkofer68,Allkofer70a,Appleton71,%
       Allkofer72,Nandi72b,Burnett73,Abdel-Monem73,Ayre73,%
       Baxendale75,Singhal83,Rastin84b,Stephens87,Grandegger93,%
       Jannakos95,Krizmanic95,Boezio00,Coutu00,Brancus00,%
       Tsuji01,LeCoultre01},
performed at different geomagnetic locations. According to CORT (see
bottom panel of Fig. \ref{f:Geo}), the sea-level charge ratio is
slightly affected by the geomagnetic field. However the dispersion
of the experimental points is too large and they do not show some
significant correlation with the geomagnetic cutoff.
%-------------------------------------------------------------------%
\begin{figure}[hbt]
\centering
\ifpdf
\includegraphics[width=8.3cm]{amh.pdf} % for PDFLaTeX
\else
\includegraphics[width=8.3cm]{amh.eps} % for LaTeX
\fi
\protect\caption{Differential momentum spectra of muons at sea level
                 for several zenith angles and angular bins. The data
                 points are from AMH spectrometer \citep{Green79}.
                 The filled areas display the expected variations of
                 the muon fluxes inside the bins
                 \protect$\Delta\theta$ indicated at the left.
                 The solid curves are calculated for average zenith
                 angles \protect$\langle\theta\rangle$. All the data
                 are scaled with the factors shown in the parentheses
                 at the right.
\label{f:AMH}}
\end{figure}
%-------------------------------------------------------------------%
%-------------------------------------------------------------------%
\begin{figure*}[bt]
\centering
\ifpdf
\includegraphics[width=\gw]{f30.pdf}~% for PDFLaTeX
\includegraphics[width=\gw]{f75.pdf} % for PDFLaTeX
\else
\includegraphics[width=\gw]{f30.eps}~% for LaTeX
\includegraphics[width=\gw]{f75.eps} % for LaTeX
\fi
\protect\caption{Differential momentum spectra of muons at sea level
                 for $\theta=30^\circ$ and $75^\circ$. The data
                 points are from \citep{Kellogg78,Jokisch79,Tsuji98}.
                 The filled areas display expected variations of the
                 muon fluxes inside the indicated angular bins.
                 The solid curves are for the corresponding
                 average zenith angles.
\label{f:F30-75}}
\end{figure*}
%-------------------------------------------------------------------%

Figure \ref{f:AMH} shows a comparison with the sea-level data from
the AMH magnetic spectrometer (Texas A \& M and University of
Houston Collaboration) measured within different zenith-angle bins
$\Delta\theta$ \citep{Green79}. The spectra for average zenith
angles and expected variations of the spectra inside the
angular bins are shown by solid lines and filed areas, respectively.
One can see a good or at least qualitative agreement with the data
for $p\lesssim100$ GeV/c and $\theta\lesssim80^\circ$. At higher
momenta and at large zenith angles the situation is spoiled.
However, a comparison between the AMH and world survey data
for vertical (see Fig.~\ref{f:FslSurvey}) suggests that there
is some systematic bias in the AMH experiment above 100 GeV/c.
The abnormal scatter of points in the near-horizontal bin
is indicative of a flaw in the large-angle measurements.

Figure \ref{f:F30-75} displays a comparison between the differential
momentum spectra of muons at ground level calculated with CORT and
the data from Brookhaven magnetic spectrometer \citep{Kellogg78},
Kiel-DESY muon spectrometer located at Hamburg \citep{Jokisch79},
and Okayama altazimuthal counter cosmic-ray telescope with a magnet
spectrometer \citep{Tsuji98}.
The filled areas display the expected variations of the muon
fluxes inside the angular bins $\Delta\theta=(26-34)^\circ$
(left panel) and $(69-81)^\circ$ (right panel).
The solid curves are for the corresponding average zenith angles
of $30^\circ$ and $75^\circ$. The calculations presented in Fig.
\ref{f:F30-75} are done without taking into account the geomagnetic
effects (both the primary spectrum cutoff and muon deflection).
These are not important under the conditions of Brookhaven and
Kiel-DESY experiments but requisite for Okayama site.%
\footnote{The Okayama telescope location is
          ($34^\circ40'$\,N,\,$133^\circ56'$\,E) and vertical
          geomagnetic cutoff rigidity is about 12.4 GV.}

Generally the geomagnetic effects must decrease the low-momentum part
of the expected muon flux (see top panel of Fig. \ref{f:Geo}). The
calculated spectra fit rather well the precise Brookhaven and
Kiel-DESY data but are in rather bad agreement (especially for
$\theta=30^\circ$) with the Okayama telescope data of poorer
statistics; the geomagnetic corrections can only aggravate
the disagreement. One can conclude from Fig. \ref{f:F30-75} that
there is a systematic flaw in the average zenith angles measured
in the Okayama experiment.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\protect\section{Numerical results for neutrinos}
\label{sec:Neutrinos}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Let us briefly sketch some results on atmospheric neutrinos (AN).
Due to geomagnetic effects, the low-energy AN spectra and angular
distributions are quite different for different sites of the globe.
Figures \ref{f:ANFluxLE} and \ref{f:ANRatioLE} display our
predictions for ten underground neutrino laboratories listed in
Table \ref{t:Labs}.
%-------------------------------------------------------------------%
\begin{table}[htb]
\centering
\protect\caption{List of some underground laboratories.
\label{t:Labs}}
\vspace{3mm}
\ifpdf
\includegraphics[width=\linewidth]{labs.pdf} % for PDFLaTeX
\else
\includegraphics[width=\linewidth]{labs.eps} % for LaTeX
\fi
\end{table}
%-------------------------------------------------------------------%

Figure \ref{f:ANFluxLE} shows the $\nu_e$, $\overline{\nu}_e$,
$\nu_\mu$ and $\overline{\nu}_\mu$ energy spectra averaged over all
zenith and azimuth angles. The ratios of the AN fluxes averaged over
the lower and upper semispheres (``up-to-down'' ratios) are shown in
Fig. \ref{f:ANRatioLE}. As a result of geomagnetic effects, both the
spectra and the up-to-down ratios for the following 6 groups of the
labs: SOUDAN + SNO + IMB, HPW, NUSEX + Fr\'ejus, Gran Sasso + Baksan,
Kamioka, and KGF are quite distinct for energies below a few GeV.

Figure \ref{f:K-GS} depicts the zenith-angle distributions of
$\nu_e$, $\overline{\nu}_e$, $\nu_\mu$ and $\overline{\nu}_\mu$
for Kamioka and Gran Sasso labs calculated with CORT using its
``standard''(KM+SS) model of particle interaction
(see sec. \ref{sec:interactions}) and also the TARGET model for
$\pi/K$ meson production (including the superposition model for
collisions of nuclei) used by Bartol group
\citep{Barr89,Agrawal96,Lipari98}.
The TARGET model for nucleon production is not included into the
CORT options yet. Instead the KM+SS model is used everywhere.
Below, we refer to CORT switched into the
TARGET meson production model as ``CORT+TARGET''.
For comparison, the result of the recent 3D calculation by
\citet{Battistoni00} based on the FLUKA 3D Monte Carlo simulation
package is also shown. It allows to ``highlight'' the 3D effects
which are very dependent on neutrino energy and direction of
arrival.
%-------------------------------------------------------------------%
\begin{figure}[thb]
\centering
\ifpdf
\includegraphics[width=\linewidth]{a.pdf} % for PDFLaTeX
\else
\includegraphics[width=\linewidth]{a.eps} % for LaTeX
\fi
\protect\caption{Scaled $4\pi$ averaged fluxes of $\nu_e$,
                 $\overline{\nu}_e$, $\nu_\mu$, and
                 $\overline{\nu}_\mu$ for ten underground
                 laboratories (see table \ref{t:Labs} for notation).
\label{f:ANFluxLE}}
\end{figure}
%-------------------------------------------------------------------%
%-------------------------------------------------------------------%
\begin{figure}[thb]
\centering
\ifpdf
\includegraphics[width=\linewidth]{r.pdf} % for PDFLaTeX
\else
\includegraphics[width=\linewidth]{r.eps} % for LaTeX
\fi
\protect\caption{Up-to-Down ratios of the $\nu_e$,
                 $\overline{\nu}_e$, $\nu_\mu$, and
                 $\overline{\nu}_\mu$ fluxes for ten underground
                 laboratories (see table \ref{t:Labs} for notation).
\label{f:ANRatioLE}}
\end{figure}
%-------------------------------------------------------------------%

A more quantitative comparison is given by Table \ref{t:FluxRatios}
which tabulates the ratios of the$\nu_e$, $\overline{\nu}_e$,
$\nu_\mu$, and $\overline{\nu}_\mu$ fluxes calculated with
CORT+TARGET and FLUKA to those with the ``standard'' CORT.
The fluxes are averaged over all directions and over several
energy bins. Table \ref{t:FluxRatios} also shows the
so-called flavor ratio%
\footnote{This quantity roughly represents the ratio of $e$
          like to $\mu$ like single-ring contained events measured
          in water Cherenkov detectors and the ``showers-to-tracks''
          ratio measured in iron detectors.}
\[
R_\nu=\left(\nu_e  +\tfrac{1}{3}\overline{\nu}_e  \right)/
      \left(\nu_\mu+\tfrac{1}{3}\overline{\nu}_\mu\right),
\]
evaluated with CORT, CORT+TARGET, and FLUKA for the same energy
bins. % and averaged in the same way.

By using Fig. \ref{f:K-GS} and Table \ref{t:FluxRatios}, one can
conclude that the present calculations with CORT for neutrino
energies below 1-2 GeV lead to the AN fluxes which are essentially
{\em lower} than those obtained by \citet{Barr89,Agrawal96,Lipari98}
and by \citet{Honda90,Honda95,Honda96} and those are in current use
for many analyses of sub-GeV and multi-GeV $\nu$ induced events in
underground detectors.%
\footnote{A comparison between the earlier AN flux calculations
          has been discussed in detail by \citet{Gaisser96}.}
On the other hand, the new AN fluxes are rather close to the
those obtained with the earlier versions of CORT
\citep{Bugaev87a,Bugaev87b,Bugaev89a,Bugaev90}. The main
differences (for neutrino energies above 1-2 GeV) are due to the
new data for the primary cosmic-ray spectrum and
composition (see Section \ref{sec:Primaries}).


%-------------------------------------------------------------------%
\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
\begin{figure*}[t]
\centering
\ifpdf
\includegraphics[width=\w]{e1k.pdf}~% for PDFLaTeX
\includegraphics[width=\w]{e2k.pdf}~% for PDFLaTeX
\includegraphics[width=\w]{e1g.pdf}~% for PDFLaTeX
\includegraphics[width=\w]{e2g.pdf} % for PDFLaTeX
\includegraphics[width=\w]{m1k.pdf}~% for PDFLaTeX
\includegraphics[width=\w]{m2k.pdf}~% for PDFLaTeX
\includegraphics[width=\w]{m1g.pdf}~% for PDFLaTeX
\includegraphics[width=\w]{m2g.pdf} % for PDFLaTeX
\else
\includegraphics[width=\w]{e1k.eps}~% for LaTeX
\includegraphics[width=\w]{e2k.eps}~% for LaTeX
\includegraphics[width=\w]{e1g.eps}~% for LaTeX
\includegraphics[width=\w]{e2g.eps} % for LaTeX
\includegraphics[width=\w]{m1k.eps}~% for LaTeX
\includegraphics[width=\w]{m2k.eps}~% for LaTeX
\includegraphics[width=\w]{m1g.eps}~% for LaTeX
\includegraphics[width=\w]{m2g.eps} % for LaTeX
\fi
\protect\caption{Zenith-angle distributions of $\nu_e$,
                 $\overline{\nu}_e$, $\nu_\mu$, and
                 $\overline{\nu}_\mu$ for Kamioka and Gran Sasso.
                 Solid curves represent the result of CORT obtained
                 with its ``standard'' interaction model.
                 The dashed curves are the result of CORT obtained
                 with the TARGET model for meson production and
                 superposition model for collisions of nuclei.
                 The circles are for the result of calculation
                 by \citet{Battistoni00} based on FLUKA 3D code.
                 The distributions are averaged over the azimuth
                 angle and over the energy bins indicated near the
                 curves.
\label{f:K-GS}}
\end{figure*}
%-------------------------------------------------------------------%
%-------------------------------------------------------------------%
\begin{table*}[b]
\centering
\protect\caption{The ratios of $4\pi$ averaged AN fluxes obtained
                 with CORT+TARGET and FLUKA to those with CORT and
                 neutrino flavor ratios $R_\nu$ calculated with CORT,
                 CORT+TARGET, and FLUKA for Kamioka.}
\label{t:FluxRatios}
\vspace{4mm}
\begin{tabular}{|c|cccc|cccc|}
\hline
               & \multicolumn{8}{c|}%
{{\bf\Red Average fluxes normalized to CORT\Black}}                 \\
\hline
$\Delta E_\nu$ & \multicolumn{4}{c|}{CORT+TARGET}
               & \multicolumn{4}{c|}{FLUKA 3D}                      \\
  (GeV)        &
 $\nu_e$    & $\overline{\nu}_e$ & $\nu_\mu$ & $\overline{\nu}_\mu$ &
 $\nu_e$    & $\overline{\nu}_e$ & $\nu_\mu$ & $\overline{\nu}_\mu$ \\
\hline      &      &      &      &      &      &      &      &      \\
  0.1--0.2  & 1.64 & 1.78 & 1.72 & 1.73 & 1.25 & 1.45 & 1.41 & 1.40 \\
  0.2--0.3  & 1.51 & 1.68 & 1.60 & 1.61 & 1.23 & 1.38 & 1.34 & 1.33 \\
  0.3--0.5  & 1.42 & 1.61 & 1.48 & 1.50 & 1.17 & 1.28 & 1.22 & 1.22 \\
  0.5--0.7  & 1.34 & 1.50 & 1.36 & 1.38 & 1.10 & 1.19 & 1.11 & 1.12 \\
  0.7--1.0  & 1.28 & 1.41 & 1.27 & 1.30 & 1.05 & 1.12 & 1.04 & 1.05 \\
  1.0--2.0  & 1.20 & 1.30 & 1.17 & 1.20 & 1.02 & 1.05 & 0.99 & 1.01 \\
  2.0--3.0  & 1.10 & 1.17 & 1.03 & 1.08 & 0.98 & 1.03 & 0.96 & 0.97 \\
  3.0--5.0  & 1.02 & 1.07 & 0.95 & 1.01 & 0.99 & 1.02 & 0.95 & 0.97 \\
  5.0--7.0  & 0.95 & 0.99 & 0.91 & 0.96 & 1.01 & 1.05 & 0.94 & 0.99 \\
  7.0--10.0 & 0.91 & 0.95 & 0.89 & 0.93 & 1.03 & 1.09 & 0.95 & 1.04 \\
 10.0--20.0 & 0.87 & 0.91 & 0.88 & 0.91 & 1.07 & 1.12 & 0.96 & 1.02 \\
 20.0--30.0 & 0.83 & 0.87 & 0.87 & 0.89 & 1.10 & 1.15 & 0.94 & 1.02 \\
            &      &      &      &      &      &      &      &      \\
\hline
\end{tabular}
\hspace{3mm}
\begin{tabular}{|c|c|c|c|}
\hline
& \multicolumn{3}{c|}{{\bf\Red Flavor ratios $R_\nu$\Black}} \\
\hline
$\Delta E_\nu$ & CORT & CORT+TARGET & FLUKA 3D \\
   (GeV)       &      &             &          \\
\hline         &      &             &          \\
  0.1--0.2     & 0.52 &    0.51     &  0.48    \\
  0.2--0.3     & 0.52 &    0.50     &  0.49    \\
  0.3--0.5     & 0.51 &    0.50     &  0.50    \\
  0.5--0.7     & 0.50 &    0.50     &  0.50    \\
  0.7--1.0     & 0.49 &    0.50     &  0.50    \\
  1.0--2.0     & 0.47 &    0.49     &  0.49    \\
  2.0--3.0     & 0.44 &    0.47     &  0.45    \\
  3.0--5.0     & 0.40 &    0.43     &  0.42    \\
  5.0--7.0     & 0.36 &    0.37     &  0.38    \\
  7.0--10.0    & 0.32 &    0.32     &  0.34    \\
 10.0--20.0    & 0.26 &    0.26     &  0.29    \\
 20.0--30.0    & 0.20 &    0.19     &  0.23    \\
               &      &             &          \\
\hline
\end{tabular}
\end{table*}

\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

The AN angular distributions calculated with CORT and with FLUKA
3D code \citep{Battistoni00,Battistoni01b} are in good agreement
for any zenith angle at energies above 0.7--0.8 GeV. On the
other hand, CORT is systematically lower than FLUKA for
$E_\nu\lesssim1$~GeV. At the energy $E_\nu=0.5$~GeV the FLUKA
$4\pi$ averaged $\nu_e$, $\overline{\nu}_e$, $\nu_\mu$ and
$\overline{\nu}_\mu$ fluxes exceed the corresponding CORT results
by about 15\%, 27\%, 18\% and 18\%, respectively. This discrepancy
is only partially due to 3D effects (Fig. \ref{f:K-GS}), which can
account for an increase $\lesssim10\%$, and it is most probably
related to differences between the hadronic interaction models
adopted in FLUKA and CORT codes.

It is pertinent to note that several recent 3D Monte Carlo
calculations of the AN flux
\citep{Honda01a,Honda01b,Liu01,Plyaskin01,Tserkovnyak01,Wentz01a}
confirmed the sharp enhancement of low-energy neutrino intensities
for near-horizontal directions predicted with FLUKA. However
the quantitative characteristics of this enhancement are very
model-dependent and vary over a broad range from one calculation
to another.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% We also note that the tendency has been for almost each group
% to decrease the AN fluxes.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\protect\subsection*{\em Prompt neutrinos}

Figures \ref{f:Nu_e} and \ref{f:Nu_m} collect the differential
energy spectra of downward going atmospheric neutrinos calculated
within a wide energy range (from 50 MeV to about 20 EeV) for 11
zenith angles. Figures show the ``conventional'' neutrino
contribution (originated from decay of pions, kaons and muons)
and the total AN spectra which include the ``prompt'' neutrino
contribution originated from semileptonic decays of charmed hadrons
(mainly $D^\pm$, $D^0$, $\overline{D}^0$ mesons and $\Lambda_c^+$
hyperons).%
\footnote{We do not discuss here the contributions from $D_s$ mesons
          and heavier states (such as $b\overline{b}$,
          $t\overline{t}$, $W$, $Z$) which become important sources
          of atmospheric (and extraterrestrial) prompt $\tau$ neutrinos
          at very high energies \citep{Pasquali99b,Athar01}.}
The prompt neutrino contribution must dominate at very high energies.
However the charm hadroproduction cross sections are very
model-dependent and cannot be unambiguously predicted for lack of a
generally accepted model. As a result the prompt neutrino
contribution and even the energies above which the prompt muon and
electron neutrinos become dominant are very uncertain as yet
\citep{Bugaev98,Pasquali98,Pasquali99a,Gelmini00a,Gelmini00b,%
Volkova01,Costa01a,Costa01b}.
The results shown in Figs. \ref{f:Nu_e} and \ref{f:Nu_m} are obtained
by using the two phenomenological approaches to the charm production
problem: the quark-gluon string model (QGSM) and recombination
quark-parton model (RQPM). The prompt muon fluxes predicted by QGSM
and RQPM are both consistent with the current deep underground data
and may be considered as the safe lower and upper limits for the
prompt muon contributions.
The basic assumptions of these models and the aspects of atmospheric
cascade calculations were described by \citet{Bugaev89b}.%
\footnote{See also Refs. \citep{Bugaev98,Naumov98a} for more details
          and related references.}

Calculations of conventional AN fluxes for energies below
$\sim100$ GeV are done with CORT for Kamioka site, while the
high-energy part of the conventional neutrino spectra as well as
the prompt neutrino contributions are obtained by using the results
of Ref. \citep{Naumov98b}. Needless to say the high-energy AN fluxes
are independent of location. The high-energy calculation takes into
account many ``thin'' effects, like $K_{e3}$ and $K_{\mu3}$ form
factors and $K_S^0$ semileptonic decays, meson regeneration and
charge exchange through reactions
$\pi^\pm+\mathrm{Air}\to\pi^{\pm(\mp)}+X$,
$\pi^\pm+\mathrm{Air}\to K^{\pm(0)}+X$, etc., as well as through
$K_{\pi2}$ and $K_{\pi3}$ decays.
The low- and high-energy parts of the spectra are smoothly merged
by using a polynomial (over $\log(E_\nu)$) clutching functions.
The results of these calculations for all underground detectors
are now embedded in a simple FORTRAN code ``CORTout''. It is based
on a two-dimensional spline interpolation over the detailed tables
and it may be useful for a fast evaluating the energy spectra and
zenith-angle distributions of conventional and prompt atmospheric
neutrinos within the energy range from 0.05 to $\sim10^{10}$ GeV.
The code is available upon request from the author.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\protect\section{Conclusions}
\label{sec:Conclusions}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Let us briefly re-state the main conclusions which follow from the
outcome of \citet{Fiorentini01a,Fiorentini01b} and from the present
study.

The results of extensive calculations performed with the new version
of code CORT by using the updated KM+SS hadronic interaction model
and the BESS+JACEE fit for the spectrum and composition of primary
cosmic rays are in substantial agreement with a quite representative
set of data on muon momentum spectra and charge ratios measured at
different atmospheric depths, zenith angles, and geomagnetic
locations.
This provides a solid evidence for the validity of our description
of hadronic interactions and shower development.
The low-energy atmospheric neutrino fluxes calculated with CORT for
many geomagnetic locations are systematically lower than those used
in the current analyses of the data from underground neutrino
experiments, while the neutrino flavor ratio and ``up-to-down''
asymmetry are essentially the same as in all recent calculations.
It is difficult to increase the AN flux without spoiling the
agreement with the current data on hadronic interactions, primary
spectrum and muon fluxes.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\protect\section*{Acknowledgments}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{acknowledgements}
It is pleasure to thank Gianni Fiorentini and Francesco Villante
for their cooperation and for many constructive discussions. I thank
Kostya Kuzmin and Segey Sinegovsky for significant contribution.
I am thankful to Mirko Boezio, Patricia Hansen, Pierre Le Coultre,
and Tomonori Wada for providing me with the data from
CAPRICE\,94/98, L3+C, and Okayama experiments before publication.
I also thank Thomas Hebbeker and Charles Timmermans for their
kind assistance in finding and understanding many experimental
data. This research was supported in part by the Ministry of
Education of the Russian Federation, the program ``Universities of
Russia'' and by the Florence L3+C group.
\end{acknowledgements}


\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
\begin{figure*}[t]
\centering
\ifpdf
\includegraphics[width=0.88\linewidth]{e1t.pdf} % for PDFLaTeX
\includegraphics[width=0.88\linewidth]{e2t.pdf} % for PDFLaTeX
\else
\includegraphics[width=0.88\linewidth]{e1t.eps} % for LaTeX
\includegraphics[width=0.88\linewidth]{e2t.eps} % for LaTeX
\fi
\protect\caption{Differential energy spectra of electron neutrinos
                 and antineutrinos for 11 zenith angles $\theta$.
                 Low-energy range is for Kamioka site. Contributions
                 from charm decay are calculated using recombination
                 quark-parton model (RQPM) and quark-gluon string
                 model (QGSM). At high energies, from smallest to
                 largest fluxes, $\cos\theta$ varies from 0 to 1
                 with an increment of 0.1 for each group of curves.
\label{f:Nu_e}}
\end{figure*}
\begin{figure*}[t]
\centering
\ifpdf
\includegraphics[width=0.88\linewidth]{m1t.pdf} % for PDFLaTeX
\includegraphics[width=0.88\linewidth]{m2t.pdf} % for PDFLaTeX
\else
\includegraphics[width=0.88\linewidth]{m1t.eps} % for LaTeX
\includegraphics[width=0.88\linewidth]{m2t.eps} % for LaTeX
\fi
\protect\caption{Differential energy spectra of muon neutrinos
                 and antineutrinos for 11 zenith angles $\theta$.
                 Low-energy range is for Kamioka site. Contributions
                 from charm decay are calculated using recombination
                 quark-parton model (RQPM) and quark-gluon string
                 model (QGSM). At high energies, from smallest to
                 largest fluxes, $\cos\theta$ varies from 0 to 1
                 with an increment of 0.1 for each group of curves.
\label{f:Nu_m}}
\end{figure*}


\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%\vspace*{3mm}
\begin{thebibliography}{999}
%\vspace*{3mm}

%-------------------------------------------------------------------%
 \bibitem[Abdel-Monem {\em et al.}(1973)]{Abdel-Monem73}
          Abdel-Monem, M.S. {\em et al.},
%%%%%%%%% Abdel-Monem, M.S., Benbrook, J.R., Osborne, A.R., and
%%%%%%%%% Sheldon, W.R.,
          in: Proc. of the 13th Intern. Cosmic Ray Conf., Denver,
          1973,
%%%%%%%%% ed. by Chasson, R.L. (University of Denver, Colorado
%%%%%%%%% Associated University Press, Boulder, CO, 1973),
          Vol.\,{\bf3}, p.\,1811. %pp.\,1811--?
%%CITATION = NONE;%%
%--------------------------------------------------------------------%
 \bibitem[Agrawal {\em et al.}(1996)]{Agrawal96}
          Agrawal, V. {\em et al.},
%%%%%%%%%``Atmospheric neutrino flux above 1 GeV'',
%%%%%%%%% Agrawal, V., Gaisser, T.K., Lipari, P., and Stanev, T.
          Phys. Rev. D\,{\bf53}, 1314, 1996. %1314--1323
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Aharonian {\em et al.}(1999)]{Aharonian99}
          Aharonian, F. {\em et al.},
%%%%%%%%% Aharonian, F., Akhperjanian, A.G., Barrio, J.A.,
%%%%%%%%% Belgarian, A.S., {\em et al.} (HEGRA Collaboration),
%%%%%%%%%``The cosmic ray proton spectrum determined with the imaging
%%%%%%%%% Atmospheric Cherenkov-technique'',
          Phys. Rev. D\,{\bf59}, 092003, 1999.
%%%%%%%%% (astro-.
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Alcaraz {\em et al.}(2000a)]{Alcaraz00a}
          Alcaraz, J. {\em et al.}
%%%%%%%%% Alcaraz, J., Alvisi, D., Ambrosi, B., Anderhub, H.,
%%%%%%%%% {\em et al.} (AMS Collaboration),
%%%%%%%%%``Protons in near Earth orbit'',
          Phys. Lett. B\,{\bf472}, 215, 2000. %215--226
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Alcaraz {\em et al.}(2000b)]{Alcaraz00b}
          Alcaraz, J. {\em et al.}
%%%%%%%%% Alcaraz, J., Alvisi, D., Ambrosi, B., Anderhub, H.,
%%%%%%%%% {\em et al.} (AMS Collaboration),
%%%%%%%%%``Cosmic protons'',
          Phys. Lett. B\,{\bf490}, 27, 2000. %27-35
%%CITATION = PHLTA,B490,27;%%
%-------------------------------------------------------------------%
 \bibitem[Alcaraz {\em et al.}(2000c)]{Alcaraz00c}
          Alcaraz, J. {\em et al.}
%%%%%%%%% Alcaraz, J., Alvisi, D., Ambrosi, B., Anderhub, H.,
%%%%%%%%% {\em et al.} (AMS Collaboration),
%%%%%%%%%``Helium in near Earth orbit'',
          Phys. Lett. B\,{\bf494}, 193, 2000. %193-202
%%CITATION = PHLTA,B494,193;%%
%-------------------------------------------------------------------%
 \bibitem[Allkofer {\em et al.}(1968)]{Allkofer68}
          Allkofer, O.C., Andresen, R.D., and Dau, W.D.,
%%%%%%%%%``The muon spectra near the geomagnetic equator'',
%%%%%%%%% in: Proc. of the 10th Intern. Cosmic Ray Conf., Calgary,
%%%%%%%%% 1967, ed. by Wilson, M.D. (National Research Council of
%%%%%%%%% Canada, Ottawa),
          Can. J. Phys. {\bf46}, S301, 1968. %S301--S305
%%CITATION = CJPHA,46,S301;%%
%-------------------------------------------------------------------%
 \bibitem[Allkofer \& Clausen(1970)]{Allkofer70a}
          Allkofer, O.C. and Clausen, K.,
%%%%%%%%% in: Proc. of the 11th Intern. Cosmic Ray Conf., Budapest,
%%%%%%%%% 1969, ed. by Gemsey, T. {\em et al.},
          Acta. Phys. Acad. Sci. Hung. {\bf29} (Suppl. 2), 698,
          1970. %698--?
%%CITATION = APAHA,29,698;%%
%-------------------------------------------------------------------%
 \bibitem[Allkofer {\em et al.}(1970)]{Allkofer70b}
          Allkofer, O.C., Dau, W.D., and Jokisch, H.,
%%%%%%%%%``The absolute vertical muon intensity of cosmic rays at
%%%%%%%%% 1 GeV/c at sea level'',
          Phys. Lett. {\bf31}\,B, 606, 1970; %606-608
%%CITATION = PHLTA,B31,606;%%
          {\em ibid.}, {\bf36}\,B, 428, 1971 (Erratum).
%%CITATION = PHLTA,B36,428;%%
%-------------------------------------------------------------------%
 \bibitem[Allkofer {\em et al.}(1971)]{Allkofer71}
          Allkofer, O.C., Carstensen, K., and Dau, W.D.,
%%%%%%%%%``The absolute cosmic ray muon spectrum at sea level'',
          Phys. Lett. {\bf36}\,B, 425, 1971. %425--427
%%CITATION = PHLTA,B36,425;%%
%-------------------------------------------------------------------%
 \bibitem[Allkofer \& Dau(1972)]{Allkofer72}
          Allkofer, O.C. and Dau, W.D.,
%%%%%%%%%``The muon charge ratio at sea level in the low momentum
%%%%%%%%% region'',
          Phys. Lett. {\bf38}\,B, 439, 1972. %439--440
%%CITATION = PHLTA,B38,439;%%
%-------------------------------------------------------------------%
 \bibitem[Allkofer {\em et al.}(1975)]{Allkofer75}
          Allkofer, O.C., Clausen, K., and Dau, W.D.,
%%%%%%%%%``The low-momentum muon spectrum near equator'',
          Lett. Nuovo Cim. {\bf12}, 107, 1975. %107--110
%%CITATION = NCLTA,12,107;%%
%-------------------------------------------------------------------%
 \bibitem[Appleton {\em et al.}(1971)]{Appleton71}
          Appleton, I.C., Hogue, M.T., and Rastin, B.C.,
%%%%%%%%%``A study of the muon momentum spectrum and
%%%%%%%%% positive--negative ratio at sea--level'',
          Nucl. Phys. B\,{\bf26}, 365, 1971. %365--389
%%CITATION = NUPHA,B26,365;%%
%--------------------------------------------------------------------%
 \bibitem[Asakimori {\em et al.}(1998)]{Asakimori98}
          Asakimori, K. {\em et al.},
%%%%%%%%% Asakimori, K., Burnett, T.H., Cherry, M.L., Chevli, K.,
%%%%%%%%% {\em et al.} (JACEE Collaboration),
%%%%%%%%%``Cosmic-ray proton and helium spectra: results from the
%%%%%%%%% JACEE experiment'',
          Astrophys. J. {\bf502}, 278, 1998. %278--283
%%CITATION = ASJOA,502,278;%%
%-------------------------------------------------------------------%
 \bibitem[Athar {\em et al.}(2001)]{Athar01}
          Athar, H. {\em et al.},
%%%%%%%%% Athar, H., Cheung, K.M., Lin, G.L., and Tseng, J.J.,
%%%%%%%%%``High-energy cosmic-ray tau neutrino flux from the
%%%%%%%%% Milky Way,'' %(v.1)
%%%%%%%%%``Comparison of high-energy galactic and atmospheric tau
%%%%%%%%% neutrino flux,'' %(v.2)
          .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Aurela \& Wolfendale(1967)]{Aurela67}
          Aurela, A.M. and Wolfendale, A.W.,
          Annales Acad. Sci. Fenn. A\,{\bf6}, Vol.\,{\bf227}, 1,
          1967.
%%CITATION = AAFPA,227,1;%%
%-------------------------------------------------------------------%
 \bibitem[Ayre {\em et al.}(1973)]{Ayre73}
          Ayre, C.A. {\em et al.},
%%%%%%%%% Ayre, C.A., Baxendale, J.M., Daniel, B.J., Hume, C.J.,
%%%%%%%%% Nandi, B.C., Thompson, M.G., Wells, S.C., Whalley, M.R.,
%%%%%%%%% and Wolfendale, A.W.,
%%%%%%%%%``The muon charge ratio in the near vertical direction in
%%%%%%%%% the range 10 GeV/c -- 450 GeV/c'',
          in: Proc. of the 13th Intern. Cosmic Ray Conf., Denver,
          1973,
%%%%%%%%% ed. by Chasson, R.L. (University of Denver, Colorado
%%%%%%%%% Associated University Press, Boulder, CO, 1973),
          Vol.\,{\bf3}, p.\,1822. %pp.\,1822--1827
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Ayre {\em et al.}(1975)]{Ayre75}
          Ayre, C.A. {\em et al.},
%%%%%%%%% Ayre, C.A., Baxendale, J.M., Hume, C.J., {\em et al.},
%%%%%%%%%``Precise measurement of the vertical muon spectrum in the
%%%%%%%%% range 20--500 GeV/c'',
          J. Phys. G: Nucl. Phys. {\bf1}, 584, 1975. %584--600
%%CITATION = JPHGB,1,584;%%
%-------------------------------------------------------------------%
 \bibitem[Baber {\em et al.}(1968a)]{Baber68a}
          Baber, S.R., Nash, W.F., and Rastin, B.C.,
%%%%%%%%%``The momentum spectrum of muons at sea level in the range
%%%%%%%%% 3--1000 GeV/c'',
          Nucl. Phys. B\,{\bf4}, 539, 1968. %539--548
%%CITATION = NUPHA,B4,539;%%
%-------------------------------------------------------------------%
 \bibitem[Baber {\em et al.}(1968b)]{Baber68b}
          Baber, S.R., Nash, W.F., and Rastin, B.C.,
%%%%%%%%%``The positive-negative ratio of muons at sea level in the
%%%%%%%%% range 3--360 GeV'',
          Nucl. Phys. B\,{\bf4}, 549, 1968. %549--558
%%CITATION = NUPHA,B4,549;%%
%-------------------------------------------------------------------%
 \bibitem[Barr {\em et al.}(1989)]{Barr89}
          Barr, G., Gaisser, T.K., and Stanev, T.,
%%%%%%%%%``Flux of atmospheric neutrinos'',
          Phys. Rev. D\,{\bf39}, 3532, 1989. %3532--3534
%%CITATION = PHRVA,D39,3532;%%
%-------------------------------------------------------------------%
 \bibitem[Baschiera {\em et al.}(1979)]{Baschiera79}
          Baschiera, B. {\em et al.},
%%%%%%%%% Baschiera, B., Basini, G., Bilokon, H., Dettorre
%%%%%%%%% Piazzoli, B., Mannocchi, G., Castagnoli, C., and
%%%%%%%%% Picchi, P.,
%%%%%%%%%``Integral and differential absolute intensity measurements
%%%%%%%%% of cosmic-ray muons below 1 GeV'',
          Nuovo Cim. C\,{\bf2}, 473, 1979. %473--487
%%CITATION = NUCIA,C2,473;%%
%-------------------------------------------------------------------%
 \bibitem[Bateman {\em et al.}(1971)]{Bateman71}
          Bateman, B.J. {\em et al.},
%%%%%%%%% Bateman, B.J., Cantrell, W.G., Durda, D.R., Duller, N.M.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Absolute measurement of the vertical cosmic ray muon
%%%%%%%%% intensity at 3--50 GeV/c near sea level'',
          Phys. Lett. {\bf36}\,B, 144, 1971. %144--148
%%CITATION = PHLTA,B36,144;%%
%-------------------------------------------------------------------%
 \bibitem[Battistoni {\em et al.}(2000)]{Battistoni00}
          Battistoni, G. {\em et al.},
%%%%%%%%% Battistoni, G., Ferrari, A., Lipari, P., Montaruli, T.,
%%%%%%%%% Sala, P.R., and Rancati, T.,
%%%%%%%%%``A 3-dimensional calculation of atmospheric neutrino
%%%%%%%%% flux'',
          Astropart. Phys. {\bf12}, 315, 2000. %315--333
%%%%%%%%% .
          For the current FLUKA results see URL \FLUKA-3D.
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Battistoni {\em et al.}(2001a)]{Battistoni01a}
          Battistoni, G. {\em et al.},
%%%%%%%%% Battistoni, G., Ferrari, A., Montaruli, T., and Sala, P.R.,
%%%%%%%%%``Comparison of the FLUKA calculations with CAPRICE94 data
%%%%%%%%% on muons in  atmosphere'',
           % (Submitted to Astropart. Phys.).
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Battistoni {\em et al.}(2001b)]{Battistoni01b}
          Battistoni, G. {\em et al.},
%%%%%%%%% Battistoni, G., Ferrari, A., Montaruli, T., and Sala, P.R.,
%%%%%%%%%``Progresses in the validation of the FLUKA atmospheric
%%%%%%%%% $\nu$ flux calculations'',
          .
%%CITATION = ;%%
%--------------------------------------------------------------------%
 \bibitem[Baxendale {\em et al.}(1975)]{Baxendale75}
          Baxendale, J.M., Hume, C.J., and Thompson, M.G.,
%%%%%%%%%``Precise measurement of the sea level muon charge ratio'',
          J. Phys. G: Nucl. Phys. {\bf1}, 781, 1975. %781--788
%%CITATION = JPHGB,1,781;%%
%--------------------------------------------------------------------%
 \bibitem[Beatty {\em et al.}(1993)]{Beatty93}
          Beatty, J.J. {\em et al.},
%%%%%%%%% Beatty, J.J., Ficenec, D.J., Tobias, S., Mitchell, J.W.,
%%%%%%%%% {\em et al.}, % [SMILI]
%%%%%%%%%``The cosmic-ray ${}^3$He/${}^4$He ratio from 100 to 1600
%%%%%%%%% MeV amu${}^{-1}$",
          Astrophys. J. {\bf413}, 268, 1993. %268--280
%%CITATION = ASJOA,518,457;%%
%-------------------------------------------------------------------%
 \bibitem[Bellotti {\em et al.}(1999)]{Bellotti99}
          Bellotti, R. {\em et al.},
%%%%%%%%% Bellotti, R., Cafagna, F., Circella, M., De Marzo, C.N.,
%%%%%%%%% {\em et al.} (WiZard/MASS2 Collaboration),
%%%%%%%%%``Balloon measurements of cosmic ray muon spectra in the
%%%%%%%%% atmosphere along with those of primary protons and helium
%%%%%%%%% nuclei over midlatitude'',
          Phys. Rev. D\,{\bf60}, 052002, 1999.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Bialas {\em et al.}(1976)]{Bialas76}
          Bialas, A., Bleszy\'nski, M., and Czy\.z, W.,
%%%%%%%%%``Multiplicity distributions in nucleus-nucleus collisions
%%%%%%%%% at high-energies'',
          Nucl. Phys. B\,{\bf111}, 461, 1976. %461--476
%%CITATION = NUPHA,B111,461;%%
%-------------------------------------------------------------------%
 \bibitem[Boezio {\em et al.}(1999a)]{Boezio99a}
          Boezio, M. {\em et al.},
%%%%%%%%% Boezio, M., Carlson, P., Francke, T., Weber, N.,
%%%%%%%%% {\em et al.} (WiZard/CAPRICE Collaboration),
%%%%%%%%%``The cosmic ray proton and helium spectra between 0.4 and
%%%%%%%%% 200 GV'',
          Astrophys. J. {\bf518}, 457, 1999. %457--472
%%CITATION = ASJOA,518,457;%%
%-------------------------------------------------------------------%
 \bibitem[Boezio {\em et al.}(1999)]{Boezio99b}
          Boezio, M. {\em et al.},
%%%%%%%%%``New measurement of the flux of atmospheric muons'',
          Phys. Rev. Lett. {\bf82}, 4757, 1999. %4757--4760
%%CITATION = PRLTA,82,4757;%%
%-------------------------------------------------------------------%
 \bibitem[Boezio {\em et al.}(2000)]{Boezio00}
          Boezio, M. {\em et al.},
%%%%%%%%%``Measurement of the flux of atmospheric muons with the
%%%%%%%%% CAPRICE\,94 apparatus'',
          Phys. Rev. D\,{\bf62}, 032007, 2000;
%%%%%%%%% .
          Boesio, M., private communication (September, 2001).
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Brancus {\em et al.}(2000)]{Brancus00}
          Brancus, I.M. {\em et al.},
%%%%%%%%% Brancus, I.M., Vulpescu, B., Bercuci, A., Badea, A.F.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Measurements of the charge ratio of atmospheric muons'',
          Acta Phys. Polon. B\,{\bf31}, 465, 2000. %465-470
%%CITATION = APPOA,B31,465;%%
%-------------------------------------------------------------------%
 \bibitem[Buckley {\em et al.}(1994)]{Buckley94}
          Buckley, J. {\em et al.},
%%%%%%%%% Buckley, J., Dwyer, J., Muller, D., Swordy, S., and
%%%%%%%%% Tang, K.K. (RICH-I Collaboration)
%%%%%%%%%``A new measurement of the flux of the light cosmic-ray
%%%%%%%%% nuclei at high energies'',
          Astrophys. J. {\bf429}, 736, 1994. %736-747
%%CITATION = ASJOA,429,736;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1984)]{Bugaev84}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``The nucleon component of cosmic radiation in the
%%%%%%%%% atmosphere at intermediate energies'',
          Investigations on Geomagnetism, Aeronomy and Solar Physics,
          {\bf73}, 198, 1984. %198--211
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1985)]{Bugaev85}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``Cosmic-ray nucleons in the atmosphere'', [$\Pi$--0385]
%%%%%%%%%``Geomagnetic effect and spectra of secondary cosmic
%%%%%%%%% nucleons'', [$\Pi$--0401]
          Inst. for Nuclear Research, USSR Acad. of Sci.,
          Reports $\Pi$--0385 and $\Pi$--0401, Moscow, 1985.
%%%%%%%%% (unpublished).
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1986)]{Bugaev86}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``Low and medium energy neutrinos in the atmosphere'',
          Izv. Akad. Nauk SSSR, Ser. Fiz. {\bf50}, 2239, 1986
          %2239--2241
          [Bull. Russ. Acad. Sci. Phys. {\bf50,\,N.\,11}, 156, 1986].
          %156--158
%%CITATION = BUPSA,50N11,156;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1987a)]{Bugaev87a}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``Cosmic-ray muons and neutrinos at low and intermediate
%%%%%%%%% energies'',
          Yad. Fiz. {\bf45}, 1380, 1987 %1380--1391
          [Sov. J. Nucl. Phys. {\bf45}, 857, 1987]. %857--864
%%CITATION = SJNCA,45,857;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1987b)]{Bugaev87b}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``Geomagnetic effect and spectra of cosmic-ray muons and
%%%%%%%%% neutrinos'',
          Inst. for Nuclear Research, USSR Acad. of Sci.,
          Report $\Pi$-0537, Moscow, 1987.
%%%%%%%%% (unpublished).
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1988)]{Bugaev88}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``Influence of geomagnetic field on low energy atmospheric
%%%%%%%%% neutrino fluxes'',
          in: Proc. of the 2nd Intern. Symposium ``Underground
          Physics'87'', Baksan Valley, 1987, ed. by Domogatsky, G.V.
          {\em et al.} (Moscow, ``Nauka'', 1988) p.\,255.
          %pp.\,255--258
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1989)]{Bugaev89a}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``On the interpretation of the Kamiokande neutrino
%%%%%%%%% experiment'',
          Phys. Lett. B\,{\bf232}, 391, 1989. %391--397
%%CITATION = PHLTA,B232,391;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev {\em et al.}(1989)]{Bugaev89b}
          Bugaev, E.V. {\em et al.},
%%%%%%%%% Bugaev, E.V., Naumov, V.A., Sinegovsky, S.I., and
%%%%%%%%% Zaslavskaya, E.S.,
%%%%%%%%%``Prompt Leptons In Cosmic Rays'',
         Nuovo Cim. C.\,{\bf12}, 41, 1989. %41--73
%%CITATION = NUCIA,C12,41;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev and Naumov(1990)]{Bugaev90}
          Bugaev, E.V. and Naumov, V.A.,
%%%%%%%%%``On the interpretation of the Kamiokande neutrino
%%%%%%%%% experiment'',
          Yad. Fiz. {\bf51}, 774, 1990 %774--776
          [Sov. J. Nucl. Phys. {\bf51}, 493, 1990]. %493--494
%%CITATION = SJNCA,51,493;%%
%-------------------------------------------------------------------%
 \bibitem[Bugaev {\em et al.}(1998)]{Bugaev98}
          Bugaev, E.V. {\em et al.},
%%%%%%%%% Bugaev, E.V., Misaki, A., Naumov, V.A., Sinegovskaya,
%%%%%%%%% T.S., Sinegovsky, S.I., and Takahashi, N.,
%%%%%%%%%``Atmospheric muon flux at sea level, underground, and
%%%%%%%%% underwater'',
          Phys. Rev. D\,{\bf58}, 054001, 1998
          (see also  for more details).
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Burnett {\em et al.}(1973)]{Burnett73}
          Burnett, T.H. {\em et al.},
%%%%%%%%% Burnett, T.H., LaMay, L.J., Masek, G.E., and Maung, T.,
%%%%%%%%%``Measurement of cosmic-ray muon charge ratio at sea level
%%%%%%%%% between energies 10 and 15000 GeV'',
          Phys. Rev. Lett. {\bf30}, 937, 1973. %937--940
%%CITATION = PRLTA,30,937;%%
%-------------------------------------------------------------------%
 \bibitem[Caro {\em et al.}(1950)]{Caro50}
          Caro, D.E., Parry, J.K., and Rathgeber, H.D.,
          Nature {\bf165}, 689, 1950. %689--?
%%CITATION = NATUA,165,689;%%
%-------------------------------------------------------------------%
 \bibitem[Coates \&  Nash(1962)]{Coates62}
          Coates, D.W. and Nash, W.F.,
          Austral. J. Phys. {\bf15}, 420, 1962. %420--?
%%CITATION = AUJPA,15,420;%%
%-------------------------------------------------------------------%
 \bibitem[Conversi(1950)]{Conversi50}
          Conversi, M.,
%%%%%%%%%``Experiments on cosmic-ray mesons and protons at several
%%%%%%%%% altitudes and latitude'',
          Phys. Rev. {\bf79}, 749, 1950. %749--767
%%%%%%%%% Actually the point at E_kin=400 MeV are from a paper by
%%%%%%%%% M. Conversi and F. Nappo (measurements at 1947) cited
%%%%%%%%% by Conversi in p.758.
%%CITATION = PHRVA,79,749;%%
%-------------------------------------------------------------------%
\bibitem[Costa(2001)]{Costa01a}
         Costa, C.G.S.,
%%%%%%%%%``The prompt lepton cookbook,''
         Astropart. Phys. {\bf16}, 193, 2001. %193--204
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
\bibitem[Costa {\em et al.}(2001)]{Costa01b}
         Costa, C.G.S., Halzen, F., and Salles, C.,
%%%%%%%%%``The prompt TeV-PeV atmospheric neutrino window,''
         .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Coutu {\em et al.}(2000)]{Coutu00}
          Coutu, S. {\em et al.},
%%%%%%%%% Coutu, S., Beatti, J.J., Du Vernois, M.A., Barwik, S.W.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Energy spectra, altitude profiles and charge ratios of
%%%%%%%%% atmospheric muons'',
          Phys. Rev. D\,{\bf62}, 032001, 2000.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Diehl {\em et al.}(1997)]{Diehl97}
          Diehl, E. {\em et al.},
%%%%%%%%% Diehl, E., Ellithorpe, D., M\"{u}ller, D., and Swordy, S.
%%%%%%%%% (RICH-II Collaboration),
          in: Proc. of the 25th Intern. Cosmic Ray Conf., Durban,
%%%%%%%%% South Africa, July 30 - August 6,
          1997,
%%%%%%%%% ed. by Potgieter, M.S., Raubenheimer, B.C., and
%%%%%%%%% van der Walt, D.J. (Wesprint, Potchefstroom, Space Research
%%%%%%%%% Unit, 1997)
          Vol.\,{\bf3}, pp.\,405. %pp.\,405--408.
          No result presented in this paper. The data are taken from
          von Rosenvinge, T.T.
%%%%%%%%%``Low-energy cosmic rays: their composition and nuclear
%%%%%%%%% cross sections'',
          {\em ibid.}, Vol.\,{\bf8}, p.237. %pp.\,237--256
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Dorman {\em et al.}(1971)]{Dorman71}
          Dorman, L.I., Smirnov, V.S., and Tiasto, M.I.,
          {\em Cosmic rays in the magnetic field of the Earth},
          ``Nauka'', Moscow, 1971.
%%CITATION = NONE;%%
%--------------------------------------------------------------------%
 \bibitem[Dorman(1972)]{Dorman72}
          Dorman, L.I.,
          {\em Meteorological effects of cosmic rays},
          ``Nauka'', Moscow, 1972.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Engel {\em et al.}(2001)]{Engel01}
          Engel, R., Gaisser, T.K., and Stanev, T.,
%%%%%%%%%``The flux of atmospheric muons'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,1029. %pp.\,1029--1032
%%CITATION = NONE;%%
%--------------------------------------------------------------------%
 \bibitem[Filosofo(1954)]{Filosofo54}
          Filosofo, I., Pohl, E., and Pohl-R\"uling, J.,
%%%%%%%%%``On the positive excess of mesons of high energy'',
          Nuovo Cim. {\bf12}, 809, 1954. %809--812
%%CITATION = NUCIA,12,809;%%
%-------------------------------------------------------------------%
 \bibitem[Fiorentini {\em et al.}(2001a)]{Fiorentini01a}
          Fiorentini, G., Naumov, V.A., and Villante, F.L.,
%%%%%%%%%``Atmospheric neutrino flux supported by recent muon
%%%%%%%%% experiments'',
          Phys. Lett. B\,{\bf510}, 173, 2001. %173--188
%%%%%%%%% (see also .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Fiorentini {\em et al.}(2001b)]{Fiorentini01b}
          Fiorentini, G., Naumov, V.A., and Villante, F.L.,
%%%%%%%%%``Atmospheric neutrino flux and muon data'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,1218 %pp.\,1218--1221
          .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Fukui(1955)]{Fukui55}
          Fukui, S.,
          J. Phys. Soc. Jap. {\bf10}, 735, 1955. %735--?
%%CITATION = JUPSA,10,735;%%
%--------------------------------------------------------------------%
 \bibitem[Gaisser {\em et al.}(1996)]{Gaisser96}
          Gaisser, T.K. {\em et al.},
%%%%%%%%% Gaisser, T.K., Honda, M., Kasahara, K., Lee, H.,
%%%%%%%%% Midorikawa, S., Naumov, V.A., and Stanev, T.,
%%%%%%%%%``Comparison of atmospheric neutrino flux calculations at
%%%%%%%%% low energies'',
          Phys. Rev. D\,{\bf54}, 5578, 1996. %5578--5584.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Gelmini {\em et al.}(2000a)]{Gelmini00a}
          Gelmini, G., Gondolo, P., and Varieschi, G.,
%%%%%%%%%``Prompt atmospheric neutrinos and muons: NLO vs LO {QCD}
%%%%%%%%% predictions,''
          Phys. Rev. D\,{\bf61}, 036005, 2000.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Gelmini {\em et al.}(2000b)]{Gelmini00b}
          Gelmini, G., Gondolo, P., and Varieschi, G.,
%%%%%%%%%``Prompt atmospheric neutrinos and muons: Dependence on
%%%%%%%%% the gluon  distribution function,''
          Phys. Rev. D\,{\bf61}, 056011, 2000.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Grandegger(1993)]{Grandegger93}
          Grandegger, W.,
          KfK Report 5122, Kernforschungszentrum Karlsruhe, 1993.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Green {\em et al.}(1979)]{Green79}
          Green, P.J. {\em et al.},
%%%%%%%%% Green, P.J., Duller, N.M., Magnuson, C.E., Choate, L.M.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Absolute intensities of medium-energy muons in the
%%%%%%%%% vertical and at zenith angles $55^\circ-85^\circ$'',
          Phys. Rev. D\,{\bf20}, 1598, 1979. %1598--1607
%%CITATION = PHRVA,D20,1598;%%
%-------------------------------------------------------------------%
 \bibitem[Grieder(2001)]{Grieder01}
          Grieder, P.K.F., {\em Cosmic rays at Earth},
          North-Holland, Elsevier Science Publishers B.\,V.,
          Amsterdam, 2001.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Hansen {\em et al.}(2001)]{Hansen01}
          Hansen, P. {\em et al.},
%%%%%%%%% Hansen, P., Ambriola, M., Bartalucci, S., Bellotti, R.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``A new measurement of muon spectra in the atmosphere'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf1}, p.\,921; %pp.\,921--924
          Hansen, P., private communication (September, 2001).
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Hayman \& Wolfendale(1962a)]{Hayman62a}
          Hayman, P.J. and Wolfendale, A.W.,
%%%%%%%%%``The momentum spectrum of cosmic ray muons near sea level
%%%%%%%%% in the momentum range 5-1000 GeV/c'',
          Proc. Phys. Soc. (London) A\,{\bf80}, 710, 1962. %710--728
%%CITATION = PPSOA,A80,710;%%
%-------------------------------------------------------------------%
 \bibitem[Hayman \& Wolfendale(1962b)]{Hayman62b}
          Hayman, P.J. and Wolfendale, A.W.,
          Nature {\bf195}, 166, 1962. %166--?
%%CITATION = NATUA,195,166;%%
%-------------------------------------------------------------------%
 \bibitem[Hebbeker \& Timmermans(2001)]{Hebbeker01}
          Hebbeker, T. and Timmermans, C.,
%%%%%%%%%``A compilation of high-energy atmospheric muon data at
%%%%%%%%% sea level'',
           (submitted to Astropart. Phys.).
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Holmes {\em et al.}(1961)]{Holmes61}
          Holmes, J.E.R., Owen, B.G., and Rodgers, A.L.,
%%%%%%%%%``Measurements of the momentum spectrum and charge ratio
%%%%%%%%% of cosmic ray muons at sea level in the momentum range
%%%%%%%%% 10 GeV/c -- 1000 GeV/c'',
          Proc. Phys. Soc. (London) A\,{\bf78}, 505, 1961. %505--515
%%CITATION = PPSOA,A78,505;%%
%-------------------------------------------------------------------%
 \bibitem[Honda {\em et al.}(1990)]{Honda90}
          Honda, M. {\em et al.},
%%%%%%%%% Honda, M., Kasahara, K., Hidaka, K., and Midorikawa, S.,
%%%%%%%%%``Atmospheric neutrino fluxes'',
          Phys. Lett. B\,{\bf248}, 193, 1990. %193--198
%%CITATION = PHLTA,B248,193;%%
%-------------------------------------------------------------------%
 \bibitem[Honda {\em et al.}(1995)]{Honda95}
          Honda, M. {\em et al.},
%%%%%%%%% Honda, M., Kajita, T., Kasahara, K., and Midorikawa, S.,
%%%%%%%%%``Calculation of the flux of atmospheric neutrinos'',
          Phys. Rev. D\,{\bf52}, 4985, 1995. %4985--5005
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Honda {\em et al.}(1996)]{Honda96}
          Honda, M. {\em et al.},
%%%%%%%%% Honda, M., Kajita, T., Kasahara, K., and Midorikawa, S.,
%%%%%%%%%``Atmospheric neutrinos'',
          Prog. Theor. Phys. Suppl. {\bf123}, 483, 1996. %483--490
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Honda {\em et al.}(2001a)]{Honda01a}
          Honda, M. {\em et al.},
%%%%%%%%% Honda, M., Kajita, T., Kasahara, K., and Midorikawa, S.,
%%%%%%%%%``Comparison of 3-dimensional and 1-dimensional schemes in
%%%%%%%%% the calculation of atmospheric neutrinos'',
          Phys. Rev. D\,{\bf64}, 053011, 2001.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Honda {\em et al.}(2001b)]{Honda01b}
          Honda, M. {\em et al.},
%%%%%%%%% Honda, M., Kajita, T., Kasahara, K., and Midorikawa, S.,
%%%%%%%%%``Calculation of the atmospheric neutrino flux improved
%%%%%%%%% with recent cosmic ray observations'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,1162. %pp.\,1162--1165
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Ichimura {\em et al.}(1993)]{Ichimura93}
          Ichimura, M. {\em et al.},
%%%%%%%%% Ichimura, M., Kogawa, M., Kuramata, S., Mito, H.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Observation of heavy cosmic-ray primaries over the wide
%%%%%%%%% energy range from $\sim100$ GeV/particle to $\sim100$
%%%%%%%%% TeV/particle: Is the celebrated ``knee'' actually so
%%%%%%%%% prominent?'',
          Phys. Rev. D\,{\bf48}, 1949, 1993. %949--1975
%%CITATION = PHRVA,D48,1949;%%
%-------------------------------------------------------------------%
 \bibitem[Jannakos(1995)]{Jannakos95}
          Jannakos, T.E.,
%%%%%%%%%``Muon capture reactions on C-12 and C-13 in the KARMEN
%%%%%%%%% detector'',
%%%%%%%%%``Myoneinfang-Reaktionen an 12C und 13C im KARMEN
%%%%%%%%% Detektor'',
          FZKA Report 5520, Forschungszentrum Karlsruhe, 1995.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Jokisch {\em et al.}(1979)]{Jokisch79}
          Jokisch, H. {\em et al.}
%%%%%%%%% Jokisch, H., Carstensen, K., Dau, W.D., Meyer, H.J., and
%%%%%%%%% Allkofer, O.C.,
%%%%%%%%%``Cosmic-ray muon spectrum up to 1 TeV at $75^\circ$ zenith
%%%%%%%%% angle'',
          Phys. Rev. D\,{\bf19}, 1368, 1979. %1368--1372
%%CITATION = PHRVA,D19,1368;%%
%--------------------------------------------------------------------%
 \bibitem[Kalinovsky {\em et al.}(1985)]{Kalinovsky85}
          Kalinovsky, A.N., Mokhov, N.V., and Nikitin, Yu.P.,
          {\em Transport of high-energy particles through matter},
          ``Energoatomizdat'', Moscow, 1985.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Kawaguchi {\em et al.}(1965)]{Kawaguchi65}
          Kawaguchi, S. {\em et al.},
%%%%%%%%% Kawaguchi, S., Sakai, T., Oda, H., Ueno, H., and
%%%%%%%%% Kamiya, Y.,
          in: Proc. of the 9th Intern. Cosmic Ray Conf., London,
          1965,
%%%%%%%%% (The Institute of Physics and The Physical Society,
%%%%%%%%% London, 1965),
          Vol.\,{\bf2}, p.\,941. %pp.\,941--?
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Kellogg {\em et al.}(1978)]{Kellogg78}
          Kellogg, R.G., Kasha, H., and Larsen, R.C.,
%%%%%%%%%``Momentum spectra, charge ratio, and zenith-angle
%%%%%%%%% dependence of cosmic-ray muons'',
          Phys. Rev. D\,{\bf17}, 98, 1978. %98-113
%%CITATION = PHRVA,D17,98;%%
%--------------------------------------------------------------------%
 \bibitem[Kimel' \& Mokhov(1974)]{Kimel'74}
          Kimel', L.R. and Mokhov, N.V.,
%%%%%%%%%``Distributions of particles in the diapason of
%%%%%%%%% $10^{-2} -10^{12}$ eV, initiated in dense media by
%%%%%%%%% high-energy hadrons'',
          Izv. Vuz. Fiz. {\bf10}, 17, 1974. %17--23
%%CITATION = IVUFA,10,17;%%
%--------------------------------------------------------------------%
 \bibitem[Kimel' \& Mokhov(1975)]{Kimel'75}
          Kimel', L.R. and Mokhov, N.V.,
%%%%%%%%%``Differential cross sections of hadron-nucleus interactions
%%%%%%%%% and some results of calculations of inter-nuclei cascades'',
          Problems of Dosimetry and Radiation Protection, {\bf14}, 41,
          1975. %41--44
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Kowalski \& Bartke(1981)]{Kowalski81}
          Kowalski, M. and Bartke, J.,
%%%%%%%%%``Probability distribution of the number of interacting
%%%%%%%%% nucleons of the incident nucleus in collisions of
%%%%%%%%% relativistic nuclei'',
          Acta Phys. Polon. B\,{\bf12}, 759, 1981. %759--?.
%%CITATION = APPOA,B12,759;%%
%-------------------------------------------------------------------%
 \bibitem[Kremer {\em et al.}(1999)]{Kremer99}
          Kremer, J. {\em et al.},
%%%%%%%%% Kremer, J., Boezio, M.,  Ambriola, M.L., Barbiellini, G.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Measurements of ground-level muons at two geomagnetic
%%%%%%%%% locations'',
          Phys. Rev. Lett. {\bf83}, 4241, 1999. %4241--4244
%%CITATION = PRLTA,83,4241;%%
%-------------------------------------------------------------------%
 \bibitem[Krizmanic {\em et al.}(1995)]{Krizmanic95}
          Krizmanic, J.F. {\em et al.},
%%%%%%%%% Krizmanic, J.F., Davis, A.J., Barbier, L.M., Christian,
%%%%%%%%% E.R., {\em et al.},
%%%%%%%%%``Measurement of the cosmic ray induced muon flux through
%%%%%%%%% the atmosphere using IMAX'',
          in: Proc. of the 24th Intern. Cosmic Ray Conf., Rome,
          1995,
%%%%%%%%% (ARTI Grafiche Editoriali SRL, Urdino, 1995),
          Vol.\,{\bf1}, p.\,593. %pp.\,593--596
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Ladron de Guevara(2001)]{Ladron01}
          Ladron de Guevara, P. (for the L3 Collaboration),
%%%%%%%%%``First results from the L3+C experiment at CERN'',
          Nucl. Phys. B (Proc. Suppl.) {\bf95}, 237, 2001. %237--240
%%CITATION = NUPHZ,95,237;%%
%-------------------------------------------------------------------%
 \bibitem[Le Coultre(2001)]{LeCoultre01}
          Le Coultre, P. (for the L3+C Collaboration),
%%%%%%%%%``Momentum dependence of the vertical muon flux and charge
%%%%%%%%% ratio measured by the L3+C experiment'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,974; %pp.\,974--976
          Le Coultre, P., private communication (August, 2001)
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Lipari {\em et al.}(1998)]{Lipari98}
          Lipari, P., Gaisser, T.K., and Stanev, T.,
%%%%%%%%%``Geomagnetic effects on atmospheric neutrinos'',
          Phys. Rev. D\,{\bf58}, 073003, 1998.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Liu {\em et al.}(2001)]{Liu01}
          Liu, Y., Derome, L., and Bu\'enerd, M.,
%%%%%%%%%``Atmospheric muon and neutrino flux from 3D simulation'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,1033. %pp.\,1033--1036
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Mason(1972)]{Mason72}
          Mason, G.M.,
%%%%%%%%%``Interstellar propagation of galactic cosmic-ray nuclei
%%%%%%%%% $1 \leq Z \leq8$ in the energy range 100--1000 MeV per
%%%%%%%%% nucleon'', %[IMP-5, 1969]
          Astrophys. J. {\bf171}, 139, 1972. %139--161
%%CITATION = ASJOA,171,139;%%
%-------------------------------------------------------------------%
 \bibitem[Menn {\em et al.}(2000)]{Menn00}
          Menn, W. {\em et al.},
%%%%%%%%% Menn, W., Hof, M., Reimer, O., Simon, M., {\em et al.}
%%%%%%%%% (IMAX Collaboration),
          Astrophys. J. {\bf533}, 281, 2000. %281--297
%%CITATION = ASJOA,533,281;%%
%-------------------------------------------------------------------%
 \bibitem[Moroney \& Parry(1954)]{Moroney54}
          Moroney, J.R. and Parry, J.K.,
          Austral. J. Phys. {\bf7}, 424, 1954. %424--?
%%CITATION = AUJPA,7,424;%%
%-------------------------------------------------------------------%
 \bibitem[Motoki {\em et al.}(2001)]{Motoki01}
          Motoki, M. {\em et al.},
%%%%%%%%% Motoki, M., Sanuki, T., Orito, S., Abe, K. {\em et al.}
%%%%%%%%% (BESS Collaboration),
%%%%%%%%%``Precise measurement of atmospheric muon fluxes at sea
%%%%%%%%% level'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf1}, p.\,927. %pp.\,927--930
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Nandi \& Sinha(1972a)]{Nandi72a}
          Nandi, B.C. and Sinha, M.S.,
%%%%%%%%%``The momentum spectrum of muons at sea level in the range
%%%%%%%%% 5--1200 GeV/c'',
          J. Phys. A: Gen. Phys. {\bf5}, 1384, 1972. %1384--1394
%%CITATION = JPAGB,5,1384;%%
%-------------------------------------------------------------------%
 \bibitem[Nandi \& Sinha(1972b)]{Nandi72b}
          Nandi, B.C. and Sinha, M.S.,
%%%%%%%%%``The charge ratio of muons at sea level in the range 5--600
%%%%%%%%% GeV/c'',
          Nucl. Phys. B\,{\bf40}, 289, 1972. %289--297
%%CITATION = NUPHA,B40,289;%%
%-------------------------------------------------------------------%%
 \bibitem[Naumov(1984)]{Naumov84}
          Naumov, V.A.,
%%%%%%%%%``A method for calculating geomagnetic corrections to
%%%%%%%%% differential energy spectra of cosmic-ray nucleons in the
%%%%%%%%% atmosphere'',
          Investigations on Geomagnetism, Aeronomy and Solar Physics,
          {\bf69}, 82, 1984. %82--94
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Naumov(1993)]{Naumov93}
          Naumov, V.A.,
%%%%%%%%%``Cosmic-ray neutrinos at low and intermediate energies''
          in: Proc. of the Intern. Workshop on $\nu_{\mu}/\nu_e$
          Problem in Atmospheric Neutrinos, Gran Sasso, 1993, eds.
          Berezinsky, V.S. and Fiorentini, G. (LNGS, L'Aquila, 1993),
          p.\,25. %pp.\,25--70.
%%CITATION = NONE;%%
%--------------------------------------------------------------------%
 \bibitem[Naumov(1998)]{Naumov98a}
          Naumov, V.A.,
%%%%%%%%%``Physics of atmospheric neutrinos (introductory overview),''
          in: Proc. the Baikal School on Fundamental Physics
          ``Astrophysics and Microworld Physics'', Irkutsk, 1998,
          eds. Naumov, V.A., Parfenov, Yu.V, and Sinegovsky, S.I.,
          (Irkutsk State University, Irkutsk, 1998), p.\,67.
          %pp.\,67--85.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Naumov(2001)]{Naumov01b}
          Naumov, V.A.,
          {\em Atmospheric muons and neutrinos}, see URL of this
          Workshop, Program and Transparencies
          (Talks and Posters): \\ \Workshop.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Naumov \& Sinegovskaya(2000b)]{Naumov00b}
          Naumov, V.A. and Sinegovskaya, T.S.,
%%%%%%%%%``Simple method for solving transport equations describing
%%%%%%%%% the propagation of cosmic-ray nucleons in the atmosphere'',
          Yad. Fiz. {\bf63}, 2020, 2000 %2020--2028
          [Phys. Atom. Nucl. {\bf63}, 1927, 2000]. %1927--1935.
%%CITATION = PANUE,63,1927;%%
%-------------------------------------------------------------------%
 \bibitem[Naumov \& Sinegovskaya(2001)]{Naumov01a}
          Naumov, V.A. and Sinegovskaya, T.S.,
%%%%%%%%%``Atmospheric proton and neutron spectra at energies above
%%%%%%%%% 1 GeV'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf1}, p.\,4173 %pp.\,4173--4176
          .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Naumov {\em et al.}(1998)]{Naumov98b}
          Naumov, V.A., Sinegovskaya, T.S., and Sinegovsky, S.I.,
%%%%%%%%%``The $K_{\ell3}$ form factors and atmospheric neutrino
%%%%%%%%% flavor ratio at high  energies'',
          Nuovo Cim. A\,{\bf111}, 129, 1998. %129--148
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Naumov {\em et al.}(2000a)]{Naumov00a}
          Naumov, V.A., Sinegovskaya, T.S., and Sinegovsky, S.I.,
%%%%%%%%%``Potential of deep-underwater neutrino telescopes for
%%%%%%%%% recording muons from charm decay'',
          Yad. Fiz. {\bf63}, 2016, 2000 %2016--2019
          [Phys. Atom. Nucl. {\bf63}, 1923, 2000]. %1923--1926
%%CITATION = PANUE,63,1923;%%
%--------------------------------------------------------------------%
 \bibitem[Naumov {\em et al.}(1994)]{Naumov94}
          Naumov, V.A., Sinegovsky, S.I., and Bugaev, E.V.,
%%%%%%%%%``A new method for calculating the energy spectrum of cosmic
%%%%%%%%% ray muons under thick layers of matter'',
          Yad. Fiz. {\bf57}, 439, 1994 %439--451
          [Phys. Atom. Nucl. {\bf57}, 412, 1994]; %412--424
%%CITATION = PANUE,57,412;%%
          .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Owen \& Wilson(1949)]{Owen49}
          Owen, B.G. and Wilson, J.G.,
%%%%%%%%%``On the ratio of positive to negative particles in the
%%%%%%%%% vertical cosmic-ray beam at sea level'',
          Proc. Phys. Soc. (London) A\,{\bf62}, 601, 1949. %601--607
%%CITATION = PPSOA,A62,601;%%
%-------------------------------------------------------------------%
 \bibitem[Owen \& Wilson(1951)]{Owen51}
          Owen, B.G. and Wilson, J.G.,
%%%%%%%%%``Further measurements of the charge ratio of $\mu$-mesons
%%%%%%%%% at sea level'',
          Proc. Phys. Soc. (London) A\,{\bf64}, 417, 1951. %417--424
%%CITATION = PPSOA,A64,417;%%
%-------------------------------------------------------------------%
 \bibitem[Owen \& Wilson(1955)]{Owen55}
          Owen, B.G. and Wilson, J.G.,
%%%%%%%%%``Measurements of the momentum spectrum of $\mu$-mesons at
%%%%%%%%% sea level. I: The momentum range $5\times10^8$ eV/c --
%%%%%%%%% $2\times10^{10}$ eV/c'',
          Proc. Phys. Soc. (London) A\,{\bf68}, 409, 1955. %409--418
%%CITATION = PPSOA,A68,409;%%
%-------------------------------------------------------------------%
 \bibitem[Pak {\em et al.}(1961)]{Pak61}
          Pak, W. {\em et al.},
%%%%%%%%% Pak, W., Ozaki, S., Roe, B.P., and Greisen, K.,
%%%%%%%%%``Momentum spectrum of muons'',
          Phys. Rev. {\bf121}, 905, 1961. %905--907
%%CITATION = PHRVA,121,905;%%
%-------------------------------------------------------------------%
 \bibitem[Pasquali {\em et al.}(1998)]{Pasquali98}
          Pasquali, L., Reno, M.H., and Sarcevic, I.,
%%%%%%%%%``Secondary decays in atmospheric charm contributions to
%%%%%%%%% the flux of  muons and muon neutrinos,''
          Astropart. Phys. {\bf9}, 193, 1998. %193--202
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Pasquali {\em et al.}(1999)]{Pasquali99a}
          Pasquali, L., Reno, M.H., and Sarcevic, I.,
%%%%%%%%%``Lepton fluxes from atmospheric charm,''
          Phys. Rev. D\,{\bf59}, 034020, 1999.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Pasquali \& Reno(1999)]{Pasquali99b}
          Pasquali, L. and Reno, M.H.,
%%%%%%%%%``Tau neutrino fluxes from atmospheric charm,''
          Phys. Rev. D\,{\bf59}, 093003, 1999.
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Petersen(2001)]{Petersen01}
          Petersen, B. (for the L3 Collaboration),
%%%%%%%%%``First results from the L3+C experiment at CERN'',
          Nuovo Cim. C\,{\bf24}, 751, 2001. %751--754
%%CITATION = NUCIA,24C,751;%%
%-------------------------------------------------------------------%
 \bibitem[Pine {\em et al.}(1959)]{Pine59}
          Pine, J., Davidson, R.J., and Greisen, K.,
          Nuovo Cim. {\bf14}, 1181, 1959. %1181--?
%%CITATION = NUCIA,14,1181;%%
%-------------------------------------------------------------------%
 \bibitem[Plyaskin(2001)]{Plyaskin01}
          Plyaskin, V.,  % (for the AMS Collaboration),
%%%%%%%%%``Calculation of atmospheric neutrino flux'',
           Phys. Lett. B\,{\bf516}, 213, 2001. %213-235
%%%%%%%%% .
%%CITATION = ;%%
%--------------------------------------------------------------------%
 \bibitem[Poirier {\em et al.}(2001)]{Poirier01}
          Poirier, J., Roesler, S., and Fasso, A.,
%%%%%%%%%``Distributions of secondary muons at sea level from cosmic
%%%%%%%%% gamma rays below 10 TeV'',
           (also SLAC-PUB-8945).
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Rastin {\em et al.}(1965)]{Rastin65}
          Rastin, B.C. {\em et al.},
%%%%%%%%% Rastin, B.C., Baber, S.R., Bull, R.M., and Nash, W.F.,
          in: Proc. of the 9th Intern. Cosmic Ray Conf., London,
          1965,
%%%%%%%%% (The Inst. of Physics and The Physical Society, London,
%%%%%%%%% 1965),
           Vol.\,{\bf2}, p.\,981. %pp.\,981--?
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Rastin(1984a)]{Rastin84a}
          Rastin, B.C.,
%%%%%%%%%``An accurate measurement of the sea-level muon spectrum
%%%%%%%%% within the range 4 to 4000 GeV/c'',
          J. Phys G: Nucl. Phys. {\bf10}, 1609, 1984. %1609--1628
%%CITATION = JPHGB,10,1609;%%
%-------------------------------------------------------------------%
 \bibitem[Rastin(1984b)]{Rastin84b}
          Rastin, B.C.,
%%%%%%%%%``A study of the muon charge ratio at sea level within
%%%%%%%%% the momentum range 4 to 2000 GeV/c'',
          J. Phys G: Nucl. Phys. {\bf10}, 1629, 1984. %1629--1638
%%CITATION = JPHGB,10,1629;%%
%-------------------------------------------------------------------%
 \bibitem[Ryan {\em et al.}(1972)]{Ryan72}
          Ryan, M.J., Ormes, J.F., and Balasubrahmanyan, V.K.,
%%%%%%%%%``Cosmic-ray proton and helium spectra above 50 GeV'',
          Phys. Rev. Lett. {\bf28}, 985, 1972; %985--988
%%CITATION = PRLTA,28,985;%%
          {\em ibid.} {\bf28}, 1497, 1972 (Erratum).
%%CITATION = PRLTA,28,1497;%%
%-------------------------------------------------------------------%
 \bibitem[Sanuki {\em et al.}(2000)]{Sanuki00}
          Sanuki, T. {\em et al.},
%%%%%%%%% Sanuki, T., Motoki, M., Matsumoto, H., Seo, E.S.,
%%%%%%%%% {\em et al.} (BESS Collaboration),
%%%%%%%%%``Precise measurement of cosmic ray proton and helium
%%%%%%%%% spectra with the BESS spectrometer'',
          Astrophys. J. {\bf545}, 1135, 2000. %1135--1142
%%%%%%%%% .
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Sanuki(2001)]{Sanuki01a}
          Sanuki, T. (for the BESS collaboration),
%%%%%%%%%``BESS data on primary cosmic rays and muons'',
          Nucl. Phys. B (Proc. Suppl.) {\bf100}, 121, 2001. %121-123
%%CITATION = NUPHZ,100,121;%%
%-------------------------------------------------------------------%
 \bibitem[Sanuki {\em et al.}(2001)]{Sanuki01b}
          Sanuki, T. {\em et al.},
%%%%%%%%% Sanuki, T., Yamamoto, Y., Motoki, M., Abe, K. {\em et al.}
%%%%%%%%% (BESS Collaboration),
%%%%%%%%%``Atmospheric muons at various altitudes'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf1}, p.\,950. %pp.\,950--953
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Seo {\em et al.}(1991)]{Seo91}
          Seo, E.S. {\em et al.},
%%%%%%%%% Seo, E.S., Ormes, J.F., Streitmatter, R.E., Stochaj, S.J.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Measurement of cosmic-ray proton and helium spectra
%%%%%%%%% during the 1987 solar minimum'',
          Astrophys. J. {\bf378}, 763, 1991. %763-772
%%CITATION = ASJOA,378,763;%%
%-------------------------------------------------------------------%
 \bibitem[Serov \& Sychev(1973)]{Serov73}
          Serov, A.Ya. and Sychev, B.S.,
%%%%%%%%%``Phenomenological formulas for double-differential
%%%%%%%%% distributions of secondary hadrons, generated in
%%%%%%%%% nucleon-nucleus collisions at high energies'',
          Transactions of the Moscow Radio-Engineering Institute,
          {\bf11}, 173, 1973. %173--187
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Sinegovskaya and Sinegovsky(2001)]{Sinegovskaya01}
          Sinegovskaya, T.S. and Sinegovsky, S.I.,
%%%%%%%%%``Prompt muon contribution to the flux underwater'',
          Phys. Rev. D\,{\bf63}, 096004, 2001.
%%CITATION = ;%%
%-------------------------------------------------------------------%
 \bibitem[Singhal(1983)]{Singhal83}
          Singhal, K.P.,
          in: Proc. of the 18th Intern. Cosmic Ray Conf., Bangalore,
          1983,
%%%%%%%%% ed. by Durgaprasad, N {\em et al.} (Tata Institute of
%%%%%%%%% Fundamental Research, Colaba, Bombay, India, 1983).
          Vol.\,{\bf7}, p.\,27. %pp.\,27--?
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Smart \& Shea(1994)]{Smart94}
          Smart, D.F. and Shea, M.A.,
%%%%%%%%%``Geomagnetic cutoffs: a review for space dosimetry
%%%%%%%%% applications'',
          Adv. Space Res. {\bf14}, 787, 1994. %787--797
%%CITATION = ASRSD,14,787;%%
%--------------------------------------------------------------------%
 \bibitem[Stephens \& Golden(1987)]{Stephens87}
          Stephens, S.A. and Golden, R.L.,
%%%%%%%%%``Determination of the energy spectrum and charge
%%%%%%%%% composition of muons at sea-level'',
          in: Proc. of the 20th Intern. Cosmic Ray Conf., Moscow,
          1987,
%%%%%%%%% ed. by Kozyarivsky, V.A. {\em et al.} (Nauka, Moscow,
%%%%%%%%% 1987),
          Vol.\,{\bf6}, p.\,173. %pp.\,173--176
%%CITATION = NONE;%%
%--------------------------------------------------------------------%
 \bibitem[Sychev(1999)]{Sychev99}
         Sychev, B.S.,
         {\em Cross sections for interactions of high-energy hadrons
         with atomic nuclei},
         Moscow Radio-Engineering Inst. of the Russian Acad. of Sci.,
         Moscow, 1999.
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Thompson {\em et al.}(1977)]{Thompson77}
          Thompson, M.G. {\em et al.},
%%%%%%%%% Thompson, M.G., Thornley, R., Whalley, M.R., and
%%%%%%%%% Wolfendale, A.W.,
          in: Proc. of the 15th Intern. Cosmic Ray Conf., Plovdiv,
          1977,
%%%%%%%%% ed. by Christov, C.Ya. {\em et al.} (Institute for
%%%%%%%%% Nuclear Research, Bulgarian Academy of Sciences, Sofia,
%%%%%%%%% 1977),
          Vol.\,{\bf6}, p.\,21. %pp.\,21--25
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Tserkovnyak {\em et al.}(2001)]{Tserkovnyak01}
          Tserkovnyak, Y. {\em et al.},
%%%%%%%%% Tserkovnyak, Y., Komar, R.J., Nally, C.W., and
%%%%%%%%% Waltham, C.E.,
%%%%%%%%%``A three-dimensional calculation of atmospheric neutrino
%%%%%%%%% fluxes'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,1196. %pp.\,1196--1199
%%CITATION = NONE;%%
%-------------------------------------------------------------------%
 \bibitem[Tsuji {\em et al.}(1998)]{Tsuji98}
          Tsuji, S. {\em et al.},
%%%%%%%%% Tsuji, S., Katayama, T., Okei, K., Wada, T., {\em et al.},
%%%%%%%%%``Measurements of muons at sea level'',
          J. Phys. G: Nucl. Part. Phys. {\bf24}, 1805,
          1998. %1805--1822
%%CITATION = JPHGB,24,1805;%%
%--------------------------------------------------------------------%
 \bibitem[Tsuji {\em et al.}(2001)]{Tsuji01}
          Tsuji, S. {\em et al.},
%%%%%%%%% Tsuji, S., Katayama, T., Kawakami, N., Nakatsuka, T.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``Atmospheric muon measurements II: Zenith angular
%%%%%%%%% dependence'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,935; %pp.\,935--938
          Wada, T., private communication (July, 2001).
%%CITATION = NONE;%%
%--------------------------------------------------------------------%
 \bibitem[Vall {\em et al.}(1986)]{Vall86}
          Vall, A.N., Naumov, V.A., and Sinegovsky, S.I.,
%%%%%%%%%``The hadronic component of high-energy cosmic rays and the
%%%%%%%%% growth of the inelastic cross sections'',
          Yad. Fiz. {\bf44}, 1240, 1986 %1240--1250
          [Sov. J. Nucl. Phys. {\bf44}, 806, 1986]. %806--812
%%CITATION = SJNCA,44,806;%%
%-------------------------------------------------------------------%
 \bibitem[Volkova \& Zatsepin(2001)]{Volkova01}
          Volkova, L.V and Zatsepin, G.T.,
%%%%%%%%%``Fluxes of cosmic ray muons and atmospheric neutrinos at
%%%%%%%%% high energies,''
          Yad. Fiz. {\bf64}, 313, 2001 %313--322
          [Phys. Atom. Nucl. {\bf64}, 266, 2001]. %266--274
%%CITATION = PANUE,64,266;%%
%-------------------------------------------------------------------%
 \bibitem[Vulpescu {\em et al.}(1998)]{Vulpescu98}
          Vulpescu, B. {\em et al.},
%%%%%%%%% Vulpescu, B., Wentz, J., Brancus, I.M., Rebel, H.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``A compact detector for the measurement of the cosmic-ray
%%%%%%%%% muon charge ratio'',
          Nucl. Instr. Meth. A\,{\bf414}, 205, 1998. %205--217
%%CITATION = NUIMA,A414,205;%%
%-------------------------------------------------------------------%
 \bibitem[Vulpescu {\em et al.}(2001)]{Vulpescu01}
          Vulpescu, B. {\em et al.},
%%%%%%%%% Vulpescu, B., Wentz, J., Brancus, I.M., Badea, A.F.,
%%%%%%%%% {\em et al.},
%%%%%%%%%``The charge ratio of atmospheric muons below 1.0 GeV/c by
%%%%%%%%% measuring the lifetime of muonic atoms in aluminium,''
          J. Phys. G\,{\bf27}, 977, 2001. %977--991
%%CITATION = JPHGB,G27,977;%%
%-------------------------------------------------------------------%
 \bibitem[Wentz {\em et al.}(2001a)]{Wentz01a}
          Wentz, J. {\em et al.},
%%%%%%%%% Wentz, J., Bercuci, A., Heck, D., Mathes, H.J.,
%%%%%%%%% Oehlschl\"ager, J., Rebel, H., and Vulpescu, B.,
%%%%%%%%%``Simulation of atmospheric neutrino fluxes with CORSIKA'',
          in: Proc. of the 27th Intern. Cosmic Ray Conf., Hamburg,
          2001, Vol.\,{\bf3}, p.\,1167. %pp.\,1167--1170
%%CITATION = NONE;%%
%--------------------------------------------------------------------%
 \bibitem[Wentz {\em et al.}(2001b)]{Wentz01b}
          Wentz, J. {\em et al.},
%%%%%%%%% Wentz, J., Mathes, H.J., Rebel, H., Badea, A.F., Bercuci,
%%%%%%%%% A., Bozdog, H., Brancus, I.M., Petcu, M., and Vulpescu, B.,
%%%%%%%%%``The relevance of the muon charge ratio for the atmospheric
%%%%%%%%% neutrino anomaly,''
          J. Phys. G\,{\bf27}, 1699, 2001. %1699--1708
%%CITATION = JPHGB,G27,1699;%%
%--------------------------------------------------------------------%
 \bibitem[Wiebel-Sooth {\em et al.}(1998)]{Wiebel-Sooth98}
          Wiebel-Sooth, B., Biermann, P.L., and Meyer, H.,
%%%%%%%%%``Cosmic rays. VII. Individual element spectra: prediction
%%%%%%%%% and data'',
          Astron. Astrophys. {\bf330}, 389, 1998. %389--398
%%CITATION = AAEJA,330,389;%%
%-------------------------------------------------------------------%

\end{thebibliography}

\end{document}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% A D D I T I O N A L    R E F E R E N C E S                        %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%-------------------------------------------------------------------%
 \bibitem[De Pascale {\em et al.}(1993)]{DePascale93}
          De Pascale, M.P. {\em et al.},
%%%%%%%%% De Pascale, M.P., Morselli, A., Picozza, P., Golden, R.L.,
%%%%%%%%% {\em et al.},
          J. Geophys. Res. {\bf98}, No.\,A3, 3501, 1993. %3501--3507
%%CITATION = JGREA,98,3501;%%
%-------------------------------------------------------------------%


\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
\begin{figure*}[t]
\centering
\ifpdf
\includegraphics[width=8.4cm]{e1.pdf}~% for PDFLaTeX
\includegraphics[width=8.4cm]{e2.pdf} % for PDFLaTeX
\includegraphics[width=8.4cm]{m1.pdf}~% for PDFLaTeX
\includegraphics[width=8.4cm]{m2.pdf} % for PDFLaTeX
\else
\includegraphics[width=8.0cm]{e1.eps}~% for LaTeX
\includegraphics[width=8.0cm]{e2.eps} % for LaTeX
\includegraphics[width=8.0cm]{m1.eps}~% for LaTeX
\includegraphics[width=8.0cm]{m2.eps} % for LaTeX
\fi
\protect\caption{Energy spectra of $\nu_e$, $\overline{\nu}_e$,
                 $\nu_\mu$, and $\overline{\nu}_\mu$ for 11
                 zenith angles (low-energy range is for Kamioka).
\label{f:Nu_Conv}}
\end{figure*}
\clearpage % DO NOT DELETE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!









CITATION = ;
CITATION = ;
CITATION = ;
CITATION = ;
CITATION = ;
CITATION = NUCIA,C12,41;
CITATION = ;
CITATION = ;
CITATION = ;
CITATION = PANUE,64,266;
CITATION = JPHGB,G27,977;
CITATION = JPHGB,G27,1699;

