#include <DQMPFCandidateAnalyzer.h>

Inheritance diagram for DQMPFCandidateAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Get the analysis. More...

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Inizialize parameters for histo binning. More...

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	Initialize run-based parameters. More...

	DQMPFCandidateAnalyzer (const edm::ParameterSet &)
	Constructor. More...

	~DQMPFCandidateAnalyzer () override
	Destructor. More...

Public Member Functions inherited from DQMOneEDAnalyzer<>
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Private Member Functions
static bool	jetSortingRule (reco::Jet x, reco::Jet y)

Private Attributes
bool	bypassAllDCSChecks_

bool	bypassAllPVChecks_

std::string	candidateType_

edm::ParameterSet	cleaningParameters_

std::vector< int >	countsPFCand_

std::vector< int >	countsPFCandRECO_

JetMETDQMDCSFilter *	DCSFilter_

std::vector< edm::ParameterSet >	diagnosticsParameters_

std::vector< double >	etaMaxPFCand_

std::vector< double >	etaMaxPFCandRECO_

std::vector< double >	etaMinPFCand_

std::vector< double >	etaMinPFCandRECO_

HLTConfigProvider	FilterhltConfig_

bool	hbhenoifilterdecision

edm::InputTag	hbheNoiseFilterResultTag_

edm::EDGetTokenT< bool >	hbheNoiseFilterResultToken_

std::string	HBHENoiseStringMiniAOD

double	hcalMin_

bool	isMiniAO_

bool	isMiniAOD_

int	LSBegin_

int	LSEnd_

MonitorElement *	m_HOverTrackP_Barrel_hPt_10_20

MonitorElement *	m_HOverTrackP_Barrel_hPt_1_10

MonitorElement *	m_HOverTrackP_Barrel_hPt_20_50

MonitorElement *	m_HOverTrackP_Barrel_hPt_50

MonitorElement *	m_HOverTrackP_EndCap_hPt_10_20

MonitorElement *	m_HOverTrackP_EndCap_hPt_1_10

MonitorElement *	m_HOverTrackP_EndCap_hPt_20_50

MonitorElement *	m_HOverTrackP_EndCap_hPt_50

MonitorElement *	m_HOverTrackP_trackPtVsEta

MonitorElement *	m_HOverTrackPVsEta_hPt_10_20

MonitorElement *	m_HOverTrackPVsEta_hPt_1_10

MonitorElement *	m_HOverTrackPVsEta_hPt_20_50

MonitorElement *	m_HOverTrackPVsEta_hPt_50

MonitorElement *	m_HOverTrackPVsTrackP_Barrel

MonitorElement *	m_HOverTrackPVsTrackP_EndCap

MonitorElement *	m_HOverTrackPVsTrackPt_Barrel

MonitorElement *	m_HOverTrackPVsTrackPt_EndCap

std::map< std::string, MonitorElement * >	map_of_MEs

edm::InputTag	METFilterMiniAODLabel2_

edm::InputTag	METFilterMiniAODLabel_

edm::EDGetTokenT< edm::TriggerResults >	METFilterMiniAODToken2_

edm::EDGetTokenT< edm::TriggerResults >	METFilterMiniAODToken_

int	miniaodfilterdec

int	miniaodfilterindex

edm::InputTag	mInputCollection_

MonitorElement *	mProfileIsoPFChHad_EcalOccupancyCentral

MonitorElement *	mProfileIsoPFChHad_EcalOccupancyEndcap

MonitorElement *	mProfileIsoPFChHad_EMPtCentral

MonitorElement *	mProfileIsoPFChHad_EMPtEndcap

MonitorElement *	mProfileIsoPFChHad_HadPtCentral

MonitorElement *	mProfileIsoPFChHad_HadPtEndcap

MonitorElement *	mProfileIsoPFChHad_HcalOccupancyCentral

MonitorElement *	mProfileIsoPFChHad_HcalOccupancyEndcap

MonitorElement *	mProfileIsoPFChHad_TrackOccupancy

MonitorElement *	mProfileIsoPFChHad_TrackPt

std::vector< MonitorElement * >	multiplicityPFCand_

std::vector< std::string >	multiplicityPFCand_name_

std::vector< std::string >	multiplicityPFCand_nameRECO_

std::vector< MonitorElement * >	multiplicityPFCandRECO_

int	numPV_

std::vector< MonitorElement * >	occupancyPFCand_

std::vector< std::string >	occupancyPFCand_name_

std::vector< std::string >	occupancyPFCand_name_puppiNolepWeight_

std::vector< std::string >	occupancyPFCand_nameRECO_

std::vector< MonitorElement * >	occupancyPFCand_puppiNolepWeight_

std::vector< MonitorElement * >	occupancyPFCandRECO_

edm::EDGetTokenT< std::vector< pat::PackedCandidate > >	pflowPackedToken_

edm::EDGetTokenT< std::vector< reco::PFCandidate > >	pflowToken_

double	ptMinCand_

std::vector< MonitorElement * >	ptPFCand_

std::vector< std::string >	ptPFCand_name_

std::vector< std::string >	ptPFCand_name_puppiNolepWeight_

std::vector< std::string >	ptPFCand_nameRECO_

std::vector< MonitorElement * >	ptPFCand_puppiNolepWeight_

std::vector< MonitorElement * >	ptPFCandRECO_

edm::InputTag	theTriggerResultsLabel_

std::vector< int >	typePFCand_

std::vector< int >	typePFCandRECO_

int	verbose_

edm::InputTag	vertexTag_

edm::EDGetTokenT< std::vector< reco::Vertex > >	vertexToken_

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer<>
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from DQMOneEDAnalyzer<>
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Protected Attributes inherited from DQMOneEDAnalyzer<>
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

DQM jetMET analysis monitoring for PFCandidates

    Jan. '16: by

    M. Artur Weber

Definition at line 60 of file DQMPFCandidateAnalyzer.h.

Constructor & Destructor Documentation

◆ DQMPFCandidateAnalyzer()

DQMPFCandidateAnalyzer::DQMPFCandidateAnalyzer ( const edm::ParameterSet & pSet )

Constructor.

Definition at line 60 of file DQMPFCandidateAnalyzer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   miniaodfilterdec = -1;
 
   candidateType_ = pSet.getUntrackedParameter<std::string>("CandType");
   //here only choice between miniaod or reco
 
   LSBegin_ = pSet.getParameter<int>("LSBegin");
   LSEnd_ = pSet.getParameter<int>("LSEnd");
 
   isMiniAOD_ = (std::string("Packed") == candidateType_);
 
   mInputCollection_ = pSet.getParameter<edm::InputTag>("PFCandidateLabel");
 
   if (isMiniAOD_) {
     pflowPackedToken_ = consumes<std::vector<pat::PackedCandidate> >(mInputCollection_);
   } else {
     pflowToken_ = consumes<std::vector<reco::PFCandidate> >(mInputCollection_);
   }
 
   miniaodfilterdec = -1;
 
   // Smallest track pt
   ptMinCand_ = pSet.getParameter<double>("ptMinCand");
 
   // Smallest raw HCAL energy linked to the track
   hcalMin_ = pSet.getParameter<double>("hcalMin");
 
   diagnosticsParameters_ = pSet.getParameter<std::vector<edm::ParameterSet> >("METDiagonisticsParameters");
 
   edm::ConsumesCollector iC = consumesCollector();
   //DCS
   DCSFilter_ = new JetMETDQMDCSFilter(pSet.getParameter<ParameterSet>("DCSFilter"), iC);
   if (isMiniAOD_) {
     METFilterMiniAODLabel_ = pSet.getParameter<edm::InputTag>("FilterResultsLabelMiniAOD");
     METFilterMiniAODToken_ = consumes<edm::TriggerResults>(METFilterMiniAODLabel_);
 
     METFilterMiniAODLabel2_ = pSet.getParameter<edm::InputTag>("FilterResultsLabelMiniAOD2");
     METFilterMiniAODToken2_ = consumes<edm::TriggerResults>(METFilterMiniAODLabel2_);
 
     HBHENoiseStringMiniAOD = pSet.getParameter<std::string>("HBHENoiseLabelMiniAOD");
   }
 
   if (!isMiniAOD_) {
     hbheNoiseFilterResultTag_ = pSet.getParameter<edm::InputTag>("HBHENoiseFilterResultLabel");
     hbheNoiseFilterResultToken_ = consumes<bool>(hbheNoiseFilterResultTag_);
   }
   //jet cleanup parameters
   cleaningParameters_ = pSet.getParameter<ParameterSet>("CleaningParameters");
 
   //Vertex requirements
   bypassAllPVChecks_ = cleaningParameters_.getParameter<bool>("bypassAllPVChecks");
   bypassAllDCSChecks_ = cleaningParameters_.getParameter<bool>("bypassAllDCSChecks");
   vertexTag_ = cleaningParameters_.getParameter<edm::InputTag>("vertexCollection");
   vertexToken_ = consumes<std::vector<reco::Vertex> >(edm::InputTag(vertexTag_));
 
   verbose_ = pSet.getParameter<int>("verbose");
 }

◆ ~DQMPFCandidateAnalyzer()

DQMPFCandidateAnalyzer::~DQMPFCandidateAnalyzer ( )

override

Destructor.

Definition at line 122 of file DQMPFCandidateAnalyzer.cc.

