CMS 3D CMS Logo
List of all members | Public Member Functions | Static Private Member Functions | Private Attributes
DQMPFCandidateAnalyzer Class Reference

#include <DQMPFCandidateAnalyzer.h>

Inheritance diagram for DQMPFCandidateAnalyzer:
DQMOneEDAnalyzer<> edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

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 final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~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 &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > 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 & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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
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<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 declare what type of product will make and with which optional label More...
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
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)
 
template<Transition B>
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 (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes () noexcept
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag) noexcept
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
- 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.

63 {
64  miniaodfilterdec = -1;
65 
66  candidateType_ = pSet.getUntrackedParameter<std::string>("CandType");
67  //here only choice between miniaod or reco
68 
69  LSBegin_ = pSet.getParameter<int>("LSBegin");
70  LSEnd_ = pSet.getParameter<int>("LSEnd");
71 
72  isMiniAOD_ = (std::string("Packed") == candidateType_);
73 
74  mInputCollection_ = pSet.getParameter<edm::InputTag>("PFCandidateLabel");
75 
76  if (isMiniAOD_) {
77  pflowPackedToken_ = consumes<std::vector<pat::PackedCandidate> >(mInputCollection_);
78  } else {
79  pflowToken_ = consumes<std::vector<reco::PFCandidate> >(mInputCollection_);
80  }
81 
82  miniaodfilterdec = -1;
83 
84  // Smallest track pt
85  ptMinCand_ = pSet.getParameter<double>("ptMinCand");
86 
87  // Smallest raw HCAL energy linked to the track
88  hcalMin_ = pSet.getParameter<double>("hcalMin");
89 
90  diagnosticsParameters_ = pSet.getParameter<std::vector<edm::ParameterSet> >("METDiagonisticsParameters");
91 
92  edm::ConsumesCollector iC = consumesCollector();
93  //DCS
94  DCSFilter_ = new JetMETDQMDCSFilter(pSet.getParameter<ParameterSet>("DCSFilter"), iC);
95  if (isMiniAOD_) {
96  METFilterMiniAODLabel_ = pSet.getParameter<edm::InputTag>("FilterResultsLabelMiniAOD");
97  METFilterMiniAODToken_ = consumes<edm::TriggerResults>(METFilterMiniAODLabel_);
98 
99  METFilterMiniAODLabel2_ = pSet.getParameter<edm::InputTag>("FilterResultsLabelMiniAOD2");
100  METFilterMiniAODToken2_ = consumes<edm::TriggerResults>(METFilterMiniAODLabel2_);
101 
102  HBHENoiseStringMiniAOD = pSet.getParameter<std::string>("HBHENoiseLabelMiniAOD");
103  }
104 
105  if (!isMiniAOD_) {
106  hbheNoiseFilterResultTag_ = pSet.getParameter<edm::InputTag>("HBHENoiseFilterResultLabel");
107  hbheNoiseFilterResultToken_ = consumes<bool>(hbheNoiseFilterResultTag_);
108  }
109  //jet cleanup parameters
110  cleaningParameters_ = pSet.getParameter<ParameterSet>("CleaningParameters");
111 
112  //Vertex requirements
113  bypassAllPVChecks_ = cleaningParameters_.getParameter<bool>("bypassAllPVChecks");
114  bypassAllDCSChecks_ = cleaningParameters_.getParameter<bool>("bypassAllDCSChecks");
115  vertexTag_ = cleaningParameters_.getParameter<edm::InputTag>("vertexCollection");
116  vertexToken_ = consumes<std::vector<reco::Vertex> >(edm::InputTag(vertexTag_));
117 
118  verbose_ = pSet.getParameter<int>("verbose");
119 }

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

◆ ~DQMPFCandidateAnalyzer()

DQMPFCandidateAnalyzer::~DQMPFCandidateAnalyzer ( )
override

Destructor.

Definition at line 122 of file DQMPFCandidateAnalyzer.cc.

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

References LogTrace.

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.

