CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Types | Private Member Functions | Private Attributes
pat::PATMuonProducer Class Reference

class definition More...

#include "PhysicsTools/PatAlgos/interface/PATMuonProducer.h"

Inheritance diagram for pat::PATMuonProducer:
edm::stream::EDProducer< edm::GlobalCache< PATMuonHeavyObjectCache > >

Public Member Functions

 PATMuonProducer (const edm::ParameterSet &iConfig, PATMuonHeavyObjectCache const *)
 default constructir More...
 
void produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
 everything that needs to be done during the event loop More...
 
 ~PATMuonProducer () override
 default destructur More...
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< PATMuonHeavyObjectCache > >
 EDProducer ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 description of config file parameters More...
 
static void globalEndJob (PATMuonHeavyObjectCache *)
 
static std::unique_ptr< PATMuonHeavyObjectCache > initializeGlobalCache (const edm::ParameterSet &iConfig)
 

Private Types

typedef std::vector< edm::Handle< edm::Association< reco::GenParticleCollection > > > GenAssociations
 
typedef std::vector< edm::Handle< edm::ValueMap< IsoDeposit > > > IsoDepositMaps
 
typedef std::pair< pat::IsolationKeys, edm::InputTag > IsolationLabel
 
typedef std::vector< IsolationLabel > IsolationLabels
 
typedef std::vector< edm::Handle< edm::ValueMap< double > > > IsolationValueMaps
 
typedef edm::RefToBase< reco::Muon > MuonBaseRef
 typedefs for convenience More...
 

Private Member Functions

void embedHighLevel (pat::Muon &aMuon, reco::TrackRef track, reco::TransientTrack &tt, reco::Vertex &primaryVertex, bool primaryVertexIsValid, reco::BeamSpot &beamspot, bool beamspotIsValid)
 
void fillHltTriggerInfo (pat::Muon &muon, edm::Handle< std::vector< pat::TriggerObjectStandAlone >> &triggerObjects, const edm::TriggerNames &names, const std::vector< std::string > &collection_names)
 
void fillL1TriggerInfo (pat::Muon &muon, edm::Handle< std::vector< pat::TriggerObjectStandAlone >> &triggerObjects, const edm::TriggerNames &names, const edm::ESHandle< GlobalTrackingGeometry > &geometry)
 
void fillMuon (Muon &aMuon, const MuonBaseRef &muonRef, const reco::CandidateBaseRef &baseRef, const GenAssociations &genMatches, const IsoDepositMaps &deposits, const IsolationValueMaps &isolationValues) const
 common muon filling, for both the standard and PF2PAT case More...
 
std::optional< GlobalPoint > getMuonDirection (const reco::MuonChamberMatch &chamberMatch, const edm::ESHandle< GlobalTrackingGeometry > &geometry, const DetId &chamberId)
 
double getRelMiniIsoPUCorrected (const pat::Muon &muon, double rho, const std::vector< double > &area)
 
bool isChargedHadron (long pdgid)
 
bool isNeutralHadron (long pdgid)
 
bool isPhoton (long pdgid)
 
double puppiCombinedIsolation (const pat::Muon &muon, const pat::PackedCandidateCollection *pc)
 
template<typename T >
void readIsolationLabels (const edm::ParameterSet &iConfig, const char *psetName, IsolationLabels &labels, std::vector< edm::EDGetTokenT< edm::ValueMap< T >>> &tokens)
 
double relMiniIsoPUCorrected (const pat::Muon &aMuon, double rho)
 
void setMuonMiniIso (pat::Muon &aMuon, const pat::PackedCandidateCollection *pc)
 

Private Attributes

bool addEfficiencies_
 add efficiencies to the muon (this will be data members of th muon even w/o embedding) More...
 
bool addGenMatch_
 add generator match information More...
 
bool addInverseBeta_
 add combined inverse beta measurement into the muon More...
 
bool addPuppiIsolation_
 add puppi isolation More...
 
bool addResolutions_
 add resolutions to the muon (this will be data members of th muon even w/o embedding) More...
 
bool addTriggerMatching_
 Trigger. More...
 
edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
 input source of the primary vertex/beamspot More...
 
edm::EDGetTokenT< edm::ValueMap< reco::MuonMETCorrectionData > > caloMETMuonCorrsToken_
 source of caloMET muon corrections More...
 
bool computeMiniIso_
 
bool computeMuonMVA_
 standard muon selectors More...
 
bool computePuppiCombinedIso_
 
bool computeSoftMuonMVA_
 
std::vector< double > effectiveAreaVec_
 
pat::helper::EfficiencyLoader efficiencyLoader_
 helper class to add efficiencies to the muon More...
 
bool embedBestTrack_
 embed the track from best muon measurement (global pflow) More...
 
bool embedCaloMETMuonCorrs_
 embed muon MET correction info for caloMET into the muon More...
 
bool embedCombinedMuon_
 embed track of the combined fit into the muon More...
 
bool embedDytMuon_
 embed track from DYT muon fit into the muon More...
 
bool embedGenMatch_
 embed the gen match information into the muon More...
 
bool embedHighLevelSelection_
 embed high level selection variables More...
 
bool embedPFCandidate_
 embed pfCandidates into the muon More...
 
bool embedPfEcalEnergy_
 add ecal PF energy More...
 
bool embedPickyMuon_
 embed track from picky muon fit into the muon More...
 
bool embedStandAloneMuon_
 embed track from muon system into the muon More...
 
bool embedTcMETMuonCorrs_
 embed muon MET correction info for tcMET into the muon More...
 
bool embedTpfmsMuon_
 embed track from tpfms muon fit into the muon More...
 
bool embedTrack_
 embed the track from inner tracker into the muon More...
 
bool embedTunePBestTrack_
 embed the track from best muon measurement (muon only) More...
 
bool forceEmbedBestTrack_
 force separate embed of the best track even if already embedded More...
 
std::vector< edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > > > genMatchTokens_
 input tags for generator match information More...
 
bool getdBFromTrack_
 switch on reading the dB information from the track More...
 
std::vector< std::string > hltCollectionFilters_
 
IsolationLabels isoDepositLabels_
 input source for isoDeposits More...
 
std::vector< edm::EDGetTokenT< edm::ValueMap< IsoDeposit > > > isoDepositTokens_
 
IsolationLabels isolationValueLabels_
 input source isolation value maps More...
 
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueTokens_
 
pat::helper::MultiIsolator isolator_
 helper class to add userdefined isolation values to the muon More...
 
pat::helper::MultiIsolator::IsolationValuePairs isolatorTmpStorage_
 isolation value pair for temporary storage before being folded into the muon More...
 
std::vector< double > miniIsoParams_
 
edm::EDGetTokenT< edm::ValueMap< reco::MuonTimeExtra > > muonTimeExtraToken_
 input tag for reading inverse beta More...
 
edm::EDGetTokenT< edm::View< reco::Muon > > muonToken_
 input source More...
 
edm::EDGetTokenT< reco::JetTagCollection > mvaBTagCollectionTag_
 
edm::EDGetTokenT< reco::JetCorrector > mvaL1Corrector_
 
edm::EDGetTokenT< reco::JetCorrector > mvaL1L2L3ResCorrector_
 
bool mvaUseJec_
 
edm::EDGetTokenT< pat::PackedCandidateCollection > pcToken_
 
edm::EDGetTokenT< reco::PFCandidateCollection > pfMuonToken_
 input source pfCandidates that will be to be transformed into pat::Muons, when using PF2PAT More...
 
GreaterByPt< Muon > pTComparator_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_charged_hadrons_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_neutral_hadrons_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_photons_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_charged_hadrons_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_neutral_hadrons_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_photons_
 
edm::EDGetTokenT< std::vector< reco::Vertex > > pvToken_
 input source of the primary vertex More...
 
bool recomputeBasicSelectors_
 
double relMiniIsoPUCorrected_
 
pat::helper::KinResolutionsLoader resolutionLoader_
 helper class to add resolutions to the muon More...
 
edm::EDGetTokenT< double > rho_
 
edm::EDGetTokenT< edm::ValueMap< reco::MuonSimInfo > > simInfo_
 MC info. More...
 
edm::EDGetTokenT< edm::ValueMap< reco::MuonMETCorrectionData > > tcMETMuonCorrsToken_
 source of tcMET muon corrections More...
 
edm::EDGetTokenT< std::vector< pat::TriggerObjectStandAlone > > triggerObjects_
 
edm::EDGetTokenT< edm::TriggerResults > triggerResults_
 
bool useParticleFlow_
 switch to use particle flow (PF2PAT) or not More...
 
pat::PATUserDataHelper< pat::Muon > userDataHelper_
 helper class to add userData to the muon More...
 
bool useUserData_
 add user data to the muon (this will be data members of th muon even w/o embedding) More...
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< PATMuonHeavyObjectCache > >
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache > LuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCache > RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 

Detailed Description

class definition

Produces pat::Muon's.

The PATMuonProducer produces analysis-level pat::Muon's starting from a collection of objects of reco::Muon.

Author
Steven Lowette, Roger Wolf
Version
Id
PATMuonProducer.h,v 1.29 2012/08/22 15:02:52 bellan Exp

Definition at line 66 of file PATMuonProducer.h.

Member Typedef Documentation

typedef std::vector<edm::Handle<edm::Association<reco::GenParticleCollection> > > pat::PATMuonProducer::GenAssociations
private

Definition at line 87 of file PATMuonProducer.h.

typedef std::vector<edm::Handle<edm::ValueMap<IsoDeposit> > > pat::PATMuonProducer::IsoDepositMaps
private

Definition at line 88 of file PATMuonProducer.h.

typedef std::pair<pat::IsolationKeys, edm::InputTag> pat::PATMuonProducer::IsolationLabel
private

Definition at line 90 of file PATMuonProducer.h.

typedef std::vector<IsolationLabel> pat::PATMuonProducer::IsolationLabels
private

Definition at line 91 of file PATMuonProducer.h.

typedef std::vector<edm::Handle<edm::ValueMap<double> > > pat::PATMuonProducer::IsolationValueMaps
private

Definition at line 89 of file PATMuonProducer.h.

typedef edm::RefToBase<reco::Muon> pat::PATMuonProducer::MuonBaseRef
private

typedefs for convenience

Definition at line 86 of file PATMuonProducer.h.

Constructor & Destructor Documentation

PATMuonProducer::PATMuonProducer ( const edm::ParameterSet &  iConfig,
PATMuonHeavyObjectCache const *   
)
explicit

default constructir

Definition at line 75 of file PATMuonProducer.cc.

References addEfficiencies_, addGenMatch_, addInverseBeta_, addPuppiIsolation_, addResolutions_, addTriggerMatching_, beamLineToken_, caloMETMuonCorrsToken_, computeMiniIso_, computeMuonMVA_, computePuppiCombinedIso_, computeSoftMuonMVA_, effectiveAreaVec_, efficiencyLoader_, embedBestTrack_, embedCaloMETMuonCorrs_, embedCombinedMuon_, embedDytMuon_, embedGenMatch_, embedHighLevelSelection_, embedPFCandidate_, embedPfEcalEnergy_, embedPickyMuon_, embedStandAloneMuon_, embedTcMETMuonCorrs_, embedTpfmsMuon_, embedTrack_, embedTunePBestTrack_, Exception, edm::ParameterSet::existsAs(), forceEmbedBestTrack_, genMatchTokens_, getdBFromTrack_, edm::ParameterSet::getParameter(), hltCollectionFilters_, isoDepositLabels_, isoDepositTokens_, isolationValueLabels_, isolationValueTokens_, miniIsoParams_, muonTimeExtraToken_, muonToken_, mvaBTagCollectionTag_, mvaL1Corrector_, mvaL1L2L3ResCorrector_, mvaUseJec_, pcToken_, pfMuonToken_, PUPPIIsolation_charged_hadrons_, PUPPIIsolation_neutral_hadrons_, PUPPIIsolation_photons_, PUPPINoLeptonsIsolation_charged_hadrons_, PUPPINoLeptonsIsolation_neutral_hadrons_, PUPPINoLeptonsIsolation_photons_, pvToken_, readIsolationLabels(), recomputeBasicSelectors_, resolutionLoader_, rho_, simInfo_, GlobalPosition_Frontier_DevDB_cff::tag, tcMETMuonCorrsToken_, triggerObjects_, triggerResults_, useParticleFlow_, userDataHelper_, useUserData_, and edm::vector_transform().

76  : relMiniIsoPUCorrected_(0),
77  useUserData_(iConfig.exists("userData")),
78  computeMuonMVA_(false),
79  computeSoftMuonMVA_(false),
80  recomputeBasicSelectors_(false),
81  mvaUseJec_(false),
82  isolator_(iConfig.exists("userIsolation") ? iConfig.getParameter<edm::ParameterSet>("userIsolation")
83  : edm::ParameterSet(),
84  consumesCollector(),
85  false) {
86  // input source
87  muonToken_ = consumes<edm::View<reco::Muon>>(iConfig.getParameter<edm::InputTag>("muonSource"));
88  // embedding of tracks
89  embedBestTrack_ = iConfig.getParameter<bool>("embedMuonBestTrack");
90  embedTunePBestTrack_ = iConfig.getParameter<bool>("embedTunePMuonBestTrack");
91  forceEmbedBestTrack_ = iConfig.getParameter<bool>("forceBestTrackEmbedding");
92  embedTrack_ = iConfig.getParameter<bool>("embedTrack");
93  embedCombinedMuon_ = iConfig.getParameter<bool>("embedCombinedMuon");
94  embedStandAloneMuon_ = iConfig.getParameter<bool>("embedStandAloneMuon");
95  // embedding of muon MET correction information
96  embedCaloMETMuonCorrs_ = iConfig.getParameter<bool>("embedCaloMETMuonCorrs");
97  embedTcMETMuonCorrs_ = iConfig.getParameter<bool>("embedTcMETMuonCorrs");
98  caloMETMuonCorrsToken_ =
99  mayConsume<edm::ValueMap<reco::MuonMETCorrectionData>>(iConfig.getParameter<edm::InputTag>("caloMETMuonCorrs"));
100  tcMETMuonCorrsToken_ =
101  mayConsume<edm::ValueMap<reco::MuonMETCorrectionData>>(iConfig.getParameter<edm::InputTag>("tcMETMuonCorrs"));
102  // pflow specific configurables
103  useParticleFlow_ = iConfig.getParameter<bool>("useParticleFlow");
104  embedPFCandidate_ = iConfig.getParameter<bool>("embedPFCandidate");
105  pfMuonToken_ = mayConsume<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfMuonSource"));
106  embedPfEcalEnergy_ = iConfig.getParameter<bool>("embedPfEcalEnergy");
107  // embedding of tracks from TeV refit
108  embedPickyMuon_ = iConfig.getParameter<bool>("embedPickyMuon");
109  embedTpfmsMuon_ = iConfig.getParameter<bool>("embedTpfmsMuon");
110  embedDytMuon_ = iConfig.getParameter<bool>("embedDytMuon");
111  // embedding of inverse beta variable information
112  addInverseBeta_ = iConfig.getParameter<bool>("addInverseBeta");
113  if (addInverseBeta_) {
114  muonTimeExtraToken_ =
115  consumes<edm::ValueMap<reco::MuonTimeExtra>>(iConfig.getParameter<edm::InputTag>("sourceMuonTimeExtra"));
116  }
117  // Switch to get the dB from the track instead of using IPTools
118  getdBFromTrack_ = iConfig.getParameter<bool>("getdBFromTrack");
119  // Monte Carlo matching
120  addGenMatch_ = iConfig.getParameter<bool>("addGenMatch");
121  if (addGenMatch_) {
122  embedGenMatch_ = iConfig.getParameter<bool>("embedGenMatch");
123  if (iConfig.existsAs<edm::InputTag>("genParticleMatch")) {
124  genMatchTokens_.push_back(consumes<edm::Association<reco::GenParticleCollection>>(
125  iConfig.getParameter<edm::InputTag>("genParticleMatch")));
126  } else {
127  genMatchTokens_ = edm::vector_transform(
128  iConfig.getParameter<std::vector<edm::InputTag>>("genParticleMatch"),
129  [this](edm::InputTag const& tag) { return consumes<edm::Association<reco::GenParticleCollection>>(tag); });
130  }
131  }
132  // efficiencies
133  addEfficiencies_ = iConfig.getParameter<bool>("addEfficiencies");
134  if (addEfficiencies_) {
135  efficiencyLoader_ =
136  pat::helper::EfficiencyLoader(iConfig.getParameter<edm::ParameterSet>("efficiencies"), consumesCollector());
137  }
138  // resolutions
139  addResolutions_ = iConfig.getParameter<bool>("addResolutions");
140  if (addResolutions_) {
141  resolutionLoader_ = pat::helper::KinResolutionsLoader(iConfig.getParameter<edm::ParameterSet>("resolutions"));
142  }
143  // puppi
144  addPuppiIsolation_ = iConfig.getParameter<bool>("addPuppiIsolation");
145  if (addPuppiIsolation_) {
146  PUPPIIsolation_charged_hadrons_ =
147  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationChargedHadrons"));
148  PUPPIIsolation_neutral_hadrons_ =
149  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationNeutralHadrons"));
150  PUPPIIsolation_photons_ =
151  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationPhotons"));
152  //puppiNoLeptons
153  PUPPINoLeptonsIsolation_charged_hadrons_ =
154  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiNoLeptonsIsolationChargedHadrons"));
155  PUPPINoLeptonsIsolation_neutral_hadrons_ =
156  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiNoLeptonsIsolationNeutralHadrons"));
157  PUPPINoLeptonsIsolation_photons_ =
158  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiNoLeptonsIsolationPhotons"));
159  }
160  // read isoDeposit labels, for direct embedding
161  readIsolationLabels(iConfig, "isoDeposits", isoDepositLabels_, isoDepositTokens_);
162  // read isolation value labels, for direct embedding
163  readIsolationLabels(iConfig, "isolationValues", isolationValueLabels_, isolationValueTokens_);
164  // check to see if the user wants to add user data
165  if (useUserData_) {
166  userDataHelper_ = PATUserDataHelper<Muon>(iConfig.getParameter<edm::ParameterSet>("userData"), consumesCollector());
167  }
168  // embed high level selection variables
169  embedHighLevelSelection_ = iConfig.getParameter<bool>("embedHighLevelSelection");
170  if (embedHighLevelSelection_) {
171  beamLineToken_ = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamLineSrc"));
172  pvToken_ = consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("pvSrc"));
173  }
174 
175  //for mini-isolation calculation
176  computeMiniIso_ = iConfig.getParameter<bool>("computeMiniIso");
177 
178  computePuppiCombinedIso_ = iConfig.getParameter<bool>("computePuppiCombinedIso");
179 
180  effectiveAreaVec_ = iConfig.getParameter<std::vector<double>>("effectiveAreaVec");
181 
182  miniIsoParams_ = iConfig.getParameter<std::vector<double>>("miniIsoParams");
183  if (computeMiniIso_ && miniIsoParams_.size() != 9) {
184  throw cms::Exception("ParameterError") << "miniIsoParams must have exactly 9 elements.\n";
185  }
186  if (computeMiniIso_ || computePuppiCombinedIso_)
187  pcToken_ = consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfCandsForMiniIso"));
188 
189  // standard selectors
190  recomputeBasicSelectors_ = iConfig.getParameter<bool>("recomputeBasicSelectors");
191  computeMuonMVA_ = iConfig.getParameter<bool>("computeMuonMVA");
192  if (computeMuonMVA_ and not computeMiniIso_)
193  throw cms::Exception("ConfigurationError") << "MiniIso is needed for Muon MVA calculation.\n";
194 
195  if (computeMuonMVA_) {
196  // pfCombinedInclusiveSecondaryVertexV2BJetTags
197  mvaBTagCollectionTag_ = consumes<reco::JetTagCollection>(iConfig.getParameter<edm::InputTag>("mvaJetTag"));
198  mvaL1Corrector_ = consumes<reco::JetCorrector>(iConfig.getParameter<edm::InputTag>("mvaL1Corrector"));
199  mvaL1L2L3ResCorrector_ = consumes<reco::JetCorrector>(iConfig.getParameter<edm::InputTag>("mvaL1L2L3ResCorrector"));
200  rho_ = consumes<double>(iConfig.getParameter<edm::InputTag>("rho"));
201  mvaUseJec_ = iConfig.getParameter<bool>("mvaUseJec");
202  }
203 
204  computeSoftMuonMVA_ = iConfig.getParameter<bool>("computeSoftMuonMVA");
205 
206  // MC info
207  simInfo_ = consumes<edm::ValueMap<reco::MuonSimInfo>>(iConfig.getParameter<edm::InputTag>("muonSimInfo"));
208 
209  addTriggerMatching_ = iConfig.getParameter<bool>("addTriggerMatching");
210  if (addTriggerMatching_) {
211  triggerObjects_ =
212  consumes<std::vector<pat::TriggerObjectStandAlone>>(iConfig.getParameter<edm::InputTag>("triggerObjects"));
213  triggerResults_ = consumes<edm::TriggerResults>(iConfig.getParameter<edm::InputTag>("triggerResults"));
214  }
215  hltCollectionFilters_ = iConfig.getParameter<std::vector<std::string>>("hltCollectionFilters");
216 
217  // produces vector of muons
218  produces<std::vector<Muon>>();
219 }
pat::PATMuonProducer::embedTpfmsMuon_
bool embedTpfmsMuon_
embed track from tpfms muon fit into the muon
Definition: PATMuonProducer.h:173
pat::PATMuonProducer::PUPPIIsolation_charged_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_charged_hadrons_
Definition: PATMuonProducer.h:219
pat::PATMuonProducer::useUserData_
bool useUserData_
add user data to the muon (this will be data members of th muon even w/o embedding) ...
Definition: PATMuonProducer.h:213
pat::PATMuonProducer::readIsolationLabels
void readIsolationLabels(const edm::ParameterSet &iConfig, const char *psetName, IsolationLabels &labels, std::vector< edm::EDGetTokenT< edm::ValueMap< T >>> &tokens)
Definition: PATMuonProducer.h:261
pat::PATMuonProducer::PUPPIIsolation_photons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_photons_
Definition: PATMuonProducer.h:221
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
pat::PATUserDataHelper
Assists in assimilating all pat::UserData into pat objects.
Definition: PATUserDataHelper.h:53
pat::PATMuonProducer::addPuppiIsolation_
bool addPuppiIsolation_
add puppi isolation
Definition: PATMuonProducer.h:217
GlobalPosition_Frontier_DevDB_cff.tag
tag
Definition: GlobalPosition_Frontier_DevDB_cff.py:11
pat::PATMuonProducer::tcMETMuonCorrsToken_
edm::EDGetTokenT< edm::ValueMap< reco::MuonMETCorrectionData > > tcMETMuonCorrsToken_
source of tcMET muon corrections
Definition: PATMuonProducer.h:169
pat::PATMuonProducer::embedTcMETMuonCorrs_
bool embedTcMETMuonCorrs_
embed muon MET correction info for tcMET into the muon
Definition: PATMuonProducer.h:167
edm::ParameterSet::existsAs
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:161
pat::PATMuonProducer::pvToken_
edm::EDGetTokenT< std::vector< reco::Vertex > > pvToken_
input source of the primary vertex
Definition: PATMuonProducer.h:203
pat::PATMuonProducer::muonTimeExtraToken_
edm::EDGetTokenT< edm::ValueMap< reco::MuonTimeExtra > > muonTimeExtraToken_
input tag for reading inverse beta
Definition: PATMuonProducer.h:181
pat::PATMuonProducer::PUPPINoLeptonsIsolation_photons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_photons_
Definition: PATMuonProducer.h:225
pat::PATMuonProducer::rho_
edm::EDGetTokenT< double > rho_
Definition: PATMuonProducer.h:234
pat::PATMuonProducer::muonToken_
edm::EDGetTokenT< edm::View< reco::Muon > > muonToken_
input source
Definition: PATMuonProducer.h:140
pat::PATMuonProducer::addEfficiencies_
bool addEfficiencies_
add efficiencies to the muon (this will be data members of th muon even w/o embedding) ...
Definition: PATMuonProducer.h:211
pat::PATMuonProducer::caloMETMuonCorrsToken_
edm::EDGetTokenT< edm::ValueMap< reco::MuonMETCorrectionData > > caloMETMuonCorrsToken_
source of caloMET muon corrections
Definition: PATMuonProducer.h:165
pat::PATMuonProducer::embedCaloMETMuonCorrs_
bool embedCaloMETMuonCorrs_
embed muon MET correction info for caloMET into the muon
Definition: PATMuonProducer.h:163
pat::PATMuonProducer::simInfo_
edm::EDGetTokenT< edm::ValueMap< reco::MuonSimInfo > > simInfo_
MC info.
Definition: PATMuonProducer.h:249
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:674
pat::PATMuonProducer::addTriggerMatching_
bool addTriggerMatching_
Trigger.
Definition: PATMuonProducer.h:252
pat::PATMuonProducer::embedBestTrack_
bool embedBestTrack_
embed the track from best muon measurement (global pflow)
Definition: PATMuonProducer.h:151
pat::PATMuonProducer::mvaL1Corrector_
edm::EDGetTokenT< reco::JetCorrector > mvaL1Corrector_
Definition: PATMuonProducer.h:232
pat::PATMuonProducer::addResolutions_
bool addResolutions_
add resolutions to the muon (this will be data members of th muon even w/o embedding) ...
Definition: PATMuonProducer.h:189
pat::PATMuonProducer::PUPPINoLeptonsIsolation_charged_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_charged_hadrons_
Definition: PATMuonProducer.h:223
pat::PATMuonProducer::miniIsoParams_
std::vector< double > miniIsoParams_
Definition: PATMuonProducer.h:147
pat::PATMuonProducer::effectiveAreaVec_
std::vector< double > effectiveAreaVec_
Definition: PATMuonProducer.h:146
pat::PATMuonProducer::PUPPINoLeptonsIsolation_neutral_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_neutral_hadrons_
Definition: PATMuonProducer.h:224
pat::PATMuonProducer::mvaBTagCollectionTag_
edm::EDGetTokenT< reco::JetTagCollection > mvaBTagCollectionTag_
Definition: PATMuonProducer.h:231
pat::PATMuonProducer::pfMuonToken_
edm::EDGetTokenT< reco::PFCandidateCollection > pfMuonToken_
input source pfCandidates that will be to be transformed into pat::Muons, when using PF2PAT ...
Definition: PATMuonProducer.h:195
pat::PATMuonProducer::isolator_
pat::helper::MultiIsolator isolator_
helper class to add userdefined isolation values to the muon
Definition: PATMuonProducer.h:240
pat::PATMuonProducer::isoDepositTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< IsoDeposit > > > isoDepositTokens_
Definition: PATMuonProducer.h:206
pat::PATMuonProducer::beamLineToken_
edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
input source of the primary vertex/beamspot
Definition: PATMuonProducer.h:201
pat::PATMuonProducer::addInverseBeta_
bool addInverseBeta_
add combined inverse beta measurement into the muon
Definition: PATMuonProducer.h:177
pat::PATMuonProducer::useParticleFlow_
bool useParticleFlow_
switch to use particle flow (PF2PAT) or not
Definition: PATMuonProducer.h:193
pat::PATMuonProducer::userDataHelper_
pat::PATUserDataHelper< pat::Muon > userDataHelper_
helper class to add userData to the muon
Definition: PATMuonProducer.h:246
pat::PATMuonProducer::embedStandAloneMuon_
bool embedStandAloneMuon_
embed track from muon system into the muon
Definition: PATMuonProducer.h:159
pat::PATMuonProducer::triggerResults_
edm::EDGetTokenT< edm::TriggerResults > triggerResults_
Definition: PATMuonProducer.h:254
pat::PATMuonProducer::embedTrack_
bool embedTrack_
embed the track from inner tracker into the muon
Definition: PATMuonProducer.h:157
pat::PATMuonProducer::addGenMatch_
bool addGenMatch_
add generator match information
Definition: PATMuonProducer.h:183
pat::PATMuonProducer::embedPfEcalEnergy_
bool embedPfEcalEnergy_
add ecal PF energy
Definition: PATMuonProducer.h:215
pat::PATMuonProducer::getdBFromTrack_
bool getdBFromTrack_
switch on reading the dB information from the track
Definition: PATMuonProducer.h:179
pat::PATMuonProducer::pcToken_
edm::EDGetTokenT< pat::PackedCandidateCollection > pcToken_
Definition: PATMuonProducer.h:143
pat::PATMuonProducer::hltCollectionFilters_
std::vector< std::string > hltCollectionFilters_
Definition: PATMuonProducer.h:255
edm::vector_transform
auto vector_transform(std::vector< InputType > const &input, Function predicate) -> std::vector< typename std::remove_cv< typename std::remove_reference< decltype(predicate(input.front()))>::type >::type >
Definition: transform.h:11
pat::PATMuonProducer::embedTunePBestTrack_
bool embedTunePBestTrack_
embed the track from best muon measurement (muon only)
Definition: PATMuonProducer.h:153
pat::PATMuonProducer::mvaL1L2L3ResCorrector_
edm::EDGetTokenT< reco::JetCorrector > mvaL1L2L3ResCorrector_
Definition: PATMuonProducer.h:233
pat::PATMuonProducer::efficiencyLoader_
pat::helper::EfficiencyLoader efficiencyLoader_
helper class to add efficiencies to the muon
Definition: PATMuonProducer.h:244
pat::PATMuonProducer::embedPickyMuon_
bool embedPickyMuon_
embed track from picky muon fit into the muon
Definition: PATMuonProducer.h:171
pat::PATMuonProducer::embedDytMuon_
bool embedDytMuon_
embed track from DYT muon fit into the muon
Definition: PATMuonProducer.h:175
pat::PATMuonProducer::triggerObjects_
edm::EDGetTokenT< std::vector< pat::TriggerObjectStandAlone > > triggerObjects_
Definition: PATMuonProducer.h:253
pat::PATMuonProducer::embedGenMatch_
bool embedGenMatch_
embed the gen match information into the muon
Definition: PATMuonProducer.h:187
pat::PATMuonProducer::isoDepositLabels_
IsolationLabels isoDepositLabels_
input source for isoDeposits
Definition: PATMuonProducer.h:205
pat::PATMuonProducer::genMatchTokens_
std::vector< edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > > > genMatchTokens_
input tags for generator match information
Definition: PATMuonProducer.h:185
pat::PATMuonProducer::isolationValueLabels_
IsolationLabels isolationValueLabels_
input source isolation value maps
Definition: PATMuonProducer.h:208
pat::PATMuonProducer::forceEmbedBestTrack_
bool forceEmbedBestTrack_
force separate embed of the best track even if already embedded
Definition: PATMuonProducer.h:155
pat::PATMuonProducer::embedPFCandidate_
bool embedPFCandidate_
embed pfCandidates into the muon
Definition: PATMuonProducer.h:197
pat::PATMuonProducer::isolationValueTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueTokens_
Definition: PATMuonProducer.h:209
pat::PATMuonProducer::computeMuonMVA_
bool computeMuonMVA_
standard muon selectors
Definition: PATMuonProducer.h:227
pat::PATMuonProducer::embedCombinedMuon_
bool embedCombinedMuon_
embed track of the combined fit into the muon
Definition: PATMuonProducer.h:161
pat::PATMuonProducer::resolutionLoader_
pat::helper::KinResolutionsLoader resolutionLoader_
helper class to add resolutions to the muon
Definition: PATMuonProducer.h:191
pat::PATMuonProducer::PUPPIIsolation_neutral_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_neutral_hadrons_
Definition: PATMuonProducer.h:220
pat::PATMuonProducer::embedHighLevelSelection_
bool embedHighLevelSelection_
embed high level selection variables
Definition: PATMuonProducer.h:199
PATMuonProducer::~PATMuonProducer ( )
override

default destructur

Definition at line 221 of file PATMuonProducer.cc.

221 {}

Member Function Documentation

void PATMuonProducer::embedHighLevel ( pat::Muon &  aMuon,
reco::TrackRef  track,
reco::TransientTrack &  tt,
reco::Vertex &  primaryVertex,
bool  primaryVertexIsValid,
reco::BeamSpot &  beamspot,
bool  beamspotIsValid 
)
private

Definition at line 1095 of file PATMuonProducer.cc.

References pat::Muon::BS2D, pat::Muon::BS3D, reco::Vertex::covariance(), DEFINE_FWK_MODULE, getdBFromTrack_, reco::BeamSpot::position(), reco::Vertex::position(), pat::Muon::PV2D, pat::Muon::PV3D, pat::Muon::PVDZ, mps_fire::result, reco::BeamSpot::rotatedCovariance3D(), pat::Muon::setDB(), IPTools::signedImpactParameter3D(), IPTools::signedTransverseImpactParameter(), and reco::Vertex::zError().

Referenced by produce().

1101  {
1102  // Generic variable to store measurements
1103  std::pair<bool,Measurement1D> result;
1104  double d0_corr;
1105  double d0_err;
1106 
1107  // Correct to PV
1108 
1109  // PV2D
1110  if (getdBFromTrack_){
1111  aMuon.setDB(track->dxy(primaryVertex.position()),
1112  track->dxyError(primaryVertex.position(), primaryVertex.covariance()),
1113  pat::Muon::PV2D);
1114  } else{
1115  result =
1116  IPTools::signedTransverseImpactParameter(tt, GlobalVector(track->px(), track->py(), track->pz()), primaryVertex);
1117  d0_corr = result.second.value();
1118  d0_err = primaryVertexIsValid ? result.second.error() : -1.0;
1119  aMuon.setDB( d0_corr, d0_err, pat::Muon::PV2D);
1120  }
1121 
1122  // PV3D
1123  result = IPTools::signedImpactParameter3D(tt, GlobalVector(track->px(), track->py(), track->pz()), primaryVertex);
1124  d0_corr = result.second.value();
1125  d0_err = primaryVertexIsValid ? result.second.error() : -1.0;
1126  aMuon.setDB(d0_corr, d0_err, pat::Muon::PV3D);
1127 
1128  // Correct to beam spot
1129 
1130  // BS2D
1131  // make a fake vertex out of beam spot
1132  reco::Vertex vBeamspot(beamspot.position(), beamspot.rotatedCovariance3D());
1133  if (getdBFromTrack_){
1134  aMuon.setDB(track->dxy(beamspot), track->dxyError(beamspot), pat::Muon::BS2D);
1135  } else{
1136  result = IPTools::signedTransverseImpactParameter(tt, GlobalVector(track->px(), track->py(), track->pz()), vBeamspot);
1137  d0_corr = result.second.value();
1138  d0_err = beamspotIsValid ? result.second.error() : -1.0;
1139  aMuon.setDB( d0_corr, d0_err, pat::Muon::BS2D);
1140  }
1141 
1142  // BS3D
1143  result = IPTools::signedImpactParameter3D(tt, GlobalVector(track->px(), track->py(), track->pz()), vBeamspot);
1144  d0_corr = result.second.value();
1145  d0_err = beamspotIsValid ? result.second.error() : -1.0;
1146  aMuon.setDB(d0_corr, d0_err, pat::Muon::BS3D);
1147 
1148  // PVDZ
1149  aMuon.setDB(
1150  track->dz(primaryVertex.position()), std::hypot(track->dzError(), primaryVertex.zError()), pat::Muon::PVDZ);
1151 }
mps_fire.result
result
Definition: mps_fire.py:291
reco::Vertex::zError
double zError() const
error on z
Definition: Vertex.h:123
IPTools::signedTransverseImpactParameter
std::pair< bool, Measurement1D > signedTransverseImpactParameter(const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)
Definition: IPTools.cc:50
IPTools::signedImpactParameter3D
std::pair< bool, Measurement1D > signedImpactParameter3D(const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)
Definition: IPTools.cc:71
reco::Vertex::covariance
double covariance(int i, int j) const
(i, j)-th element of error matrix, i, j = 0, ... 2
Definition: Vertex.h:130
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:109
pat::Muon::setDB
void setDB(double dB, double edB, IPTYPE type=PV2D)
Definition: Muon.h:245
pat::PATMuonProducer::getdBFromTrack_
bool getdBFromTrack_
switch on reading the dB information from the track
Definition: PATMuonProducer.h:179
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
reco::BeamSpot::rotatedCovariance3D
Covariance3DMatrix rotatedCovariance3D() const
Definition: BeamSpot.cc:78
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10
void PATMuonProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

description of config file parameters

Definition at line 947 of file PATMuonProducer.cc.

References edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addNode(), edm::ParameterSetDescription::addOptional(), pat::helper::KinResolutionsLoader::fillDescription(), pat::PATUserDataHelper< ObjectType >::fillDescription(), edm::ParameterSetDescription::ifValue(), or, edm::ParameterSetDescription::setAllowAnything(), edm::ParameterSetDescription::setComment(), and edm::ParameterDescriptionNode::setComment().

947  {
948  edm::ParameterSetDescription iDesc;
949  iDesc.setComment("PAT muon producer module");
950 
951  // input source
952  iDesc.add<edm::InputTag>("muonSource", edm::InputTag("no default"))->setComment("input collection");
953 
954  // embedding
955  iDesc.add<bool>("embedMuonBestTrack", true)->setComment("embed muon best track (global pflow)");
956  iDesc.add<bool>("embedTunePMuonBestTrack", true)->setComment("embed muon best track (muon only)");
957  iDesc.add<bool>("forceBestTrackEmbedding", true)
958  ->setComment(
959  "force embedding separately the best tracks even if they're already embedded e.g. as tracker or global "
960  "tracks");
961  iDesc.add<bool>("embedTrack", true)->setComment("embed external track");
962  iDesc.add<bool>("embedStandAloneMuon", true)->setComment("embed external stand-alone muon");
963  iDesc.add<bool>("embedCombinedMuon", false)->setComment("embed external combined muon");
964  iDesc.add<bool>("embedPickyMuon", false)->setComment("embed external picky track");
965  iDesc.add<bool>("embedTpfmsMuon", false)->setComment("embed external tpfms track");
966  iDesc.add<bool>("embedDytMuon", false)->setComment("embed external dyt track ");
967 
968  // embedding of MET muon corrections
969  iDesc.add<bool>("embedCaloMETMuonCorrs", true)->setComment("whether to add MET muon correction for caloMET or not");
970  iDesc.add<edm::InputTag>("caloMETMuonCorrs", edm::InputTag("muonMETValueMapProducer", "muCorrData"))
971  ->setComment("source of MET muon corrections for caloMET");
972  iDesc.add<bool>("embedTcMETMuonCorrs", true)->setComment("whether to add MET muon correction for tcMET or not");
973  iDesc.add<edm::InputTag>("tcMETMuonCorrs", edm::InputTag("muonTCMETValueMapProducer", "muCorrData"))
974  ->setComment("source of MET muon corrections for tcMET");
975 
976  // pf specific parameters
977  iDesc.add<edm::InputTag>("pfMuonSource", edm::InputTag("pfMuons"))->setComment("particle flow input collection");
978  iDesc.add<bool>("useParticleFlow", false)->setComment("whether to use particle flow or not");
979  iDesc.add<bool>("embedPFCandidate", false)->setComment("embed external particle flow object");
980  iDesc.add<bool>("embedPfEcalEnergy", true)->setComment("add ecal energy as reconstructed by PF");
981 
982  // inverse beta computation
983  iDesc.add<bool>("addInverseBeta", true)->setComment("add combined inverse beta");
984  iDesc.add<edm::InputTag>("sourceInverseBeta", edm::InputTag("muons", "combined"))
985  ->setComment("source of inverse beta values");
986 
987  // switch to get the IP from the best track instead of running IPTools
988  iDesc.add<bool>("getdBFromTrack", false)->setComment("switch IP2D computation to use the best track one");
989 
990  // MC matching configurables
991  iDesc.add<bool>("addGenMatch", true)->setComment("add MC matching");
992  iDesc.add<bool>("embedGenMatch", false)->setComment("embed MC matched MC information");
993  std::vector<edm::InputTag> emptySourceVector;
994  iDesc
995  .addNode(edm::ParameterDescription<edm::InputTag>("genParticleMatch", edm::InputTag(), true) xor
996  edm::ParameterDescription<std::vector<edm::InputTag>>("genParticleMatch", emptySourceVector, true))
997  ->setComment("input with MC match information");
998 
999  // mini-iso
1000  iDesc.add<bool>("computeMiniIso", false)->setComment("whether or not to compute and store electron mini-isolation");
1001  iDesc.add<bool>("computePuppiCombinedIso", false)
1002  ->setComment("whether or not to compute and store puppi combined isolation");
1003 
1004  iDesc.add<edm::InputTag>("pfCandsForMiniIso", edm::InputTag("packedPFCandidates"))
1005  ->setComment("collection to use to compute mini-iso");
1006  iDesc.add<std::vector<double>>("miniIsoParams", std::vector<double>())
1007  ->setComment("mini-iso parameters to use for muons");
1008 
1009  iDesc.add<bool>("addTriggerMatching", false)->setComment("add L1 and HLT matching to offline muon");
1010 
1011  pat::helper::KinResolutionsLoader::fillDescription(iDesc);
1012 
1013  // IsoDeposit configurables
1014  edm::ParameterSetDescription isoDepositsPSet;
1015  isoDepositsPSet.addOptional<edm::InputTag>("tracker");
1016  isoDepositsPSet.addOptional<edm::InputTag>("ecal");
1017  isoDepositsPSet.addOptional<edm::InputTag>("hcal");
1018  isoDepositsPSet.addOptional<edm::InputTag>("particle");
1019  isoDepositsPSet.addOptional<edm::InputTag>("pfChargedHadrons");
1020  isoDepositsPSet.addOptional<edm::InputTag>("pfChargedAll");
1021  isoDepositsPSet.addOptional<edm::InputTag>("pfPUChargedHadrons");
1022  isoDepositsPSet.addOptional<edm::InputTag>("pfNeutralHadrons");
1023  isoDepositsPSet.addOptional<edm::InputTag>("pfPhotons");
1024  isoDepositsPSet.addOptional<std::vector<edm::InputTag>>("user");
1025  iDesc.addOptional("isoDeposits", isoDepositsPSet);
1026 
1027  // isolation values configurables
1028  edm::ParameterSetDescription isolationValuesPSet;
1029  isolationValuesPSet.addOptional<edm::InputTag>("tracker");
1030  isolationValuesPSet.addOptional<edm::InputTag>("ecal");
1031  isolationValuesPSet.addOptional<edm::InputTag>("hcal");
1032  isolationValuesPSet.addOptional<edm::InputTag>("particle");
1033  isolationValuesPSet.addOptional<edm::InputTag>("pfChargedHadrons");
1034  isolationValuesPSet.addOptional<edm::InputTag>("pfChargedAll");
1035  isolationValuesPSet.addOptional<edm::InputTag>("pfPUChargedHadrons");
1036  isolationValuesPSet.addOptional<edm::InputTag>("pfNeutralHadrons");
1037  isolationValuesPSet.addOptional<edm::InputTag>("pfPhotons");
1038  iDesc.addOptional("isolationValues", isolationValuesPSet);
1039 
1040  iDesc.ifValue(edm::ParameterDescription<bool>("addPuppiIsolation", false, true),
1041  true >> (edm::ParameterDescription<edm::InputTag>(
1042  "puppiIsolationChargedHadrons",
1043  edm::InputTag("muonPUPPIIsolation", "h+-DR030-ThresholdVeto000-ConeVeto000"),
1044  true) and
1045  edm::ParameterDescription<edm::InputTag>(
1046  "puppiIsolationNeutralHadrons",
1047  edm::InputTag("muonPUPPIIsolation", "h0-DR030-ThresholdVeto000-ConeVeto001"),
1048  true) and
1049  edm::ParameterDescription<edm::InputTag>(
1050  "puppiIsolationPhotons",
1051  edm::InputTag("muonPUPPIIsolation", "gamma-DR030-ThresholdVeto000-ConeVeto001"),
1052  true) and
1053  edm::ParameterDescription<edm::InputTag>(
1054  "puppiNoLeptonsIsolationChargedHadrons",
1055  edm::InputTag("muonPUPPINoLeptonsIsolation", "h+-DR030-ThresholdVeto000-ConeVeto000"),
1056  true) and
1057  edm::ParameterDescription<edm::InputTag>(
1058  "puppiNoLeptonsIsolationNeutralHadrons",
1059  edm::InputTag("muonPUPPINoLeptonsIsolation", "h0-DR030-ThresholdVeto000-ConeVeto001"),
1060  true) and
1061  edm::ParameterDescription<edm::InputTag>(
1062  "puppiNoLeptonsIsolationPhotons",
1063  edm::InputTag("muonPUPPINoLeptonsIsolation", "gamma-DR030-ThresholdVeto000-ConeVeto001"),
1064  true)) or
1065  false >> edm::EmptyGroupDescription());
1066 
1067  // Efficiency configurables
1068  edm::ParameterSetDescription efficienciesPSet;
1069  efficienciesPSet.setAllowAnything(); // TODO: the pat helper needs to implement a description.
1070  iDesc.add("efficiencies", efficienciesPSet);
1071  iDesc.add<bool>("addEfficiencies", false);
1072 
1073  // Check to see if the user wants to add user data
1074  edm::ParameterSetDescription userDataPSet;
1075  PATUserDataHelper<Muon>::fillDescription(userDataPSet);
1076  iDesc.addOptional("userData", userDataPSet);
1077 
1078  edm::ParameterSetDescription isolationPSet;
1079  isolationPSet.setAllowAnything(); // TODO: the pat helper needs to implement a description.
1080  iDesc.add("userIsolation", isolationPSet);
1081 
1082  iDesc.add<bool>("embedHighLevelSelection", true)->setComment("embed high level selection");
1083  edm::ParameterSetDescription highLevelPSet;
1084  highLevelPSet.setAllowAnything();
1085  iDesc.addNode(edm::ParameterDescription<edm::InputTag>("beamLineSrc", edm::InputTag(), true))
1086  ->setComment("input with high level selection");
1087  iDesc.addNode(edm::ParameterDescription<edm::InputTag>("pvSrc", edm::InputTag(), true))
1088  ->setComment("input with high level selection");
1089 
1090  //descriptions.add("PATMuonProducer", iDesc);
1091 }
edm::ParameterDescriptionNode::setComment
void setComment(std::string const &value)
Definition: ParameterDescriptionNode.cc:101
edm::ParameterSetDescription::ifValue
ParameterDescriptionNode * ifValue(ParameterDescription< T > const &switchParameter, std::unique_ptr< ParameterDescriptionCases< T >> cases)
Definition: ParameterSetDescription.h:220
edm::ParameterSetDescription::addOptional
ParameterDescriptionBase * addOptional(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:105
edm::ParameterSetDescription::setAllowAnything
void setAllowAnything()
allow any parameter label/value pairs
Definition: ParameterSetDescription.cc:37
edm::ParameterSetDescription::addNode
ParameterDescriptionNode * addNode(ParameterDescriptionNode const &node)
Definition: ParameterSetDescription.cc:41
pat::PATUserDataHelper::fillDescription
static void fillDescription(edm::ParameterSetDescription &iDesc)
Definition: PATUserDataHelper.h:153
edm::ParameterSetDescription::setComment
void setComment(std::string const &value)
Definition: ParameterSetDescription.cc:33
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
pat::helper::KinResolutionsLoader::fillDescription
static void fillDescription(edm::ParameterSetDescription &iDesc)
Method for documentation and validation of PSet.
Definition: KinResolutionsLoader.cc:36
void PATMuonProducer::fillHltTriggerInfo ( pat::Muon &  muon,
edm::Handle< std::vector< pat::TriggerObjectStandAlone >> &  triggerObjects,
const edm::TriggerNames &  names,
const std::vector< std::string > &  collection_names 
)
private

Definition at line 286 of file PATMuonProducer.cc.

References pat::PATObject< ObjectType >::addTriggerObjectMatch(), boostedElectronIsolation_cff::deltaR, metsig::muon, dataset::name, hgcalPlots::obj, trigger::TriggerMuon, TriggerAnalyzer::triggerObjects, and pat::TriggerObjectStandAlone::unpackPathNames().

Referenced by produce().

289  {
290  // WARNING: in a case of close-by muons the dR matching may select both muons.
291  // It's better to select the best match for a given collection.
292  for (const auto& triggerObject : *triggerObjects) {
293  if (triggerObject.hasTriggerObjectType(trigger::TriggerMuon)) {
294  bool keepIt = false;
295  for (const auto& name : collection_filter_names) {
296  if (triggerObject.hasCollection(name)) {
297  keepIt = true;
298  break;
299  }
300  }
301  if (not keepIt)
302  continue;
303  if (deltaR(triggerObject.p4(), muon) > 0.1)
304  continue;
305  pat::TriggerObjectStandAlone obj(triggerObject);
306  obj.unpackPathNames(names);
307  muon.addTriggerObjectMatch(obj);
308  }
309  }
310 }
boostedElectronIsolation_cff.deltaR
deltaR
Definition: boostedElectronIsolation_cff.py:77
pat::PATObject::addTriggerObjectMatch
void addTriggerObjectMatch(const TriggerObjectStandAlone &trigObj)
add a trigger match
Definition: PATObject.h:204
pat::TriggerObjectStandAlone
Analysis-level trigger object class (stand-alone)
Definition: TriggerObjectStandAlone.h:32
void PATMuonProducer::fillL1TriggerInfo ( pat::Muon &  muon,
edm::Handle< std::vector< pat::TriggerObjectStandAlone >> &  triggerObjects,
const edm::TriggerNames &  names,
const edm::ESHandle< GlobalTrackingGeometry > &  geometry 
)
private

Definition at line 234 of file PATMuonProducer.cc.

References funct::abs(), pat::PATObject< ObjectType >::addTriggerObjectMatch(), MuonSubdetId::CSC, hiPixelPairStep_cff::deltaPhi, boostedElectronIsolation_cff::deltaR, MuonSubdetId::DT, getMuonDirection(), reco::Muon::matches(), hgcalPlots::obj, trigger::TriggerL1Mu, TriggerAnalyzer::triggerObjects, and pat::TriggerObjectStandAlone::unpackPathNames().

Referenced by produce().

237  {
238  // L1 trigger object parameters are defined at MB2/ME2. Use the muon
239  // chamber matching information to get the local direction of the
240  // muon trajectory and convert it to a global direction to match the
241  // trigger objects
242 
243  std::optional<GlobalPoint> muonPosition;
244  // Loop over chambers
245  // initialize muonPosition with any available match, just in case
246  // the second station is missing - it's better folling back to
247  // dR matching at IP
248  for (const auto& chamberMatch : aMuon.matches()) {
249  if (chamberMatch.id.subdetId() == MuonSubdetId::DT) {
250  DTChamberId detId(chamberMatch.id.rawId());
251  if (abs(detId.station()) > 3)
252  continue;
253  muonPosition = getMuonDirection(chamberMatch, geometry, detId);
254  if (abs(detId.station()) == 2)
255  break;
256  }
257  if (chamberMatch.id.subdetId() == MuonSubdetId::CSC) {
258  CSCDetId detId(chamberMatch.id.rawId());
259  if (abs(detId.station()) > 3)
260  continue;
261  muonPosition = getMuonDirection(chamberMatch, geometry, detId);
262  if (abs(detId.station()) == 2)
263  break;
264  }
265  }
266  if (not muonPosition)
267  return;
268  for (const auto& triggerObject : *triggerObjects) {
269  if (triggerObject.hasTriggerObjectType(trigger::TriggerL1Mu)) {
270  if (fabs(triggerObject.eta()) < 0.001) {
271  // L1 is defined in X-Y plain
272  if (deltaPhi(triggerObject.phi(), muonPosition->phi()) > 0.1)
273  continue;
274  } else {
275  // 3D L1
276  if (deltaR(triggerObject.p4(), *muonPosition) > 0.15)
277  continue;
278  }
279  pat::TriggerObjectStandAlone obj(triggerObject);
280  obj.unpackPathNames(names);
281  aMuon.addTriggerObjectMatch(obj);
282  }
283  }
284 }
trigger::TriggerL1Mu
enum start value shifted to 81 so as to avoid clashes with PDG codes
Definition: TriggerTypeDefs.h:31
boostedElectronIsolation_cff.deltaR
deltaR
Definition: boostedElectronIsolation_cff.py:77
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
hiPixelPairStep_cff.deltaPhi
deltaPhi
Definition: hiPixelPairStep_cff.py:87
MuonSubdetId::DT
static constexpr int DT
Definition: MuonSubdetId.h:12
MuonSubdetId::CSC
static constexpr int CSC
Definition: MuonSubdetId.h:13
pat::PATMuonProducer::getMuonDirection
std::optional< GlobalPoint > getMuonDirection(const reco::MuonChamberMatch &chamberMatch, const edm::ESHandle< GlobalTrackingGeometry > &geometry, const DetId &chamberId)
Definition: PATMuonProducer.cc:223
pat::TriggerObjectStandAlone
Analysis-level trigger object class (stand-alone)
Definition: TriggerObjectStandAlone.h:32
void PATMuonProducer::fillMuon ( Muon &  aMuon,
const MuonBaseRef &  muonRef,
const reco::CandidateBaseRef &  baseRef,
const GenAssociations &  genMatches,
const IsoDepositMaps &  deposits,
const IsolationValueMaps &  isolationValues 
) const
private

common muon filling, for both the standard and PF2PAT case

Definition at line 785 of file PATMuonProducer.cc.

References addGenMatch_, pat::PATObject< ObjectType >::addGenParticleRef(), edm::contains(), reco::Muon::DYT, efficiencyLoader_, embedBestTrack_, pat::Muon::embedCombinedMuon(), embedCombinedMuon_, pat::Muon::embedDytMuon(), embedDytMuon_, embedGenMatch_, pat::PATObject< ObjectType >::embedGenParticle(), pat::Muon::embedMuonBestTrack(), pat::Muon::embedPickyMuon(), embedPickyMuon_, pat::Muon::embedStandAloneMuon(), embedStandAloneMuon_, pat::Muon::embedTpfmsMuon(), embedTpfmsMuon_, pat::Muon::embedTrack(), embedTrack_, embedTunePBestTrack_, pat::Muon::embedTunePMuonBestTrack(), pat::helper::EfficiencyLoader::enabled(), pat::helper::KinResolutionsLoader::enabled(), plotBeamSpotDB::first, forceEmbedBestTrack_, mps_fire::i, edm::RefToBase< T >::id(), reco::Muon::isAValidMuonTrack(), reco::Muon::isGlobalMuon(), isoDepositLabels_, isolationValueLabels_, gen::n, pat::Muon::pfCandidateRef(), reco::Muon::Picky, resolutionLoader_, pat::helper::EfficiencyLoader::setEfficiencies(), pat::Lepton< LeptonType >::setIsoDeposit(), pat::Lepton< LeptonType >::setIsolation(), reco::LeafCandidate::setP4(), pat::helper::KinResolutionsLoader::setResolutions(), source, reco::Muon::TPFMS, and useParticleFlow_.

Referenced by produce().

790  {
791  // in the particle flow algorithm,
792  // the muon momentum is recomputed.
793  // the new value is stored as the momentum of the
794  // resulting PFCandidate of type Muon, and choosen
795  // as the pat::Muon momentum
796  if (useParticleFlow_)
797  aMuon.setP4(aMuon.pfCandidateRef()->p4());
798  if (embedTrack_)
799  aMuon.embedTrack();
800  if (embedStandAloneMuon_)
801  aMuon.embedStandAloneMuon();
802  if (embedCombinedMuon_)
803  aMuon.embedCombinedMuon();
804 
805  // embed the TeV refit track refs (only available for globalMuons)
806  if (aMuon.isGlobalMuon()) {
807  if (embedPickyMuon_ && aMuon.isAValidMuonTrack(reco::Muon::Picky))
808  aMuon.embedPickyMuon();
809  if (embedTpfmsMuon_ && aMuon.isAValidMuonTrack(reco::Muon::TPFMS))
810  aMuon.embedTpfmsMuon();
811  if (embedDytMuon_ && aMuon.isAValidMuonTrack(reco::Muon::DYT))
812  aMuon.embedDytMuon();
813  }
814 
815  // embed best tracks (at the end, so unless forceEmbedBestTrack_ is true we can save some space not embedding them twice)
816  if (embedBestTrack_)
817  aMuon.embedMuonBestTrack(forceEmbedBestTrack_);
818  if (embedTunePBestTrack_)
819  aMuon.embedTunePMuonBestTrack(forceEmbedBestTrack_);
820 
821  // store the match to the generated final state muons
822  if (addGenMatch_) {
823  for (size_t i = 0, n = genMatches.size(); i < n; ++i) {
824  reco::GenParticleRef genMuon = (*genMatches[i])[baseRef];
825  aMuon.addGenParticleRef(genMuon);
826  }
827  if (embedGenMatch_)
828  aMuon.embedGenParticle();
829  }
830  if (efficiencyLoader_.enabled()) {
831  efficiencyLoader_.setEfficiencies(aMuon, muonRef);
832  }
833 
834  for (size_t j = 0, nd = deposits.size(); j < nd; ++j) {
835  if (useParticleFlow_) {
836  if (deposits[j]->contains(baseRef.id())) {
837  aMuon.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[baseRef]);
838  } else if (deposits[j]->contains(muonRef.id())) {
839  aMuon.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[muonRef]);
840  } else {
841  reco::CandidatePtr source = aMuon.pfCandidateRef()->sourceCandidatePtr(0);
842  aMuon.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[source]);
843  }
844  } else {
845  aMuon.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[muonRef]);
846  }
847  }
848 
849  for (size_t j = 0; j < isolationValues.size(); ++j) {
850  if (useParticleFlow_) {
851  if (isolationValues[j]->contains(baseRef.id())) {
852  aMuon.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[baseRef]);
853  } else if (isolationValues[j]->contains(muonRef.id())) {
854  aMuon.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[muonRef]);
855  } else {
856  reco::CandidatePtr source = aMuon.pfCandidateRef()->sourceCandidatePtr(0);
857  aMuon.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[source]);
858  }
859  } else {
860  aMuon.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[muonRef]);
861  }
862  }
863 
864  if (resolutionLoader_.enabled()) {
865  resolutionLoader_.setResolutions(aMuon);
866  }
867 }
pat::PATMuonProducer::embedTpfmsMuon_
bool embedTpfmsMuon_
embed track from tpfms muon fit into the muon
Definition: PATMuonProducer.h:173
pat::helper::EfficiencyLoader::enabled
bool enabled() const
&#39;true&#39; if this there is at least one efficiency configured
Definition: EfficiencyLoader.h:25
pat::Muon::embedDytMuon
void embedDytMuon()
embed reference to the above dyt Track
pat::Muon::embedTpfmsMuon
void embedTpfmsMuon()
embed reference to the above tpfms Track
mps_fire.i
i
Definition: mps_fire.py:338
pat::Lepton::setIsolation
void setIsolation(IsolationKeys key, float value)
Definition: Lepton.h:99
edm::contains
bool contains(EventRange const &lh, EventID const &rh)
Definition: EventRange.cc:38
pat::Muon::embedCombinedMuon
void embedCombinedMuon()
set reference to Track reconstructed in both tracked and muon detector (reimplemented from reco::Muon...
pat::Muon::embedTunePMuonBestTrack
void embedTunePMuonBestTrack(bool force=false)
pat::Muon::embedMuonBestTrack
void embedMuonBestTrack(bool force=false)
pat::PATMuonProducer::embedBestTrack_
bool embedBestTrack_
embed the track from best muon measurement (global pflow)
Definition: PATMuonProducer.h:151
pat::helper::KinResolutionsLoader::setResolutions
void setResolutions(pat::PATObject< T > &obj) const
Sets the efficiencies for this object, using the reference to the original objects.
Definition: KinResolutionsLoader.h:49
edm::RefToBase::id
ProductID id() const
Definition: RefToBase.h:242
pat::Muon::embedStandAloneMuon
void embedStandAloneMuon()
set reference to Track reconstructed in the muon detector only (reimplemented from reco::Muon) ...
pat::helper::KinResolutionsLoader::enabled
bool enabled() const
&#39;true&#39; if this there is at least one efficiency configured
Definition: KinResolutionsLoader.h:27
pat::Lepton::setIsoDeposit
void setIsoDeposit(IsolationKeys key, const IsoDeposit &dep)
Sets the IsoDeposit associated with some key; if it is already existent, it is overwritten.
Definition: Lepton.h:176
pat::Muon::embedTrack
void embedTrack()
set reference to Track reconstructed in the tracker only (reimplemented from reco::Muon) ...
pat::PATMuonProducer::useParticleFlow_
bool useParticleFlow_
switch to use particle flow (PF2PAT) or not
Definition: PATMuonProducer.h:193
pat::PATMuonProducer::embedStandAloneMuon_
bool embedStandAloneMuon_
embed track from muon system into the muon
Definition: PATMuonProducer.h:159
pat::PATMuonProducer::embedTrack_
bool embedTrack_
embed the track from inner tracker into the muon
Definition: PATMuonProducer.h:157
pat::PATMuonProducer::addGenMatch_
bool addGenMatch_
add generator match information
Definition: PATMuonProducer.h:183
reco::Muon::isGlobalMuon
bool isGlobalMuon() const override
Definition: Muon.h:291
pat::PATObject::embedGenParticle
void embedGenParticle()
Definition: PATObject.h:694
CandIsolatorFromDeposits_cfi.deposits
deposits
Definition: CandIsolatorFromDeposits_cfi.py:4
pat::Muon::pfCandidateRef
reco::PFCandidateRef pfCandidateRef() const
reco::Muon::isAValidMuonTrack
bool isAValidMuonTrack(const MuonTrackType &type) const
pat::PATObject::addGenParticleRef
void addGenParticleRef(const reco::GenParticleRef &ref)
Definition: PATObject.h:678
pat::PATMuonProducer::embedTunePBestTrack_
bool embedTunePBestTrack_
embed the track from best muon measurement (muon only)
Definition: PATMuonProducer.h:153
electronProducer_cff.isolationValues
isolationValues
Definition: electronProducer_cff.py:26
pat::helper::EfficiencyLoader::setEfficiencies
void setEfficiencies(pat::PATObject< T > &obj, const R &originalRef) const
Sets the efficiencies for this object, using the reference to the original objects.
Definition: EfficiencyLoader.h:42
pat::PATMuonProducer::efficiencyLoader_
pat::helper::EfficiencyLoader efficiencyLoader_
helper class to add efficiencies to the muon
Definition: PATMuonProducer.h:244
pat::PATMuonProducer::embedPickyMuon_
bool embedPickyMuon_
embed track from picky muon fit into the muon
Definition: PATMuonProducer.h:171
pat::PATMuonProducer::embedDytMuon_
bool embedDytMuon_
embed track from DYT muon fit into the muon
Definition: PATMuonProducer.h:175
pat::PATMuonProducer::embedGenMatch_
bool embedGenMatch_
embed the gen match information into the muon
Definition: PATMuonProducer.h:187
pat::PATMuonProducer::isoDepositLabels_
IsolationLabels isoDepositLabels_
input source for isoDeposits
Definition: PATMuonProducer.h:205
pat::PATMuonProducer::isolationValueLabels_
IsolationLabels isolationValueLabels_
input source isolation value maps
Definition: PATMuonProducer.h:208
pat::Muon::embedPickyMuon
void embedPickyMuon()
embed reference to the above picky Track
pat::PATMuonProducer::forceEmbedBestTrack_
bool forceEmbedBestTrack_
force separate embed of the best track even if already embedded
Definition: PATMuonProducer.h:155
pat::PATMuonProducer::embedCombinedMuon_
bool embedCombinedMuon_
embed track of the combined fit into the muon
Definition: PATMuonProducer.h:161
pat::PATMuonProducer::resolutionLoader_
pat::helper::KinResolutionsLoader resolutionLoader_
helper class to add resolutions to the muon
Definition: PATMuonProducer.h:191
source
static std::string const source
Definition: EdmProvDump.cc:47
reco::LeafCandidate::setP4
void setP4(const LorentzVector &p4) final
set 4-momentum
Definition: LeafCandidate.h:143
std::optional< GlobalPoint > PATMuonProducer::getMuonDirection ( const reco::MuonChamberMatch &  chamberMatch,
const edm::ESHandle< GlobalTrackingGeometry > &  geometry,
const DetId &  chamberId 
)
private

Definition at line 223 of file PATMuonProducer.cc.

References GlobalTrackingGeometry::idToDet(), GeomDet::toGlobal(), reco::MuonChamberMatch::x, and reco::MuonChamberMatch::y.

Referenced by fillL1TriggerInfo().

225  {
226  const GeomDet* chamberGeometry = geometry->idToDet(chamberId);
227  if (chamberGeometry) {
228  LocalPoint localPosition(chamberMatch.x, chamberMatch.y, 0);
229  return std::optional<GlobalPoint>(std::in_place, chamberGeometry->toGlobal(localPosition));
230  }
231  return std::optional<GlobalPoint>();
232 }
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:54
GlobalTrackingGeometry::idToDet
const GeomDet * idToDet(DetId) const override
Definition: GlobalTrackingGeometry.cc:44
double PATMuonProducer::getRelMiniIsoPUCorrected ( const pat::Muon &  muon,
double  rho,
const std::vector< double > &  area 
)
private

Definition at line 884 of file PATMuonProducer.cc.

References custom_jme_cff::area, MTVHistoProducerAlgoForTrackerBlock_cfi::maxdr, MTVHistoProducerAlgoForTrackerBlock_cfi::mindr, pat::miniIsoDr(), miniIsoParams_, pat::Lepton< LeptonType >::miniPFIsolation(), pat::muonRelMiniIsoPUCorrected(), and reco::LeafCandidate::polarP4().

Referenced by produce().

884  {
885  double mindr(miniIsoParams_[0]);
886  double maxdr(miniIsoParams_[1]);
887  double kt_scale(miniIsoParams_[2]);
888  double drcut = pat::miniIsoDr(muon.polarP4(), mindr, maxdr, kt_scale);
889  return pat::muonRelMiniIsoPUCorrected(muon.miniPFIsolation(), muon.polarP4(), drcut, rho, area);
890 }
reco::LeafCandidate::polarP4
const PolarLorentzVector & polarP4() const final
four-momentum Lorentz vector
Definition: LeafCandidate.h:101
MTVHistoProducerAlgoForTrackerBlock_cfi.mindr
mindr
Definition: MTVHistoProducerAlgoForTrackerBlock_cfi.py:79
pat::PATMuonProducer::miniIsoParams_
std::vector< double > miniIsoParams_
Definition: PATMuonProducer.h:147
pat::Lepton::miniPFIsolation
const PFIsolation & miniPFIsolation() const
Definition: Lepton.h:196
pat::muonRelMiniIsoPUCorrected
double muonRelMiniIsoPUCorrected(const PFIsolation &iso, const reco::Candidate::PolarLorentzVector &p4, double dr, double rho, const std::vector< double > &area)
Definition: MiniIsolation.cc:60
MTVHistoProducerAlgoForTrackerBlock_cfi.maxdr
maxdr
Definition: MTVHistoProducerAlgoForTrackerBlock_cfi.py:80
custom_jme_cff.area
area
Definition: custom_jme_cff.py:133
pat::miniIsoDr
float miniIsoDr(const reco::Candidate::PolarLorentzVector &p4, float mindr, float maxdr, float kt_scale)
Definition: MiniIsolation.cc:15
static void pat::PATMuonProducer::globalEndJob ( PATMuonHeavyObjectCache *  )
inlinestatic

Definition at line 77 of file PATMuonProducer.h.

References fillDescriptions(), and iEvent.

77 {}
static std::unique_ptr<PATMuonHeavyObjectCache> pat::PATMuonProducer::initializeGlobalCache ( const edm::ParameterSet &  iConfig)
inlinestatic

Definition at line 73 of file PATMuonProducer.h.

73  {
74  return std::make_unique<PATMuonHeavyObjectCache>(iConfig);
75  }
bool PATMuonProducer::isChargedHadron ( long  pdgid)
private

Definition at line 942 of file PATMuonProducer.cc.

References funct::abs().

Referenced by puppiCombinedIsolation().

942 { return std::abs(pdgid) == 211; }
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
BPhysicsValidation_cfi.pdgid
pdgid
Definition: BPhysicsValidation_cfi.py:8
bool PATMuonProducer::isNeutralHadron ( long  pdgid)
private

Definition at line 940 of file PATMuonProducer.cc.

References funct::abs().

Referenced by puppiCombinedIsolation().

940 { return std::abs(pdgid) == 130; }
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
BPhysicsValidation_cfi.pdgid
pdgid
Definition: BPhysicsValidation_cfi.py:8
bool PATMuonProducer::isPhoton ( long  pdgid)
private

Definition at line 944 of file PATMuonProducer.cc.

Referenced by puppiCombinedIsolation().

944 { return pdgid == 22; }
BPhysicsValidation_cfi.pdgid
pdgid
Definition: BPhysicsValidation_cfi.py:8
void PATMuonProducer::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
override

everything that needs to be done during the event loop

Definition at line 312 of file PATMuonProducer.cc.

References addGenMatch_, addInverseBeta_, addPuppiIsolation_, addTriggerMatching_, beamLineToken_, ecalDrivenElectronSeedsParameters_cff::beamSpot, pat::helper::MultiIsolator::beginEvent(), TransientTrackBuilder::build(), muonProducer_cfi::caloMETMuonCorrs, caloMETMuonCorrsToken_, computeMiniIso_, computeMuonMVA_, computePuppiCombinedIso_, computeSoftMuonMVA_, CandIsolatorFromDeposits_cfi::deposits, PVValHelper::dz, reco::PFCandidate::ecalEnergy(), effectiveAreaVec_, efficiencyLoader_, pat::Muon::embedCaloMETMuonCorrs(), embedCaloMETMuonCorrs_, embedCombinedMuon_, embedGenMatch_, embedHighLevel(), embedHighLevelSelection_, pat::Muon::embedPFCandidate(), embedPFCandidate_, embedPfEcalEnergy_, pat::Muon::embedTcMETMuonCorrs(), embedTcMETMuonCorrs_, pat::helper::EfficiencyLoader::enabled(), pat::helper::KinResolutionsLoader::enabled(), pat::helper::MultiIsolator::enabled(), pat::helper::MultiIsolator::endEvent(), pat::helper::MultiIsolator::fill(), fillHltTriggerInfo(), fillL1TriggerInfo(), fillMuon(), genMatchTokens_, geometry, edm::EventSetup::get(), edm::Event::getByToken(), getRelMiniIsoPUCorrected(), cutBasedMuonId_MuonPOG_V0_cff::globalTrack, hltCollectionFilters_, mps_fire::i, training_settings::idx, pat::Muon::initSimInfo(), susybsm::HSCParticleType::innerTrack, edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), isoDepositTokens_, electronProducer_cff::isolationValues, isolationValueTokens_, isolator_, isolatorTmpStorage_, edm::EventBase::isRealData(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), reco::Muon::LowPtMvaLoose, reco::Muon::LowPtMvaMedium, muon::makeSelectorBitset(), reco::Muon::MiniIsoLoose, reco::Muon::MiniIsoMedium, reco::Muon::MiniIsoTight, reco::Muon::MiniIsoVeryTight, eostools::move(), reco::Muon::MultiIsoMedium, reco::PFCandidate::muonRef(), extraflags_cff::muons, muonTimeExtraToken_, muonToken_, DetachedQuadStep_cff::mva, mvaBTagCollectionTag_, muonProducer_cfi::mvaL1Corrector, mvaL1Corrector_, muonProducer_cfi::mvaL1L2L3ResCorrector, mvaL1L2L3ResCorrector_, reco::Muon::MvaLoose, reco::Muon::MvaMedium, reco::Muon::MvaTight, mvaUseJec_, reco::Muon::MvaVTight, reco::Muon::MvaVVTight, pat::helper::EfficiencyLoader::newEvent(), pat::helper::KinResolutionsLoader::newEvent(), nhits, hemisphereProducer_cfi::patMuons, pcToken_, pfMuonToken_, reco::Vertex::position(), impactParameterTagInfos_cfi::primaryVertex, edm::Handle< T >::product(), pTComparator_, puppiCombinedIsolation(), PUPPIIsolation_charged_hadrons_, PUPPIIsolation_neutral_hadrons_, PUPPIIsolation_photons_, reco::Muon::PuppiIsoLoose, reco::Muon::PuppiIsoMedium, reco::Muon::PuppiIsoTight, PUPPINoLeptonsIsolation_charged_hadrons_, PUPPINoLeptonsIsolation_neutral_hadrons_, PUPPINoLeptonsIsolation_photons_, edm::Event::put(), MetAnalyzer::pv(), pat::Muon::PV2D, pat::Muon::PV3D, pvToken_, pat::Muon::readTimeExtra(), recomputeBasicSelectors_, resolutionLoader_, rho_, edm::Event::run(), pat::Muon::setExtSimType(), pat::Lepton< LeptonType >::setIsolation(), pat::Muon::setIsolationPUPPI(), pat::Muon::setIsolationPUPPINoLeptons(), setMuonMiniIso(), pat::Muon::setNormChi2(), pat::Muon::setNumberOfValidHits(), pat::Muon::setPFCandidateRef(), pat::Muon::setPfEcalEnergy(), pat::Muon::setSimBX(), pat::Muon::setSimEta(), pat::Muon::setSimFlavour(), pat::Muon::setSimHeaviestMotherFlavour(), pat::Muon::setSimMatchQuality(), pat::Muon::setSimMotherPdgId(), pat::Muon::setSimPdgId(), pat::Muon::setSimPhi(), pat::Muon::setSimProdRho(), pat::Muon::setSimProdZ(), pat::Muon::setSimPt(), pat::Muon::setSimTpEvent(), pat::Muon::setSimType(), simInfo_, reco::Muon::SoftMvaId, muonProducer_cfi::tcMETMuonCorrs, tcMETMuonCorrsToken_, edm::Event::triggerNames(), TriggerAnalyzer::triggerObjects, triggerObjects_, triggerResults, triggerResults_, groupFilesInBlocks::tt, useParticleFlow_, userDataHelper_, and useUserData_.

