CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
TracksterToSimTracksterAssociatorProducer Class Reference

#include <TracksterToSimTracksterAssociatorProducer.h>

Inheritance diagram for TracksterToSimTracksterAssociatorProducer:
edm::global::EDProducer<> edm::global::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 TracksterToSimTracksterAssociatorProducer (const edm::ParameterSet &)
 
 ~TracksterToSimTracksterAssociatorProducer () 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 noexcept final
 
bool wantsGlobalRuns () const noexcept final
 
bool wantsInputProcessBlocks () const noexcept final
 
bool wantsProcessBlocks () const noexcept final
 
bool wantsStreamLuminosityBlocks () const noexcept final
 
bool wantsStreamRuns () const noexcept 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
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > 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
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices 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 Member Functions

void produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override
 

Private Attributes

edm::EDGetTokenT< std::vector< reco::CaloCluster > > layerClustersCollectionToken_
 
edm::EDGetTokenT< ticl::AssociationMap< ticl::mapWithFraction, std::vector< reco::CaloCluster >, std::vector< ticl::Trackster > > > LayerClusterToSimTracksterMapToken_
 
edm::EDGetTokenT< ticl::AssociationMap< ticl::mapWithFraction, std::vector< reco::CaloCluster >, std::vector< ticl::Trackster > > > LayerClusterToTracksterMapToken_
 
edm::EDGetTokenT< std::vector< ticl::Trackster > > recoTracksterCollectionToken_
 
edm::EDGetTokenT< std::vector< ticl::Trackster > > simTracksterCollectionToken_
 

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 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

Definition at line 23 of file TracksterToSimTracksterAssociatorProducer.h.

Constructor & Destructor Documentation

◆ TracksterToSimTracksterAssociatorProducer()

TracksterToSimTracksterAssociatorProducer::TracksterToSimTracksterAssociatorProducer ( const edm::ParameterSet &  pset)
explicit

Definition at line 6 of file TracksterToSimTracksterAssociatorProducer.cc.

References edm::ProductRegistryHelper::produces().

7  : recoTracksterCollectionToken_(
8  consumes<std::vector<ticl::Trackster>>(pset.getParameter<edm::InputTag>("tracksters"))),
9  simTracksterCollectionToken_(
10  consumes<std::vector<ticl::Trackster>>(pset.getParameter<edm::InputTag>("simTracksters"))),
11  layerClustersCollectionToken_(
12  consumes<std::vector<reco::CaloCluster>>(pset.getParameter<edm::InputTag>("layerClusters"))),
13  LayerClusterToTracksterMapToken_(
14  consumes<
15  ticl::AssociationMap<ticl::mapWithFraction, std::vector<reco::CaloCluster>, std::vector<ticl::Trackster>>>(
16  pset.getParameter<edm::InputTag>("tracksterMap"))),
17  LayerClusterToSimTracksterMapToken_(
18  consumes<
19  ticl::AssociationMap<ticl::mapWithFraction, std::vector<reco::CaloCluster>, std::vector<ticl::Trackster>>>(
20  pset.getParameter<edm::InputTag>("simTracksterMap"))) {
21  produces<
22  ticl::AssociationMap<ticl::mapWithFractionAndScore, std::vector<ticl::Trackster>, std::vector<ticl::Trackster>>>(
23  "tracksterToSimTracksterMap");
24  produces<
25  ticl::AssociationMap<ticl::mapWithFractionAndScore, std::vector<ticl::Trackster>, std::vector<ticl::Trackster>>>(
26  "simTracksterToTracksterMap");
27 }
TracksterToSimTracksterAssociatorProducer::layerClustersCollectionToken_
edm::EDGetTokenT< std::vector< reco::CaloCluster > > layerClustersCollectionToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:35
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:160
TracksterToSimTracksterAssociatorProducer::simTracksterCollectionToken_
edm::EDGetTokenT< std::vector< ticl::Trackster > > simTracksterCollectionToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:34
ticl::mapWithFraction
std::vector< std::vector< std::pair< unsigned int, float > >> mapWithFraction
Definition: TICLAssociationMap.h:19
TracksterToSimTracksterAssociatorProducer::LayerClusterToTracksterMapToken_
edm::EDGetTokenT< ticl::AssociationMap< ticl::mapWithFraction, std::vector< reco::CaloCluster >, std::vector< ticl::Trackster > > > LayerClusterToTracksterMapToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:38
TracksterToSimTracksterAssociatorProducer::LayerClusterToSimTracksterMapToken_
edm::EDGetTokenT< ticl::AssociationMap< ticl::mapWithFraction, std::vector< reco::CaloCluster >, std::vector< ticl::Trackster > > > LayerClusterToSimTracksterMapToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:41
TracksterToSimTracksterAssociatorProducer::recoTracksterCollectionToken_
edm::EDGetTokenT< std::vector< ticl::Trackster > > recoTracksterCollectionToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:33

