/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/Validation/RecoTau/plugins/DQMHistPlotter.cc

Go to the documentation of this file.

#include "Validation/RecoTau/plugins/DQMHistPlotter.h"

#include "Validation/RecoTau/plugins/dqmAuxFunctions.h"

// framework & common header files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

//DQM services
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include <TCanvas.h>
#include <TPad.h>
#include <TPostScript.h>
#include <TStyle.h>
#include <TROOT.h>
#include <TMath.h>

#include <iostream>

//defaults for cfgEntryProcess
const double defaultScaleFactor = 1.;
const std::string type_smMC = "smMC";
const std::string type_bsmMC = "bsmMC";
const std::string type_smSumMC = "smSumMC";
const std::string type_Data = "Data";

//defaults for cfgEntryAxisX
const double defaultMinX = -1.;
const double defaultMaxX = -1.;
const double defaultXaxisTitleOffset = 1.0;
const double defaultXaxisTitleSize = 0.05;

//defaults for cfgEntryAxisY
const double defaultMinY_linear = 0.;
const double defaultMinY_log = 1.e-2;
const double defaultMaxY_linear = -1.;
const double defaultMaxY_log = -1.;
const std::string yScale_linear = "linear";
const std::string yScale_log = "log";
const std::string defaultYscale = yScale_linear;
const double defaultYaxisTitleOffset = 1.0;
const double defaultYaxisTitleSize = 0.05;
const double defaultYaxisMaximumScaleFactor_linear = 1.6;
const double defaultYaxisMaximumScaleFactor_log = 5.e+2;

double TauDQMHistPlotter::cfgEntryAxisY::yAxisNorm_ = 0.;

// defaults for cfgEntryLegend
const double defaultLegendPosX = 0.50;
const double defaultLegendPosY = 0.55;
const double defaultLegendSizeX =  0.39;
const double defaultLegendSizeY =  0.34;
const std::string defaultLegendHeader = "";
const std::string defaultLegendOptions = "brNDC";
const int defaultLegendBorderSize = 0;
const int defaultLegendFillColor = 0;

// defaults for cfgEntryLabel
const double defaultLabelPosX = 0.66;
const double defaultLabelPosY = 0.82;
const double defaultLabelSizeX = 0.26;
const double defaultLabelSizeY = 0.10;
const std::string defaultLabelOptions = "brNDC";
const int defaultLabelBorderSize = 0;
const int defaultLabelFillColor = 0;
const int defaultLabelTextColor = 1;
const double defaultLabelTextSize = 0.05;
const int defaultLabelTextAlign = 22; // horizontally and vertically centered, see documentation of TAttText
const double defaultLabelTextAngle = 0.;

// defaults for cfgEntryDrawOption
const int defaultMarkerColor = 1;
const int defaultMarkerSize = 1;
const int defaultMarkerStyle = 2;
const int defaultLineColor = 0;
const int defaultLineStyle = 1;
const int defaultLineWidth = 2;
const int defaultFillColor = 0;
const int defaultFillStyle = 1001;
const std::string defaultDrawOption = "";
const std::string defaultDrawOptionLegend = "lpf";

const std::string drawOption_eBand = "eBand";

// defaults for cfgEntryDrawJob
const bool defaultDoOverlayPlots = false;

// global defaults
const int defaultCanvasSizeX = 800;
const int defaultCanvasSizeY = 600;

const std::string drawOptionSeparator = "#.#";

const int verbosity = 0;

template <class T>
void checkCfgDef(const std::string& cfgEntryName, std::map<std::string, T>& def, int& errorFlag,
                 const std::string& defType, const std::string& drawJobName)
{
  if ( def.find(cfgEntryName) == def.end() ) {
    edm::LogError ("checkCfgDef") << " " << defType << " = " << cfgEntryName << " undefined, needed by drawJob = " << drawJobName << " !!";
    errorFlag = 1;
  }
}

template <class T>
void checkCfgDefs(const std::vector<std::string>& cfgEntryNames, std::map<std::string, T>& def, int& errorFlag,
                  const std::string& defType, const std::string& drawJobName)
{
  for ( std::vector<std::string>::const_iterator cfgEntryName = cfgEntryNames.begin();
        cfgEntryName != cfgEntryNames.end(); ++cfgEntryName ) {
    checkCfgDef(*cfgEntryName, def, errorFlag, defType, drawJobName);
  }
}