312  {
313  // get the tracking Geometry
314  edm::ESHandle<GlobalTrackingGeometry> geometry;
315  iSetup.get<GlobalTrackingGeometryRecord>().get(geometry);
316  if (!geometry.isValid())
317  throw cms::Exception("FatalError") << "Unable to find GlobalTrackingGeometryRecord in event!\n";
318 
319  // switch off embedding (in unschedules mode)
320  if (iEvent.isRealData()) {
321  addGenMatch_ = false;
322  embedGenMatch_ = false;
323  }
324 
325  edm::Handle<edm::View<reco::Muon>> muons;
326  iEvent.getByToken(muonToken_, muons);
327 
328  edm::Handle<pat::PackedCandidateCollection> pc;
329  if (computeMiniIso_ || computePuppiCombinedIso_)
330  iEvent.getByToken(pcToken_, pc);
331 
332  // get the ESHandle for the transient track builder,
333  // if needed for high level selection embedding
334  edm::ESHandle<TransientTrackBuilder> trackBuilder;
335 
336  if (isolator_.enabled())
337  isolator_.beginEvent(iEvent, iSetup);
338  if (efficiencyLoader_.enabled())
339  efficiencyLoader_.newEvent(iEvent);
340  if (resolutionLoader_.enabled())
341  resolutionLoader_.newEvent(iEvent, iSetup);
342 
343  IsoDepositMaps deposits(isoDepositTokens_.size());
344  for (size_t j = 0; j < isoDepositTokens_.size(); ++j) {
345  iEvent.getByToken(isoDepositTokens_[j], deposits[j]);
346  }
347 
348  IsolationValueMaps isolationValues(isolationValueTokens_.size());
349  for (size_t j = 0; j < isolationValueTokens_.size(); ++j) {
350  iEvent.getByToken(isolationValueTokens_[j], isolationValues[j]);
351  }
352 
353  //value maps for puppi isolation
354  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_charged_hadrons;
355  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_neutral_hadrons;
356  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_photons;
357  //value maps for puppiNoLeptons isolation
358  edm::Handle<edm::ValueMap<float>> PUPPINoLeptonsIsolation_charged_hadrons;
359  edm::Handle<edm::ValueMap<float>> PUPPINoLeptonsIsolation_neutral_hadrons;
360  edm::Handle<edm::ValueMap<float>> PUPPINoLeptonsIsolation_photons;
361  if (addPuppiIsolation_) {
362  //puppi
363  iEvent.getByToken(PUPPIIsolation_charged_hadrons_, PUPPIIsolation_charged_hadrons);
364  iEvent.getByToken(PUPPIIsolation_neutral_hadrons_, PUPPIIsolation_neutral_hadrons);
365  iEvent.getByToken(PUPPIIsolation_photons_, PUPPIIsolation_photons);
366  //puppiNoLeptons
367  iEvent.getByToken(PUPPINoLeptonsIsolation_charged_hadrons_, PUPPINoLeptonsIsolation_charged_hadrons);
368  iEvent.getByToken(PUPPINoLeptonsIsolation_neutral_hadrons_, PUPPINoLeptonsIsolation_neutral_hadrons);
369  iEvent.getByToken(PUPPINoLeptonsIsolation_photons_, PUPPINoLeptonsIsolation_photons);
370  }
371 
372  // inputs for muon mva
373  edm::Handle<reco::JetTagCollection> mvaBTagCollectionTag;
374  edm::Handle<reco::JetCorrector> mvaL1Corrector;
375  edm::Handle<reco::JetCorrector> mvaL1L2L3ResCorrector;
376  if (computeMuonMVA_) {
377  iEvent.getByToken(mvaBTagCollectionTag_, mvaBTagCollectionTag);
378  iEvent.getByToken(mvaL1Corrector_, mvaL1Corrector);
379  iEvent.getByToken(mvaL1L2L3ResCorrector_, mvaL1L2L3ResCorrector);
380  }
381 
382  // prepare the MC genMatchTokens_
383  GenAssociations genMatches(genMatchTokens_.size());
384  if (addGenMatch_) {
385  for (size_t j = 0, nd = genMatchTokens_.size(); j < nd; ++j) {
386  iEvent.getByToken(genMatchTokens_[j], genMatches[j]);
387  }
388  }
389 
390  // prepare the high level selection: needs beamline
391  // OR primary vertex, depending on user selection
392  reco::Vertex primaryVertex;
393  reco::BeamSpot beamSpot;
394  bool beamSpotIsValid = false;
395  bool primaryVertexIsValid = false;
396  if (embedHighLevelSelection_) {
397  // get the beamspot
398  edm::Handle<reco::BeamSpot> beamSpotHandle;
399  iEvent.getByToken(beamLineToken_, beamSpotHandle);
400 
401  // get the primary vertex
402  edm::Handle<std::vector<reco::Vertex>> pvHandle;
403  iEvent.getByToken(pvToken_, pvHandle);
404 
405  if (beamSpotHandle.isValid()) {
406  beamSpot = *beamSpotHandle;
407  beamSpotIsValid = true;
408  } else {
409  edm::LogError("DataNotAvailable") << "No beam spot available from EventSetup, not adding high level selection \n";
410  }
411  if (pvHandle.isValid() && !pvHandle->empty()) {
412  primaryVertex = pvHandle->at(0);
413  primaryVertexIsValid = true;
414  } else {
415  edm::LogError("DataNotAvailable")
416  << "No primary vertex available from EventSetup, not adding high level selection \n";
417  }
418  // this is needed by the IPTools methods from the tracking group
419  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", trackBuilder);
420  }
421 
422  // MC info
423  edm::Handle<edm::ValueMap<reco::MuonSimInfo>> simInfo;
424  bool simInfoIsAvailalbe = iEvent.getByToken(simInfo_, simInfo);
425 
426  // this will be the new object collection
427  std::vector<Muon>* patMuons = new std::vector<Muon>();
428 
429  edm::Handle<reco::PFCandidateCollection> pfMuons;
430  if (useParticleFlow_) {
431  // get the PFCandidates of type muons
432  iEvent.getByToken(pfMuonToken_, pfMuons);
433 
434  unsigned index = 0;
435  for (reco::PFCandidateConstIterator i = pfMuons->begin(); i != pfMuons->end(); ++i, ++index) {
436  const reco::PFCandidate& pfmu = *i;
437  //const reco::IsolaPFCandidate& pfmu = *i;
438  const reco::MuonRef& muonRef = pfmu.muonRef();
439  assert(muonRef.isNonnull());
440 
441  MuonBaseRef muonBaseRef(muonRef);
442  Muon aMuon(muonBaseRef);
443 
444  if (useUserData_) {
445  userDataHelper_.add(aMuon, iEvent, iSetup);
446  }
447 
448  // embed high level selection
449  if (embedHighLevelSelection_) {
450  // get the tracks
451  reco::TrackRef innerTrack = muonBaseRef->innerTrack();
452  reco::TrackRef globalTrack = muonBaseRef->globalTrack();
453  reco::TrackRef bestTrack = muonBaseRef->muonBestTrack();
454  reco::TrackRef chosenTrack = innerTrack;
455  // Make sure the collection it points to is there
456  if (bestTrack.isNonnull() && bestTrack.isAvailable())
457  chosenTrack = bestTrack;
458 
459  if (chosenTrack.isNonnull() && chosenTrack.isAvailable()) {
460  unsigned int nhits = chosenTrack->numberOfValidHits(); // ????
461  aMuon.setNumberOfValidHits(nhits);
462 
463  reco::TransientTrack tt = trackBuilder->build(chosenTrack);
464  embedHighLevel(aMuon, chosenTrack, tt, primaryVertex, primaryVertexIsValid, beamSpot, beamSpotIsValid);
465  }
466 
467  if (globalTrack.isNonnull() && globalTrack.isAvailable() && !embedCombinedMuon_) {
468  double norm_chi2 = globalTrack->chi2() / globalTrack->ndof();
469  aMuon.setNormChi2(norm_chi2);
470  }
471  }
472  reco::PFCandidateRef pfRef(pfMuons, index);
473  //reco::PFCandidatePtr ptrToMother(pfMuons,index);
474  reco::CandidateBaseRef pfBaseRef(pfRef);
475 
476  aMuon.setPFCandidateRef(pfRef);
477  if (embedPFCandidate_)
478  aMuon.embedPFCandidate();
479  fillMuon(aMuon, muonBaseRef, pfBaseRef, genMatches, deposits, isolationValues);
480 
481  if (computeMiniIso_)
482  setMuonMiniIso(aMuon, pc.product());
483 
484  if (addPuppiIsolation_) {
485  aMuon.setIsolationPUPPI((*PUPPIIsolation_charged_hadrons)[muonBaseRef],
486  (*PUPPIIsolation_neutral_hadrons)[muonBaseRef],
487  (*PUPPIIsolation_photons)[muonBaseRef]);
488 
489  aMuon.setIsolationPUPPINoLeptons((*PUPPINoLeptonsIsolation_charged_hadrons)[muonBaseRef],
490  (*PUPPINoLeptonsIsolation_neutral_hadrons)[muonBaseRef],
491  (*PUPPINoLeptonsIsolation_photons)[muonBaseRef]);
492  } else {
493  aMuon.setIsolationPUPPI(-999., -999., -999.);
494  aMuon.setIsolationPUPPINoLeptons(-999., -999., -999.);
495  }
496 
497  if (embedPfEcalEnergy_) {
498  aMuon.setPfEcalEnergy(pfmu.ecalEnergy());
499  }
500 
501  patMuons->push_back(aMuon);
502  }
503  } else {
504  edm::Handle<edm::View<reco::Muon>> muons;
505  iEvent.getByToken(muonToken_, muons);
506 
507  // embedding of muon MET corrections
508  edm::Handle<edm::ValueMap<reco::MuonMETCorrectionData>> caloMETMuonCorrs;
509  //edm::ValueMap<reco::MuonMETCorrectionData> caloMETmuCorValueMap;
510  if (embedCaloMETMuonCorrs_) {
511  iEvent.getByToken(caloMETMuonCorrsToken_, caloMETMuonCorrs);
512  //caloMETmuCorValueMap = *caloMETmuCorValueMap_h;
513  }
514  edm::Handle<edm::ValueMap<reco::MuonMETCorrectionData>> tcMETMuonCorrs;
515  //edm::ValueMap<reco::MuonMETCorrectionData> tcMETmuCorValueMap;
516  if (embedTcMETMuonCorrs_) {
517  iEvent.getByToken(tcMETMuonCorrsToken_, tcMETMuonCorrs);
518  //tcMETmuCorValueMap = *tcMETmuCorValueMap_h;
519  }
520 
521  if (embedPfEcalEnergy_) {
522  // get the PFCandidates of type muons
523  iEvent.getByToken(pfMuonToken_, pfMuons);
524  }
525 
526  edm::Handle<edm::ValueMap<reco::MuonTimeExtra>> muonsTimeExtra;
527  if (addInverseBeta_) {
528  // get MuonTimerExtra value map
529  iEvent.getByToken(muonTimeExtraToken_, muonsTimeExtra);
530  }
531 
532  for (edm::View<reco::Muon>::const_iterator itMuon = muons->begin(); itMuon != muons->end(); ++itMuon) {
533  // construct the Muon from the ref -> save ref to original object
534  unsigned int idx = itMuon - muons->begin();
535  MuonBaseRef muonRef = muons->refAt(idx);
536  reco::CandidateBaseRef muonBaseRef(muonRef);
537 
538  Muon aMuon(muonRef);
539  fillMuon(aMuon, muonRef, muonBaseRef, genMatches, deposits, isolationValues);
540  if (computeMiniIso_)
541  setMuonMiniIso(aMuon, pc.product());
542  if (addPuppiIsolation_) {
543  aMuon.setIsolationPUPPI((*PUPPIIsolation_charged_hadrons)[muonRef],
544  (*PUPPIIsolation_neutral_hadrons)[muonRef],
545  (*PUPPIIsolation_photons)[muonRef]);
546  aMuon.setIsolationPUPPINoLeptons((*PUPPINoLeptonsIsolation_charged_hadrons)[muonRef],
547  (*PUPPINoLeptonsIsolation_neutral_hadrons)[muonRef],
548  (*PUPPINoLeptonsIsolation_photons)[muonRef]);
549  } else {
550  aMuon.setIsolationPUPPI(-999., -999., -999.);
551  aMuon.setIsolationPUPPINoLeptons(-999., -999., -999.);
552  }
553 
554  // Isolation
555  if (isolator_.enabled()) {
556  //reco::CandidatePtr mother = ptrToMother->sourceCandidatePtr(0);
557  isolator_.fill(*muons, idx, isolatorTmpStorage_);
558  typedef pat::helper::MultiIsolator::IsolationValuePairs IsolationValuePairs;
559  // better to loop backwards, so the vector is resized less times
560  for (IsolationValuePairs::const_reverse_iterator it = isolatorTmpStorage_.rbegin(),
561  ed = isolatorTmpStorage_.rend();
562  it != ed;
563  ++it) {
564  aMuon.setIsolation(it->first, it->second);
565  }
566  }
567 
568  // for (size_t j = 0, nd = deposits.size(); j < nd; ++j) {
569  // aMuon.setIsoDeposit(isoDepositLabels_[j].first,
570  // (*deposits[j])[muonRef]);
571  // }
572 
573  // add sel to selected
574  edm::Ptr<reco::Muon> muonsPtr = muons->ptrAt(idx);
575  if (useUserData_) {
576  userDataHelper_.add(aMuon, iEvent, iSetup);
577  }
578 
579  // embed high level selection
580  if (embedHighLevelSelection_) {
581  // get the tracks
582  reco::TrackRef innerTrack = itMuon->innerTrack();
583  reco::TrackRef globalTrack = itMuon->globalTrack();
584  reco::TrackRef bestTrack = itMuon->muonBestTrack();
585  reco::TrackRef chosenTrack = innerTrack;
586  // Make sure the collection it points to is there
587  if (bestTrack.isNonnull() && bestTrack.isAvailable())
588  chosenTrack = bestTrack;
589  if (chosenTrack.isNonnull() && chosenTrack.isAvailable()) {
590  unsigned int nhits = chosenTrack->numberOfValidHits(); // ????
591  aMuon.setNumberOfValidHits(nhits);
592 
593  reco::TransientTrack tt = trackBuilder->build(chosenTrack);
594  embedHighLevel(aMuon, chosenTrack, tt, primaryVertex, primaryVertexIsValid, beamSpot, beamSpotIsValid);
595  }
596 
597  if (globalTrack.isNonnull() && globalTrack.isAvailable()) {
598  double norm_chi2 = globalTrack->chi2() / globalTrack->ndof();
599  aMuon.setNormChi2(norm_chi2);
600  }
601  }
602 
603  // embed MET muon corrections
604  if (embedCaloMETMuonCorrs_)
605  aMuon.embedCaloMETMuonCorrs((*caloMETMuonCorrs)[muonRef]);
606  if (embedTcMETMuonCorrs_)
607  aMuon.embedTcMETMuonCorrs((*tcMETMuonCorrs)[muonRef]);
608 
609  if (embedPfEcalEnergy_) {
610  aMuon.setPfEcalEnergy(-99.0);
611  for (const reco::PFCandidate& pfmu : *pfMuons) {
612  if (pfmu.muonRef().isNonnull()) {
613  if (pfmu.muonRef().id() != muonRef.id())
614  throw cms::Exception("Configuration")
615  << "Muon reference within PF candidates does not point to the muon collection." << std::endl;
616  if (pfmu.muonRef().key() == muonRef.key()) {
617  aMuon.setPfEcalEnergy(pfmu.ecalEnergy());
618  }
619  }
620  }
621  }
622  if (addInverseBeta_) {
623  aMuon.readTimeExtra((*muonsTimeExtra)[muonRef]);
624  }
625  // MC info
626  aMuon.initSimInfo();
627  if (simInfoIsAvailalbe) {
628  const auto& msi = (*simInfo)[muonBaseRef];
629  aMuon.setSimType(msi.primaryClass);
630  aMuon.setExtSimType(msi.extendedClass);
631  aMuon.setSimFlavour(msi.flavour);
632  aMuon.setSimHeaviestMotherFlavour(msi.heaviestMotherFlavour);
633  aMuon.setSimPdgId(msi.pdgId);
634  aMuon.setSimMotherPdgId(msi.motherPdgId);
635  aMuon.setSimBX(msi.tpBX);
636  aMuon.setSimTpEvent(msi.tpEvent);
637  aMuon.setSimProdRho(msi.vertex.Rho());
638  aMuon.setSimProdZ(msi.vertex.Z());
639  aMuon.setSimPt(msi.p4.pt());
640  aMuon.setSimEta(msi.p4.eta());
641  aMuon.setSimPhi(msi.p4.phi());
642  aMuon.setSimMatchQuality(msi.tpAssoQuality);
643  }
644  patMuons->push_back(aMuon);
645  }
646  }
647 
648  // sort muons in pt
649  std::sort(patMuons->begin(), patMuons->end(), pTComparator_);
650 
651  // Store standard muon selection decisions and jet related
652  // quantaties.
653  // Need a separate loop over muons to have all inputs properly
654  // computed and stored in the object.
655  edm::Handle<double> rho;
656  if (computeMuonMVA_)
657  iEvent.getByToken(rho_, rho);
658  const reco::Vertex* pv(nullptr);
659  if (primaryVertexIsValid)
660  pv = &primaryVertex;
661 
662  edm::Handle<std::vector<pat::TriggerObjectStandAlone>> triggerObjects;
663  edm::Handle<edm::TriggerResults> triggerResults;
664  bool triggerObjectsAvailable = false;
665  bool triggerResultsAvailable = false;
666  if (addTriggerMatching_) {
667  triggerObjectsAvailable = iEvent.getByToken(triggerObjects_, triggerObjects);
668  triggerResultsAvailable = iEvent.getByToken(triggerResults_, triggerResults);
669  }
670 
671  for (auto& muon : *patMuons) {
672  // trigger info
673  if (addTriggerMatching_ and triggerObjectsAvailable and triggerResultsAvailable) {
674  const edm::TriggerNames& triggerNames(iEvent.triggerNames(*triggerResults));
675  fillL1TriggerInfo(muon, triggerObjects, triggerNames, geometry);
676  fillHltTriggerInfo(muon, triggerObjects, triggerNames, hltCollectionFilters_);
677  }
678 
679  if (recomputeBasicSelectors_) {
680  muon.setSelectors(0);
681  bool isRun2016BCDEF = (272728 <= iEvent.run() && iEvent.run() <= 278808);
682  muon.setSelectors(muon::makeSelectorBitset(muon, pv, isRun2016BCDEF));
683  }
684  float miniIsoValue = -1;
685  if (computeMiniIso_) {
686  // MiniIsolation working points
687 
688  miniIsoValue = getRelMiniIsoPUCorrected(muon, *rho, effectiveAreaVec_);
689 
690  muon.setSelector(reco::Muon::MiniIsoLoose, miniIsoValue < 0.40);
691  muon.setSelector(reco::Muon::MiniIsoMedium, miniIsoValue < 0.20);
692  muon.setSelector(reco::Muon::MiniIsoTight, miniIsoValue < 0.10);
693  muon.setSelector(reco::Muon::MiniIsoVeryTight, miniIsoValue < 0.05);
694  }
695 
696  double puppiCombinedIsolationPAT = -1;
697  if (computePuppiCombinedIso_) {
698  puppiCombinedIsolationPAT = puppiCombinedIsolation(muon, pc.product());
699  muon.setSelector(reco::Muon::PuppiIsoLoose, puppiCombinedIsolationPAT < 0.27);
700  muon.setSelector(reco::Muon::PuppiIsoMedium, puppiCombinedIsolationPAT < 0.22);
701  muon.setSelector(reco::Muon::PuppiIsoTight, puppiCombinedIsolationPAT < 0.12);
702  }
703 
704  float jetPtRatio = 0.0;
705  float jetPtRel = 0.0;
706  float mva = 0.0;
707  float mva_lowpt = 0.0;
708  if (computeMuonMVA_ && primaryVertexIsValid && computeMiniIso_) {
709  if (mvaUseJec_) {
710  mva = globalCache()->muonMvaEstimator()->computeMva(muon,
711  primaryVertex,
712  *(mvaBTagCollectionTag.product()),
713  jetPtRatio,
714  jetPtRel,
715  miniIsoValue,
716  &*mvaL1Corrector,
717  &*mvaL1L2L3ResCorrector);
718  mva_lowpt = globalCache()->muonLowPtMvaEstimator()->computeMva(muon,
719  primaryVertex,
720  *(mvaBTagCollectionTag.product()),
721  jetPtRatio,
722  jetPtRel,
723  miniIsoValue,
724  &*mvaL1Corrector,
725  &*mvaL1L2L3ResCorrector);
726 
727  } else {
728  mva = globalCache()->muonMvaEstimator()->computeMva(
729  muon, primaryVertex, *(mvaBTagCollectionTag.product()), jetPtRatio, jetPtRel, miniIsoValue);
730  mva_lowpt = globalCache()->muonLowPtMvaEstimator()->computeMva(
731  muon, primaryVertex, *(mvaBTagCollectionTag.product()), jetPtRatio, jetPtRel, miniIsoValue);
732  }
733 
734  muon.setMvaValue(mva);
735  muon.setLowPtMvaValue(mva_lowpt);
736  muon.setJetPtRatio(jetPtRatio);
737  muon.setJetPtRel(jetPtRel);
738 
739  // multi-isolation
740  if (computeMiniIso_) {
741  muon.setSelector(reco::Muon::MultiIsoMedium,
742  miniIsoValue < 0.11 && (muon.jetPtRatio() > 0.74 || muon.jetPtRel() > 6.8));
743  }
744 
745  // MVA working points
746  // https://twiki.cern.ch/twiki/bin/viewauth/CMS/LeptonMVA
747  double dB2D = fabs(muon.dB(pat::Muon::PV2D));
748  double dB3D = fabs(muon.dB(pat::Muon::PV3D));
749  double edB3D = fabs(muon.edB(pat::Muon::PV3D));
750  double sip3D = edB3D > 0 ? dB3D / edB3D : 0.0;
751  double dz = fabs(muon.muonBestTrack()->dz(primaryVertex.position()));
752 
753  // muon preselection
754  if (muon.pt() > 5 and muon.isLooseMuon() and muon.passed(reco::Muon::MiniIsoLoose) and sip3D < 8.0 and
755  dB2D < 0.05 and dz < 0.1) {
756  muon.setSelector(reco::Muon::MvaLoose, muon.mvaValue() > -0.60);
757  muon.setSelector(reco::Muon::MvaMedium, muon.mvaValue() > -0.20);
758  muon.setSelector(reco::Muon::MvaTight, muon.mvaValue() > 0.15);
759  muon.setSelector(reco::Muon::MvaVTight, muon.mvaValue() > 0.45);
760  muon.setSelector(reco::Muon::MvaVVTight, muon.mvaValue() > 0.9);
761  }
762  if (muon.pt() > 5 and muon.isLooseMuon() and sip3D < 4 and dB2D < 0.5 and dz < 1) {
763  muon.setSelector(reco::Muon::LowPtMvaLoose, muon.lowptMvaValue() > -0.60);
764  muon.setSelector(reco::Muon::LowPtMvaMedium, muon.lowptMvaValue() > -0.20);
765  }
766  }
767 
768  //SOFT MVA
769  if (computeSoftMuonMVA_) {
770  float mva = globalCache()->softMuonMvaEstimator()->computeMva(muon);
771  muon.setSoftMvaValue(mva);
772  //preselection in SoftMuonMvaEstimator.cc
773  muon.setSelector(reco::Muon::SoftMvaId, muon.softMvaValue() > 0.58); //WP choose for bmm4
774  }
775  }
776 
777  // put products in Event
778  std::unique_ptr<std::vector<Muon>> ptr(patMuons);
779  iEvent.put(std::move(ptr));
780 
781  if (isolator_.enabled())
782  isolator_.endEvent();
783 }
pat::PATMuonProducer::PUPPIIsolation_charged_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_charged_hadrons_
Definition: PATMuonProducer.h:219
pat::helper::EfficiencyLoader::enabled
bool enabled() const
&#39;true&#39; if this there is at least one efficiency configured
Definition: EfficiencyLoader.h:25
pat::PATMuonProducer::useUserData_
bool useUserData_
add user data to the muon (this will be data members of th muon even w/o embedding) ...
Definition: PATMuonProducer.h:213
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:575
pat::PATMuonProducer::PUPPIIsolation_photons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_photons_
Definition: PATMuonProducer.h:221
reco::PFCandidate::ecalEnergy
double ecalEnergy() const
return corrected Ecal energy
Definition: PFCandidate.h:222
pat::helper::EfficiencyLoader::newEvent
void newEvent(const edm::Event &event)
To be called for each new event, reads in the ValueMaps for efficiencies.
Definition: EfficiencyLoader.cc:23
muonProducer_cfi.caloMETMuonCorrs
caloMETMuonCorrs
Definition: muonProducer_cfi.py:48
pat::PATMuonProducer::addPuppiIsolation_
bool addPuppiIsolation_
add puppi isolation
Definition: PATMuonProducer.h:217
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:125
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:251
pat::PATMuonProducer::tcMETMuonCorrsToken_
edm::EDGetTokenT< edm::ValueMap< reco::MuonMETCorrectionData > > tcMETMuonCorrsToken_
source of tcMET muon corrections
Definition: PATMuonProducer.h:169
mps_fire.i
i
Definition: mps_fire.py:338
pat::PATMuonProducer::embedTcMETMuonCorrs_
bool embedTcMETMuonCorrs_
embed muon MET correction info for tcMET into the muon
Definition: PATMuonProducer.h:167
pat::PATMuonProducer::pvToken_
edm::EDGetTokenT< std::vector< reco::Vertex > > pvToken_
input source of the primary vertex
Definition: PATMuonProducer.h:203
pat::PATMuonProducer::muonTimeExtraToken_
edm::EDGetTokenT< edm::ValueMap< reco::MuonTimeExtra > > muonTimeExtraToken_
input tag for reading inverse beta
Definition: PATMuonProducer.h:181
pat::PATMuonProducer::PUPPINoLeptonsIsolation_photons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_photons_
Definition: PATMuonProducer.h:225
pat::PATMuonProducer::rho_
edm::EDGetTokenT< double > rho_
Definition: PATMuonProducer.h:234
pat::PATMuonProducer::muonToken_
edm::EDGetTokenT< edm::View< reco::Muon > > muonToken_
input source
Definition: PATMuonProducer.h:140
pat::PATMuonProducer::caloMETMuonCorrsToken_
edm::EDGetTokenT< edm::ValueMap< reco::MuonMETCorrectionData > > caloMETMuonCorrsToken_
source of caloMET muon corrections
Definition: PATMuonProducer.h:165
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
hemisphereProducer_cfi.patMuons
patMuons
Definition: hemisphereProducer_cfi.py:6
pat::PATMuonProducer::embedCaloMETMuonCorrs_
bool embedCaloMETMuonCorrs_
embed muon MET correction info for caloMET into the muon
Definition: PATMuonProducer.h:163
pat::PATMuonProducer::simInfo_
edm::EDGetTokenT< edm::ValueMap< reco::MuonSimInfo > > simInfo_
MC info.
Definition: PATMuonProducer.h:249
pat::PATMuonProducer::setMuonMiniIso
void setMuonMiniIso(pat::Muon &aMuon, const pat::PackedCandidateCollection *pc)
Definition: PATMuonProducer.cc:869
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:20
pat::PATMuonProducer::addTriggerMatching_
bool addTriggerMatching_
Trigger.
Definition: PATMuonProducer.h:252
pat::PATMuonProducer::mvaL1Corrector_
edm::EDGetTokenT< reco::JetCorrector > mvaL1Corrector_
Definition: PATMuonProducer.h:232
pat::PATMuonProducer::PUPPINoLeptonsIsolation_charged_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_charged_hadrons_
Definition: PATMuonProducer.h:223
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:263
reco::PFCandidateConstIterator
PFCandidateCollection::const_iterator PFCandidateConstIterator
iterator
Definition: PFCandidateFwd.h:18
edm::EventBase::isRealData
bool isRealData() const
Definition: EventBase.h:62
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:109
pat::PATMuonProducer::embedHighLevel
void embedHighLevel(pat::Muon &aMuon, reco::TrackRef track, reco::TransientTrack &tt, reco::Vertex &primaryVertex, bool primaryVertexIsValid, reco::BeamSpot &beamspot, bool beamspotIsValid)
Definition: PATMuonProducer.cc:1095
cutBasedMuonId_MuonPOG_V0_cff.globalTrack
globalTrack
Definition: cutBasedMuonId_MuonPOG_V0_cff.py:135
pat::PATMuonProducer::effectiveAreaVec_
std::vector< double > effectiveAreaVec_
Definition: PATMuonProducer.h:146
ecalDrivenElectronSeedsParameters_cff.beamSpot
beamSpot
Definition: ecalDrivenElectronSeedsParameters_cff.py:15
pat::PATMuonProducer::PUPPINoLeptonsIsolation_neutral_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_neutral_hadrons_
Definition: PATMuonProducer.h:224
edm::Ref::id
ProductID id() const
Accessor for product ID.
Definition: Ref.h:257
pat::PATMuonProducer::mvaBTagCollectionTag_
edm::EDGetTokenT< reco::JetTagCollection > mvaBTagCollectionTag_
Definition: PATMuonProducer.h:231
pat::PATMuonProducer::pfMuonToken_
edm::EDGetTokenT< reco::PFCandidateCollection > pfMuonToken_
input source pfCandidates that will be to be transformed into pat::Muons, when using PF2PAT ...
Definition: PATMuonProducer.h:195
pat::PATMuonProducer::isolator_
pat::helper::MultiIsolator isolator_
helper class to add userdefined isolation values to the muon
Definition: PATMuonProducer.h:240
muonProducer_cfi.mvaL1Corrector
mvaL1Corrector
Definition: muonProducer_cfi.py:123
pat::PATMuonProducer::isoDepositTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< IsoDeposit > > > isoDepositTokens_
Definition: PATMuonProducer.h:206
pat::helper::KinResolutionsLoader::enabled
bool enabled() const
&#39;true&#39; if this there is at least one efficiency configured
Definition: KinResolutionsLoader.h:27
pat::PATMuonProducer::beamLineToken_
edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
input source of the primary vertex/beamspot
Definition: PATMuonProducer.h:201
pat::PATMuonProducer::addInverseBeta_
bool addInverseBeta_
add combined inverse beta measurement into the muon
Definition: PATMuonProducer.h:177
pat::PATMuonProducer::useParticleFlow_
bool useParticleFlow_
switch to use particle flow (PF2PAT) or not
Definition: PATMuonProducer.h:193
pat::PATMuonProducer::userDataHelper_
pat::PATUserDataHelper< pat::Muon > userDataHelper_
helper class to add userData to the muon
Definition: PATMuonProducer.h:246
pat::helper::MultiIsolator::enabled
bool enabled() const
True if it has a non null configuration.
Definition: MultiIsolator.h:50
pat::PATMuonProducer::IsoDepositMaps
std::vector< edm::Handle< edm::ValueMap< IsoDeposit > > > IsoDepositMaps
Definition: PATMuonProducer.h:88
pat::PATMuonProducer::triggerResults_
edm::EDGetTokenT< edm::TriggerResults > triggerResults_
Definition: PATMuonProducer.h:254
pat::helper::MultiIsolator::beginEvent
void beginEvent(const edm::Event &event, const edm::EventSetup &eventSetup)
Definition: MultiIsolator.cc:66
pat::PATMuonProducer::pTComparator_
GreaterByPt< Muon > pTComparator_
Definition: PATMuonProducer.h:238
pat::PATMuonProducer::addGenMatch_
bool addGenMatch_
add generator match information
Definition: PATMuonProducer.h:183
Muon
Definition: Muon.py:1
pat::PATMuonProducer::fillL1TriggerInfo
void fillL1TriggerInfo(pat::Muon &muon, edm::Handle< std::vector< pat::TriggerObjectStandAlone >> &triggerObjects, const edm::TriggerNames &names, const edm::ESHandle< GlobalTrackingGeometry > &geometry)
Definition: PATMuonProducer.cc:234
pat::PATMuonProducer::embedPfEcalEnergy_
bool embedPfEcalEnergy_
add ecal PF energy
Definition: PATMuonProducer.h:215
pat::helper::KinResolutionsLoader::newEvent
void newEvent(const edm::Event &event, const edm::EventSetup &setup)
To be called for each new event, reads in the EventSetup object.
Definition: KinResolutionsLoader.cc:28
DetachedQuadStep_cff.mva
mva
Definition: DetachedQuadStep_cff.py:211
CandIsolatorFromDeposits_cfi.deposits
deposits
Definition: CandIsolatorFromDeposits_cfi.py:4
edm::Event::run
RunNumber_t run() const
Definition: Event.h:101
MetAnalyzer.pv
def pv(vc)
Definition: MetAnalyzer.py:7
muonProducer_cfi.mvaL1L2L3ResCorrector
mvaL1L2L3ResCorrector
Definition: muonProducer_cfi.py:124
pat::PATMuonProducer::pcToken_
edm::EDGetTokenT< pat::PackedCandidateCollection > pcToken_
Definition: PATMuonProducer.h:143
extraflags_cff.muons
muons
Definition: extraflags_cff.py:17
triggerResults
static std::string const triggerResults
Definition: EdmProvDump.cc:45
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
pat::PATMuonProducer::puppiCombinedIsolation
double puppiCombinedIsolation(const pat::Muon &muon, const pat::PackedCandidateCollection *pc)
Definition: PATMuonProducer.cc:892
pat::PATMuonProducer::hltCollectionFilters_
std::vector< std::string > hltCollectionFilters_
Definition: PATMuonProducer.h:255
reco::PFCandidate::muonRef
reco::MuonRef muonRef() const
Definition: PFCandidate.cc:459
pat::PATMuonProducer::fillHltTriggerInfo
void fillHltTriggerInfo(pat::Muon &muon, edm::Handle< std::vector< pat::TriggerObjectStandAlone >> &triggerObjects, const edm::TriggerNames &names, const std::vector< std::string > &collection_names)
Definition: PATMuonProducer.cc:286
muon::makeSelectorBitset
reco::Muon::Selector makeSelectorBitset(reco::Muon const &muon, reco::Vertex const *vertex=0, bool run2016_hip_mitigation=false)
Definition: MuonSelectors.cc:1002
impactParameterTagInfos_cfi.primaryVertex
primaryVertex
Definition: impactParameterTagInfos_cfi.py:5
electronProducer_cff.isolationValues
isolationValues
Definition: electronProducer_cff.py:26
edm::Handle::product
T const * product() const
Definition: Handle.h:74
pat::PATMuonProducer::getRelMiniIsoPUCorrected
double getRelMiniIsoPUCorrected(const pat::Muon &muon, double rho, const std::vector< double > &area)
Definition: PATMuonProducer.cc:884
training_settings.idx
idx
Definition: training_settings.py:16
pat::PATMuonProducer::GenAssociations
std::vector< edm::Handle< edm::Association< reco::GenParticleCollection > > > GenAssociations
Definition: PATMuonProducer.h:87
pat::PATMuonProducer::mvaL1L2L3ResCorrector_
edm::EDGetTokenT< reco::JetCorrector > mvaL1L2L3ResCorrector_
Definition: PATMuonProducer.h:233
pat::PATMuonProducer::fillMuon
void fillMuon(Muon &aMuon, const MuonBaseRef &muonRef, const reco::CandidateBaseRef &baseRef, const GenAssociations &genMatches, const IsoDepositMaps &deposits, const IsolationValueMaps &isolationValues) const
common muon filling, for both the standard and PF2PAT case
Definition: PATMuonProducer.cc:785
geometry
ESHandle< TrackerGeometry > geometry
Definition: TkLasBeamFitter.cc:176
pat::helper::MultiIsolator::IsolationValuePairs
std::vector< std::pair< pat::IsolationKeys, float > > IsolationValuePairs
Definition: MultiIsolator.h:16
reco::PFCandidate
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:40
pat::PATMuonProducer::IsolationValueMaps
std::vector< edm::Handle< edm::ValueMap< double > > > IsolationValueMaps
Definition: PATMuonProducer.h:89
pat::PATMuonProducer::efficiencyLoader_
pat::helper::EfficiencyLoader efficiencyLoader_
helper class to add efficiencies to the muon
Definition: PATMuonProducer.h:244
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
pat::PATMuonProducer::isolatorTmpStorage_
pat::helper::MultiIsolator::IsolationValuePairs isolatorTmpStorage_
isolation value pair for temporary storage before being folded into the muon
Definition: PATMuonProducer.h:242
pat::PATMuonProducer::triggerObjects_
edm::EDGetTokenT< std::vector< pat::TriggerObjectStandAlone > > triggerObjects_
Definition: PATMuonProducer.h:253
pat::PATMuonProducer::embedGenMatch_
bool embedGenMatch_
embed the gen match information into the muon
Definition: PATMuonProducer.h:187
pat::PATMuonProducer::genMatchTokens_
std::vector< edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > > > genMatchTokens_
input tags for generator match information
Definition: PATMuonProducer.h:185
edm::ProductID::id
ProductIndex id() const
Definition: ProductID.h:38
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:44
pat::PATMuonProducer::MuonBaseRef
edm::RefToBase< reco::Muon > MuonBaseRef
typedefs for convenience
Definition: PATMuonProducer.h:86
pat::PATMuonProducer::embedPFCandidate_
bool embedPFCandidate_
embed pfCandidates into the muon
Definition: PATMuonProducer.h:197
pat::PATMuonProducer::isolationValueTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueTokens_
Definition: PATMuonProducer.h:209
pat::PATMuonProducer::computeMuonMVA_
bool computeMuonMVA_
standard muon selectors
Definition: PATMuonProducer.h:227
pat::PATMuonProducer::embedCombinedMuon_
bool embedCombinedMuon_
embed track of the combined fit into the muon
Definition: PATMuonProducer.h:161
muonProducer_cfi.tcMETMuonCorrs
tcMETMuonCorrs
Definition: muonProducer_cfi.py:51
pat::PATMuonProducer::resolutionLoader_
pat::helper::KinResolutionsLoader resolutionLoader_
helper class to add resolutions to the muon
Definition: PATMuonProducer.h:191
edm::Event::triggerNames
edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const override
Definition: Event.cc:256
pat::PATMuonProducer::PUPPIIsolation_neutral_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_neutral_hadrons_
Definition: PATMuonProducer.h:220
eostools.move
def move(src, dest)
Definition: eostools.py:511
pat::helper::MultiIsolator::fill
void fill(const edm::View< T > &coll, int idx, IsolationValuePairs &isolations) const
Definition: MultiIsolator.h:82
pat::PATMuonProducer::embedHighLevelSelection_
bool embedHighLevelSelection_
embed high level selection variables
Definition: PATMuonProducer.h:199
double PATMuonProducer::puppiCombinedIsolation ( const pat::Muon &  muon,
const pat::PackedCandidateCollection *  pc 
)
private

