CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
PatTriggerAnalyzer Class Reference
Inheritance diagram for PatTriggerAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 PatTriggerAnalyzer (const edm::ParameterSet &iConfig)
 default constructor More...
 
 ~PatTriggerAnalyzer () override
 default destructor More...
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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)
 

Private Member Functions

void analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup) override
 everythin that needs to be done during the event loop More...
 
void beginJob () override
 everythin that needs to be done before the event loop More...
 
void endJob () override
 everythin that needs to be done after the event loop More...
 

Private Attributes

std::map< std::string, TH1D * > histos1D_
 histogram management More...
 
std::map< std::string, TH2D * > histos2D_
 
unsigned maxID_
 maximal id for meanPt plot More...
 
unsigned minID_
 minimal id for meanPt plot More...
 
std::string muonMatch_
 input for trigger match objects More...
 
edm::EDGetTokenT< pat::MuonCollectionmuonsToken_
 input for muons More...
 
std::map< unsigned, unsigned > sumN_
 internals for meanPt histogram calculation More...
 
std::map< unsigned, double > sumPt_
 
edm::InputTag trigger_
 input for patTrigger More...
 
edm::EDGetTokenT< pat::TriggerEventtriggerEventToken_
 input for patTriggerEvent More...
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsInputProcessBlocks ()
 
static bool wantsProcessBlocks ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- 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)
 

Detailed Description

Definition at line 15 of file PatTriggerAnalyzer.cc.

Constructor & Destructor Documentation

◆ PatTriggerAnalyzer()

PatTriggerAnalyzer::PatTriggerAnalyzer ( const edm::ParameterSet iConfig)
explicit

default constructor

Definition at line 60 of file PatTriggerAnalyzer.cc.

61  : // pat::Trigger
62  trigger_(iConfig.getParameter<edm::InputTag>("trigger")),
63  // pat::TriggerEvent
64  triggerEventToken_(consumes<TriggerEvent>(iConfig.getParameter<edm::InputTag>("triggerEvent"))),
65  // muon input collection
66  muonsToken_(consumes<MuonCollection>(iConfig.getParameter<edm::InputTag>("muons"))),
67  // muon match objects
68  muonMatch_(iConfig.getParameter<std::string>("muonMatch")),
69  // minimal id for of all trigger objects
70  minID_(iConfig.getParameter<unsigned>("minID")),
71  // maximal id for of all trigger objects
72  maxID_(iConfig.getParameter<unsigned>("maxID")),
73  histos1D_(),
74  histos2D_() {}

◆ ~PatTriggerAnalyzer()

PatTriggerAnalyzer::~PatTriggerAnalyzer ( )
override

default destructor

Definition at line 76 of file PatTriggerAnalyzer.cc.

76 {}

Member Function Documentation

◆ analyze()

void PatTriggerAnalyzer::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

everythin that needs to be done during the event loop

Implements edm::EDAnalyzer.

Definition at line 113 of file PatTriggerAnalyzer.cc.

113  {
114  // PAT trigger event
117 
118  // PAT object collection
120  iEvent.getByToken(muonsToken_, muons);
121 
122  // PAT trigger helper for trigger matching information
123  const helper::TriggerMatchHelper matchHelper;
124 
125  /*
126  kinematics comparison
127  */
128 
129  // loop over muon references (PAT muons have been used in the matcher in task 3)
130  for (size_t iMuon = 0; iMuon < muons->size(); ++iMuon) {
131  // we need all these ingedients to recieve matched trigger objects from the matchHelper
132  const TriggerObjectRef trigRef(matchHelper.triggerMatchObject(muons, iMuon, muonMatch_, iEvent, *triggerEvent));
133  // finally we can fill the histograms
134  if (trigRef.isAvailable() && trigRef.isNonnull()) { // check references (necessary!)
135  histos2D_["ptTrigCand"]->Fill(muons->at(iMuon).pt(), trigRef->pt());
136  histos2D_["etaTrigCand"]->Fill(muons->at(iMuon).eta(), trigRef->eta());
137  histos2D_["phiTrigCand"]->Fill(muons->at(iMuon).phi(), trigRef->phi());
138  }
139  }
140 
141  /*
142  turn-on curve
143  */
144 
145  // get the trigger objects corresponding to the used matching (HLT muons)
146  const TriggerObjectRefVector trigRefs(triggerEvent->objects(trigger::TriggerMuon));
147  // loop over selected trigger objects
148  for (TriggerObjectRefVector::const_iterator iTrig = trigRefs.begin(); iTrig != trigRefs.end(); ++iTrig) {
149  // get all matched candidates for the trigger object
150  const reco::CandidateBaseRefVector candRefs(
151  matchHelper.triggerMatchCandidates((*iTrig), muonMatch_, iEvent, *triggerEvent));
152  if (candRefs.empty())
153  continue;
154  // fill the histogram...
155  // (only for the first match, since we resolved ambiguities in the matching configuration,
156  // so that we have one at maximum per trigger object)
157  reco::CandidateBaseRef muonRef(candRefs.at(0));
158  if (muonRef.isAvailable() && muonRef.isNonnull()) {
159  histos1D_["turnOn"]->Fill(muonRef->pt());
160  }
161  }
162 
163  /*
164  mean pt
165  */
166 
167  // loop over all trigger match objects from minID to maxID; have
168  // a look to DataFormats/HLTReco/interface/TriggerTypeDefs.h to
169  // know more about the available trrigger object id's
170  for (unsigned id = minID_; id <= maxID_; ++id) {
171  // vector of all objects for a given object id
172  const TriggerObjectRefVector objRefs(triggerEvent->objects(id));
173  // buffer the number of objects
174  sumN_[id] += objRefs.size();
175  // iterate the objects and buffer the pt of the objects
176  for (TriggerObjectRefVector::const_iterator iRef = objRefs.begin(); iRef != objRefs.end(); ++iRef) {
177  sumPt_[id] += (*iRef)->pt();
178  }
179  }
180 }

