CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
EgammaHLTGsfTrackVarProducer Class Reference
Inheritance diagram for EgammaHLTGsfTrackVarProducer:
edm::global::EDProducer<> edm::global::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 EgammaHLTGsfTrackVarProducer (const edm::ParameterSet &)
 
void produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override
 
- Public Member Functions inherited from edm::global::EDProducer<>
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
EDProducer & operator= (const EDProducer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () 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
 
std::vector< bool > const & recordProvenanceList () 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)
 
TypeLabelList const & typeLabelList () const
 used by the fwk to register the list of products of this module More...
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () 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
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::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

const edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > chi2MapPutToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dEtaMapPutToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dEtaSeedMapPutToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dPhiMapPutToken_
 
const edm::EDGetTokenT< reco::ElectronCollection > electronToken_
 
const edm::EDGetTokenT< reco::GsfTrackCollection > gsfTrackToken_
 
const int lowerTrackNrToRemoveCut_
 
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magneticFieldToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > missingHitsMapPutToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > nLayerITMapPutToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > oneOverESeedMinusOneOverPMapPutToken_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > oneOverESuperMinusOneOverPMapPutToken_
 
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandToken_
 
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > trackerGeometryToken_
 
const int upperTrackNrToRemoveCut_
 
const bool useDefaultValuesForBarrel_
 
const bool useDefaultValuesForEndcap_
 
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > validHitsMapPutToken_
 

Additional Inherited Members

- Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
template<typename T >
using BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >
 
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::ProducerBase
template<Transition Tr = Transition::Event>
auto produces (std::string instanceName) noexcept
 declare what type of product will make and with which optional label More...
 
template<Transition B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<BranchType B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
template<Transition Tr = Transition::Event>
auto produces () noexcept
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename 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<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Author
Roberto Covarelli (CERN)
Id
EgammaHLTGsfTrackVarProducer.cc,v 1.1 2012/01/23 12:56:38 sharper Exp

Definition at line 41 of file EgammaHLTGsfTrackVarProducer.cc.

Constructor & Destructor Documentation

◆ EgammaHLTGsfTrackVarProducer()

EgammaHLTGsfTrackVarProducer::EgammaHLTGsfTrackVarProducer ( const edm::ParameterSet &  config)
explicit

Definition at line 72 of file EgammaHLTGsfTrackVarProducer.cc.

73  : recoEcalCandToken_(
74  consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"))),
75  electronToken_{consumes<reco::ElectronCollection>(config.getParameter<edm::InputTag>("inputCollection"))},
76  gsfTrackToken_{consumes<reco::GsfTrackCollection>(config.getParameter<edm::InputTag>("inputCollection"))},
77  beamSpotToken_{consumes<reco::BeamSpot>(config.getParameter<edm::InputTag>("beamSpotProducer"))},
78  magneticFieldToken_{esConsumes()},
79  trackerGeometryToken_{esConsumes()},
80  upperTrackNrToRemoveCut_{config.getParameter<int>("upperTrackNrToRemoveCut")},
81  lowerTrackNrToRemoveCut_{config.getParameter<int>("lowerTrackNrToRemoveCut")},
82  useDefaultValuesForBarrel_{config.getParameter<bool>("useDefaultValuesForBarrel")},
83  useDefaultValuesForEndcap_{config.getParameter<bool>("useDefaultValuesForEndcap")},
84  dEtaMapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("Deta").setBranchAlias("deta")},
85  dEtaSeedMapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("DetaSeed").setBranchAlias("detaseed")},
86  dPhiMapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("Dphi").setBranchAlias("dphi")},
87  oneOverESuperMinusOneOverPMapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("OneOESuperMinusOneOP")},
88  oneOverESeedMinusOneOverPMapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("OneOESeedMinusOneOP")},
89  missingHitsMapPutToken_{
90  produces<reco::RecoEcalCandidateIsolationMap>("MissingHits").setBranchAlias("missinghits")},
91  validHitsMapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("ValidHits").setBranchAlias("validhits")},
92  nLayerITMapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("NLayerIT").setBranchAlias("nlayerit")},
93  chi2MapPutToken_{produces<reco::RecoEcalCandidateIsolationMap>("Chi2").setBranchAlias("chi2")} {}
EgammaHLTGsfTrackVarProducer::electronToken_
const edm::EDGetTokenT< reco::ElectronCollection > electronToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:49
EgammaHLTGsfTrackVarProducer::recoEcalCandToken_
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:48
config
Definition: config.py:1
EgammaHLTGsfTrackVarProducer::oneOverESuperMinusOneOverPMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > oneOverESuperMinusOneOverPMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:64
EgammaHLTGsfTrackVarProducer::trackerGeometryToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > trackerGeometryToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:54
EgammaHLTGsfTrackVarProducer::chi2MapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > chi2MapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:69
EgammaHLTGsfTrackVarProducer::beamSpotToken_
const edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:51
EgammaHLTGsfTrackVarProducer::validHitsMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > validHitsMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:67
EgammaHLTGsfTrackVarProducer::dEtaMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dEtaMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:61
EgammaHLTGsfTrackVarProducer::magneticFieldToken_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magneticFieldToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:53
EgammaHLTGsfTrackVarProducer::oneOverESeedMinusOneOverPMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > oneOverESeedMinusOneOverPMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:65
EgammaHLTGsfTrackVarProducer::dPhiMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dPhiMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:63
EgammaHLTGsfTrackVarProducer::nLayerITMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > nLayerITMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:68
EgammaHLTGsfTrackVarProducer::dEtaSeedMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dEtaSeedMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:62
EgammaHLTGsfTrackVarProducer::gsfTrackToken_
const edm::EDGetTokenT< reco::GsfTrackCollection > gsfTrackToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:50
EgammaHLTGsfTrackVarProducer::missingHitsMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > missingHitsMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:66

