MagneticFieldModel.h

Go to the documentation of this file.

#ifndef _fwk_MagneticFieldModel_h_
#define _fwk_MagneticFieldModel_h_

//#include <utl/AugerException.h>
#include <utl/Vector.h>
#include <utl/Point.h>
#include <utl/AugerUnits.h>
#include <utl/UTCDateTime.h>

#include <fwk/CoordinateSystemRegistry.h>
#include <fwk/LocalCoordinateSystem.h>

#include "ParametricGeoMagneticField.h"
#include "ParametricGeoMagChange.h"

namespace fwk {

class MagneticFieldModel {

 public:
    static MagneticFieldModel& instance();

    static utl::Vector GetMagneticFieldVector(const utl::Point &position,
                                              const utl::TimeStamp &time);

    static double GetDeclination(const utl::Point &position, 
                                 const utl::TimeStamp &time);

    static double GetInclination(const utl::Point &position, 
                                 const utl::TimeStamp &time);

    static double GetStrength(const utl::Point &position,
                              const utl::TimeStamp &time) 
    {return GetMagneticFieldVector(position, time).GetMag();}


    ParametricGeoMagneticField ModelForYear(double year);
    ParametricGeoMagneticField ModelForTime(const utl::UTCDate &date);
    ParametricGeoMagneticField ModelForTime(const utl::TimeStamp &time)
    { 
      const utl::UTCDateTime date(time);
      return ModelForTime(date);
    }

 private:
   MagneticFieldModel();
  void Init();
  
 protected:
  std::vector <ParametricGeoMagneticField> fModels;
  ParametricGeoMagChange *fChangeModel;
  
  double decimalYear(int year, int month, int day);
};

} // fwk

#endif