CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
EgammaHLTPFChargedIsolationProducer Class Reference

#include <EgammaHLTPFChargedIsolationProducer.h>

Inheritance diagram for EgammaHLTPFChargedIsolationProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 EgammaHLTPFChargedIsolationProducer (const edm::ParameterSet &)
 
void produce (edm::Event &, const edm::EventSetup &) override
 
 ~EgammaHLTPFChargedIsolationProducer ()
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducer ()
 
- 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 ()
 
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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
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
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

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

Private Attributes

const edm::EDGetTokenT< reco::BeamSpotbeamSpotProducer_
 
const double drMax_
 
const double drVetoBarrel_
 
const double drVetoEndcap_
 
const double dxyMax_
 
const double dzMax_
 
edm::EDGetTokenT< reco::ElectronCollectionelectronProducer_
 
const edm::EDGetTokenT< reco::PFCandidateCollectionpfCandidateProducer_
 
const int pfToUse_
 
const double ptMin_
 
edm::EDGetTokenT< reco::RecoEcalCandidateCollectionrecoEcalCandidateProducer_
 
const bool useGsfTrack_
 
const bool useSCRefs_
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
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 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

Definition at line 31 of file EgammaHLTPFChargedIsolationProducer.h.

Constructor & Destructor Documentation

EgammaHLTPFChargedIsolationProducer::EgammaHLTPFChargedIsolationProducer ( const edm::ParameterSet config)
explicit
Author
Matteo Sani (UCSD)

Definition at line 20 of file EgammaHLTPFChargedIsolationProducer.cc.

References electronProducer_, edm::ParameterSet::getParameter(), recoEcalCandidateProducer_, and useSCRefs_.

20  :
21  pfCandidateProducer_(consumes<reco::PFCandidateCollection>(config.getParameter<edm::InputTag>("pfCandidatesProducer"))),
22  beamSpotProducer_ (consumes<reco::BeamSpot>(config.getParameter<edm::InputTag>("beamSpotProducer"))),
23  useGsfTrack_ (config.getParameter<bool>("useGsfTrack")),
24  useSCRefs_ (config.getParameter<bool>("useSCRefs")),
25  drMax_ (config.getParameter<double>("drMax")),
26  drVetoBarrel_ (config.getParameter<double>("drVetoBarrel")),
27  drVetoEndcap_ (config.getParameter<double>("drVetoEndcap")),
28  ptMin_ (config.getParameter<double>("ptMin")),
29  dzMax_ (config.getParameter<double>("dzMax")),
30  dxyMax_ (config.getParameter<double>("dxyMax")),
31  pfToUse_ (config.getParameter<int>("pfCandidateType")) {
32 
33  if(useSCRefs_) {
34  recoEcalCandidateProducer_ = consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"));
35  produces < reco::RecoEcalCandidateIsolationMap >();
36  } else {
37  electronProducer_ = consumes<reco::ElectronCollection>(config.getParameter<edm::InputTag>("electronProducer"));
38  produces < reco::ElectronIsolationMap >();
39  }
40 }
T getParameter(std::string const &) const
const edm::EDGetTokenT< reco::PFCandidateCollection > pfCandidateProducer_
edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandidateProducer_
const edm::EDGetTokenT< reco::BeamSpot > beamSpotProducer_
edm::EDGetTokenT< reco::ElectronCollection > electronProducer_
EgammaHLTPFChargedIsolationProducer::~EgammaHLTPFChargedIsolationProducer ( )
inline

Definition at line 34 of file EgammaHLTPFChargedIsolationProducer.h.

34 {};

Member Function Documentation

void EgammaHLTPFChargedIsolationProducer::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 42 of file EgammaHLTPFChargedIsolationProducer.cc.

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