Member Function Documentation

◆ fillDescriptions()

void EgammaHLTGsfTrackVarProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 95 of file EgammaHLTGsfTrackVarProducer.cc.

References edm::ConfigurationDescriptions::add(), and submitPVResolutionJobs::desc.

95  {
96  edm::ParameterSetDescription desc;
97  desc.add<edm::InputTag>(("recoEcalCandidateProducer"), edm::InputTag("hltRecoEcalSuperClusterActivityCandidate"));
98  desc.add<edm::InputTag>(("inputCollection"), edm::InputTag("hltActivityElectronGsfTracks"));
99  desc.add<edm::InputTag>(("beamSpotProducer"), edm::InputTag("hltOnlineBeamSpot"));
100  desc.add<int>(("upperTrackNrToRemoveCut"), 9999);
101  desc.add<int>(("lowerTrackNrToRemoveCut"), -1);
102  desc.add<bool>(("useDefaultValuesForBarrel"), false);
103  desc.add<bool>(("useDefaultValuesForEndcap"), false);
104 
105  descriptions.add("hltEgammaHLTGsfTrackVarProducer", desc);
106 }
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57

◆ produce()

void EgammaHLTGsfTrackVarProducer::produce ( edm::StreamID  ,
edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
) const
overridevirtual

Implements edm::global::EDProducerBase.

Definition at line 107 of file EgammaHLTGsfTrackVarProducer.cc.

References funct::abs(), AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, beamSpotToken_, chi2MapPutToken_, dEtaMapPutToken_, dEtaSeedMapPutToken_, dPhiMapPutToken_, electronToken_, edm::EventSetup::getData(), ecalDrivenGsfElectronCoresHGC_cff::gsfTracks, gsfTrackToken_, iEvent, MultiTrajectoryStateTransform::innerStateOnSurface(), edm::AssociationMap< Tag >::insert(), edm::Ref< C, T, F >::isNonnull(), HLT_2022v15_cff::magneticField, magneticFieldToken_, missingHitsMapPutToken_, nLayerITMapPutToken_, TrackCollections2monitor_cff::normalizedChi2, TrackCollections2monitor_cff::numberOfValidHits, oneOverESeedMinusOneOverPMapPutToken_, oneOverESuperMinusOneOverPMapPutToken_, multiTrajectoryStateMode::positionFromModeCartesian(), recoEcalCandToken_, trackerGeometryToken_, upperTrackNrToRemoveCut_, useDefaultValuesForBarrel_, useDefaultValuesForEndcap_, and validHitsMapPutToken_.

