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 19 of file LeptonUpdater.cc.

Constructor & Destructor Documentation

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

Definition at line 23 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_.

23  :
24  src_(consumes<std::vector<T>>(iConfig.getParameter<edm::InputTag>("src"))),
25  vertices_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("vertices"))),
26  beamLineToken_(consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamspot"))),
27  computeMiniIso_(iConfig.getParameter<bool>("computeMiniIso")),
28  fixDxySign_(iConfig.getParameter<bool>("fixDxySign"))
29  {
30  //for mini-isolation calculation
31  if (computeMiniIso_) {
32  readMiniIsoParams(iConfig);
33  pcToken_ = consumes<pat::PackedCandidateCollection >(iConfig.getParameter<edm::InputTag>("pfCandsForMiniIso"));
34  }
36  if (typeid(T) == typeid(pat::Muon)) recomputeMuonBasicSelectors_ = iConfig.getParameter<bool>("recomputeMuonBasicSelectors");
37  produces<std::vector<T>>();
38  }
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 40 of file LeptonUpdater.cc.

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

40 {}

Member Function Documentation

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

Definition at line 44 of file LeptonUpdater.cc.

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

44  {
46  desc.add<edm::InputTag>("src")->setComment("Lepton collection");
47  desc.add<edm::InputTag>("vertices")->setComment("Vertex collection");
48  desc.add<edm::InputTag>("beamspot", edm::InputTag("offlineBeamSpot"))->setComment("Beam spot");
49  desc.add<bool>("computeMiniIso", false)->setComment("Recompute miniIsolation");
50  desc.add<bool>("fixDxySign", false)->setComment("Fix the IP sign");
51  desc.addOptional<edm::InputTag>("pfCandsForMiniIso", edm::InputTag("packedPFCandidates"))->setComment("PackedCandidate collection used for miniIso");
52  if (typeid(T) == typeid(pat::Muon)) {
53  desc.add<bool>("recomputeMuonBasicSelectors",false)->setComment("Recompute basic cut-based muon selector flags");
54  desc.addOptional<std::vector<double>>("miniIsoParams")->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")->setComment("Parameters used for miniIso in the barrel (as in PATElectronProducer)");
58  desc.addOptional<std::vector<double>>("miniIsoParamsE")->setComment("Parameters used for miniIso in the endcap (as in PATElectronProducer)");
59  descriptions.add("electronsUpdated", desc);
60  }
61  }
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 107 of file LeptonUpdater.cc.

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

107  {
108  return miniIsoParams_[lep.isEE()];
109  }
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 122 of file LeptonUpdater.cc.

References funct::abs(), pat::LeptonUpdater< T >::beamLineToken_, ecalDrivenElectronSeedsParameters_cff::beamSpot, pat::LeptonUpdater< T >::computeMiniIso_, 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().

122  {
124  iEvent.getByToken(src_, src);
125 
127  iEvent.getByToken(vertices_, vertices);
128  const reco::Vertex & pv = vertices->front();
129 
131  if(computeMiniIso_) iEvent.getByToken(pcToken_, pc);
132 
133  edm::Handle<reco::BeamSpot> beamSpotHandle;
134  iEvent.getByToken(beamLineToken_, beamSpotHandle);
136  bool beamSpotIsValid = false;
137  if (beamSpotHandle.isValid()) {
138  beamSpot = *beamSpotHandle;
139  beamSpotIsValid = true;
140  } else {
141  edm::LogError("DataNotAvailable") << "No beam spot available \n";
142  }
143 
144  std::unique_ptr<std::vector<T>> out(new std::vector<T>(*src));
145 
146  const bool do_hip_mitigation_2016 = recomputeMuonBasicSelectors_ && (272728 <= iEvent.run() && iEvent.run() <= 278808);
147 
148  for (unsigned int i = 0, n = src->size(); i < n; ++i) {
149  T & lep = (*out)[i];
150  setDZ(lep, pv);
151  if (computeMiniIso_) {
152  const auto & params = miniIsoParams(lep);
153  pat::PFIsolation miniiso = pat::getMiniPFIsolation(pc.product(), lep.polarP4(),
154  params[0], params[1], params[2],
155  params[3], params[4], params[5],
156  params[6], params[7], params[8]);
157  lep.setMiniPFIsolation(miniiso);
158  }
159  if (recomputeMuonBasicSelectors_) recomputeMuonBasicSelectors(lep,pv,do_hip_mitigation_2016);
160  //Fixing the sign of impact parameters
161  if (fixDxySign_) {
162  float signPV = 1.;
163  float signBS = 1.;
164  if (beamSpotIsValid) {
165  signBS = copysign(1., lep.bestTrack()->dxy(beamSpot));
166  }
167  signPV = copysign(1., lep.bestTrack()->dxy(pv.position()));
168  lep.setDB(abs(lep.dB(T::PV2D)) * signPV, lep.edB(T::PV2D), T::PV2D);
169  lep.setDB(abs(lep.dB(T::BS2D)) * signBS, lep.edB(T::BS2D), T::BS2D);
170  }
171  }
172 
173  iEvent.put(std::move(out));
174 }
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:125
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
const Point & position() const
position
Definition: Vertex.h:109
edm::EDGetTokenT< std::vector< T > > src_
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
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
template<typename T >
void pat::LeptonUpdater< T >::readMiniIsoParams ( const edm::ParameterSet iConfig)
inline

Definition at line 65 of file LeptonUpdater.cc.

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

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

65  {
66  miniIsoParams_[0] = iConfig.getParameter<std::vector<double> >("miniIsoParams");
67  if(miniIsoParams_[0].size() != 9) throw cms::Exception("ParameterError", "miniIsoParams must have exactly 9 elements.\n");
68  }
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 100 of file LeptonUpdater.cc.

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

100  {
101  miniIsoParams_[0] = iConfig.getParameter<std::vector<double> >("miniIsoParamsB");
102  miniIsoParams_[1] = iConfig.getParameter<std::vector<double> >("miniIsoParamsE");
103  if(miniIsoParams_[0].size() != 9) throw cms::Exception("ParameterError", "miniIsoParamsB must have exactly 9 elements.\n");
104  if(miniIsoParams_[1].size() != 9) throw cms::Exception("ParameterError", "miniIsoParamsE must have exactly 9 elements.\n");
105  }
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
template<>
void pat::LeptonUpdater< pat::Muon >::recomputeMuonBasicSelectors ( pat::Muon lep,
const reco::Vertex pv,
const bool  do_hip_mitigation_2016 
) const

Definition at line 115 of file LeptonUpdater.cc.

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

115  {
116  lep.setSelectors(muon::makeSelectorBitset(lep, &pv, do_hip_mitigation_2016));
117  }
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 63 of file LeptonUpdater.cc.

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

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

Definition at line 88 of file LeptonUpdater.cc.

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

88  {
89  auto track = anElectron.gsfTrack();
90  anElectron.setDB( track->dz(pv.position()), std::hypot(track->dzError(), pv.zError()), pat::Electron::PVDZ );
91  }
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 94 of file LeptonUpdater.cc.

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

94  {
95  auto track = aMuon.muonBestTrack();
96  aMuon.setDB( track->dz(pv.position()), std::hypot(track->dzError(), pv.zError()), pat::Muon::PVDZ );
97  }
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:245
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 77 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 79 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 75 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 76 of file LeptonUpdater.cc.

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