#ifndef MARS_MAstroCatalog #define MARS_MAstroCatalog #ifndef ROOT_TVector3 #include #endif #ifndef ROOT_TList #include #endif #ifndef MARS_MGMap #include #endif class MTime; class MObservatory; class TArrayI; class TGToolTip; class MVector3 : public TVector3 { private: enum VectorType_t { kIsInvalid, kIsRaDec, kIsZdAz, kIsAltAz, kIsArbitrary }; VectorType_t fType; TString fName; public: MVector3() { fType=kIsInvalid; } MVector3(const TVector3 &v3) : TVector3(v3) { fType=kIsArbitrary; } Double_t Magnitude() const { return -2.5*TMath::Log10(Mag()); } void SetRaDec(Double_t ra, Double_t dec, Double_t mag=0) { fType = kIsRaDec; SetMagThetaPhi(pow(10, -mag/2.5), TMath::Pi()/2-dec, ra); } void SetName(const TString &str) { fName = str.Strip(TString::kBoth); } void SetZdAz(Double_t zd, Double_t az, Double_t mag=0) { fType = kIsZdAz; SetMagThetaPhi(pow(10, -mag/2.5), zd, az); } void SetAltAz(Double_t alt, Double_t az, Double_t mag=0) { fType = kIsAltAz; SetMagThetaPhi(pow(10, -mag/2.5), TMath::Pi()/2-alt, az); } const char *GetName() const { return fName; } void WriteBinary(ostream &out) const; void ReadBinary(istream &in); ClassDef(MVector3, 1) // A specialized TVector3 storing a star-name }; class MAstroCatalog : public TObject { private: Double_t fLimMag; // [1] Limiting Magnitude Double_t fRadiusFOV; // [deg] Radius of Field of View TGToolTip *fToolTip; //! The tooltip currently displayed void ShowToolTip(Int_t px, Int_t py, const char *txt); TString FindToken(TString &line, Char_t tok=','); Int_t atoi(const TSubString &sub); Float_t atof(const TSubString &sub); Int_t atoi(const TString &s); Float_t atof(const TString &s); //#if ROOT_VERSION_CODE < ROOT_VERSION(4,00,03) Bool_t fPlainScreen; //! Just a dummy!!!! ([Set,Is]Freezed) Bool_t fNoToolTips; //! Just a dummy!!!! ([Set,Is]Freezed) //#endif virtual Int_t ConvertToPad(const TVector3 &w, TVector2 &v) const; virtual void AddPrimitives(TString o); virtual void SetRangePad(Option_t *o); Int_t Convert(const TRotation &rot, TVector2 &v) const; void Draw(const TVector2 &v0, const TRotation &rot, TArrayI &dx, TArrayI &dy, Int_t stepx, Int_t stepy, Int_t type); void DrawPrimitives(Option_t *o); Bool_t DrawLine(const TVector2 &v0, Int_t dx, Int_t dy, const TRotation &rot, Int_t type); void DrawGrid(const TVector3 &v0, const TRotation &rot, Int_t type); void Paint(Option_t *o=""); Int_t DistancetoPrimitive(Int_t px, Int_t py); protected: enum { kHasChanged = BIT(15), // Display has changed kGuiActive = BIT(16), // GUI is interactive kPlainScreen = BIT(17), // View is a plain screen view kMirrorX = BIT(18), // Mirror display in X kMirrorY = BIT(19), // Mirror display in Y kNoToolTips = BIT(20), // suppress tooltips kDrawingImage = BIT(21) // just drawing into a bitmap }; MGMap fMapG; //! A map with all gui primitives and tooltips TList fList; // List of stars loaded MVector3 fRaDec; // pointing position MObservatory *fObservatory; // Possible obervatory location MTime *fTime; // Possible observation time TRotation AlignCoordinates(const TVector3 &v) const; virtual TString GetPadTitle() const; TRotation GetGrid(Bool_t local); void DrawStar(Double_t x, Double_t y, const TVector3 &v, Int_t col, const char *txt=0, Bool_t resize=kFALSE); void Update(Bool_t upd=kFALSE); void ExecuteEventKbd(Int_t keycode, Int_t keysym); void ExecuteEvent(Int_t event, Int_t mp1, Int_t mp2); char *GetObjectInfo(Int_t px, Int_t py) const; public: MAstroCatalog(); ~MAstroCatalog(); void SetTime(const MTime &time); void SetObservatory(const MObservatory &obs); void SetLimMag(Double_t mag) { fLimMag=mag; Update(); } // *MENU* *ARGS={mag=>fLimMag} void SetRadiusFOV(Double_t deg) { //const Double_t max = TestBit(kPlainScreen) ? 90 : 55; const Double_t max = TestBit(kPlainScreen) ? 180 : 90; if (deg>max) deg=max; if (deg<1) deg=1; fRadiusFOV=deg; Update(); } // *MENU* *ARGS={deg=>fRadiusFOV} void SetRadiusFOV(Double_t pixsize, Double_t w, Double_t h) { // pixsize [arcsec/pixel] // w [pixel] // h [pixel] const Double_t scale = TMath::Hypot(w, h)/2; SetRadiusFOV(scale*pixsize/3600); } // Set Radius of FOV using the pixsize [arcsec/pix], width and height [pixel] of image void SetRaDec(Double_t ra, Double_t dec) { fRaDec.SetRaDec(ra, dec, 1); Update(); } void SetRaDec(const TVector3 &v) { fRaDec=v; Update(); } void SetGuiActive(Bool_t b=kTRUE) { b ? SetBit(kGuiActive) : ResetBit(kGuiActive); } void SetPlainScreen(Bool_t b=kTRUE) { b ? SetBit(kPlainScreen) : ResetBit(kPlainScreen); Update(); } // *TOGGLE* *GETTER=IsPlainScreen Bool_t IsPlainScreen() const { return TestBit(kPlainScreen); } void SetNoToolTips(Bool_t b=kTRUE) { b ? SetBit(kNoToolTips) : ResetBit(kNoToolTips); } // *TOGGLE* *GETTER=HasNoToolTips Bool_t HasNoToolTips() const { return TestBit(kNoToolTips); } Double_t GetLimMag() const { return fLimMag; } // Get Limiting Magnitude Double_t GetRadiusFOV() const { return fRadiusFOV; } // Get maximum radius of Field Of View void Delete(Option_t *o="") { fList.Delete(); fMapG.Delete(); } // Delete list of stars Int_t ReadXephem(TString catalog = "/usr/X11R6/lib/xephem/catalogs/YBS.edb"); Int_t ReadNGC2000(TString catalog = "ngc2000.dat"); Int_t ReadBSC(TString catalog = "bsc5.dat"); Int_t ReadHeasarcPPM(TString catalog = "heasarc_ppm.tdat", TString fout=""); Int_t ReadCompressed(TString catalog); Bool_t AddObject(Float_t ra, Float_t dec, Float_t mag, TString name=""); void Print(Option_t *o="") const { fList.Print(); } // Print all stars TList *GetList() { return &fList; } // Return list of stars //void PaintImg(Int_t id=0, Option_t *o=""); void PaintImg(unsigned char *buf, int w=768, int h=576, Option_t *o=0); void Draw(Option_t *o=""); void SetDrawOption(Option_t *option="") { TObject::SetDrawOption(option); Update(kTRUE); } //*MENU* virtual void EventInfo(Int_t event, Int_t px, Int_t py, TObject *selected=0); ClassDef(MAstroCatalog, 1) // Display class for star catalogs }; #endif