◆ ~TracksterToSimTracksterAssociatorProducer()

TracksterToSimTracksterAssociatorProducer::~TracksterToSimTracksterAssociatorProducer ( )
override

Definition at line 29 of file TracksterToSimTracksterAssociatorProducer.cc.

29 {}

Member Function Documentation

◆ fillDescriptions()

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

Definition at line 192 of file TracksterToSimTracksterAssociatorProducer.cc.

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

192  {
193  edm::ParameterSetDescription desc;
194  desc.add<edm::InputTag>("tracksters", edm::InputTag("trackstersProducer"));
195  desc.add<edm::InputTag>("simTracksters", edm::InputTag("simTrackstersProducer"));
196  desc.add<edm::InputTag>("layerClusters", edm::InputTag("hgcalMergeLayerClusters"));
197  desc.add<edm::InputTag>("tracksterMap", edm::InputTag("tracksterAssociatorProducer"));
198  desc.add<edm::InputTag>("simTracksterMap", edm::InputTag("simTracksterAssociatorProducer"));
199  descriptions.add("tracksterToSimTracksterAssociatorProducer", desc);
200 }
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57

◆ produce()

void TracksterToSimTracksterAssociatorProducer::produce ( edm::StreamID  ,
edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
) const
overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 31 of file TracksterToSimTracksterAssociatorProducer.cc.

References mps_fire::i, iEvent, hltEgammaHGCALIDVarsL1Seeded_cfi::layerClusters, layerClustersCollectionToken_, LayerClusterToSimTracksterMapToken_, LayerClusterToTracksterMapToken_, SiStripPI::min, eostools::move(), recoTracksterCollectionToken_, offlineSlimmedPrimaryVertices_cfi::score, egamma::sharedEnergy(), simTracksterCollectionToken_, jetUpdater_cfi::sort, and tier0::unique().