template <class T>
const T* findCfgDef(const std::string& cfgEntryName, std::map<std::string, T>& def,
                    const std::string& defType, const std::string& drawJobName)
{
  typename std::map<std::string, T>::const_iterator it = def.find(cfgEntryName);
  if ( it != def.end() ) {
    return &(it->second);
  } else {
    edm::LogError ("findCfgDef") << " " << defType << " = " << cfgEntryName << " undefined, needed by drawJob = " << drawJobName << " !!";
    return NULL;
  } 
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

typedef std::pair<TH1*, std::string> histoDrawEntry;

void drawHistograms(const std::list<histoDrawEntry>& histogramList, bool& isFirstHistogram)
{
  for ( std::list<histoDrawEntry>::const_iterator it = histogramList.begin();
        it != histogramList.end(); ++it ) {
    std::string drawOption = ( isFirstHistogram ) ? it->second : std::string(it->second).append("same");
    it->first->Draw(drawOption.data());
    isFirstHistogram = false;
  }
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

bool find_vstring(const std::vector<std::string>& vs, const std::string& s)
{
  for ( std::vector<std::string>::const_iterator it = vs.begin();
        it != vs.end(); ++it ) {
    if ( (*it) == s ) return true;
  }
  return false;
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::cfgEntryProcess::cfgEntryProcess(const std::string& name, const edm::ParameterSet& cfg)
{
  name_ = name;

  dqmDirectory_ = cfg.getParameter<std::string>("dqmDirectory");
  
  legendEntry_ = cfg.getParameter<std::string>("legendEntry");
  legendEntryErrorBand_ = ( cfg.exists("legendEntryErrorBand") ) ? 
    cfg.getParameter<std::string>("legendEntryErrorBand") : std::string(legendEntry_).append(" Uncertainty");

  type_ = cfg.getParameter<std::string>("type");
  
  if ( verbosity ) print();
}

void TauDQMHistPlotter::cfgEntryProcess::print() const
{
  std::cout << "<TauDQMHistPlotter::cfgEntryProcess::print>:" << std::endl;
  std::cout << " name = " << name_ << std::endl;
  std::cout << " dqmDirectory = " << dqmDirectory_ << std::endl;
  std::cout << " legendEntry = " << legendEntry_ << std::endl;
  std::cout << " legendEntryErrorBand = " << legendEntryErrorBand_ << std::endl;
  std::cout << " type = " << type_ << std::endl;
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::cfgEntryAxisX::cfgEntryAxisX(const std::string& name, const edm::ParameterSet& cfg)
{
  name_ = name;

  minX_ = ( cfg.exists("minX") ) ? cfg.getParameter<double>("minX") : defaultMinX;
  maxX_ = ( cfg.exists("maxX") ) ? cfg.getParameter<double>("maxX") : defaultMaxX;
  xAxisTitle_ = cfg.getParameter<std::string>("xAxisTitle");
  xAxisTitleOffset_ = ( cfg.exists("xAxisTitleOffset") ) ? cfg.getParameter<double>("xAxisTitleOffset") : defaultXaxisTitleOffset;
  xAxisTitleSize_ = ( cfg.exists("xAxisTitleSize") ) ? cfg.getParameter<double>("xAxisTitleSize") : defaultXaxisTitleSize;

  if ( verbosity ) print();
}

void TauDQMHistPlotter::cfgEntryAxisX::print() const
{
  std::cout << "<TauDQMHistPlotter::cfgEntryAxisX::print>:" << std::endl;
  std::cout << " name = " << name_ << std::endl;
  std::cout << " minX_ = " << minX_ << std::endl;
  std::cout << " maxX_ = " << maxX_ << std::endl;
  std::cout << " xAxisTitle = " << xAxisTitle_ << std::endl;
  std::cout << " xAxisTitleOffset = " << xAxisTitleOffset_ << std::endl;
  std::cout << " xAxisTitleSize = " << xAxisTitleSize_ << std::endl;
}

void TauDQMHistPlotter::cfgEntryAxisX::applyTo(TH1* histogram) const
{
  if ( histogram ) {
    double xMin = ( minX_ != defaultMinX ) ? minX_ : histogram->GetXaxis()->GetXmin();
    double xMax = ( maxX_ != defaultMaxX ) ? maxX_ : histogram->GetXaxis()->GetXmax();
    histogram->SetAxisRange(xMin, xMax, "X");
    histogram->GetXaxis()->SetTitle(xAxisTitle_.data());
    histogram->GetXaxis()->SetTitleOffset(xAxisTitleOffset_);
    histogram->GetXaxis()->SetTitleSize(xAxisTitleSize_);
  }
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::cfgEntryAxisY::cfgEntryAxisY(const std::string& name, const edm::ParameterSet& cfg)
{
  name_ = name;

  minY_linear_ = ( cfg.exists("minY_linear") ) ? cfg.getParameter<double>("minY_linear") : defaultMinY_linear;
  minY_log_ = ( cfg.exists("minY_log") ) ? cfg.getParameter<double>("minY_log") : defaultMinY_log;
  maxY_linear_ = ( cfg.exists("maxY_linear") ) ? cfg.getParameter<double>("maxY_linear") : defaultMaxY_linear;
  maxY_log_ = ( cfg.exists("maxY_log") ) ? cfg.getParameter<double>("maxY_log") : defaultMaxY_log;
  yScale_ = ( cfg.exists("yScale") ) ? cfg.getParameter<std::string>("yScale") : defaultYscale;
  yAxisTitle_ = cfg.getParameter<std::string>("yAxisTitle");
  yAxisTitleOffset_ = ( cfg.exists("yAxisTitleOffset") ) ? cfg.getParameter<double>("yAxisTitleOffset") : defaultYaxisTitleOffset;
  yAxisTitleSize_ = ( cfg.exists("yAxisTitleSize") ) ? cfg.getParameter<double>("yAxisTitleSize") : defaultYaxisTitleSize;

  if ( verbosity ) print();
}

void TauDQMHistPlotter::cfgEntryAxisY::print() const
{
  std::cout << "<TauDQMHistPlotter::cfgEntryAxisY::print>:" << std::endl;
  std::cout << " name = " << name_ << std::endl;
  std::cout << " minY_linear = " << minY_linear_ << std::endl;
  std::cout << " minY_log = " << minY_log_ << std::endl;
  std::cout << " maxY_linear = " << maxY_linear_ << std::endl;
  std::cout << " maxY_log = " << maxY_log_ << std::endl;
  std::cout << " yScale = " << yScale_ << std::endl;
  std::cout << " yAxisTitle = " << yAxisTitle_ << std::endl;
  std::cout << " yAxisTitleOffset = " << yAxisTitleOffset_ << std::endl;
  std::cout << " yAxisTitleSize = " << yAxisTitleSize_ << std::endl;
}

void TauDQMHistPlotter::cfgEntryAxisY::applyTo(TH1* histogram) const
{
  if ( histogram ) {
    bool yLogScale = ( yScale_ == yScale_log ) ? true : false;
    double minY = ( yLogScale ) ? minY_log_ : minY_linear_;
    histogram->SetMinimum(minY);
    double maxY = ( yLogScale ) ? maxY_log_ : maxY_linear_;
    double defaultMaxY = ( yLogScale ) ? defaultMaxY_log : defaultMaxY_linear;
    if ( maxY != defaultMaxY ) {
//--- normalize y-axis range using given configuration parameter
      histogram->SetMaximum(maxY);
    } else {
//--- in case configuration parameter for y-axis range not explicitely given,
//    normalize y-axis range to maximum of any histogram included in drawJob
//    times defaultYaxisMaximumScaleFactor (apply scale factor in order to make space for legend)
      double defaultYaxisMaximumScaleFactor = ( yLogScale ) ? defaultYaxisMaximumScaleFactor_log : defaultYaxisMaximumScaleFactor_linear;
      histogram->SetMaximum(defaultYaxisMaximumScaleFactor*yAxisNorm_);
    }
    histogram->GetYaxis()->SetTitle(yAxisTitle_.data());
    histogram->GetYaxis()->SetTitleOffset(yAxisTitleOffset_);
    histogram->GetYaxis()->SetTitleSize(yAxisTitleSize_);
  }
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::cfgEntryLegend::cfgEntryLegend(const std::string& name, const edm::ParameterSet& cfg)
{
  name_ = name;

  posX_ = ( cfg.exists("posX") ) ? cfg.getParameter<double>("posX") : defaultLegendPosX;
  posY_ = ( cfg.exists("posY") ) ? cfg.getParameter<double>("posY") : defaultLegendPosY;
  sizeX_ = ( cfg.exists("sizeX") ) ? cfg.getParameter<double>("sizeX") : defaultLegendSizeX;
  sizeY_ = ( cfg.exists("sizeY") ) ? cfg.getParameter<double>("sizeY") : defaultLegendSizeY;
  header_ = ( cfg.exists("header") ) ? cfg.getParameter<std::string>("header") : defaultLegendHeader;
  option_ = ( cfg.exists("option") ) ? cfg.getParameter<std::string>("option") : defaultLegendOptions;
  borderSize_ = ( cfg.exists("borderSize") ) ? cfg.getParameter<int>("borderSize") : defaultLegendBorderSize;
  fillColor_ = ( cfg.exists("fillColor") ) ? cfg.getParameter<int>("fillColor") : defaultLegendFillColor;
 
  if ( verbosity ) print();
}

void TauDQMHistPlotter::cfgEntryLegend::print() const
{
  std::cout << "<TauDQMHistPlotter::cfgEntryLegend::print>:" << std::endl;
  std::cout << " name = " << name_ << std::endl;
  std::cout << " posX = " << posX_ << std::endl;
  std::cout << " posY = " << posY_ << std::endl;
  std::cout << " sizeX = " << sizeX_ << std::endl;
  std::cout << " sizeY = " << sizeY_ << std::endl;
  std::cout << " header = " << header_ << std::endl;
  std::cout << " option = " << option_ << std::endl;
  std::cout << " borderSize = " << borderSize_ << std::endl;
  std::cout << " fillColor = " << fillColor_ << std::endl;
}

void TauDQMHistPlotter::cfgEntryLegend::applyTo(TLegend* legend) const
{
  if ( legend ) {
    legend->SetX1(posX_);
    legend->SetY1(posY_);
    legend->SetX2(posX_ + sizeX_);
    legend->SetY2(posY_ + sizeY_);
    legend->SetHeader(header_.data());
    legend->SetOption(option_.data());
    legend->SetBorderSize(borderSize_);
    legend->SetFillColor(fillColor_);
  }
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::cfgEntryLabel::cfgEntryLabel(const std::string& name, const edm::ParameterSet& cfg)
{
  name_ = name;

  posX_ = ( cfg.exists("posX") ) ? cfg.getParameter<double>("posX") : defaultLabelPosX;
  posY_ = ( cfg.exists("posY") ) ? cfg.getParameter<double>("posY") : defaultLabelPosY;
  sizeX_ = ( cfg.exists("sizeX") ) ? cfg.getParameter<double>("sizeX") : defaultLabelSizeX;
  sizeY_ = ( cfg.exists("sizeY") ) ? cfg.getParameter<double>("sizeY") : defaultLabelSizeY;
  option_ = ( cfg.exists("option") ) ? cfg.getParameter<std::string>("option") : defaultLabelOptions;
  borderSize_ = ( cfg.exists("borderSize") ) ? cfg.getParameter<int>("borderSize") : defaultLabelBorderSize;
  fillColor_ = ( cfg.exists("fillColor") ) ? cfg.getParameter<int>("fillColor") : defaultLabelFillColor;
  textColor_ = ( cfg.exists("textColor") ) ? cfg.getParameter<int>("textColor") : defaultLabelTextColor;
  textSize_ = ( cfg.exists("textSize") ) ? cfg.getParameter<double>("textSize") : defaultLabelTextSize;
  textAlign_ = ( cfg.exists("textAlign") ) ? cfg.getParameter<int>("textAlign") : defaultLabelTextAlign;
  textAngle_ = ( cfg.exists("textAngle") ) ? cfg.getParameter<double>("textAngle") : defaultLabelTextAngle;
  text_ = cfg.getParameter<vstring>("text");
 
  if ( verbosity ) print();
}

void TauDQMHistPlotter::cfgEntryLabel::print() const
{
  std::cout << "<TauDQMHistPlotter::cfgEntryLabel::print>:" << std::endl;
  std::cout << " name = " << name_ << std::endl;
  std::cout << " posX = " << posX_ << std::endl;
  std::cout << " posY = " << posY_ << std::endl;
  std::cout << " sizeX = " << sizeX_ << std::endl;
  std::cout << " sizeY = " << sizeY_ << std::endl;
  std::cout << " option = " << option_ << std::endl;
  std::cout << " borderSize = " << borderSize_ << std::endl;
  std::cout << " fillColor = " << fillColor_ << std::endl;
  std::cout << " textColor = " << textColor_ << std::endl;
  std::cout << " textSize = " << textSize_ << std::endl;
  std::cout << " textAlign = " << textAlign_ << std::endl;
  std::cout << " textAngle = " << textAngle_ << std::endl;
  std::cout << " text = " << format_vstring(text_) << std::endl;
}

void TauDQMHistPlotter::cfgEntryLabel::applyTo(TPaveText* label) const
{
  if ( label ) {
//--- WARNING: need to call TPaveText::SetX1NDC, **not** TPaveText::SetX1 !!
//             (see documentation of base-class constructor
//               TPave::TPave(Double_t, Double_t,Double_t, Double_t, Int_t, Option_t*)
//              in TPave.cxx for details)
    label->SetX1NDC(posX_);
    label->SetY1NDC(posY_);
    label->SetX2NDC(posX_ + sizeX_);
    label->SetY2NDC(posY_ + sizeY_);
    label->SetOption(option_.data());
    label->SetBorderSize(borderSize_);
    label->SetFillColor(fillColor_);
    label->SetTextColor(textColor_);
    label->SetTextSize(textSize_);
    label->SetTextAlign(textAlign_);
    label->SetTextAngle(textAngle_);
    for ( vstring::const_iterator line = text_.begin();
          line != text_.end(); ++line ) {
      label->AddText(line->data());
    }
  }
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::cfgEntryDrawOption::cfgEntryDrawOption(const std::string& name, const edm::ParameterSet& cfg)
{
  name_ = name;

  markerColor_ = ( cfg.exists("markerColor") ) ? cfg.getParameter<int>("markerColor") : defaultMarkerColor;
  markerSize_  = ( cfg.exists("markerSize")  ) ? cfg.getParameter<double>("markerSize")  : defaultMarkerSize;
  markerStyle_ = ( cfg.exists("markerStyle") ) ? cfg.getParameter<int>("markerStyle") : defaultMarkerStyle;
  
  lineColor_ = ( cfg.exists("lineColor") ) ? cfg.getParameter<int>("lineColor") : defaultLineColor;
  lineStyle_ = ( cfg.exists("lineStyle") ) ? cfg.getParameter<int>("lineStyle") : defaultLineStyle;
  lineWidth_ = ( cfg.exists("lineWidth") ) ? cfg.getParameter<int>("lineWidth") : defaultLineWidth;
  
  fillColor_ = ( cfg.exists("fillColor") ) ? cfg.getParameter<int>("fillColor") : defaultFillColor;
  fillStyle_ = ( cfg.exists("fillStyle") ) ? cfg.getParameter<int>("fillStyle") : defaultFillStyle;
  
  drawOption_ = ( cfg.exists("drawOption") ) ? cfg.getParameter<std::string>("drawOption") : defaultDrawOption;
  drawOptionLegend_ = ( cfg.exists("drawOptionLegend") ) ? cfg.getParameter<std::string>("drawOptionLegend") : defaultDrawOptionLegend;
  
  if ( verbosity ) print();
}

TauDQMHistPlotter::cfgEntryDrawOption::cfgEntryDrawOption(const std::string& name, const cfgEntryDrawOption& blueprint)
  : name_(name),
    markerColor_(blueprint.markerColor_), markerSize_(blueprint.markerSize_), markerStyle_(blueprint.markerStyle_),
    lineColor_(blueprint.lineColor_), lineStyle_(blueprint.lineStyle_), lineWidth_(blueprint.lineWidth_),
    fillColor_(blueprint.fillColor_), fillStyle_(blueprint.fillStyle_),
    drawOption_(blueprint.drawOption_), drawOptionLegend_(blueprint.drawOptionLegend_)
{
  if ( verbosity ) print();
}

void TauDQMHistPlotter::cfgEntryDrawOption::print() const
{
  std::cout << "<TauDQMHistPlotter::cfgEntryDrawOption::print>:" << std::endl;
  std::cout << " name = " << name_ << std::endl;
  std::cout << " markerColor = " << markerColor_ << std::endl;
  std::cout << " markerSize = " << markerSize_ << std::endl;
  std::cout << " markerStyle = " << markerStyle_ << std::endl;
  std::cout << " lineColor = " << lineColor_ << std::endl;
  std::cout << " lineStyle = " << lineStyle_ << std::endl;
  std::cout << " lineWidth = " << lineWidth_ << std::endl;
  std::cout << " fillColor = " << fillColor_ << std::endl;
  std::cout << " fillStyle = " << fillStyle_ << std::endl;
  std::cout << " drawOption = " << drawOption_ << std::endl;
  std::cout << " drawOptionLegend = " << drawOptionLegend_ << std::endl;
}

void TauDQMHistPlotter::cfgEntryDrawOption::applyTo(TH1* histogram) const 
{
  if ( histogram ) {
    histogram->SetMarkerColor(markerColor_);
    histogram->SetMarkerSize(markerSize_);
    histogram->SetMarkerStyle(markerStyle_);  
    histogram->SetLineColor(lineColor_);
    histogram->SetLineStyle(lineStyle_);
    histogram->SetLineWidth(lineWidth_);   
    histogram->SetFillColor(fillColor_);
    histogram->SetFillStyle(fillStyle_);
  }
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::plotDefEntry::plotDefEntry(const std::string& dqmMonitorElement, 
                                           const std::string& drawOptionEntry,
                                           const std::string& legendEntry,
                                           const std::string& legendEntryErrorBand,
                                           const std::string& process,
                                           bool doStack)
  : dqmMonitorElement_(dqmMonitorElement), 
    drawOptionEntry_(drawOptionEntry),
    legendEntry_(legendEntry),
    legendEntryErrorBand_(legendEntryErrorBand),
    process_(process), 
    doStack_(doStack),
    isErrorBand_(false)
{
  //if ( verbosity ) print();
}

TauDQMHistPlotter::plotDefEntry::plotDefEntry(const plotDefEntry& blueprint)
  : dqmMonitorElement_(blueprint.dqmMonitorElement_), 
    drawOptionEntry_(blueprint.drawOptionEntry_),
    legendEntry_(blueprint.legendEntry_),
    legendEntryErrorBand_(blueprint.legendEntryErrorBand_),
    process_(blueprint.process_), 
    doStack_(blueprint.doStack_),
    isErrorBand_(false)
{
  //if ( verbosity ) print();
}

void TauDQMHistPlotter::plotDefEntry::print() const
{
  std::cout << "<TauDQMHistPlotter::plotDefEntry::print>:" << std::endl;
  std::cout << " dqmMonitorElement = " << dqmMonitorElement_ << std::endl;
  std::cout << " drawOptionEntry = " << drawOptionEntry_ << std::endl;
  std::cout << " legendEntry = " << legendEntry_ << std::endl;
  std::cout << " legendEntryErrorBand = " << legendEntryErrorBand_ << std::endl;
  std::cout << " process = " << process_ << std::endl;
  std::cout << " doStack = " << doStack_ << std::endl;
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::cfgEntryDrawJob::cfgEntryDrawJob(const std::string& name, 
                                                 const plotDefList& plotDefList, 
                                                 const std::string& title,
                                                 const std::string& xAxis, const std::string& yAxis, 
                                                 const std::string& legend, 
                                                 const vstring& labels)
{
  name_ = name;

  for ( plotDefList::const_iterator it = plotDefList.begin();
        it != plotDefList.end(); ++it ) {
    plots_.push_back(plotDefEntry(*it));
  }
  
  title_ = title;

  xAxis_ = xAxis;
  yAxis_ = yAxis;
  
  legend_ = legend;

  for ( vstring::const_iterator it = labels.begin();
        it != labels.end(); ++it ) {
    labels_.push_back(std::string(*it));
  }

  if ( verbosity ) print();
}

void TauDQMHistPlotter::cfgEntryDrawJob::print() const
{
  std::cout << "<TauDQMHistPlotter::cfgSetDrawJob::print>:" << std::endl;
  std::cout << " name = " << name_ << std::endl;
  std::cout << "plots = {" << std::endl;
  for ( plotDefList::const_iterator plot = plots_.begin();
        plot != plots_.end(); ++plot ) {
    plot->print();
  }
  std::cout << "}" << std::endl;
  std::cout << " title = " << title_ << std::endl;
  std::cout << " xAxis = " << xAxis_ << std::endl;
  std::cout << " yAxis = " << yAxis_ << std::endl;
  std::cout << " legend = " << legend_ << std::endl;
  std::cout << " labels = " << format_vstring(labels_) << std::endl;
}

//
//-----------------------------------------------------------------------------------------------------------------------
//

TauDQMHistPlotter::TauDQMHistPlotter(const edm::ParameterSet& cfg)
{
  if( verbosity)
    std::cout << "<TauDQMHistPlotter::TauDQMHistPlotter>:" << std::endl;

  toFile_ = cfg.getParameter<bool>("PrintToFile");
  cfgError_ = 0;

//--- configure processes  
  //std::cout << "--> configuring processes..." << std::endl;
  edm::ParameterSet cfgParSet_processes = cfg.getParameter<edm::ParameterSet>("processes");
  readCfgParameter<cfgEntryProcess>(cfgParSet_processes, processes_);
  
//--- check that process types are defined
  //std::cout << "--> checking configuration parameters..." << std::endl;

  int numProcesses_Data = 0;
  int numProcesses_sumMC = 0;
  for ( std::map<std::string, cfgEntryProcess>::const_iterator process = processes_.begin();
        process != processes_.end(); ++process ) {
    const std::string& type = process->second.type_;

    if ( !((type == type_smMC) ||
           (type == type_bsmMC) ||
           (type == type_smSumMC) ||
           (type == type_Data)) ) {
      edm::LogError ("TauDQMHistPlotter") << " Undefined process type = " << type << " !!";
      cfgError_ = 1;
    }

    if ( type == type_smSumMC ) ++numProcesses_sumMC;
    if ( type == type_Data ) ++numProcesses_Data;
  }

  if ( (numProcesses_Data > 1) || (numProcesses_sumMC > 1) ) {
    edm::LogError ("TauDQMHistPlotter") << " Cannot have more than one process of types sumMC and Data !!";
    cfgError_ = 1;
  }

//--- configure x-axes  
  //std::cout << "--> configuring x-axes..." << std::endl;
  edm::ParameterSet cfgParSet_xAxes = cfg.getParameter<edm::ParameterSet>("xAxes");
  readCfgParameter<cfgEntryAxisX>(cfgParSet_xAxes, xAxes_);

//--- configure y-axes  
  //std::cout << "--> configuring y-axes..." << std::endl;
  edm::ParameterSet cfgParSet_yAxes = cfg.getParameter<edm::ParameterSet>("yAxes");
  readCfgParameter<cfgEntryAxisY>(cfgParSet_yAxes, yAxes_);

//--- configure legends
  //std::cout << "--> configuring legends..." << std::endl;
  edm::ParameterSet cfgParSet_legends = cfg.getParameter<edm::ParameterSet>("legends");
  readCfgParameter<cfgEntryLegend>(cfgParSet_legends, legends_);

//--- configure labels
  //std::cout << "--> configuring labels..." << std::endl;
  edm::ParameterSet cfgParSet_labels = cfg.getParameter<edm::ParameterSet>("labels");
  readCfgParameter<cfgEntryLabel>(cfgParSet_labels, labels_);

//--- configure drawOptions
  //std::cout << "--> configuring drawOptions..." << std::endl;
  if ( cfg.exists("drawOptionSets") ) {
    edm::ParameterSet drawOptionSets = cfg.getParameter<edm::ParameterSet>("drawOptionSets");
    vstring drawOptionSetNames = drawOptionSets.getParameterNamesForType<edm::ParameterSet>();
    for ( vstring::const_iterator drawOptionSetName = drawOptionSetNames.begin(); 
          drawOptionSetName != drawOptionSetNames.end(); ++drawOptionSetName ) {
      edm::ParameterSet drawOptionSet = drawOptionSets.getParameter<edm::ParameterSet>(*drawOptionSetName);

      vstring drawOptionEntryNames = drawOptionSet.getParameterNamesForType<edm::ParameterSet>();
      for ( vstring::const_iterator drawOptionEntryName = drawOptionEntryNames.begin(); 
            drawOptionEntryName != drawOptionEntryNames.end(); ++drawOptionEntryName ) {
        edm::ParameterSet drawOptionEntry = drawOptionSet.getParameter<edm::ParameterSet>(*drawOptionEntryName);

        std::string drawOptionEntryName_full = std::string(*drawOptionSetName).append(drawOptionSeparator).append(*drawOptionEntryName);
        drawOptionEntries_.insert(std::pair<std::string, cfgEntryDrawOption>
                                  (drawOptionEntryName_full, cfgEntryDrawOption(drawOptionEntryName_full, drawOptionEntry)));
      }
    }
  }

  if ( cfg.exists("drawOptionEntries") ) {
    edm::ParameterSet cfgParSet_drawOptionEntries = cfg.getParameter<edm::ParameterSet>("drawOptionEntries");
    readCfgParameter<cfgEntryDrawOption>(cfgParSet_drawOptionEntries, drawOptionEntries_);
  }

//--- configure drawJobs
  //std::cout << "--> configuring drawJobs..." << std::endl;
  edm::ParameterSet drawJobs = cfg.getParameter<edm::ParameterSet>("drawJobs");
  vstring drawJobNames = drawJobs.getParameterNamesForType<edm::ParameterSet>();
  for ( vstring::const_iterator drawJobName = drawJobNames.begin(); 
        drawJobName != drawJobNames.end(); ++drawJobName ) {
    edm::ParameterSet drawJob = drawJobs.getParameter<edm::ParameterSet>(*drawJobName);

    std::map<int, plotDefList> plotDefMap;

    if ( drawJob.existsAs<edm::ParameterSet>("plots") ) { // display same monitor element for different processes
      edm::ParameterSet plots = drawJob.getParameter<edm::ParameterSet>("plots");

      vstring dqmMonitorElements = plots.getParameter<vstring>("dqmMonitorElements");
      vstring processes = plots.getParameter<vstring>("processes");

      std::string drawOptionSet = drawJob.getParameter<std::string>("drawOptionSet");
      //std::cout << "drawOptionSet = " << drawOptionSet << std::endl;

      vstring stack = ( cfg.exists("stack") ) ? drawJob.getParameter<vstring>("stack") : vstring();

      for ( vstring::const_iterator process = processes.begin();
            process != processes.end(); ++process ) {
        int index = 0;
        for ( vstring::const_iterator dqmMonitorElement = dqmMonitorElements.begin();
              dqmMonitorElement != dqmMonitorElements.end(); ++dqmMonitorElement ) {
          bool stack_dqmMonitorElement = find_vstring(stack, *process);
          std::string drawOptionEntry = std::string(drawOptionSet).append(drawOptionSeparator).append(*process);
          plotDefMap[index].push_back(plotDefEntry(*dqmMonitorElement, drawOptionEntry, "", "", *process, stack_dqmMonitorElement));
          ++index;
        }
      }
    } else { // display different monitor elements for same process
      typedef std::vector<edm::ParameterSet> vParameterSet;
      vParameterSet plots = drawJob.getParameter<vParameterSet>("plots");

      std::string process = ( drawJob.exists("process") ) ? drawJob.getParameter<std::string>("process") : "";
      //std::cout << "process (globally set) = " << process << std::endl;

      for ( vParameterSet::const_iterator plot = plots.begin(); 
            plot != plots.end(); ++plot ) {

        if ( process == "" || plot->exists("process")) {
          process = plot->getParameter<std::string>("process");
          //std::cout << "process (locally set) = " << process << std::endl;
        }

        std::string drawOptionEntry = plot->getParameter<std::string>("drawOptionEntry");
        //std::cout << "drawOptionEntry = " << drawOptionEntry << std::endl;

        std::string legendEntry = "", legendEntryErrorBand = "";
        if ( plot->exists("legendEntry") ) {
          legendEntry = plot->getParameter<std::string>("legendEntry");
          legendEntryErrorBand = ( plot->exists("legendEntryErrorBand") ) ? 
            plot->getParameter<std::string>("legendEntryErrorBand") : std::string(legendEntry).append(" Uncertainty"); 
        }
        //std::cout << "legendEntry = " << legendEntry << std::endl;
        //std::cout << "legendEntryErrorBand = " << legendEntryErrorBand << std::endl;

        vstring dqmMonitorElements = plot->getParameter<vstring>("dqmMonitorElements");
        int index = 0;
        for ( vstring::const_iterator dqmMonitorElement = dqmMonitorElements.begin();
              dqmMonitorElement != dqmMonitorElements.end(); ++dqmMonitorElement ) {
          plotDefMap[index].push_back(plotDefEntry(*dqmMonitorElement, drawOptionEntry, legendEntry, legendEntryErrorBand, process, false));
          ++index;
        }
      }
    }

//--- check that number of displayed monitor elements is the same for each plot
    unsigned numMonitorElements_ref = 0;
    bool isFirstEntry = true;
    for ( std::map<int, plotDefList>::const_iterator plot = plotDefMap.begin(); 
          plot != plotDefMap.end(); ++plot ) {
      if ( isFirstEntry ) {
        numMonitorElements_ref = plot->second.size();
        isFirstEntry = false;
      } else {
        if ( plot->second.size() != numMonitorElements_ref ) {
          edm::LogError ("TauDQMHistPlotter::TauDQMHistPlotter") << " Numbers of dqmMonitorElements must be the same for all plots" 
                                                           << " --> skipping drawJob = " << (*drawJobName) << " !!";
          cfgError_ = 1;
        }
      }
    }

//--- expand process directories in names of dqmMonitorElements
    for ( std::map<int, plotDefList>::iterator plot = plotDefMap.begin(); 
          plot != plotDefMap.end(); ++plot ) {
      for ( plotDefList::iterator entry = plot->second.begin();
            entry != plot->second.end(); ++entry ) {
        std::string dqmMonitorElement = entry->dqmMonitorElement_;
        std::string process = entry->process_;

        std::map<std::string, cfgEntryProcess>::const_iterator it = processes_.find(process);
        if ( it != processes_.end() ) {
          std::string process_dqmDirectory = it->second.dqmDirectory_;

          //std::cout << "replacing processDir = " << process_dqmDirectory << " in drawJob = " << (*drawJobName) << std::endl;

          int errorFlag = 0;
          std::string dqmMonitorElement_expanded = replace_string(dqmMonitorElement, processDirKeyword, process_dqmDirectory, 0, 1, errorFlag);
          //std::cout << " dqmMonitorElement_expanded = " << dqmMonitorElement_expanded << std::endl;

          if ( !errorFlag ) {
            entry->dqmMonitorElement_ = dqmMonitorElement_expanded;
          } else {
            cfgError_ = 1;
          }
        } else {
          edm::LogError ("TauDQMHistPlotter::TauDQMHistPlotter") << " Undefined process = " << process << " !!";
          cfgError_ = 1;
        }
      }
    }

    std::string title = ( drawJob.exists("title") ) ? drawJob.getParameter<std::string>("title") : ""; 

    std::string xAxis = drawJob.getParameter<std::string>("xAxis"); 
    std::string yAxis = drawJob.getParameter<std::string>("yAxis");
    
    std::string legend = drawJob.getParameter<std::string>("legend");
    
    vstring labels = ( drawJob.exists("labels") ) ? drawJob.getParameter<vstring>("labels") : vstring();

//--- expand parameters in names of dqmMonitorElements;
//    create drawJob objects
    for ( std::map<int, plotDefList>::iterator plot = plotDefMap.begin(); 
          plot != plotDefMap.end(); ++plot ) {
      if ( drawJob.exists("parameter") ) {
        vstring vparameter = drawJob.getParameter<vstring>("parameter");
        //std::cout << "replacing parameter = " << format_vstring(vparameter) << " in drawJob = " << (*drawJobName) << std::endl;

        for ( vstring::const_iterator parameter = vparameter.begin();
              parameter != vparameter.end(); ++parameter ) {

          plotDefList plot_expanded;

          for ( plotDefList::const_iterator entry = plot->second.begin();
                entry != plot->second.end(); ++entry ) {
            std::string dqmMonitorElement = entry->dqmMonitorElement_;
            
            int errorFlag = 0;
            std::string dqmMonitorElement_expanded = replace_string(dqmMonitorElement, parKeyword, *parameter, 1, 1, errorFlag);
            //std::cout << " dqmMonitorElement_expanded = " << dqmMonitorElement_expanded << std::endl;
            if ( !errorFlag ) {
              plot_expanded.push_back(plotDefEntry(dqmMonitorElement_expanded, entry->drawOptionEntry_, 
                                                   entry->legendEntry_, entry->legendEntryErrorBand_, entry->process_, entry->doStack_));
            } else {
              cfgError_ = 1;
            }
          }

          int errorFlag = 0;
          std::string title_expanded = replace_string(title, parKeyword, *parameter, 0, 1, errorFlag);
          //std::cout << " title_expanded = " << title_expanded << std::endl;
          std::string xAxis_expanded = replace_string(xAxis, parKeyword, *parameter, 0, 1, errorFlag);
          //std::cout << " xAxis_expanded = " << xAxis_expanded << std::endl;
          std::string yAxis_expanded = replace_string(yAxis, parKeyword, *parameter, 0, 1, errorFlag);
          //std::cout << " yAxis_expanded = " << yAxis_expanded << std::endl;
          if ( errorFlag ) cfgError_ = 1;

          drawJobs_.push_back(cfgEntryDrawJob(std::string(*drawJobName).append(*parameter),
                                              plot_expanded, title_expanded, xAxis_expanded, yAxis_expanded, legend, labels));
        }
      } else {
        drawJobs_.push_back(cfgEntryDrawJob(*drawJobName, 
                                            plot->second, title, xAxis, yAxis, legend, labels));
      }
    }
  }

//--- check that all information neccessary to process drawJob is defined;
  for ( std::list<cfgEntryDrawJob>::const_iterator drawJob = drawJobs_.begin();
        drawJob != drawJobs_.end(); ++drawJob ) {
    for ( plotDefList::const_iterator plot = drawJob->plots_.begin();
          plot != drawJob->plots_.end(); ++plot ) {
      checkCfgDef<cfgEntryDrawOption>(plot->drawOptionEntry_, drawOptionEntries_, cfgError_, "drawOptionEntry", drawJob->name_);
      checkCfgDef<cfgEntryProcess>(plot->process_, processes_, cfgError_, "process", drawJob->name_);
    }

    checkCfgDef<cfgEntryAxisX>(drawJob->xAxis_, xAxes_, cfgError_, "xAxis", drawJob->name_);
    checkCfgDef<cfgEntryAxisY>(drawJob->yAxis_, yAxes_, cfgError_, "yAxis", drawJob->name_);

    checkCfgDef<cfgEntryLegend>(drawJob->legend_, legends_, cfgError_, "legend", drawJob->name_);

    checkCfgDefs<cfgEntryLabel>(drawJob->labels_, labels_, cfgError_, "label", drawJob->name_);
  }

//--- configure canvas size
  //std::cout << "--> configuring canvas size..." << std::endl;
  canvasSizeX_ = ( cfg.exists("canvasSizeX") ) ? cfg.getParameter<int>("canvasSizeX") : defaultCanvasSizeX;
  canvasSizeY_ = ( cfg.exists("canvasSizeY") ) ? cfg.getParameter<int>("canvasSizeY") : defaultCanvasSizeY;
 
//--- configure output files
  //std::cout << "--> configuring postscript output file..." << std::endl;

  outputFilePath_ = ( cfg.exists("outputFilePath") ) ? cfg.getParameter<std::string>("outputFilePath") : "";
  if ( outputFilePath_.rbegin() != outputFilePath_.rend() ) {
    if ( (*outputFilePath_.rbegin()) == '/' ) outputFilePath_.erase(outputFilePath_.length() - 1);
  }
  //std::cout << " outputFilePath = " << outputFilePath_ << std::endl;

  outputFileName_ = ( cfg.exists("outputFileName") ) ? cfg.getParameter<std::string>("outputFileName") : "";
  //std::cout << " outputFileName = " << outputFileName_ << std::endl;

  indOutputFileName_ = ( cfg.exists("indOutputFileName") ) ? cfg.getParameter<std::string>("indOutputFileName") : "";
  if ( indOutputFileName_ != "" && indOutputFileName_.find('.') == std::string::npos ) {
    edm::LogError ("TauDQMHistPlotter") << " Failed to determine type of graphics format from indOutputFileName = " << indOutputFileName_ << " !!";
    cfgError_ = 1;
  }
  //std::cout << " indOutputFileName = " << indOutputFileName_ << std::endl;

//--- check that exactly one type of output is specified for the plots
//    (either separate graphics files displaying one plot each 
//     or postscript file displaying all plots on successive pages;
//     cannot create both types of output simultaneously,
//     as TCanvas::Print seems to interfere with TPostScript::NewPage)
  if ( outputFileName_ == "" && indOutputFileName_ == "" ) {
    edm::LogError ("TauDQMHistPlotter") << " Either outputFileName or indOutputFileName must be specified !!";
    cfgError_ = 1;
  }

  if ( outputFileName_ != "" && indOutputFileName_ != "" ) {
    edm::LogError ("TauDQMHistPlotter") << " Must not specify outputFileName and indOutputFileName simultaneously !!";
    cfgError_ = 1;
  }

  if( verbosity )
    std::cout << "done." << std::endl;
}

TauDQMHistPlotter::~TauDQMHistPlotter() 
{
// nothing to be done yet...
}

void TauDQMHistPlotter::analyze(const edm::Event&, const edm::EventSetup&)
{
// nothing to be done yet...
}

void TauDQMHistPlotter::endRun(const edm::Run& r, const edm::EventSetup& c)
{
  if( verbosity )
    std::cout << "<TauDQMHistPlotter::endJob>:" << std::endl;

//--- check that configuration parameters contain no errors
  if ( cfgError_ ) {
    edm::LogError ("endJob") << " Error in Configuration ParameterSet --> histograms will NOT be plotted !!";
    return;
  }

//--- check that DQMStore service is available
  if ( !edm::Service<DQMStore>().isAvailable() ) {
    edm::LogError ("endJob") << " Failed to access dqmStore --> histograms will NOT be plotted !!";
    return;
  }

  DQMStore& dqmStore = (*edm::Service<DQMStore>());

//--- stop ROOT from keeping references to all hsitograms
  //TH1::AddDirectory(false);

//--- stop ROOT from opening X-window for canvas output
//    (in order to be able to run in batch mode) 
  gROOT->SetBatch(true);

//--- initialize graphical output;
//    open postscript file
  TCanvas canvas("TauDQMHistPlotter","TauDQMHistPlotter", canvasSizeX_, canvasSizeY_);
  canvas.SetFillColor(10);

//--- restrict area in which histograms are drawn to quadratic TPad in the center of the TCanvas,
//    in order to make space for axis labels...
  //TPad pad("EWKTauPad", "EWKTauPad", 0.02, 0.15, 0.98, 0.85);
  //pad.SetFillColor(10);
  //pad.Draw();
  //pad.Divide(1,1);
  //pad.cd(1);

  TPostScript* ps = NULL;
  if ( outputFileName_ != "" ) {
    std::string psFileName = ( outputFilePath_ != "" ) ? std::string(outputFilePath_).append("/").append(outputFileName_) : outputFileName_;
    ps = new TPostScript(psFileName.data(), 112);
  }
    
//--- process drawJobs
  for ( std::list<cfgEntryDrawJob>::const_iterator drawJob = drawJobs_.begin(); 
        drawJob != drawJobs_.end(); ++drawJob ) {
    const std::string& drawJobName = drawJob->name_;
    if( verbosity )
      std::cout << "--> processing drawJob " << drawJobName << "..." << std::endl;

//--- prepare internally used histogram data-structures
    TH1* stackedHistogram_sum = NULL;
    std::list<TH1*> histogramsToDelete;
    std::list<plotDefEntry*> drawOptionsToDelete;

    typedef std::pair<TH1*, const plotDefEntry*> histogram_drawOption_pair;
    std::list<histogram_drawOption_pair> allHistograms;

    for ( plotDefList::const_iterator plot = drawJob->plots_.begin();
          plot != drawJob->plots_.end(); ++plot ) {

      std::string dqmMonitorElementName_full = dqmDirectoryName(std::string(dqmRootDirectory)).append(plot->dqmMonitorElement_);
      if( verbosity )
        std::cout << " dqmMonitorElementName_full = " << dqmMonitorElementName_full << std::endl;
      MonitorElement* dqmMonitorElement = dqmStore.get(dqmMonitorElementName_full);

      TH1* histogram = dqmMonitorElement->getTH1F();
      if(verbosity)
        std::cout<<"Got Histogram "<<std::endl;
      //      TH1* histogram = ( dqmMonitorElement ) ? dynamic_cast<TH1*>(dqmMonitorElement->getTH1()->Clone()) : NULL;
      //histogramsToDelete.push_back(histogram);

      if ( histogram == NULL ) {
        edm::LogError ("endJob") << " Failed to access dqmMonitorElement = " << dqmMonitorElementName_full <<","
                                 << " needed by drawJob = " << drawJobName << " --> histograms will NOT be plotted !!";
        continue;
      }

      if ( !histogram->GetSumw2N() ) histogram->Sumw2();

      const cfgEntryDrawOption* drawOptionConfig = 
        findCfgDef<cfgEntryDrawOption>(plot->drawOptionEntry_, drawOptionEntries_, "drawOptionEntry", drawJobName);
      if ( drawOptionConfig == NULL ) {
        edm::LogError ("endJob") << " Failed to access information needed by drawJob = " << drawJobName 
                                 << " --> histograms will NOT be plotted !!";
        return;
      }
      
      if ( drawOptionConfig->drawOption_ == drawOption_eBand ) {
//--- add histogram displaying central value as solid line
        TH1* histogram_centralValue = dynamic_cast<TH1*>(histogram->Clone());
        histogram_centralValue->SetName(std::string(histogram->GetName()).append("_centralValue").data());
        cfgEntryDrawOption drawOptionConfig_centralValue(*drawOptionConfig);
        drawOptionConfig_centralValue.fillColor_ = 0;
        drawOptionConfig_centralValue.fillStyle_ = 0;
        drawOptionConfig_centralValue.drawOption_ = "hist";
        drawOptionConfig_centralValue.drawOptionLegend_ = "l";
        std::string drawOptionName_centralValue = std::string(plot->drawOptionEntry_).append("_centralValue");
//--- entries in std::map need to be unique,
//    so need to check whether drawOptionEntry already exists...
        if ( drawOptionEntries_.find(drawOptionName_centralValue) == drawOptionEntries_.end() ) 
          drawOptionEntries_.insert(std::pair<std::string, cfgEntryDrawOption>
                                    (drawOptionName_centralValue, cfgEntryDrawOption(drawOptionName_centralValue, drawOptionConfig_centralValue)));
        plotDefEntry* plot_centralValue = new plotDefEntry(*plot);
        plot_centralValue->drawOptionEntry_ = drawOptionName_centralValue;
        allHistograms.push_back(histogram_drawOption_pair(histogram_centralValue, plot_centralValue));
        histogramsToDelete.push_back(histogram_centralValue);
        drawOptionsToDelete.push_back(plot_centralValue);

//--- add histogram displaying uncertainty as shaded error band
        TH1* histogram_ErrorBand = dynamic_cast<TH1*>(histogram->Clone());
        histogram_ErrorBand->SetName(std::string(histogram->GetName()).append("_ErrorBand").data());
        cfgEntryDrawOption drawOptionConfig_ErrorBand(*drawOptionConfig);
        drawOptionConfig_ErrorBand.markerColor_ = drawOptionConfig_ErrorBand.fillColor_;
        drawOptionConfig_ErrorBand.markerSize_ = 0.;
        drawOptionConfig_ErrorBand.lineColor_ = drawOptionConfig_ErrorBand.fillColor_;
        drawOptionConfig_ErrorBand.lineWidth_ = 0;
        drawOptionConfig_ErrorBand.drawOption_ = "e2";
        drawOptionConfig_ErrorBand.drawOptionLegend_ = "f";
        std::string drawOptionName_ErrorBand = std::string(plot->drawOptionEntry_).append("_ErrorBand");
//--- entries in std::map need to be unique,
//    so need to check whether drawOptionEntry already exists...
        if ( drawOptionEntries_.find(drawOptionName_ErrorBand) == drawOptionEntries_.end() ) 
          drawOptionEntries_.insert(std::pair<std::string, cfgEntryDrawOption>
                                    (drawOptionName_ErrorBand, cfgEntryDrawOption(drawOptionName_ErrorBand, drawOptionConfig_ErrorBand)));
        plotDefEntry* plot_ErrorBand = new plotDefEntry(*plot);
        plot_ErrorBand->drawOptionEntry_ = drawOptionName_ErrorBand;
        plot_ErrorBand->isErrorBand_ = true;
        allHistograms.push_back(histogram_drawOption_pair(histogram_ErrorBand, plot_ErrorBand));
        histogramsToDelete.push_back(histogram_ErrorBand);
        drawOptionsToDelete.push_back(plot_ErrorBand);
      } else if ( plot->doStack_ ) {
        TH1* stackedHistogram = dynamic_cast<TH1*>(histogram->Clone());
        if ( stackedHistogram_sum ) stackedHistogram->Add(stackedHistogram_sum);
        stackedHistogram_sum = stackedHistogram;
        histogramsToDelete.push_back(stackedHistogram);
        allHistograms.push_back(histogram_drawOption_pair(stackedHistogram, &(*plot)));
      } else {  
        allHistograms.push_back(histogram_drawOption_pair(histogram, &(*plot)));
      }
    }

//--- determine normalization of y-axis
//    (maximum of any of the histograms included in drawJob)
    double yAxisNorm = 0.;
    for ( std::list<histogram_drawOption_pair>::const_iterator it = allHistograms.begin();
          it != allHistograms.end(); ++it ) {
      yAxisNorm = TMath::Max(yAxisNorm, it->first->GetMaximum());
    }
    //std::cout << " yAxisNorm = " << yAxisNorm << std::endl;
    cfgEntryAxisY::setNorm(yAxisNorm);

//--- prepare histograms for drawing
    const cfgEntryAxisX* xAxisConfig = findCfgDef<cfgEntryAxisX>(drawJob->xAxis_, xAxes_, "xAxis", drawJobName);
    const cfgEntryAxisY* yAxisConfig = findCfgDef<cfgEntryAxisY>(drawJob->yAxis_, yAxes_, "yAxis", drawJobName);
    const cfgEntryLegend* legendConfig = findCfgDef<cfgEntryLegend>(drawJob->legend_, legends_, "legend", drawJobName);
    if ( xAxisConfig == NULL || yAxisConfig == NULL || legendConfig == NULL ) {
      edm::LogError ("endJob") << " Failed to access information needed by drawJob = " << drawJobName 
                               << " --> histograms will NOT be plotted !!";
      return;
    }

//--- WARNING: need to call 
//              TLegend::TLegend(Double_t, Double_t,Double_t, Double_t, const char* = "", Option_t* = "brNDC")
//             constructor, as TLegend::TLegend default constructor causes the created TLegend object to behave differently !!
    TLegend legend(defaultLegendPosX, defaultLegendPosY, defaultLegendPosX + defaultLegendSizeX, defaultLegendPosY + defaultLegendSizeY);
    legendConfig->applyTo(&legend);

    std::list<histoDrawEntry> smProcessHistogramList;
    std::list<histoDrawEntry> bsmProcessHistogramList;
    std::list<histoDrawEntry> smSumHistogramList;
    std::list<histoDrawEntry> smSumUncertaintyHistogramList;
    std::list<histoDrawEntry> dataHistogramList;

    for ( std::list<histogram_drawOption_pair>::const_iterator it = allHistograms.begin();
          it != allHistograms.end(); ++it ) {
      TH1* histogram = it->first;
      const plotDefEntry* drawOption = it->second;

      const cfgEntryDrawOption* drawOptionConfig = 
        findCfgDef<cfgEntryDrawOption>(drawOption->drawOptionEntry_, drawOptionEntries_, "drawOptionEntry", drawJobName);
      const cfgEntryProcess* processConfig = findCfgDef<cfgEntryProcess>(drawOption->process_, processes_, "process", drawJobName);
      if ( drawOptionConfig == NULL || processConfig == NULL ) {
        edm::LogError ("endJob") << " Failed to access information needed by drawJob = " << drawJobName 
                                 << " --> histograms will NOT be plotted !!";
        return;
      }
      
      if ( drawJob->title_ != "" ) histogram->SetTitle(drawJob->title_.data());

      xAxisConfig->applyTo(histogram);
      yAxisConfig->applyTo(histogram);

      bool yLogScale = ( yAxisConfig->yScale_ == yScale_log ) ? true : false;
      //std::cout << " yLogScale = " << yLogScale << std::endl; 
      //pad.SetLogy(yLogScale);
      canvas.SetLogy(yLogScale);

      drawOptionConfig->applyTo(histogram);
      histogram->SetStats(false);

      if ( drawOption->isErrorBand_ ) {
        smSumUncertaintyHistogramList.push_back(histoDrawEntry(histogram, drawOptionConfig->drawOption_.data()));
      } else {
        if ( processConfig->type_ == type_smMC ) {
          smProcessHistogramList.push_back(histoDrawEntry(histogram, drawOptionConfig->drawOption_.data()));
        } else if ( processConfig->type_ == type_bsmMC ) {
          bsmProcessHistogramList.push_back(histoDrawEntry(histogram, drawOptionConfig->drawOption_.data()));
        } else if ( processConfig->type_ == type_smSumMC ) {
          smSumHistogramList.push_back(histoDrawEntry(histogram, drawOptionConfig->drawOption_.data()));
        } else if ( processConfig->type_ == type_Data ) {
          dataHistogramList.push_back(histoDrawEntry(histogram, drawOptionConfig->drawOption_.data()));
        } 
      }

      std::string legendEntry, legendDrawOption;
      if ( drawOption->isErrorBand_ ) {
        legendEntry = ( drawOption->legendEntryErrorBand_ != "" ) ? drawOption->legendEntryErrorBand_ : processConfig->legendEntryErrorBand_;
        legendDrawOption = "f";
      } else {
        legendEntry = ( drawOption->legendEntry_ != "" ) ? drawOption->legendEntry_ : processConfig->legendEntry_;
        legendDrawOption = drawOptionConfig->drawOptionLegend_;
      } 

      legend.AddEntry(histogram, legendEntry.data(), legendDrawOption.data());
    }

    std::list<TPaveText> labels;
    for ( vstring::const_iterator labelName = drawJob->labels_.begin();
          labelName != drawJob->labels_.end(); ++labelName ) {
      const cfgEntryLabel* labelConfig = findCfgDef<cfgEntryLabel>(*labelName, labels_, "label", drawJobName);
      
      TPaveText label;
      labelConfig->applyTo(&label);

      labels.push_back(label);
    }

//--- draw histograms 
//   - in the order:
//    1. uncertainty on sum of all Standard Model processes
//    2. sum of all Standard Model processes
//    3. individual Standard Model processes
//    4. individual beyond the Standard Model processes
//    5. data
    bool isFirstHistogram = true;
    drawHistograms(smSumUncertaintyHistogramList, isFirstHistogram);
    drawHistograms(smSumHistogramList, isFirstHistogram);

//--- process histograms for individual Standard Model processes
//    in reverse order, so that most stacked histogram gets drawn first
    for ( std::list<histoDrawEntry>::reverse_iterator it = smProcessHistogramList.rbegin();
          it != smProcessHistogramList.rend(); ++it ) {
      std::string drawOption = ( isFirstHistogram ) ? it->second : std::string(it->second).append("same");
      it->first->Draw(drawOption.data());
      isFirstHistogram = false;
    }
        
    drawHistograms(bsmProcessHistogramList, isFirstHistogram);
    drawHistograms(dataHistogramList, isFirstHistogram);

    legend.Draw();

    for ( std::list<TPaveText>::iterator label = labels.begin();
          label != labels.end(); ++label ) {
      label->Draw();
    }

    //pad.RedrawAxis();
    
    canvas.Update();
    //pad.Update();

    if ( indOutputFileName_ != "" && toFile_) {
      int errorFlag = 0;
      std::string modIndOutputFileName = replace_string(indOutputFileName_, plotKeyword, drawJobName, 1, 1, errorFlag);
      if ( !errorFlag ) {
        std::string fullFileName = ( outputFilePath_ != "" ) ? 
          std::string(outputFilePath_).append("/").append(modIndOutputFileName) : modIndOutputFileName;
        canvas.Print(fullFileName.data());
      } else {
        edm::LogError("endJob") << " Failed to decode indOutputFileName = " << indOutputFileName_ << " --> skipping !!";
      }
    }

    if ( ps ) ps->NewPage();

//--- delete temporarily created histogram and drawOption objects
    for ( std::list<TH1*>::const_iterator histogram = histogramsToDelete.begin();
          histogram != histogramsToDelete.end(); ++histogram ) {
      delete (*histogram);
    }

    for ( std::list<plotDefEntry*>::const_iterator drawOption = drawOptionsToDelete.begin();
          drawOption != drawOptionsToDelete.end(); ++drawOption ) {
      delete (*drawOption);
    }
  }

//--- close postscript file
  canvas.Clear();
  if(verbosity)
    std::cout << "done." << std::endl;
  if ( ps ) ps->Close();
  delete ps;
}

#include "FWCore/Framework/interface/MakerMacros.h"

DEFINE_FWK_MODULE(TauDQMHistPlotter);