HLTMonElectronConsumer Class Reference

#include <DQM/HLTMonElectronConsumer/src/HLTMonElectronConsumer.cc>

Inheritance diagram for HLTMonElectronConsumer:

Public Member Functions
	HLTMonElectronConsumer (const edm::ParameterSet &)
	~HLTMonElectronConsumer ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()

Private Attributes
DQMStore *	dbe
std::string	dirname_
MonitorElement *	isocheck
std::string	isodirname_
edm::InputTag	isotag_
MonitorElement *	isototal
ofstream	logFile_
bool	monitorDaemon_
std::string	outputFile_
std::string	pixeldirname_
MonitorElement *	pixelEff
MonitorElement *	pixelhistosEt [4]
MonitorElement *	pixelhistosEta [4]
MonitorElement *	pixelhistosEtaOut [2]
MonitorElement *	pixelhistosEtOut [2]
MonitorElement *	pixelhistosPhi [4]
MonitorElement *	pixelhistosPhiOut [2]
edm::InputTag	pixeltag_
MonitorElement *	pixeltotal
MonitorElement *	trackEff

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
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

Description: This is a DQM source meant to be an example for general development of HLT DQM code. Based on the general structure used for L1TMonitor DQM sources.

Implementation: <Notes on="" implementation>="">

Definition at line 51 of file HLTMonElectronConsumer.h.

Constructor & Destructor Documentation

HLTMonElectronConsumer::HLTMonElectronConsumer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 14 of file HLTMonElectronConsumer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), edm::InputTag::label(), LogDebug, NULL, cppFunctionSkipper::operator, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  
  LogDebug("HLTMonElectronConsumer") << "constructor...." ;
  
  logFile_.open("HLTMonElectronConsumer.log");
  
  dbe = NULL;
  if (iConfig.getUntrackedParameter < bool > ("DQMStore", false)) {
    dbe = Service < DQMStore > ().operator->();
    dbe->setVerbose(0);
  }
  
  outputFile_ =
    iConfig.getUntrackedParameter <std::string>("outputFile", "");
  if (outputFile_.size() != 0) {
    LogInfo("HLTMonElectronConsumer") << "L1T Monitoring histograms will be saved to " 
                  << outputFile_ ;
  }
  else {
    outputFile_ = "L1TDQM.root";
  }
  
  bool disable =
    iConfig.getUntrackedParameter < bool > ("disableROOToutput", false);
  if (disable) {
    outputFile_ = "";
  }
  
  pixeltag_=iConfig.getParameter<edm::InputTag>("PixelTag");
  isotag_=iConfig.getParameter<edm::InputTag>("IsoTag");
  
  dirname_="HLT/HLTMonElectron/"+iConfig.getParameter<std::string>("@module_label");
  pixeldirname_="HLT/HLTMonElectron/"+pixeltag_.label();
  isodirname_="HLT/HLTMonElectron/"+isotag_.label();
  
  if (dbe != NULL) {
    dbe->setCurrentFolder(dirname_);
  }
  
}

HLTMonElectronConsumer::~HLTMonElectronConsumer

(

)

Definition at line 57 of file HLTMonElectronConsumer.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void HLTMonElectronConsumer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 72 of file HLTMonElectronConsumer.cc.

References i, j, and NULL.