42  {
44  desc.add<edm::InputTag>("electronProducer", edm::InputTag("hltEle27WP80PixelMatchElectronsL1SeededPF"));
45  desc.add<edm::InputTag>("recoEcalCandidateProducer", edm::InputTag("hltL1SeededRecoEcalCandidatePF"));
46  desc.add<edm::InputTag>("pfCandidatesProducer", edm::InputTag("hltParticleFlowReg"));
47  desc.add<edm::InputTag>("beamSpotProducer", edm::InputTag("hltOnlineBeamSpot"));
48  desc.add<bool>("useGsfTrack", false);
49  desc.add<bool>("useSCRefs", false);
50  desc.add<double>("drMax", 0.3);
51  desc.add<double>("drVetoBarrel", 0.02);
52  desc.add<double>("drVetoEndcap", 0.02);
53  desc.add<double>("ptMin", 0.0);
54  desc.add<double>("dzMax", 0.2);
55  desc.add<double>("dxyMax", 0.1);
56  desc.add<int>("pfCandidateType", 1);
57  descriptions.add(("hltEgammaHLTPFChargedIsolationProducer"), desc);
58 }
ParameterDescriptionBase * add(U const &iLabel, T const &value)
void add(std::string const &label, ParameterSetDescription const &psetDescription)
void EgammaHLTPFChargedIsolationProducer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
override

Definition at line 60 of file EgammaHLTPFChargedIsolationProducer.cc.

References beamSpotProducer_, deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, drMax_, drVetoBarrel_, drVetoEndcap_, particleFlowDisplacedVertexCandidate_cfi::dxy, dxyMax_, allConversions_cfi::dz, dzMax_, electronProducer_, edm::Event::getByToken(), mps_fire::i, edm::AssociationMap< Tag >::insert(), reco::LeafCandidate::momentum(), reco::PFCandidate::particleId(), pfCandidateProducer_, pfToUse_, reco::BeamSpot::position(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), ptMin_, edm::Event::put(), recoEcalCandidateProducer_, reco::PFCandidate::trackRef(), and useSCRefs_.

60  {
61 
65  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
66 
67  iEvent.getByToken(pfCandidateProducer_, pfHandle);
68  const reco::PFCandidateCollection* forIsolation = pfHandle.product();
69 
70  if(useSCRefs_) {
71 
72  iEvent.getByToken(recoEcalCandidateProducer_, recoEcalCandHandle);
73  reco::RecoEcalCandidateIsolationMap recoEcalCandMap(recoEcalCandHandle);
74 
75  iEvent.getByToken(beamSpotProducer_, recoBeamSpotHandle);
76  const reco::BeamSpot::Point& beamSpotPosition = recoBeamSpotHandle->position();
77 
78  float dRveto = -1;
79 
80  for(unsigned int iReco=0; iReco<recoEcalCandHandle->size(); iReco++) {
81  reco::RecoEcalCandidateRef candRef(recoEcalCandHandle, iReco);
82 
83  if (fabs(candRef->eta())<1.479)
84  dRveto = drVetoBarrel_;
85  else
86  dRveto = drVetoEndcap_;
87 
88  // Shift the RecoEcalCandidate direction vector according to the vertex
89  math::XYZVector candDirectionWrtVtx(candRef->superCluster()->x() - beamSpotPosition.x(),
90  candRef->superCluster()->y() - beamSpotPosition.y(),
91  candRef->superCluster()->z() - beamSpotPosition.z());
92 
93  float sum = 0;
94 
95  // Loop over the PFCandidates
96  for(unsigned i=0; i<forIsolation->size(); i++) {
97  const reco::PFCandidate& pfc = (*forIsolation)[i];
98 
99  //require that the PFCandidate is a charged hadron
100  if (pfc.particleId() == pfToUse_) {
101 
102  if(pfc.pt() < ptMin_) continue;
103 
104  float dz = fabs(pfc.trackRef()->dz(beamSpotPosition));
105  if(dz > dzMax_) continue;
106 
107  float dxy = fabs(pfc.trackRef()->dxy(beamSpotPosition));
108  if(fabs(dxy) > dxyMax_) continue;
109 
110  float dR = deltaR(candDirectionWrtVtx.Eta(), candDirectionWrtVtx.Phi(), pfc.momentum().Eta(), pfc.momentum().Phi());
111  if(dR > drMax_ || dR < dRveto) continue;
112 
113  sum += pfc.pt();
114  }
115  }
116 
117  recoEcalCandMap.insert(candRef, sum);
118  }
119  iEvent.put(std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandMap));
120 
121  } else {
122 
123  iEvent.getByToken(electronProducer_,electronHandle);
124  reco::ElectronIsolationMap eleMap(electronHandle);
125 
126  float dRveto = -1;
127 
128  for(unsigned int iEl=0; iEl<electronHandle->size(); iEl++) {
129  reco::ElectronRef eleRef(electronHandle, iEl);
130  //const reco::Track* eleTrk = useGsfTrack_ ? &*eleRef->gsfTrack() : &*eleRef->track();
131 
132  if (fabs(eleRef->eta())<1.479)
133  dRveto = drVetoBarrel_;
134  else
135  dRveto = drVetoEndcap_;
136 
137  float sum = 0;
138 
139  // Loop over the PFCandidates
140  for(unsigned i=0; i<forIsolation->size(); i++) {
141  const reco::PFCandidate& pfc = (*forIsolation)[i];
142 
143  //require that the PFCandidate is a charged hadron
144  if (pfc.particleId() == pfToUse_) {
145 
146  if(pfc.pt() < ptMin_) continue;
147 
148  float dz = fabs(pfc.trackRef()->dz(eleRef->vertex()));
149  if(dz > dzMax_) continue;
150 
151  float dxy = fabs(pfc.trackRef()->dxy(eleRef->vertex()));
152  if(fabs(dxy) > dxyMax_) continue;
153 
154  float dR = deltaR(eleRef->eta(), eleRef->phi(), pfc.momentum().Eta(), pfc.momentum().Phi());
155  if(dR > drMax_ || dR < dRveto) continue;
156 
157  sum += pfc.pt();
158  }
159  }
160 
161  eleMap.insert(eleRef, sum);
162  }
163  iEvent.put(std::make_unique<reco::ElectronIsolationMap>(eleMap));
164  }
165 }
virtual double pt() const final
transverse momentum
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
const edm::EDGetTokenT< reco::PFCandidateCollection > pfCandidateProducer_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
virtual Vector momentum() const final
spatial momentum vector
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:29
reco::TrackRef trackRef() const
Definition: PFCandidate.cc:438
edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandidateProducer_
const edm::EDGetTokenT< reco::BeamSpot > beamSpotProducer_
edm::EDGetTokenT< reco::ElectronCollection > electronProducer_
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
std::vector< reco::PFCandidate > PFCandidateCollection
collection of PFCandidates
T const * product() const
Definition: Handle.h:81
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:39
const Point & position() const
position
Definition: BeamSpot.h:62
virtual ParticleType particleId() const
Definition: PFCandidate.h:373