107  {
108  // Get the HLT filtered objects
109  auto recoEcalCandHandle = iEvent.getHandle(recoEcalCandToken_);
110 
111  auto const& beamSpotPosition = iEvent.get(beamSpotToken_).position();
112  auto const& magneticField = iSetup.getData(magneticFieldToken_);
113  auto const& trackerGeometry = iSetup.getData(trackerGeometryToken_);
114 
115  TransverseImpactPointExtrapolator extrapolator{
116  GsfPropagatorAdapter{AnalyticalPropagator(&magneticField, anyDirection)}};
117 
118  reco::RecoEcalCandidateIsolationMap dEtaMap(recoEcalCandHandle);
119  reco::RecoEcalCandidateIsolationMap dEtaSeedMap(recoEcalCandHandle);
120  reco::RecoEcalCandidateIsolationMap dPhiMap(recoEcalCandHandle);
121  reco::RecoEcalCandidateIsolationMap oneOverESuperMinusOneOverPMap(recoEcalCandHandle);
122  reco::RecoEcalCandidateIsolationMap oneOverESeedMinusOneOverPMap(recoEcalCandHandle);
123  reco::RecoEcalCandidateIsolationMap missingHitsMap(recoEcalCandHandle);
124  reco::RecoEcalCandidateIsolationMap validHitsMap(recoEcalCandHandle);
125  reco::RecoEcalCandidateIsolationMap nLayerITMap(recoEcalCandHandle);
126  reco::RecoEcalCandidateIsolationMap chi2Map(recoEcalCandHandle);
127 
128  for (unsigned int iRecoEcalCand = 0; iRecoEcalCand < recoEcalCandHandle->size(); ++iRecoEcalCand) {
129  reco::RecoEcalCandidateRef recoEcalCandRef(recoEcalCandHandle, iRecoEcalCand);
130 
131  const reco::SuperClusterRef scRef = recoEcalCandRef->superCluster();
132  // the idea is that we can take the tracks from properly associated
133  // electrons or just take all gsf tracks with that sc as a seed
134  std::vector<const reco::GsfTrack*> gsfTracks;
135  if (auto electronHandle = iEvent.getHandle(electronToken_)) {
136  for (auto const& ele : *electronHandle) {
137  if (ele.superCluster() == scRef) {
138  gsfTracks.push_back(ele.gsfTrack().get());
139  }
140  }
141  } else {
142  for (auto const& trk : iEvent.get(gsfTrackToken_)) {
143  auto elseed = trk.extra()->seedRef().castTo<reco::ElectronSeedRef>();
144  if (elseed->caloCluster().castTo<reco::SuperClusterRef>() == scRef) {
145  gsfTracks.push_back(&trk);
146  }
147  }
148  }
149 
150  int nLayerITValue = -1;
151  int validHitsValue = 0;
152  float chi2Value = 9999999.;
153  float missingHitsValue = 9999999;
154  float dEtaInValue = 999999;
155  float dEtaSeedInValue = 999999;
156  float dPhiInValue = 999999;
157  float oneOverESuperMinusOneOverPValue = 999999;
158  float oneOverESeedMinusOneOverPValue = 999999;
159 
160  const int nrTracks = gsfTracks.size();
161  const bool rmCutsDueToNrTracks = nrTracks <= lowerTrackNrToRemoveCut_ || nrTracks >= upperTrackNrToRemoveCut_;
162  // to use the default values, we require at least one track...
163  const bool useDefaultValues = std::abs(recoEcalCandRef->eta()) < 1.479
164  ? useDefaultValuesForBarrel_ && nrTracks >= 1
165  : useDefaultValuesForEndcap_ && nrTracks >= 1;
166 
167  if (rmCutsDueToNrTracks || useDefaultValues) {
168  nLayerITValue = 100;
169  dEtaInValue = 0;
170  dEtaSeedInValue = 0;
171  dPhiInValue = 0;
172  missingHitsValue = 0;
173  validHitsValue = 100;
174  chi2Value = 0;
175  oneOverESuperMinusOneOverPValue = 0;
176  oneOverESeedMinusOneOverPValue = 0;
177  } else {
178  for (size_t trkNr = 0; trkNr < gsfTracks.size(); trkNr++) {
179  GlobalPoint scPos(scRef->x(), scRef->y(), scRef->z());
180 
181  GlobalPoint trackExtrapToSC;
182  {
183  auto innTSOS =
184  MultiTrajectoryStateTransform::innerStateOnSurface(*gsfTracks[trkNr], trackerGeometry, &magneticField);
185  auto posTSOS = extrapolator.extrapolate(innTSOS, scPos);
186  multiTrajectoryStateMode::positionFromModeCartesian(posTSOS, trackExtrapToSC);
187  }
188 
189  EleRelPointPair scAtVtx(scRef->position(), trackExtrapToSC, beamSpotPosition);
190 
191  float trkP = gsfTracks[trkNr]->p();
192  if (scRef->energy() != 0 && trkP != 0) {
193  if (std::abs(1 / scRef->energy() - 1 / trkP) < oneOverESuperMinusOneOverPValue) {
194  oneOverESuperMinusOneOverPValue = std::abs(1 / scRef->energy() - 1 / trkP);
195  }
196  }
197  if (scRef->seed().isNonnull() && scRef->seed()->energy() != 0 && trkP != 0) {
198  if (std::abs(1 / scRef->seed()->energy() - 1 / trkP) < oneOverESeedMinusOneOverPValue) {
199  oneOverESeedMinusOneOverPValue = std::abs(1 / scRef->seed()->energy() - 1 / trkP);
200  }
201  }
202 
203  if (gsfTracks[trkNr]->missingInnerHits() < missingHitsValue) {
204  missingHitsValue = gsfTracks[trkNr]->missingInnerHits();
205  }
206 
207  // we are saving the best value, and highest value of validHits is the
208  // best
209  if (gsfTracks[trkNr]->numberOfValidHits() > validHitsValue) {
210  validHitsValue = gsfTracks[trkNr]->numberOfValidHits();
211  }
212 
213  if (gsfTracks[trkNr]->hitPattern().pixelLayersWithMeasurement() > nLayerITValue) {
214  nLayerITValue = gsfTracks[trkNr]->hitPattern().pixelLayersWithMeasurement();
215  }
216 
217  if (gsfTracks[trkNr]->normalizedChi2() < chi2Value) {
218  chi2Value = gsfTracks[trkNr]->normalizedChi2();
219  }
220 
221  if (std::abs(scAtVtx.dEta()) < dEtaInValue) {
222  // we are allowing them to come from different tracks
223  dEtaInValue = std::abs(scAtVtx.dEta());
224  }
225 
226  if (std::abs(scAtVtx.dEta()) < dEtaSeedInValue) {
227  dEtaSeedInValue = std::abs(scAtVtx.dEta() - scRef->position().eta() + scRef->seed()->position().eta());
228  }
229 
230  if (std::abs(scAtVtx.dPhi()) < dPhiInValue) {
231  // we are allowing them to come from different tracks
232  dPhiInValue = std::abs(scAtVtx.dPhi());
233  }
234  }
235  }
236 
237  dEtaMap.insert(recoEcalCandRef, dEtaInValue);
238  dEtaSeedMap.insert(recoEcalCandRef, dEtaSeedInValue);
239  dPhiMap.insert(recoEcalCandRef, dPhiInValue);
240  oneOverESuperMinusOneOverPMap.insert(recoEcalCandRef, oneOverESuperMinusOneOverPValue);
241  oneOverESeedMinusOneOverPMap.insert(recoEcalCandRef, oneOverESeedMinusOneOverPValue);
242  missingHitsMap.insert(recoEcalCandRef, missingHitsValue);
243  validHitsMap.insert(recoEcalCandRef, validHitsValue);
244  nLayerITMap.insert(recoEcalCandRef, nLayerITValue);
245  chi2Map.insert(recoEcalCandRef, chi2Value);
246  }
247 
248  iEvent.emplace(dEtaMapPutToken_, dEtaMap);
249  iEvent.emplace(dEtaSeedMapPutToken_, dEtaSeedMap);
250  iEvent.emplace(dPhiMapPutToken_, dPhiMap);
251  iEvent.emplace(oneOverESuperMinusOneOverPMapPutToken_, oneOverESuperMinusOneOverPMap);
252  iEvent.emplace(oneOverESeedMinusOneOverPMapPutToken_, oneOverESeedMinusOneOverPMap);
253  iEvent.emplace(missingHitsMapPutToken_, missingHitsMap);
254  iEvent.emplace(validHitsMapPutToken_, validHitsMap);
255  iEvent.emplace(nLayerITMapPutToken_, nLayerITMap);
256  iEvent.emplace(chi2MapPutToken_, chi2Map);
257 }
EgammaHLTGsfTrackVarProducer::electronToken_
const edm::EDGetTokenT< reco::ElectronCollection > electronToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:49
edm::EventSetup::getData
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Definition: EventSetup.h:119
HLT_2022v15_cff.magneticField
magneticField
Definition: HLT_2022v15_cff.py:487
EgammaHLTGsfTrackVarProducer::recoEcalCandToken_
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:48
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:238
TrackCollections2monitor_cff.normalizedChi2
normalizedChi2
Definition: TrackCollections2monitor_cff.py:252
EgammaHLTGsfTrackVarProducer::oneOverESuperMinusOneOverPMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > oneOverESuperMinusOneOverPMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:64
EgammaHLTGsfTrackVarProducer::trackerGeometryToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > trackerGeometryToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:54
EgammaHLTGsfTrackVarProducer::chi2MapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > chi2MapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:69
EgammaHLTGsfTrackVarProducer::beamSpotToken_
const edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:51
iEvent
int iEvent
Definition: GenABIO.cc:224
multiTrajectoryStateMode::positionFromModeCartesian
bool positionFromModeCartesian(TrajectoryStateOnSurface const &tsos, GlobalPoint &position)
Definition: MultiTrajectoryStateMode.cc:63
MultiTrajectoryStateTransform::innerStateOnSurface
TrajectoryStateOnSurface innerStateOnSurface(const reco::GsfTrack &tk) const
Definition: MultiTrajectoryStateTransform.cc:23
EgammaHLTGsfTrackVarProducer::validHitsMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > validHitsMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:67
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
EgammaHLTGsfTrackVarProducer::dEtaMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dEtaMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:61
TrackCollections2monitor_cff.numberOfValidHits
numberOfValidHits
Definition: TrackCollections2monitor_cff.py:253
AnalyticalPropagator_cfi.AnalyticalPropagator
AnalyticalPropagator
Definition: AnalyticalPropagator_cfi.py:3
EgammaHLTGsfTrackVarProducer::magneticFieldToken_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magneticFieldToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:53
ecalDrivenGsfElectronCoresHGC_cff.gsfTracks
gsfTracks
Definition: ecalDrivenGsfElectronCoresHGC_cff.py:4
EgammaHLTGsfTrackVarProducer::oneOverESeedMinusOneOverPMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > oneOverESeedMinusOneOverPMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:65
EgammaHLTGsfTrackVarProducer::dPhiMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dPhiMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:63
EgammaHLTGsfTrackVarProducer::nLayerITMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > nLayerITMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:68
EgammaHLTGsfTrackVarProducer::dEtaSeedMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > dEtaSeedMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:62
EgammaHLTGsfTrackVarProducer::gsfTrackToken_
const edm::EDGetTokenT< reco::GsfTrackCollection > gsfTrackToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:50
EgammaHLTGsfTrackVarProducer::missingHitsMapPutToken_
const edm::EDPutTokenT< reco::RecoEcalCandidateIsolationMap > missingHitsMapPutToken_
Definition: EgammaHLTGsfTrackVarProducer.cc:66

