CMS 3D CMS Logo

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

Public Member Functions

void analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
 
void beginJob () override
 
void endJob () override
 
 ExampleMuonAnalyzer (const edm::ParameterSet &pset)
 Constructor. More...
 
 ~ExampleMuonAnalyzer () override
 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
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

TH1F * h4MuInvMass
 
TH1I * hChamberMatched
 
TH1F * hEHcal
 
TH1F * hMuCaloCompatibility
 
TH1I * hMuIdAlgo
 
TH1F * hMuIso03CaloComb
 
TH1F * hMuIso03SumPt
 
TH1I * hMuonType
 
TH1F * hMuSegCompatibility
 
TH1I * hNMuons
 
TH1F * hPtRec
 
TH2F * hPtReso
 
TH1F * hPtSTATKDiff
 
edm::EDGetTokenT< pat::MuonCollectiontheMuonToken
 

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

Analyzer of the muon objects

Author
R. Bellan - CERN ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 21 of file MuonAnalyzer.cc.

Constructor & Destructor Documentation

◆ ExampleMuonAnalyzer()

ExampleMuonAnalyzer::ExampleMuonAnalyzer ( const edm::ParameterSet pset)

Constructor.

Definition at line 66 of file MuonAnalyzer.cc.

References TFileService::kSharedResource, and muonDTDigis_cfi::pset.

66  {
67  theMuonToken = consumes<pat::MuonCollection>(pset.getParameter<InputTag>("MuonCollection"));
68  usesResource(TFileService::kSharedResource);
69 }
static const std::string kSharedResource
Definition: TFileService.h:76
edm::EDGetTokenT< pat::MuonCollection > theMuonToken
Definition: MuonAnalyzer.cc:38

◆ ~ExampleMuonAnalyzer()

ExampleMuonAnalyzer::~ExampleMuonAnalyzer ( )
override

Destructor.

Definition at line 72 of file MuonAnalyzer.cc.

72 {}

Member Function Documentation

◆ analyze()

void ExampleMuonAnalyzer::analyze ( const edm::Event event,
const edm::EventSetup eventSetup 
)
overridevirtual

simple selection... Do not want to write here my super-secret Higgs analysis ;-)

Implements edm::one::EDAnalyzerBase.

Definition at line 102 of file MuonAnalyzer.cc.

References muon::All, muon::AllArbitrated, muon::AllStandAloneMuons, muon::AllTrackerMuons, change_name::diff, muon::GlobalMuonPromptTight, muon::isGoodMuon(), HLT_2024v11_cff::muon, DiMuonV_cfg::muons, muon::segmentCompatibility(), METSignificanceObjects_cfi::selectedMuons, muon::TM2DCompatibilityLoose, muon::TM2DCompatibilityTight, muon::TMLastStationLoose, muon::TMLastStationOptimizedLowPtLoose, muon::TMLastStationOptimizedLowPtTight, muon::TMLastStationTight, muon::TMOneStationLoose, muon::TMOneStationTight, and muon::TrackerMuonArbitrated.

