#include <HLTMonPhotonClient.h>

Inheritance diagram for HLTMonPhotonClient:

Public Member Functions
	HLTMonPhotonClient (const edm::ParameterSet &)

	~HLTMonPhotonClient ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

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

virtual void	beginJob ()

virtual void	endJob ()

Private Attributes
MonitorElement *	cumFilterEff

DQMStore *	dbe

std::string	dirname_

MonitorElement *	eventCounter

ofstream	logFile_

bool	monitorDaemon_

std::string	outputFile_

MonitorElement *	pixelhistosEt [4]

MonitorElement *	pixelhistosEta [4]

MonitorElement *	pixelhistosEtaOut [4]

MonitorElement *	pixelhistosEtOut [4]

MonitorElement *	pixelhistosPhi [4]

MonitorElement *	pixelhistosPhiOut [4]

MonitorElement *	relFilterEff

std::string	sourcedirname_

edm::InputTag	sourcetag_

std::vector< edm::InputTag >	theHLTCollectionLabels

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
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 30 of file HLTMonPhotonClient.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file HLTMonPhotonClient.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), LogDebug, NULL, and cppFunctionSkipper::operator.

 {
   
   LogDebug("HLTMonPhotonClient") << "constructor...." ;
   
   logFile_.open("HLTMonPhotonClient.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("HLTMonPhotonClient") << "Photon Trigger Monitoring histograms will be saved to " 
                   << outputFile_ ;
   }
   else {
     outputFile_ = "PhotonDQM.root";
   }
   
   bool disable =
     iConfig.getUntrackedParameter < bool > ("disableROOToutput", false);
   if (disable) {
     outputFile_ = "";
   }
   
   sourcetag_=iConfig.getParameter<edm::InputTag>("SourceTag");
 
   theHLTCollectionLabels = iConfig.getParameter<std::vector<edm::InputTag> >("theHLTCollectionLabels");
   
   dirname_="HLT/HLTMonPhoton/"+iConfig.getParameter<std::string>("@module_label");
   sourcedirname_="HLT/HLTMonPhoton/"+sourcetag_.label();
   
   if (dbe != NULL) {
     dbe->setCurrentFolder(dirname_);
   }
   
 }

HLTMonPhotonClient::~HLTMonPhotonClient ( )

Definition at line 57 of file HLTMonPhotonClient.cc.

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

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 72 of file HLTMonPhotonClient.cc.

References i, and j.

 {
   TH1F* refhisto = eventCounter->getTH1F();
   if(refhisto->GetBinContent(1) != 0){
     //    relFilterEff->setBinContent(6,refhisto->GetBinContent(5)/refhisto->GetBinContent(2));
     //    cumFilterEff->setBinContent(5,refhisto->GetBinContent(5)/refhisto->GetBinContent(1));
     cumFilterEff->setBinContent(4,refhisto->GetBinContent(4)/refhisto->GetBinContent(1));
     cumFilterEff->setBinContent(3,refhisto->GetBinContent(3)/refhisto->GetBinContent(1));
     cumFilterEff->setBinContent(2,refhisto->GetBinContent(2)/refhisto->GetBinContent(1));
   }else{
     relFilterEff->setBinContent(6,0);
     cumFilterEff->setBinContent(5,0);
     cumFilterEff->setBinContent(4,0);
     cumFilterEff->setBinContent(3,0);
     cumFilterEff->setBinContent(2,0);
   }
   for(int i = 0; i<4; i++){
     if(refhisto->GetBinContent(4-i) != 0){
       relFilterEff->setBinContent(5-i,refhisto->GetBinContent(5-i)/refhisto->GetBinContent(4-i));
     }else{
       relFilterEff->setBinContent(5-i,0);
     }
   }
   if(refhisto->GetBinContent(1) != 0){
     relFilterEff->setBinContent(1,refhisto->GetBinContent(1)/refhisto->GetBinContent(1));
     cumFilterEff->setBinContent(2,refhisto->GetBinContent(2)/refhisto->GetBinContent(1));
     cumFilterEff->setBinContent(1,refhisto->GetBinContent(1)/refhisto->GetBinContent(1));
   }else{
     relFilterEff->setBinContent(1,0);
     cumFilterEff->setBinContent(2,0);
     cumFilterEff->setBinContent(1,0);
   }
   
   TH1F* num;
   TH1F* denom;
 
   for(int i=0; i<3; i++){   
     denom = pixelhistosEt[0]->getTH1F();
     num = pixelhistosEt[i+1]->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.);
       }
     }
     denom = pixelhistosEta[0]->getTH1F();
     num = pixelhistosEta[i+1]->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.);
     }
     denom = pixelhistosPhi[0]->getTH1F();
     num = pixelhistosPhi[i+1]->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.);
     }
   }
 
   denom = pixelhistosEt[0]->getTH1F();
   num = pixelhistosEt[3]->getTH1F();
   for(int j=1; j <= pixelhistosEtOut[3]->getNbinsX();j++ ){
     if(denom->GetBinContent(j)!=0){
       pixelhistosEtOut[3]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
     }
     else{
       pixelhistosEtOut[3]->setBinContent(j,0.);
     }
   }
 
   denom = pixelhistosEta[0]->getTH1F();
   num = pixelhistosEta[3]->getTH1F();
   for(int j=1; j <= pixelhistosEtaOut[3]->getNbinsX();j++ ){
     if(denom->GetBinContent(j)!=0){
       pixelhistosEtaOut[3]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
     }
     else{
       pixelhistosEtaOut[3]->setBinContent(j,0.);
     }
   }
 
   denom = pixelhistosPhi[0]->getTH1F();
   num = pixelhistosPhi[3]->getTH1F();
   for(int j=1; j <= pixelhistosPhiOut[3]->getNbinsX();j++ ){
     if(denom->GetBinContent(j)!=0){
       pixelhistosPhiOut[3]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
     }
     else{
       pixelhistosPhiOut[3]->setBinContent(j,0.);
     }
   }   
   
 }