References LogTrace.

                                                 {
   delete DCSFilter_;
   LogTrace("DQMPFCandidateAnalyzer") << "[DQMPFCandidateAnalyzer] Saving the histos";
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Get the analysis.

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 666 of file DQMPFCandidateAnalyzer.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), DummyCfis::c, reco::PFBlockElementCluster::clusterRef(), gather_cfg::cout, fed_dqm_sourceclient-live_cfg::DirName, ECAL, bookConverter::elements, l1tnanotables_cff::et, HcalObjRepresent::Fill(), HCAL, mps_fire::i, iEvent, cuy::ii, edm::HandleBase::isValid(), dqmiolumiharvest::j, LogDebug, beamSpotPI::nTracks, pfLinker_cff::particleFlow, edm::Handle< T >::product(), AlCaHLTBitMon_QueryRunRegistry::string, and spclusmultinvestigator_cfi::vertexCollection.

                                                                                         {
   //Vertex information
   Handle<VertexCollection> vertexHandle;
   iEvent.getByToken(vertexToken_, vertexHandle);
 
   if (!vertexHandle.isValid()) {
     LogDebug("") << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
     if (verbose_)
       std::cout << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
   }
   numPV_ = 0;
   if (vertexHandle.isValid()) {
     VertexCollection vertexCollection = *(vertexHandle.product());
     numPV_ = vertexCollection.size();
   }
   bool bPrimaryVertex = (bypassAllPVChecks_ || (numPV_ > 0));
 
   int myLuminosityBlock;
   myLuminosityBlock = iEvent.luminosityBlock();
 
   if (myLuminosityBlock < LSBegin_)
     return;
   if (myLuminosityBlock > LSEnd_ && LSEnd_ > 0)
     return;
 
   if (verbose_)
     std::cout << "METAnalyzer analyze" << std::endl;
 
   std::string DirName = "JetMET/PFCandidates/" + mInputCollection_.label();
 
   bool hbhenoifilterdecision = true;
   if (!isMiniAOD_) {  //not checked for MiniAOD -> for miniaod decision filled as "triggerResults" bool
     edm::Handle<bool> HBHENoiseFilterResultHandle;
     iEvent.getByToken(hbheNoiseFilterResultToken_, HBHENoiseFilterResultHandle);
     if (!HBHENoiseFilterResultHandle.isValid()) {
       LogDebug("") << "METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
       if (verbose_)
         std::cout << "METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
     }
     hbhenoifilterdecision = *HBHENoiseFilterResultHandle;
   } else {  //need to check if we go for version 1 or version 2
     edm::Handle<edm::TriggerResults> metFilterResults;
     iEvent.getByToken(METFilterMiniAODToken_, metFilterResults);
     if (metFilterResults.isValid()) {
       if (miniaodfilterindex != -1) {
         hbhenoifilterdecision = metFilterResults->accept(miniaodfilterindex);
       }
     } else {
       iEvent.getByToken(METFilterMiniAODToken2_, metFilterResults);
       if (metFilterResults.isValid()) {
         if (miniaodfilterindex != -1) {
           hbhenoifilterdecision = metFilterResults->accept(miniaodfilterindex);
         }
       }
     }
   }
 
   //DCS Filter
   bool bDCSFilter = (bypassAllDCSChecks_ || DCSFilter_->filter(iEvent, iSetup));
 
   for (unsigned int i = 0; i < countsPFCand_.size(); i++) {
     countsPFCand_[i] = 0;
   }
   if (bDCSFilter && hbhenoifilterdecision && bPrimaryVertex) {
     if (isMiniAOD_) {
       edm::Handle<std::vector<pat::PackedCandidate> > packedParticleFlow;
       iEvent.getByToken(pflowPackedToken_, packedParticleFlow);
       //11, 13, 22 for el/mu/gamma, 211 chargedHadron, 130 neutralHadrons, 1 and 2 hadHF and EGammaHF
       for (unsigned int i = 0; i < packedParticleFlow->size(); i++) {
         const pat::PackedCandidate& c = packedParticleFlow->at(i);
         for (unsigned int j = 0; j < typePFCand_.size(); j++) {
           if (abs(c.pdgId()) == typePFCand_[j]) {
             //second check for endcap, if inside barrel Max and Min symmetric around 0
             if (((c.eta() > etaMinPFCand_[j]) && (c.eta() < etaMaxPFCand_[j])) ||
                 ((c.eta() > (-etaMaxPFCand_[j])) && (c.eta() < (-etaMinPFCand_[j])))) {
               countsPFCand_[j] += 1;
               ptPFCand_[j] = map_of_MEs[DirName + "/" + ptPFCand_name_[j]];
               if (ptPFCand_[j] && ptPFCand_[j]->getRootObject())
                 ptPFCand_[j]->Fill(c.eta(), c.phi(), c.pt() * c.puppiWeight());
               occupancyPFCand_[j] = map_of_MEs[DirName + "/" + occupancyPFCand_name_[j]];
               if (occupancyPFCand_[j] && occupancyPFCand_[j]->getRootObject())
                 occupancyPFCand_[j]->Fill(c.eta(), c.phi(), c.puppiWeight());
               ptPFCand_puppiNolepWeight_[j] = map_of_MEs[DirName + "/" + ptPFCand_name_puppiNolepWeight_[j]];
               if (ptPFCand_puppiNolepWeight_[j] && ptPFCand_puppiNolepWeight_[j]->getRootObject())
                 ptPFCand_puppiNolepWeight_[j]->Fill(c.eta(), c.phi(), c.pt() * c.puppiWeightNoLep());
               occupancyPFCand_puppiNolepWeight_[j] =
                   map_of_MEs[DirName + "/" + occupancyPFCand_name_puppiNolepWeight_[j]];
               if (occupancyPFCand_puppiNolepWeight_[j] && occupancyPFCand_puppiNolepWeight_[j]->getRootObject()) {
                 occupancyPFCand_puppiNolepWeight_[j]->Fill(c.eta(), c.phi(), c.puppiWeightNoLep());
               }
             }
           }
         }
       }
       for (unsigned int j = 0; j < countsPFCand_.size(); j++) {
         multiplicityPFCand_[j] = map_of_MEs[DirName + "/" + multiplicityPFCand_name_[j]];
         if (multiplicityPFCand_[j] && multiplicityPFCand_[j]->getRootObject()) {
           multiplicityPFCand_[j]->Fill(countsPFCand_[j]);
         }
       }
     } else {
       edm::Handle<std::vector<reco::PFCandidate> > particleFlow;
       iEvent.getByToken(pflowToken_, particleFlow);
       for (unsigned int i = 0; i < particleFlow->size(); i++) {
         const reco::PFCandidate& c = particleFlow->at(i);
         for (unsigned int j = 0; j < typePFCandRECO_.size(); j++) {
           if (c.particleId() == typePFCandRECO_[j]) {
             //second check for endcap, if inside barrel Max and Min symmetric around 0
             if (((c.eta() > etaMinPFCandRECO_[j]) && (c.eta() < etaMaxPFCandRECO_[j])) ||
                 ((c.eta() > (-etaMaxPFCandRECO_[j])) && (c.eta() < (-etaMinPFCandRECO_[j])))) {
               countsPFCandRECO_[j] += 1;
               ptPFCandRECO_[j] = map_of_MEs[DirName + "/" + ptPFCand_nameRECO_[j]];
               if (ptPFCandRECO_[j] && ptPFCandRECO_[j]->getRootObject())
                 ptPFCandRECO_[j]->Fill(c.eta(), c.phi(), c.pt());
               occupancyPFCandRECO_[j] = map_of_MEs[DirName + "/" + occupancyPFCand_nameRECO_[j]];
               if (occupancyPFCandRECO_[j] && occupancyPFCandRECO_[j]->getRootObject())
                 occupancyPFCandRECO_[j]->Fill(c.eta(), c.phi());
             }
             //fill quantities for isolated charged hadron quantities
             //only for charged hadrons
             if (c.particleId() == 1 && c.pt() > ptMinCand_) {
               // At least 1 GeV in HCAL
               double ecalRaw = c.rawEcalEnergy();
               double hcalRaw = c.rawHcalEnergy();
               if ((ecalRaw + hcalRaw) > hcalMin_) {
                 const PFCandidate::ElementsInBlocks& theElements = c.elementsInBlocks();
                 if (theElements.empty())
                   continue;
                 unsigned int iTrack = -999;
                 std::vector<unsigned int> iECAL;  // =999;
                 std::vector<unsigned int> iHCAL;  // =999;
                 const reco::PFBlockRef blockRef = theElements[0].first;
                 const edm::OwnVector<reco::PFBlockElement>& elements = blockRef->elements();
                 // Check that there is only one track in the block.
                 unsigned int nTracks = 0;
                 for (unsigned int iEle = 0; iEle < elements.size(); iEle++) {
                   // Find the tracks in the block
                   PFBlockElement::Type type = elements[iEle].type();
                   switch (type) {
                     case PFBlockElement::TRACK:
                       iTrack = iEle;
                       nTracks++;
                       break;
                     case PFBlockElement::ECAL:
                       iECAL.push_back(iEle);
                       break;
                     case PFBlockElement::HCAL:
                       iHCAL.push_back(iEle);
                       break;
                     default:
                       continue;
                   }
                 }
                 if (nTracks == 1) {
                   // Characteristics of the track
                   const reco::PFBlockElementTrack& et =
                       dynamic_cast<const reco::PFBlockElementTrack&>(elements[iTrack]);
                   mProfileIsoPFChHad_TrackOccupancy = map_of_MEs[DirName + "/" + "IsoPfChHad_Track_profile"];
                   if (mProfileIsoPFChHad_TrackOccupancy && mProfileIsoPFChHad_TrackOccupancy->getRootObject())
                     mProfileIsoPFChHad_TrackOccupancy->Fill(et.trackRef()->eta(), et.trackRef()->phi());
                   mProfileIsoPFChHad_TrackPt = map_of_MEs[DirName + "/" + "IsoPfChHad_TrackPt"];
                   if (mProfileIsoPFChHad_TrackPt && mProfileIsoPFChHad_TrackPt->getRootObject())
                     mProfileIsoPFChHad_TrackPt->Fill(et.trackRef()->eta(), et.trackRef()->phi(), et.trackRef()->pt());
                   if (c.rawEcalEnergy() == 0) {  //isolated hadron, nothing in ECAL
                     //right now take corrected hcalEnergy, do we want the rawHcalEnergy instead
                     m_HOverTrackP_trackPtVsEta = map_of_MEs[DirName + "/" + "HOverTrackP_trackPtVsEta"];
                     if (m_HOverTrackP_trackPtVsEta && m_HOverTrackP_trackPtVsEta->getRootObject())
                       m_HOverTrackP_trackPtVsEta->Fill(c.pt(), c.eta(), c.hcalEnergy() / et.trackRef()->p());
                     if (c.pt() > 1 && c.pt() < 10) {
                       m_HOverTrackPVsEta_hPt_1_10 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_1_10"];
                       if (m_HOverTrackPVsEta_hPt_1_10 && m_HOverTrackPVsEta_hPt_1_10->getRootObject())
                         m_HOverTrackPVsEta_hPt_1_10->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
                     } else if (c.pt() > 10 && c.pt() < 20) {
                       m_HOverTrackPVsEta_hPt_10_20 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_10_20"];
                       if (m_HOverTrackPVsEta_hPt_10_20 && m_HOverTrackPVsEta_hPt_10_20->getRootObject())
                         m_HOverTrackPVsEta_hPt_10_20->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
                     } else if (c.pt() > 20 && c.pt() < 50) {
                       m_HOverTrackPVsEta_hPt_20_50 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_20_50"];
                       if (m_HOverTrackPVsEta_hPt_20_50 && m_HOverTrackPVsEta_hPt_20_50->getRootObject())
                         m_HOverTrackPVsEta_hPt_20_50->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
                     } else if (c.pt() > 50) {
                       m_HOverTrackPVsEta_hPt_50 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_50"];
                       if (m_HOverTrackPVsEta_hPt_50 && m_HOverTrackPVsEta_hPt_50->getRootObject())
                         m_HOverTrackPVsEta_hPt_50->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
                     }
                     if (fabs(c.eta() < 1.392)) {
                       if (c.pt() > 1 && c.pt() < 10) {
                         m_HOverTrackP_Barrel_hPt_1_10 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_1_10"];
                         if (m_HOverTrackP_Barrel_hPt_1_10 && m_HOverTrackP_Barrel_hPt_1_10->getRootObject())
                           m_HOverTrackP_Barrel_hPt_1_10->Fill(c.hcalEnergy() / et.trackRef()->p());
                       } else if (c.pt() > 10 && c.pt() < 20) {
                         m_HOverTrackP_Barrel_hPt_10_20 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_10_20"];
                         if (m_HOverTrackP_Barrel_hPt_10_20 && m_HOverTrackP_Barrel_hPt_10_20->getRootObject())
                           m_HOverTrackP_Barrel_hPt_10_20->Fill(c.hcalEnergy() / et.trackRef()->p());
                       } else if (c.pt() > 20 && c.pt() < 50) {
                         m_HOverTrackP_Barrel_hPt_20_50 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_20_50"];
                         if (m_HOverTrackP_Barrel_hPt_20_50 && m_HOverTrackP_Barrel_hPt_20_50->getRootObject())
                           m_HOverTrackP_Barrel_hPt_20_50->Fill(c.hcalEnergy() / et.trackRef()->p());
                       } else if (c.pt() > 50) {
                         m_HOverTrackP_Barrel_hPt_50 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_50"];
                         if (m_HOverTrackP_Barrel_hPt_50 && m_HOverTrackP_Barrel_hPt_50->getRootObject())
                           m_HOverTrackP_Barrel_hPt_50->Fill(c.hcalEnergy() / et.trackRef()->p());
                       }
                       m_HOverTrackPVsTrackP_Barrel = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackP_Barrel"];
                       if (m_HOverTrackPVsTrackP_Barrel && m_HOverTrackPVsTrackP_Barrel->getRootObject())
                         m_HOverTrackPVsTrackP_Barrel->Fill(et.trackRef()->p(), c.hcalEnergy() / et.trackRef()->p());
                       m_HOverTrackPVsTrackPt_Barrel = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackPt_Barrel"];
                       if (m_HOverTrackPVsTrackPt_Barrel && m_HOverTrackPVsTrackPt_Barrel->getRootObject())
                         m_HOverTrackPVsTrackPt_Barrel->Fill(et.trackRef()->pt(), c.hcalEnergy() / et.trackRef()->p());
                     } else {
                       m_HOverTrackPVsTrackP_EndCap = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackP_EndCap"];
                       if (m_HOverTrackPVsTrackP_EndCap && m_HOverTrackPVsTrackP_EndCap->getRootObject())
                         m_HOverTrackPVsTrackP_EndCap->Fill(et.trackRef()->p(), c.hcalEnergy() / et.trackRef()->p());
                       m_HOverTrackPVsTrackPt_EndCap = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackPt_EndCap"];
                       if (m_HOverTrackPVsTrackPt_EndCap && m_HOverTrackPVsTrackPt_EndCap->getRootObject())
                         m_HOverTrackPVsTrackPt_EndCap->Fill(et.trackRef()->pt(), c.hcalEnergy() / et.trackRef()->p());
                       if (c.pt() > 1 && c.pt() < 10) {
                         m_HOverTrackP_EndCap_hPt_1_10 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_1_10"];
                         if (m_HOverTrackP_EndCap_hPt_1_10 && m_HOverTrackP_EndCap_hPt_1_10->getRootObject())
                           m_HOverTrackP_EndCap_hPt_1_10->Fill(c.hcalEnergy() / et.trackRef()->p());
                       } else if (c.pt() > 10 && c.pt() < 20) {
                         m_HOverTrackP_EndCap_hPt_10_20 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_10_20"];
                         if (m_HOverTrackP_EndCap_hPt_10_20 && m_HOverTrackP_EndCap_hPt_10_20->getRootObject())
                           m_HOverTrackP_EndCap_hPt_10_20->Fill(c.hcalEnergy() / et.trackRef()->p());
                       } else if (c.pt() > 20 && c.pt() < 50) {
                         m_HOverTrackP_EndCap_hPt_20_50 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_20_50"];
                         if (m_HOverTrackP_EndCap_hPt_20_50 && m_HOverTrackP_EndCap_hPt_20_50->getRootObject())
                           m_HOverTrackP_EndCap_hPt_20_50->Fill(c.hcalEnergy() / et.trackRef()->p());
                       } else if (c.pt() > 50) {
                         m_HOverTrackP_EndCap_hPt_50 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_50"];
                         if (m_HOverTrackP_EndCap_hPt_50 && m_HOverTrackP_EndCap_hPt_50->getRootObject())
                           m_HOverTrackP_EndCap_hPt_50->Fill(c.hcalEnergy() / et.trackRef()->p());
                       }
                     }
                   }
                   //ECAL element
                   for (unsigned int ii = 0; ii < iECAL.size(); ii++) {
                     const reco::PFBlockElementCluster& eecal =
                         dynamic_cast<const reco::PFBlockElementCluster&>(elements[iECAL[ii]]);
                     if (fabs(eecal.clusterRef()->eta()) < 1.479) {
                       mProfileIsoPFChHad_EcalOccupancyCentral =
                           map_of_MEs[DirName + "/" + "IsoPfChHad_ECAL_profile_central"];
                       if (mProfileIsoPFChHad_EcalOccupancyCentral &&
                           mProfileIsoPFChHad_EcalOccupancyCentral->getRootObject())
                         mProfileIsoPFChHad_EcalOccupancyCentral->Fill(eecal.clusterRef()->eta(),
                                                                       eecal.clusterRef()->phi());
                       mProfileIsoPFChHad_EMPtCentral = map_of_MEs[DirName + "/" + "IsoPfChHad_EMPt_central"];
                       if (mProfileIsoPFChHad_EMPtCentral && mProfileIsoPFChHad_EMPtCentral->getRootObject())
                         mProfileIsoPFChHad_EMPtCentral->Fill(
                             eecal.clusterRef()->eta(), eecal.clusterRef()->phi(), eecal.clusterRef()->pt());
                     } else {
                       mProfileIsoPFChHad_EcalOccupancyEndcap =
                           map_of_MEs[DirName + "/" + "IsoPfChHad_ECAL_profile_endcap"];
                       if (mProfileIsoPFChHad_EcalOccupancyEndcap &&
                           mProfileIsoPFChHad_EcalOccupancyEndcap->getRootObject())
                         mProfileIsoPFChHad_EcalOccupancyEndcap->Fill(eecal.clusterRef()->eta(),
                                                                      eecal.clusterRef()->phi());
                       mProfileIsoPFChHad_EMPtEndcap = map_of_MEs[DirName + "/" + "IsoPfChHad_EMPt_endcap"];
                       if (mProfileIsoPFChHad_EMPtEndcap && mProfileIsoPFChHad_EMPtEndcap->getRootObject())
                         mProfileIsoPFChHad_EMPtEndcap->Fill(
                             eecal.clusterRef()->eta(), eecal.clusterRef()->phi(), eecal.clusterRef()->pt());
                     }
                   }
                   //HCAL element
                   for (unsigned int ii = 0; ii < iHCAL.size(); ii++) {
                     const reco::PFBlockElementCluster& ehcal =
                         dynamic_cast<const reco::PFBlockElementCluster&>(elements[iHCAL[ii]]);
                     if (fabs(ehcal.clusterRef()->eta()) < 1.740) {
                       mProfileIsoPFChHad_HcalOccupancyCentral =
                           map_of_MEs[DirName + "/" + "IsoPfChHad_HCAL_profile_central"];
                       if (mProfileIsoPFChHad_HcalOccupancyCentral &&
                           mProfileIsoPFChHad_HcalOccupancyCentral->getRootObject())
                         mProfileIsoPFChHad_HcalOccupancyCentral->Fill(ehcal.clusterRef()->eta(),
                                                                       ehcal.clusterRef()->phi());
                       mProfileIsoPFChHad_HadPtCentral = map_of_MEs[DirName + "/" + "IsoPfChHad_HadPt_central"];
                       if (mProfileIsoPFChHad_HadPtCentral && mProfileIsoPFChHad_HadPtCentral->getRootObject())
                         mProfileIsoPFChHad_HadPtCentral->Fill(
                             ehcal.clusterRef()->eta(), ehcal.clusterRef()->phi(), ehcal.clusterRef()->pt());
                     } else {
                       mProfileIsoPFChHad_HcalOccupancyEndcap =
                           map_of_MEs[DirName + "/" + "IsoPfChHad_HCAL_profile_endcap"];
                       if (mProfileIsoPFChHad_HcalOccupancyEndcap &&
                           mProfileIsoPFChHad_HcalOccupancyEndcap->getRootObject())
                         mProfileIsoPFChHad_HcalOccupancyEndcap->Fill(ehcal.clusterRef()->eta(),
                                                                      ehcal.clusterRef()->phi());
                       mProfileIsoPFChHad_HadPtEndcap = map_of_MEs[DirName + "/" + "IsoPfChHad_HadPt_endcap"];
                       if (mProfileIsoPFChHad_HadPtEndcap && mProfileIsoPFChHad_HadPtEndcap->getRootObject())
                         mProfileIsoPFChHad_HadPtEndcap->Fill(
                             ehcal.clusterRef()->eta(), ehcal.clusterRef()->phi(), ehcal.clusterRef()->pt());
                     }
                   }
                 }
               }
             }
           }
         }
       }
       for (unsigned int j = 0; j < countsPFCandRECO_.size(); j++) {
         multiplicityPFCandRECO_[j] = map_of_MEs[DirName + "/" + multiplicityPFCand_nameRECO_[j]];
         if (multiplicityPFCandRECO_[j] && multiplicityPFCandRECO_[j]->getRootObject()) {
           multiplicityPFCandRECO_[j]->Fill(countsPFCandRECO_[j]);
         }
       }
     }  //candidate loop for both miniaod and reco
   }
 }