102  {
103  // Get the Muon collection
105  event.getByToken(theMuonToken, muons);
106 
107  // How many muons in the event?
108  hNMuons->Fill(muons->size());
109 
111 
112  // Let's look inside the muon collection.
113  for (pat::MuonCollection::const_iterator muon = muons->begin(); muon != muons->end(); ++muon) {
114  // pT spectra of muons
115  hPtRec->Fill(muon->pt());
116 
117  // what is the resolution in pt? Easy! We have the association with generated information
118  // cout<<muon->pt()<<" "<<muon->genParticle()->pt()<<endl;
119  if (muon->genLepton() != nullptr) {
120  double reso = (muon->pt() - muon->genLepton()->pt()) / muon->genLepton()->pt();
121  hPtReso->Fill(muon->genLepton()->pt(), reso);
122  }
123 
124  // What is the energy deposit in HCal?
125  if (muon->isEnergyValid())
126  hEHcal->Fill(muon->calEnergy().had);
127 
128  // Which type of muons in the collection?
129  if (muon->isStandAloneMuon())
130  if (muon->isGlobalMuon())
131  if (muon->isTrackerMuon())
132  hMuonType->Fill(1); // STA + GLB + TM
133  else
134  hMuonType->Fill(2); // STA + GLB
135  else if (muon->isTrackerMuon())
136  hMuonType->Fill(3); // STA + TM
137  else
138  hMuonType->Fill(5); // STA
139  else if (muon->isTrackerMuon())
140  hMuonType->Fill(4); // TM
141 
142  // ...mmm I want to study the relative resolution of the STA track with respect to the Tracker track.
143  // or I want to look at a track stab
144  if (muon->isGlobalMuon()) {
145  double diff = muon->innerTrack()->pt() - muon->standAloneMuon()->pt();
146  hPtSTATKDiff->Fill(diff);
147  }
148 
149  // Muon ID quantities
150 
151  // Muon in CMS are usually MIP. What is the compatibility of a muon HP using only claorimeters?
152  if (muon->isCaloCompatibilityValid())
153  hMuCaloCompatibility->Fill(muon->caloCompatibility());
154 
155  // The muon system can also be used just as only for ID. What is the compatibility of a muon HP using only calorimeters?
157 
158  // How many chambers have been associated to a muon track?
159  hChamberMatched->Fill(muon->numberOfChambers());
160  // If you look at MuonSegmentMatcher class you will see a lot of interesting quantities to look at!
161  // you can get the list of matched info using matches()
162 
163  // Muon ID selection. As described in AN-2008/098
164  if (muon::isGoodMuon(*muon, muon::All)) // dummy options - always true
165  hMuIdAlgo->Fill(0);
166  if (muon::isGoodMuon(*muon, muon::AllStandAloneMuons)) // checks isStandAloneMuon flag
167  hMuIdAlgo->Fill(1);
168  if (muon::isGoodMuon(*muon, muon::AllTrackerMuons)) // checks isTrackerMuon flag
169  hMuIdAlgo->Fill(2);
170  if (muon::isGoodMuon(*muon, muon::TrackerMuonArbitrated)) // resolve ambiguity of sharing segments
171  hMuIdAlgo->Fill(3);
172  if (muon::isGoodMuon(*muon, muon::AllArbitrated)) // all muons with the tracker muon arbitrated
173  hMuIdAlgo->Fill(4);
174  if (muon::isGoodMuon(*muon, muon::GlobalMuonPromptTight)) // global muons with tighter fit requirements
175  hMuIdAlgo->Fill(5);
176  if (muon::isGoodMuon(*muon, muon::TMLastStationLoose)) // penetration depth loose selector
177  hMuIdAlgo->Fill(6);
178  if (muon::isGoodMuon(*muon, muon::TMLastStationTight)) // penetration depth tight selector
179  hMuIdAlgo->Fill(7);
180  if (muon::isGoodMuon(*muon, muon::TM2DCompatibilityLoose)) // likelihood based loose selector
181  hMuIdAlgo->Fill(8);
182  if (muon::isGoodMuon(*muon, muon::TM2DCompatibilityTight)) // likelihood based tight selector
183  hMuIdAlgo->Fill(9);
184  if (muon::isGoodMuon(*muon, muon::TMOneStationLoose)) // require one well matched segment
185  hMuIdAlgo->Fill(10);
186  if (muon::isGoodMuon(*muon, muon::TMOneStationTight)) // require one well matched segment
187  hMuIdAlgo->Fill(11);
188  if (muon::isGoodMuon(*muon,
189  muon::TMLastStationOptimizedLowPtLoose)) // combination of TMLastStation and TMOneStation
190  hMuIdAlgo->Fill(12);
191  if (muon::isGoodMuon(*muon,
192  muon::TMLastStationOptimizedLowPtTight)) // combination of TMLastStation and TMOneStation
193  hMuIdAlgo->Fill(13);
194 
195  // Isolation variables. There are many type of isolation. You can even build your own combining the output of
196  // muon->isolationR03(). E.g.: 1.2*muon->isolationR03().emEt + 0.8*muon->isolationR03().hadEt
197  // *** WARNING *** it is just an EXAMPLE!
198  if (muon->isIsolationValid()) {
199  hMuIso03CaloComb->Fill(1.2 * muon->isolationR03().emEt + 0.8 * muon->isolationR03().hadEt);
200  hMuIso03SumPt->Fill(muon->isolationR03().sumPt);
201  }
202 
203  // OK, let see if we understood everything.
204  // Suppose we are searching for H->ZZ->4mu.
205  // In mean the 4 muons have/are:
206  // high pt (but 1 out of 4 can be at quite low pt)
207  // isolated
208  // so, we can make some requirements
209  if (muon->isolationR03().sumPt < 0.2) {
210  if (muon->isGlobalMuon() || muon::isGoodMuon(*muon, muon::TM2DCompatibilityLoose) ||
212  selectedMuons.push_back(*muon);
213  }
214  }
215 
217  if (selectedMuons.size() == 4) {
219  for (pat::MuonCollection::const_iterator muon = selectedMuons.begin(); muon != selectedMuons.end(); ++muon) {
220  p4CM = p4CM + muon->p4();
221  }
222  h4MuInvMass->Fill(p4CM.mass());
223  }
224 }
selectedMuons
====================================== Muons
muons
the two sets of parameters below are mutually exclusive, depending if RECO or ALCARECO is used the us...
Definition: DiMuonV_cfg.py:214
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
math::XYZTLorentzVector LorentzVector
Lorentz vector.
Definition: Candidate.h:36
std::vector< Muon > MuonCollection
Definition: Muon.h:35
edm::EDGetTokenT< pat::MuonCollection > theMuonToken
Definition: MuonAnalyzer.cc:38

◆ beginJob()

void ExampleMuonAnalyzer::beginJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 74 of file MuonAnalyzer.cc.

References TFileService::make().

74  {
75  // Book histograms
76  edm::Service<TFileService> fileService;
77  hPtRec = fileService->make<TH1F>("pT", "p_{T}^{rec}", 250, 0, 120);
78  hPtReso = fileService->make<TH2F>("pT_Reso", "(p_{T}^{rec}-p_{T}^{sim})/p_{T}^{sim}", 250, 0, 120, 100, -0.2, 0.2);
79  hNMuons = fileService->make<TH1I>("NMuons", "Number of muons per event", 20, 0, 20);
80 
81  hEHcal = fileService->make<TH1F>("EHCal", "Energy deposit in HCAL", 100, 0, 10);
82 
83  // ID
84  hMuonType = fileService->make<TH1I>("MuonType", "Type of Muons", 5, 1, 6);
85  hPtSTATKDiff = fileService->make<TH1F>("DiffPt_STA_TK", "p^{TK}_{T}-p^{STA}_T", 200, -50, 50);
86  hMuCaloCompatibility = fileService->make<TH1F>("CaloCompatibility", "Muon HP using Calorimeters only", 100, 0, 1);
87  hMuSegCompatibility = fileService->make<TH1F>("SegmentCompatibility", "Muon HP using segments only", 100, 0, 1);
88  hChamberMatched = fileService->make<TH1I>("NumMatchedChamber", "Number of matched chambers", 7, 0, 7);
89  hMuIdAlgo = fileService->make<TH1I>("MuonIDSelectors", "Results of muon id selectors", 13, 0, 13);
90 
91  // Isolation
92  hMuIso03SumPt = fileService->make<TH1F>("MuIso03SumPt", "Isolation #Delta(R)=0.3: SumPt", 200, 0, 10);
94  fileService->make<TH1F>("MuIso03CaloComb", "Isolation #Delta(R)=0.3: 1.2*ECAL+0.8HCAL", 200, 0, 10);
95 
96  // 4Mu invariant mass
97  h4MuInvMass = fileService->make<TH1F>("InvMass4MuSystem", "Invariant mass of the 4 muons system", 200, 0, 500);
98 }
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64

◆ endJob()

void ExampleMuonAnalyzer::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 100 of file MuonAnalyzer.cc.

100 {}

Member Data Documentation

◆ h4MuInvMass

TH1F* ExampleMuonAnalyzer::h4MuInvMass
private

Definition at line 59 of file MuonAnalyzer.cc.

◆ hChamberMatched

TH1I* ExampleMuonAnalyzer::hChamberMatched
private

Definition at line 52 of file MuonAnalyzer.cc.

◆ hEHcal

TH1F* ExampleMuonAnalyzer::hEHcal
private

Definition at line 44 of file MuonAnalyzer.cc.

◆ hMuCaloCompatibility

TH1F* ExampleMuonAnalyzer::hMuCaloCompatibility
private

Definition at line 50 of file MuonAnalyzer.cc.

◆ hMuIdAlgo

TH1I* ExampleMuonAnalyzer::hMuIdAlgo
private

Definition at line 53 of file MuonAnalyzer.cc.

◆ hMuIso03CaloComb

TH1F* ExampleMuonAnalyzer::hMuIso03CaloComb
private

Definition at line 57 of file MuonAnalyzer.cc.

◆ hMuIso03SumPt

TH1F* ExampleMuonAnalyzer::hMuIso03SumPt
private

Definition at line 56 of file MuonAnalyzer.cc.

◆ hMuonType

TH1I* ExampleMuonAnalyzer::hMuonType
private

Definition at line 46 of file MuonAnalyzer.cc.

◆ hMuSegCompatibility

TH1F* ExampleMuonAnalyzer::hMuSegCompatibility
private

Definition at line 51 of file MuonAnalyzer.cc.

◆ hNMuons

TH1I* ExampleMuonAnalyzer::hNMuons
private

Definition at line 41 of file MuonAnalyzer.cc.

◆ hPtRec

TH1F* ExampleMuonAnalyzer::hPtRec
private

Definition at line 42 of file MuonAnalyzer.cc.

◆ hPtReso

TH2F* ExampleMuonAnalyzer::hPtReso
private

Definition at line 43 of file MuonAnalyzer.cc.

◆ hPtSTATKDiff

TH1F* ExampleMuonAnalyzer::hPtSTATKDiff
private

Definition at line 47 of file MuonAnalyzer.cc.

◆ theMuonToken

edm::EDGetTokenT<pat::MuonCollection> ExampleMuonAnalyzer::theMuonToken
private

Definition at line 38 of file MuonAnalyzer.cc.