257 JParser<> zap(
"General purpose plot program for 2D ROOT objects.");
259 zap[
'f'] =
make_field(inputFile,
"<input file>:<object name>");
264 zap[
'x'] =
make_field(X,
"x-abscissa range") = JRange_t::DEFAULT_RANGE();
265 zap[
'y'] =
make_field(Y,
"y-abscissa range") = JRange_t::DEFAULT_RANGE();
266 zap[
'z'] =
make_field(Z,
"ordinate range") = JRange_t::DEFAULT_RANGE();
267 zap[
'X'] =
make_field(logx,
"logarithmic x-axis (-XX log10 axis)");
270 zap[
'P'] =
make_field(project,
"projection") =
"",
"xy",
"yx",
"xz",
"zx",
"yz",
"zy";
279 <<
"\"" <<
JName_t <<
"\" -> ROOT name; "
280 <<
"\"" <<
JTitle_t <<
"\" -> ROOT title)") =
"KM3NeT preliminary";
287 catch(
const exception &error) {
298 ((
TRootCanvas *) cv->GetCanvasImp())->
Connect(
"CloseWindow()",
"TApplication", tp,
"Terminate()");
307 gROOT->SetStyle(
"gplot");
311 cv->SetFillStyle(4000);
316 JMarkerAttributes::getInstance().setMarkerSize(
markerSize);
337 ERROR(
"File: " << input->getFullFilename() <<
" not opened." <<
endl);
347 const TString tag(key->GetName());
353 if (tag.Contains(
regexp) && isTObject(key)) {
356 title = key->GetName();
358 title = key->GetTitle();
365 TH3&
h3 =
dynamic_cast<TH3&
>(*object);
367 object =
h3.Project3D(project.c_str());
373 TH2& h2 =
dynamic_cast<TH2&
>(*object);
375 h2.SetStats(
stats != -1);
377 xmin = min(xmin, h2.GetXaxis()->GetXmin());
378 xmax = max(xmax, h2.GetXaxis()->GetXmax());
379 ymin = min(
ymin, h2.GetYaxis()->GetXmin());
380 ymax = max(
ymax, h2.GetYaxis()->GetXmax());
381 zmin = min(zmin,
logz ? h2.GetMinimum(0.0) : h2.GetMinimum());
382 zmax = max(zmax, h2.GetMaximum());
388 TGraph&
g1 =
dynamic_cast<TGraph&
>(*object);
390 for (
Int_t i = 0; i !=
g1.GetN(); ++i) {
400 if (
dynamic_cast<TGraph*
>(i->get()) !=
NULL) {
405 static_cast<TAttMarker&
>(
g1) = JMarkerAttributes::getInstance().get(
ng);
411 TGraph2D& g2 =
dynamic_cast<TGraph2D&
>(*object);
413 for (
Int_t i = 0; i != g2.GetN(); ++i) {
414 if (!
logz || g2.GetZ()[i] > 0.0) {
415 xmin = min(xmin, g2.GetX()[i]);
416 xmax = max(xmax, g2.GetX()[i]);
419 zmin = min(zmin, g2.GetZ()[i]);
420 zmax = max(zmax, g2.GetZ()[i]);
425 g2.SetMinimum(Z.getLowerLimit());
426 g2.SetMaximum(Z.getUpperLimit());
433 TF2& f2 =
dynamic_cast<TF2&
>(*object);
435 f2.SetLineColor(
kRed);
449 zmin = min(zmin, f2.GetMinimum());
450 zmax = max(zmax, f2.GetMaximum());
456 for (TString buffer[] = {
object.getLabel(), input->getFilename().c_str(),
"" }, *i = buffer; *i !=
""; ++i) {
458 *i = (*i)(
TRegexp(
"\\[.*\\]"));
460 if ((*i).Length() > 2) {
461 object.setLabel((*i)(1, (*i).Length() - 2));
465 if (
dynamic_cast<TH2*
> (
object.get()) !=
NULL ||
466 dynamic_cast<TGraph*
> (
object.get()) !=
NULL ||
467 dynamic_cast<TGraph2D*
>(
object.get()) !=
NULL ||
469 dynamic_cast<TLine*
> (
object.get()) !=
NULL ||
470 dynamic_cast<TText*
> (
object.get()) !=
NULL ||
471 dynamic_cast<TF2*
> (
object.get()) !=
NULL) {
473 DEBUG(
"Add object: " << tag <<
" with label <" <<
object.getLabel() <<
">" <<
endl);
483 ERROR(
"For other objects than 2D histograms, use JPlot1D" <<
endl);
497 if (option.find(
"COLZ") != string::npos ||
498 option.find(
"colz") != string::npos) {
499 gPad->SetRightMargin(0.20);
503 xmin = X.getLowerLimit();
504 xmax = X.getUpperLimit();
508 ymin = Y.getLowerLimit();
509 ymax = Y.getUpperLimit();
513 zmin = Z.getLowerLimit();
514 zmax = Z.getUpperLimit();
515 }
else if (zmax > zmin) {
516 setRange(zmin, zmax,
logz);
541 if (logx) {
gPad->SetLogx(); }
545 master->SetTitle(title.c_str());
547 master->GetXaxis()->SetRangeUser(xmin, xmax);
550 if (logx > 1) { setLogarithmicX(
master); }
553 if (logx > 2) {
master->GetXaxis()->SetNoExponent(); }
554 if (
logy > 2) {
master->GetYaxis()->SetNoExponent(); }
563 master->GetXaxis()->SetMoreLogLabels((logx == 1 &&
log10(xmax/xmin) < 2) ||
564 (logx > 1 && xmax-xmin < 2));
569 master->SetNdivisions(i->second, i->first.c_str());
574 master->Draw(option.c_str());
577 if (logx > 1 ||
logy > 1) {
579 if (
dynamic_cast<TH2*
> (i->get()) !=
master &&
580 dynamic_cast<TGraph2D*
>(i->get()) !=
master) {
605 if (
grid.count(
'x') ||
grid.count(
'X')) {
gPad->SetGridx(); }
606 if (
grid.count(
'y') ||
grid.count(
'Y')) {
gPad->SetGridy(); }
618 DEBUG(
"Draw " << (*i)->GetName() <<
' ' << (*i)->GetTitle() <<
endl);
620 string buffer(option);
626 TH2& h2 =
dynamic_cast<TH2&
>(*(*i));
635 TGraph&
g1 =
dynamic_cast<TGraph&
>(*(*i));
643 TGraph2D& g2 =
dynamic_cast<TGraph2D&
>(*(*i));
652 TF2& f2 =
dynamic_cast<TF2&
>(*(*i));
663 (*i)->Draw(buffer.c_str());