Member Data Documentation

◆ beamSpotToken_

const edm::EDGetTokenT<reco::BeamSpot> EgammaHLTGsfTrackVarProducer::beamSpotToken_
private

Definition at line 51 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ chi2MapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::chi2MapPutToken_
private

Definition at line 69 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ dEtaMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::dEtaMapPutToken_
private

Definition at line 61 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ dEtaSeedMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::dEtaSeedMapPutToken_
private

Definition at line 62 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ dPhiMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::dPhiMapPutToken_
private

Definition at line 63 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ electronToken_

const edm::EDGetTokenT<reco::ElectronCollection> EgammaHLTGsfTrackVarProducer::electronToken_
private

Definition at line 49 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ gsfTrackToken_

const edm::EDGetTokenT<reco::GsfTrackCollection> EgammaHLTGsfTrackVarProducer::gsfTrackToken_
private

Definition at line 50 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ lowerTrackNrToRemoveCut_

const int EgammaHLTGsfTrackVarProducer::lowerTrackNrToRemoveCut_
private

Definition at line 57 of file EgammaHLTGsfTrackVarProducer.cc.

◆ magneticFieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> EgammaHLTGsfTrackVarProducer::magneticFieldToken_
private

Definition at line 53 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ missingHitsMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::missingHitsMapPutToken_
private

