DQMHistNormalizer Class Reference

Inheritance diagram for DQMHistNormalizer:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override
	DQMHistNormalizer (const edm::ParameterSet &)
virtual void	endJob () override
virtual void	endRun (const edm::Run &r, const edm::EventSetup &c) override
virtual	~DQMHistNormalizer ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
lat::Regexp *	buildRegex (const string &expr)

Private Attributes
vector< string >	plotNamesToNormalize_
string	reference_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Class to produce efficiency histograms by dividing nominator by denominator histograms

Author

Christian Veelken, UC Davis

Definition at line 36 of file DQMHistNormalizer.cc.

Constructor & Destructor Documentation

DQMHistNormalizer::DQMHistNormalizer ( const edm::ParameterSet &  cfg )

explicit

Definition at line 52 of file DQMHistNormalizer.cc.

                                                              :
  plotNamesToNormalize_(cfg.getParameter< std::vector<string> >("plotNamesToNormalize")),
  reference_(cfg.getParameter< string >("reference"))
{
  //std::cout << "<DQMHistNormalizer::DQMHistNormalizer>:" << std::endl;
}

DQMHistNormalizer::~DQMHistNormalizer ( )

virtual

Definition at line 59 of file DQMHistNormalizer.cc.

{
//--- nothing to be done yet
}

Member Function Documentation

void DQMHistNormalizer::analyze ( const edm::Event &  , const edm::EventSetup &    )

overridevirtual

Definition at line 64 of file DQMHistNormalizer.cc.

{
//--- nothing to be done yet
}

lat::Regexp * DQMHistNormalizer::buildRegex ( const string &  expr )

private

Definition at line 69 of file DQMHistNormalizer.cc.

References MillePedeFileConverter_cfg::e, raiseDQMError(), and corrVsCorr::rx.

Referenced by endRun().

{
  lat::Regexp* rx = 0;
  try
  {
    rx = new lat::Regexp(expr, 0, lat::Regexp::Wildcard);
    rx->study();
  }
  catch (lat::Error &e)
  {
    raiseDQMError("DQMStore", "Invalid regular expression '%s': %s",
          expr.c_str(), e.explain().c_str());
  }
  return rx;
}

virtual void DQMHistNormalizer::endJob ( void  )

inlineoverridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 43 of file DQMHistNormalizer.cc.

References EnergyCorrector::c, Ecal2004TBTDCRanges_v1_cff::endRun, and alignCSCRings::r.

{}

void DQMHistNormalizer::endRun ( const edm::Run &  r, const edm::EventSetup &  c  )

overridevirtual

Definition at line 85 of file DQMHistNormalizer.cc.

References buildRegex(), DEFINE_FWK_MODULE, dir, fastjetJetProducer_validation_cfg::DQMStore, create_public_lumi_plots::hist, plotNamesToNormalize_, and reference_.

{
  //std::cout << "<DQMHistNormalizer::endJob>:" << std::endl;

//--- 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>());

  vector<MonitorElement *> allOurMEs = dqmStore.getAllContents("RecoTauV/");
  lat::Regexp* refregex = buildRegex("*RecoTauV/*/" + reference_);
  vector<lat::Regexp*> toNormRegex;
  for ( std::vector<string>::const_iterator toNorm = plotNamesToNormalize_.begin(); toNorm != plotNamesToNormalize_.end(); ++toNorm )
    toNormRegex.push_back( buildRegex("*RecoTauV/*/" + *toNorm) );

  map<string, MonitorElement *> refsMap;
  vector<MonitorElement *> toNormElements;

  for(vector<MonitorElement *>::const_iterator element = allOurMEs.begin(); element != allOurMEs.end(); ++element){
    string pathname = (*element)->getFullname();
    //cout << pathname << endl;
    //Matches reference
    if(refregex->match(pathname)){
      //cout << "Matched to ref" << endl;
      string dir = pathname.substr(0, pathname.rfind("/"));
      if(refsMap.find(dir) != refsMap.end()){
    edm::LogInfo("DQMHistNormalizer")<<"DQMHistNormalizer::endRun: Warning! found multiple normalizing references for dir: "<<dir<<"!";
    edm::LogInfo("DQMHistNormalizer")<<"     " << (refsMap[dir])->getFullname();
    edm::LogInfo("DQMHistNormalizer")<<"     " << pathname;
      }
      else{
    refsMap[dir] = *element;
      }
    }

    //Matches targets
    for(vector<lat::Regexp*>::const_iterator reg = toNormRegex.begin(); reg != toNormRegex.end(); ++reg){
      if((*reg)->match(pathname)){
    //cout << "Matched to target" << endl;
    toNormElements.push_back(*element);
    //cout << "Filled the collection" << endl;
      }
    }
  }

  delete refregex;
  for(vector<lat::Regexp*>::const_iterator reg = toNormRegex.begin(); reg != toNormRegex.end(); ++reg)
    delete *reg;

  for(vector<MonitorElement *>::const_iterator matchingElement = toNormElements.begin(); matchingElement != toNormElements.end(); ++matchingElement){
    string meName = (*matchingElement)->getFullname();
    string dir = meName.substr(0, meName.rfind("/"));
    
    if(refsMap.find(dir) == refsMap.end()){
      edm::LogInfo("DQMHistNormalizer")<<"DQMHistNormalizer::endRun: Error! normalizing references for "<<meName<<" not found! Skipping...";
      continue;
    }
      
    float norm = refsMap[dir]->getTH1()->GetEntries();
    TH1* hist = (*matchingElement)->getTH1();
    if ( norm != 0. ) {
      if( !hist->GetSumw2N() ) hist->Sumw2();
      hist->Scale(1/norm);//use option "width" to divide the bin contents and errors by the bin width?
    }else{
      edm::LogInfo("DQMHistNormalizer")<<"DQMHistNormalizer::endRun: Error! Normalization failed in "<<hist->GetTitle()<<"!";
    }
      
  }//    for(vector<MonitorElement *>::const_iterator matchingElement = matchingElemts.begin(); matchingElement = matchingElemts.end(); ++matchingElement)
}

Member Data Documentation

vector<string> DQMHistNormalizer::plotNamesToNormalize_

private

Definition at line 48 of file DQMHistNormalizer.cc.

Referenced by endRun().

string DQMHistNormalizer::reference_

private

Definition at line 49 of file DQMHistNormalizer.cc.

Referenced by endRun().