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 const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
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 wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 15 of file PatTriggerAnalyzer.cc.

Constructor & Destructor Documentation

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_() {}
T getParameter(std::string const &) const
std::map< std::string, TH2D * > histos2D_
std::map< std::string, TH1D * > histos1D_
histogram management
unsigned minID_
minimal id for meanPt plot
edm::EDGetTokenT< pat::MuonCollection > muonsToken_
input for muons
edm::InputTag trigger_
input for patTrigger
unsigned maxID_
maximal id for meanPt plot
std::string muonMatch_
input for trigger match objects
edm::EDGetTokenT< pat::TriggerEvent > triggerEventToken_
input for patTriggerEvent
PatTriggerAnalyzer::~PatTriggerAnalyzer ( )
override

default destructor

Definition at line 76 of file PatTriggerAnalyzer.cc.

76 {}

Member Function Documentation

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

everythin that needs to be done during the event loop

Definition at line 113 of file PatTriggerAnalyzer.cc.

References edm::Event::getByToken(), histos1D_, histos2D_, triggerObjects_cff::id, maxID_, minID_, muonMatch_, PDWG_BPHSkim_cff::muons, muonsToken_, pat::TriggerEvent::objects(), sumN_, sumPt_, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEventToken_, and trigger::TriggerMuon.

113  {
114  // PAT trigger event
116  iEvent.getByToken(triggerEventToken_, triggerEvent);
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 }
std::map< unsigned, unsigned > sumN_
internals for meanPt histogram calculation
std::map< std::string, TH2D * > histos2D_
std::map< unsigned, double > sumPt_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
std::map< std::string, TH1D * > histos1D_
histogram management
const TriggerObjectCollection * objects() const
Definition: TriggerEvent.h:294
unsigned minID_
minimal id for meanPt plot
edm::EDGetTokenT< pat::MuonCollection > muonsToken_
input for muons
unsigned maxID_
maximal id for meanPt plot
std::string muonMatch_
input for trigger match objects
edm::EDGetTokenT< pat::TriggerEvent > triggerEventToken_
input for patTriggerEvent
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.

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

78  {
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 }
std::map< unsigned, unsigned > sumN_
internals for meanPt histogram calculation
const double Pi
std::map< std::string, TH2D * > histos2D_
std::map< unsigned, double > sumPt_
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
std::map< std::string, TH1D * > histos1D_
histogram management
unsigned minID_
minimal id for meanPt plot
edm::Service< TFileService > fileService
unsigned maxID_
maximal id for meanPt plot
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.

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

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 }
std::map< unsigned, unsigned > sumN_
internals for meanPt histogram calculation
std::map< unsigned, double > sumPt_
std::map< std::string, TH1D * > histos1D_
histogram management
unsigned minID_
minimal id for meanPt plot
unsigned maxID_
maximal id for meanPt plot

Member Data Documentation

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

histogram management

Definition at line 44 of file PatTriggerAnalyzer.cc.

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

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

Definition at line 45 of file PatTriggerAnalyzer.cc.

Referenced by analyze(), and beginJob().

unsigned PatTriggerAnalyzer::maxID_
private

maximal id for meanPt plot

Definition at line 41 of file PatTriggerAnalyzer.cc.

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

unsigned PatTriggerAnalyzer::minID_
private

minimal id for meanPt plot

Definition at line 39 of file PatTriggerAnalyzer.cc.

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

std::string PatTriggerAnalyzer::muonMatch_
private

input for trigger match objects

Definition at line 37 of file PatTriggerAnalyzer.cc.

Referenced by analyze().

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

input for muons

Definition at line 35 of file PatTriggerAnalyzer.cc.

Referenced by analyze().

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

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

Definition at line 49 of file PatTriggerAnalyzer.cc.

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

edm::InputTag PatTriggerAnalyzer::trigger_
private

input for patTrigger

Definition at line 31 of file PatTriggerAnalyzer.cc.

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

input for patTriggerEvent

Definition at line 33 of file PatTriggerAnalyzer.cc.

Referenced by analyze().