◆ bookHistograms()

void DQMPFCandidateAnalyzer::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &
	)

overridevirtual

Inizialize parameters for histo binning.

Implements DQMOneEDAnalyzer<>.

Definition at line 128 of file DQMPFCandidateAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), fed_dqm_sourceclient-live_cfg::DirName, JetChargeProducer_cfi::exp, nano_mu_digi_cff::float, mps_fire::i, dqm-mbProfile::log, M_PI, LepHTMonitor_cff::nbins_eta, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and findQualityFiles::v.

                                                                                                               {
   ibooker.setCurrentFolder("JetMET/PFCandidates/" + mInputCollection_.label());
   std::string DirName = "JetMET/PFCandidates/" + mInputCollection_.label();
 
   if (!isMiniAOD_) {
     if (!occupancyPFCandRECO_.empty())
       occupancyPFCandRECO_.clear();
     if (!occupancyPFCand_nameRECO_.empty())
       occupancyPFCand_nameRECO_.clear();
     if (!etaMinPFCandRECO_.empty())
       etaMinPFCandRECO_.clear();
     if (!etaMaxPFCandRECO_.empty())
       etaMaxPFCandRECO_.clear();
     if (!typePFCandRECO_.empty())
       typePFCandRECO_.clear();
     if (!countsPFCandRECO_.empty())
       countsPFCandRECO_.clear();
     if (!ptPFCandRECO_.empty())
       ptPFCandRECO_.clear();
     if (!ptPFCand_nameRECO_.empty())
       ptPFCand_nameRECO_.clear();
     if (!multiplicityPFCandRECO_.empty())
       multiplicityPFCandRECO_.clear();
     if (!multiplicityPFCand_nameRECO_.empty())
       multiplicityPFCand_nameRECO_.clear();
     for (std::vector<edm::ParameterSet>::const_iterator v = diagnosticsParameters_.begin();
          v != diagnosticsParameters_.end();
          v++) {
       int etaNBinsPFCand = v->getParameter<int>("etaNBins");
       double etaMinPFCand = v->getParameter<double>("etaMin");
       double etaMaxPFCand = v->getParameter<double>("etaMax");
       int phiNBinsPFCand = v->getParameter<int>("phiNBins");
       double phiMinPFCand = v->getParameter<double>("phiMin");
       double phiMaxPFCand = v->getParameter<double>("phiMax");
       int nMinPFCand = v->getParameter<int>("nMin");
       int nMaxPFCand = v->getParameter<int>("nMax");
       int nbinsPFCand = v->getParameter<double>("nbins");
       etaMinPFCandRECO_.push_back(etaMinPFCand);
       etaMaxPFCandRECO_.push_back(etaMaxPFCand);
       typePFCandRECO_.push_back(v->getParameter<int>("type"));
       countsPFCandRECO_.push_back(0);
       multiplicityPFCandRECO_.push_back(
           ibooker.book1D(std::string(v->getParameter<std::string>("name")).append("_multiplicity_").c_str(),
                          std::string(v->getParameter<std::string>("name")) + "multiplicity",
                          nbinsPFCand,
                          nMinPFCand,
                          nMaxPFCand));
       multiplicityPFCand_nameRECO_.push_back(
           std::string(v->getParameter<std::string>("name")).append("_multiplicity_"));
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + multiplicityPFCand_nameRECO_[multiplicityPFCand_nameRECO_.size() - 1],
           multiplicityPFCandRECO_[multiplicityPFCandRECO_.size() - 1]));
 
       //push back names first, we need to create histograms with the name and fill it for endcap plots later
       occupancyPFCand_nameRECO_.push_back(std::string(v->getParameter<std::string>("name")).append("_occupancy_"));
 
       ptPFCand_nameRECO_.push_back(std::string(v->getParameter<std::string>("name")).append("_pt_"));
       //special booking for endcap plots, merge plots for eminus and eplus into one plot, using variable binning
       //barrel plots have eta-boundaries on minus and plus side
       //parameters start on minus side
       if (etaMinPFCand * etaMaxPFCand < 0) {  //barrel plots, plot only in barrel region
         occupancyPFCandRECO_.push_back(
             ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_occupancy_").c_str(),
                            std::string(v->getParameter<std::string>("name")) + "occupancy",
                            etaNBinsPFCand,
                            etaMinPFCand,
                            etaMaxPFCand,
                            phiNBinsPFCand,
                            phiMinPFCand,
                            phiMaxPFCand));
         ptPFCandRECO_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_").c_str(),
                                                std::string(v->getParameter<std::string>("name")) + "pt",
                                                etaNBinsPFCand,
                                                etaMinPFCand,
                                                etaMaxPFCand,
                                                phiNBinsPFCand,
                                                phiMinPFCand,
                                                phiMaxPFCand));
       } else {  //endcap or forward plots,
         const int nbins_eta_endcap = 2 * (etaNBinsPFCand + 1);
         double eta_limits_endcap[nbins_eta_endcap];
         for (int i = 0; i < nbins_eta_endcap; i++) {
           if (i < (etaNBinsPFCand + 1)) {
             eta_limits_endcap[i] = etaMinPFCand + i * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
           } else {
             eta_limits_endcap[i] =
                 -etaMaxPFCand + (i - (etaNBinsPFCand + 1)) * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
           }
         }
         TH2F* hist_temp_occup = new TH2F((occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size() - 1]).c_str(),
                                          "occupancy",
                                          nbins_eta_endcap - 1,
                                          eta_limits_endcap,
                                          phiNBinsPFCand,
                                          phiMinPFCand,
                                          phiMaxPFCand);
         occupancyPFCandRECO_.push_back(
             ibooker.book2D(occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size() - 1], hist_temp_occup));
         TH2F* hist_temp_pt = new TH2F((ptPFCand_nameRECO_[ptPFCand_nameRECO_.size() - 1]).c_str(),
                                       "pt",
                                       nbins_eta_endcap - 1,
                                       eta_limits_endcap,
                                       phiNBinsPFCand,
                                       phiMinPFCand,
                                       phiMaxPFCand);
         ptPFCandRECO_.push_back(ibooker.book2D(ptPFCand_nameRECO_[ptPFCand_nameRECO_.size() - 1], hist_temp_pt));
       }
 
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size() - 1],
           occupancyPFCandRECO_[occupancyPFCandRECO_.size() - 1]));
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + ptPFCand_nameRECO_[ptPFCand_nameRECO_.size() - 1], ptPFCandRECO_[ptPFCandRECO_.size() - 1]));
     }
 
     mProfileIsoPFChHad_TrackOccupancy = ibooker.book2D(
         "IsoPfChHad_Track_profile", "Isolated PFChHadron Tracker_occupancy", 108, -2.7, 2.7, 160, -M_PI, M_PI);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_Track_profile",
                                                               mProfileIsoPFChHad_TrackOccupancy));
     mProfileIsoPFChHad_TrackPt =
         ibooker.book2D("IsoPfChHad_TrackPt", "Isolated PFChHadron TrackPt", 108, -2.7, 2.7, 160, -M_PI, M_PI);
     map_of_MEs.insert(
         std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_TrackPt", mProfileIsoPFChHad_TrackPt));
 
     mProfileIsoPFChHad_EcalOccupancyCentral = ibooker.book2D("IsoPfChHad_ECAL_profile_central",
                                                              "IsolatedPFChHa ECAL occupancy (Barrel)",
                                                              180,
                                                              -1.479,
                                                              1.479,
                                                              125,
                                                              -M_PI,
                                                              M_PI);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_ECAL_profile_central",
                                                               mProfileIsoPFChHad_EcalOccupancyCentral));
     mProfileIsoPFChHad_EMPtCentral = ibooker.book2D(
         "IsoPfChHad_EMPt_central", "Isolated PFChHadron HadPt_central", 180, -1.479, 1.479, 360, -M_PI, M_PI);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_EMPt_central",
                                                               mProfileIsoPFChHad_EMPtCentral));
 
     mProfileIsoPFChHad_EcalOccupancyEndcap = ibooker.book2D(
         "IsoPfChHad_ECAL_profile_endcap", "IsolatedPFChHa ECAL occupancy (Endcap)", 110, -2.75, 2.75, 125, -M_PI, M_PI);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_ECAL_profile_endcap",
                                                               mProfileIsoPFChHad_EcalOccupancyEndcap));
     mProfileIsoPFChHad_EMPtEndcap =
         ibooker.book2D("IsoPfChHad_EMPt_endcap", "Isolated PFChHadron EMPt_endcap", 110, -2.75, 2.75, 125, -M_PI, M_PI);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_EMPt_endcap",
                                                               mProfileIsoPFChHad_EMPtEndcap));
 
     const int nbins_eta = 16;
 
     double eta_limits[nbins_eta] = {-2.650,
                                     -2.500,
                                     -2.322,
                                     -2.172,
                                     -2.043,
                                     -1.930,
                                     -1.830,
                                     -1.740,
                                     1.740,
                                     1.830,
                                     1.930,
                                     2.043,
                                     2.172,
                                     2.3122,
                                     2.500,
                                     2.650};
 
     TH2F* hist_temp_HCAL = new TH2F("IsoPfChHad_HCAL_profile_endcap",
                                     "IsolatedPFChHa HCAL occupancy (outer endcap)",
                                     nbins_eta - 1,
                                     eta_limits,
                                     36,
                                     -M_PI,
                                     M_PI);
     TH2F* hist_tempPt_HCAL = (TH2F*)hist_temp_HCAL->Clone("Isolated PFCHHadron HadPt (outer endcap)");
 
     mProfileIsoPFChHad_HcalOccupancyCentral = ibooker.book2D("IsoPfChHad_HCAL_profile_central",
                                                              "IsolatedPFChHa HCAL occupancy (Central Part)",
                                                              40,
                                                              -1.740,
                                                              1.740,
                                                              72,
                                                              -M_PI,
                                                              M_PI);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HCAL_profile_central",
                                                               mProfileIsoPFChHad_HcalOccupancyCentral));
     mProfileIsoPFChHad_HadPtCentral = ibooker.book2D(
         "IsoPfChHad_HadPt_central", "Isolated PFChHadron HadPt_central", 40, -1.740, 1.740, 72, -M_PI, M_PI);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HadPt_central",
                                                               mProfileIsoPFChHad_HadPtCentral));
 
     mProfileIsoPFChHad_HcalOccupancyEndcap = ibooker.book2D("IsoPfChHad_HCAL_profile_endcap", hist_temp_HCAL);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HCAL_profile_endcap",
                                                               mProfileIsoPFChHad_HcalOccupancyEndcap));
     mProfileIsoPFChHad_HadPtEndcap = ibooker.book2D("IsoPfChHad_HadPt_endcap", hist_tempPt_HCAL);
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HadPt_endcap",
                                                               mProfileIsoPFChHad_HadPtEndcap));
 
     //actual HCAL segmentation in pseudorapidity -> reduce by a factor of two
     //const int nbins_eta_hcal_total=54;
     //double eta_limits_hcal_total[nbins_eta_hcal_total]=
     //   {-2.650,-2.500,-2.322,-2.172,-2.043,-1.930,-1.830,-1.740,-1.653,-1.566,-1.479,-1.392,-1.305,
     //   -1.218,-1.131,-1.044,-0.957,-0.870,-0.783,-0.696,-0.609,-0.522,-0.435,-0.348,-0.261,-0.174,-0.087,0.0,
     //   0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609, 0.783, 0.870, 0.957, 1.044, 1.131, 1.218
     //   1.305, 1.392, 1.479, 1.566, 1.653, 1.740, 1.830, 1.930, 2.043, 2.172, 2.322, 2.500, 2.650}
     //
 
     const int nbins_eta_hcal_total = 28;
     double eta_limits_hcal_total[nbins_eta_hcal_total] = {
         -2.650, -2.322, -2.043, -1.830, -1.653, -1.479, -1.305, -1.131, -0.957, -0.783, -0.609, -0.435, -0.261, -0.087,
         0.087,  0.261,  0.435,  0.609,  0.783,  0.957,  1.131,  1.305,  1.479,  1.653,  1.830,  2.043,  2.322,  2.650};
     float eta_limits_hcal_total_f[nbins_eta_hcal_total];
     float log_bin_spacing = log(200.) / 40.;
     const int nbins_pt_total_hcal = 41;
     double pt_limits_hcal[nbins_pt_total_hcal];
     float pt_limits_hcal_f[nbins_pt_total_hcal];
     for (int i = 0; i < nbins_pt_total_hcal; i++) {
       pt_limits_hcal[i] = exp(i * log_bin_spacing);
       pt_limits_hcal_f[i] = exp(i * log_bin_spacing);
     }
     for (int i = 0; i < nbins_eta_hcal_total; i++) {
       eta_limits_hcal_total_f[i] = (float)eta_limits_hcal_total[i];
     }
     m_HOverTrackP_trackPtVsEta = ibooker.book2D("HOverTrackP_trackPtVsEta",
                                                 "HOverTrackP_trackPtVsEta",
                                                 nbins_pt_total_hcal - 1,
                                                 pt_limits_hcal_f,
                                                 nbins_eta_hcal_total - 1,
                                                 eta_limits_hcal_total_f);
     m_HOverTrackPVsTrackP_Barrel = ibooker.bookProfile(
         "HOverTrackPVsTrackP_Barrel", "HOverTrackPVsTrackP_Barrel", nbins_pt_total_hcal - 1, pt_limits_hcal, 0, 4, " ");
     m_HOverTrackPVsTrackP_EndCap = ibooker.bookProfile(
         "HOverTrackPVsTrackP_EndCap", "HOverTrackPVsTrackP_EndCap", nbins_pt_total_hcal - 1, pt_limits_hcal, 0, 4, " ");
     m_HOverTrackPVsTrackPt_Barrel = ibooker.bookProfile("HOverTrackPVsTrackPt_Barrel",
                                                         "HOverTrackPVsTrackPt_Barrel",
                                                         nbins_pt_total_hcal - 1,
                                                         pt_limits_hcal,
                                                         0,
                                                         4,
                                                         " ");
     m_HOverTrackPVsTrackPt_EndCap = ibooker.bookProfile("HOverTrackPVsTrackPt_EndCap",
                                                         "HOverTrackPVsTrackPt_EndCap",
                                                         nbins_pt_total_hcal - 1,
                                                         pt_limits_hcal,
                                                         0,
                                                         4,
                                                         " ");
 
     m_HOverTrackPVsEta_hPt_1_10 = ibooker.bookProfile("HOverTrackPVsEta_hPt_1_10",
                                                       "HOverTrackPVsEta, 1<hPt<10 GeV",
                                                       nbins_eta_hcal_total - 1,
                                                       eta_limits_hcal_total,
                                                       0,
                                                       4,
                                                       " ");
     m_HOverTrackPVsEta_hPt_10_20 = ibooker.bookProfile("HOverTrackPVsEta_hPt_10_20",
                                                        "HOverTrackPVsEta, 10<hPt<20 GeV",
                                                        nbins_eta_hcal_total - 1,
                                                        eta_limits_hcal_total,
                                                        0,
                                                        4,
                                                        " ");
     m_HOverTrackPVsEta_hPt_20_50 = ibooker.bookProfile("HOverTrackPVsEta_hPt_20_50",
                                                        "HOverTrackPVsEta, 20<hPt<50 GeV",
                                                        nbins_eta_hcal_total - 1,
                                                        eta_limits_hcal_total,
                                                        0,
                                                        4,
                                                        " ");
     m_HOverTrackPVsEta_hPt_50 = ibooker.bookProfile("HOverTrackPVsEta_hPt_50",
                                                     "HOverTrackPVsEta, hPt>50 GeV",
                                                     nbins_eta_hcal_total - 1,
                                                     eta_limits_hcal_total,
                                                     0,
                                                     4,
                                                     " ");
 
     m_HOverTrackP_Barrel_hPt_1_10 =
         ibooker.book1D("HOverTrackP_Barrel_hPt_1_10", "HOverTrackP_B, 1<hPt<10 GeV", 50, 0, 4);
     m_HOverTrackP_Barrel_hPt_10_20 =
         ibooker.book1D("HOverTrackP_Barrel_hPt_10_20", "HOverTrackP_B, 10<hPt<20 GeV", 50, 0, 4);
     m_HOverTrackP_Barrel_hPt_20_50 =
         ibooker.book1D("HOverTrackP_Barrel_hPt_20_50", "HOverTrackP_B, 20<hPt<50 GeV", 50, 0, 4);
     m_HOverTrackP_Barrel_hPt_50 = ibooker.book1D("HOverTrackP_Barrel_hPt_50", "HOverTrackP_B, hPt>50 GeV", 50, 0, 4);
 
     m_HOverTrackP_EndCap_hPt_1_10 =
         ibooker.book1D("HOverTrackP_EndCap_hPt_1_10", "HOverTrackP_E, 1<hPt<10 GeV", 50, 0, 4);
     m_HOverTrackP_EndCap_hPt_10_20 =
         ibooker.book1D("HOverTrackP_EndCap_hPt_10_20", "HOverTrackP_E, 10<hPt<20 GeV", 50, 0, 4);
     m_HOverTrackP_EndCap_hPt_20_50 =
         ibooker.book1D("HOverTrackP_EndCap_hPt_20_50", "HOverTrackP_E, 20<hPt<50 GeV", 50, 0, 4);
     m_HOverTrackP_EndCap_hPt_50 = ibooker.book1D("HOverTrackP_EndCap_hPt_50", "HOverTrackP_E, hPt>50 GeV", 50, 0, 4);
 
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_trackPtVsEta",
                                                               m_HOverTrackP_trackPtVsEta));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackP_Barrel",
                                                               m_HOverTrackPVsTrackP_Barrel));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackP_EndCap",
                                                               m_HOverTrackPVsTrackP_EndCap));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackPt_Barrel",
                                                               m_HOverTrackPVsTrackPt_Barrel));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackPt_EndCap",
                                                               m_HOverTrackPVsTrackPt_EndCap));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_1_10",
                                                               m_HOverTrackPVsEta_hPt_1_10));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_10_20",
                                                               m_HOverTrackPVsEta_hPt_10_20));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_20_50",
                                                               m_HOverTrackPVsEta_hPt_20_50));
     map_of_MEs.insert(
         std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_50", m_HOverTrackPVsEta_hPt_50));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_1_10",
                                                               m_HOverTrackP_Barrel_hPt_1_10));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_10_20",
                                                               m_HOverTrackP_Barrel_hPt_10_20));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_20_50",
                                                               m_HOverTrackP_Barrel_hPt_20_50));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_50",
                                                               m_HOverTrackP_Barrel_hPt_50));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_1_10",
                                                               m_HOverTrackP_EndCap_hPt_1_10));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_10_20",
                                                               m_HOverTrackP_EndCap_hPt_10_20));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_20_50",
                                                               m_HOverTrackP_EndCap_hPt_20_50));
     map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_50",
                                                               m_HOverTrackP_EndCap_hPt_50));
   } else {  //MiniAOD workflow
     if (!occupancyPFCand_.empty())
       occupancyPFCand_.clear();
     if (!occupancyPFCand_name_.empty())
       occupancyPFCand_name_.clear();
     if (!occupancyPFCand_puppiNolepWeight_.empty())
       occupancyPFCand_puppiNolepWeight_.clear();
     if (!occupancyPFCand_name_puppiNolepWeight_.empty())
       occupancyPFCand_name_puppiNolepWeight_.clear();
     if (!etaMinPFCand_.empty())
       etaMinPFCand_.clear();
     if (!etaMaxPFCand_.empty())
       etaMaxPFCand_.clear();
     if (!typePFCand_.empty())
       typePFCand_.clear();
     if (!countsPFCand_.empty())
       countsPFCand_.clear();
     if (!ptPFCand_.empty())
       ptPFCand_.clear();
     if (!ptPFCand_name_.empty())
       ptPFCand_name_.clear();
     if (!ptPFCand_puppiNolepWeight_.empty())
       ptPFCand_puppiNolepWeight_.clear();
     if (!ptPFCand_name_puppiNolepWeight_.empty())
       ptPFCand_name_puppiNolepWeight_.clear();
     if (!multiplicityPFCand_.empty())
       multiplicityPFCand_.clear();
     if (!multiplicityPFCand_name_.empty())
       multiplicityPFCand_name_.clear();
     for (std::vector<edm::ParameterSet>::const_iterator v = diagnosticsParameters_.begin();
          v != diagnosticsParameters_.end();
          v++) {
       int etaNBinsPFCand = v->getParameter<int>("etaNBins");
       double etaMinPFCand = v->getParameter<double>("etaMin");
       double etaMaxPFCand = v->getParameter<double>("etaMax");
       int phiNBinsPFCand = v->getParameter<int>("phiNBins");
       double phiMinPFCand = v->getParameter<double>("phiMin");
       double phiMaxPFCand = v->getParameter<double>("phiMax");
       int nMinPFCand = v->getParameter<int>("nMin");
       int nMaxPFCand = v->getParameter<int>("nMax");
       int nbinsPFCand = v->getParameter<double>("nbins");
 
       // etaNBins_.push_back(etaNBins);
       etaMinPFCand_.push_back(etaMinPFCand);
       etaMaxPFCand_.push_back(etaMaxPFCand);
       typePFCand_.push_back(v->getParameter<int>("type"));
       countsPFCand_.push_back(0);
 
       multiplicityPFCand_.push_back(
           ibooker.book1D(std::string(v->getParameter<std::string>("name")).append("_multiplicity_").c_str(),
                          std::string(v->getParameter<std::string>("name")) + "multiplicity",
                          nbinsPFCand,
                          nMinPFCand,
                          nMaxPFCand));
       multiplicityPFCand_name_.push_back(std::string(v->getParameter<std::string>("name")).append("_multiplicity_"));
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + multiplicityPFCand_name_[multiplicityPFCand_name_.size() - 1],
           multiplicityPFCand_[multiplicityPFCand_.size() - 1]));
       //push back names first, we need to create histograms with the name and fill it for endcap plots later
       occupancyPFCand_name_.push_back(
           std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiWeight_"));
       ptPFCand_name_.push_back(std::string(v->getParameter<std::string>("name")).append("_pt_puppiWeight_"));
       //push back names first, we need to create histograms with the name and fill it for endcap plots later
       occupancyPFCand_name_puppiNolepWeight_.push_back(
           std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiNolepWeight_"));
       ptPFCand_name_puppiNolepWeight_.push_back(
           std::string(v->getParameter<std::string>("name")).append("_pt_puppiNolepWeight_"));
       //special booking for endcap plots, merge plots for eminus and eplus into one plot, using variable binning
       //barrel plots have eta-boundaries on minus and plus side
       //parameters start on minus side
       if (etaMinPFCand * etaMaxPFCand < 0) {  //barrel plots, plot only in barrel region
         occupancyPFCand_.push_back(
             ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiWeight_").c_str(),
                            std::string(v->getParameter<std::string>("name")) + "occupancy_puppiWeight_",
                            etaNBinsPFCand,
                            etaMinPFCand,
                            etaMaxPFCand,
                            phiNBinsPFCand,
                            phiMinPFCand,
                            phiMaxPFCand));
         ptPFCand_.push_back(
             ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_puppiWeight_").c_str(),
                            std::string(v->getParameter<std::string>("name")) + "pt_puppiWeight_",
                            etaNBinsPFCand,
                            etaMinPFCand,
                            etaMaxPFCand,
                            phiNBinsPFCand,
                            phiMinPFCand,
                            phiMaxPFCand));
         occupancyPFCand_puppiNolepWeight_.push_back(ibooker.book2D(
             std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiNolepWeight_").c_str(),
             std::string(v->getParameter<std::string>("name")) + "occupancy_puppiNolepWeight",
             etaNBinsPFCand,
             etaMinPFCand,
             etaMaxPFCand,
             phiNBinsPFCand,
             phiMinPFCand,
             phiMaxPFCand));
         ptPFCand_puppiNolepWeight_.push_back(
             ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_puppiNolepWeight_").c_str(),
                            std::string(v->getParameter<std::string>("name")) + "pt_puppiNolepWeight",
                            etaNBinsPFCand,
                            etaMinPFCand,
                            etaMaxPFCand,
                            phiNBinsPFCand,
                            phiMinPFCand,
                            phiMaxPFCand));
       } else {  //endcap or forward plots,
         const int nbins_eta_endcap = 2 * (etaNBinsPFCand + 1);
         double eta_limits_endcap[nbins_eta_endcap];
         for (int i = 0; i < nbins_eta_endcap; i++) {
           if (i < (etaNBinsPFCand + 1)) {
             eta_limits_endcap[i] = etaMinPFCand + i * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
           } else {
             eta_limits_endcap[i] =
                 -etaMaxPFCand + (i - (etaNBinsPFCand + 1)) * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
           }
         }
         TH2F* hist_temp_occup = new TH2F((occupancyPFCand_name_[occupancyPFCand_name_.size() - 1]).c_str(),
                                          (occupancyPFCand_name_[occupancyPFCand_name_.size() - 1]).c_str(),
                                          nbins_eta_endcap - 1,
                                          eta_limits_endcap,
                                          phiNBinsPFCand,
                                          phiMinPFCand,
                                          phiMaxPFCand);
         occupancyPFCand_.push_back(
             ibooker.book2D(occupancyPFCand_name_[occupancyPFCand_name_.size() - 1], hist_temp_occup));
         TH2F* hist_temp_pt = new TH2F((ptPFCand_name_[ptPFCand_name_.size() - 1]).c_str(),
                                       (ptPFCand_name_[ptPFCand_name_.size() - 1]).c_str(),
                                       nbins_eta_endcap - 1,
                                       eta_limits_endcap,
                                       phiNBinsPFCand,
                                       phiMinPFCand,
                                       phiMaxPFCand);
         ptPFCand_.push_back(ibooker.book2D(ptPFCand_name_[ptPFCand_name_.size() - 1], hist_temp_pt));
         TH2F* hist_temp_occup_puppiNolepWeight = new TH2F(
             (occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
             (occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
             nbins_eta_endcap - 1,
             eta_limits_endcap,
             phiNBinsPFCand,
             phiMinPFCand,
             phiMaxPFCand);
         occupancyPFCand_puppiNolepWeight_.push_back(
             ibooker.book2D(occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1],
                            hist_temp_occup_puppiNolepWeight));
         TH2F* hist_temp_pt_puppiNolepWeight =
             new TH2F((ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
                      (ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
                      nbins_eta_endcap - 1,
                      eta_limits_endcap,
                      phiNBinsPFCand,
                      phiMinPFCand,
                      phiMaxPFCand);
         ptPFCand_puppiNolepWeight_.push_back(
             ibooker.book2D(ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1],
                            hist_temp_pt_puppiNolepWeight));
       }
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1],
           occupancyPFCand_puppiNolepWeight_[occupancyPFCand_puppiNolepWeight_.size() - 1]));
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1],
           ptPFCand_puppiNolepWeight_[ptPFCand_puppiNolepWeight_.size() - 1]));
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + occupancyPFCand_name_[occupancyPFCand_name_.size() - 1],
           occupancyPFCand_[occupancyPFCand_.size() - 1]));
       map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
           DirName + "/" + ptPFCand_name_[ptPFCand_name_.size() - 1], ptPFCand_[ptPFCand_.size() - 1]));
     }
   }
 }