Definition at line 66 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ nLayerITMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::nLayerITMapPutToken_
private

Definition at line 68 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ oneOverESeedMinusOneOverPMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::oneOverESeedMinusOneOverPMapPutToken_
private

Definition at line 65 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ oneOverESuperMinusOneOverPMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::oneOverESuperMinusOneOverPMapPutToken_
private

Definition at line 64 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ recoEcalCandToken_

const edm::EDGetTokenT<reco::RecoEcalCandidateCollection> EgammaHLTGsfTrackVarProducer::recoEcalCandToken_
private

Definition at line 48 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ trackerGeometryToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> EgammaHLTGsfTrackVarProducer::trackerGeometryToken_
private

Definition at line 54 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ upperTrackNrToRemoveCut_

const int EgammaHLTGsfTrackVarProducer::upperTrackNrToRemoveCut_
private

Definition at line 56 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ useDefaultValuesForBarrel_

const bool EgammaHLTGsfTrackVarProducer::useDefaultValuesForBarrel_
private

Definition at line 58 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ useDefaultValuesForEndcap_

const bool EgammaHLTGsfTrackVarProducer::useDefaultValuesForEndcap_
private

Definition at line 59 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

◆ validHitsMapPutToken_

const edm::EDPutTokenT<reco::RecoEcalCandidateIsolationMap> EgammaHLTGsfTrackVarProducer::validHitsMapPutToken_
private

Definition at line 67 of file EgammaHLTGsfTrackVarProducer.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.14