MultiClustersFromTrackstersProducer Class Reference

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
	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
const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 18 of file MultiClustersFromTrackstersProducer.cc.

Constructor & Destructor Documentation

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

inlineoverride

Definition at line 21 of file MultiClustersFromTrackstersProducer.cc.

{}

Member Function Documentation

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

static

Definition at line 39 of file MultiClustersFromTrackstersProducer.cc.

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

                                                                                                     {
  // 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);
}

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

override

Definition at line 48 of file MultiClustersFromTrackstersProducer.cc.

References reco::CaloCluster::addHitAndFraction(), SplitLinear::begin, counter, dataset::end, HLT_FULL_cff::fraction, edm::Event::getByToken(), reco::CaloCluster::hgcal_em, mps_fire::i, layer_clusters_token_, eostools::move(), reco::HGCalMultiCluster::push_back(), edm::PtrVector< T >::push_back(), edm::Event::put(), reco::CaloCluster::setAlgoId(), reco::CaloCluster::setCorrectedEnergy(), reco::CaloCluster::setEnergy(), reco::CaloCluster::setPosition(), reco::PFCluster::setTime(), groupFilesInBlocks::temp, tracksters_token_, HcalDetIdTransform::transform(), hgcalPerformanceValidation::val, and histoStyle::weight.

                                                                                        {
  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));
}

Member Data Documentation

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

private

Definition at line 27 of file MultiClustersFromTrackstersProducer.cc.

Referenced by produce().

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

private

Definition at line 28 of file MultiClustersFromTrackstersProducer.cc.

Referenced by produce().