◆ dqmBeginRun()

void DQMPFCandidateAnalyzer::dqmBeginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Initialize run-based parameters.

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 629 of file DQMPFCandidateAnalyzer.cc.

References newFWLiteAna::found, mps_fire::i, cond::persistency::search(), electronDataDiscovery::search2(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                         {
   miniaodfilterindex = -1;
 
   if (isMiniAOD_) {
     bool changed_filter = true;
     if (FilterhltConfig_.init(iRun, iSetup, METFilterMiniAODLabel_.process(), changed_filter)) {
       miniaodfilterdec = 0;
       for (unsigned int i = 0; i < FilterhltConfig_.size(); i++) {
         std::string search = FilterhltConfig_.triggerName(i).substr(
             5);  //actual label of filter, the first 5 items are Flag_, so stripped off
         std::string search2 =
             HBHENoiseStringMiniAOD;  //all filters end with DQM, which is not in the flag --> ONLY not for HBHEFilters
         std::size_t found = search2.find(search);
         if (found != std::string::npos) {
           miniaodfilterindex = i;
         }
       }
     } else if (FilterhltConfig_.init(iRun, iSetup, METFilterMiniAODLabel2_.process(), changed_filter)) {
       miniaodfilterdec = 1;
       for (unsigned int i = 0; i < FilterhltConfig_.size(); i++) {
         std::string search = FilterhltConfig_.triggerName(i).substr(
             5);  //actual label of filter, the first 5 items are Flag_, so stripped off
         std::string search2 =
             HBHENoiseStringMiniAOD;  //all filters end with DQM, which is not in the flag --> ONLY not for HBHEFilters
         std::size_t found = search2.find(search);
         if (found != std::string::npos) {
           miniaodfilterindex = i;
         }
       }
     } else {
       edm::LogWarning("MiniAOD MET Filter HLT OBject version")
           << "nothing found with both RECO and reRECO label" << std::endl;
     }
   }
 }

◆ jetSortingRule()

static bool DQMPFCandidateAnalyzer::jetSortingRule	(	reco::Jet	x,
		reco::Jet	y
	)

inlinestaticprivate

Definition at line 79 of file DQMPFCandidateAnalyzer.h.

References x, and y.

79 { return x.pt() > y.pt(); }

DDAxes::y

DDAxes::x

Member Data Documentation

◆ bypassAllDCSChecks_

bool DQMPFCandidateAnalyzer::bypassAllDCSChecks_

private

Definition at line 175 of file DQMPFCandidateAnalyzer.h.

◆ bypassAllPVChecks_

bool DQMPFCandidateAnalyzer::bypassAllPVChecks_

private

Definition at line 174 of file DQMPFCandidateAnalyzer.h.

◆ candidateType_

std::string DQMPFCandidateAnalyzer::candidateType_

private

Definition at line 106 of file DQMPFCandidateAnalyzer.h.

◆ cleaningParameters_

edm::ParameterSet DQMPFCandidateAnalyzer::cleaningParameters_

private

Definition at line 113 of file DQMPFCandidateAnalyzer.h.

◆ countsPFCand_

std::vector<int> DQMPFCandidateAnalyzer::countsPFCand_

private

Definition at line 160 of file DQMPFCandidateAnalyzer.h.

◆ countsPFCandRECO_

std::vector<int> DQMPFCandidateAnalyzer::countsPFCandRECO_

private

Definition at line 166 of file DQMPFCandidateAnalyzer.h.

◆ DCSFilter_

JetMETDQMDCSFilter* DQMPFCandidateAnalyzer::DCSFilter_

private

Definition at line 111 of file DQMPFCandidateAnalyzer.h.

◆ diagnosticsParameters_

std::vector<edm::ParameterSet> DQMPFCandidateAnalyzer::diagnosticsParameters_

private

Definition at line 114 of file DQMPFCandidateAnalyzer.h.

◆ etaMaxPFCand_

std::vector<double> DQMPFCandidateAnalyzer::etaMaxPFCand_

private

Definition at line 159 of file DQMPFCandidateAnalyzer.h.

◆ etaMaxPFCandRECO_

std::vector<double> DQMPFCandidateAnalyzer::etaMaxPFCandRECO_

private

Definition at line 165 of file DQMPFCandidateAnalyzer.h.

◆ etaMinPFCand_

std::vector<double> DQMPFCandidateAnalyzer::etaMinPFCand_

private

Definition at line 159 of file DQMPFCandidateAnalyzer.h.

◆ etaMinPFCandRECO_

std::vector<double> DQMPFCandidateAnalyzer::etaMinPFCandRECO_

private

Definition at line 165 of file DQMPFCandidateAnalyzer.h.

◆ FilterhltConfig_

HLTConfigProvider DQMPFCandidateAnalyzer::FilterhltConfig_

private

Definition at line 92 of file DQMPFCandidateAnalyzer.h.

◆ hbhenoifilterdecision

bool DQMPFCandidateAnalyzer::hbhenoifilterdecision

private

Definition at line 98 of file DQMPFCandidateAnalyzer.h.

◆ hbheNoiseFilterResultTag_

edm::InputTag DQMPFCandidateAnalyzer::hbheNoiseFilterResultTag_

private

Definition at line 88 of file DQMPFCandidateAnalyzer.h.

◆ hbheNoiseFilterResultToken_

edm::EDGetTokenT<bool> DQMPFCandidateAnalyzer::hbheNoiseFilterResultToken_

private

Definition at line 89 of file DQMPFCandidateAnalyzer.h.

◆ HBHENoiseStringMiniAOD

std::string DQMPFCandidateAnalyzer::HBHENoiseStringMiniAOD

private

Definition at line 90 of file DQMPFCandidateAnalyzer.h.

◆ hcalMin_

double DQMPFCandidateAnalyzer::hcalMin_

private

Definition at line 118 of file DQMPFCandidateAnalyzer.h.

◆ isMiniAO_

bool DQMPFCandidateAnalyzer::isMiniAO_

private

Definition at line 108 of file DQMPFCandidateAnalyzer.h.

◆ isMiniAOD_

bool DQMPFCandidateAnalyzer::isMiniAOD_

private

Definition at line 179 of file DQMPFCandidateAnalyzer.h.

◆ LSBegin_

int DQMPFCandidateAnalyzer::LSBegin_

private

Definition at line 171 of file DQMPFCandidateAnalyzer.h.

◆ LSEnd_

int DQMPFCandidateAnalyzer::LSEnd_

private

Definition at line 172 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_Barrel_hPt_10_20

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_10_20

private

Definition at line 133 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_Barrel_hPt_1_10

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_1_10

private

Definition at line 132 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_Barrel_hPt_20_50

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_20_50

private

Definition at line 134 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_Barrel_hPt_50

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_50

private

Definition at line 135 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_EndCap_hPt_10_20

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_10_20

private

Definition at line 138 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_EndCap_hPt_1_10

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_1_10

private

Definition at line 137 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_EndCap_hPt_20_50

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_20_50

private

Definition at line 139 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_EndCap_hPt_50

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_50

private

Definition at line 140 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackP_trackPtVsEta

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackP_trackPtVsEta

private

Definition at line 120 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsEta_hPt_10_20

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_10_20

private

Definition at line 128 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsEta_hPt_1_10

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_1_10

private

Definition at line 127 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsEta_hPt_20_50

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_20_50

private

Definition at line 129 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsEta_hPt_50

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_50

private

Definition at line 130 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsTrackP_Barrel

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackP_Barrel

private

Definition at line 121 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsTrackP_EndCap

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackP_EndCap

private

Definition at line 122 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsTrackPt_Barrel

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackPt_Barrel

private

Definition at line 124 of file DQMPFCandidateAnalyzer.h.

◆ m_HOverTrackPVsTrackPt_EndCap

MonitorElement* DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackPt_EndCap

private

Definition at line 125 of file DQMPFCandidateAnalyzer.h.

◆ map_of_MEs

std::map<std::string, MonitorElement*> DQMPFCandidateAnalyzer::map_of_MEs

private

Definition at line 177 of file DQMPFCandidateAnalyzer.h.

◆ METFilterMiniAODLabel2_

edm::InputTag DQMPFCandidateAnalyzer::METFilterMiniAODLabel2_

private

Definition at line 95 of file DQMPFCandidateAnalyzer.h.

◆ METFilterMiniAODLabel_

edm::InputTag DQMPFCandidateAnalyzer::METFilterMiniAODLabel_

private

Definition at line 93 of file DQMPFCandidateAnalyzer.h.

◆ METFilterMiniAODToken2_

edm::EDGetTokenT<edm::TriggerResults> DQMPFCandidateAnalyzer::METFilterMiniAODToken2_

private

Definition at line 96 of file DQMPFCandidateAnalyzer.h.

◆ METFilterMiniAODToken_

edm::EDGetTokenT<edm::TriggerResults> DQMPFCandidateAnalyzer::METFilterMiniAODToken_

private

Definition at line 94 of file DQMPFCandidateAnalyzer.h.

◆ miniaodfilterdec

int DQMPFCandidateAnalyzer::miniaodfilterdec

private

Definition at line 101 of file DQMPFCandidateAnalyzer.h.

◆ miniaodfilterindex

int DQMPFCandidateAnalyzer::miniaodfilterindex

private

Definition at line 99 of file DQMPFCandidateAnalyzer.h.

◆ mInputCollection_

edm::InputTag DQMPFCandidateAnalyzer::mInputCollection_

private

Definition at line 86 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_EcalOccupancyCentral

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EcalOccupancyCentral

private

Definition at line 150 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_EcalOccupancyEndcap

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EcalOccupancyEndcap

private

Definition at line 151 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_EMPtCentral

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EMPtCentral

private

Definition at line 144 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_EMPtEndcap

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EMPtEndcap

private

Definition at line 145 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_HadPtCentral

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HadPtCentral

private

Definition at line 142 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_HadPtEndcap

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HadPtEndcap

private

Definition at line 143 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_HcalOccupancyCentral

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HcalOccupancyCentral

private

Definition at line 148 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_HcalOccupancyEndcap

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HcalOccupancyEndcap

private

Definition at line 149 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_TrackOccupancy

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_TrackOccupancy

private

Definition at line 152 of file DQMPFCandidateAnalyzer.h.

◆ mProfileIsoPFChHad_TrackPt

MonitorElement* DQMPFCandidateAnalyzer::mProfileIsoPFChHad_TrackPt

private

Definition at line 146 of file DQMPFCandidateAnalyzer.h.

◆ multiplicityPFCand_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::multiplicityPFCand_

private

Definition at line 155 of file DQMPFCandidateAnalyzer.h.

◆ multiplicityPFCand_name_

std::vector<std::string> DQMPFCandidateAnalyzer::multiplicityPFCand_name_

private

Definition at line 156 of file DQMPFCandidateAnalyzer.h.

◆ multiplicityPFCand_nameRECO_

std::vector<std::string> DQMPFCandidateAnalyzer::multiplicityPFCand_nameRECO_

private

Definition at line 164 of file DQMPFCandidateAnalyzer.h.

◆ multiplicityPFCandRECO_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::multiplicityPFCandRECO_

private

Definition at line 163 of file DQMPFCandidateAnalyzer.h.

◆ numPV_

int DQMPFCandidateAnalyzer::numPV_

private

Definition at line 168 of file DQMPFCandidateAnalyzer.h.

◆ occupancyPFCand_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::occupancyPFCand_

private

Definition at line 155 of file DQMPFCandidateAnalyzer.h.

◆ occupancyPFCand_name_

std::vector<std::string> DQMPFCandidateAnalyzer::occupancyPFCand_name_

private

Definition at line 156 of file DQMPFCandidateAnalyzer.h.

◆ occupancyPFCand_name_puppiNolepWeight_

std::vector<std::string> DQMPFCandidateAnalyzer::occupancyPFCand_name_puppiNolepWeight_

private

Definition at line 158 of file DQMPFCandidateAnalyzer.h.

◆ occupancyPFCand_nameRECO_

std::vector<std::string> DQMPFCandidateAnalyzer::occupancyPFCand_nameRECO_

private

Definition at line 164 of file DQMPFCandidateAnalyzer.h.

◆ occupancyPFCand_puppiNolepWeight_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::occupancyPFCand_puppiNolepWeight_

private

Definition at line 157 of file DQMPFCandidateAnalyzer.h.

◆ occupancyPFCandRECO_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::occupancyPFCandRECO_

private

Definition at line 163 of file DQMPFCandidateAnalyzer.h.

◆ pflowPackedToken_

edm::EDGetTokenT<std::vector<pat::PackedCandidate> > DQMPFCandidateAnalyzer::pflowPackedToken_

private

Definition at line 104 of file DQMPFCandidateAnalyzer.h.

◆ pflowToken_

edm::EDGetTokenT<std::vector<reco::PFCandidate> > DQMPFCandidateAnalyzer::pflowToken_

private

Definition at line 103 of file DQMPFCandidateAnalyzer.h.

◆ ptMinCand_

double DQMPFCandidateAnalyzer::ptMinCand_

private

Definition at line 116 of file DQMPFCandidateAnalyzer.h.

◆ ptPFCand_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::ptPFCand_

private

Definition at line 155 of file DQMPFCandidateAnalyzer.h.

◆ ptPFCand_name_

std::vector<std::string> DQMPFCandidateAnalyzer::ptPFCand_name_

private

Definition at line 156 of file DQMPFCandidateAnalyzer.h.

◆ ptPFCand_name_puppiNolepWeight_

std::vector<std::string> DQMPFCandidateAnalyzer::ptPFCand_name_puppiNolepWeight_

private

Definition at line 158 of file DQMPFCandidateAnalyzer.h.

◆ ptPFCand_nameRECO_

std::vector<std::string> DQMPFCandidateAnalyzer::ptPFCand_nameRECO_

private

Definition at line 164 of file DQMPFCandidateAnalyzer.h.

◆ ptPFCand_puppiNolepWeight_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::ptPFCand_puppiNolepWeight_

private

Definition at line 157 of file DQMPFCandidateAnalyzer.h.

◆ ptPFCandRECO_

std::vector<MonitorElement*> DQMPFCandidateAnalyzer::ptPFCandRECO_

private

Definition at line 163 of file DQMPFCandidateAnalyzer.h.

◆ theTriggerResultsLabel_

edm::InputTag DQMPFCandidateAnalyzer::theTriggerResultsLabel_

private

Definition at line 87 of file DQMPFCandidateAnalyzer.h.

◆ typePFCand_

std::vector<int> DQMPFCandidateAnalyzer::typePFCand_

private

Definition at line 160 of file DQMPFCandidateAnalyzer.h.

◆ typePFCandRECO_

std::vector<int> DQMPFCandidateAnalyzer::typePFCandRECO_

private

Definition at line 166 of file DQMPFCandidateAnalyzer.h.

◆ verbose_

int DQMPFCandidateAnalyzer::verbose_

private

Definition at line 169 of file DQMPFCandidateAnalyzer.h.

◆ vertexTag_

edm::InputTag DQMPFCandidateAnalyzer::vertexTag_

private

Definition at line 83 of file DQMPFCandidateAnalyzer.h.

◆ vertexToken_

edm::EDGetTokenT<std::vector<reco::Vertex> > DQMPFCandidateAnalyzer::vertexToken_

private

Definition at line 84 of file DQMPFCandidateAnalyzer.h.

Public Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DQMPFCandidateAnalyzer()

◆ ~DQMPFCandidateAnalyzer()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

◆ jetSortingRule()

Member Data Documentation

◆ bypassAllDCSChecks_

◆ bypassAllPVChecks_

◆ candidateType_

◆ cleaningParameters_

◆ countsPFCand_

◆ countsPFCandRECO_

◆ DCSFilter_

◆ diagnosticsParameters_

◆ etaMaxPFCand_

◆ etaMaxPFCandRECO_

◆ etaMinPFCand_

◆ etaMinPFCandRECO_

◆ FilterhltConfig_

◆ hbhenoifilterdecision

◆ hbheNoiseFilterResultTag_

◆ hbheNoiseFilterResultToken_

◆ HBHENoiseStringMiniAOD

◆ hcalMin_

◆ isMiniAO_

◆ isMiniAOD_

◆ LSBegin_

◆ LSEnd_

◆ m_HOverTrackP_Barrel_hPt_10_20

◆ m_HOverTrackP_Barrel_hPt_1_10

◆ m_HOverTrackP_Barrel_hPt_20_50

◆ m_HOverTrackP_Barrel_hPt_50

◆ m_HOverTrackP_EndCap_hPt_10_20

◆ m_HOverTrackP_EndCap_hPt_1_10

◆ m_HOverTrackP_EndCap_hPt_20_50

◆ m_HOverTrackP_EndCap_hPt_50

◆ m_HOverTrackP_trackPtVsEta

◆ m_HOverTrackPVsEta_hPt_10_20

◆ m_HOverTrackPVsEta_hPt_1_10

◆ m_HOverTrackPVsEta_hPt_20_50

◆ m_HOverTrackPVsEta_hPt_50

◆ m_HOverTrackPVsTrackP_Barrel

◆ m_HOverTrackPVsTrackP_EndCap

◆ m_HOverTrackPVsTrackPt_Barrel

◆ m_HOverTrackPVsTrackPt_EndCap

◆ map_of_MEs

◆ METFilterMiniAODLabel2_

◆ METFilterMiniAODLabel_

◆ METFilterMiniAODToken2_

◆ METFilterMiniAODToken_

◆ miniaodfilterdec

◆ miniaodfilterindex

◆ mInputCollection_

◆ mProfileIsoPFChHad_EcalOccupancyCentral

◆ mProfileIsoPFChHad_EcalOccupancyEndcap

◆ mProfileIsoPFChHad_EMPtCentral

◆ mProfileIsoPFChHad_EMPtEndcap

◆ mProfileIsoPFChHad_HadPtCentral

◆ mProfileIsoPFChHad_HadPtEndcap

◆ mProfileIsoPFChHad_HcalOccupancyCentral

◆ mProfileIsoPFChHad_HcalOccupancyEndcap

◆ mProfileIsoPFChHad_TrackOccupancy

◆ mProfileIsoPFChHad_TrackPt

◆ multiplicityPFCand_

◆ multiplicityPFCand_name_

◆ multiplicityPFCand_nameRECO_

◆ multiplicityPFCandRECO_

◆ numPV_

◆ occupancyPFCand_

◆ occupancyPFCand_name_

◆ occupancyPFCand_name_puppiNolepWeight_

◆ occupancyPFCand_nameRECO_