CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
pat::LeptonUpdater< T > Class Template Reference
Inheritance diagram for pat::LeptonUpdater< T >:
edm::global::EDProducer<> edm::global::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 LeptonUpdater (const edm::ParameterSet &iConfig)
 
const std::vector< double > & miniIsoParams (const T &lep) const
 
template<>
const std::vector< double > & miniIsoParams (const pat::Electron &lep) const
 
void produce (edm::StreamID, edm::Event &, edm::EventSetup const &) const override
 
void readMiniIsoParams (const edm::ParameterSet &iConfig)
 
template<>
void readMiniIsoParams (const edm::ParameterSet &iConfig)
 
void recomputeMuonBasicSelectors (T &, const reco::Vertex &, const bool) const
 
template<>
void recomputeMuonBasicSelectors (pat::Muon &lep, const reco::Vertex &pv, const bool do_hip_mitigation_2016) const
 
void setDZ (T &lep, const reco::Vertex &pv) const
 
template<>
void setDZ (pat::Electron &anElectron, const reco::Vertex &pv) const
 
template<>
void setDZ (pat::Muon &aMuon, const reco::Vertex &pv) const
 
 ~LeptonUpdater () override
 
- Public Member Functions inherited from edm::global::EDProducer<>
 EDProducer ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsStreamLuminosityBlocks () const final
 
bool wantsStreamRuns () const final
 
- Public Member Functions inherited from edm::global::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Attributes

edm::EDGetTokenT< reco::BeamSpotbeamLineToken_
 
bool computeMiniIso_
 
bool fixDxySign_
 
std::vector< double > miniIsoParams_ [2]
 
edm::EDGetTokenT< pat::PackedCandidateCollectionpcToken_
 
bool recomputeMuonBasicSelectors_
 
edm::EDGetTokenT< std::vector< T > > src_
 
edm::EDGetTokenT< std::vector< reco::Vertex > > vertices_
 

Additional Inherited Members

- Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

template<typename T>
class pat::LeptonUpdater< T >

Definition at line 18 of file LeptonUpdater.cc.

Constructor & Destructor Documentation

template<typename T >
pat::LeptonUpdater< T >::LeptonUpdater ( const edm::ParameterSet iConfig)
inlineexplicit

Definition at line 20 of file LeptonUpdater.cc.

References pat::LeptonUpdater< T >::computeMiniIso_, edm::ParameterSet::getParameter(), pat::LeptonUpdater< T >::pcToken_, pat::LeptonUpdater< T >::readMiniIsoParams(), and pat::LeptonUpdater< T >::recomputeMuonBasicSelectors_.