33  {
34  edm::Handle<std::vector<ticl::Trackster>> recoTrackstersHandle;
35  iEvent.getByToken(recoTracksterCollectionToken_, recoTrackstersHandle);
36  const auto& recoTracksters = *recoTrackstersHandle;
37 
38  edm::Handle<std::vector<ticl::Trackster>> simTrackstersHandle;
39  iEvent.getByToken(simTracksterCollectionToken_, simTrackstersHandle);
40  const auto& simTracksters = *simTrackstersHandle;
41 
42  edm::Handle<std::vector<reco::CaloCluster>> layerClustersHandle;
43  iEvent.getByToken(layerClustersCollectionToken_, layerClustersHandle);
44  const auto& layerClusters = *layerClustersHandle;
45 
46  edm::Handle<ticl::AssociationMap<ticl::mapWithFraction, std::vector<reco::CaloCluster>, std::vector<ticl::Trackster>>>
47  layerClusterToTracksterMapHandle;
48  iEvent.getByToken(LayerClusterToTracksterMapToken_, layerClusterToTracksterMapHandle);
49  const auto& layerClusterToTracksterMap = *layerClusterToTracksterMapHandle;
50 
51  edm::Handle<ticl::AssociationMap<ticl::mapWithFraction, std::vector<reco::CaloCluster>, std::vector<ticl::Trackster>>>
52  layerClusterToSimTracksterMapHandle;
53  iEvent.getByToken(LayerClusterToSimTracksterMapToken_, layerClusterToSimTracksterMapHandle);
54  const auto& layerClusterToSimTracksterMap = *layerClusterToSimTracksterMapHandle;
55 
56  auto tracksterToSimTracksterMap = std::make_unique<
57  ticl::AssociationMap<ticl::mapWithFractionAndScore, std::vector<ticl::Trackster>, std::vector<ticl::Trackster>>>(
58  recoTrackstersHandle, simTrackstersHandle, iEvent);
59  auto simTracksterToTracksterMap = std::make_unique<
60  ticl::AssociationMap<ticl::mapWithFractionAndScore, std::vector<ticl::Trackster>, std::vector<ticl::Trackster>>>(
61  simTrackstersHandle, recoTrackstersHandle, iEvent);
62 
63  for (unsigned int tracksterIndex = 0; tracksterIndex < recoTracksters.size(); ++tracksterIndex) {
64  const auto& recoTrackster = recoTracksters[tracksterIndex];
65  edm::Ref<std::vector<ticl::Trackster>> recoTracksterRef(recoTrackstersHandle, tracksterIndex);
66  const auto& layerClustersIds = recoTrackster.vertices();
67  float recoToSimScoresDenominator = 0.f;
68  ticl::mapWithFraction layerClusterToAssociatedSimTracksterMap(layerClustersIds.size());
69  std::vector<unsigned int> associatedSimTracksterIndices;
70  for (unsigned int i = 0; i < layerClustersIds.size(); ++i) {
71  unsigned int layerClusterId = layerClustersIds[i];
72  const auto& layerCluster = layerClusters[layerClusterId];
73  float recoFraction = 1.f / recoTrackster.vertex_multiplicity(i);
74  float squaredRecoFraction = recoFraction * recoFraction;
75  float squaredLayerClusterEnergy = layerCluster.energy() * layerCluster.energy();
76  recoToSimScoresDenominator += squaredLayerClusterEnergy * squaredRecoFraction;
77  const auto& simTracksterVec = layerClusterToSimTracksterMap[layerClusterId];
78  for (const auto& [simTracksterIndex, simSharedEnergy] : simTracksterVec) {
79  layerClusterToAssociatedSimTracksterMap[i].push_back({simTracksterIndex, simSharedEnergy});
80  associatedSimTracksterIndices.push_back(simTracksterIndex);
81  }
82  }
83 
84  // Keep only unique associatedSimTracksterIndices
85  std::sort(associatedSimTracksterIndices.begin(), associatedSimTracksterIndices.end());
86  associatedSimTracksterIndices.erase(
87  std::unique(associatedSimTracksterIndices.begin(), associatedSimTracksterIndices.end()),
88  associatedSimTracksterIndices.end());
89 
90  // Add missing sim tracksters with 0 shared energy to layerClusterToAssociatedSimTracksterMap
91  for (unsigned int i = 0; i < layerClustersIds.size(); ++i) {
92  unsigned int layerClusterId = layerClustersIds[i];
93  const auto& simTracksterVec = layerClusterToSimTracksterMap[layerClusterId];
94  for (unsigned int simTracksterIndex : associatedSimTracksterIndices) {
95  if (std::find_if(simTracksterVec.begin(), simTracksterVec.end(), [simTracksterIndex](const auto& pair) {
96  return pair.first == simTracksterIndex;
97  }) == simTracksterVec.end()) {
98  layerClusterToAssociatedSimTracksterMap[i].push_back({simTracksterIndex, 0.f});
99  }
100  }
101  }
102 
103  const float invDenominator = 1.f / recoToSimScoresDenominator;
104 
105  for (unsigned int i = 0; i < layerClustersIds.size(); ++i) {
106  unsigned int layerClusterId = layerClustersIds[i];
107  const auto& layerCluster = layerClusters[layerClusterId];
108  float recoFraction = 1.f / recoTrackster.vertex_multiplicity(i);
109  float squaredRecoFraction = recoFraction * recoFraction;
110  float squaredLayerClusterEnergy = layerCluster.energy() * layerCluster.energy();
111  float recoSharedEnergy = layerCluster.energy() * recoFraction;
112  float invLayerClusterEnergy = 1.f / layerCluster.energy();
113  const auto& simTracksterVec = layerClusterToAssociatedSimTracksterMap[i];
114  for (const auto& [simTracksterIndex, simSharedEnergy] : simTracksterVec) {
115  edm::Ref<std::vector<ticl::Trackster>> simTracksterRef(simTrackstersHandle, simTracksterIndex);
116  float sharedEnergy = std::min(simSharedEnergy, recoSharedEnergy);
117  float simFraction = simSharedEnergy * invLayerClusterEnergy;
118  float score = invDenominator *
119  std::min(squaredRecoFraction, (recoFraction - simFraction) * (recoFraction - simFraction)) *
120  squaredLayerClusterEnergy;
121  tracksterToSimTracksterMap->insert(recoTracksterRef, simTracksterRef, sharedEnergy, score);
122  }
123  }
124  }
125 
126  for (unsigned int tracksterIndex = 0; tracksterIndex < simTracksters.size(); ++tracksterIndex) {
127  const auto& simTrackster = simTracksters[tracksterIndex];
128  edm::Ref<std::vector<ticl::Trackster>> simTracksterRef(simTrackstersHandle, tracksterIndex);
129  const auto& layerClustersIds = simTrackster.vertices();
130  float simToRecoScoresDenominator = 0.f;
131  ticl::mapWithFraction layerClusterToAssociatedTracksterMap(layerClustersIds.size());
132  std::vector<unsigned int> associatedRecoTracksterIndices;
133  for (unsigned int i = 0; i < layerClustersIds.size(); ++i) {
134  unsigned int layerClusterId = layerClustersIds[i];
135  const auto& layerCluster = layerClusters[layerClusterId];
136  float simFraction = 1.f / simTrackster.vertex_multiplicity(i);
137  float squaredSimFraction = simFraction * simFraction;
138  float squaredLayerClusterEnergy = layerCluster.energy() * layerCluster.energy();
139  simToRecoScoresDenominator += squaredLayerClusterEnergy * squaredSimFraction;
140  const auto& recoTracksterVec = layerClusterToTracksterMap[layerClusterId];
141  for (const auto& [recoTracksterIndex, recoSharedEnergy] : recoTracksterVec) {
142  layerClusterToAssociatedTracksterMap[i].push_back({recoTracksterIndex, recoSharedEnergy});
143  associatedRecoTracksterIndices.push_back(recoTracksterIndex);
144  }
145  }
146  // keep only unique associatedRecoTracksterIndices
147  std::sort(associatedRecoTracksterIndices.begin(), associatedRecoTracksterIndices.end());
148  associatedRecoTracksterIndices.erase(
149  std::unique(associatedRecoTracksterIndices.begin(), associatedRecoTracksterIndices.end()),
150  associatedRecoTracksterIndices.end());
151  // for each layer cluster, loop over associatedRecoTracksterIndices and add the missing reco tracksters with 0 shared energy to layerClusterToAssociatedTracksterMap
152  for (unsigned int i = 0; i < layerClustersIds.size(); ++i) {
153  unsigned int layerClusterId = layerClustersIds[i];
154  const auto& recoTracksterVec = layerClusterToTracksterMap[layerClusterId];
155  for (unsigned int recoTracksterIndex : associatedRecoTracksterIndices) {
156  if (std::find_if(recoTracksterVec.begin(), recoTracksterVec.end(), [recoTracksterIndex](const auto& pair) {
157  return pair.first == recoTracksterIndex;
158  }) == recoTracksterVec.end()) {
159  layerClusterToAssociatedTracksterMap[i].push_back({recoTracksterIndex, 0.f});
160  }
161  }
162  }
163 
164  const float invDenominator = 1.f / simToRecoScoresDenominator;
165 
166  for (unsigned int i = 0; i < layerClustersIds.size(); ++i) {
167  unsigned int layerClusterId = layerClustersIds[i];
168  const auto& layerCluster = layerClusters[layerClusterId];
169  float simFraction = 1.f / simTrackster.vertex_multiplicity(i);
170  float squaredSimFraction = simFraction * simFraction;
171  float squaredLayerClusterEnergy = layerCluster.energy() * layerCluster.energy();
172  float simSharedEnergy = layerCluster.energy() * simFraction;
173  float invLayerClusterEnergy = 1.f / layerCluster.energy();
174  const auto& recoTracksterVec = layerClusterToAssociatedTracksterMap[i];
175  for (const auto& [recoTracksterIndex, recoSharedEnergy] : recoTracksterVec) {
176  edm::Ref<std::vector<ticl::Trackster>> recoTracksterRef(recoTrackstersHandle, recoTracksterIndex);
177  float sharedEnergy = std::min(recoSharedEnergy, simSharedEnergy);
178  float recoFraction = recoSharedEnergy * invLayerClusterEnergy;
179  float score = invDenominator *
180  std::min(squaredSimFraction, (simFraction - recoFraction) * (simFraction - recoFraction)) *
181  squaredLayerClusterEnergy;
182  simTracksterToTracksterMap->insert(tracksterIndex, recoTracksterIndex, sharedEnergy, score);
183  }
184  }
185  }
186  tracksterToSimTracksterMap->sort(true);
187  simTracksterToTracksterMap->sort(true);
188  iEvent.put(std::move(tracksterToSimTracksterMap), "tracksterToSimTracksterMap");
189  iEvent.put(std::move(simTracksterToTracksterMap), "simTracksterToTracksterMap");
190 }
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:30
mps_fire.i
i
Definition: mps_fire.py:429
offlineSlimmedPrimaryVertices_cfi.score
score
Definition: offlineSlimmedPrimaryVertices_cfi.py:6
TracksterToSimTracksterAssociatorProducer::layerClustersCollectionToken_
edm::EDGetTokenT< std::vector< reco::CaloCluster > > layerClustersCollectionToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:35
hltEgammaHGCALIDVarsL1Seeded_cfi.layerClusters
layerClusters
Definition: hltEgammaHGCALIDVarsL1Seeded_cfi.py:5
TracksterToSimTracksterAssociatorProducer::simTracksterCollectionToken_
edm::EDGetTokenT< std::vector< ticl::Trackster > > simTracksterCollectionToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:34
ticl::mapWithFraction
std::vector< std::vector< std::pair< unsigned int, float > >> mapWithFraction
Definition: TICLAssociationMap.h:19
iEvent
int iEvent
Definition: GenABIO.cc:224
tier0.unique
def unique(seq, keepstr=True)
Definition: tier0.py:24
TracksterToSimTracksterAssociatorProducer::LayerClusterToTracksterMapToken_
edm::EDGetTokenT< ticl::AssociationMap< ticl::mapWithFraction, std::vector< reco::CaloCluster >, std::vector< ticl::Trackster > > > LayerClusterToTracksterMapToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:38
TracksterToSimTracksterAssociatorProducer::LayerClusterToSimTracksterMapToken_
edm::EDGetTokenT< ticl::AssociationMap< ticl::mapWithFraction, std::vector< reco::CaloCluster >, std::vector< ticl::Trackster > > > LayerClusterToSimTracksterMapToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:41
egamma::sharedEnergy
float sharedEnergy(reco::CaloCluster const &clu1, reco::CaloCluster const &clu2, EcalRecHitCollection const &barrelRecHits, EcalRecHitCollection const &endcapRecHits)
Definition: EgAmbiguityTools.cc:140
TracksterToSimTracksterAssociatorProducer::recoTracksterCollectionToken_
edm::EDGetTokenT< std::vector< ticl::Trackster > > recoTracksterCollectionToken_
Definition: TracksterToSimTracksterAssociatorProducer.h:33
eostools.move
def move(src, dest)
Definition: eostools.py:511

Member Data Documentation

◆ layerClustersCollectionToken_

edm::EDGetTokenT<std::vector<reco::CaloCluster> > TracksterToSimTracksterAssociatorProducer::layerClustersCollectionToken_
private

Definition at line 35 of file TracksterToSimTracksterAssociatorProducer.h.

Referenced by produce().

◆ LayerClusterToSimTracksterMapToken_

edm::EDGetTokenT< ticl::AssociationMap<ticl::mapWithFraction, std::vector<reco::CaloCluster>, std::vector<ticl::Trackster> > > TracksterToSimTracksterAssociatorProducer::LayerClusterToSimTracksterMapToken_
private

Definition at line 41 of file TracksterToSimTracksterAssociatorProducer.h.

Referenced by produce().

◆ LayerClusterToTracksterMapToken_

edm::EDGetTokenT< ticl::AssociationMap<ticl::mapWithFraction, std::vector<reco::CaloCluster>, std::vector<ticl::Trackster> > > TracksterToSimTracksterAssociatorProducer::LayerClusterToTracksterMapToken_
private

Definition at line 38 of file TracksterToSimTracksterAssociatorProducer.h.

Referenced by produce().

◆ recoTracksterCollectionToken_

edm::EDGetTokenT<std::vector<ticl::Trackster> > TracksterToSimTracksterAssociatorProducer::recoTracksterCollectionToken_
private

Definition at line 33 of file TracksterToSimTracksterAssociatorProducer.h.

Referenced by produce().

◆ simTracksterCollectionToken_

edm::EDGetTokenT<std::vector<ticl::Trackster> > TracksterToSimTracksterAssociatorProducer::simTracksterCollectionToken_
private

Definition at line 34 of file TracksterToSimTracksterAssociatorProducer.h.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.14