Context Navigation

source: trunk/MagicSoft/Mars/datacenter/macros/plotdb.C@ 8553

Visit:

Last change on this file since 8553 was 8478, checked in by tbretz, 18 years ago
* empty log message *
File size: 24.7 KB

Line
1	/* ======================================================================== *\
2	! $Name: not supported by cvs2svn $:$Id: plotdb.C,v 1.34 2007-05-09 12:15:51 tbretz Exp $
3	! --------------------------------------------------------------------------
4	!
5	! *
6	! * This file is part of MARS, the MAGIC Analysis and Reconstruction
7	! * Software. It is distributed to you in the hope that it can be a useful
8	! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
9	! * It is distributed WITHOUT ANY WARRANTY.
10	! *
11	! * Permission to use, copy, modify and distribute this software and its
12	! * documentation for any purpose is hereby granted without fee,
13	! * provided that the above copyright notice appear in all copies and
14	! * that both that copyright notice and this permission notice appear
15	! * in supporting documentation. It is provided "as is" without express
16	! * or implied warranty.
17	! *
18	!
19	!
20	! Author(s): Thomas Bretz, 05/2005 <mailto:tbretz@astro.uni-wuerzburg.de>
21	! Author(s): Daniela Dorner, 05/2005 <mailto:dorner@astro.uni-wuerzburg.de>
22	!
23	! Copyright: MAGIC Software Development, 2000-2006
24	!
25	!
26	\* ======================================================================== */
27
28	/////////////////////////////////////////////////////////////////////////////
29	//
30	// plotdb.C
31	// ========
32	//
33	// This macro is used to read quality parameters from the DB and plot them.
34	//
35	// The parameters are from the following files:
36	// calib*.root:mean conversion factor, mean arrival time, rms arrival time
37	// (each parameter for inner and outer camera)
38	// signal*.root: mean pedestal rms (for inner and outer camera)
39	// star*.root: PSF, # of Muons, Effective OnTime, Muon rate,
40	// Ratio MC/Data(MuonSize) and mean number of islands
41	//
42	// In the DB these values are stored in the tables Calibration and Star.
43	//
44	// Usage:
45	// .x plotdb.C --> all values in the DB are plotted
46	// You can chose are certain period:
47	// .x plotdb.C(25) --> all values from period 25 are plotted
48	// or a time period from a certain date to a certain date
49	// .x plotdb.C("2004-11-14 00:00:00", "2005-02-28 00:00:00")
50	// --> all values from 14.11.2004 0h to 28.2.2005 0h are plotted
51	// or all data, but with dataset data highlighted
52	// .x plotdb.C("dataset.txt")
53	// --> the sequences defined in dataset.txt are highlighted (blue:on, red:off)
54	// --> You can also add a dataset-name as last argument to one of the
55	// calls above
56	//
57	// Make sure, that database and password are corretly set in a resource
58	// file called sql.rc and the resource file is found.
59	//
60	/////////////////////////////////////////////////////////////////////////////
61	#include <iostream>
62	#include <iomanip>
63
64	#include <TH1.h>
65	#include <TEnv.h>
66	#include <TPad.h>
67	#include <TLine.h>
68	#include <TText.h>
69	#include <TFrame.h>
70	#include <TStyle.h>
71	#include <TCanvas.h>
72	#include <TPRegexp.h>
73	#include <TSQLRow.h>
74	#include <TSQLResult.h>
75	#include <TGraphErrors.h>
76
77	#include "MTime.h"
78	#include "MAstro.h"
79	#include "MDataSet.h"
80	#include "MSQLMagic.h"
81	#include "MStatusDisplay.h"
82
83	class MPlot : public MParContainer
84	{
85	public:
86	// Possible constants to group-by (average) over a certain period
87	enum GroupBy_t
88	{
89	kNone,
90	kGroupByPrimary,
91	kGroupByHour,
92	kGroupByNight,
93	kGroupByWeek,
94	kGroupByMonth,
95	kGroupBySeason,
96	kGroupByYear
97	};
98
99	private:
100	MSQLMagic &fServer; // Reference to the sql-server class
101
102	MDataSet *fDataSet; // A possible dtaset to highlite single points
103
104	TString fPrimaryDate; // The name of the data we plot
105	TString fPrimaryNumber; // The corresponding name for the key number
106	TString fSecondary; // The value versus which the second plot is made
107
108	TString fRequestFrom; // Start of a requested date range
109	TString fRequestTo; // End of a requested date range
110	Int_t fRequestPeriod; // A possible requested period
111
112	Float_t fPlotMin;
113	Float_t fPlotMax;
114
115	Float_t fHistMin;
116	Float_t fHistMax;
117
118	TString fDescription; // The description (title) of the plot
119	TString fNameTab; // The name of the tab in the display
120
121	TString fCondition; // An additional condition added to the query
122	GroupBy_t fGroupBy; // A possible Group-By flag
123
124	// --------------------------------------------------------------------------
125	//
126	// Function to plot the result of the query
127	//
128	void PlotTable(TSQLResult &res, TString name, Float_t fmin, Float_t fmax, Float_t resolution)
129	{
130	// Enable all otions in the statistics box
131	gStyle->SetOptStat(111111);
132
133	TSQLRow *row;
134
135	// Create TGraph objects
136	TGraph &gt = res.GetFieldCount()>4 ? new TGraphErrors : new TGraph;
137	gt.SetNameTitle(name, Form("%s vs Time", name.Data()));
138	gt.SetMarkerStyle(kFullDotMedium);
139
140	TGraph gz;
141	gz.SetNameTitle(name, Form("%s vs <Zd>", name.Data()));
142	gz.SetMarkerStyle(kFullDotMedium);
143
144	TGraph gt0, gt1;
145	gt0.SetMarkerColor(kRed);
146	gt1.SetMarkerColor(kBlue);
147	gt0.SetMarkerStyle(kFullDotLarge);
148	gt1.SetMarkerStyle(kFullDotLarge);
149
150	TGraph gz0, gz1;
151	gz0.SetMarkerColor(kRed);
152	gz1.SetMarkerColor(kBlue);
153	gz0.SetMarkerStyle(kFullDotLarge);
154	gz1.SetMarkerStyle(kFullDotLarge);
155
156	Int_t first = -1;
157	Int_t last = -1;
158
159	// Loop over the data
160	while ((row=res.Next()))
161	{
162	// Get all fields of this row
163	const char date = (row)[0];
164	const char zd = (row)[1];
165	const char val = (row)[2];
166	const char snum = res.GetFieldCount()>3 ? (row)[3] : 0;
167	const char verr = res.GetFieldCount()>4 ? (row)[5] : 0;
168	if (!date \|\| !val \|\| !zd)
169	continue;
170
171	// check if date is valid
172	MTime t(date);
173	if (!t.SetSqlDateTime(date))
174	continue;
175
176	// check if it belongs to the requested MAGIC period
177	if (fRequestPeriod>0 && MAstro::GetMagicPeriod(t.GetMjd())!=fRequestPeriod)
178	continue;
179
180	// Get axis range
181	if (first<0)
182	first = TMath::Nint(TMath::Floor(t.GetMjd()));
183	last = TMath::Nint(TMath::Ceil(t.GetMjd()));
184
185	// Convert a possible key number into a integer
186	UInt_t seq = snum ? atoi(snum) : 0;
187
188	// convert primary and secondary value into floats
189	Float_t value = atof(val);
190	Float_t zenith = atof(zd);
191
192	// If a datset is given add the point to the special TGraphs
193	// used for highliting these dates
194	if (fDataSet)
195	{
196	if (fDataSet->HasOnSequence(seq))
197	{
198	gt1.SetPoint(gt1.GetN(), t.GetAxisTime(), value);
199	gz1.SetPoint(gz1.GetN(), zenith, value);
200	}
201
202	if (fDataSet->HasOffSequence(seq))
203	{
204	gt0.SetPoint(gt0.GetN(), t.GetAxisTime(), value);
205	gz0.SetPoint(gz0.GetN(), zenith, value);
206	}
207	}
208
209	// Add Data to TGraph
210	gt.SetPoint(gt.GetN(), t.GetAxisTime(), value);
211	gz.SetPoint(gz.GetN(), zenith, value);
212
213	// Set error-bar, if one
214	if (verr)
215	static_cast<TGraphErrors&>(gt).SetPointError(gt.GetN()-1, 0, atof(verr));
216	}
217
218	// If this is done earlier the plots remain empty since root 5.12/00
219	if (fmax>fmin)
220	{
221	gt.SetMinimum(fmin);
222	gt.SetMaximum(fmax);
223	gz.SetMinimum(fmin);
224	gz.SetMaximum(fmax);
225	}
226
227	gROOT->SetSelectedPad(0);
228
229	// Create a TCanvas or open a new tab
230	TString title = fNameTab.IsNull() ? name(name.First('.')+2, name.Length()) : fNameTab;
231	TCanvas &c = fDisplay ? fDisplay->AddTab(title) : *new TCanvas;
232	// Set fillcolor, remove border and divide pad
233	c.SetFillColor(kWhite);
234	c.SetBorderMode(0);
235	c.Divide(1,2);
236
237	// Output mean and rms to console
238	cerr << setprecision(4) << setw(10) << title << ": ";
239	if (gt.GetN()==0)
240	{
241	cerr << " <empty>" << endl;
242	return;
243	}
244	cerr << setw(8) << gt.GetMean(2) << "+-" << setw(8) << gt.GetRMS(2) << " ";
245	if (gt0.GetN()>0 \|\| gt1.GetN()>0)
246	{
247	cerr << setw(8) << gt1.GetMean(2) << "+-" << setw(8) << gt1.GetRMS(2) << " ";
248	cerr << setw(8) << gt0.GetMean(2) << "+-" << setw(8) << gt0.GetRMS(2);
249	}
250	cerr << endl;
251
252	TVirtualPad *pad = gPad;
253
254	// draw contants of pad 2 (counting starts at 0)
255	pad->cd(2);
256	gPad->SetBorderMode(0);
257	gPad->SetFrameBorderMode(0);
258	gPad->SetGridy();
259
260	gPad->SetLeftMargin(0.06);
261	gPad->SetRightMargin(0.06);
262	gPad->SetBottomMargin(0.08);
263
264	// format axis
265	TH1 *h = gt.GetHistogram();
266
267	h->SetXTitle("Time");
268	h->SetYTitle(name);
269	h->GetXaxis()->SetTimeDisplay(1);
270	h->GetYaxis()->SetTitleOffset(0.8);
271	h->GetXaxis()->SetTitleOffset(1.0);
272	h->GetXaxis()->SetLabelOffset(0.01);
273
274	// draw TGraph
275	gt.DrawClone("AP");
276	if (gt0.GetN()>0)
277	gt0.DrawClone("P");
278	if (gt1.GetN()>0)
279	gt1.DrawClone("P");
280
281	// Add lines and text showing the MAGIC periods
282	TLine l;
283	TText t;
284	Int_t num=0;
285	l.SetLineStyle(kDotted);
286	l.SetLineColor(kBlue);
287	t.SetTextColor(kBlue);
288	l.SetLineWidth(1);
289	t.SetTextSize(h->GetXaxis()->GetLabelSize());
290	t.SetTextAlign(21);
291	Int_t p0 = MAstro::GetMagicPeriod(first);
292	for (Int_t p = first; p<last; p++)
293	{
294	Int_t p1 = MAstro::GetMagicPeriod(p);
295	if (p1!=p0)
296	{
297	l.DrawLine(MTime(p).GetAxisTime(), h->GetMinimum(), MTime(p).GetAxisTime(), h->GetMaximum());
298	t.DrawText(MTime(p+15).GetAxisTime(), h->GetMaximum(), Form("%d", p1));
299	num++;
300	}
301	p0 = p1;
302	}
303	if (num<4)
304	gPad->SetGridx();
305
306	const Double_t min = fHistMin>fHistMax ? h->GetMinimum()-resolution/2 : fHistMin;
307	const Double_t max = fHistMin>fHistMax ? h->GetMaximum()+resolution/2 : fHistMax;
308
309	// Use this to save the pad with the time development to a file
310	//gPad->SaveAs(Form("plotdb-%s.eps", title.Data()));
311
312	// Go back to first (upper) pad, format it and divide it again
313	pad->cd(1);
314	gPad->SetBorderMode(0);
315	gPad->SetFrameBorderMode(0);
316	gPad->Divide(2,1);
317
318	TVirtualPad *pad2 = gPad;
319
320	// format left pad
321	pad2->cd(1);
322	gPad->SetBorderMode(0);
323	gPad->SetFrameBorderMode(0);
324	gPad->SetGridx();
325	gPad->SetGridy();
326
327	// Create histogram
328	const Int_t n = resolution>0 ? TMath::Nint((max-min)/resolution) : 50;
329
330	TH1F hist("Hist", Form("Distribution of %s", fDescription.IsNull() ? name.Data() : fDescription.Data()), n, min, max);
331	hist.SetDirectory(0);
332
333	// Fill data into histogra,
334	for (int i=0; i<gt.GetN(); i++)
335	hist.Fill(gt.GetY()[i]);
336
337	// Format histogram
338	if (fDescription.IsNull())
339	hist.SetXTitle(name);
340	hist.SetYTitle("Counts");
341
342	// plot histogram
343	hist.DrawCopy("");
344
345	// format right pad
346	pad2->cd(2);
347	gPad->SetBorderMode(0);
348	gPad->SetFrameBorderMode(0);
349	gPad->SetGridy();
350
351	// format graph
352	TH1 *h2 = gz.GetHistogram();
353
354	h2->SetXTitle("Zd");
355	h2->SetYTitle(name);
356
357	// draw graph
358	gz.DrawClone("AP");
359
360	if (gz0.GetN()>0)
361	gz0.DrawClone("P");
362	if (gz1.GetN()>0)
363	gz1.DrawClone("P");
364	}
365
366	public:
367	MPlot(MSQLMagic &server) : fServer(server), fDataSet(NULL),
368	fRequestPeriod(-1), fPlotMin(0), fPlotMax(-1), fHistMin(0), fHistMax(-1), fGroupBy(kNone)
369	{
370	}
371	~MPlot()
372	{
373	if (fDataSet)
374	delete fDataSet;
375	}
376	void SetDataSet(const TString filename)
377	{
378	if (fDataSet)
379	{
380	delete fDataSet;
381	fDataSet = NULL;
382	}
383	if (!filename.IsNull())
384	fDataSet = new MDataSet(filename);
385	}
386	void SetPlotRange(Float_t min, Float_t max, Int_t n=5) { fPlotMin = min; fPlotMax = max; }
387	void SetHistRange(Float_t min, Float_t max) { fHistMin = min; fHistMax = max; }
388	void SetRequestRange(const char from="", const char to="") { fRequestFrom = from; fRequestTo = to; }
389	void SetRequestPeriod(Int_t n=-1) { fRequestPeriod = n; }
390	void SetCondition(const char *cond="") { fCondition = cond; }
391	void SetDescription(const char d, const char t=0) { fDescription = d; fNameTab = t; }
392	void SetGroupBy(GroupBy_t b=kGroupByWeek) { fGroupBy=b; }
393	void SetPrimaryDate(const char *ts) { fPrimaryDate=ts; }
394	void SetPrimaryNumber(const char *ts) { fPrimaryNumber=ts; }
395	void SetSecondary(const char *ts) { fSecondary=ts; }
396
397	Bool_t Plot(const char *value, Float_t min=0, Float_t max=-1, Float_t resolution=0)
398	{
399	TString named = fPrimaryDate;
400	TString named2 = fSecondary;
401	TString namev = value;
402
403	TString tablev = namev(0, namev.First('.'));
404	TString valuev = namev(namev.First('.')+1, namev.Length());
405
406	TString tabled = named(0, named.First('.'));
407	TString valued = named(named.First('.')+1, named.Length());
408
409	TString query="SELECT ";
410	switch (fGroupBy)
411	{
412	case kNone:
413	case kGroupByPrimary:
414	query += valued;
415	break;
416	case kGroupByHour:
417	query += Form("DATE_FORMAT(%s, '%%Y-%%m-%%d %%H:30:00') AS %s ", fPrimaryDate.Data(), valued.Data());
418	break;
419	case kGroupByNight:
420	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-%%d 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data());
421	break;
422	case kGroupByWeek:
423	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%x%%v') AS %s ", fPrimaryDate.Data(), valued.Data());
424	break;
425	case kGroupByMonth:
426	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data());
427	break;
428	case kGroupBySeason:
429	//query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data());
430	break;
431	case kGroupByYear:
432	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-08-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data());
433	break;
434	}
435
436	if (fGroupBy==kNone)
437	{
438	query += ", ";
439	query += fSecondary;
440	query += ", ";
441	query += value;
442	query += ", ";
443	query += fPrimaryNumber;
444	query += " ";
445	}
446	else
447	{
448	query += ", AVG(";
449	query += fSecondary;
450	query += "), AVG(";
451	query += value;
452	query += "), ";
453	query += fPrimaryNumber;
454	query += ", STD(";
455	query += fSecondary;
456	query += "), STD(";
457	query += value;
458	query += ") ";
459	}
460
461	query += Form("FROM %s ", tabled.Data());
462
463	const Bool_t interval = !fRequestFrom.IsNull() && !fRequestTo.IsNull();
464
465	TString where(fCondition);
466	if (!fDataSet && !interval && tablev=="Star")
467	{
468	if (!where.IsNull())
469	where += " AND ";
470	where += "Star.fMuonNumber>300 ";
471	}
472
473	if (interval)
474	{
475	if (!where.IsNull())
476	where += " AND ";
477	where += Form("%s BETWEEN '%s' AND '%s' ",
478	fPrimaryDate.Data(), fRequestFrom.Data(), fRequestTo.Data());
479	}
480
481	// ------------------------------
482
483	query += fServer.GetJoins(tabled, query+" "+where);
484
485	if (!where.IsNull())
486	{
487	query += "WHERE ";
488	query += where;
489	}
490
491	if (fGroupBy!=kNone)
492	{
493	query += Form("GROUP BY %s ", valued.Data());
494	//query += Form(" HAVING COUNT(%s)=(COUNT(*)+1)/2 ", valuev.Data());
495	}
496	query += Form("ORDER BY %s ", valued.Data());
497
498
499	// ------------------------------
500
501	TSQLResult *res = fServer.Query(query);
502	if (!res)
503	{
504	cout << "ERROR - Query failed: " << query << endl;
505	return kFALSE;
506	}
507
508	if (max>min)
509	PlotTable(*res, namev, min, max, resolution);
510	else
511	PlotTable(*res, namev, fPlotMin, fPlotMax, resolution);
512
513
514	delete res;
515	return kTRUE;
516	}
517	};
518
519	void plotall(MPlot &plot)
520	{
521	//plot.SetGroupBy(MPlot::kGroupByWeek);
522
523	plot.SetPrimaryDate("Sequences.fRunStart");
524	plot.SetPrimaryNumber("Sequences.fSequenceFirst");
525	plot.SetSecondary("(Sequences.fZenithDistanceMin+Sequences.fZenithDistanceMax)/2");
526
527	//inner camera
528	//from calib*.root
529	plot.SetDescription("Conversion Factor inner Camera;C_{I} [phe/fadc cnts]", "ConvI");
530	plot.Plot("Calibration.fConvFactorInner", 0, 0.7, 0.002);
531	plot.SetDescription("Mean Arrival Time inner Camera;T_{I} [sl]", "ArrTmI");
532	plot.Plot("Calibration.fArrTimeMeanInner", 0, 40.0, 0.1);
533	plot.SetDescription("RMS Arrival Time inner Camera;\\sigma_{T,I} [sl]", "RmsArrTmI");
534	plot.Plot("Calibration.fArrTimeRmsInner", 0, 3.0, 0.01);
535	plot.SetDescription("Number of unsuitable pixels inner Camera;N{I}", "UnsuitI");
536	plot.Plot("Calibration.fUnsuitableInner", 0, 75, 1);
537	plot.SetDescription("Number of unsuitable pixels >50%;N", "Unsuit50");
538	plot.Plot("Calibration.fUnsuitable50", 0, 75, 1);
539	plot.SetDescription("Number of unsuitable pixels >1%;N", "Unsuit01");
540	plot.Plot("Calibration.fUnsuitable01", 0, 75, 1);
541
542	//from signal*.root
543	plot.SetDescription("Mean Pedestal RMS inner Camera;\\sigma_{P,I} [phe]", "PedRmsI");
544	plot.Plot("Calibration.fMeanPedRmsInner", 0, 3.5, 0.05);
545	plot.SetDescription("Mean Signal inner Camera;S_{I} [phe]", "SignalI");
546	plot.Plot("Calibration.fMeanSignalInner", 0, 7.0, 0.05);
547	plot.SetDescription("Mean PulsePosCheck (falling edge) inner camera;T [sl]", "ChkPos");
548	plot.Plot("Calibration.fPulsePosCheckMean", 1, 40.0, 0.1);
549	plot.SetDescription("Rms PulsePosCheck (falling edge) inner camera;T [sl]", "ChkRms");
550	plot.Plot("Calibration.fPulsePosCheckRms", 0, 5.0, 0.1);
551	plot.SetDescription("Mean calibrated PulsePos;T", "PulPos");
552	plot.Plot("Calibration.fPulsePosMean", 1, 40.0, 0.1);
553	plot.SetDescription("Rms calibrated PulsePos;T", "PulRms");
554	plot.Plot("Calibration.fPulsePosRms", 0, 3.0, 0.02);
555
556	plot.SetDescription("Hi-/Lo-Gain offset;", "PulOff");
557	plot.Plot("Calibration.fPulsePosOffMed", -0.33, 0.33, 0.01);
558	plot.SetDescription("Hi-/Lo-Gain ratio;", "HiLoRatio");
559	plot.Plot("Calibration.fHiLoGainRatioMed", 10, 12.5, 0.05);
560
561	//plot.SetDescription("Pulse Variance;", "PulVar");
562	//plot.Plot("Calibration.fPulsePosVar", 0, 0.03, 0.001);
563
564	//from star*.root
565	//muon
566	plot.SetDescription("Point Spred Function;PSF [mm]");
567	plot.Plot("Star.fPSF", 0, 30, 0.5);
568	plot.SetDescription("Muon Calibration Ratio Data/MC;r [1]", "MuonCal");
569	plot.Plot("Star.fRatio", 0, 200, 0.5);
570	plot.SetDescription("Muon Rate after Muon Cuts;R [Hz]");
571	plot.Plot("Star.fMuonRate", 0, 2.0, 0.05);
572	//quality
573	plot.SetDescription("Camera Inhomogeneity;\\sigma [%]", "Inhom");
574	plot.Plot("Star.fInhomogeneity", 0, 100, 1);
575	plot.SetDescription("Camera Spark Rate;R [Hz]", "Sparks");
576	plot.Plot("Star.fSparkRate", 0.075, 2.425, 0.05);
577	//imgpar
578	plot.SetDescription("Mean Number of Islands after cleaning;N [#]", "NumIsl");
579	plot.Plot("Star.fMeanNumberIslands", 0.5, 4.5, 0.01);
580	plot.SetDescription("Measures effective on time;T_{eff} [s]", "EffOn");
581	plot.Plot("Star.fEffOnTime", 0, 10000, 150);
582	plot.SetDescription("Relative effective on time;T_{eff}/T_{obs} [ratio]", "RelTime");
583	plot.Plot("Star.fEffOnTime/Sequences.fRunTime", 0.006, 1.506, 0.01);
584	plot.SetDescription("Datarate [Hz]", "Rate");
585	plot.Plot("Star.fDataRate", 0, 600, 10);
586	plot.SetDescription("Maximum Humidity [%]", "Hum");
587	plot.Plot("Star.fMaxHumidity", 0, 100, 1);
588
589	//muon
590	//plot.SetDescription("Number of Muons after Muon Cuts;N [#]");
591	//plot.Plot("Star.fMuonNumber", 0, 10000, 100);
592
593	// starguider
594	plot.SetDescription("Median No. Stars recognized by the starguider;N_{0}", "StarsMed");
595	plot.Plot("Star.fNumStarsMed", 0, 100, 1);
596	plot.SetDescription("RMS No. Stars recognized by the starguider;\\sigma_{N_{0}}", "StarsRMS");
597	plot.Plot("Star.fNumStarsRMS", 0, 25, 1);
598	plot.SetDescription("Median No. Stars correlated by the starguider;N", "CorMed");
599	plot.Plot("Star.fNumStarsCorMed", 0, 100, 1);
600	plot.SetDescription("RMS No. Stars correlated by the starguider;\\sigma_{N}", "CorRMS");
601	plot.Plot("Star.fNumStarsCorRMS", 0, 25, 1);
602	plot.SetDescription("Relative number of correlated stars;N/N_{0} [%]", "StarsRel");
603	plot.Plot("Star.fNumStarsCorMed/Star.fNumStarsMed*100", 0, 100, 10);
604	plot.SetDescription("Median skbrightess measured by the starguider;B [au]", "BrightMed");
605	plot.Plot("Star.fBrightnessMed", 0, 111, 1);
606	plot.SetDescription("RMS skybrightess measured by the starguider;\\sigma_{B} [au]", "BrightRMS");
607	plot.Plot("Star.fBrightnessRMS", 0, 64, 1);
608
609	//outer camera
610	//from calib*.root
611	plot.SetDescription("Conversion Factor outer Camera;C_{O} [phe/fadc cnts]", "ConvO");
612	plot.Plot("Calibration.fConvFactorOuter", 0, 3.0, 0.01);
613	plot.SetDescription("Mean Arrival Time outer Camera;T_{O} [sl]", "ArrTmO");
614	plot.Plot("Calibration.fArrTimeMeanOuter", 0, 45, 0.1);
615	plot.SetDescription("RMS Arrival Time outer Camera;\\sigma_{T,O} [sl]", "RmsArrTmO");
616	plot.Plot("Calibration.fArrTimeRmsOuter", 0, 3.0, 0.01);
617	plot.SetDescription("Number of unsuitable pixels outer Camera;N{O}", "UnsuitO");
618	plot.Plot("Calibration.fUnsuitableOuter", 0, 25, 1);
619	//from signal*.root
620	plot.SetDescription("Mean Pedestal RMS outer Camera;\\sigma_{P,O} [phe]", "PedRmsO");
621	plot.Plot("Calibration.fMeanPedRmsOuter", 0, 4.0, 0.05);
622	plot.SetDescription("Mean Signal outer Camera;S_{O} [phe]", "SignalO");
623	plot.Plot("Calibration.fMeanSignalOuter", 0, 4.0, 0.05);
624	}
625
626	int plotdb(TString from, TString to, const char *dataset=0)
627	{
628	TEnv env("sql.rc");
629
630	MSQLMagic serv(env);
631	if (!serv.IsConnected())
632	{
633	cout << "ERROR - Connection to database failed." << endl;
634	return 0;
635	}
636
637	cout << "plotdb" << endl;
638	cout << "------" << endl;
639	cout << endl;
640	cout << "Connected to " << serv.GetName() << endl;
641	cout << endl;
642
643	MStatusDisplay *d = new MStatusDisplay;
644	d->SetWindowName(serv.GetName());
645	d->SetTitle(serv.GetName());
646
647	MPlot plot(serv);
648	plot.SetDataSet(dataset);
649	plot.SetDisplay(d);
650	plot.SetRequestRange(from, to);
651	plotall(plot);
652	d->SaveAsRoot("plotdb.root");
653	d->SaveAsPS("plotdb.ps");
654
655	return 1;
656	}
657
658	int plotdb(const char *ds)
659	{
660	TEnv env("sql.rc");
661
662	MSQLMagic serv(env);
663	if (!serv.IsConnected())
664	{
665	cout << "ERROR - Connection to database failed." << endl;
666	return 0;
667	}
668
669	cout << "plotdb" << endl;
670	cout << "------" << endl;
671	cout << endl;
672	cout << "Connected to " << serv.GetName() << endl;
673	cout << endl;
674
675	MStatusDisplay *d = new MStatusDisplay;
676	d->SetWindowName(serv.GetName());
677	d->SetTitle(serv.GetName());
678
679	MPlot plot(serv);
680	plot.SetDataSet(ds);
681	plot.SetDisplay(d);
682	plot.SetRequestRange("", "");
683	plotall(plot);
684	d->SaveAsRoot("plotdb.root");
685	d->SaveAsPS("plotdb.ps");
686
687	return 1;
688	}
689
690	int plotdb(Int_t period, const char *dataset="")
691	{
692	TEnv env("sql.rc");
693
694	MSQLMagic serv(env);
695	if (!serv.IsConnected())
696	{
697	cout << "ERROR - Connection to database failed." << endl;
698	return 0;
699	}
700
701	cout << "plotdb" << endl;
702	cout << "------" << endl;
703	cout << endl;
704	cout << "Connected to " << serv.GetName() << endl;
705	cout << endl;
706
707	MStatusDisplay *d = new MStatusDisplay;
708	d->SetWindowName(serv.GetName());
709	d->SetTitle(serv.GetName());
710
711	MPlot plot(serv);
712	plot.SetDataSet(dataset);
713	plot.SetDisplay(d);
714	plot.SetRequestPeriod(period);
715	plotall(plot);
716	d->SaveAsRoot("plotdb.root");
717	d->SaveAsPS("plotdb.ps");
718
719	return 1;
720	}
721
722	int plotdb()
723	{
724	return plotdb("", "");
725	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats: