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
std::string	label_

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 ( const edm::ParameterSet & ps )

Definition at line 34 of file MultiClustersFromTrackstersProducer.cc.

References label_.

     : label_(ps.getParameter<std::string>("label")),
       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>>(label_);
 }

MultiClustersFromTrackstersProducer::~MultiClustersFromTrackstersProducer ( )

inlineoverride

Definition at line 21 of file MultiClustersFromTrackstersProducer.cc.

References fillDescriptions(), and produce().

21 {}

Member Function Documentation

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

static

Definition at line 41 of file MultiClustersFromTrackstersProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~MultiClustersFromTrackstersProducer().

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

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

override

Definition at line 51 of file MultiClustersFromTrackstersProducer.cc.

References begin, end, HLT_2018_cff::fraction, edm::Event::getByToken(), reco::CaloCluster::hgcal_em, mps_fire::i, training_settings::idx, label_, layer_clusters_token_, eostools::move(), edm::PtrVector< T >::push_back(), groupFilesInBlocks::temp, tracksters_token_, HcalDetIdTransform::transform(), heppy_batch::val, and mps_merge::weight.

Referenced by ~MultiClustersFromTrackstersProducer().

                                                                                         {
   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) {
     std::array<double, 3> baricenter{{0., 0., 0.}};
     double total_weight = 0.;
     reco::HGCalMultiCluster temp;
     int counter = 0;
     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);
     multiclusters->push_back(temp);
   });
 
   evt.put(std::move(multiclusters), label_);
 }