666  {
667  //Vertex information
668  Handle<VertexCollection> vertexHandle;
669  iEvent.getByToken(vertexToken_, vertexHandle);
670 
671  if (!vertexHandle.isValid()) {
672  LogDebug("") << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
673  if (verbose_)
674  std::cout << "CaloMETAnalyzer: Could not find vertex collection" << std::endl;
675  }
676  numPV_ = 0;
677  if (vertexHandle.isValid()) {
678  VertexCollection vertexCollection = *(vertexHandle.product());
679  numPV_ = vertexCollection.size();
680  }
681  bool bPrimaryVertex = (bypassAllPVChecks_ || (numPV_ > 0));
682 
683  int myLuminosityBlock;
684  myLuminosityBlock = iEvent.luminosityBlock();
685 
686  if (myLuminosityBlock < LSBegin_)
687  return;
688  if (myLuminosityBlock > LSEnd_ && LSEnd_ > 0)
689  return;
690 
691  if (verbose_)
692  std::cout << "METAnalyzer analyze" << std::endl;
693 
694  std::string DirName = "JetMET/PFCandidates/" + mInputCollection_.label();
695 
696  bool hbhenoifilterdecision = true;
697  if (!isMiniAOD_) { //not checked for MiniAOD -> for miniaod decision filled as "triggerResults" bool
698  edm::Handle<bool> HBHENoiseFilterResultHandle;
699  iEvent.getByToken(hbheNoiseFilterResultToken_, HBHENoiseFilterResultHandle);
700  if (!HBHENoiseFilterResultHandle.isValid()) {
701  LogDebug("") << "METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
702  if (verbose_)
703  std::cout << "METAnalyzer: Could not find HBHENoiseFilterResult" << std::endl;
704  }
705  hbhenoifilterdecision = *HBHENoiseFilterResultHandle;
706  } else { //need to check if we go for version 1 or version 2
707  edm::Handle<edm::TriggerResults> metFilterResults;
708  iEvent.getByToken(METFilterMiniAODToken_, metFilterResults);
709  if (metFilterResults.isValid()) {
710  if (miniaodfilterindex != -1) {
711  hbhenoifilterdecision = metFilterResults->accept(miniaodfilterindex);
712  }
713  } else {
714  iEvent.getByToken(METFilterMiniAODToken2_, metFilterResults);
715  if (metFilterResults.isValid()) {
716  if (miniaodfilterindex != -1) {
717  hbhenoifilterdecision = metFilterResults->accept(miniaodfilterindex);
718  }
719  }
720  }
721  }
722 
723  //DCS Filter
724  bool bDCSFilter = (bypassAllDCSChecks_ || DCSFilter_->filter(iEvent, iSetup));
725 
726  for (unsigned int i = 0; i < countsPFCand_.size(); i++) {
727  countsPFCand_[i] = 0;
728  }
729  if (bDCSFilter && hbhenoifilterdecision && bPrimaryVertex) {
730  if (isMiniAOD_) {
731  edm::Handle<std::vector<pat::PackedCandidate> > packedParticleFlow;
732  iEvent.getByToken(pflowPackedToken_, packedParticleFlow);
733  //11, 13, 22 for el/mu/gamma, 211 chargedHadron, 130 neutralHadrons, 1 and 2 hadHF and EGammaHF
734  for (unsigned int i = 0; i < packedParticleFlow->size(); i++) {
735  const pat::PackedCandidate& c = packedParticleFlow->at(i);
736  for (unsigned int j = 0; j < typePFCand_.size(); j++) {
737  if (abs(c.pdgId()) == typePFCand_[j]) {
738  //second check for endcap, if inside barrel Max and Min symmetric around 0
739  if (((c.eta() > etaMinPFCand_[j]) && (c.eta() < etaMaxPFCand_[j])) ||
740  ((c.eta() > (-etaMaxPFCand_[j])) && (c.eta() < (-etaMinPFCand_[j])))) {
741  countsPFCand_[j] += 1;
742  ptPFCand_[j] = map_of_MEs[DirName + "/" + ptPFCand_name_[j]];
743  if (ptPFCand_[j] && ptPFCand_[j]->getRootObject())
744  ptPFCand_[j]->Fill(c.eta(), c.phi(), c.pt() * c.puppiWeight());
745  occupancyPFCand_[j] = map_of_MEs[DirName + "/" + occupancyPFCand_name_[j]];
746  if (occupancyPFCand_[j] && occupancyPFCand_[j]->getRootObject())
747  occupancyPFCand_[j]->Fill(c.eta(), c.phi(), c.puppiWeight());
748  ptPFCand_puppiNolepWeight_[j] = map_of_MEs[DirName + "/" + ptPFCand_name_puppiNolepWeight_[j]];
749  if (ptPFCand_puppiNolepWeight_[j] && ptPFCand_puppiNolepWeight_[j]->getRootObject())
750  ptPFCand_puppiNolepWeight_[j]->Fill(c.eta(), c.phi(), c.pt() * c.puppiWeightNoLep());
751  occupancyPFCand_puppiNolepWeight_[j] =
752  map_of_MEs[DirName + "/" + occupancyPFCand_name_puppiNolepWeight_[j]];
753  if (occupancyPFCand_puppiNolepWeight_[j] && occupancyPFCand_puppiNolepWeight_[j]->getRootObject()) {
754  occupancyPFCand_puppiNolepWeight_[j]->Fill(c.eta(), c.phi(), c.puppiWeightNoLep());
755  }
756  }
757  }
758  }
759  }
760  for (unsigned int j = 0; j < countsPFCand_.size(); j++) {
761  multiplicityPFCand_[j] = map_of_MEs[DirName + "/" + multiplicityPFCand_name_[j]];
762  if (multiplicityPFCand_[j] && multiplicityPFCand_[j]->getRootObject()) {
763  multiplicityPFCand_[j]->Fill(countsPFCand_[j]);
764  }
765  }
766  } else {
767  edm::Handle<std::vector<reco::PFCandidate> > particleFlow;
768  iEvent.getByToken(pflowToken_, particleFlow);
769  for (unsigned int i = 0; i < particleFlow->size(); i++) {
770  const reco::PFCandidate& c = particleFlow->at(i);
771  for (unsigned int j = 0; j < typePFCandRECO_.size(); j++) {
772  if (c.particleId() == typePFCandRECO_[j]) {
773  //second check for endcap, if inside barrel Max and Min symmetric around 0
774  if (((c.eta() > etaMinPFCandRECO_[j]) && (c.eta() < etaMaxPFCandRECO_[j])) ||
775  ((c.eta() > (-etaMaxPFCandRECO_[j])) && (c.eta() < (-etaMinPFCandRECO_[j])))) {
776  countsPFCandRECO_[j] += 1;
777  ptPFCandRECO_[j] = map_of_MEs[DirName + "/" + ptPFCand_nameRECO_[j]];
778  if (ptPFCandRECO_[j] && ptPFCandRECO_[j]->getRootObject())
779  ptPFCandRECO_[j]->Fill(c.eta(), c.phi(), c.pt());
780  occupancyPFCandRECO_[j] = map_of_MEs[DirName + "/" + occupancyPFCand_nameRECO_[j]];
781  if (occupancyPFCandRECO_[j] && occupancyPFCandRECO_[j]->getRootObject())
782  occupancyPFCandRECO_[j]->Fill(c.eta(), c.phi());
783  }
784  //fill quantities for isolated charged hadron quantities
785  //only for charged hadrons
786  if (c.particleId() == 1 && c.pt() > ptMinCand_) {
787  // At least 1 GeV in HCAL
788  double ecalRaw = c.rawEcalEnergy();
789  double hcalRaw = c.rawHcalEnergy();
790  if ((ecalRaw + hcalRaw) > hcalMin_) {
791  const PFCandidate::ElementsInBlocks& theElements = c.elementsInBlocks();
792  if (theElements.empty())
793  continue;
794  unsigned int iTrack = -999;
795  std::vector<unsigned int> iECAL; // =999;
796  std::vector<unsigned int> iHCAL; // =999;
797  const reco::PFBlockRef blockRef = theElements[0].first;
798  const edm::OwnVector<reco::PFBlockElement>& elements = blockRef->elements();
799  // Check that there is only one track in the block.
800  unsigned int nTracks = 0;
801  for (unsigned int iEle = 0; iEle < elements.size(); iEle++) {
802  // Find the tracks in the block
803  PFBlockElement::Type type = elements[iEle].type();
804  switch (type) {
805  case PFBlockElement::TRACK:
806  iTrack = iEle;
807  nTracks++;
808  break;
809  case PFBlockElement::ECAL:
810  iECAL.push_back(iEle);
811  break;
812  case PFBlockElement::HCAL:
813  iHCAL.push_back(iEle);
814  break;
815  default:
816  continue;
817  }
818  }
819  if (nTracks == 1) {
820  // Characteristics of the track
821  const reco::PFBlockElementTrack& et =
822  dynamic_cast<const reco::PFBlockElementTrack&>(elements[iTrack]);
823  mProfileIsoPFChHad_TrackOccupancy = map_of_MEs[DirName + "/" + "IsoPfChHad_Track_profile"];
824  if (mProfileIsoPFChHad_TrackOccupancy && mProfileIsoPFChHad_TrackOccupancy->getRootObject())
825  mProfileIsoPFChHad_TrackOccupancy->Fill(et.trackRef()->eta(), et.trackRef()->phi());
826  mProfileIsoPFChHad_TrackPt = map_of_MEs[DirName + "/" + "IsoPfChHad_TrackPt"];
827  if (mProfileIsoPFChHad_TrackPt && mProfileIsoPFChHad_TrackPt->getRootObject())
828  mProfileIsoPFChHad_TrackPt->Fill(et.trackRef()->eta(), et.trackRef()->phi(), et.trackRef()->pt());
829  if (c.rawEcalEnergy() == 0) { //isolated hadron, nothing in ECAL
830  //right now take corrected hcalEnergy, do we want the rawHcalEnergy instead
831  m_HOverTrackP_trackPtVsEta = map_of_MEs[DirName + "/" + "HOverTrackP_trackPtVsEta"];
832  if (m_HOverTrackP_trackPtVsEta && m_HOverTrackP_trackPtVsEta->getRootObject())
833  m_HOverTrackP_trackPtVsEta->Fill(c.pt(), c.eta(), c.hcalEnergy() / et.trackRef()->p());
834  if (c.pt() > 1 && c.pt() < 10) {
835  m_HOverTrackPVsEta_hPt_1_10 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_1_10"];
836  if (m_HOverTrackPVsEta_hPt_1_10 && m_HOverTrackPVsEta_hPt_1_10->getRootObject())
837  m_HOverTrackPVsEta_hPt_1_10->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
838  } else if (c.pt() > 10 && c.pt() < 20) {
839  m_HOverTrackPVsEta_hPt_10_20 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_10_20"];
840  if (m_HOverTrackPVsEta_hPt_10_20 && m_HOverTrackPVsEta_hPt_10_20->getRootObject())
841  m_HOverTrackPVsEta_hPt_10_20->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
842  } else if (c.pt() > 20 && c.pt() < 50) {
843  m_HOverTrackPVsEta_hPt_20_50 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_20_50"];
844  if (m_HOverTrackPVsEta_hPt_20_50 && m_HOverTrackPVsEta_hPt_20_50->getRootObject())
845  m_HOverTrackPVsEta_hPt_20_50->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
846  } else if (c.pt() > 50) {
847  m_HOverTrackPVsEta_hPt_50 = map_of_MEs[DirName + "/" + "HOverTrackPVsEta_hPt_50"];
848  if (m_HOverTrackPVsEta_hPt_50 && m_HOverTrackPVsEta_hPt_50->getRootObject())
849  m_HOverTrackPVsEta_hPt_50->Fill(c.eta(), c.hcalEnergy() / et.trackRef()->p());
850  }
851  if (fabs(c.eta() < 1.392)) {
852  if (c.pt() > 1 && c.pt() < 10) {
853  m_HOverTrackP_Barrel_hPt_1_10 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_1_10"];
854  if (m_HOverTrackP_Barrel_hPt_1_10 && m_HOverTrackP_Barrel_hPt_1_10->getRootObject())
855  m_HOverTrackP_Barrel_hPt_1_10->Fill(c.hcalEnergy() / et.trackRef()->p());
856  } else if (c.pt() > 10 && c.pt() < 20) {
857  m_HOverTrackP_Barrel_hPt_10_20 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_10_20"];
858  if (m_HOverTrackP_Barrel_hPt_10_20 && m_HOverTrackP_Barrel_hPt_10_20->getRootObject())
859  m_HOverTrackP_Barrel_hPt_10_20->Fill(c.hcalEnergy() / et.trackRef()->p());
860  } else if (c.pt() > 20 && c.pt() < 50) {
861  m_HOverTrackP_Barrel_hPt_20_50 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_20_50"];
862  if (m_HOverTrackP_Barrel_hPt_20_50 && m_HOverTrackP_Barrel_hPt_20_50->getRootObject())
863  m_HOverTrackP_Barrel_hPt_20_50->Fill(c.hcalEnergy() / et.trackRef()->p());
864  } else if (c.pt() > 50) {
865  m_HOverTrackP_Barrel_hPt_50 = map_of_MEs[DirName + "/" + "HOverTrackP_Barrel_hPt_50"];
866  if (m_HOverTrackP_Barrel_hPt_50 && m_HOverTrackP_Barrel_hPt_50->getRootObject())
867  m_HOverTrackP_Barrel_hPt_50->Fill(c.hcalEnergy() / et.trackRef()->p());
868  }
869  m_HOverTrackPVsTrackP_Barrel = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackP_Barrel"];
870  if (m_HOverTrackPVsTrackP_Barrel && m_HOverTrackPVsTrackP_Barrel->getRootObject())
871  m_HOverTrackPVsTrackP_Barrel->Fill(et.trackRef()->p(), c.hcalEnergy() / et.trackRef()->p());
872  m_HOverTrackPVsTrackPt_Barrel = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackPt_Barrel"];
873  if (m_HOverTrackPVsTrackPt_Barrel && m_HOverTrackPVsTrackPt_Barrel->getRootObject())
874  m_HOverTrackPVsTrackPt_Barrel->Fill(et.trackRef()->pt(), c.hcalEnergy() / et.trackRef()->p());
875  } else {
876  m_HOverTrackPVsTrackP_EndCap = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackP_EndCap"];
877  if (m_HOverTrackPVsTrackP_EndCap && m_HOverTrackPVsTrackP_EndCap->getRootObject())
878  m_HOverTrackPVsTrackP_EndCap->Fill(et.trackRef()->p(), c.hcalEnergy() / et.trackRef()->p());
879  m_HOverTrackPVsTrackPt_EndCap = map_of_MEs[DirName + "/" + "HOverTrackPVsTrackPt_EndCap"];
880  if (m_HOverTrackPVsTrackPt_EndCap && m_HOverTrackPVsTrackPt_EndCap->getRootObject())
881  m_HOverTrackPVsTrackPt_EndCap->Fill(et.trackRef()->pt(), c.hcalEnergy() / et.trackRef()->p());
882  if (c.pt() > 1 && c.pt() < 10) {
883  m_HOverTrackP_EndCap_hPt_1_10 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_1_10"];
884  if (m_HOverTrackP_EndCap_hPt_1_10 && m_HOverTrackP_EndCap_hPt_1_10->getRootObject())
885  m_HOverTrackP_EndCap_hPt_1_10->Fill(c.hcalEnergy() / et.trackRef()->p());
886  } else if (c.pt() > 10 && c.pt() < 20) {
887  m_HOverTrackP_EndCap_hPt_10_20 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_10_20"];
888  if (m_HOverTrackP_EndCap_hPt_10_20 && m_HOverTrackP_EndCap_hPt_10_20->getRootObject())
889  m_HOverTrackP_EndCap_hPt_10_20->Fill(c.hcalEnergy() / et.trackRef()->p());
890  } else if (c.pt() > 20 && c.pt() < 50) {
891  m_HOverTrackP_EndCap_hPt_20_50 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_20_50"];
892  if (m_HOverTrackP_EndCap_hPt_20_50 && m_HOverTrackP_EndCap_hPt_20_50->getRootObject())
893  m_HOverTrackP_EndCap_hPt_20_50->Fill(c.hcalEnergy() / et.trackRef()->p());
894  } else if (c.pt() > 50) {
895  m_HOverTrackP_EndCap_hPt_50 = map_of_MEs[DirName + "/" + "HOverTrackP_EndCap_hPt_50"];
896  if (m_HOverTrackP_EndCap_hPt_50 && m_HOverTrackP_EndCap_hPt_50->getRootObject())
897  m_HOverTrackP_EndCap_hPt_50->Fill(c.hcalEnergy() / et.trackRef()->p());
898  }
899  }
900  }
901  //ECAL element
902  for (unsigned int ii = 0; ii < iECAL.size(); ii++) {
903  const reco::PFBlockElementCluster& eecal =
904  dynamic_cast<const reco::PFBlockElementCluster&>(elements[iECAL[ii]]);
905  if (fabs(eecal.clusterRef()->eta()) < 1.479) {
906  mProfileIsoPFChHad_EcalOccupancyCentral =
907  map_of_MEs[DirName + "/" + "IsoPfChHad_ECAL_profile_central"];
908  if (mProfileIsoPFChHad_EcalOccupancyCentral &&
909  mProfileIsoPFChHad_EcalOccupancyCentral->getRootObject())
910  mProfileIsoPFChHad_EcalOccupancyCentral->Fill(eecal.clusterRef()->eta(),
911  eecal.clusterRef()->phi());
912  mProfileIsoPFChHad_EMPtCentral = map_of_MEs[DirName + "/" + "IsoPfChHad_EMPt_central"];
913  if (mProfileIsoPFChHad_EMPtCentral && mProfileIsoPFChHad_EMPtCentral->getRootObject())
914  mProfileIsoPFChHad_EMPtCentral->Fill(
915  eecal.clusterRef()->eta(), eecal.clusterRef()->phi(), eecal.clusterRef()->pt());
916  } else {
917  mProfileIsoPFChHad_EcalOccupancyEndcap =
918  map_of_MEs[DirName + "/" + "IsoPfChHad_ECAL_profile_endcap"];
919  if (mProfileIsoPFChHad_EcalOccupancyEndcap &&
920  mProfileIsoPFChHad_EcalOccupancyEndcap->getRootObject())
921  mProfileIsoPFChHad_EcalOccupancyEndcap->Fill(eecal.clusterRef()->eta(),
922  eecal.clusterRef()->phi());
923  mProfileIsoPFChHad_EMPtEndcap = map_of_MEs[DirName + "/" + "IsoPfChHad_EMPt_endcap"];
924  if (mProfileIsoPFChHad_EMPtEndcap && mProfileIsoPFChHad_EMPtEndcap->getRootObject())
925  mProfileIsoPFChHad_EMPtEndcap->Fill(
926  eecal.clusterRef()->eta(), eecal.clusterRef()->phi(), eecal.clusterRef()->pt());
927  }
928  }
929  //HCAL element
930  for (unsigned int ii = 0; ii < iHCAL.size(); ii++) {
931  const reco::PFBlockElementCluster& ehcal =
932  dynamic_cast<const reco::PFBlockElementCluster&>(elements[iHCAL[ii]]);
933  if (fabs(ehcal.clusterRef()->eta()) < 1.740) {
934  mProfileIsoPFChHad_HcalOccupancyCentral =
935  map_of_MEs[DirName + "/" + "IsoPfChHad_HCAL_profile_central"];
936  if (mProfileIsoPFChHad_HcalOccupancyCentral &&
937  mProfileIsoPFChHad_HcalOccupancyCentral->getRootObject())
938  mProfileIsoPFChHad_HcalOccupancyCentral->Fill(ehcal.clusterRef()->eta(),
939  ehcal.clusterRef()->phi());
940  mProfileIsoPFChHad_HadPtCentral = map_of_MEs[DirName + "/" + "IsoPfChHad_HadPt_central"];
941  if (mProfileIsoPFChHad_HadPtCentral && mProfileIsoPFChHad_HadPtCentral->getRootObject())
942  mProfileIsoPFChHad_HadPtCentral->Fill(
943  ehcal.clusterRef()->eta(), ehcal.clusterRef()->phi(), ehcal.clusterRef()->pt());
944  } else {
945  mProfileIsoPFChHad_HcalOccupancyEndcap =
946  map_of_MEs[DirName + "/" + "IsoPfChHad_HCAL_profile_endcap"];
947  if (mProfileIsoPFChHad_HcalOccupancyEndcap &&
948  mProfileIsoPFChHad_HcalOccupancyEndcap->getRootObject())
949  mProfileIsoPFChHad_HcalOccupancyEndcap->Fill(ehcal.clusterRef()->eta(),
950  ehcal.clusterRef()->phi());
951  mProfileIsoPFChHad_HadPtEndcap = map_of_MEs[DirName + "/" + "IsoPfChHad_HadPt_endcap"];
952  if (mProfileIsoPFChHad_HadPtEndcap && mProfileIsoPFChHad_HadPtEndcap->getRootObject())
953  mProfileIsoPFChHad_HadPtEndcap->Fill(
954  ehcal.clusterRef()->eta(), ehcal.clusterRef()->phi(), ehcal.clusterRef()->pt());
955  }
956  }
957  }
958  }
959  }
960  }
961  }
962  }
963  for (unsigned int j = 0; j < countsPFCandRECO_.size(); j++) {
964  multiplicityPFCandRECO_[j] = map_of_MEs[DirName + "/" + multiplicityPFCand_nameRECO_[j]];
965  if (multiplicityPFCandRECO_[j] && multiplicityPFCandRECO_[j]->getRootObject()) {
966  multiplicityPFCandRECO_[j]->Fill(countsPFCandRECO_[j]);
967  }
968  }
969  } //candidate loop for both miniaod and reco
970  }
971 }

