/**
   \file
   Support utilities for the AIRES shower file. Shamelessly copied from 
   the CORSIKA implementation

   \author Troy Porter
   \version $Id: AiresUtilities.h 14717 2009-09-17 20:24:36Z lukas $
   \date 06 May 2004
   \ingroup aires
*/

#ifndef _io_AiresUtilities_h__
#define _io_AiresUtilities_h_

static const char CVSId_io_AiresUtilities[] = 

"$Id: AiresUtilities.h 14717 2009-09-17 20:24:36Z lukas $";

#include <utl/GeometryUtilities.h>
#include <utl/AugerUnits.h>

#include <iostream>

namespace io {

  namespace Aires {

    // This value just initialises the variable. The declination proper is 
    // set via the value taken from the Aires shower file header record. See
    // AiresShowerFile::Open() for further details.
    static double kMagneticFieldDeclination = 0.0; 
    
    /// Convert AIRES particle code to PDG
    int AiresToPDG(int theAiresCode);

    /// Returns the azimuth rotated from AIRES's system to Auger standard
    inline double AiresAzimuthToAuger(const double airesAzimuth) { 
      
      using utl::deg;
      using utl::GeometryUtilities::NormalizeAngle;

      return NormalizeAngle(airesAzimuth + 90*deg - io::Aires::kMagneticFieldDeclination); 

    }

  } // namespace AIRES

} // io


#endif 

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: