dqmAuxFunctions.cc

Go to the documentation of this file.

#include "Validation/RecoTau/plugins/dqmAuxFunctions.h"
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "DQMServices/Core/interface/DQMStore.h"
 
#include <TH1.h>
 
#include <iostream>
 
std::string replace_string(const std::string& src,
                           const std::string& keyword,
                           const std::string& parameter,
                           unsigned minReplacements,
                           unsigned maxReplacements,
                           int& errorFlag) {
  std::string modSrc = src;
  unsigned numReplacements = 0;
  while (modSrc.find(keyword) != std::string::npos) {
    modSrc.replace(modSrc.find(keyword), keyword.length(), parameter);
    ++numReplacements;
  }
  if ((numReplacements < minReplacements) || (numReplacements > maxReplacements)) {
    edm::LogError("replace_string") << " Failed to replace parameter = " << parameter << " in src = " << src << ","
                                    << " numReplacements = " << numReplacements << " (min = " << minReplacements
                                    << ", max = " << maxReplacements << ") !!";
    errorFlag = 1;
  }
  return modSrc;
}
 
//
//-----------------------------------------------------------------------------------------------------------------------
//
 
std::string format_vstring(const std::vector<std::string>& vs) {
  std::ostringstream os;
 
  os << "{ ";
 
  unsigned numEntries = vs.size();
  for (unsigned iEntry = 0; iEntry < numEntries; ++iEntry) {
    os << vs[iEntry];
    if (iEntry < (numEntries - 1))
      os << ", ";
  }
 
  os << " }";
 
  return os.str();
}
 
//
//-----------------------------------------------------------------------------------------------------------------------
//
 
std::string dqmDirectoryName(const std::string& directory) {
  std::string dirName = directory;
  //if ( dirName == "" || dirName.find_last_of(dqmSeparator) != (dirName.length() - 1) )  dirName.append(dqmSeparator);
  //--- add tailing '/'
  if (!dirName.empty() && dirName.find_last_of(dqmSeparator) != (dirName.length() - 1))
    dirName.append(dqmSeparator);
  return dirName;
}
 
std::string dqmSubDirectoryName_merged(const std::string& directory, const std::string& subdirectory) {
  std::string subDirName = subdirectory;
  //--- remove characters specifying directory part from name of subdirectory
  if (subDirName.find(directory) <= 1)
    subDirName.replace(subDirName.find(directory), directory.length(), "");
  //--- remove tailing '/'s
  while (subDirName.find(dqmSeparator) == 0)
    subDirName.replace(subDirName.find(dqmSeparator), dqmSeparator.length(), "");
  return subDirName;
}
 
//
//-----------------------------------------------------------------------------------------------------------------------
//
 
void dqmRegisterHistogram(dqm::legacy::DQMStore& dqmStore, TH1* histogram, const std::string& name) {
  //std::cout << "<dqmRegisterHistogram>:" << std::endl;
  //histogram->SetName(std::string(histogram->GetName()).append("_copied").data());
  histogram->SetName(histogram->GetName());
  if (TH1F* h = dynamic_cast<TH1F*>(histogram)) {
    //std::cout << " --> calling DQMStore::book1D" << std::endl;
    dqmStore.book1D(name, h);
  } else if (TH1S* h = dynamic_cast<TH1S*>(histogram)) {
    //std::cout << " --> calling DQMStore::book1@" << std::endl;
    dqmStore.book1S(name, h);
  } else if (TH2F* h = dynamic_cast<TH2F*>(histogram)) {
    //std::cout << " --> calling DQMStore::book2D" << std::endl;
    dqmStore.book2D(name, h);
  } else if (TH2S* h = dynamic_cast<TH2S*>(histogram)) {
    //std::cout << " --> calling DQMStore::book2S" << std::endl;
    dqmStore.book2S(name, h);
  } else if (TH3F* h = dynamic_cast<TH3F*>(histogram)) {
    //std::cout << " --> calling DQMStore::book3D" << std::endl;
    dqmStore.book3D(name, h);
  } else if (TProfile* h = dynamic_cast<TProfile*>(histogram)) {
    //std::cout << " --> calling DQMStore::bookProfile" << std::endl;
    dqmStore.bookProfile(name, h);
  } else if (TProfile2D* h = dynamic_cast<TProfile2D*>(histogram)) {
    //std::cout << " --> calling DQMStore::bookProfile2D" << std::endl;
    dqmStore.bookProfile2D(name, h);
  }
}
 
void dqmCopyRecursively(dqm::legacy::DQMStore& dqmStore,
                        const std::string& inputDirectory,
                        const std::string& outputDirectory,
                        double scaleFactor,
                        int mode,
                        bool rmInputDirectory) {
  //std::cout << "<copyRecursively>:" << std::endl;
  //std::cout << " inputDirectory = " << inputDirectory << std::endl;
  //std::cout << " outputDirectory = " << outputDirectory << std::endl;
  //std::cout << " rmInputDirectory = " << rmInputDirectory << std::endl;
 
  //--- copy all monitor elements in current inputDirectory to the outputDirectory
  dqmStore.setCurrentFolder(inputDirectory);
  std::vector<std::string> meNames = dqmStore.getMEs();
  for (std::vector<std::string>::const_iterator meName = meNames.begin(); meName != meNames.end(); ++meName) {
    std::string meName_full = dqmDirectoryName(inputDirectory).append(*meName);
    //std::cout << " meName_full = " <<  meName_full << std::endl;
 
    dqmStore.setCurrentFolder(inputDirectory);
    dqm::legacy::MonitorElement* meInput = dqmStore.get(meName_full);
    //std::cout << " meInput = " << meInput << std::endl;
    if (!meInput) {
      edm::LogError("copyRecursively") << " Failed to access meName = " << (*meName) << " in DQMStore"
                                       << " --> skipping !!";
      continue;
    }
 
    TH1* histogram = meInput->getTH1();
    //std::cout << " histogram = " << histogram << std::endl;
    if (!histogram) {
      edm::LogError("copyRecursively") << " Failed to access histogram associated to meName = " << (*meName)
                                       << " in DQMStore"
                                       << " --> skipping !!";
      continue;
    }
 
    std::unique_ptr<TH1> clone(dynamic_cast<TH1*>(histogram->Clone()));
    clone->Scale(scaleFactor);
 
    dqmStore.setCurrentFolder(outputDirectory);
    dqm::legacy::MonitorElement* meOutput = dqmStore.get(dqmDirectoryName(outputDirectory).append(*meName));
    //std::cout << " meOutput = " << meOutput << std::endl;
    //--- check if outputHistogram does already exist
    if (meOutput) {
      switch (mode) {
        case 1:  // print error message
          edm::LogError("copyRecursively")
              << " meName = " << (*meName) << " already exists in outputDirectory = " << outputDirectory
              << " --> skipping !!";
          break;
        case 2:  // overwrite outputHistogram
          dqmRegisterHistogram(dqmStore, clone.release(), *meName);
          break;
        case 3:  // add histogram to outputHistogram
          meOutput->getTH1()->Add(clone.get(), scaleFactor);
      }
    } else {
      dqmRegisterHistogram(dqmStore, clone.release(), *meName);
    }
  }
 
  //--- call function recursively for all sub-directories
  dqmStore.setCurrentFolder(inputDirectory);
  std::vector<std::string> dirNames = dqmStore.getSubdirs();
  for (std::vector<std::string>::const_iterator dirName = dirNames.begin(); dirName != dirNames.end(); ++dirName) {
    std::string subDirName = dqmSubDirectoryName_merged(inputDirectory, *dirName);
    //std::cout << " subDirName = " << subDirName << std::endl;
 
    std::string inputDirName_full = dqmDirectoryName(inputDirectory).append(subDirName);
    //std::cout << " inputDirName_full = " << inputDirName_full << std::endl;
 
    std::string outputDirName_full = dqmDirectoryName(outputDirectory).append(subDirName);
    //std::cout << " outputDirName_full = " << outputDirName_full << std::endl;
 
    dqmCopyRecursively(dqmStore, inputDirName_full, outputDirName_full, scaleFactor, mode, rmInputDirectory);
  }
}
 
//
//-----------------------------------------------------------------------------------------------------------------------
//
 
void separateHistogramFromDirectoryName(const std::string& histogramAndDirectoryName,
                                        std::string& histogramName,
                                        std::string& directoryName) {
  //std::cout << "<separateHistogramFromDirectoryName>:" << std::endl;
 
  std::string tempName = histogramAndDirectoryName;
 
  //--- remove DQM root directory from histogram name
  std::string::size_type dqmRootDirectoryPos = tempName.find(dqmRootDirectory);
  if (dqmRootDirectoryPos != std::string::npos) {
    tempName.replace(dqmRootDirectoryPos, dqmRootDirectory.size(), "");
  }
 
  //std::cout << " tempName = " << tempName << std::endl;
 
  //--- extract directory from histogram name
  std::string::size_type lastPos;
  std::string::size_type nextPos = tempName.find(dqmSeparator);
  do {
    lastPos = nextPos;
    nextPos = tempName.find(dqmSeparator, lastPos + 1);
  } while (nextPos != std::string::npos);
 
  histogramName = (lastPos != std::string::npos) ? std::string(tempName, lastPos + 1, tempName.length()) : tempName;
  directoryName = (lastPos != std::string::npos) ? std::string(tempName, 0, lastPos) : "";
 
  //std::cout << " histogramName = " << histogramName << std::endl;
  //std::cout << " directoryName = " << directoryName << std::endl;
}