References funct::abs(), edm::HLTGlobalStatus::accept(), c, reco::PFBlockElementCluster::clusterRef(), gather_cfg::cout, fed_dqm_sourceclient-live_cfg::DirName, ECAL, bookConverter::elements, EgHLTOffHistBins_cfi::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.

◆ 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.

128  {
129  ibooker.setCurrentFolder("JetMET/PFCandidates/" + mInputCollection_.label());
130  std::string DirName = "JetMET/PFCandidates/" + mInputCollection_.label();
131 
132  if (!isMiniAOD_) {
133  if (!occupancyPFCandRECO_.empty())
134  occupancyPFCandRECO_.clear();
135  if (!occupancyPFCand_nameRECO_.empty())
136  occupancyPFCand_nameRECO_.clear();
137  if (!etaMinPFCandRECO_.empty())
138  etaMinPFCandRECO_.clear();
139  if (!etaMaxPFCandRECO_.empty())
140  etaMaxPFCandRECO_.clear();
141  if (!typePFCandRECO_.empty())
142  typePFCandRECO_.clear();
143  if (!countsPFCandRECO_.empty())
144  countsPFCandRECO_.clear();
145  if (!ptPFCandRECO_.empty())
146  ptPFCandRECO_.clear();
147  if (!ptPFCand_nameRECO_.empty())
148  ptPFCand_nameRECO_.clear();
149  if (!multiplicityPFCandRECO_.empty())
150  multiplicityPFCandRECO_.clear();
151  if (!multiplicityPFCand_nameRECO_.empty())
152  multiplicityPFCand_nameRECO_.clear();
153  for (std::vector<edm::ParameterSet>::const_iterator v = diagnosticsParameters_.begin();
154  v != diagnosticsParameters_.end();
155  v++) {
156  int etaNBinsPFCand = v->getParameter<int>("etaNBins");
157  double etaMinPFCand = v->getParameter<double>("etaMin");
158  double etaMaxPFCand = v->getParameter<double>("etaMax");
159  int phiNBinsPFCand = v->getParameter<int>("phiNBins");
160  double phiMinPFCand = v->getParameter<double>("phiMin");
161  double phiMaxPFCand = v->getParameter<double>("phiMax");
162  int nMinPFCand = v->getParameter<int>("nMin");
163  int nMaxPFCand = v->getParameter<int>("nMax");
164  int nbinsPFCand = v->getParameter<double>("nbins");
165  etaMinPFCandRECO_.push_back(etaMinPFCand);
166  etaMaxPFCandRECO_.push_back(etaMaxPFCand);
167  typePFCandRECO_.push_back(v->getParameter<int>("type"));
168  countsPFCandRECO_.push_back(0);
169  multiplicityPFCandRECO_.push_back(
170  ibooker.book1D(std::string(v->getParameter<std::string>("name")).append("_multiplicity_").c_str(),
171  std::string(v->getParameter<std::string>("name")) + "multiplicity",
172  nbinsPFCand,
173  nMinPFCand,
174  nMaxPFCand));
175  multiplicityPFCand_nameRECO_.push_back(
176  std::string(v->getParameter<std::string>("name")).append("_multiplicity_"));
177  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
178  DirName + "/" + multiplicityPFCand_nameRECO_[multiplicityPFCand_nameRECO_.size() - 1],
179  multiplicityPFCandRECO_[multiplicityPFCandRECO_.size() - 1]));
180 
181  //push back names first, we need to create histograms with the name and fill it for endcap plots later
182  occupancyPFCand_nameRECO_.push_back(std::string(v->getParameter<std::string>("name")).append("_occupancy_"));
183 
184  ptPFCand_nameRECO_.push_back(std::string(v->getParameter<std::string>("name")).append("_pt_"));
185  //special booking for endcap plots, merge plots for eminus and eplus into one plot, using variable binning
186  //barrel plots have eta-boundaries on minus and plus side
187  //parameters start on minus side
188  if (etaMinPFCand * etaMaxPFCand < 0) { //barrel plots, plot only in barrel region
189  occupancyPFCandRECO_.push_back(
190  ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_occupancy_").c_str(),
191  std::string(v->getParameter<std::string>("name")) + "occupancy",
192  etaNBinsPFCand,
193  etaMinPFCand,
194  etaMaxPFCand,
195  phiNBinsPFCand,
196  phiMinPFCand,
197  phiMaxPFCand));
198  ptPFCandRECO_.push_back(ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_").c_str(),
199  std::string(v->getParameter<std::string>("name")) + "pt",
200  etaNBinsPFCand,
201  etaMinPFCand,
202  etaMaxPFCand,
203  phiNBinsPFCand,
204  phiMinPFCand,
205  phiMaxPFCand));
206  } else { //endcap or forward plots,
207  const int nbins_eta_endcap = 2 * (etaNBinsPFCand + 1);
208  double eta_limits_endcap[nbins_eta_endcap];
209  for (int i = 0; i < nbins_eta_endcap; i++) {
210  if (i < (etaNBinsPFCand + 1)) {
211  eta_limits_endcap[i] = etaMinPFCand + i * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
212  } else {
213  eta_limits_endcap[i] =
214  -etaMaxPFCand + (i - (etaNBinsPFCand + 1)) * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
215  }
216  }
217  TH2F* hist_temp_occup = new TH2F((occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size() - 1]).c_str(),
218  "occupancy",
219  nbins_eta_endcap - 1,
220  eta_limits_endcap,
221  phiNBinsPFCand,
222  phiMinPFCand,
223  phiMaxPFCand);
224  occupancyPFCandRECO_.push_back(
225  ibooker.book2D(occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size() - 1], hist_temp_occup));
226  TH2F* hist_temp_pt = new TH2F((ptPFCand_nameRECO_[ptPFCand_nameRECO_.size() - 1]).c_str(),
227  "pt",
228  nbins_eta_endcap - 1,
229  eta_limits_endcap,
230  phiNBinsPFCand,
231  phiMinPFCand,
232  phiMaxPFCand);
233  ptPFCandRECO_.push_back(ibooker.book2D(ptPFCand_nameRECO_[ptPFCand_nameRECO_.size() - 1], hist_temp_pt));
234  }
235 
236  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
237  DirName + "/" + occupancyPFCand_nameRECO_[occupancyPFCand_nameRECO_.size() - 1],
238  occupancyPFCandRECO_[occupancyPFCandRECO_.size() - 1]));
239  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
240  DirName + "/" + ptPFCand_nameRECO_[ptPFCand_nameRECO_.size() - 1], ptPFCandRECO_[ptPFCandRECO_.size() - 1]));
241  }
242 
243  mProfileIsoPFChHad_TrackOccupancy = ibooker.book2D(
244  "IsoPfChHad_Track_profile", "Isolated PFChHadron Tracker_occupancy", 108, -2.7, 2.7, 160, -M_PI, M_PI);
245  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_Track_profile",
246  mProfileIsoPFChHad_TrackOccupancy));
247  mProfileIsoPFChHad_TrackPt =
248  ibooker.book2D("IsoPfChHad_TrackPt", "Isolated PFChHadron TrackPt", 108, -2.7, 2.7, 160, -M_PI, M_PI);
249  map_of_MEs.insert(
250  std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_TrackPt", mProfileIsoPFChHad_TrackPt));
251 
252  mProfileIsoPFChHad_EcalOccupancyCentral = ibooker.book2D("IsoPfChHad_ECAL_profile_central",
253  "IsolatedPFChHa ECAL occupancy (Barrel)",
254  180,
255  -1.479,
256  1.479,
257  125,
258  -M_PI,
259  M_PI);
260  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_ECAL_profile_central",
261  mProfileIsoPFChHad_EcalOccupancyCentral));
262  mProfileIsoPFChHad_EMPtCentral = ibooker.book2D(
263  "IsoPfChHad_EMPt_central", "Isolated PFChHadron HadPt_central", 180, -1.479, 1.479, 360, -M_PI, M_PI);
264  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_EMPt_central",
265  mProfileIsoPFChHad_EMPtCentral));
266 
267  mProfileIsoPFChHad_EcalOccupancyEndcap = ibooker.book2D(
268  "IsoPfChHad_ECAL_profile_endcap", "IsolatedPFChHa ECAL occupancy (Endcap)", 110, -2.75, 2.75, 125, -M_PI, M_PI);
269  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_ECAL_profile_endcap",
270  mProfileIsoPFChHad_EcalOccupancyEndcap));
271  mProfileIsoPFChHad_EMPtEndcap =
272  ibooker.book2D("IsoPfChHad_EMPt_endcap", "Isolated PFChHadron EMPt_endcap", 110, -2.75, 2.75, 125, -M_PI, M_PI);
273  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_EMPt_endcap",
274  mProfileIsoPFChHad_EMPtEndcap));
275 
276  const int nbins_eta = 16;
277 
278  double eta_limits[nbins_eta] = {-2.650,
279  -2.500,
280  -2.322,
281  -2.172,
282  -2.043,
283  -1.930,
284  -1.830,
285  -1.740,
286  1.740,
287  1.830,
288  1.930,
289  2.043,
290  2.172,
291  2.3122,
292  2.500,
293  2.650};
294 
295  TH2F* hist_temp_HCAL = new TH2F("IsoPfChHad_HCAL_profile_endcap",
296  "IsolatedPFChHa HCAL occupancy (outer endcap)",
297  nbins_eta - 1,
298  eta_limits,
299  36,
300  -M_PI,
301  M_PI);
302  TH2F* hist_tempPt_HCAL = (TH2F*)hist_temp_HCAL->Clone("Isolated PFCHHadron HadPt (outer endcap)");
303 
304  mProfileIsoPFChHad_HcalOccupancyCentral = ibooker.book2D("IsoPfChHad_HCAL_profile_central",
305  "IsolatedPFChHa HCAL occupancy (Central Part)",
306  40,
307  -1.740,
308  1.740,
309  72,
310  -M_PI,
311  M_PI);
312  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HCAL_profile_central",
313  mProfileIsoPFChHad_HcalOccupancyCentral));
314  mProfileIsoPFChHad_HadPtCentral = ibooker.book2D(
315  "IsoPfChHad_HadPt_central", "Isolated PFChHadron HadPt_central", 40, -1.740, 1.740, 72, -M_PI, M_PI);
316  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HadPt_central",
317  mProfileIsoPFChHad_HadPtCentral));
318 
319  mProfileIsoPFChHad_HcalOccupancyEndcap = ibooker.book2D("IsoPfChHad_HCAL_profile_endcap", hist_temp_HCAL);
320  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HCAL_profile_endcap",
321  mProfileIsoPFChHad_HcalOccupancyEndcap));
322  mProfileIsoPFChHad_HadPtEndcap = ibooker.book2D("IsoPfChHad_HadPt_endcap", hist_tempPt_HCAL);
323  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "IsoPfChHad_HadPt_endcap",
324  mProfileIsoPFChHad_HadPtEndcap));
325 
326  //actual HCAL segmentation in pseudorapidity -> reduce by a factor of two
327  //const int nbins_eta_hcal_total=54;
328  //double eta_limits_hcal_total[nbins_eta_hcal_total]=
329  // {-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,
330  // -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,
331  // 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
332  // 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}
333  //
334 
335  const int nbins_eta_hcal_total = 28;
336  double eta_limits_hcal_total[nbins_eta_hcal_total] = {
337  -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,
338  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};
339  float eta_limits_hcal_total_f[nbins_eta_hcal_total];
340  float log_bin_spacing = log(200.) / 40.;
341  const int nbins_pt_total_hcal = 41;
342  double pt_limits_hcal[nbins_pt_total_hcal];
343  float pt_limits_hcal_f[nbins_pt_total_hcal];
344  for (int i = 0; i < nbins_pt_total_hcal; i++) {
345  pt_limits_hcal[i] = exp(i * log_bin_spacing);
346  pt_limits_hcal_f[i] = exp(i * log_bin_spacing);
347  }
348  for (int i = 0; i < nbins_eta_hcal_total; i++) {
349  eta_limits_hcal_total_f[i] = (float)eta_limits_hcal_total[i];
350  }
351  m_HOverTrackP_trackPtVsEta = ibooker.book2D("HOverTrackP_trackPtVsEta",
352  "HOverTrackP_trackPtVsEta",
353  nbins_pt_total_hcal - 1,
354  pt_limits_hcal_f,
355  nbins_eta_hcal_total - 1,
356  eta_limits_hcal_total_f);
357  m_HOverTrackPVsTrackP_Barrel = ibooker.bookProfile(
358  "HOverTrackPVsTrackP_Barrel", "HOverTrackPVsTrackP_Barrel", nbins_pt_total_hcal - 1, pt_limits_hcal, 0, 4, " ");
359  m_HOverTrackPVsTrackP_EndCap = ibooker.bookProfile(
360  "HOverTrackPVsTrackP_EndCap", "HOverTrackPVsTrackP_EndCap", nbins_pt_total_hcal - 1, pt_limits_hcal, 0, 4, " ");
361  m_HOverTrackPVsTrackPt_Barrel = ibooker.bookProfile("HOverTrackPVsTrackPt_Barrel",
362  "HOverTrackPVsTrackPt_Barrel",
363  nbins_pt_total_hcal - 1,
364  pt_limits_hcal,
365  0,
366  4,
367  " ");
368  m_HOverTrackPVsTrackPt_EndCap = ibooker.bookProfile("HOverTrackPVsTrackPt_EndCap",
369  "HOverTrackPVsTrackPt_EndCap",
370  nbins_pt_total_hcal - 1,
371  pt_limits_hcal,
372  0,
373  4,
374  " ");
375 
376  m_HOverTrackPVsEta_hPt_1_10 = ibooker.bookProfile("HOverTrackPVsEta_hPt_1_10",
377  "HOverTrackPVsEta, 1<hPt<10 GeV",
378  nbins_eta_hcal_total - 1,
379  eta_limits_hcal_total,
380  0,
381  4,
382  " ");
383  m_HOverTrackPVsEta_hPt_10_20 = ibooker.bookProfile("HOverTrackPVsEta_hPt_10_20",
384  "HOverTrackPVsEta, 10<hPt<20 GeV",
385  nbins_eta_hcal_total - 1,
386  eta_limits_hcal_total,
387  0,
388  4,
389  " ");
390  m_HOverTrackPVsEta_hPt_20_50 = ibooker.bookProfile("HOverTrackPVsEta_hPt_20_50",
391  "HOverTrackPVsEta, 20<hPt<50 GeV",
392  nbins_eta_hcal_total - 1,
393  eta_limits_hcal_total,
394  0,
395  4,
396  " ");
397  m_HOverTrackPVsEta_hPt_50 = ibooker.bookProfile("HOverTrackPVsEta_hPt_50",
398  "HOverTrackPVsEta, hPt>50 GeV",
399  nbins_eta_hcal_total - 1,
400  eta_limits_hcal_total,
401  0,
402  4,
403  " ");
404 
405  m_HOverTrackP_Barrel_hPt_1_10 =
406  ibooker.book1D("HOverTrackP_Barrel_hPt_1_10", "HOverTrackP_B, 1<hPt<10 GeV", 50, 0, 4);
407  m_HOverTrackP_Barrel_hPt_10_20 =
408  ibooker.book1D("HOverTrackP_Barrel_hPt_10_20", "HOverTrackP_B, 10<hPt<20 GeV", 50, 0, 4);
409  m_HOverTrackP_Barrel_hPt_20_50 =
410  ibooker.book1D("HOverTrackP_Barrel_hPt_20_50", "HOverTrackP_B, 20<hPt<50 GeV", 50, 0, 4);
411  m_HOverTrackP_Barrel_hPt_50 = ibooker.book1D("HOverTrackP_Barrel_hPt_50", "HOverTrackP_B, hPt>50 GeV", 50, 0, 4);
412 
413  m_HOverTrackP_EndCap_hPt_1_10 =
414  ibooker.book1D("HOverTrackP_EndCap_hPt_1_10", "HOverTrackP_E, 1<hPt<10 GeV", 50, 0, 4);
415  m_HOverTrackP_EndCap_hPt_10_20 =
416  ibooker.book1D("HOverTrackP_EndCap_hPt_10_20", "HOverTrackP_E, 10<hPt<20 GeV", 50, 0, 4);
417  m_HOverTrackP_EndCap_hPt_20_50 =
418  ibooker.book1D("HOverTrackP_EndCap_hPt_20_50", "HOverTrackP_E, 20<hPt<50 GeV", 50, 0, 4);
419  m_HOverTrackP_EndCap_hPt_50 = ibooker.book1D("HOverTrackP_EndCap_hPt_50", "HOverTrackP_E, hPt>50 GeV", 50, 0, 4);
420 
421  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_trackPtVsEta",
422  m_HOverTrackP_trackPtVsEta));
423  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackP_Barrel",
424  m_HOverTrackPVsTrackP_Barrel));
425  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackP_EndCap",
426  m_HOverTrackPVsTrackP_EndCap));
427  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackPt_Barrel",
428  m_HOverTrackPVsTrackPt_Barrel));
429  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsTrackPt_EndCap",
430  m_HOverTrackPVsTrackPt_EndCap));
431  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_1_10",
432  m_HOverTrackPVsEta_hPt_1_10));
433  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_10_20",
434  m_HOverTrackPVsEta_hPt_10_20));
435  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_20_50",
436  m_HOverTrackPVsEta_hPt_20_50));
437  map_of_MEs.insert(
438  std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackPVsEta_hPt_50", m_HOverTrackPVsEta_hPt_50));
439  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_1_10",
440  m_HOverTrackP_Barrel_hPt_1_10));
441  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_10_20",
442  m_HOverTrackP_Barrel_hPt_10_20));
443  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_20_50",
444  m_HOverTrackP_Barrel_hPt_20_50));
445  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_Barrel_hPt_50",
446  m_HOverTrackP_Barrel_hPt_50));
447  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_1_10",
448  m_HOverTrackP_EndCap_hPt_1_10));
449  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_10_20",
450  m_HOverTrackP_EndCap_hPt_10_20));
451  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_20_50",
452  m_HOverTrackP_EndCap_hPt_20_50));
453  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(DirName + "/" + "HOverTrackP_EndCap_hPt_50",
454  m_HOverTrackP_EndCap_hPt_50));
455  } else { //MiniAOD workflow
456  if (!occupancyPFCand_.empty())
457  occupancyPFCand_.clear();
458  if (!occupancyPFCand_name_.empty())
459  occupancyPFCand_name_.clear();
460  if (!occupancyPFCand_puppiNolepWeight_.empty())
461  occupancyPFCand_puppiNolepWeight_.clear();
462  if (!occupancyPFCand_name_puppiNolepWeight_.empty())
463  occupancyPFCand_name_puppiNolepWeight_.clear();
464  if (!etaMinPFCand_.empty())
465  etaMinPFCand_.clear();
466  if (!etaMaxPFCand_.empty())
467  etaMaxPFCand_.clear();
468  if (!typePFCand_.empty())
469  typePFCand_.clear();
470  if (!countsPFCand_.empty())
471  countsPFCand_.clear();
472  if (!ptPFCand_.empty())
473  ptPFCand_.clear();
474  if (!ptPFCand_name_.empty())
475  ptPFCand_name_.clear();
476  if (!ptPFCand_puppiNolepWeight_.empty())
477  ptPFCand_puppiNolepWeight_.clear();
478  if (!ptPFCand_name_puppiNolepWeight_.empty())
479  ptPFCand_name_puppiNolepWeight_.clear();
480  if (!multiplicityPFCand_.empty())
481  multiplicityPFCand_.clear();
482  if (!multiplicityPFCand_name_.empty())
483  multiplicityPFCand_name_.clear();
484  for (std::vector<edm::ParameterSet>::const_iterator v = diagnosticsParameters_.begin();
485  v != diagnosticsParameters_.end();
486  v++) {
487  int etaNBinsPFCand = v->getParameter<int>("etaNBins");
488  double etaMinPFCand = v->getParameter<double>("etaMin");
489  double etaMaxPFCand = v->getParameter<double>("etaMax");
490  int phiNBinsPFCand = v->getParameter<int>("phiNBins");
491  double phiMinPFCand = v->getParameter<double>("phiMin");
492  double phiMaxPFCand = v->getParameter<double>("phiMax");
493  int nMinPFCand = v->getParameter<int>("nMin");
494  int nMaxPFCand = v->getParameter<int>("nMax");
495  int nbinsPFCand = v->getParameter<double>("nbins");
496 
497  // etaNBins_.push_back(etaNBins);
498  etaMinPFCand_.push_back(etaMinPFCand);
499  etaMaxPFCand_.push_back(etaMaxPFCand);
500  typePFCand_.push_back(v->getParameter<int>("type"));
501  countsPFCand_.push_back(0);
502 
503  multiplicityPFCand_.push_back(
504  ibooker.book1D(std::string(v->getParameter<std::string>("name")).append("_multiplicity_").c_str(),
505  std::string(v->getParameter<std::string>("name")) + "multiplicity",
506  nbinsPFCand,
507  nMinPFCand,
508  nMaxPFCand));
509  multiplicityPFCand_name_.push_back(std::string(v->getParameter<std::string>("name")).append("_multiplicity_"));
510  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
511  DirName + "/" + multiplicityPFCand_name_[multiplicityPFCand_name_.size() - 1],
512  multiplicityPFCand_[multiplicityPFCand_.size() - 1]));
513  //push back names first, we need to create histograms with the name and fill it for endcap plots later
514  occupancyPFCand_name_.push_back(
515  std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiWeight_"));
516  ptPFCand_name_.push_back(std::string(v->getParameter<std::string>("name")).append("_pt_puppiWeight_"));
517  //push back names first, we need to create histograms with the name and fill it for endcap plots later
518  occupancyPFCand_name_puppiNolepWeight_.push_back(
519  std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiNolepWeight_"));
520  ptPFCand_name_puppiNolepWeight_.push_back(
521  std::string(v->getParameter<std::string>("name")).append("_pt_puppiNolepWeight_"));
522  //special booking for endcap plots, merge plots for eminus and eplus into one plot, using variable binning
523  //barrel plots have eta-boundaries on minus and plus side
524  //parameters start on minus side
525  if (etaMinPFCand * etaMaxPFCand < 0) { //barrel plots, plot only in barrel region
526  occupancyPFCand_.push_back(
527  ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiWeight_").c_str(),
528  std::string(v->getParameter<std::string>("name")) + "occupancy_puppiWeight_",
529  etaNBinsPFCand,
530  etaMinPFCand,
531  etaMaxPFCand,
532  phiNBinsPFCand,
533  phiMinPFCand,
534  phiMaxPFCand));
535  ptPFCand_.push_back(
536  ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_puppiWeight_").c_str(),
537  std::string(v->getParameter<std::string>("name")) + "pt_puppiWeight_",
538  etaNBinsPFCand,
539  etaMinPFCand,
540  etaMaxPFCand,
541  phiNBinsPFCand,
542  phiMinPFCand,
543  phiMaxPFCand));
544  occupancyPFCand_puppiNolepWeight_.push_back(ibooker.book2D(
545  std::string(v->getParameter<std::string>("name")).append("_occupancy_puppiNolepWeight_").c_str(),
546  std::string(v->getParameter<std::string>("name")) + "occupancy_puppiNolepWeight",
547  etaNBinsPFCand,
548  etaMinPFCand,
549  etaMaxPFCand,
550  phiNBinsPFCand,
551  phiMinPFCand,
552  phiMaxPFCand));
553  ptPFCand_puppiNolepWeight_.push_back(
554  ibooker.book2D(std::string(v->getParameter<std::string>("name")).append("_pt_puppiNolepWeight_").c_str(),
555  std::string(v->getParameter<std::string>("name")) + "pt_puppiNolepWeight",
556  etaNBinsPFCand,
557  etaMinPFCand,
558  etaMaxPFCand,
559  phiNBinsPFCand,
560  phiMinPFCand,
561  phiMaxPFCand));
562  } else { //endcap or forward plots,
563  const int nbins_eta_endcap = 2 * (etaNBinsPFCand + 1);
564  double eta_limits_endcap[nbins_eta_endcap];
565  for (int i = 0; i < nbins_eta_endcap; i++) {
566  if (i < (etaNBinsPFCand + 1)) {
567  eta_limits_endcap[i] = etaMinPFCand + i * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
568  } else {
569  eta_limits_endcap[i] =
570  -etaMaxPFCand + (i - (etaNBinsPFCand + 1)) * (etaMaxPFCand - etaMinPFCand) / (double)etaNBinsPFCand;
571  }
572  }
573  TH2F* hist_temp_occup = new TH2F((occupancyPFCand_name_[occupancyPFCand_name_.size() - 1]).c_str(),
574  (occupancyPFCand_name_[occupancyPFCand_name_.size() - 1]).c_str(),
575  nbins_eta_endcap - 1,
576  eta_limits_endcap,
577  phiNBinsPFCand,
578  phiMinPFCand,
579  phiMaxPFCand);
580  occupancyPFCand_.push_back(
581  ibooker.book2D(occupancyPFCand_name_[occupancyPFCand_name_.size() - 1], hist_temp_occup));
582  TH2F* hist_temp_pt = new TH2F((ptPFCand_name_[ptPFCand_name_.size() - 1]).c_str(),
583  (ptPFCand_name_[ptPFCand_name_.size() - 1]).c_str(),
584  nbins_eta_endcap - 1,
585  eta_limits_endcap,
586  phiNBinsPFCand,
587  phiMinPFCand,
588  phiMaxPFCand);
589  ptPFCand_.push_back(ibooker.book2D(ptPFCand_name_[ptPFCand_name_.size() - 1], hist_temp_pt));
590  TH2F* hist_temp_occup_puppiNolepWeight = new TH2F(
591  (occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
592  (occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
593  nbins_eta_endcap - 1,
594  eta_limits_endcap,
595  phiNBinsPFCand,
596  phiMinPFCand,
597  phiMaxPFCand);
598  occupancyPFCand_puppiNolepWeight_.push_back(
599  ibooker.book2D(occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1],
600  hist_temp_occup_puppiNolepWeight));
601  TH2F* hist_temp_pt_puppiNolepWeight =
602  new TH2F((ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
603  (ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1]).c_str(),
604  nbins_eta_endcap - 1,
605  eta_limits_endcap,
606  phiNBinsPFCand,
607  phiMinPFCand,
608  phiMaxPFCand);
609  ptPFCand_puppiNolepWeight_.push_back(
610  ibooker.book2D(ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1],
611  hist_temp_pt_puppiNolepWeight));
612  }
613  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
614  DirName + "/" + occupancyPFCand_name_puppiNolepWeight_[occupancyPFCand_name_puppiNolepWeight_.size() - 1],
615  occupancyPFCand_puppiNolepWeight_[occupancyPFCand_puppiNolepWeight_.size() - 1]));
616  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
617  DirName + "/" + ptPFCand_name_puppiNolepWeight_[ptPFCand_name_puppiNolepWeight_.size() - 1],
618  ptPFCand_puppiNolepWeight_[ptPFCand_puppiNolepWeight_.size() - 1]));
619  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
620  DirName + "/" + occupancyPFCand_name_[occupancyPFCand_name_.size() - 1],
621  occupancyPFCand_[occupancyPFCand_.size() - 1]));
622  map_of_MEs.insert(std::pair<std::string, MonitorElement*>(
623  DirName + "/" + ptPFCand_name_[ptPFCand_name_.size() - 1], ptPFCand_[ptPFCand_.size() - 1]));
624  }
625  }
626 }

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

◆ 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.

629  {
630  miniaodfilterindex = -1;
631 
632  if (isMiniAOD_) {
633  bool changed_filter = true;
634  if (FilterhltConfig_.init(iRun, iSetup, METFilterMiniAODLabel_.process(), changed_filter)) {
635  miniaodfilterdec = 0;
636  for (unsigned int i = 0; i < FilterhltConfig_.size(); i++) {
637  std::string search = FilterhltConfig_.triggerName(i).substr(
638  5); //actual label of filter, the first 5 items are Flag_, so stripped off
639  std::string search2 =
640  HBHENoiseStringMiniAOD; //all filters end with DQM, which is not in the flag --> ONLY not for HBHEFilters
641  std::size_t found = search2.find(search);
642  if (found != std::string::npos) {
643  miniaodfilterindex = i;
644  }
645  }
646  } else if (FilterhltConfig_.init(iRun, iSetup, METFilterMiniAODLabel2_.process(), changed_filter)) {
647  miniaodfilterdec = 1;
648  for (unsigned int i = 0; i < FilterhltConfig_.size(); i++) {
649  std::string search = FilterhltConfig_.triggerName(i).substr(
650  5); //actual label of filter, the first 5 items are Flag_, so stripped off
651  std::string search2 =
652  HBHENoiseStringMiniAOD; //all filters end with DQM, which is not in the flag --> ONLY not for HBHEFilters
653  std::size_t found = search2.find(search);
654  if (found != std::string::npos) {
655  miniaodfilterindex = i;
656  }
657  }
658  } else {
659  edm::LogWarning("MiniAOD MET Filter HLT OBject version")
660  << "nothing found with both RECO and reRECO label" << std::endl;
661  }
662  }
663 }

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

◆ jetSortingRule()

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

Definition at line 79 of file DQMPFCandidateAnalyzer.h.

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

References x, and y.

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.

DQMPFCandidateAnalyzer::METFilterMiniAODToken2_
edm::EDGetTokenT< edm::TriggerResults > METFilterMiniAODToken2_
Definition: DQMPFCandidateAnalyzer.h:96
JetMETDQMDCSFilter
Definition: JetMETDQMDCSFilter.h:24
DQMPFCandidateAnalyzer::miniaodfilterdec
int miniaodfilterdec
Definition: DQMPFCandidateAnalyzer.h:101
DDAxes::y
DQMPFCandidateAnalyzer::mInputCollection_
edm::InputTag mInputCollection_
Definition: DQMPFCandidateAnalyzer.h:86
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_TrackPt
MonitorElement * mProfileIsoPFChHad_TrackPt
Definition: DQMPFCandidateAnalyzer.h:146
mps_fire.i
i
Definition: mps_fire.py:428
DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_1_10
MonitorElement * m_HOverTrackP_EndCap_hPt_1_10
Definition: DQMPFCandidateAnalyzer.h:137
dqmMemoryStats.float
float
Definition: dqmMemoryStats.py:127
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EMPtCentral
MonitorElement * mProfileIsoPFChHad_EMPtCentral
Definition: DQMPFCandidateAnalyzer.h:144
JetMETDQMDCSFilter::filter
bool filter(const edm::Event &evt, const edm::EventSetup &es)
Definition: JetMETDQMDCSFilter.cc:63
DQMPFCandidateAnalyzer::occupancyPFCandRECO_
std::vector< MonitorElement * > occupancyPFCandRECO_
Definition: DQMPFCandidateAnalyzer.h:163
edm::Handle::product
T const * product() const
Definition: Handle.h:70
DQMPFCandidateAnalyzer::numPV_
int numPV_
Definition: DQMPFCandidateAnalyzer.h:168
DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_1_10
MonitorElement * m_HOverTrackPVsEta_hPt_1_10
Definition: DQMPFCandidateAnalyzer.h:127
DQMPFCandidateAnalyzer::ptMinCand_
double ptMinCand_
Definition: DQMPFCandidateAnalyzer.h:116
DQMPFCandidateAnalyzer::pflowToken_
edm::EDGetTokenT< std::vector< reco::PFCandidate > > pflowToken_
Definition: DQMPFCandidateAnalyzer.h:103
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
gather_cfg.cout
cout
Definition: gather_cfg.py:144
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89301
DQMPFCandidateAnalyzer::etaMaxPFCand_
std::vector< double > etaMaxPFCand_
Definition: DQMPFCandidateAnalyzer.h:159
DQMPFCandidateAnalyzer::occupancyPFCand_
std::vector< MonitorElement * > occupancyPFCand_
Definition: DQMPFCandidateAnalyzer.h:155
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EcalOccupancyCentral
MonitorElement * mProfileIsoPFChHad_EcalOccupancyCentral
Definition: DQMPFCandidateAnalyzer.h:150
DQMPFCandidateAnalyzer::DCSFilter_
JetMETDQMDCSFilter * DCSFilter_
Definition: DQMPFCandidateAnalyzer.h:111
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:47
edm::InputTag::process
std::string const & process() const
Definition: InputTag.h:40
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DQMPFCandidateAnalyzer::etaMinPFCand_
std::vector< double > etaMinPFCand_
Definition: DQMPFCandidateAnalyzer.h:159
DDAxes::x
DQMPFCandidateAnalyzer::miniaodfilterindex
int miniaodfilterindex
Definition: DQMPFCandidateAnalyzer.h:99
findQualityFiles.v
v
Definition: findQualityFiles.py:179
DQMPFCandidateAnalyzer::LSBegin_
int LSBegin_
Definition: DQMPFCandidateAnalyzer.h:171
DQMPFCandidateAnalyzer::etaMaxPFCandRECO_
std::vector< double > etaMaxPFCandRECO_
Definition: DQMPFCandidateAnalyzer.h:165
BeamSpotPI::nTracks
Definition: BeamSpotPayloadInspectorHelper.h:42
newFWLiteAna.found
found
Definition: newFWLiteAna.py:118
edm::Handle< VertexCollection >
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
DQMPFCandidateAnalyzer::occupancyPFCand_nameRECO_
std::vector< std::string > occupancyPFCand_nameRECO_
Definition: DQMPFCandidateAnalyzer.h:164
DQMPFCandidateAnalyzer::map_of_MEs
std::map< std::string, MonitorElement * > map_of_MEs
Definition: DQMPFCandidateAnalyzer.h:177
edm::Ref< PFBlockCollection >
reco::PFBlockElement::Type
Type
Definition: PFBlockElement.h:30
DQMPFCandidateAnalyzer::occupancyPFCand_puppiNolepWeight_
std::vector< MonitorElement * > occupancyPFCand_puppiNolepWeight_
Definition: DQMPFCandidateAnalyzer.h:157
DQMPFCandidateAnalyzer::ptPFCand_name_
std::vector< std::string > ptPFCand_name_
Definition: DQMPFCandidateAnalyzer.h:156
DQMPFCandidateAnalyzer::occupancyPFCand_name_puppiNolepWeight_
std::vector< std::string > occupancyPFCand_name_puppiNolepWeight_
Definition: DQMPFCandidateAnalyzer.h:158
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
DQMPFCandidateAnalyzer::ptPFCand_nameRECO_
std::vector< std::string > ptPFCand_nameRECO_
Definition: DQMPFCandidateAnalyzer.h:164
ECAL
Definition: HCALResponse.h:21
DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_20_50
MonitorElement * m_HOverTrackP_Barrel_hPt_20_50
Definition: DQMPFCandidateAnalyzer.h:134
DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackP_EndCap
MonitorElement * m_HOverTrackPVsTrackP_EndCap
Definition: DQMPFCandidateAnalyzer.h:122
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HadPtCentral
MonitorElement * mProfileIsoPFChHad_HadPtCentral
Definition: DQMPFCandidateAnalyzer.h:142
DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackPt_EndCap
MonitorElement * m_HOverTrackPVsTrackPt_EndCap
Definition: DQMPFCandidateAnalyzer.h:125
DQMPFCandidateAnalyzer::FilterhltConfig_
HLTConfigProvider FilterhltConfig_
Definition: DQMPFCandidateAnalyzer.h:92
DQMPFCandidateAnalyzer::bypassAllDCSChecks_
bool bypassAllDCSChecks_
Definition: DQMPFCandidateAnalyzer.h:175
DQMPFCandidateAnalyzer::ptPFCand_puppiNolepWeight_
std::vector< MonitorElement * > ptPFCand_puppiNolepWeight_
Definition: DQMPFCandidateAnalyzer.h:157
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
DQMPFCandidateAnalyzer::multiplicityPFCand_nameRECO_
std::vector< std::string > multiplicityPFCand_nameRECO_
Definition: DQMPFCandidateAnalyzer.h:164
DQMPFCandidateAnalyzer::ptPFCand_name_puppiNolepWeight_
std::vector< std::string > ptPFCand_name_puppiNolepWeight_
Definition: DQMPFCandidateAnalyzer.h:158
dqm::implementation::IBooker::bookProfile
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int, double lowY, double highY, char const *option="s", FUNC onbooking=NOOP())
Definition: DQMStore.h:322
DQMPFCandidateAnalyzer::bypassAllPVChecks_
bool bypassAllPVChecks_
Definition: DQMPFCandidateAnalyzer.h:174
DQMPFCandidateAnalyzer::multiplicityPFCand_
std::vector< MonitorElement * > multiplicityPFCand_
Definition: DQMPFCandidateAnalyzer.h:155
DQMPFCandidateAnalyzer::HBHENoiseStringMiniAOD
std::string HBHENoiseStringMiniAOD
Definition: DQMPFCandidateAnalyzer.h:90
pfLinker_cff.particleFlow
particleFlow
Definition: pfLinker_cff.py:5
cond::persistency::search
std::vector< T >::const_iterator search(const cond::Time_t &val, const std::vector< T > &container)
Definition: IOVProxy.cc:21
HCAL
Definition: HCALResponse.h:21
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_TrackOccupancy
MonitorElement * mProfileIsoPFChHad_TrackOccupancy
Definition: DQMPFCandidateAnalyzer.h:152
reco::PFCandidate::ElementsInBlocks
std::vector< ElementInBlock > ElementsInBlocks
Definition: PFCandidate.h:379
electronDataDiscovery.search2
def search2()
Definition: electronDataDiscovery.py:280
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:63
DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_20_50
MonitorElement * m_HOverTrackPVsEta_hPt_20_50
Definition: DQMPFCandidateAnalyzer.h:129
LepHTMonitor_cff.nbins_eta
nbins_eta
Definition: LepHTMonitor_cff.py:41
DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_50
MonitorElement * m_HOverTrackP_Barrel_hPt_50
Definition: DQMPFCandidateAnalyzer.h:135
DQMPFCandidateAnalyzer::ptPFCandRECO_
std::vector< MonitorElement * > ptPFCandRECO_
Definition: DQMPFCandidateAnalyzer.h:163
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233
edm::ParameterSet
Definition: ParameterSet.h:47
DQMPFCandidateAnalyzer::cleaningParameters_
edm::ParameterSet cleaningParameters_
Definition: DQMPFCandidateAnalyzer.h:113
DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_10_20
MonitorElement * m_HOverTrackPVsEta_hPt_10_20
Definition: DQMPFCandidateAnalyzer.h:128
DQMPFCandidateAnalyzer::typePFCandRECO_
std::vector< int > typePFCandRECO_
Definition: DQMPFCandidateAnalyzer.h:166
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_20_50
MonitorElement * m_HOverTrackP_EndCap_hPt_20_50
Definition: DQMPFCandidateAnalyzer.h:139
pat::PackedCandidate
Definition: PackedCandidate.h:22
DQMPFCandidateAnalyzer::multiplicityPFCandRECO_
std::vector< MonitorElement * > multiplicityPFCandRECO_
Definition: DQMPFCandidateAnalyzer.h:163
DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_10_20
MonitorElement * m_HOverTrackP_EndCap_hPt_10_20
Definition: DQMPFCandidateAnalyzer.h:138
DQMPFCandidateAnalyzer::hcalMin_
double hcalMin_
Definition: DQMPFCandidateAnalyzer.h:118
iEvent
int iEvent
Definition: GenABIO.cc:224
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
DQMPFCandidateAnalyzer::LSEnd_
int LSEnd_
Definition: DQMPFCandidateAnalyzer.h:172
DQMPFCandidateAnalyzer::METFilterMiniAODToken_
edm::EDGetTokenT< edm::TriggerResults > METFilterMiniAODToken_
Definition: DQMPFCandidateAnalyzer.h:94
DQMPFCandidateAnalyzer::isMiniAOD_
bool isMiniAOD_
Definition: DQMPFCandidateAnalyzer.h:179
EgHLTOffHistBins_cfi.et
et
Definition: EgHLTOffHistBins_cfi.py:8
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EMPtEndcap
MonitorElement * mProfileIsoPFChHad_EMPtEndcap
Definition: DQMPFCandidateAnalyzer.h:145
DQMPFCandidateAnalyzer::m_HOverTrackPVsEta_hPt_50
MonitorElement * m_HOverTrackPVsEta_hPt_50
Definition: DQMPFCandidateAnalyzer.h:130
DQMPFCandidateAnalyzer::hbheNoiseFilterResultToken_
edm::EDGetTokenT< bool > hbheNoiseFilterResultToken_
Definition: DQMPFCandidateAnalyzer.h:89
DQMPFCandidateAnalyzer::m_HOverTrackP_EndCap_hPt_50
MonitorElement * m_HOverTrackP_EndCap_hPt_50
Definition: DQMPFCandidateAnalyzer.h:140
DQMPFCandidateAnalyzer::etaMinPFCandRECO_
std::vector< double > etaMinPFCandRECO_
Definition: DQMPFCandidateAnalyzer.h:165
DQMPFCandidateAnalyzer::typePFCand_
std::vector< int > typePFCand_
Definition: DQMPFCandidateAnalyzer.h:160
DQMPFCandidateAnalyzer::hbheNoiseFilterResultTag_
edm::InputTag hbheNoiseFilterResultTag_
Definition: DQMPFCandidateAnalyzer.h:88
DQMPFCandidateAnalyzer::countsPFCandRECO_
std::vector< int > countsPFCandRECO_
Definition: DQMPFCandidateAnalyzer.h:166
DQMPFCandidateAnalyzer::diagnosticsParameters_
std::vector< edm::ParameterSet > diagnosticsParameters_
Definition: DQMPFCandidateAnalyzer.h:114
DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackPt_Barrel
MonitorElement * m_HOverTrackPVsTrackPt_Barrel
Definition: DQMPFCandidateAnalyzer.h:124
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DQMPFCandidateAnalyzer::candidateType_
std::string candidateType_
Definition: DQMPFCandidateAnalyzer.h:106
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:1053
dqm::impl::MonitorElement::getRootObject
virtual const TObject * getRootObject() const
Definition: MonitorElement.cc:967
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HcalOccupancyEndcap
MonitorElement * mProfileIsoPFChHad_HcalOccupancyEndcap
Definition: DQMPFCandidateAnalyzer.h:149
bookConverter.elements
elements
Definition: bookConverter.py:147
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HadPtEndcap
MonitorElement * mProfileIsoPFChHad_HadPtEndcap
Definition: DQMPFCandidateAnalyzer.h:143
DQMPFCandidateAnalyzer::verbose_
int verbose_
Definition: DQMPFCandidateAnalyzer.h:169
DQMPFCandidateAnalyzer::vertexToken_
edm::EDGetTokenT< std::vector< reco::Vertex > > vertexToken_
Definition: DQMPFCandidateAnalyzer.h:84
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d'tor
Definition: HLTConfigProvider.cc:36
DQMPFCandidateAnalyzer::multiplicityPFCand_name_
std::vector< std::string > multiplicityPFCand_name_
Definition: DQMPFCandidateAnalyzer.h:156
spclusmultinvestigator_cfi.vertexCollection
vertexCollection
Definition: spclusmultinvestigator_cfi.py:4
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:177
reco::PFBlockElementCluster
Cluster Element.
Definition: PFBlockElementCluster.h:16
DQMPFCandidateAnalyzer::ptPFCand_
std::vector< MonitorElement * > ptPFCand_
Definition: DQMPFCandidateAnalyzer.h:155
reco::PFBlockElementTrack
Track Element.
Definition: PFBlockElementTrack.h:17
DQMPFCandidateAnalyzer::METFilterMiniAODLabel2_
edm::InputTag METFilterMiniAODLabel2_
Definition: DQMPFCandidateAnalyzer.h:95
fed_dqm_sourceclient-live_cfg.DirName
DirName
Definition: fed_dqm_sourceclient-live_cfg.py:56
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
reco::PFBlockElementCluster::clusterRef
const PFClusterRef & clusterRef() const override
Definition: PFBlockElementCluster.h:29
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
edm::HLTGlobalStatus::accept
bool accept() const
Has at least one path accepted the event?
Definition: HLTGlobalStatus.h:49
reco::PFCandidate
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:41
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:70
DQMPFCandidateAnalyzer::hbhenoifilterdecision
bool hbhenoifilterdecision
Definition: DQMPFCandidateAnalyzer.h:98
DQMPFCandidateAnalyzer::m_HOverTrackP_trackPtVsEta
MonitorElement * m_HOverTrackP_trackPtVsEta
Definition: DQMPFCandidateAnalyzer.h:120
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:234
c
auto & c
Definition: CAHitNtupletGeneratorKernelsImpl.h:56
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
DQMPFCandidateAnalyzer::m_HOverTrackPVsTrackP_Barrel
MonitorElement * m_HOverTrackPVsTrackP_Barrel
Definition: DQMPFCandidateAnalyzer.h:121
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_EcalOccupancyEndcap
MonitorElement * mProfileIsoPFChHad_EcalOccupancyEndcap
Definition: DQMPFCandidateAnalyzer.h:151
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
DQMPFCandidateAnalyzer::vertexTag_
edm::InputTag vertexTag_
Definition: DQMPFCandidateAnalyzer.h:83
DQMPFCandidateAnalyzer::countsPFCand_
std::vector< int > countsPFCand_
Definition: DQMPFCandidateAnalyzer.h:160
cuy.ii
ii
Definition: cuy.py:589
DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_10_20
MonitorElement * m_HOverTrackP_Barrel_hPt_10_20
Definition: DQMPFCandidateAnalyzer.h:133
DQMPFCandidateAnalyzer::METFilterMiniAODLabel_
edm::InputTag METFilterMiniAODLabel_
Definition: DQMPFCandidateAnalyzer.h:93
edm::InputTag
Definition: InputTag.h:15
edm::ConsumesCollector
Definition: ConsumesCollector.h:45
DQMPFCandidateAnalyzer::m_HOverTrackP_Barrel_hPt_1_10
MonitorElement * m_HOverTrackP_Barrel_hPt_1_10
Definition: DQMPFCandidateAnalyzer.h:132
DQMPFCandidateAnalyzer::occupancyPFCand_name_
std::vector< std::string > occupancyPFCand_name_
Definition: DQMPFCandidateAnalyzer.h:156
DQMPFCandidateAnalyzer::pflowPackedToken_
edm::EDGetTokenT< std::vector< pat::PackedCandidate > > pflowPackedToken_
Definition: DQMPFCandidateAnalyzer.h:104
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
edm::OwnVector< reco::PFBlockElement >
DQMPFCandidateAnalyzer::mProfileIsoPFChHad_HcalOccupancyCentral
MonitorElement * mProfileIsoPFChHad_HcalOccupancyCentral
Definition: DQMPFCandidateAnalyzer.h:148
Generated for CMSSW Reference Manual by doxygen 1.8.16