Inheritance diagram for MultiClustersFromTrackstersProducer:

Public Member Functions
	MultiClustersFromTrackstersProducer (const edm::ParameterSet &)

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

	~MultiClustersFromTrackstersProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
edm::EDGetTokenT< std::vector< reco::CaloCluster > >	layer_clusters_token_

edm::EDGetTokenT< std::vector< ticl::Trackster > >	tracksters_token_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 18 of file MultiClustersFromTrackstersProducer.cc.

Constructor & Destructor Documentation

◆ MultiClustersFromTrackstersProducer()

MultiClustersFromTrackstersProducer::MultiClustersFromTrackstersProducer ( const edm::ParameterSet & ps )

Definition at line 33 of file MultiClustersFromTrackstersProducer.cc.

     : layer_clusters_token_(consumes<std::vector<reco::CaloCluster>>(ps.getParameter<edm::InputTag>("LayerClusters"))),
       tracksters_token_(consumes<std::vector<ticl::Trackster>>(ps.getParameter<edm::InputTag>("Tracksters"))) {
   produces<std::vector<reco::HGCalMultiCluster>>();
 }

◆ ~MultiClustersFromTrackstersProducer()

MultiClustersFromTrackstersProducer::~MultiClustersFromTrackstersProducer ( )

inlineoverride

Definition at line 21 of file MultiClustersFromTrackstersProducer.cc.

21 {}

Member Function Documentation

◆ fillDescriptions()

void MultiClustersFromTrackstersProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 39 of file MultiClustersFromTrackstersProducer.cc.

                                                                                                      {
   // hgcalMultiClusters
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("Tracksters", edm::InputTag("Tracksters", "TrackstersByCA"));
   desc.add<edm::InputTag>("LayerClusters", edm::InputTag("hgcalLayerClusters"));
   desc.addUntracked<unsigned int>("verbosity", 3);
   descriptions.add("multiClustersFromTrackstersProducer", desc);
 }

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

◆ produce()

void MultiClustersFromTrackstersProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 48 of file MultiClustersFromTrackstersProducer.cc.

                                                                                         {
   auto multiclusters = std::make_unique<std::vector<reco::HGCalMultiCluster>>();
   edm::Handle<std::vector<ticl::Trackster>> tracksterHandle;
   evt.getByToken(tracksters_token_, tracksterHandle);
  
   edm::Handle<std::vector<reco::CaloCluster>> layer_clustersHandle;
   evt.getByToken(layer_clusters_token_, layer_clustersHandle);
  
   auto const& tracksters = *tracksterHandle;
   auto const& layerClusters = *layer_clustersHandle;
  
   edm::PtrVector<reco::BasicCluster> clusterPtrs;
   for (unsigned i = 0; i < layerClusters.size(); ++i) {
     edm::Ptr<reco::BasicCluster> ptr(layer_clustersHandle, i);
     clusterPtrs.push_back(ptr);
   }
  
   std::for_each(std::begin(tracksters), std::end(tracksters), [&](auto const& trackster) {
     // Create an empty multicluster if the trackster has no layer clusters.
     // This could happen when a seed leads to no trackster and a dummy one is produced.
  
     std::array<double, 3> baricenter{{0., 0., 0.}};
     double total_weight = 0.;
     reco::HGCalMultiCluster temp;
     int counter = 0;
     if (!trackster.vertices().empty()) {
       std::for_each(std::begin(trackster.vertices()), std::end(trackster.vertices()), [&](unsigned int idx) {
         temp.push_back(clusterPtrs[idx]);
         auto fraction = 1.f / trackster.vertex_multiplicity(counter++);
         for (auto const& cell : clusterPtrs[idx]->hitsAndFractions()) {
           temp.addHitAndFraction(cell.first, cell.second * fraction);
         }
         auto weight = clusterPtrs[idx]->energy() * fraction;
         total_weight += weight;
         baricenter[0] += clusterPtrs[idx]->x() * weight;
         baricenter[1] += clusterPtrs[idx]->y() * weight;
         baricenter[2] += clusterPtrs[idx]->z() * weight;
       });
       std::transform(
           std::begin(baricenter), std::end(baricenter), std::begin(baricenter), [&total_weight](double val) -> double {
             return val / total_weight;
           });
     }
     temp.setEnergy(total_weight);
     temp.setCorrectedEnergy(total_weight);
     temp.setPosition(math::XYZPoint(baricenter[0], baricenter[1], baricenter[2]));
     temp.setAlgoId(reco::CaloCluster::hgcal_em);
     temp.setTime(trackster.time(), trackster.timeError());
     multiclusters->push_back(temp);
   });
  
   evt.put(std::move(multiclusters));
 }

References mps_fire::end, HLT_FULL_cff::fraction, edm::Event::getByToken(), reco::CaloCluster::hgcal_em, mps_fire::i, heavyIonCSV_trainingSettings::idx, layer_clusters_token_, HLTEgPhaseIITestSequence_cff::layerClusters, eostools::move(), edm::PtrVector< T >::push_back(), edm::Event::put(), groupFilesInBlocks::temp, tracksters_token_, HcalDetIdTransform::transform(), heppy_batch::val, and mps_merge::weight.

Member Data Documentation

◆ layer_clusters_token_

edm::EDGetTokenT<std::vector<reco::CaloCluster> > MultiClustersFromTrackstersProducer::layer_clusters_token_

private

Definition at line 27 of file MultiClustersFromTrackstersProducer.cc.

Referenced by produce().

◆ tracksters_token_

edm::EDGetTokenT<std::vector<ticl::Trackster> > MultiClustersFromTrackstersProducer::tracksters_token_