#ifndef MARS_MObservatory #define MARS_MObservatory #ifndef MARS_MParContainer #include "MParContainer.h" #endif class MObservatory : public MParContainer { public: enum LocationName_t { kMagic1, kWuerzburgCity }; private: LocationName_t fObservatoryKey; //! TString fObservatoryName; //! Name of the observatory Double_t fLongitude; //! [rad] Longitude of observatory (+ east) Double_t fLatitude; //! [rad] Latitude of observatory (+ north) Double_t fSinLatitude; //! Sin component for faster access Double_t fCosLatitude; //! Cos component for faster access Double_t fHeight; //! [m] height of observatory void Init(const char *name, const char *title); public: MObservatory(const char *name=NULL, const char *title=NULL); MObservatory(LocationName_t key, const char *name=NULL, const char *title=NULL); void SetLocation(LocationName_t name); void Print(Option_t *o=0) const; const TString &GetObservatoryName() const { return fObservatoryName; } Double_t GetLatitudeDeg() const { return fLatitude*kRad2Deg; } //[deg] Double_t GetLongitudeDeg() const { return fLongitude*kRad2Deg; } //[deg] Double_t GetLatitudeRad() const { return fLatitude; } //[rad] Double_t GetLongitudeRad() const { return fLongitude; } //[rad] Double_t GetPhi() const { return fLatitude; } //[rad] Double_t GetElong() const { return fLongitude; } //[rad] Double_t GetSinPhi() const { return fSinLatitude; } Double_t GetCosPhi() const { return fCosLatitude; } Double_t GetHeight() const { return fHeight; } void RotationAngle(Double_t theta, Double_t phi, Double_t &sin, Double_t &cos) const; Double_t RotationAngle(Double_t theta, Double_t phi) const; LocationName_t GetObservatoryKey() const { return fObservatoryKey; } ClassDef(MObservatory, 0) // class storing observatory locations }; #endif