{
  //total pixelmatch efficiencies from summary histo
  if(pixeltotal!=NULL){
    LogInfo("HLTMonElectronConsumer") << "  pixelhisto  " <<  pixeltotal->getBinContent(1);
    if(pixeltotal->getBinContent(1)!=0)
      pixelEff->Fill(pixeltotal->getBinContent(2)/pixeltotal->getBinContent(1));
    else
      pixelEff->Fill(0.);
    
    if(pixeltotal->getBinContent(3)!=0)
      trackEff->Fill(pixeltotal->getBinContent(4)/pixeltotal->getBinContent(3));
    else
      trackEff->Fill(0.);

    // efficiency as kinematic function
    for(int i =0; i<2 ;i++){
      TH1F* num;
      TH1F* denom;
      num=pixelhistosEt[2*i+1]->getTH1F();
      denom=pixelhistosEt[2*i]->getTH1F();
      for(int j=1; j <= pixelhistosEtOut[i]->getNbinsX();j++ ){
    if(denom->GetBinContent(j)!=0)
      pixelhistosEtOut[i]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
    else
      pixelhistosEtOut[i]->setBinContent(j,0.);
      }
      num=pixelhistosEta[2*i+1]->getTH1F();
      denom=pixelhistosEta[2*i]->getTH1F();
      for(int j=1; j <= pixelhistosEtaOut[i]->getNbinsX();j++ ){
    if(denom->GetBinContent(j)!=0)
      pixelhistosEtaOut[i]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
    else
      pixelhistosEtaOut[i]->setBinContent(j,0.);
      }
      num=pixelhistosPhi[2*i+1]->getTH1F();
      denom=pixelhistosPhi[2*i]->getTH1F();
      for(int j=1; j <= pixelhistosPhiOut[i]->getNbinsX();j++ ){
    if(denom->GetBinContent(j)!=0)
      pixelhistosPhiOut[i]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
    else
      pixelhistosPhiOut[i]->setBinContent(j,0.);
      }
    }
  }else
    LogInfo("HLTMonElectronConsumer") << " empty pixelhisto  " ;

  if(isototal!=NULL){
    TH1F* refhist = isototal->getTH1F();
    for(int i =1; i<= refhist->GetNbinsX();i++){      
      if(refhist->GetMaximum(i)!=0)
    isocheck->setBinContent(i,refhist->GetBinContent(i)/refhist->GetMaximum());
      else
    isocheck->setBinContent(i,0.);
    } 
  }else
    LogInfo("HLTMonElectronConsumer") << " empty isohisto  " ;
  


}

void HLTMonElectronConsumer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 137 of file HLTMonElectronConsumer.cc.

References DQMStore::book1D(), DQMStore::bookFloat(), DQMStore::get(), MonitorElement::getTH1F(), i, cppFunctionSkipper::operator, DQMStore::rmdir(), DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
  
  DQMStore *dbe = 0;
  dbe = Service < DQMStore > ().operator->();
  
  if (dbe) {
    dbe->setCurrentFolder(dirname_);
    dbe->rmdir(dirname_);
  }
  
  
  if (dbe) {
    dbe->setCurrentFolder(dirname_);
    

    // load pixel MEs
    std::string tmpname = pixeldirname_ + "/total eff";
    LogInfo("HLTMonElectronConsumer") << " reading histo: "  << tmpname;    
    pixeltotal=dbe->get(tmpname);
    TH1F* refhist;
    if(pixeltotal!=0){
      for(int i = 0; i<4; i++){
    LogInfo("HLTMonElectronConsumer") << "loop iteration: "<<i ;
    refhist=pixeltotal->getTH1F();
        LogInfo("HLTMonElectronConsumer") << "retrieving: " <<  pixeldirname_ + "/" + refhist->GetXaxis()->GetBinLabel(i+1) + "eta";
    tmpname = pixeldirname_ + "/" + refhist->GetXaxis()->GetBinLabel(i+1) + "eta";
    pixelhistosEta[i]=dbe->get(tmpname);
    tmpname = pixeldirname_ + "/" + refhist->GetXaxis()->GetBinLabel(i+1) + "et";
    pixelhistosEt[i]=dbe->get(tmpname);
    tmpname = pixeldirname_ + "/" + refhist->GetXaxis()->GetBinLabel(i+1) + "phi";
    pixelhistosPhi[i]=dbe->get(tmpname);
      }
      LogInfo("HLTMonElectronConsumer") << "Et ";
      refhist = pixelhistosEt[0]->getTH1F();
      pixelhistosEtOut[0]  =dbe->book1D("pixel eff et","pixel eff et",refhist->GetNbinsX(),refhist->GetXaxis()->GetXmin(),refhist->GetXaxis()->GetXmax());
      pixelhistosEtOut[1]  =dbe->book1D("track eff et","track eff et",refhist->GetNbinsX(),refhist->GetXaxis()->GetXmin(),refhist->GetXaxis()->GetXmax());
      LogInfo("HLTMonElectronConsumer") << "Eta ";
      refhist = pixelhistosEta[0]->getTH1F();
      pixelhistosEtaOut[0] =dbe->book1D("pixel eff eta","pixel eff eta",refhist->GetNbinsX(),refhist->GetXaxis()->GetXmin(),refhist->GetXaxis()->GetXmax());
      pixelhistosEtaOut[1] =dbe->book1D("track eff eta","track eff eta",refhist->GetNbinsX(),refhist->GetXaxis()->GetXmin(),refhist->GetXaxis()->GetXmax());
      LogInfo("HLTMonElectronConsumer") << "Phi ";
      refhist = pixelhistosPhi[0]->getTH1F();
      pixelhistosPhiOut[0] =dbe->book1D("pixel eff phi","pixel eff phi",refhist->GetNbinsX(),refhist->GetXaxis()->GetXmin(),refhist->GetXaxis()->GetXmax());
      pixelhistosPhiOut[1] =dbe->book1D("track eff phi","track eff phi",refhist->GetNbinsX(),refhist->GetXaxis()->GetXmin(),refhist->GetXaxis()->GetXmax());
    }else
      LogInfo("HLTMonElectronConsumer") << "pixelhisto doesn't exist during beginJob" ;
      
    pixelEff=dbe->bookFloat("total pixelmatch");
    trackEff=dbe->bookFloat("total trackmatch");

    LogInfo("HLTMonElectronConsumer") << "writing: "  << pixelEff->getPathname();


    tmpname = isodirname_ + "/total eff";
    LogInfo("HLTMonElectronConsumer") << " reading histo: "  << tmpname;    
    isototal=dbe->get(tmpname);
    refhist = isototal->getTH1F();
    isocheck = dbe->book1D("consistency check","consistency check",refhist->GetNbinsX(),refhist->GetXaxis()->GetXmin(),refhist->GetXaxis()->GetXmax());
    

  } // end "if(dbe)"
}

