Index: trunk/MagicSoft/Mars/mpointing/MHPointing.cc =================================================================== --- trunk/MagicSoft/Mars/mpointing/MHPointing.cc (revision 7198) +++ trunk/MagicSoft/Mars/mpointing/MHPointing.cc (revision 7198) @@ -0,0 +1,292 @@ +/* ======================================================================== *\ +! +! * +! * This file is part of MARS, the MAGIC Analysis and Reconstruction +! * Software. It is distributed to you in the hope that it can be a useful +! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. +! * It is distributed WITHOUT ANY WARRANTY. +! * +! * Permission to use, copy, modify and distribute this software and its +! * documentation for any purpose is hereby granted without fee, +! * provided that the above copyright notice appear in all copies and +! * that both that copyright notice and this permission notice appear +! * in supporting documentation. It is provided "as is" without express +! * or implied warranty. +! * +! +! +! Author(s): Thomas Bretz, 05/2005 +! +! Copyright: MAGIC Software Development, 2000-2005 +! +! +\* ======================================================================== */ + +///////////////////////////////////////////////////////////////////////////// +// +// MHPointing +// +// Display drive information +// +//////////////////////////////////////////////////////////////////////////// +#include "MHPointing.h" + +#include +#include +#include + +#include "MLog.h" +#include "MLogManip.h" + +#include "MParList.h" + +#include "MTime.h" + +#include "MReportDrive.h" +#include "MReportStarguider.h" + +ClassImp(MHPointing); + +using namespace std; + +void MHPointing::ResetGraph(TGraph &g) const +{ + g.Set(1); + g.SetPoint(0, 0, 0); // Dummy Point + g.SetEditable(); // Used as flag: First point? yes/no +} + +void MHPointing::InitGraph(TGraph &g) const +{ + ResetGraph(g); + g.SetMarkerStyle(kFullDotMedium); +} + +void MHPointing::AddPoint(TGraph &g, Double_t x, Double_t y) const +{ + if (g.IsEditable()) + { + g.RemovePoint(0); + g.SetEditable(kFALSE); + } + + g.SetPoint(g.GetN(), x, y); +} + +// -------------------------------------------------------------------------- +// +// Setup histograms +// +MHPointing::MHPointing(const char *name, const char *title) +: /*fTime(0),*/ fReportCosy(0), fReportSG(0) +{ + fName = name ? name : "MHPointing"; + fTitle = title ? title : "Graphs for rate data"; + + // Init Graphs + fDevTimeSG.SetNameTitle("DevSG", "Drive (black) and starguider (blue) deviation"); + fDevTimeCosy.SetNameTitle("DevCosy", "Cosy deviation"); + fBrightness.SetNameTitle("Brigtness", "Arbitrary Sky Brightness"); + fNumStars.SetNameTitle("NumStars", "Number of stars identified by starguider"); + fDevZd.SetNameTitle("DevZd", "Starguider Deviation Zd (blue), Az(black)"); + fDevAz.SetNameTitle("DevAz", "Starguider Deviation Az"); + + InitGraph(fDevTimeSG); + InitGraph(fDevTimeCosy); + InitGraph(fBrightness); + InitGraph(fNumStars); + InitGraph(fDevZd); + InitGraph(fDevAz); + + fDevTimeSG.SetMinimum(0); + fDevTimeCosy.SetMinimum(0); + fBrightness.SetMinimum(0); + fNumStars.SetMinimum(0); + + fDevTimeSG.SetMarkerColor(kBlue); + fDevZd.SetMarkerColor(kBlue); +} + +// -------------------------------------------------------------------------- +// +// Setup the Binning for the histograms automatically if the correct +// instances of MBinning +// +Bool_t MHPointing::SetupFill(const MParList *plist) +{ + //fTime = (MTime*)plist->FindObject("MTime"); + //if (!fTime) + //{ + // *fLog << warn << "MTime not found... abort." << endl; + // return kFALSE; + //} + fReportSG = (MReportStarguider*)plist->FindObject("MReportStarguider"); + if (!fReportSG) + { + *fLog << warn << "MReportStarguider not found... abort." << endl; + return kFALSE; + } + fReportCosy = (MReportDrive*)plist->FindObject("MReportDrive"); + if (!fReportCosy) + { + *fLog << warn << "MReportDrive not found... abort." << endl; + return kFALSE; + } + + // Reset Graphs + ResetGraph(fDevTimeSG); + ResetGraph(fDevTimeCosy); + + return kTRUE; +} + +// -------------------------------------------------------------------------- +// +// Fill the histograms with data from a MMuonCalibPar and +// MMuonSearchPar container. +// +Bool_t MHPointing::Fill(const MParContainer *par, const Stat_t w) +{ + const MTime *t = dynamic_cast(par); + if (!t) + { + *fLog << err <<"MHPointing::Fill - ERROR: MTime not given as argument... abort." << endl; + return kERROR; + } + + const Double_t tm = t->GetAxisTime(); + + if (t->GetName()==(TString)"MTimeStarguider") + { + AddPoint(fDevTimeSG, tm, fReportSG->GetDevAbs()); + AddPoint(fBrightness, tm, fReportSG->GetSkyBrightness()); + AddPoint(fNumStars, tm, fReportSG->GetNumIdentifiedStars()); + AddPoint(fDevZd, tm, fReportSG->GetDevZd()); + AddPoint(fDevAz, tm, fReportSG->GetDevAz()); + return kTRUE; + } + + if (t->GetName()==(TString)"MTimeDrive") + { + AddPoint(fDevTimeCosy, tm, fReportCosy->GetAbsError()*60); + return kTRUE; + } + + return kTRUE; +} + +void MHPointing::DrawGraph(TGraph &g, const char *y) const +{ + TH1 *h = g.GetHistogram(); + if (h) + { + TAxis *axe = h->GetXaxis(); + axe->SetLabelSize(0.033); + axe->SetTimeFormat("%H:%M:%S %F1995-01-01 00:00:00 GMT"); + axe->SetTimeDisplay(1); + axe->SetTitle("Time"); + if (y) + h->SetYTitle(y); + } +} +/* +// -------------------------------------------------------------------------- +// +// Update position of an axis on the right side of the histogram +// +void MHPointing::UpdateRightAxis(TGraph &g) const +{ + TH1 &h = *g.GetHistogram(); + + const Double_t max = h.GetMaximum(); + if (max==0) + return; + + TGaxis *axis = (TGaxis*)gPad->FindObject("RightAxis"); + if (!axis) + return; + + axis->SetX1(g.GetXaxis()->GetXmax()); + axis->SetX2(g.GetXaxis()->GetXmax()); + axis->SetY1(gPad->GetUymin()); + axis->SetY2(gPad->GetUymax()); + axis->SetWmax(max); +} + +// -------------------------------------------------------------------------- +// +// Draw an axis on the right side of the histogram +// + +void MHPointing::DrawRightAxis(const char *title) const +{ + TGaxis *axis = new TGaxis(gPad->GetUxmax(), gPad->GetUymin(), + gPad->GetUxmax(), gPad->GetUymax(), + 0, 1, 510, "+L"); + axis->SetName("RightAxis"); + axis->SetTitle(title); + axis->SetTitleOffset(0.9); + axis->SetTextColor(kRed); + axis->SetBit(kCanDelete); + axis->Draw(); +} +*/ +// -------------------------------------------------------------------------- +// +// This displays the TGraph like you expect it to be (eg. time on the axis) +// It should also make sure, that the displayed time always is UTC and +// not local time... +// +void MHPointing::Draw(Option_t *opt) +{ + TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this); + pad->SetBorderMode(0); + + AppendPad(); + + pad->Divide(2,2); + + pad->cd(1); + gPad->SetBorderMode(0); + gPad->SetGridx(); + gPad->SetGridy(); + fDevTimeSG.Draw("AP"); + fDevTimeCosy.Draw("P"); + + pad->cd(2); + gPad->SetBorderMode(0); + gPad->SetGridx(); + gPad->SetGridy(); + fBrightness.Draw("AP"); + + pad->cd(3); + gPad->SetBorderMode(0); + gPad->SetGridx(); + gPad->SetGridy(); + fDevZd.Draw("AP"); + fDevAz.Draw("P"); + + pad->cd(4); + gPad->SetBorderMode(0); + gPad->SetGridx(); + gPad->SetGridy(); + fNumStars.Draw("AP"); +} + +void MHPointing::Paint(Option_t *o) +{ + DrawGraph(fDevTimeSG, "\\Delta [arcmin]"); + DrawGraph(fDevTimeCosy, "\\Delta [arcmin]"); + DrawGraph(fBrightness, "Brightness [au]"); + DrawGraph(fNumStars, "N"); + DrawGraph(fDevZd, "\\Delta [arcmin]"); + DrawGraph(fDevAz, "\\Delta [arcmin]"); + /* + gPad->cd(1); + + if (gPad) + { + fHumidity.GetHistogram()->GetYaxis()->SetTitleColor(kBlue); + UpdateRightAxis(fTemperature); + }*/ +} Index: trunk/MagicSoft/Mars/mpointing/MHPointing.h =================================================================== --- trunk/MagicSoft/Mars/mpointing/MHPointing.h (revision 7198) +++ trunk/MagicSoft/Mars/mpointing/MHPointing.h (revision 7198) @@ -0,0 +1,52 @@ +#ifndef MARS_MHPointing +#define MARS_MHPointing + +#ifndef MARS_MH +#include "MH.h" +#endif +#ifndef ROOT_TGraph +#include +#endif + +class MTime; +class MReportDrive; +class MReportStarguider; + +class MHPointing : public MH +{ +private: + MReportDrive *fReportCosy; //! + MReportStarguider *fReportSG; //! + + TGraph fDevTimeSG; // Starguider deviation versus time + TGraph fDevTimeCosy; // Drive deviation versus time + + TGraph fBrightness; // Arbitrary sky brightness + + TGraph fNumStars; // Number of stars identified by starguider + + TGraph fDevZd; // Starguider deviation Zd + TGraph fDevAz; // Starguider deviation Az + + void ResetGraph(TGraph &g) const; + void InitGraph(TGraph &g) const; + void AddPoint(TGraph &g, Double_t x, Double_t y) const; + void DrawGraph(TGraph &g, const char *y=0) const; + //void UpdateRightAxis(TGraph &g) const; + //void DrawRightAxis(const char *title) const; + +public: + MHPointing(const char *name=NULL, const char *title=NULL); + + Bool_t SetupFill(const MParList *plist); + Bool_t Fill(const MParContainer *par, const Stat_t w=1); + + void Draw(Option_t *opt=""); + void Paint(Option_t *opt=""); + + ClassDef(MHPointing, 1) // Histogram to display tracking/pointing information +}; + +#endif + + Index: trunk/MagicSoft/Mars/mpointing/Makefile =================================================================== --- trunk/MagicSoft/Mars/mpointing/Makefile (revision 7196) +++ trunk/MagicSoft/Mars/mpointing/Makefile (revision 7198) @@ -26,4 +26,5 @@ MPointingPosCalc.cc \ MPointingPosInterpolate.cc \ + MHPointing.cc \ MSrcPosCam.cc \ MSrcPosCalc.cc \