21  : src_(consumes<std::vector<T>>(iConfig.getParameter<edm::InputTag>("src"))),
22  vertices_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("vertices"))),
23  beamLineToken_(consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamspot"))),
24  computeMiniIso_(iConfig.getParameter<bool>("computeMiniIso")),
25  fixDxySign_(iConfig.getParameter<bool>("fixDxySign")) {
26  //for mini-isolation calculation
27  if (computeMiniIso_) {
28  readMiniIsoParams(iConfig);
29  pcToken_ = consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfCandsForMiniIso"));
30  }
32  if (typeid(T) == typeid(pat::Muon))
33  recomputeMuonBasicSelectors_ = iConfig.getParameter<bool>("recomputeMuonBasicSelectors");
34  produces<std::vector<T>>();
35  }
T getParameter(std::string const &) const
edm::EDGetTokenT< std::vector< T > > src_
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
bool recomputeMuonBasicSelectors_
edm::EDGetTokenT< pat::PackedCandidateCollection > pcToken_
void readMiniIsoParams(const edm::ParameterSet &iConfig)
long double T
Analysis-level muon class.
Definition: Muon.h:51
edm::EDGetTokenT< std::vector< reco::Vertex > > vertices_
template<typename T >
pat::LeptonUpdater< T >::~LeptonUpdater ( )
inlineoverride

Definition at line 37 of file LeptonUpdater.cc.

References pat::LeptonUpdater< T >::produce().

37 {}

Member Function Documentation

template<typename T >
static void pat::LeptonUpdater< T >::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
inlinestatic

Definition at line 41 of file LeptonUpdater.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and edm::ParameterSetDescription::addOptional().

41  {
43  desc.add<edm::InputTag>("src")->setComment("Lepton collection");
44  desc.add<edm::InputTag>("vertices")->setComment("Vertex collection");
45  desc.add<edm::InputTag>("beamspot", edm::InputTag("offlineBeamSpot"))->setComment("Beam spot");
46  desc.add<bool>("computeMiniIso", false)->setComment("Recompute miniIsolation");
47  desc.add<bool>("fixDxySign", false)->setComment("Fix the IP sign");
48  desc.addOptional<edm::InputTag>("pfCandsForMiniIso", edm::InputTag("packedPFCandidates"))
49  ->setComment("PackedCandidate collection used for miniIso");
50  if (typeid(T) == typeid(pat::Muon)) {
51  desc.add<bool>("recomputeMuonBasicSelectors", false)
52  ->setComment("Recompute basic cut-based muon selector flags");
53  desc.addOptional<std::vector<double>>("miniIsoParams")
54  ->setComment("Parameters used for miniIso (as in PATMuonProducer)");
55  descriptions.add("muonsUpdated", desc);
56  } else if (typeid(T) == typeid(pat::Electron)) {
57  desc.addOptional<std::vector<double>>("miniIsoParamsB")
58  ->setComment("Parameters used for miniIso in the barrel (as in PATElectronProducer)");
59  desc.addOptional<std::vector<double>>("miniIsoParamsE")
60  ->setComment("Parameters used for miniIso in the endcap (as in PATElectronProducer)");
61  descriptions.add("electronsUpdated", desc);
62  }
63  }
ParameterDescriptionBase * addOptional(U const &iLabel, T const &value)
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Analysis-level electron class.
Definition: Electron.h:52
void add(std::string const &label, ParameterSetDescription const &psetDescription)
long double T
Analysis-level muon class.
Definition: Muon.h:51
template<typename T >
const std::vector<double>& pat::LeptonUpdater< T >::miniIsoParams ( const T lep) const
inline
template<>
const std::vector< double > & pat::LeptonUpdater< pat::Electron >::miniIsoParams ( const pat::Electron lep) const

Definition at line 112 of file LeptonUpdater.cc.

References reco::GsfElectron::isEE(), and pat::LeptonUpdater< T >::miniIsoParams_.

112  {
113  return miniIsoParams_[lep.isEE()];
114  }
bool isEE() const
Definition: GsfElectron.h:357
std::vector< double > miniIsoParams_[2]
template<typename T >
void pat::LeptonUpdater< T >::produce ( edm::StreamID  ,
edm::Event iEvent,
edm::EventSetup const &   
) const
overridevirtual

Implements edm::global::EDProducerBase.

Definition at line 131 of file LeptonUpdater.cc.

References funct::abs(), pat::LeptonUpdater< T >::beamLineToken_, ecalDrivenElectronSeedsParameters_cff::beamSpot, pat::LeptonUpdater< T >::computeMiniIso_, constexpr, pat::LeptonUpdater< T >::fixDxySign_, edm::Event::getByToken(), pat::getMiniPFIsolation(), mps_fire::i, edm::HandleBase::isValid(), pat::LeptonUpdater< T >::miniIsoParams(), eostools::move(), gen::n, MillePedeFileConverter_cfg::out, pat::LeptonUpdater< T >::pcToken_, reco::Vertex::position(), edm::Handle< T >::product(), edm::Event::put(), MetAnalyzer::pv(), pat::LeptonUpdater< T >::recomputeMuonBasicSelectors(), pat::LeptonUpdater< T >::recomputeMuonBasicSelectors_, edm::Event::run(), pat::LeptonUpdater< T >::setDZ(), TrackRefitter_38T_cff::src, pat::LeptonUpdater< T >::src_, electrons_cff::vertices, and pat::LeptonUpdater< T >::vertices_.

Referenced by pat::LeptonUpdater< T >::~LeptonUpdater().

131  {
133  iEvent.getByToken(src_, src);
134 
136  iEvent.getByToken(vertices_, vertices);
137  const reco::Vertex &pv = vertices->front();
138 
140  if (computeMiniIso_)
141  iEvent.getByToken(pcToken_, pc);
142 
143  edm::Handle<reco::BeamSpot> beamSpotHandle;
144  iEvent.getByToken(beamLineToken_, beamSpotHandle);
146  bool beamSpotIsValid = false;
147  if (beamSpotHandle.isValid()) {
148  beamSpot = *beamSpotHandle;
149  beamSpotIsValid = true;
150  } else {
151  edm::LogError("DataNotAvailable") << "No beam spot available \n";
152  }
153 
154  std::unique_ptr<std::vector<T>> out(new std::vector<T>(*src));
155 
156  const bool do_hip_mitigation_2016 =
157  recomputeMuonBasicSelectors_ && (272728 <= iEvent.run() && iEvent.run() <= 278808);
158 
159  for (unsigned int i = 0, n = src->size(); i < n; ++i) {
160  T &lep = (*out)[i];
161  setDZ(lep, pv);
162  if (computeMiniIso_) {
163  const auto &params = miniIsoParams(lep);
165  lep.polarP4(),
166  params[0],
167  params[1],
168  params[2],
169  params[3],
170  params[4],
171  params[5],
172  params[6],
173  params[7],
174  params[8]);
175  lep.setMiniPFIsolation(miniiso);
176  }
178  recomputeMuonBasicSelectors(lep, pv, do_hip_mitigation_2016);
179  //Fixing the sign of impact parameters
180  if (fixDxySign_) {
181  float signPV = 1.;
182  float signBS = 1.;
183  if (beamSpotIsValid) {
184  if constexpr (std::is_same_v<T, pat::Electron>)
185  signBS = copysign(1., lep.gsfTrack()->dxy(beamSpot));
186  else
187  signBS = copysign(1., lep.bestTrack()->dxy(beamSpot));
188  }
189  if constexpr (std::is_same_v<T, pat::Electron>)
190  signPV = copysign(1., lep.gsfTrack()->dxy(pv.position()));
191  else
192  signPV = copysign(1., lep.bestTrack()->dxy(pv.position()));
193  lep.setDB(abs(lep.dB(T::PV2D)) * signPV, lep.edB(T::PV2D), T::PV2D);
194  lep.setDB(abs(lep.dB(T::BS2D)) * signBS, lep.edB(T::BS2D), T::BS2D);
195  }
196  }
197 
198  iEvent.put(std::move(out));
199 }
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
edm::EDGetTokenT< std::vector< T > > src_
Definition: HeavyIon.h:7
void put(edm::Event &evt, double value, const char *instanceName)
edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
bool recomputeMuonBasicSelectors_
RunNumber_t run() const
Definition: Event.h:101
def pv(vc)
Definition: MetAnalyzer.py:7
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
void setDZ(T &lep, const reco::Vertex &pv) const
bool isValid() const
Definition: HandleBase.h:74
void recomputeMuonBasicSelectors(T &, const reco::Vertex &, const bool) const
edm::EDGetTokenT< pat::PackedCandidateCollection > pcToken_
T const * product() const
Definition: Handle.h:74
const std::vector< double > & miniIsoParams(const T &lep) const
static int position[264][3]
Definition: ReadPGInfo.cc:509
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)
long double T
edm::EDGetTokenT< std::vector< reco::Vertex > > vertices_
def move(src, dest)
Definition: eostools.py:511
#define constexpr
template<typename T >
void pat::LeptonUpdater< T >::readMiniIsoParams ( const edm::ParameterSet iConfig)
inline