void HLTMonElectronConsumer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 203 of file HLTMonElectronConsumer.cc.

                               {

//     std::cout << "HLTMonElectronConsumer: end job...." << std::endl;
 
   if (outputFile_.size() != 0 && dbe)
     dbe->save(outputFile_);
 
   return;
}

Member Data Documentation

DQMStore* HLTMonElectronConsumer::dbe

private

Definition at line 64 of file HLTMonElectronConsumer.h.

std::string HLTMonElectronConsumer::dirname_

private

Definition at line 81 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::isocheck

private

Definition at line 70 of file HLTMonElectronConsumer.h.

std::string HLTMonElectronConsumer::isodirname_

private

Definition at line 83 of file HLTMonElectronConsumer.h.

edm::InputTag HLTMonElectronConsumer::isotag_

private

Definition at line 67 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::isototal

private

Definition at line 69 of file HLTMonElectronConsumer.h.

ofstream HLTMonElectronConsumer::logFile_

private

Definition at line 85 of file HLTMonElectronConsumer.h.

bool HLTMonElectronConsumer::monitorDaemon_

private

Definition at line 84 of file HLTMonElectronConsumer.h.

std::string HLTMonElectronConsumer::outputFile_

private

Definition at line 86 of file HLTMonElectronConsumer.h.

std::string HLTMonElectronConsumer::pixeldirname_

private

Definition at line 82 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixelEff

private

Definition at line 78 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixelhistosEt[4]

private

Definition at line 71 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixelhistosEta[4]

private

Definition at line 72 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixelhistosEtaOut[2]

private

Definition at line 75 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixelhistosEtOut[2]

private

Definition at line 74 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixelhistosPhi[4]

private

Definition at line 73 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixelhistosPhiOut[2]

private

Definition at line 76 of file HLTMonElectronConsumer.h.

edm::InputTag HLTMonElectronConsumer::pixeltag_

private

Definition at line 66 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::pixeltotal

private

Definition at line 77 of file HLTMonElectronConsumer.h.

MonitorElement* HLTMonElectronConsumer::trackEff

private

Definition at line 79 of file HLTMonElectronConsumer.h.