On the parametrization of atmospheric muon angular flux underwater
  The analytical expression for angular integral flux of atmospheric muons in matter with the explicit relation of its parameters with those of the sea level spectrum is obtained. The fitting formula for the sea level muon spectrum at different zenith angles for spherical atmosphere is proposed. The concrete calculations for pure water are presented. Fluctuations of muon energy losses are taken into account by means of parametrized correction factor calculated using survival probabilities resulted from Monte Carlo simulations. Parametrizations of all continuous energy losses are obtained with using the most recent expressions for muon interaction cross-sections. The corresponding parametrization errors and field of method application are comprehensively discussed. The proposed formulae could be useful primarily for experimentalists processing data of arrays located deep under water or under ice. 
S.I.Klimushin, E.V.Bugaev, I.A.Sokalski (Institute for Nuclear   Research, Russian Academy of Sciences, Moscow)