Context Navigation

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

Visit:

Last change on this file since 8993 was 8988, checked in by tbretz, 16 years ago
* empty log message *
File size: 27.3 KB

Line
1	/* ======================================================================== *\
2	! $Name: not supported by cvs2svn $:$Id: plotdb.C,v 1.48 2008-06-30 09:36:34 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-2008
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	// Create TGraph objects
134	TGraph &gt = res.GetFieldCount()>4 ? new TGraphErrors : new TGraph;
135	gt.SetNameTitle(Form("%s_time", res.GetFieldName(2)), Form("%s vs %s", res.GetFieldName(2), res.GetFieldName(0)));
136	gt.SetMarkerStyle(kFullDotMedium);
137
138	TGraph gz;
139	gz.SetNameTitle(res.GetFieldName(2), Form("%s vs %s", res.GetFieldName(2), res.GetFieldName(1)));
140	gz.SetMarkerStyle(kFullDotMedium);
141
142	TGraph gt0, gt1;
143	gt0.SetMarkerColor(kRed);
144	gt1.SetMarkerColor(kBlue);
145	gt0.SetMarkerStyle(kFullDotLarge);
146	gt1.SetMarkerStyle(kFullDotLarge);
147
148	TGraph gz0, gz1;
149	gz0.SetMarkerColor(kRed);
150	gz1.SetMarkerColor(kBlue);
151	gz0.SetMarkerStyle(kFullDotLarge);
152	gz1.SetMarkerStyle(kFullDotLarge);
153
154	Int_t first = -1;
155	Int_t last = -1;
156
157	// Loop over the data
158	TSQLRow *row = 0;
159	while ((row=res.Next()))
160	{
161	// Get all fields of this row
162	const char date = (row)[0];
163	const char zd = (row)[1];
164	const char val = (row)[2];
165	const char snum = res.GetFieldCount()>3 ? (row)[3] : 0;
166	const char verr = res.GetFieldCount()>4 ? (row)[5] : 0;
167
168	delete row;
169
170	if (!date \|\| !val \|\| !zd)
171	continue;
172
173	// check if date is valid
174	MTime t(date);
175	if (!t.SetSqlDateTime(date))
176	continue;
177
178	// check if it belongs to the requested MAGIC period
179	if (fRequestPeriod>0 && MAstro::GetMagicPeriod(t.GetMjd())!=fRequestPeriod)
180	continue;
181
182	// Get axis range
183	if (first<0)
184	first = TMath::Nint(TMath::Floor(t.GetMjd()));
185	last = TMath::Nint(TMath::Ceil(t.GetMjd()));
186
187	// Convert a possible key number into a integer
188	UInt_t seq = snum ? atoi(snum) : 0;
189
190	// convert primary and secondary value into floats
191	Float_t value = atof(val);
192	Float_t zenith = atof(zd);
193
194	// If a datset is given add the point to the special TGraphs
195	// used for highliting these dates
196	if (fDataSet)
197	{
198	if (fDataSet->HasOnSequence(seq))
199	{
200	gt1.SetPoint(gt1.GetN(), t.GetAxisTime(), value);
201	gz1.SetPoint(gz1.GetN(), zenith, value);
202	}
203
204	if (fDataSet->HasOffSequence(seq))
205	{
206	gt0.SetPoint(gt0.GetN(), t.GetAxisTime(), value);
207	gz0.SetPoint(gz0.GetN(), zenith, value);
208	}
209	}
210
211	// Add Data to TGraph
212	gt.SetPoint(gt.GetN(), t.GetAxisTime(), value);
213	gz.SetPoint(gz.GetN(), zenith, value);
214
215	// Set error-bar, if one
216	if (verr)
217	static_cast<TGraphErrors&>(gt).SetPointError(gt.GetN()-1, 0, atof(verr));
218	}
219
220	// If this is done earlier the plots remain empty since root 5.12/00
221	if (fmax>fmin)
222	{
223	gt.SetMinimum(fmin);
224	gt.SetMaximum(fmax);
225	gz.SetMinimum(fmin);
226	gz.SetMaximum(fmax);
227	}
228
229	gROOT->SetSelectedPad(0);
230
231	TString cname = fNameTab.IsNull() ? name(name.First('.')+2, name.Length()) : fNameTab;
232	TString ctitle = fDescription.IsNull() ? name.Data() : fDescription.Data();
233	if (ctitle.First(';')>0)
234	ctitle.Remove(ctitle.First(';'), ctitle.Length());
235
236	// Create a TCanvas or open a new tab
237	TCanvas &c = fDisplay ? fDisplay->AddTab(cname, ctitle) : *new TCanvas(cname, ctitle);
238	// Set fillcolor, remove border and divide pad
239	c.SetFillColor(kWhite);
240	c.SetFrameBorderMode(0);
241	c.SetBorderMode(0);
242	c.Divide(1,2);
243
244	// Output mean and rms to console
245	cerr << setprecision(4) << setw(10) << cname << ": ";
246	if (gt.GetN()==0)
247	{
248	cerr << " <empty>" << endl;
249	return;
250	}
251	cerr << setw(8) << gt.GetMean(2) << "+-" << setw(8) << gt.GetRMS(2) << " ";
252	if (gt0.GetN()>0 \|\| gt1.GetN()>0)
253	{
254	cerr << setw(8) << gt1.GetMean(2) << "+-" << setw(8) << gt1.GetRMS(2) << " ";
255	cerr << setw(8) << gt0.GetMean(2) << "+-" << setw(8) << gt0.GetRMS(2);
256	}
257	cerr << endl;
258
259	// format axis
260	TH1 *h = gt.GetHistogram();
261
262	const Double_t min = fHistMin>fHistMax ? h->GetMinimum()-resolution/2 : fHistMin;
263	const Double_t max = fHistMin>fHistMax ? h->GetMaximum()+resolution/2 : fHistMax;
264
265	// Create histogram
266	const Int_t n = resolution>0 ? TMath::Nint((max-min)/resolution) : 50;
267
268	TH1F hist("Hist", Form("Distribution of %s", fDescription.IsNull() ? name.Data() : fDescription.Data()), n, min, max);
269	hist.SetDirectory(0);
270
271	// Fill data into histogra,
272	for (int i=0; i<gt.GetN(); i++)
273	hist.Fill(gt.GetY()[i]);
274
275	// Format histogram
276	if (fDescription.IsNull())
277	hist.SetXTitle(name);
278	hist.SetYTitle("Counts");
279
280	TVirtualPad *pad = gPad;
281
282	// draw contants of pad 2 (counting starts at 0)
283	pad->cd(2);
284	gPad->SetBorderMode(0);
285	gPad->SetFrameBorderMode(0);
286	gPad->SetGridy();
287
288	gPad->SetLeftMargin(0.06);
289	gPad->SetRightMargin(0.06);
290	gPad->SetBottomMargin(0.08);
291
292	h->SetXTitle("Time");
293	h->SetYTitle(hist.GetXaxis()->GetTitle());
294	h->GetXaxis()->SetTimeDisplay(1);
295	h->GetYaxis()->SetTitleOffset(0.8);
296	h->GetXaxis()->SetTitleOffset(1.0);
297	h->GetXaxis()->SetLabelOffset(0.01);
298
299	// draw TGraph
300	gROOT->SetSelectedPad(0);
301	gt.DrawClone("AP");
302	if (gt0.GetN()>0)
303	gt0.DrawClone("P");
304	if (gt1.GetN()>0)
305	gt1.DrawClone("P");
306
307	// Add lines and text showing the MAGIC periods
308	TLine l;
309	TText t;
310	Int_t num=0;
311	l.SetLineStyle(kDotted);
312	l.SetLineColor(kBlue);
313	t.SetTextColor(kBlue);
314	l.SetLineWidth(1);
315	t.SetTextSize(h->GetXaxis()->GetLabelSize());
316	t.SetTextAlign(21);
317	Int_t p0 = MAstro::GetMagicPeriod(first);
318	for (Int_t p = first; p<last; p++)
319	{
320	Int_t p1 = MAstro::GetMagicPeriod(p);
321	if (p1!=p0)
322	{
323	l.DrawLine(MTime(p).GetAxisTime(), h->GetMinimum(), MTime(p).GetAxisTime(), h->GetMaximum());
324	t.DrawText(MTime(p+15).GetAxisTime(), h->GetMaximum(), Form("%d", p1));
325	num++;
326	}
327	p0 = p1;
328	}
329	if (num<4)
330	gPad->SetGridx();
331
332	// Go back to first (upper) pad, format it and divide it again
333	pad->cd(1);
334	gPad->SetBorderMode(0);
335	gPad->SetFrameBorderMode(0);
336	gPad->Divide(2,1);
337
338	TVirtualPad *pad2 = gPad;
339
340	// format left pad
341	pad2->cd(1);
342	gPad->SetBorderMode(0);
343	gPad->SetFrameBorderMode(0);
344	gPad->SetGridx();
345	gPad->SetGridy();
346
347	// plot histogram
348	hist.DrawCopy("");
349
350	// format right pad
351	pad2->cd(2);
352	gPad->SetBorderMode(0);
353	gPad->SetFrameBorderMode(0);
354	gPad->SetGridy();
355
356	// draw graph
357	gROOT->SetSelectedPad(0);
358	gz.GetXaxis()->SetTitle(res.GetFieldName(1));
359	gz.GetYaxis()->SetTitle(hist.GetXaxis()->GetTitle());
360	gz.DrawClone("AP");
361
362	if (gz0.GetN()>0)
363	gz0.DrawClone("P");
364	if (gz1.GetN()>0)
365	gz1.DrawClone("P");
366	}
367
368	public:
369	MPlot(MSQLMagic &server) : fServer(server), fDataSet(NULL),
370	fRequestPeriod(-1), fPlotMin(0), fPlotMax(-1), fHistMin(0), fHistMax(-1), fGroupBy(kNone)
371	{
372	}
373	~MPlot()
374	{
375	if (fDataSet)
376	delete fDataSet;
377	}
378	void SetDataSet(const TString filename)
379	{
380	if (fDataSet)
381	{
382	delete fDataSet;
383	fDataSet = NULL;
384	}
385	if (!filename.IsNull())
386	fDataSet = new MDataSet(filename);
387	}
388	void SetPlotRange(Float_t min, Float_t max/, Int_t n=5/) { fPlotMin = min; fPlotMax = max; }
389	void SetHistRange(Float_t min, Float_t max) { fHistMin = min; fHistMax = max; }
390	void SetRequestRange(const char from="", const char to="") { fRequestFrom = from; fRequestTo = to; }
391	void SetRequestPeriod(Int_t n=-1) { fRequestPeriod = n; }
392	void SetCondition(const char *cond="") { fCondition = cond; }
393	void SetDescription(const char d, const char t=0) { fDescription = d; fNameTab = t; }
394	void SetTabName(const char *t) { fNameTab = t; }
395	void SetGroupBy(GroupBy_t b=kGroupByWeek) { fGroupBy=b; }
396	void SetPrimaryDate(const char *ts) { fPrimaryDate=ts; }
397	void SetPrimaryNumber(const char *ts) { fPrimaryNumber=ts; }
398	void SetSecondary(const char *ts) { fSecondary=ts; }
399
400	MStatusDisplay *GetDisplay() { return fDisplay; }
401
402	Bool_t Plot(const char *value, Float_t min=0, Float_t max=-1, Float_t resolution=0)
403	{
404	TString named = fPrimaryDate;
405	TString named2 = fSecondary;
406	TString namev = value;
407
408	TString tablev = namev(0, namev.First('.'));
409	TString valuev = namev(namev.First('.')+1, namev.Length());
410
411	TString tabled = named(0, named.First('.'));
412	TString valued = named(named.First('.')+1, named.Length());
413
414	TString query="SELECT ";
415	switch (fGroupBy)
416	{
417	case kNone:
418	case kGroupByPrimary:
419	query += Form("%s AS %s", valued.Data(), valued.Data()+1);
420	break;
421	case kGroupByHour:
422	query += Form("DATE_FORMAT(%s, '%%Y-%%m-%%d %%H:30:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1);
423	break;
424	case kGroupByNight:
425	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-%%d 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1);
426	break;
427	case kGroupByWeek:
428	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%x%%v') AS %s ", fPrimaryDate.Data(), valued.Data()+1);
429	break;
430	case kGroupByMonth:
431	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1);
432	break;
433	case kGroupBySeason:
434	//query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-%%m-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1);
435	break;
436	case kGroupByYear:
437	query += Form("DATE_FORMAT(ADDDATE(%s,Interval 12 hour), '%%Y-08-15 00:00:00') AS %s ", fPrimaryDate.Data(), valued.Data()+1);
438	break;
439	}
440
441	if (fGroupBy==kNone)
442	{
443	query += ", ";
444	query += fSecondary;
445	query += ", ";
446	query += value;
447	query += ", ";
448	query += fPrimaryNumber;
449	query += " ";
450	}
451	else
452	{
453	query += ", AVG(";
454	query += fSecondary;
455	query += "), AVG(";
456	query += value;
457	query += "), ";
458	query += fPrimaryNumber;
459	query += ", STD(";
460	query += fSecondary;
461	query += "), STD(";
462	query += value;
463	query += ") ";
464	}
465
466	query += Form("FROM %s ", tabled.Data());
467
468	const Bool_t interval = !fRequestFrom.IsNull() && !fRequestTo.IsNull();
469
470	TString where(fCondition);
471	/*
472	if (!fDataSet && !interval && tablev=="Star")
473	{
474	if (!where.IsNull())
475	where += " AND ";
476	where += "Star.fMuonNumber>300 ";
477	}
478	*/
479	if (interval)
480	{
481	if (!where.IsNull())
482	where += " AND ";
483	where += Form("%s BETWEEN '%s' AND '%s' ",
484	fPrimaryDate.Data(), fRequestFrom.Data(), fRequestTo.Data());
485	}
486
487	// ------------------------------
488
489	query += fServer.GetJoins(tabled, query+" "+where);
490
491	if (!where.IsNull())
492	{
493	query += "WHERE ";
494	query += where;
495	query += " ";
496	}
497
498	if (fGroupBy!=kNone)
499	{
500	query += Form("GROUP BY %s ", valued.Data()+1);
501	//query += Form(" HAVING COUNT(%s)=(COUNT(*)+1)/2 ", valuev.Data());
502	}
503	query += Form("ORDER BY %s ", valued.Data()+1);
504
505	// ------------------------------
506
507	TSQLResult *res = fServer.Query(query);
508	if (!res)
509	{
510	cout << "ERROR - Query failed: " << query << endl;
511	return kFALSE;
512	}
513
514	if (max>min)
515	PlotTable(*res, namev, min, max, resolution);
516	else
517	PlotTable(*res, namev, fPlotMin, fPlotMax, resolution);
518
519
520	delete res;
521	return kTRUE;
522	}
523	};
524
525	void plotalldb(MPlot &plot)
526	{
527	//plot.SetGroupBy(MPlot::kGroupByNight);
528
529	plot.SetPrimaryDate("Sequences.fRunStart");
530	plot.SetPrimaryNumber("Sequences.fSequenceFirst");
531	plot.SetSecondary("(Sequences.fZenithDistanceMin+Sequences.fZenithDistanceMax)/2 as '<Zd>'");
532
533	//inner camera
534	//from calib*.root
535	plot.SetDescription("Conversion Factor inner Camera;C_{I} [phe/fadc cnts]", "ConvI");
536	plot.Plot("Calibration.fConvFactorInner", 0, 0.7, 0.002);
537	plot.SetDescription("Mean Arrival Time inner Camera;T_{I} [sl]", "ArrTmI");
538	plot.Plot("Calibration.fArrTimeMeanInner", 0, 40.0, 0.1);
539	plot.SetDescription("RMS Arrival Time inner Camera;\\sigma_{T,I} [sl]", "RmsArrTmI");
540	plot.Plot("Calibration.fArrTimeRmsInner", 0, 4.5, 0.01);
541	plot.SetDescription("Number of unsuitable pixels inner Camera;N_{I}", "UnsuitI");
542	plot.Plot("Calibration.fUnsuitableInner", 0, 75, 1);
543	plot.SetDescription("Number of unsuitable pixels >50%;N", "Unsuit50");
544	plot.Plot("Calibration.fUnsuitable50", 0, 75, 1);
545	plot.SetDescription("Number of unsuitable pixels >1%;N", "Unsuit01");
546	plot.Plot("Calibration.fUnsuitable01", 0, 200, 5);
547
548	//from signal*.root
549	plot.SetDescription("Mean Pedestal RMS inner Camera;\\sigma_{P,I} [phe]", "PedRmsI");
550	plot.Plot("Calibration.fMeanPedRmsInner", 0, 3.5, 0.05);
551	plot.SetDescription("Mean Signal inner Camera;S_{I} [phe]", "SignalI");
552	plot.Plot("Calibration.fMeanSignalInner", 0, 7.0, 0.05);
553
554	plot.SetDescription("Mean PulsePosCheck (maximum slice) inner camera;T [sl]", "ChkPos");
555	plot.Plot("Calibration.fPulsePosCheckMean", 1, 50.0, 0.1);
556	plot.SetDescription("Rms PulsePosCheck (maximum slice) inner camera;T [sl]", "ChkRms");
557	plot.Plot("Calibration.fPulsePosCheckRms", 0, 12.0, 0.1);
558	plot.SetDescription("Mean calibrated PulsePos (as extracted);T [ns]", "PulPos");
559	plot.Plot("Calibration.fPulsePosMean", 1, 40.0, 0.1);
560	plot.SetDescription("Rms calibrated PulsePos (as extracted);T [ns]", "PulRms");
561	plot.Plot("Calibration.fPulsePosRms", 0, 3.0, 0.02);
562
563	plot.SetDescription("Ratio of accepted calibration events in cal-run;f [%]", "RatioCal");
564	plot.Plot("Calibration.fRatioCalEvents", 0, 105, 5);
565
566	plot.SetDescription("Average rate of events with lvl1 trigger;R [Hz]", "RateTrig");
567	plot.Plot("Calibration.fRateTrigEvts", 0, 500, 5);
568	plot.SetDescription("Average rate of events with only Sum trigger;R [Hz]", "RateSum");
569	plot.Plot("Calibration.fRateSumEvts", 0, 1500, 5);
570	//plot.SetDescription("Ratio of only Sum to Lvl1 triggers;R [Hz]", "RateRatio");
571	//plot.Plot("100*Calibration.fRateSumEvts/Calibration.fRateTrigEvts", 0, 500, 1);
572	plot.SetDescription("Average rate of events with calibration trigger;R [Hz]", "RateCal");
573	plot.Plot("Calibration.fRateCalEvts", 0, 75, 1);
574	plot.SetDescription("Average rate of events with pedestal trigger;R [Hz]", "RatePed");
575	plot.Plot("Calibration.fRatePedEvts", 0, 50, 1);
576	plot.SetDescription("Average rate of events with ped+cosmics trigger;R [Hz]", "RatePT");
577	plot.Plot("Calibration.fRatePedTrigEvts", 0, 4.0, 0.1);
578	plot.SetDescription("Average rate of events without trigger pattern;R [Hz]", "Rate0");
579	plot.Plot("Calibration.fRateNullEvts", 0, 3.5, 0.1);
580	plot.SetDescription("Average rate of unknown trigger pattern;R [Hz]", "RateUnknown");
581	plot.Plot("Calibration.fRateUnknownEvts", 0, 3.5, 0.1);
582
583	plot.SetDescription("Hi-/Lo-Gain offset;", "PulOff");
584	plot.Plot("Calibration.fPulsePosOffMed", -0.33, 0.5, 0.01);
585	plot.SetDescription("Hi-/Lo-Gain ratio;", "HiLoRatio");
586	plot.Plot("Calibration.fHiLoGainRatioMed", 10, 15, 0.05);
587
588	//plot.SetDescription("Pulse Variance;", "PulVar");
589	//plot.Plot("Calibration.fPulsePosVar", 0, 0.03, 0.001);
590
591
592	//from star*.root
593	//muon
594	plot.SetCondition("Star.fMuonNumber>300");
595	plot.SetDescription("Point Spred Function;PSF [mm]");
596	plot.Plot("Star.fPSF", 0, 30, 0.5);
597	plot.SetDescription("Muon Calibration Ratio Data/MC;r [1]", "MuonCal");
598	plot.Plot("Star.fRatio", 0, 200, 0.5);
599	plot.SetDescription("Muon Rate after Muon Cuts;R [Hz]");
600	plot.Plot("Star.fMuonRate", 0, 2.0, 0.05);
601	plot.SetCondition();
602
603	//quality
604	plot.SetDescription("Datarate [Hz]", "Rate");
605	plot.Plot("Star.fDataRate", 0, 600, 10);
606	plot.SetDescription("Camera Inhomogeneity;\\sigma [%]", "Inhom");
607	plot.Plot("Star.fInhomogeneity", 0, 100, 1);
608	plot.SetDescription("Camera Spark Rate;R [Hz]", "Sparks");
609	plot.Plot("Star.fSparkRate", 0.075, 2.425, 0.05);
610	plot.SetDescription("Relative effective on time;T_{eff}/T_{obs} [ratio]", "RelTime");
611	plot.Plot("Star.fEffOnTime/Sequences.fRunTime", 0.006, 1.506, 0.01);
612	//imgpar
613	plot.SetDescription("Mean Number of Islands after cleaning;N [#]", "NumIsl");
614	plot.Plot("Star.fMeanNumberIslands", 0.5, 4.5, 0.01);
615	//weather
616	plot.SetDescription("Average Temperature [\\circ C];T [\\circ C]", "Temp");
617	plot.Plot("Star.fAvgTemperature", -5, 25, 1);
618	plot.SetDescription("Average Cloudiness [%]", "AvgClouds");
619	plot.Plot("Star.fAvgCloudiness", 0, 100, 1);
620	plot.SetDescription("RMS Cloudiness [%]", "RmsClouds");
621	plot.Plot("Star.fRmsCloudiness", 0, 30, 1);
622	plot.SetDescription("Sky Temperature [K]", "SkyTemp");
623	plot.Plot("Star.fAvgTempSky", 180, 280, 1);
624	plot.SetDescription("Maximum Humidity [%]", "MaxHum");
625	plot.Plot("Star.fMaxHumidity", 0, 100, 1);
626	plot.SetDescription("Average Humidity [%]", "AvgHum");
627	plot.Plot("Star.fAvgHumidity", 0, 100, 1);
628	plot.SetDescription("Average Wind Speed [km/h];v [km/h]", "Wind");
629	plot.Plot("Star.fAvgWindSpeed", 0, 50, 1);
630
631	//muon
632	//plot.SetDescription("Number of Muons after Muon Cuts;N [#]");
633	//plot.Plot("Star.fMuonNumber", 0, 10000, 100);
634
635	// starguider
636	plot.SetDescription("Median No. Stars recognized by the starguider;N_{0}", "StarsMed");
637	plot.Plot("Star.fNumStarsMed", 0, 100, 1);
638	plot.SetDescription("RMS No. Stars recognized by the starguider;\\sigma_{N_{0}}", "StarsRMS");
639	plot.Plot("Star.fNumStarsRMS", 0, 25, 1);
640	plot.SetDescription("Median No. Stars correlated by the starguider;N", "CorMed");
641	plot.Plot("Star.fNumStarsCorMed", 0, 100, 1);
642	plot.SetDescription("RMS No. Stars correlated by the starguider;\\sigma_{N}", "CorRMS");
643	plot.Plot("Star.fNumStarsCorRMS", 0, 25, 1);
644	plot.SetDescription("Relative number of correlated stars;N/N_{0} [%]", "StarsRel");
645	plot.Plot("Star.fNumStarsCorMed/Star.fNumStarsMed*100", 0, 100, 10);
646	plot.SetDescription("Median skbrightess measured by the starguider;B [au]", "BrightMed");
647	plot.Plot("Star.fBrightnessMed", 0, 111, 1);
648	plot.SetDescription("RMS skybrightess measured by the starguider;\\sigma_{B} [au]", "BrightRMS");
649	plot.Plot("Star.fBrightnessRMS", 0, 64, 1);
650
651	//outer camera
652	//from calib*.root
653	plot.SetDescription("Conversion Factor outer Camera;C_{O} [phe/fadc cnts]", "ConvO");
654	plot.Plot("Calibration.fConvFactorOuter", 0, 3.0, 0.01);
655	plot.SetDescription("Mean Arrival Time outer Camera;T_{O} [sl]", "ArrTmO");
656	plot.Plot("Calibration.fArrTimeMeanOuter", 0, 45, 0.1);
657	plot.SetDescription("RMS Arrival Time outer Camera;\\sigma_{T,O} [sl]", "RmsArrTmO");
658	plot.Plot("Calibration.fArrTimeRmsOuter", 0, 4.5, 0.01);
659	plot.SetDescription("Number of unsuitable pixels outer Camera;N_{O}", "UnsuitO");
660	plot.Plot("Calibration.fUnsuitableOuter", 0, 25, 1);
661	//from signal*.root
662	plot.SetDescription("Mean Pedestal RMS outer Camera;\\sigma_{P,O} [phe]", "PedRmsO");
663	plot.Plot("Calibration.fMeanPedRmsOuter", 0, 4.0, 0.05);
664	plot.SetDescription("Mean Signal outer Camera;S_{O} [phe]", "SignalO");
665	plot.Plot("Calibration.fMeanSignalOuter", 0, 4.0, 0.05);
666	}
667
668	int plotdb(TString from, TString to, const char *dataset=0)
669	{
670	TEnv env("sql.rc");
671
672	MSQLMagic serv(env);
673	if (!serv.IsConnected())
674	{
675	cout << "ERROR - Connection to database failed." << endl;
676	return 0;
677	}
678
679	cout << "plotdb" << endl;
680	cout << "------" << endl;
681	cout << endl;
682	cout << "Connected to " << serv.GetName() << endl;
683	cout << endl;
684
685	MStatusDisplay *d = new MStatusDisplay;
686	d->SetWindowName(serv.GetName());
687	d->SetTitle(serv.GetName());
688
689	MPlot plot(serv);
690	plot.SetDataSet(dataset);
691	plot.SetDisplay(d);
692	plot.SetRequestRange(from, to);
693	plotalldb(plot);
694	d->SaveAsRoot("plotdb.root");
695	d->SaveAsPS("plotdb.ps");
696
697	return 1;
698	}
699
700	int plotdb(const char *ds)
701	{
702	TEnv env("sql.rc");
703
704	MSQLMagic serv(env);
705	if (!serv.IsConnected())
706	{
707	cout << "ERROR - Connection to database failed." << endl;
708	return 0;
709	}
710
711	cout << "plotdb" << endl;
712	cout << "------" << endl;
713	cout << endl;
714	cout << "Connected to " << serv.GetName() << endl;
715	cout << endl;
716
717	MStatusDisplay *d = new MStatusDisplay;
718	d->SetWindowName(serv.GetName());
719	d->SetTitle(serv.GetName());
720
721	MPlot plot(serv);
722	plot.SetDataSet(ds);
723	plot.SetDisplay(d);
724	plot.SetRequestRange("", "");
725	plotalldb(plot);
726	d->SaveAsRoot("plotdb.root");
727	d->SaveAsPS("plotdb.ps");
728
729	return 1;
730	}
731
732	int plotdb(Int_t period, const char *dataset="")
733	{
734	TEnv env("sql.rc");
735
736	MSQLMagic serv(env);
737	if (!serv.IsConnected())
738	{
739	cout << "ERROR - Connection to database failed." << endl;
740	return 0;
741	}
742
743	cout << "plotdb" << endl;
744	cout << "------" << endl;
745	cout << endl;
746	cout << "Connected to " << serv.GetName() << endl;
747	cout << endl;
748
749	MStatusDisplay *d = new MStatusDisplay;
750	d->SetWindowName(serv.GetName());
751	d->SetTitle(serv.GetName());
752
753	MPlot plot(serv);
754	plot.SetDataSet(dataset);
755	plot.SetDisplay(d);
756	plot.SetRequestPeriod(period);
757	plotalldb(plot);
758	d->SaveAsRoot("plotdb.root");
759	d->SaveAsPS("plotdb.ps");
760
761	return 1;
762	}
763
764	int plotdb()
765	{
766	return plotdb("", "");
767	}

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

Download in other formats: