Basic2DGenericTopoClusterizer Class Reference

Inheritance diagram for Basic2DGenericTopoClusterizer:

Public Member Functions
	Basic2DGenericTopoClusterizer (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
	Basic2DGenericTopoClusterizer (const B2DGT &)=delete
void	buildClusters (const edm::Handle< reco::PFRecHitCollection > &, const std::vector< bool > &, const std::vector< bool > &, reco::PFClusterCollection &) override
B2DGT &	operator= (const B2DGT &)=delete
	~Basic2DGenericTopoClusterizer () override=default
Public Member Functions inherited from InitialClusteringStepBase
	InitialClusteringStepBase (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
	InitialClusteringStepBase (const ICSB &)=delete
std::ostream &	operator<< (std::ostream &o) const
ICSB &	operator= (const ICSB &)=delete
void	reset ()
virtual void	update (const edm::EventSetup &)
virtual void	updateEvent (const edm::Event &)
virtual	~InitialClusteringStepBase ()=default

Private Types
typedef Basic2DGenericTopoClusterizer	B2DGT

Private Member Functions
void	buildTopoCluster (const edm::Handle< reco::PFRecHitCollection > &, const std::vector< bool > &, unsigned int, std::vector< bool > &, reco::PFCluster &)

Private Attributes
const bool	_useCornerCells

Additional Inherited Members
Protected Types inherited from InitialClusteringStepBase
typedef std::tuple< std::vector< int >, std::vector< double >, std::vector< double > >	I3tuple
Protected Member Functions inherited from InitialClusteringStepBase
reco::PFRecHitRef	makeRefhit (const edm::Handle< reco::PFRecHitCollection > &h, const unsigned i) const
Protected Attributes inherited from InitialClusteringStepBase
const std::unordered_map< std::string, int >	_layerMap
unsigned	_nClustersFound
unsigned	_nSeeds
std::unordered_map< int, I3tuple >	_thresholds

Detailed Description

Definition at line 5 of file Basic2DGenericTopoClusterizer.cc.

Member Typedef Documentation

B2DGT

typedef Basic2DGenericTopoClusterizer Basic2DGenericTopoClusterizer::B2DGT

private

Definition at line 6 of file Basic2DGenericTopoClusterizer.cc.

Constructor & Destructor Documentation

Basic2DGenericTopoClusterizer() [1/2]

Basic2DGenericTopoClusterizer::Basic2DGenericTopoClusterizer ( const edm::ParameterSet &  conf, edm::ConsumesCollector &  cc  )

inline

Definition at line 9 of file Basic2DGenericTopoClusterizer.cc.

      : InitialClusteringStepBase(conf, cc), _useCornerCells(conf.getParameter<bool>("useCornerCells")) {}

~Basic2DGenericTopoClusterizer()

Basic2DGenericTopoClusterizer::~Basic2DGenericTopoClusterizer ( )

overridedefault

Basic2DGenericTopoClusterizer() [2/2]

Basic2DGenericTopoClusterizer::Basic2DGenericTopoClusterizer ( const B2DGT &  )

delete

Member Function Documentation

buildClusters()

void Basic2DGenericTopoClusterizer::buildClusters ( const edm::Handle< reco::PFRecHitCollection > &  input, const std::vector< bool > &  rechitMask, const std::vector< bool > &  seedable, reco::PFClusterCollection &  output  )

overridevirtual

Implements InitialClusteringStepBase.

Definition at line 43 of file Basic2DGenericTopoClusterizer.cc.

References buildTopoCluster(), hfClusterShapes_cfi::hits, mps_fire::i, input, dqmiolumiharvest::j, fileCollector::seed, HLT_FULL_cff::seeds, jetUpdater_cfi::sort, and groupFilesInBlocks::temp.

                                                                                   {
  auto const& hits = *input;
  std::vector<bool> used(hits.size(), false);
  std::vector<unsigned int> seeds;

  // get the seeds and sort them descending in energy
  seeds.reserve(hits.size());
  for (unsigned int i = 0; i < hits.size(); ++i) {
    if (!rechitMask[i] || !seedable[i] || used[i])
      continue;
    seeds.emplace_back(i);
  }
  // maxHeap would be better
  std::sort(
      seeds.begin(), seeds.end(), [&](unsigned int i, unsigned int j) { return hits[i].energy() > hits[j].energy(); });

  reco::PFCluster temp;
  for (auto seed : seeds) {
    if (!rechitMask[seed] || !seedable[seed] || used[seed])
      continue;
    temp.reset();
    buildTopoCluster(input, rechitMask, seed, used, temp);
    if (!temp.recHitFractions().empty())
      output.push_back(temp);
  }
}

buildTopoCluster()

void Basic2DGenericTopoClusterizer::buildTopoCluster ( const edm::Handle< reco::PFRecHitCollection > &  input, const std::vector< bool > &  rechitMask, unsigned int  kcell, std::vector< bool > &  used, reco::PFCluster &  topocluster  )

private

Definition at line 73 of file Basic2DGenericTopoClusterizer.cc.

References InitialClusteringStepBase::_thresholds, _useCornerCells, funct::abs(), reco::PFCluster::addRecHitFraction(), hcalRecHitTable_cff::depth, PFLayer::HCAL_BARREL1, PFLayer::HCAL_BARREL2, PFLayer::HCAL_ENDCAP, input, createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, LOGDRESSED, InitialClusteringStepBase::makeRefhit(), and particleFlowZeroSuppressionECAL_cff::thresholds.

Referenced by buildClusters().

                                                                                 {
  auto const& cell = (*input)[kcell];
  int cell_layer = (int)cell.layer();
  if (cell_layer == PFLayer::HCAL_BARREL2 && std::abs(cell.positionREP().eta()) > 0.34) {
    cell_layer *= 100;
  }

  auto const& thresholds = _thresholds.find(cell_layer)->second;
  double thresholdE = 0.;
  double thresholdPT2 = 0.;

  for (unsigned int j = 0; j < (std::get<1>(thresholds)).size(); ++j) {
    int depth = std::get<0>(thresholds)[j];

    if ((cell_layer == PFLayer::HCAL_BARREL1 && cell.depth() == depth) ||
        (cell_layer == PFLayer::HCAL_ENDCAP && cell.depth() == depth) ||
        (cell_layer != PFLayer::HCAL_BARREL1 && cell_layer != PFLayer::HCAL_ENDCAP)) {
      thresholdE = std::get<1>(thresholds)[j];
      thresholdPT2 = std::get<2>(thresholds)[j];
    }
  }

  if (cell.energy() < thresholdE || cell.pt2() < thresholdPT2) {
    LOGDRESSED("GenericTopoCluster::buildTopoCluster()")
        << "RecHit " << cell.detId() << " with enegy " << cell.energy() << " GeV was rejected!." << std::endl;
    return;
  }

  auto k = kcell;
  used[k] = true;
  auto ref = makeRefhit(input, k);
  topocluster.addRecHitFraction(reco::PFRecHitFraction(ref, 1.0));

  auto const& neighbours = (_useCornerCells ? cell.neighbours8() : cell.neighbours4());

  for (auto nb : neighbours) {
    if (used[nb] || !rechitMask[nb]) {
      LOGDRESSED("GenericTopoCluster::buildTopoCluster()")
          << "  RecHit " << cell.detId() << "\'s"
          << " neighbor RecHit " << input->at(nb).detId() << " with enegy " << input->at(nb).energy()
          << " GeV was rejected!"
          << " Reasons : " << used[nb] << " (used) " << !rechitMask[nb] << " (masked)." << std::endl;
      continue;
    }
    buildTopoCluster(input, rechitMask, nb, used, topocluster);
  }
}

operator=()

B2DGT& Basic2DGenericTopoClusterizer::operator= ( const B2DGT &  )

delete

Member Data Documentation

_useCornerCells

const bool Basic2DGenericTopoClusterizer::_useCornerCells

private

Definition at line 21 of file Basic2DGenericTopoClusterizer.cc.

Referenced by buildTopoCluster().