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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), and HLT_2018_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 begin, end, HLT_2018_cff::fraction, edm::Event::getByToken(), reco::CaloCluster::hgcal_em, mps_fire::i, charmTagsComputerCvsB_cfi::idx, layer_clusters_token_, hltEgammaHGCALIDVarsL1Seeded_cfi::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_