Index: trunk/MagicSoft/Mars/Changelog =================================================================== --- trunk/MagicSoft/Mars/Changelog (revision 8993) +++ trunk/MagicSoft/Mars/Changelog (revision 8994) @@ -29,4 +29,9 @@ - check for the combination of telescope-, run- and file-number to be unique. + + * datacenter/macros/plotrundb.C: + - use MPlot from plotdb + - removed obsolete DataCheck plots + - plot only range of dates which makes sense Index: trunk/MagicSoft/Mars/datacenter/macros/plotrundb.C =================================================================== --- trunk/MagicSoft/Mars/datacenter/macros/plotrundb.C (revision 8993) +++ trunk/MagicSoft/Mars/datacenter/macros/plotrundb.C (revision 8994) @@ -1,4 +1,4 @@ /* ======================================================================== *\ -! $Name: not supported by cvs2svn $:$Id: plotrundb.C,v 1.4 2007-10-07 21:12:15 tbretz Exp $ +! $Name: not supported by cvs2svn $:$Id: plotrundb.C,v 1.5 2008-07-03 16:29:48 tbretz Exp $ ! -------------------------------------------------------------------------- ! @@ -21,5 +21,5 @@ ! Author(s): Daniela Dorner, 05/2005 ! -! Copyright: MAGIC Software Development, 2000-2006 +! Copyright: MAGIC Software Development, 2000-2008 ! ! @@ -59,460 +59,5 @@ // ///////////////////////////////////////////////////////////////////////////// -#include -#include - -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include - -#include "MTime.h" -#include "MAstro.h" -#include "MDataSet.h" -#include "MSQLMagic.h" -#include "MStatusDisplay.h" - -class MPlot : public MParContainer -{ -public: - // Possible constants to group-by (average) over a certain period - enum GroupBy_t - { - kNone, - kGroupByPrimary, - kGroupByHour, - kGroupByNight, - kGroupByWeek, - kGroupByMonth, - kGroupBySeason, - kGroupByYear - }; - -private: - MSQLMagic &fServer; // Reference to the sql-server class - - MDataSet *fDataSet; // A possible dtaset to highlite single points - - TString fPrimaryDate; // The name of the data we plot - TString fPrimaryNumber; // The corresponding name for the key number - TString fSecondary; // The value versus which the second plot is made - - TString fRequestFrom; // Start of a requested date range - TString fRequestTo; // End of a requested date range - Int_t fRequestPeriod; // A possible requested period - - Float_t fPlotMin; - Float_t fPlotMax; - - Float_t fHistMin; - Float_t fHistMax; - - TString fDescription; // The description (title) of the plot - TString fNameTab; // The name of the tab in the display - - TString fCondition; // An additional condition added to the query - GroupBy_t fGroupBy; // A possible Group-By flag - - // -------------------------------------------------------------------------- - // - // Function to plot the result of the query - // - void PlotTable(TSQLResult &res, TString name, Float_t fmin, Float_t fmax, Float_t resolution) - { - // Enable all otions in the statistics box - gStyle->SetOptStat(111111); - - // Create TGraph objects - TGraph > = res.GetFieldCount()>4 ? *new TGraphErrors : *new TGraph; - gt.SetNameTitle(name, Form("%s vs Time", name.Data())); - gt.SetMarkerStyle(kFullDotMedium); - - TGraph gz; - gz.SetNameTitle(name, Form("%s vs ", name.Data())); - gz.SetMarkerStyle(kFullDotMedium); - - TGraph gt0, gt1; - gt0.SetMarkerColor(kRed); - gt1.SetMarkerColor(kBlue); - gt0.SetMarkerStyle(kFullDotLarge); - gt1.SetMarkerStyle(kFullDotLarge); - - TGraph gz0, gz1; - gz0.SetMarkerColor(kRed); - gz1.SetMarkerColor(kBlue); - gz0.SetMarkerStyle(kFullDotLarge); - gz1.SetMarkerStyle(kFullDotLarge); - - Int_t first = -1; - Int_t last = -1; - - // Loop over the data - TSQLRow *row=0; - while ((row=res.Next())) - { - // Get all fields of this row - const char *date = (*row)[0]; - const char *zd = (*row)[1]; - const char *val = (*row)[2]; - const char *snum = res.GetFieldCount()>3 ? (*row)[3] : 0; - const char *verr = res.GetFieldCount()>4 ? (*row)[5] : 0; - - delete row; - - if (!date || !val || !zd) - continue; - - // check if date is valid - MTime t(date); - if (!t.SetSqlDateTime(date)) - continue; - - // check if it belongs to the requested MAGIC period - if (fRequestPeriod>0 && MAstro::GetMagicPeriod(t.GetMjd())!=fRequestPeriod) - continue; - - // Get axis range - if (first<0) - first = TMath::Nint(TMath::Floor(t.GetMjd())); - last = TMath::Nint(TMath::Ceil(t.GetMjd())); - - // Convert a possible key number into a integer - UInt_t seq = snum ? atoi(snum) : 0; - - // convert primary and secondary value into floats - Float_t value = atof(val); - Float_t zenith = atof(zd); - - // If a datset is given add the point to the special TGraphs - // used for highliting these dates - if (fDataSet) - { - if (fDataSet->HasOnSequence(seq)) - { - gt1.SetPoint(gt1.GetN(), t.GetAxisTime(), value); - gz1.SetPoint(gz1.GetN(), zenith, value); - } - - if (fDataSet->HasOffSequence(seq)) - { - gt0.SetPoint(gt0.GetN(), t.GetAxisTime(), value); - gz0.SetPoint(gz0.GetN(), zenith, value); - } - } - - // Add Data to TGraph - gt.SetPoint(gt.GetN(), t.GetAxisTime(), value); - gz.SetPoint(gz.GetN(), zenith, value); - - // Set error-bar, if one - if (verr) - static_cast(gt).SetPointError(gt.GetN()-1, 0, atof(verr)); - } - - // If this is done earlier the plots remain empty since root 5.12/00 - if (fmax>fmin) - { - gt.SetMinimum(fmin); - gt.SetMaximum(fmax); - gz.SetMinimum(fmin); - gz.SetMaximum(fmax); - } - - gROOT->SetSelectedPad(0); - - // Create a TCanvas or open a new tab - TString title = fNameTab.IsNull() ? name(name.First('.')+2, name.Length()) : fNameTab; - TCanvas &c = fDisplay ? fDisplay->AddTab(title) : *new TCanvas; - // Set fillcolor, remove border and divide pad - c.SetFillColor(kWhite); - c.SetBorderMode(0); - c.Divide(1,2); - - // Output mean and rms to console - cerr << setprecision(4) << setw(10) << title << ": "; - if (gt.GetN()==0) - { - cerr << " " << endl; - return; - } - cerr << setw(8) << gt.GetMean(2) << "+-" << setw(8) << gt.GetRMS(2) << " "; - if (gt0.GetN()>0 || gt1.GetN()>0) - { - cerr << setw(8) << gt1.GetMean(2) << "+-" << setw(8) << gt1.GetRMS(2) << " "; - cerr << setw(8) << gt0.GetMean(2) << "+-" << setw(8) << gt0.GetRMS(2); - } - cerr << endl; - - TVirtualPad *pad = gPad; - - // draw contants of pad 2 (counting starts at 0) - pad->cd(2); - gPad->SetBorderMode(0); - gPad->SetFrameBorderMode(0); - gPad->SetGridy(); - - gPad->SetLeftMargin(0.06); - gPad->SetRightMargin(0.06); - gPad->SetBottomMargin(0.08); - - // format axis - TH1 *h = gt.GetHistogram(); - - h->SetXTitle("Time"); - h->SetYTitle(name); - h->GetXaxis()->SetTimeDisplay(1); - h->GetYaxis()->SetTitleOffset(0.8); - h->GetXaxis()->SetTitleOffset(1.0); - h->GetXaxis()->SetLabelOffset(0.01); - - // draw TGraph - gt.DrawClone("AP"); - if (gt0.GetN()>0) - gt0.DrawClone("P"); - if (gt1.GetN()>0) - gt1.DrawClone("P"); - - // Add lines and text showing the MAGIC periods - TLine l; - TText t; - Int_t num=0; - l.SetLineStyle(kDotted); - l.SetLineColor(kBlue); - t.SetTextColor(kBlue); - l.SetLineWidth(1); - t.SetTextSize(h->GetXaxis()->GetLabelSize()); - t.SetTextAlign(21); - Int_t p0 = MAstro::GetMagicPeriod(first); - for (Int_t p = first; pGetMinimum(), MTime(p).GetAxisTime(), h->GetMaximum()); - t.DrawText(MTime(p+15).GetAxisTime(), h->GetMaximum(), Form("%d", p1)); - num++; - } - p0 = p1; - } - if (num<4) - gPad->SetGridx(); - - const Double_t min = fHistMin>fHistMax ? h->GetMinimum()-resolution/2 : fHistMin; - const Double_t max = fHistMin>fHistMax ? h->GetMaximum()+resolution/2 : fHistMax; - - // Use this to save the pad with the time development to a file - //gPad->SaveAs(Form("plotdb-%s.eps", title.Data())); - - // Go back to first (upper) pad, format it and divide it again - pad->cd(1); - gPad->SetBorderMode(0); - gPad->SetFrameBorderMode(0); - gPad->Divide(2,1); - - TVirtualPad *pad2 = gPad; - - // format left pad - pad2->cd(1); - gPad->SetBorderMode(0); - gPad->SetFrameBorderMode(0); - gPad->SetGridx(); - gPad->SetGridy(); - - // Create histogram - const Int_t n = resolution>0 ? TMath::Nint((max-min)/resolution) : 50; - - TH1F hist("Hist", Form("Distribution of %s", fDescription.IsNull() ? name.Data() : fDescription.Data()), n, min, max); - hist.SetDirectory(0); - - // Fill data into histogra, - for (int i=0; icd(2); - gPad->SetBorderMode(0); - gPad->SetFrameBorderMode(0); - gPad->SetGridy(); - - // format graph - TH1 *h2 = gz.GetHistogram(); - - h2->SetXTitle("Zd"); - h2->SetYTitle(name); - - // draw graph - gz.DrawClone("AP"); - - if (gz0.GetN()>0) - gz0.DrawClone("P"); - if (gz1.GetN()>0) - gz1.DrawClone("P"); - } - -public: - MPlot(MSQLMagic &server) : fServer(server), fDataSet(NULL), - fRequestPeriod(-1), fPlotMin(0), fPlotMax(-1), fHistMin(0), fHistMax(-1), fGroupBy(kNone) - { - } - ~MPlot() - { - if (fDataSet) - delete fDataSet; - } - void SetDataSet(const TString filename) - { - if (fDataSet) - { - delete fDataSet; - fDataSet = NULL; - } - if (!filename.IsNull()) - fDataSet = new MDataSet(filename); - } - void SetPlotRange(Float_t min, Float_t max, Int_t n=5) { fPlotMin = min; fPlotMax = max; } - void SetHistRange(Float_t min, Float_t max) { fHistMin = min; fHistMax = max; } - void SetRequestRange(const char *from="", const char *to="") { fRequestFrom = from; fRequestTo = to; } - void SetRequestPeriod(Int_t n=-1) { fRequestPeriod = n; } - void SetCondition(const char *cond="") { fCondition = cond; } - void SetDescription(const char *d, const char *t=0) { fDescription = d; fNameTab = t; } - void SetGroupBy(GroupBy_t b=kGroupByWeek) { fGroupBy=b; } - void SetPrimaryDate(const char *ts) { fPrimaryDate=ts; } - void SetPrimaryNumber(const char *ts) { fPrimaryNumber=ts; } - void SetSecondary(const char *ts) { fSecondary=ts; } - - Bool_t Plot(const char *value, Float_t min=0, Float_t max=-1, Float_t resolution=0) - { - TString named = fPrimaryDate; - TString named2 = fSecondary; - TString namev = value; - - TString tablev = namev(0, namev.First('.')); - TString valuev = namev(namev.First('.')+1, namev.Length()); - - TString tabled = named(0, named.First('.')); - TString valued = named(named.First('.')+1, named.Length()); - - TString query="SELECT "; - switch (fGroupBy) - { - case kNone: - case kGroupByPrimary: - query += Form("%s AS %s", valued.Data(), valued.Data()+1); - break; - case kGroupByHour: - query += Form("DATE_FORMAT(%s, '%%Y-%%m-%%d %%H:30:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1); - break; - case kGroupByNight: - query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-%%d 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1); - break; - case kGroupByWeek: - query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%x%%v') AS %s ", fPrimaryDate.Data(), valued.Data()+1); - break; - case kGroupByMonth: - query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1); - break; - case kGroupBySeason: - //query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1); - break; - case kGroupByYear: - query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-08-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1); - break; - } - - if (fGroupBy==kNone) - { - query += ", "; - query += fSecondary; - query += ", "; - query += value; - query += ", "; - query += fPrimaryNumber; - query += " "; - } - else - { - query += ", AVG("; - query += fSecondary; - query += "), AVG("; - query += value; - query += "), "; - query += fPrimaryNumber; - query += ", STD("; - query += fSecondary; - query += "), STD("; - query += value; - query += ") "; - } - - query += Form("FROM %s ", tabled.Data()); - - const Bool_t interval = !fRequestFrom.IsNull() && !fRequestTo.IsNull(); - - TString where(fCondition); - if (!where.IsNull()) - where += " AND "; - where += "RunData.fExcludedFDAKEY=1 "; - - if (interval) - { - if (!where.IsNull()) - where += " AND "; - where += Form("%s BETWEEN '%s' AND '%s' ", - fPrimaryDate.Data(), fRequestFrom.Data(), fRequestTo.Data()); - } - - // ------------------------------ - - query += fServer.GetJoins(tabled, query+" "+where); - - if (!where.IsNull()) - { - query += "WHERE "; - query += where; - } - - if (fGroupBy!=kNone) - { - query += Form("GROUP BY %s ", valued.Data()+1); - //query += Form(" HAVING COUNT(%s)=(COUNT(*)+1)/2 ", valuev.Data()); - } - query += Form("ORDER BY %s ", valued.Data()+1); - - - // ------------------------------ - - TSQLResult *res = fServer.Query(query); - if (!res) - { - cout << "ERROR - Query failed: " << query << endl; - return kFALSE; - } - - if (max>min) - PlotTable(*res, namev, min, max, resolution); - else - PlotTable(*res, namev, fPlotMin, fPlotMax, resolution); - - - delete res; - return kTRUE; - } -}; +#include "plotdb.C" void plotall(MPlot &plot) @@ -523,4 +68,10 @@ plot.SetPrimaryNumber("RunData.fRunNumber"); plot.SetSecondary("RunData.fZenithDistance"); + + MTime t(-1); + TString cond = "RunData.fExcludedFDAKEY=1 AND RunData.fRunStart BETWEEN '2004-01-01' AND '"; + cond += t.GetSqlDateTime(); + cond += "'"; + plot.SetCondition(cond); //inner camera @@ -536,5 +87,5 @@ plot.SetDescription("L2 Trigger rate before prescaler;R_{L2,U}", "Unpresc"); plot.Plot("RunData.fL2RateUnpresc", -0.5, 999.5, 5); - +/* //from signal*.root plot.SetDescription("Signal Position; [sl]", "PosMean"); @@ -555,4 +106,5 @@ plot.SetDescription("Interleaved Signal Height FWHM;\\sigma_{H} [counts]", "ISigFWHM"); plot.Plot("DataCheck.fHeightFWHMInterlaced", -0.5, 155.5, 1); +*/ }