Jpp  15.0.2
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JPlot1D.cc
Go to the documentation of this file.
1 #include <string>
2 #include <iostream>
3 #include <iomanip>
4 #include <vector>
5 #include <set>
6 #include <cmath>
7 
8 #include "TROOT.h"
9 #include "TFile.h"
10 #include "TClass.h"
11 #include "TApplication.h"
12 #include "TCanvas.h"
13 #include "TKey.h"
14 #include "TStyle.h"
15 #include "TAttMarker.h"
16 #include "TAttLine.h"
17 #include "TH1.h"
18 #include "TH2.h"
19 #include "TH3.h"
20 #include "TF1.h"
21 #include "TF2.h"
22 #include "THStack.h"
23 #include "TGraph.h"
24 #include "TGraphErrors.h"
25 #include "TMultiGraph.h"
26 #include "TProfile.h"
27 #include "TEllipse.h"
28 #include "TMarker.h"
29 #include "TLine.h"
30 #include "TLegend.h"
31 #include "TString.h"
32 #include "TRegexp.h"
33 #include "TText.h"
34 
35 #include "JTools/JRange.hh"
36 #include "JLang/JSinglePointer.hh"
37 #include "JROOT/JStyle.hh"
38 #include "JROOT/JCanvas.hh"
39 #include "JROOT/JMarkerAttributes.hh"
40 #include "JROOT/JLineAttributes.hh"
41 #include "JROOT/JLegend.hh"
42 #include "JGizmo/JRootObjectID.hh"
43 #include "JGizmo/JRootObject.hh"
44 #include "JGizmo/JGizmoToolkit.hh"
45 
46 #include "Jeep/JPrint.hh"
47 #include "Jeep/JParser.hh"
48 #include "Jeep/JMessage.hh"
49 
50 namespace {
51 
52  const char* const JName_t = "?"; //!< Draw histogram name as title of plot
53  const char* const JTitle_t = "%"; //!< Draw histogram title as title of plot
54 }
55 
56 
57 /**
58  * \file
59  * General purpose plot program for 1D ROOT objects.
60  * The option <tt>-f</tt> corresponds to <tt><file name>:<object name></tt>.
61  * \author mdejong
62  */
63 int main(int argc, char **argv)
64 {
65  using namespace std;
66  using namespace JPP;
67 
68  typedef JRange<double> JRange_t;
69 
70  vector<JRootObjectID> inputFile;
71  string outputFile;
72  JCanvas canvas;
73  int stats;
74  string legend;
75  JRange_t X;
76  JRange_t Y;
77  JRange_t Z;
78  JCounter logx;
79  bool logy;
80  bool logz;
81  char project;
82  string xLabel;
83  string yLabel;
84  JCounter drawLine;
85  bool fillArea;
86  int lineWidth;
87  double markerSize;
88  string option;
89  set<char> grid;
90  bool batch;
91  string title;
92  pair<string, int> Ndivisions;
93  int group;
94  int debug;
95  string xTimeFormat;
96 
97  try {
98 
99  JParser<> zap("General purpose plot program for 1D ROOT objects.");
100 
101  zap['f'] = make_field(inputFile, "<input file>:<object name>");
102  zap['o'] = make_field(outputFile, "graphics output") = "";
103  zap['w'] = make_field(canvas, "size of canvas <nx>x<ny> [pixels]") = JCanvas(500, 500);
104  zap['s'] = make_field(stats) = -1;
105  zap['L'] = make_field(legend, "position legend e.g. TR") = "", "TL", "TR", "BR", "BL";
106  zap['x'] = make_field(X, "abscissa range") = JRange_t();
107  zap['y'] = make_field(Y, "ordinate range") = JRange_t();
108  zap['z'] = make_field(Z, "ordinate range of projection)") = JRange_t();
109  zap['X'] = make_field(logx, "logarithmic x-axis (-XX log10 axis)");
110  zap['Y'] = make_field(logy, "logarithmic y-axis");
111  zap['Z'] = make_field(logz, "logarithmic y-axis; after projection");
112  zap['P'] = make_field(project, "projection") = '\0', 'x', 'X', 'y', 'Y';
113  zap['>'] = make_field(xLabel, "x-axis label") = "";
114  zap['^'] = make_field(yLabel, "y-axis label") = "";
115  zap['C'] = make_field(drawLine);
116  zap['F'] = make_field(fillArea);
117  zap['l'] = make_field(lineWidth, "line width") = 2;
118  zap['S'] = make_field(markerSize, "marker size") = 1.0;
119  zap['O'] = make_field(option, "plotting option") = "";
120  zap['G'] = make_field(grid, "grid lines [X][Y]") = JPARSER::initialised();
121  zap['B'] = make_field(batch, "batch processing");
122  zap['T'] = make_field(title, "title") = "KM3NeT preliminary";
123  zap['N'] = make_field(Ndivisions, "axis divisioning (e.g. \"X 505\")") = JPARSER::initialised();
124  zap['g'] = make_field(group, "group colour codes of objects") = 1;
125  zap['t'] = make_field(xTimeFormat, "set time format for x-axis, e.g. \%d\\/\%m\\/\\%y%F1970-01-01 00:00:00") = "";
126  zap['d'] = make_field(debug) = 0;
127 
128  zap(argc, argv);
129  }
130  catch(const exception &error) {
131  FATAL(error.what() << endl);
132  }
133 
134 
135  gROOT->SetBatch(batch);
136 
137  TApplication* tp = new TApplication("user", NULL, NULL);
138  TCanvas* cv = new TCanvas("c1", "c1", canvas.x, canvas.y);
139 
140  JSinglePointer<TStyle> gStyle(new JStyle("gplot", cv->GetWw(), cv->GetWh()));
141 
142  gROOT->SetStyle("gplot");
143  gROOT->ForceStyle();
144 
145 
146  cv->SetFillStyle(4000);
147  cv->SetFillColor(kWhite);
148  cv->Divide(1,1);
149  cv->cd(1);
150 
151 
152  JMarkerAttributes::getInstance().setMarkerSize(markerSize);
153  JLineAttributes ::getInstance().setLineWidth (lineWidth);
154 
155 
156  Double_t xmin = numeric_limits<double>::max();
157  Double_t xmax = numeric_limits<double>::lowest();
158 
159  Double_t ymin = numeric_limits<double>::max();
160  Double_t ymax = numeric_limits<double>::lowest();
161 
162 
163  vector<JRootObject> listOfObjects;
164 
165  const bool px = (project == 'x' || project == 'X'); // projection on x-axis of (2|3)D histogram
166  const bool py = (project == 'y' || project == 'Y'); // projection on y-axis of (2|3)D histogram
167  const bool pz = (project == 'z' || project == 'Z'); // projection on z-axis of 3D histogram
168 
169  logy = (logy || logz);
170 
171  if (px) {
172  swap(Y, Z); // Y becomes range in TH2::ProjectionX() and Z becomes y-axis range
173  }
174 
175  if (py) {
176  swap(X, Z); // X becomes range in TH2::ProjectionY()
177  swap(Y, X); // Y becomes x-axis range and Z becomes y-axis range
178  }
179 
180  TH1* master = NULL;
181 
182  for (vector<JRootObjectID>::const_iterator input = inputFile.begin(); input != inputFile.end(); ++input) {
183 
184  DEBUG("Input: " << *input << endl);
185 
186  TDirectory* dir = getDirectory(*input);
187 
188  if (dir == NULL) {
189  ERROR("File: " << input->getFullFilename() << " not opened." << endl);
190  continue;
191  }
192 
193  const TRegexp regexp(input->getObjectName());
194 
195  TIter iter(dir->GetListOfKeys());
196 
197  for (TKey* key; (key = (TKey*) iter.Next()) != NULL; ) {
198 
199  const TString tag(key->GetName());
200 
201  DEBUG("Key: " << tag << " match = " << tag.Contains(regexp) << endl);
202 
203  // option match
204 
205  if (tag.Contains(regexp)) {
206 
207  if (title == JName_t) {
208  title = key->GetName();
209  } else if (title == JTitle_t) {
210  title = key->GetTitle();
211  }
212 
213  JRootObject object(key->ReadObj());
214 
215  TAttMarker marker = JMarkerAttributes::getInstance().get(0);
216  TAttLine line = JLineAttributes ::getInstance().get(0);
217 
218  if (group > 1)
219  marker.SetMarkerColor(JMarkerAttributes::getInstance().get(listOfObjects.size()/group).GetMarkerColor());
220  else
221  marker = JMarkerAttributes::getInstance().get(listOfObjects.size());
222 
223  if (drawLine == 1)
224  line = JLineAttributes::getInstance().get(listOfObjects.size());
225  else
226  line.SetLineColor(marker.GetMarkerColor());
227 
228  // projections
229 
230  try {
231 
232  TProfile& h1 = dynamic_cast<TProfile&>(*object);
233 
234  object = h1.ProjectionX();
235  }
236  catch(exception&) {}
237 
238  try {
239 
240  TH2& h2 = dynamic_cast<TH2&>(*object);
241 
242  if (px) {
243 
244  if (Z != JRange_t())
245  object = h2.ProjectionX(MAKE_CSTRING(h2.GetName() << "_px" << LABEL_TERMINATOR << listOfObjects.size()),
246  h2.GetYaxis()->FindBin(Z.getLowerLimit()),
247  h2.GetYaxis()->FindBin(Z.getUpperLimit()) - 1);
248  else
249  object = h2.ProjectionX(MAKE_CSTRING(h2.GetName() << "_px" << LABEL_TERMINATOR << listOfObjects.size()),
250  1,
251  h2.GetYaxis()->GetNbins());
252 
253  } else if (py) {
254 
255  if (Z != JRange_t())
256  object = h2.ProjectionY(MAKE_CSTRING(h2.GetName() << "_py" << LABEL_TERMINATOR << listOfObjects.size()),
257  h2.GetXaxis()->FindBin(Z.getLowerLimit()),
258  h2.GetXaxis()->FindBin(Z.getUpperLimit()) - 1);
259  else
260  object = h2.ProjectionY(MAKE_CSTRING(h2.GetName() << "_py" << LABEL_TERMINATOR << listOfObjects.size()),
261  1,
262  h2.GetXaxis()->GetNbins());
263 
264  } else {
265 
266  ERROR("For 2D histograms, use option option -P for projections or use JPlot2D" << endl);
267 
268  continue;
269  }
270  }
271  catch(exception&) {}
272 
273  try {
274 
275  TH3& h3 = dynamic_cast<TH3&>(*object);
276 
277  if (px) {
278 
279  object = h3.ProjectionX(MAKE_CSTRING(h3.GetName() << "_px" << LABEL_TERMINATOR << listOfObjects.size()));
280 
281  } else if (py) {
282 
283  object = h3.ProjectionY(MAKE_CSTRING(h3.GetName() << "_py" << LABEL_TERMINATOR << listOfObjects.size()));
284 
285  } else if (pz) {
286 
287  object = h3.ProjectionZ(MAKE_CSTRING(h3.GetName() << "_pz" << LABEL_TERMINATOR << listOfObjects.size()));
288 
289  } else {
290 
291  ERROR("For 3D histograms, use option option -P for projections or use JPlot2D -P <projection>" << endl);
292 
293  continue;
294  }
295  }
296  catch(exception&) {}
297 
298  // colouring
299 
300  try {
301  if (dynamic_cast<TMarker*>(object.get()) == NULL) {
302  dynamic_cast<TAttMarker&>(*object) = marker;
303  } }
304  catch(exception&) {}
305 
306  try {
307  if (dynamic_cast<TLine*>(object.get()) == NULL) {
308  dynamic_cast<TAttLine&> (*object) = line;
309  }
310  }
311  catch(exception&) {}
312 
313  if (fillArea) {
314 
315  try {
316 
317  TAttFill& fill = dynamic_cast<TAttFill&>(*object);
318 
319  fill.SetFillColor(marker.GetMarkerColor());
320  }
321  catch(exception&) {}
322  }
323 
324  // set errors
325 
326  if (drawLine) {
327 
328  try {
329 
330  TH1& h1 = dynamic_cast<TH1&>(*object);
331 
332  for (int i = 1; i <= h1.GetNbinsX(); ++i) {
333  h1.SetBinError(i, 0.0);
334  }
335  }
336  catch(exception&) {}
337 
338  try {
339 
340  TGraphErrors& g1 = dynamic_cast<TGraphErrors&>(*object);
341 
342  for (Int_t i = 0; i != g1.GetN(); ++i) {
343  g1.GetEX()[i] = 0.0;
344  g1.GetEY()[i] = 0.0;
345  }
346  }
347  catch(exception&) {}
348  }
349 
350  // min-max
351 
352  try {
353 
354  TH1& h1 = dynamic_cast<TH1&>(*object);
355 
356  h1.SetStats(stats != -1);
357 
358  xmin = min(xmin, h1.GetXaxis()->GetXmin());
359  xmax = max(xmax, h1.GetXaxis()->GetXmax());
360  ymin = min(ymin, logy ? h1.GetMinimum(0.0) : h1.GetMinimum());
361  ymax = max(ymax, h1.GetMaximum());
362  }
363  catch(exception&) {}
364 
365  try {
366 
367  TGraph& g1 = dynamic_cast<TGraph&>(*object);
368 
369  for (Int_t i = 0; i != g1.GetN(); ++i) {
370 
371  xmin = min(xmin, g1.GetX()[i]);
372  xmax = max(xmax, g1.GetX()[i]);
373 
374  if (!logy || g1.GetY()[i] > 0.0) {
375  ymin = min(ymin, g1.GetY()[i]);
376  ymax = max(ymax, g1.GetY()[i]);
377  }
378  }
379  }
380  catch(exception&) {}
381 
382  try {
383 
384  TGraphErrors& g1 = dynamic_cast<TGraphErrors&>(*object);
385 
386  for (Int_t i = 0; i != g1.GetN(); ++i) {
387  if (!logy || g1.GetY()[i] - g1.GetEY()[i] > 0.0) { ymin = min(ymin, g1.GetY()[i] - g1.GetEY()[i]); }
388  if (!logy || g1.GetY()[i] + g1.GetEY()[i] > 0.0) { ymax = max(ymax, g1.GetY()[i] + g1.GetEY()[i]); }
389  }
390  }
391  catch(exception&) {}
392 
393  try {
394 
395  TMultiGraph& m1 = dynamic_cast<TMultiGraph&>(*object);
396 
397  for (TIter i1(m1.GetListOfGraphs()); TGraph* g1 = dynamic_cast<TGraph*>(i1()); ) {
398 
399  for (Int_t i = 0; i != g1->GetN(); ++i) {
400 
401  xmin = min(xmin, g1->GetX()[i]);
402  xmax = max(xmax, g1->GetX()[i]);
403 
404  if (!logy || g1->GetY()[i] > 0.0) {
405  ymin = min(ymin, g1->GetY()[i]);
406  ymax = max(ymax, g1->GetY()[i]);
407  }
408  }
409  }
410  }
411  catch(exception&) {}
412 
413  try {
414 
415  TF2& f2 = dynamic_cast<TF2&>(*object);
416  TF1* f1 = NULL;
417 
418  TString formula = f2.GetExpFormula();
419  TString _z_ = TString::Format("%f", 0.5 * (Z.getLowerLimit() + Z.getUpperLimit()));
420 
421  double __xmin;
422  double __xmax;
423  double __ymin;
424  double __ymax;
425 
426  f2.GetRange(__xmin, __ymin, __xmax, __ymax);
427 
428  if (px) {
429 
430  formula.ReplaceAll("y", _z_);
431 
432  f1 = new TF1(MAKE_CSTRING(f2.GetName() << "_px" << LABEL_TERMINATOR << listOfObjects.size()), formula);
433 
434  f1->SetRange(__xmin, __xmax);
435 
436  } else if (py) {
437 
438  formula.ReplaceAll("x", _z_);
439  formula.ReplaceAll("y", "x");
440 
441  f1 = new TF1(MAKE_CSTRING(f2.GetName() << "_py" << LABEL_TERMINATOR << listOfObjects.size()), formula);
442 
443  f1->SetRange(__ymin, __ymax);
444 
445  } else {
446 
447  ERROR("For 2D functions, use option option -P for projections or use JPlot2D" << endl);
448 
449  continue;
450  }
451 
452  DEBUG("TF1: " << f1->GetExpFormula() << endl);
453 
454  f1->SetParameters(f2.GetParameters());
455 
456  object = f1;
457  }
458  catch(exception&) {}
459 
460  try {
461 
462  TF1& f1 = dynamic_cast<TF1&>(*object);
463 
464  double __xmin;
465  double __xmax;
466 
467  f1.GetRange(__xmin, __xmax);
468 
469  xmin = min(xmin, __xmin);
470  xmax = max(xmax, __xmax);
471  ymin = min(ymin, f1.GetMinimum());
472  ymax = max(ymax, f1.GetMaximum());
473  }
474  catch(exception&) {}
475 
476  try {
477 
478  THStack& hs = dynamic_cast<THStack&>(*object);
479 
480  NOTICE("THStack" << endl);
481 
482  TIterator* iterator = hs.GetHists()->MakeIterator();
483 
484  for (size_t index = 1; TObject* i = iterator->Next(); ++index) {
485 
486  TH1& h1 = dynamic_cast<TH1&>(*i);
487 
488  NOTICE("TH1[" << index << "] " << h1.GetName() << endl);
489 
490  xmin = min(xmin, h1.GetXaxis()->GetXmin());
491  xmax = max(xmax, h1.GetXaxis()->GetXmax());
492 
493  ymin = min(ymin, logy ? h1.GetMinimum(0.0) : h1.GetMinimum());
494  ymax = max(ymax, h1.GetMaximum());
495 
496  h1.SetLineWidth(1);
497  h1.SetLineColor(kBlack);
498 
499  h1.SetFillColor(JMarkerAttributes::getInstance().get(index).GetMarkerColor());
500  }
501  }
502  catch(exception&) {}
503 
504  try {
505 
506  TLine& h1 = dynamic_cast<TLine&>(*object);
507 
508  xmin = min(xmin, h1.GetX1());
509  xmax = max(xmax, h1.GetX2());
510  ymin = min(ymin, h1.GetY1());
511  ymax = max(ymax, h1.GetY2());
512  }
513  catch(exception&) {}
514 
515  // label
516 
517  for (TString buffer[] = { object.getLabel(), input->getFilename().c_str(), "" }, *i = buffer; *i != ""; ++i) {
518 
519  *i = (*i)(TRegexp("\\[.*\\]"));
520 
521  DEBUG("Label: <" << *i << ">" << endl);
522 
523  if (i->Length() > 2) {
524  object.setLabel((*i)(1, i->Length() - 2));
525  }
526  }
527 
528  DEBUG("Add object: " << tag << " with label " << object.getLabel() << endl);
529 
530  if (master == NULL) {
531  master = dynamic_cast<TH1*>(object.get());
532  }
533 
534  listOfObjects.push_back(object);
535  }
536  }
537  }
538 
539  if (listOfObjects.empty()) {
540  ERROR("Nothing to draw." << endl);
541  }
542 
543  for (vector<JRootObject>::iterator i = listOfObjects.begin(); i != listOfObjects.end(); ++i) {
544 
545  // set line attributes of associated functions
546  // for single objects, change colour if same range else change style
547  // for multiple objecs, keep colour and change style
548 
549  TH1* h1 = dynamic_cast<TH1*> (i->get());
550  TGraph* g1 = dynamic_cast<TGraph*> (i->get());
551  TMultiGraph* m1 = dynamic_cast<TMultiGraph*>(i->get());
552  TAttLine* ls = dynamic_cast<TAttLine*> (i->get());
553 
554  TList list;
555  TIterator* iterator = NULL;
556 
557  if (h1 != NULL) {
558 
559  iterator = h1->GetListOfFunctions()->MakeIterator();
560 
561  } else if (g1 != NULL) {
562 
563  iterator = g1->GetListOfFunctions()->MakeIterator();
564 
565  } else if (m1 != NULL) {
566 
567  for (TIter i1(m1->GetListOfGraphs()); TGraph* gi = dynamic_cast<TGraph*>(i1()); ) {
568  for (TIter i2(gi->GetListOfFunctions()); TF1* fi = dynamic_cast<TF1*>(i2()); ) {
569  list.Add(fi);
570  }
571  }
572 
573  iterator = list.MakeIterator();
574  }
575 
576  if (iterator != NULL) {
577 
578  Double_t x1[] = { numeric_limits<Double_t>::max(), numeric_limits<Double_t>::max() };
579  Double_t x2[] = { numeric_limits<Double_t>::lowest(), numeric_limits<Double_t>::lowest() };
580 
581  for (int nf = 0, ns = 0, nc = 1; TF1* f1 = (TF1*) iterator->Next(); ++nf) {
582 
583  f1->GetRange(x1[1], x2[1]);
584  f1->SetNpx(1000);
585 
586  dynamic_cast<TAttLine&>(*f1) = JLineAttributes::getInstance().get(0);
587 
588  if (listOfObjects.size() == 1) {
589 
590  if (x1[0] == x1[1] &&
591  x2[0] == x2[1])
592  ++nc; // change colour
593  else if (nf != 0)
594  ++ns; // change style
595 
596  f1->SetLineStyle(JLineAttributes ::getInstance().get(ns).GetLineStyle());
597  f1->SetLineColor(JMarkerAttributes::getInstance().get(nc).GetMarkerColor());
598 
599  } else {
600 
601  // keep colour of base object and accordingly modify line style.
602 
603  f1->SetLineColor(ls->GetLineColor());
604  f1->SetLineStyle(JLineAttributes::getInstance().get(ns++).GetLineStyle());
605  }
606 
607  x1[0] = x1[1];
608  x2[0] = x2[1];
609 
610  // set limits
611 
612  double __xmin;
613  double __xmax;
614 
615  f1->GetRange(__xmin, __xmax);
616 
617  ymin = min(ymin, f1->GetMinimum(__xmin, __xmax));
618  ymax = max(ymax, f1->GetMaximum(__xmin, __xmax));
619  }
620  }
621  }
622 
623  // plot frame
624 
625  if (X != JRange_t()) {
626  xmin = X.getLowerLimit();
627  xmax = X.getUpperLimit();
628  }
629 
630  if (Y != JRange_t()) {
631  ymin = Y.getLowerLimit();
632  ymax = Y.getUpperLimit();
633  } else {
634  setRange(ymin, ymax, logy);
635  }
636 
637  cv->cd(1);
638 
639  if (!listOfObjects.empty()) {
640 
641  if (master == NULL || master->GetXaxis()->GetXmin() > xmin || master->GetXaxis()->GetXmax() < xmax) {
642 
643  if (X != JRange_t()) {
644 
645  master = new TH1D("__H__", NULL, 100, X.getLowerLimit(), X.getUpperLimit());
646 
647  master->SetStats(kFALSE);
648 
649  } else if (xmin < xmax) {
650 
651  master = new TH1D("__H__", NULL, 100, xmin, xmax);
652 
653  master->SetStats(kFALSE);
654  }
655  }
656  }
657 
658  if (master == NULL) {
659 
660  TText* p = new TText(0.5, 0.5, MAKE_CSTRING("No data"));
661 
662  p->SetTextAlign(21);
663  p->SetTextAngle(45);
664  p->Draw();
665 
666  } else {
667 
668  if (logx) { gPad->SetLogx(); }
669  if (logy) { gPad->SetLogy(); }
670 
671  master->GetXaxis()->SetRangeUser(xmin, xmax);
672 
673  if (logx > 1) {
674  setLogarithmicX(master);
675  }
676 
677  master->SetTitle(title.c_str());
678 
679  master->SetMinimum(ymin);
680  master->SetMaximum(ymax);
681 
682  if (xLabel != "") { master->GetXaxis()->SetTitle(xLabel.c_str()); master->GetXaxis()->CenterTitle(true); }
683  if (yLabel != "") { master->GetYaxis()->SetTitle(yLabel.c_str()); master->GetYaxis()->CenterTitle(true); }
684 
685  master->GetXaxis()->SetMoreLogLabels((logx == 1 && log10(xmax/xmin) < 2) ||
686  (logx > 1 && (xmax-xmin) < 2));
687 
688  master->GetYaxis()->SetMoreLogLabels( logy && log10(ymax/ymin) < 2);
689  master->GetYaxis()->SetNoExponent ( logy && log10(ymax/ymin) < 2);
690 
691  if (Ndivisions.first != "") {
692  master->SetNdivisions(Ndivisions.second, Ndivisions.first.c_str());
693  }
694 
695  if (xTimeFormat != "") {
696 
697  master->GetXaxis()->SetTimeDisplay(1);
698 
699  if (xTimeFormat == "utc") {
700  master->GetXaxis()->SetTimeFormat("#splitline{}{#splitline{%d-%m-%y}{ %H:%M}}");
701  master->GetXaxis()->SetTimeOffset(0.0, "gmt");
702  } else if (xTimeFormat == "UTC") {
703  master->GetXaxis()->SetTimeFormat("%d-%m-%y");
704  master->GetXaxis()->SetTimeOffset(0.0, "gmt");
705  } else {
706  master->GetXaxis()->SetTimeFormat(xTimeFormat.c_str());
707  }
708  }
709 
710  master->Draw(option.c_str());
711  }
712 
713  if (logx > 1) {
714  for (vector<JRootObject>::iterator i = listOfObjects.begin(); i != listOfObjects.end(); ++i) {
715  if (dynamic_cast<TH1*>(i->get()) != master) {
716  setLogarithmicX(dynamic_cast<TH1*> (i->get()));
717  setLogarithmicX(dynamic_cast<TGraph*> (i->get()));
718  setLogarithmicX(dynamic_cast<TMultiGraph*>(i->get()));
719  setLogarithmicX(dynamic_cast<TF1*> (i->get()));
720  }
721  }
722  }
723 
724  if (grid.count('x') || grid.count('X')) { gPad->SetGridx(); }
725  if (grid.count('y') || grid.count('Y')) { gPad->SetGridy(); }
726 
727  if (stats != -1)
728  gStyle->SetOptStat(stats);
729  else
730  gStyle->SetOptFit(kFALSE);
731 
732 
733  for (vector<JRootObject>::const_iterator i = listOfObjects.begin(); i != listOfObjects.end(); ++i) {
734 
735  DEBUG("Draw " << (*i)->GetName() << ' ' << (*i)->GetTitle() << endl);
736 
737  string buffer(option);
738 
739  //if (!dynamic_cast<THStack*>(i->get())) {
740  // buffer += "SAMES";
741  //}
742 
743  buffer += "SAME";
744 
745  TF1* f1 = dynamic_cast<TF1*> (i->get());
746  TGraph* g1 = dynamic_cast<TGraph*> (i->get());
747  TMultiGraph* q1 = dynamic_cast<TMultiGraph*>(i->get());
748 
749  if (f1 != NULL) {
750  f1->SetNpx(1000);
751  }
752 
753  if (g1 != NULL) {
754  if (g1->GetN() > 1 && drawLine)
755  buffer += "L"; // drawing cut line
756  else
757  buffer += "P"; // drawing point(s)
758  }
759 
760  if (q1 != NULL) {
761 
762  for (TIter i1(q1->GetListOfGraphs()); TGraph* gi = dynamic_cast<TGraph*>(i1()); ) {
763 
764  string zbuf = buffer;
765 
766  if (gi->GetN() > 1 && drawLine)
767  zbuf += "L"; // drawing cut line
768  else
769  zbuf += "P"; // drawing point(s)
770 
771  gi->Draw(zbuf.c_str());
772  }
773 
774  } else {
775 
776  (*i).Draw(buffer.c_str());
777  }
778  }
779 
780  //gPad->RedrawAxis();
781 
782  if (legend != "") {
783 
784  Ssiz_t height = listOfObjects.size();
785  Ssiz_t width = 1;
786 
787  for (vector<JRootObject>::const_iterator i = listOfObjects.begin(); i != listOfObjects.end(); ++i) {
788  width = max(width, i->getLabel().Length());
789  }
790 
791  TLegend* lg = getLegend(width, height, legend);
792 
793  for (vector<JRootObject>::const_iterator i = listOfObjects.begin(); i != listOfObjects.end(); ++i) {
794  lg->AddEntry(*i, " " + i->getLabel(), isTAttLine(*i) ? "L" : "LPF");
795  }
796 
797  lg->Draw();
798  }
799 
800  cv->Update();
801 
802  if (outputFile != "") {
803  cv->SaveAs(outputFile.c_str());
804  }
805 
806  if (!batch) {
807  tp->Run();
808  }
809 
810  return (master != NULL ? 0 : 1);
811 }
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1500
master
then echo Test string reversed by master(hit< return > to continue)." $DIR/JProcess -c "$DIR/JEcho" -rC fi if (( 1 ))
main
int main(int argc, char *argv[])
Definition: Main.cc:15
JGIZMO::setLogarithmicX
void setLogarithmicX(TF1 *f1)
Make parameter x of function logarithmic (e.g. after filling with log10()).
Definition: JGizmoToolkit.hh:517
JDETECTOR::getLabel
std::string getLabel(const JLocation &location)
Get module label for monitoring and other applications.
Definition: JLocation.hh:246
TObject
Definition: JRoot.hh:19
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition: JPrint.hh:151
h1
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:71
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition: JParser.hh:66
Y
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
Definition: detector-XY:fit1d.sh:33
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JGIZMO::isTAttLine
bool isTAttLine(const TObject *object)
Get drawing option of TH1.
Definition: JGizmoToolkit.hh:184
JPrint.hh
I/O formatting auxiliaries.
JLANG::getInstance
T & getInstance(const T &object)
Get static instance from temporary object.
Definition: JObject.hh:75
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1961
NOTICE
#define NOTICE(A)
Definition: JMessage.hh:64
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
X
then break fi done getCenter read X Y Z let X
Definition: acoustics-fit.sh:387
debug
int debug
debug level
Definition: JSirene.cc:63
JMessage.hh
General purpose messaging.
JGIZMO::LABEL_TERMINATOR
static const char LABEL_TERMINATOR
label terminator
Definition: JRootObject.hh:23
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
JGIZMO::setRange
void setRange(double &xmin, double &xmax, const bool logx)
Set axis range.
Definition: JGizmoToolkit.hh:942
JRange.hh
Auxiliary class to define a range between two values.
JParser.hh
Utility class to parse command line options.
JROOT::getLegend
TLegend * getLegend(const Int_t width, const Int_t height, const std::string option="TR")
Get legend.
Definition: JLegend.hh:28
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
JPARSER::JCounter
Auxiliary class to handle multiple boolean-like I/O.
Definition: JParser.hh:221
JGIZMO::getDirectory
TDirectory * getDirectory(const JRootObjectID &id)
Get TDirectory pointer.
Definition: JGizmoToolkit.hh:131
Z
do set_variable MODULE getModule a $WORKDIR detector_a datx L $STRING JEditDetector a $WORKDIR detector_a datx M $MODULE setz o $WORKDIR detector_a datx JEditDetector a $WORKDIR detector_b datx M $MODULE setz o $WORKDIR detector_b datx done echo Output stored at $WORKDIR detector_a datx and $WORKDIR tripod_a txt JDrawDetector2D a $WORKDIR detector_a datx a $WORKDIR detector_b datx L BL o detector $FORMAT $BATCH JDrawDetector2D T $WORKDIR tripod_a txt T $WORKDIR tripod_b txt L BL o tripod $FORMAT $BATCH JCompareDetector a $WORKDIR detector_a datx b $WORKDIR detector_b datx o $WORKDIR abc root &dev null for KEY in X Y Z
Definition: JAlignDetector.sh:81
g1
Double_t g1(const Double_t x)
Function.
Definition: JQuantiles.cc:25
Generated by   doxygen 1.8.5