References histos1D_, histos2D_, triggerObjects_cff::id, iEvent, maxID_, minID_, muonMatch_, PDWG_BPHSkim_cff::muons, muonsToken_, sumN_, sumPt_, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEventToken_, and trigger::TriggerMuon.

◆ beginJob()

void PatTriggerAnalyzer::beginJob ( void  )
overrideprivatevirtual

everythin that needs to be done before the event loop

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file PatTriggerAnalyzer.cc.

78  {
79  edm::Service<TFileService> fileService;
80 
81  // pt correlation plot
82  histos2D_["ptTrigCand"] =
83  fileService->make<TH2D>("ptTrigCand", "Object vs. candidate p_{T} (GeV)", 60, 0., 300., 60, 0., 300.);
84  histos2D_["ptTrigCand"]->SetXTitle("candidate p_{T} (GeV)");
85  histos2D_["ptTrigCand"]->SetYTitle("object p_{T} (GeV)");
86  // eta correlation plot
87  histos2D_["etaTrigCand"] =
88  fileService->make<TH2D>("etaTrigCand", "Object vs. candidate #eta", 50, -2.5, 2.5, 50, -2.5, 2.5);
89  histos2D_["etaTrigCand"]->SetXTitle("candidate #eta");
90  histos2D_["etaTrigCand"]->SetYTitle("object #eta");
91  // phi correlation plot
92  histos2D_["phiTrigCand"] = fileService->make<TH2D>(
93  "phiTrigCand", "Object vs. candidate #phi", 60, -TMath::Pi(), TMath::Pi(), 60, -TMath::Pi(), TMath::Pi());
94  histos2D_["phiTrigCand"]->SetXTitle("candidate #phi");
95  histos2D_["phiTrigCand"]->SetYTitle("object #phi");
96  // turn-on curves
97  histos1D_["turnOn"] = fileService->make<TH1D>("turnOn", "p_{T} (GeV) of matched candidate", 10, 0., 50.);
98  histos1D_["turnOn"]->SetXTitle("candidate p_{T} (GeV)");
99  histos1D_["turnOn"]->SetYTitle("# of objects");
100  // mean pt for all trigger objects
101  histos1D_["ptMean"] = fileService->make<TH1D>(
102  "ptMean", "Mean p_{T} (GeV) per trigger object type", maxID_ - minID_ + 1, minID_ - 0.5, maxID_ + 0.5);
103  histos1D_["ptMean"]->SetXTitle("trigger object type");
104  histos1D_["ptMean"]->SetYTitle("mean p_{T} (GeV)");
105 
106  // initialize counters for mean pt calculation
107  for (unsigned id = minID_; id <= maxID_; ++id) {
108  sumN_[id] = 0;
109  sumPt_[id] = 0.;
110  }
111 }

References histos1D_, histos2D_, triggerObjects_cff::id, TFileService::make(), maxID_, minID_, Pi, sumN_, and sumPt_.

◆ endJob()

void PatTriggerAnalyzer::endJob ( void  )
overrideprivatevirtual

everythin that needs to be done after the event loop

Reimplemented from edm::EDAnalyzer.

Definition at line 182 of file PatTriggerAnalyzer.cc.