Member Data Documentation

const edm::EDGetTokenT<reco::BeamSpot> EgammaHLTPFChargedIsolationProducer::beamSpotProducer_
private

Definition at line 44 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

const double EgammaHLTPFChargedIsolationProducer::drMax_
private

Definition at line 49 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

const double EgammaHLTPFChargedIsolationProducer::drVetoBarrel_
private

Definition at line 50 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

const double EgammaHLTPFChargedIsolationProducer::drVetoEndcap_
private

Definition at line 51 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

const double EgammaHLTPFChargedIsolationProducer::dxyMax_
private

Definition at line 54 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

const double EgammaHLTPFChargedIsolationProducer::dzMax_
private

Definition at line 53 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

edm::EDGetTokenT<reco::ElectronCollection> EgammaHLTPFChargedIsolationProducer::electronProducer_
private
const edm::EDGetTokenT<reco::PFCandidateCollection> EgammaHLTPFChargedIsolationProducer::pfCandidateProducer_
private

Definition at line 43 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

const int EgammaHLTPFChargedIsolationProducer::pfToUse_
private

Definition at line 55 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

const double EgammaHLTPFChargedIsolationProducer::ptMin_
private

Definition at line 52 of file EgammaHLTPFChargedIsolationProducer.h.

Referenced by produce().

edm::EDGetTokenT<reco::RecoEcalCandidateCollection> EgammaHLTPFChargedIsolationProducer::recoEcalCandidateProducer_
private
const bool EgammaHLTPFChargedIsolationProducer::useGsfTrack_
private

Definition at line 46 of file EgammaHLTPFChargedIsolationProducer.h.

const bool EgammaHLTPFChargedIsolationProducer::useSCRefs_
private