Definition at line 892 of file PATMuonProducer.cc.

References funct::abs(), reco::deltaPhi(), reco::deltaR2(), reco::LeafCandidate::eta(), isChargedHadron(), isNeutralHadron(), isPhoton(), LogTrace, dqm::qstatus::OTHER, objects.autophobj::particleType, reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().

Referenced by produce().

892  {
893  double dR_threshold = 0.4;
894  double dR2_threshold = dR_threshold * dR_threshold;
895  double mix_fraction = 0.5;
896  enum particleType { CH = 0, NH = 1, PH = 2, OTHER = 100000 };
897  double val_PuppiWithLep = 0.0;
898  double val_PuppiWithoutLep = 0.0;
899 
900  for (const auto& cand : *pc) { //pat::pat::PackedCandidate loop start
901 
902  const particleType pType =
903  isChargedHadron(cand.pdgId()) ? CH : isNeutralHadron(cand.pdgId()) ? NH : isPhoton(cand.pdgId()) ? PH : OTHER;
904  if (pType == OTHER) {
905  if (cand.pdgId() != 1 && cand.pdgId() != 2 && abs(cand.pdgId()) != 11 && abs(cand.pdgId()) != 13) {
906  LogTrace("PATMuonProducer") << "candidate with PDGID = " << cand.pdgId()
907  << " is not CH/NH/PH/e/mu or 1/2 (and this is removed from isolation calculation)"
908  << std::endl;
909  }
910  continue;
911  }
912  double d_eta = std::abs(cand.eta() - muon.eta());
913  if (d_eta > dR_threshold)
914  continue;
915 
916  double d_phi = std::abs(reco::deltaPhi(cand.phi(), muon.phi()));
917  if (d_phi > dR_threshold)
918  continue;
919 
920  double dR2 = reco::deltaR2(cand, muon);
921  if (dR2 > dR2_threshold)
922  continue;
923  if (pType == CH && dR2 < 0.0001 * 0.0001)
924  continue;
925  if (pType == NH && dR2 < 0.01 * 0.01)
926  continue;
927  if (pType == PH && dR2 < 0.01 * 0.01)
928  continue;
929  val_PuppiWithLep += cand.pt() * cand.puppiWeight();
930  val_PuppiWithoutLep += cand.pt() * cand.puppiWeightNoLep();
931 
932  } //pat::pat::PackedCandidate loop end
933 
934  double reliso_Puppi_withLep = val_PuppiWithLep / muon.pt();
935  double reliso_Puppi_withoutlep = val_PuppiWithoutLep / muon.pt();
936  double reliso_Puppi_combined = mix_fraction * reliso_Puppi_withLep + (1.0 - mix_fraction) * reliso_Puppi_withoutlep;
937  return reliso_Puppi_combined;
938 }
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:22
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:137
dqm::qstatus::OTHER
static const int OTHER
Definition: DQMDefinitions.h:12
pat::PATMuonProducer::isChargedHadron
bool isChargedHadron(long pdgid)
Definition: PATMuonProducer.cc:942
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:671
reco::deltaR2
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16
pat::PATMuonProducer::isPhoton
bool isPhoton(long pdgid)
Definition: PATMuonProducer.cc:944
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:133
pat::PATMuonProducer::isNeutralHadron
bool isNeutralHadron(long pdgid)
Definition: PATMuonProducer.cc:940
template<typename T >
void pat::PATMuonProducer::readIsolationLabels ( const edm::ParameterSet &  iConfig,
const char *  psetName,
IsolationLabels &  labels,
std::vector< edm::EDGetTokenT< edm::ValueMap< T >>> &  tokens 
)
private