void HLTMonPhotonClient::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 175 of file HLTMonPhotonClient.cc.

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

 {
   
   DQMStore *dbe = 0;
   dbe = Service < DQMStore > ().operator->();
   
   if (dbe) {
     dbe->setCurrentFolder(dirname_);
     dbe->rmdir(dirname_);
   }
     
   if (dbe) {
     dbe->setCurrentFolder(dirname_);
 
     std::string tmpname = "";
       
     for(int i = 0; i<4; i++){
       LogInfo("HLTMonPhotonClient") << "loop iteration: "<<i ;
       tmpname = sourcedirname_ + "/" + theHLTCollectionLabels[i].label() + " Eta Dist";
       pixelhistosEta[i]=dbe->get(tmpname);
       tmpname = sourcedirname_ + "/" + theHLTCollectionLabels[i].label() + " Et Dist";
       pixelhistosEt[i]=dbe->get(tmpname);
       tmpname = sourcedirname_ + "/" + theHLTCollectionLabels[i].label() + " Phi Dist";
       pixelhistosPhi[i]=dbe->get(tmpname);
     }
 
     TH1F* refhist;
     int nBins;
     double xMin;
     double xMax;
 
     refhist = pixelhistosEt[0]->getTH1F();
     nBins = refhist->GetNbinsX();
     xMin = refhist->GetXaxis()->GetXmin();
     xMax = refhist->GetXaxis()->GetXmax();
     pixelhistosEtOut[0]  =dbe->book1D("EcalIsol Eff vs Et","EcalIsol Eff vs Et",nBins,xMin,xMax);
     pixelhistosEtOut[1]  =dbe->book1D("HcalIsol Eff vs Et","HcalIsol Eff vs Et",nBins,xMin,xMax);
     pixelhistosEtOut[2]  =dbe->book1D("Track Isol Eff vs Et","Track Isol Eff vs Et",nBins,xMin,xMax);
     pixelhistosEtOut[3]  =dbe->book1D("Total Eff vs Et","Total Eff vs Et",nBins,xMin,xMax);
 
     refhist = pixelhistosEta[0]->getTH1F();
     nBins = refhist->GetNbinsX();
     xMin = refhist->GetXaxis()->GetXmin();
     xMax = refhist->GetXaxis()->GetXmax();
     pixelhistosEtaOut[0] =dbe->book1D("EcalIsol Eff vs Eta","EcalIsol Eff vs Eta",nBins,xMin,xMax);
     pixelhistosEtaOut[1] =dbe->book1D("HcalIsol Eff vs Eta","HcalIsol Eff vs Eta",nBins,xMin,xMax);
     pixelhistosEtaOut[2] =dbe->book1D("Track Isol Eff vs Eta","Track Isol Eff vs Eta",nBins,xMin,xMax);
     pixelhistosEtaOut[3] =dbe->book1D("Total Eff vs Eta","Total Eff vs Eta",nBins,xMin,xMax);
 
     refhist = pixelhistosPhi[0]->getTH1F();
     nBins = refhist->GetNbinsX();
     xMin = refhist->GetXaxis()->GetXmin();
     xMax = refhist->GetXaxis()->GetXmax();
     pixelhistosPhiOut[0] =dbe->book1D("EcalIsol Eff vs Phi","EcalIsol Eff vs Phi",nBins,xMin,xMax);
     pixelhistosPhiOut[1] =dbe->book1D("HcalIsol Eff vs Phi","HcalIsol Eff vs Phi",nBins,xMin,xMax);
     pixelhistosPhiOut[2] =dbe->book1D("Track Isol Eff vs Phi","Track Isol Eff vs Phi",nBins,xMin,xMax);
     pixelhistosPhiOut[3] =dbe->book1D("Total Eff vs Phi","Total Eff vs Phi",nBins,xMin,xMax);
 
     tmpname = sourcedirname_ + "/Evts Passing Filters";    
     eventCounter = dbe->get(tmpname);
 
     relFilterEff = dbe->book1D("Relative Filter Effs","Relative Filter Effs",5,0,5);
     relFilterEff->setBinLabel(1,"HLTIsoEtFilterEff");
     relFilterEff->setBinLabel(2,"HLTIsoEcalIsolEff");
     relFilterEff->setBinLabel(3,"HLTIsoHcalIsolEff");
     relFilterEff->setBinLabel(4,"HLTIsoTrackIsolEff");
     //    relFilterEff->setBinLabel(5,"FIXME");
     //    relFilterEff->setBinLabel(6,"FIXME");
 
     cumFilterEff = dbe->book1D("Cumulative Filter Effs","Cumulative Filter Effs",5,0,5);
     cumFilterEff->setBinLabel(1,"HLTIsoEtFilterEff");
     cumFilterEff->setBinLabel(2,"HLTIsoEcalIsolEff");
     cumFilterEff->setBinLabel(3,"HLTIsoHcalIsolEff");
     cumFilterEff->setBinLabel(4,"HLTIsoTrackIsolEff");
     //    cumFilterEff->setBinLabel(5,"TrackIsolEff");
 
   } // end "if(dbe)"
 }