Definition at line 67 of file LeptonUpdater.cc.

References edm::ParameterSet::getParameter(), pat::LeptonUpdater< T >::miniIsoParams_, and findQualityFiles::size.

Referenced by pat::LeptonUpdater< T >::LeptonUpdater().

67  {
68  miniIsoParams_[0] = iConfig.getParameter<std::vector<double>>("miniIsoParams");
69  if (miniIsoParams_[0].size() != 9)
70  throw cms::Exception("ParameterError", "miniIsoParams must have exactly 9 elements.\n");
71  }
size
Write out results.
T getParameter(std::string const &) const
std::vector< double > miniIsoParams_[2]
template<>
void pat::LeptonUpdater< pat::Electron >::readMiniIsoParams ( const edm::ParameterSet iConfig)

Definition at line 103 of file LeptonUpdater.cc.

References edm::ParameterSet::getParameter(), pat::LeptonUpdater< T >::miniIsoParams_, and findQualityFiles::size.

103  {
104  miniIsoParams_[0] = iConfig.getParameter<std::vector<double>>("miniIsoParamsB");
105  miniIsoParams_[1] = iConfig.getParameter<std::vector<double>>("miniIsoParamsE");
106  if (miniIsoParams_[0].size() != 9)
107  throw cms::Exception("ParameterError", "miniIsoParamsB must have exactly 9 elements.\n");
108  if (miniIsoParams_[1].size() != 9)
109  throw cms::Exception("ParameterError", "miniIsoParamsE must have exactly 9 elements.\n");
110  }
size
Write out results.
T getParameter(std::string const &) const
std::vector< double > miniIsoParams_[2]
template<typename T >
void pat::LeptonUpdater< T >::recomputeMuonBasicSelectors ( T lep,
const reco::Vertex pv,
const bool  do_hip_mitigation_2016 
) const

Definition at line 117 of file LeptonUpdater.cc.

Referenced by pat::LeptonUpdater< T >::miniIsoParams(), and pat::LeptonUpdater< T >::produce().

119  {}
template<>
void pat::LeptonUpdater< pat::Muon >::recomputeMuonBasicSelectors ( pat::Muon lep,
const reco::Vertex pv,
const bool  do_hip_mitigation_2016 
) const

Definition at line 122 of file LeptonUpdater.cc.

References muon::makeSelectorBitset(), and reco::Muon::setSelectors().

124  {
125  lep.setSelectors(muon::makeSelectorBitset(lep, &pv, do_hip_mitigation_2016));
126  }
void setSelectors(uint64_t selectors)
Definition: Muon.h:232
reco::Muon::Selector makeSelectorBitset(reco::Muon const &muon, reco::Vertex const *vertex=0, bool run2016_hip_mitigation=false)
template<typename T >
void pat::LeptonUpdater< T >::setDZ ( T lep,
const reco::Vertex pv 
) const
inline

Definition at line 65 of file LeptonUpdater.cc.

Referenced by pat::LeptonUpdater< T >::produce().

65 {}
template<>
void pat::LeptonUpdater< pat::Electron >::setDZ ( pat::Electron anElectron,
const reco::Vertex pv 
) const

Definition at line 91 of file LeptonUpdater.cc.

References pat::Electron::gsfTrack(), reco::Vertex::position(), pat::Electron::PVDZ, pat::Electron::setDB(), HiIsolationCommonParameters_cff::track, and reco::Vertex::zError().

91  {
92  auto track = anElectron.gsfTrack();
93  anElectron.setDB(track->dz(pv.position()), std::hypot(track->dzError(), pv.zError()), pat::Electron::PVDZ);
94  }
double zError() const
error on z
Definition: Vertex.h:123
const Point & position() const
position
Definition: Vertex.h:109
void setDB(double dB, double edB, IPTYPE type)
Set impact parameter of a certain type and its uncertainty.
reco::GsfTrackRef gsfTrack() const override
override the reco::GsfElectron::gsfTrack method, to access the internal storage of the supercluster ...
template<>
void pat::LeptonUpdater< pat::Muon >::setDZ ( pat::Muon aMuon,
const reco::Vertex pv 
) const

Definition at line 97 of file LeptonUpdater.cc.

References pat::Muon::muonBestTrack(), reco::Vertex::position(), pat::Muon::PVDZ, pat::Muon::setDB(), HiIsolationCommonParameters_cff::track, and reco::Vertex::zError().

97  {
98  auto track = aMuon.muonBestTrack();
99  aMuon.setDB(track->dz(pv.position()), std::hypot(track->dzError(), pv.zError()), pat::Muon::PVDZ);
100  }
double zError() const
error on z
Definition: Vertex.h:123
const Point & position() const
position
Definition: Vertex.h:109
void setDB(double dB, double edB, IPTYPE type=PV2D)
Definition: Muon.h:247
reco::TrackRef muonBestTrack() const override
Track selected to be the best measurement of the muon parameters (including PFlow global information)...

Member Data Documentation

template<typename T >
edm::EDGetTokenT<reco::BeamSpot> pat::LeptonUpdater< T >::beamLineToken_
private

Definition at line 80 of file LeptonUpdater.cc.

Referenced by pat::LeptonUpdater< T >::produce().

template<typename T >
bool pat::LeptonUpdater< T >::computeMiniIso_
private
template<typename T >
bool pat::LeptonUpdater< T >::fixDxySign_
private

Definition at line 82 of file LeptonUpdater.cc.

Referenced by pat::LeptonUpdater< T >::produce().

template<typename T >
std::vector<double> pat::LeptonUpdater< T >::miniIsoParams_[2]
private
template<typename T >
edm::EDGetTokenT<pat::PackedCandidateCollection> pat::LeptonUpdater< T >::pcToken_
private
template<typename T >
bool pat::LeptonUpdater< T >::recomputeMuonBasicSelectors_
private
template<typename T >
edm::EDGetTokenT<std::vector<T> > pat::LeptonUpdater< T >::src_
private

Definition at line 78 of file LeptonUpdater.cc.

Referenced by pat::LeptonUpdater< T >::produce().

template<typename T >
edm::EDGetTokenT<std::vector<reco::Vertex> > pat::LeptonUpdater< T >::vertices_
private

Definition at line 79 of file LeptonUpdater.cc.

Referenced by pat::LeptonUpdater< T >::produce().