#include <PFJetDQMAnalyzer.h>

Inheritance diagram for PFJetDQMAnalyzer:

Public Member Functions
	PFJetDQMAnalyzer (const edm::ParameterSet &parameterSet)

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
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

void	beginJob ()

void	endJob ()

void	storeBadEvents (edm::Event const &, float &val)

Private Attributes
std::string	benchmarkLabel_

edm::InputTag	inputLabel_

edm::InputTag	matchLabel_

edm::EDGetTokenT< edm::View < reco::Jet > >	myJet_

edm::EDGetTokenT< edm::View < reco::Jet > >	myMatchedJet_

int	nBadEvents_

PFJetMonitor	pfJetMonitor_

edm::ParameterSet	pSet_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

Definition at line 14 of file PFJetDQMAnalyzer.h.

Constructor & Destructor Documentation

PFJetDQMAnalyzer::PFJetDQMAnalyzer ( const edm::ParameterSet & parameterSet )

Definition at line 22 of file PFJetDQMAnalyzer.cc.

References benchmarkLabel_, edm::ParameterSet::getParameter(), inputLabel_, matchLabel_, myJet_, myMatchedJet_, edm::parameterSet(), pfJetMonitor_, pSet_, PFJetMonitor::setParameters(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   pSet_                = parameterSet;
   inputLabel_          = pSet_.getParameter<edm::InputTag>("InputCollection");
   matchLabel_          = pSet_.getParameter<edm::InputTag>("MatchCollection");
   benchmarkLabel_      = pSet_.getParameter<std::string>("BenchmarkLabel"); 
 
   pfJetMonitor_.setParameters(parameterSet);  // set parameters for booking histograms and validating jet 
   
   myJet_ = consumes< edm::View<reco::Jet> >(inputLabel_);
   myMatchedJet_ = consumes< edm::View<reco::Jet> >(matchLabel_);
 }

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 51 of file PFJetDQMAnalyzer.cc.

References PFJetMonitor::fill(), edm::Event::getByToken(), edm::ParameterSet::getParameter(), edm::HandleBase::isValid(), myJet_, myMatchedJet_, nBadEvents_, pfJetMonitor_, pSet_, and storeBadEvents().

                                                 {
 
   edm::Handle< edm::View<reco::Jet> > jetCollection;
   iEvent.getByToken(myJet_, jetCollection);   
   
   edm::Handle< edm::View<reco::Jet> > matchedJetCollection; 
   iEvent.getByToken(myMatchedJet_, matchedJetCollection);
 
   float maxRes = 0.0;
   float minRes = 99.99;
   float jetpT = 0.0;
   if (jetCollection.isValid() && matchedJetCollection.isValid()) {
     //pfJetMonitor_.fill( *jetCollection, *matchedJetCollection, minRes, maxRes);  // match collections and fill pt eta phi and charge histos for candidate jet, fill delta_x_VS_y histos for matched couples, book and fill delta_frac_VS_frac histos for matched couples
     pfJetMonitor_.fill( *jetCollection, *matchedJetCollection, minRes, maxRes, jetpT, pSet_);  // match collections and fill pt eta phi and charge histos for candidate jet, fill delta_x_VS_y histos for matched couples, book and fill delta_frac_VS_frac histos for matched couples
     edm::ParameterSet skimPS = pSet_.getParameter<edm::ParameterSet>("SkimParameter");
     if ( (skimPS.getParameter<bool>("switchOn")) &&  
          (nBadEvents_ <= skimPS.getParameter<int32_t>("maximumNumberToBeStored")) ) {
       if (jetpT > skimPS.getParameter<double>("minimumJetpT")) { 
     if ( minRes < skimPS.getParameter<double>("lowerCutOffOnResolution")) {
       storeBadEvents(iEvent,minRes);
       nBadEvents_++;
     } else if (maxRes > skimPS.getParameter<double>("upperCutOffOnResolution")) {
       storeBadEvents(iEvent,maxRes);
       nBadEvents_++;
     }
       } // minimum jet pT check
     }
   }
 }

void PFJetDQMAnalyzer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 38 of file PFJetDQMAnalyzer.cc.

References benchmarkLabel_, Benchmark::DQM_, nBadEvents_, cppFunctionSkipper::operator, cmsHarvester::path, pfJetMonitor_, pSet_, DQMStore::setCurrentFolder(), PFJetMonitor::setup(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                 {
 
   Benchmark::DQM_ = edm::Service<DQMStore>().operator->();
   // part of the following could be put in the base class
   std::string path = "ParticleFlow/" + benchmarkLabel_;
   Benchmark::DQM_->setCurrentFolder(path.c_str());
   edm::LogInfo("PFJetDQMAnalyzer") << " PFJetDQMAnalyzer::beginJob " << "Histogram Folder path set to "<< path;
   pfJetMonitor_.setup(pSet_);  // booking histograms of type delta_frac_VS_frac from PFJetMonitor, pt_ eta_ phi_ and charge_ from CandidateBenchmark,  delta_x_VS_y from MatchCandidateBenchmark
   nBadEvents_ = 0;
 }

void PFJetDQMAnalyzer::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 98 of file PFJetDQMAnalyzer.cc.

98 {

99 }

void PFJetDQMAnalyzer::storeBadEvents	(	edm::Event const &	iEvent,
		float &	val
	)

private

Definition at line 81 of file PFJetDQMAnalyzer.cc.

References benchmarkLabel_, DQMStore::bookFloat(), Benchmark::DQM_, edm::EventID::event(), MonitorElement::Fill(), DQMStore::get(), edm::EventBase::id(), edm::EventID::luminosityBlock(), cmsHarvester::path, MonitorElement::Reset(), edm::EventID::run(), DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by analyze().

                                                                         {
   unsigned int runNb  = iEvent.id().run();
   unsigned int evtNb  = iEvent.id().event();
   unsigned int lumiNb = iEvent.id().luminosityBlock();
   
   std::string path = "ParticleFlow/" + benchmarkLabel_ + "/BadEvents";
   Benchmark::DQM_->setCurrentFolder(path.c_str());
   std::ostringstream eventid_str;
   eventid_str << runNb << "_"<< evtNb << "_" << lumiNb;
   MonitorElement* me = Benchmark::DQM_->get(path + "/" + eventid_str.str());
   if (me) me->Reset();
   else me = Benchmark::DQM_->bookFloat(eventid_str.str());
   me->Fill(val);  
 }