void HLTMonPhotonClient::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 256 of file HLTMonPhotonClient.cc.

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

Member Data Documentation

MonitorElement* HLTMonPhotonClient::cumFilterEff

private

Definition at line 55 of file HLTMonPhotonClient.h.

DQMStore* HLTMonPhotonClient::dbe

private

Definition at line 43 of file HLTMonPhotonClient.h.

std::string HLTMonPhotonClient::dirname_

private

Definition at line 57 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::eventCounter

private

Definition at line 53 of file HLTMonPhotonClient.h.

ofstream HLTMonPhotonClient::logFile_

private

Definition at line 60 of file HLTMonPhotonClient.h.

bool HLTMonPhotonClient::monitorDaemon_

private

Definition at line 59 of file HLTMonPhotonClient.h.

std::string HLTMonPhotonClient::outputFile_

private

Definition at line 61 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::pixelhistosEt[4]

private

Definition at line 47 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::pixelhistosEta[4]

private

Definition at line 48 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::pixelhistosEtaOut[4]

private

Definition at line 51 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::pixelhistosEtOut[4]

private

Definition at line 50 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::pixelhistosPhi[4]

private

Definition at line 49 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::pixelhistosPhiOut[4]

private

Definition at line 52 of file HLTMonPhotonClient.h.

MonitorElement* HLTMonPhotonClient::relFilterEff

private

Definition at line 54 of file HLTMonPhotonClient.h.

std::string HLTMonPhotonClient::sourcedirname_

private

Definition at line 58 of file HLTMonPhotonClient.h.

edm::InputTag HLTMonPhotonClient::sourcetag_

private

Definition at line 45 of file HLTMonPhotonClient.h.

std::vector<edm::InputTag> HLTMonPhotonClient::theHLTCollectionLabels

private

Definition at line 62 of file HLTMonPhotonClient.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation