CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
PatTriggerAnalyzer Class Reference
Inheritance diagram for PatTriggerAnalyzer:
edm::one::EDAnalyzer< edm::one::SharedResources > edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

 PatTriggerAnalyzer (const edm::ParameterSet &iConfig)
 default constructor More...
 
 ~PatTriggerAnalyzer () override
 default destructor More...
 
- Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
const EDAnalyzeroperator= (const EDAnalyzer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () 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
 
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 const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
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::MuonCollection
muonsToken_
 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::TriggerEvent
triggerEventToken_
 input for patTriggerEvent More...
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< Bconsumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename 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<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

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.

References TFileService::kSharedResource.

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_() {
75  usesResource(TFileService::kSharedResource);
76 }
static const std::string kSharedResource
Definition: TFileService.h:76
std::map< std::string, TH2D * > histos2D_
std::map< std::string, TH1D * > histos1D_
histogram management
unsigned minID_
minimal id for meanPt plot
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
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 78 of file PatTriggerAnalyzer.cc.

78 {}

Member Function Documentation

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

everythin that needs to be done during the event loop

Implements edm::one::EDAnalyzerBase.

Definition at line 115 of file PatTriggerAnalyzer.cc.

References edm::Event::getByToken(), histos1D_, histos2D_, gpuClustering::id, maxID_, minID_, muonMatch_, patZpeak::muons, muonsToken_, sumN_, sumPt_, triggerEventToken_, and trigger::TriggerMuon.

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

Definition at line 80 of file PatTriggerAnalyzer.cc.

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

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

Definition at line 184 of file PatTriggerAnalyzer.cc.

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

184  {
185  // normalize the entries for the mean pt plot
186  for (unsigned id = minID_; id <= maxID_; ++id) {
187  if (sumN_[id] != 0)
188  histos1D_["ptMean"]->Fill(id, sumPt_[id] / sumN_[id]);
189  }
190 }
std::map< unsigned, unsigned > sumN_
internals for meanPt histogram calculation
uint16_t *__restrict__ id
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().