fill label vector from the contents of the parameter set, for the embedding of isoDeposits or userIsolation values

Definition at line 261 of file PATMuonProducer.h.

References pat::EcalIso, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), pat::HcalIso, crabWrapper::key, label, pat::PfAllParticleIso, pat::PfChargedAllIso, pat::PfChargedHadronIso, pat::PfGammaIso, pat::PfNeutralHadronIso, pat::PfPUChargedHadronIso, pat::TrackIso, pat::UserBaseIso, and edm::vector_transform().

Referenced by PATMuonProducer().

264  {
265  labels.clear();
266 
267  if (iConfig.exists(psetName)) {
268  edm::ParameterSet depconf = iConfig.getParameter<edm::ParameterSet>(psetName);
269 
270  if (depconf.exists("tracker"))
271  labels.push_back(std::make_pair(pat::TrackIso, depconf.getParameter<edm::InputTag>("tracker")));
272  if (depconf.exists("ecal"))
273  labels.push_back(std::make_pair(pat::EcalIso, depconf.getParameter<edm::InputTag>("ecal")));
274  if (depconf.exists("hcal"))
275  labels.push_back(std::make_pair(pat::HcalIso, depconf.getParameter<edm::InputTag>("hcal")));
276  if (depconf.exists("pfAllParticles")) {
277  labels.push_back(std::make_pair(pat::PfAllParticleIso, depconf.getParameter<edm::InputTag>("pfAllParticles")));
278  }
279  if (depconf.exists("pfChargedHadrons")) {
280  labels.push_back(
281  std::make_pair(pat::PfChargedHadronIso, depconf.getParameter<edm::InputTag>("pfChargedHadrons")));
282  }
283  if (depconf.exists("pfChargedAll")) {
284  labels.push_back(std::make_pair(pat::PfChargedAllIso, depconf.getParameter<edm::InputTag>("pfChargedAll")));
285  }
286  if (depconf.exists("pfPUChargedHadrons")) {
287  labels.push_back(
288  std::make_pair(pat::PfPUChargedHadronIso, depconf.getParameter<edm::InputTag>("pfPUChargedHadrons")));
289  }
290  if (depconf.exists("pfNeutralHadrons")) {
291  labels.push_back(
292  std::make_pair(pat::PfNeutralHadronIso, depconf.getParameter<edm::InputTag>("pfNeutralHadrons")));
293  }
294  if (depconf.exists("pfPhotons")) {
295  labels.push_back(std::make_pair(pat::PfGammaIso, depconf.getParameter<edm::InputTag>("pfPhotons")));
296  }
297  if (depconf.exists("user")) {
298  std::vector<edm::InputTag> userdeps = depconf.getParameter<std::vector<edm::InputTag>>("user");
299  std::vector<edm::InputTag>::const_iterator it = userdeps.begin(), ed = userdeps.end();
300  int key = pat::IsolationKeys::UserBaseIso;
301  for (; it != ed; ++it, ++key) {
302  labels.push_back(std::make_pair(pat::IsolationKeys(key), *it));
303  }
304  tokens = edm::vector_transform(
305  labels, [this](IsolationLabel const& label) { return consumes<edm::ValueMap<T>>(label.second); });
306  }
307  }
308  tokens = edm::vector_transform(labels, [this](pat::PATMuonProducer::IsolationLabel const& label) {
309  return consumes<edm::ValueMap<T>>(label.second);
310  });
311 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:674
pat::IsolationKeys
IsolationKeys
Enum defining isolation keys.
Definition: Isolation.h:9
label
char const * label
Definition: PFTauDecayModeTools.cc:11
pat::PATMuonProducer::IsolationLabel
std::pair< pat::IsolationKeys, edm::InputTag > IsolationLabel
Definition: PATMuonProducer.h:90
edm::vector_transform
auto vector_transform(std::vector< InputType > const &input, Function predicate) -> std::vector< typename std::remove_cv< typename std::remove_reference< decltype(predicate(input.front()))>::type >::type >
Definition: transform.h:11
double pat::PATMuonProducer::relMiniIsoPUCorrected ( const pat::Muon &  aMuon,
double  rho 
)
private
void PATMuonProducer::setMuonMiniIso ( pat::Muon &  aMuon,
const pat::PackedCandidateCollection *  pc 
)
private

Definition at line 869 of file PATMuonProducer.cc.

References pat::getMiniPFIsolation(), miniIsoParams_, reco::LeafCandidate::polarP4(), and pat::Lepton< LeptonType >::setMiniPFIsolation().

Referenced by produce().

869  {
870  pat::PFIsolation miniiso = pat::getMiniPFIsolation(pc,
871  aMuon.polarP4(),
872  miniIsoParams_[0],
873  miniIsoParams_[1],
874  miniIsoParams_[2],
875  miniIsoParams_[3],
876  miniIsoParams_[4],
877  miniIsoParams_[5],
878  miniIsoParams_[6],
879  miniIsoParams_[7],
880  miniIsoParams_[8]);
881  aMuon.setMiniPFIsolation(miniiso);
882 }
reco::LeafCandidate::polarP4
const PolarLorentzVector & polarP4() const final
four-momentum Lorentz vector
Definition: LeafCandidate.h:101
pat::PATMuonProducer::miniIsoParams_
std::vector< double > miniIsoParams_
Definition: PATMuonProducer.h:147
pat::Lepton::setMiniPFIsolation
void setMiniPFIsolation(PFIsolation const &iso)
Definition: Lepton.h:197
pat::getMiniPFIsolation
PFIsolation getMiniPFIsolation(const pat::PackedCandidateCollection *pfcands, const reco::Candidate::PolarLorentzVector &p4, float mindr=0.05, float maxdr=0.2, float kt_scale=10.0, float ptthresh=0.5, float deadcone_ch=0.0001, float deadcone_pu=0.01, float deadcone_ph=0.01, float deadcone_nh=0.01, float dZ_cut=0.0)
Definition: MiniIsolation.cc:19

Member Data Documentation

bool pat::PATMuonProducer::addEfficiencies_
private

add efficiencies to the muon (this will be data members of th muon even w/o embedding)

Definition at line 211 of file PATMuonProducer.h.

Referenced by PATMuonProducer().

bool pat::PATMuonProducer::addGenMatch_
private

add generator match information

Definition at line 183 of file PATMuonProducer.h.

Referenced by fillMuon(), PATMuonProducer(), and produce().

bool pat::PATMuonProducer::addInverseBeta_
private

add combined inverse beta measurement into the muon

Definition at line 177 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::addPuppiIsolation_
private

add puppi isolation

Definition at line 217 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::addResolutions_
private

add resolutions to the muon (this will be data members of th muon even w/o embedding)

Definition at line 189 of file PATMuonProducer.h.

Referenced by PATMuonProducer().

bool pat::PATMuonProducer::addTriggerMatching_
private

Trigger.

Definition at line 252 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<reco::BeamSpot> pat::PATMuonProducer::beamLineToken_
private

input source of the primary vertex/beamspot

Definition at line 201 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<reco::MuonMETCorrectionData> > pat::PATMuonProducer::caloMETMuonCorrsToken_
private

source of caloMET muon corrections

Definition at line 165 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::computeMiniIso_
private

Definition at line 144 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::computeMuonMVA_
private

standard muon selectors

Definition at line 227 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::computePuppiCombinedIso_
private

Definition at line 145 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::computeSoftMuonMVA_
private

Definition at line 228 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

std::vector<double> pat::PATMuonProducer::effectiveAreaVec_
private

Definition at line 146 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

pat::helper::EfficiencyLoader pat::PATMuonProducer::efficiencyLoader_
private

helper class to add efficiencies to the muon

Definition at line 244 of file PATMuonProducer.h.

Referenced by fillMuon(), PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedBestTrack_
private

embed the track from best muon measurement (global pflow)

Definition at line 151 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

bool pat::PATMuonProducer::embedCaloMETMuonCorrs_
private

embed muon MET correction info for caloMET into the muon

Definition at line 163 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedCombinedMuon_
private

embed track of the combined fit into the muon

Definition at line 161 of file PATMuonProducer.h.

Referenced by fillMuon(), PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedDytMuon_
private

embed track from DYT muon fit into the muon

Definition at line 175 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

bool pat::PATMuonProducer::embedGenMatch_
private

embed the gen match information into the muon

Definition at line 187 of file PATMuonProducer.h.

Referenced by fillMuon(), PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedHighLevelSelection_
private

embed high level selection variables

Definition at line 199 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedPFCandidate_
private

embed pfCandidates into the muon

Definition at line 197 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedPfEcalEnergy_
private

add ecal PF energy

Definition at line 215 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedPickyMuon_
private

embed track from picky muon fit into the muon

Definition at line 171 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

bool pat::PATMuonProducer::embedStandAloneMuon_
private

embed track from muon system into the muon

Definition at line 159 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

bool pat::PATMuonProducer::embedTcMETMuonCorrs_
private

embed muon MET correction info for tcMET into the muon

Definition at line 167 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::embedTpfmsMuon_
private

embed track from tpfms muon fit into the muon

Definition at line 173 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

bool pat::PATMuonProducer::embedTrack_
private

embed the track from inner tracker into the muon

Definition at line 157 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

bool pat::PATMuonProducer::embedTunePBestTrack_
private

embed the track from best muon measurement (muon only)

Definition at line 153 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

bool pat::PATMuonProducer::forceEmbedBestTrack_
private

force separate embed of the best track even if already embedded

Definition at line 155 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

std::vector<edm::EDGetTokenT<edm::Association<reco::GenParticleCollection> > > pat::PATMuonProducer::genMatchTokens_
private

input tags for generator match information

Definition at line 185 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::getdBFromTrack_
private

switch on reading the dB information from the track

Definition at line 179 of file PATMuonProducer.h.

Referenced by embedHighLevel(), and PATMuonProducer().

std::vector<std::string> pat::PATMuonProducer::hltCollectionFilters_
private

Definition at line 255 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

IsolationLabels pat::PATMuonProducer::isoDepositLabels_
private

input source for isoDeposits

Definition at line 205 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

std::vector<edm::EDGetTokenT<edm::ValueMap<IsoDeposit> > > pat::PATMuonProducer::isoDepositTokens_
private

Definition at line 206 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

IsolationLabels pat::PATMuonProducer::isolationValueLabels_
private

input source isolation value maps

Definition at line 208 of file PATMuonProducer.h.

Referenced by fillMuon(), and PATMuonProducer().

std::vector<edm::EDGetTokenT<edm::ValueMap<double> > > pat::PATMuonProducer::isolationValueTokens_
private

Definition at line 209 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

pat::helper::MultiIsolator pat::PATMuonProducer::isolator_
private

helper class to add userdefined isolation values to the muon

Definition at line 240 of file PATMuonProducer.h.

Referenced by produce().

pat::helper::MultiIsolator::IsolationValuePairs pat::PATMuonProducer::isolatorTmpStorage_
private

isolation value pair for temporary storage before being folded into the muon

Definition at line 242 of file PATMuonProducer.h.

Referenced by produce().

std::vector<double> pat::PATMuonProducer::miniIsoParams_
private

Definition at line 147 of file PATMuonProducer.h.

Referenced by getRelMiniIsoPUCorrected(), PATMuonProducer(), and setMuonMiniIso().

edm::EDGetTokenT<edm::ValueMap<reco::MuonTimeExtra> > pat::PATMuonProducer::muonTimeExtraToken_
private

input tag for reading inverse beta

Definition at line 181 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::View<reco::Muon> > pat::PATMuonProducer::muonToken_
private

input source

Definition at line 140 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<reco::JetTagCollection> pat::PATMuonProducer::mvaBTagCollectionTag_
private

Definition at line 231 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<reco::JetCorrector> pat::PATMuonProducer::mvaL1Corrector_
private

Definition at line 232 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<reco::JetCorrector> pat::PATMuonProducer::mvaL1L2L3ResCorrector_
private

Definition at line 233 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::mvaUseJec_
private

Definition at line 230 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<pat::PackedCandidateCollection> pat::PATMuonProducer::pcToken_
private

Definition at line 143 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<reco::PFCandidateCollection> pat::PATMuonProducer::pfMuonToken_
private

input source pfCandidates that will be to be transformed into pat::Muons, when using PF2PAT

Definition at line 195 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

GreaterByPt<Muon> pat::PATMuonProducer::pTComparator_
private

— tools — comparator for pt ordering

Definition at line 238 of file PATMuonProducer.h.

Referenced by produce().

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATMuonProducer::PUPPIIsolation_charged_hadrons_
private

Definition at line 219 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATMuonProducer::PUPPIIsolation_neutral_hadrons_
private

Definition at line 220 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATMuonProducer::PUPPIIsolation_photons_
private

Definition at line 221 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATMuonProducer::PUPPINoLeptonsIsolation_charged_hadrons_
private

Definition at line 223 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATMuonProducer::PUPPINoLeptonsIsolation_neutral_hadrons_
private

Definition at line 224 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATMuonProducer::PUPPINoLeptonsIsolation_photons_
private

Definition at line 225 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<std::vector<reco::Vertex> > pat::PATMuonProducer::pvToken_
private

input source of the primary vertex

Definition at line 203 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::recomputeBasicSelectors_
private

Definition at line 229 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

double pat::PATMuonProducer::relMiniIsoPUCorrected_
private

Definition at line 148 of file PATMuonProducer.h.

pat::helper::KinResolutionsLoader pat::PATMuonProducer::resolutionLoader_
private

helper class to add resolutions to the muon

Definition at line 191 of file PATMuonProducer.h.

Referenced by fillMuon(), PATMuonProducer(), and produce().

edm::EDGetTokenT<double> pat::PATMuonProducer::rho_
private

Definition at line 234 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<reco::MuonSimInfo> > pat::PATMuonProducer::simInfo_
private

MC info.

Definition at line 249 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<reco::MuonMETCorrectionData> > pat::PATMuonProducer::tcMETMuonCorrsToken_
private

source of tcMET muon corrections

Definition at line 169 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<std::vector<pat::TriggerObjectStandAlone> > pat::PATMuonProducer::triggerObjects_
private

Definition at line 253 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

edm::EDGetTokenT<edm::TriggerResults> pat::PATMuonProducer::triggerResults_
private

Definition at line 254 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::useParticleFlow_
private

switch to use particle flow (PF2PAT) or not

Definition at line 193 of file PATMuonProducer.h.

Referenced by fillMuon(), PATMuonProducer(), and produce().

pat::PATUserDataHelper<pat::Muon> pat::PATMuonProducer::userDataHelper_
private

helper class to add userData to the muon

Definition at line 246 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

bool pat::PATMuonProducer::useUserData_
private

add user data to the muon (this will be data members of th muon even w/o embedding)

Definition at line 213 of file PATMuonProducer.h.

Referenced by PATMuonProducer(), and produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11