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