182  {
183  // normalize the entries for the mean pt plot
184  for (unsigned id = minID_; id <= maxID_; ++id) {
185  if (sumN_[id] != 0)
186  histos1D_["ptMean"]->Fill(id, sumPt_[id] / sumN_[id]);
187  }
188 }

References histos1D_, triggerObjects_cff::id, maxID_, minID_, sumN_, and sumPt_.

Member Data Documentation

◆ histos1D_

std::map<std::string, TH1D*> PatTriggerAnalyzer::histos1D_
private

histogram management

Definition at line 44 of file PatTriggerAnalyzer.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ histos2D_

std::map<std::string, TH2D*> PatTriggerAnalyzer::histos2D_
private

Definition at line 45 of file PatTriggerAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ maxID_

unsigned PatTriggerAnalyzer::maxID_
private

maximal id for meanPt plot

Definition at line 41 of file PatTriggerAnalyzer.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ minID_

unsigned PatTriggerAnalyzer::minID_
private

minimal id for meanPt plot

Definition at line 39 of file PatTriggerAnalyzer.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ muonMatch_

std::string PatTriggerAnalyzer::muonMatch_
private

input for trigger match objects

Definition at line 37 of file PatTriggerAnalyzer.cc.

Referenced by analyze().

◆ muonsToken_

edm::EDGetTokenT<pat::MuonCollection> PatTriggerAnalyzer::muonsToken_
private

input for muons

Definition at line 35 of file PatTriggerAnalyzer.cc.

Referenced by analyze().

◆ sumN_

std::map<unsigned, unsigned> PatTriggerAnalyzer::sumN_
private

internals for meanPt histogram calculation

Definition at line 48 of file PatTriggerAnalyzer.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ sumPt_

std::map<unsigned, double> PatTriggerAnalyzer::sumPt_
private

Definition at line 49 of file PatTriggerAnalyzer.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ trigger_

edm::InputTag PatTriggerAnalyzer::trigger_
private

input for patTrigger

Definition at line 31 of file PatTriggerAnalyzer.cc.

◆ triggerEventToken_

edm::EDGetTokenT<pat::TriggerEvent> PatTriggerAnalyzer::triggerEventToken_
private

input for patTriggerEvent

Definition at line 33 of file PatTriggerAnalyzer.cc.

Referenced by analyze().

PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
PatTriggerAnalyzer::histos1D_
std::map< std::string, TH1D * > histos1D_
histogram management
Definition: PatTriggerAnalyzer.cc:44
PatTriggerAnalyzer::sumN_
std::map< unsigned, unsigned > sumN_
internals for meanPt histogram calculation
Definition: PatTriggerAnalyzer.cc:48
PatTriggerAnalyzer::sumPt_
std::map< unsigned, double > sumPt_
Definition: PatTriggerAnalyzer.cc:49
edm::RefVector
Definition: EDProductfwd.h:27
edm::Handle
Definition: AssociativeIterator.h:50
edm::Ref
Definition: AssociativeIterator.h:58
trigger::TriggerMuon
Definition: TriggerTypeDefs.h:79
PatTriggerAnalyzer::histos2D_
std::map< std::string, TH2D * > histos2D_
Definition: PatTriggerAnalyzer.cc:45
PatTriggerAnalyzer::muonMatch_
std::string muonMatch_
input for trigger match objects
Definition: PatTriggerAnalyzer.cc:37
PatTriggerAnalyzer::trigger_
edm::InputTag trigger_
input for patTrigger
Definition: PatTriggerAnalyzer.cc:31
edm::Service< TFileService >
iEvent
int iEvent
Definition: GenABIO.cc:224
PDWG_DiPhoton_SD_cff.triggerEvent
triggerEvent
Definition: PDWG_DiPhoton_SD_cff.py:39
edm::RefToBaseVector
Definition: EDProductfwd.h:29
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
PatTriggerAnalyzer::triggerEventToken_
edm::EDGetTokenT< pat::TriggerEvent > triggerEventToken_
input for patTriggerEvent
Definition: PatTriggerAnalyzer.cc:33
PatTriggerAnalyzer::maxID_
unsigned maxID_
maximal id for meanPt plot
Definition: PatTriggerAnalyzer.cc:41
PatTriggerAnalyzer::muonsToken_
edm::EDGetTokenT< pat::MuonCollection > muonsToken_
input for muons
Definition: PatTriggerAnalyzer.cc:35
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:29
edm::RefVectorIterator
Definition: EDProductfwd.h:33
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::RefToBase< Candidate >
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
edm::InputTag
Definition: InputTag.h:15
PatTriggerAnalyzer::minID_
unsigned minID_
minimal id for meanPt plot
Definition: PatTriggerAnalyzer.cc:39