Basic2DGenericTopoClusterizer Class Reference

#include <Basic2DGenericTopoClusterizer.h>

Inheritance diagram for Basic2DGenericTopoClusterizer:

Public Member Functions
	Basic2DGenericTopoClusterizer (const edm::ParameterSet &conf)
	Basic2DGenericTopoClusterizer (const B2DGT &)=delete
void	buildClusters (const edm::Handle< reco::PFRecHitCollection > &, const std::vector< bool > &, const std::vector< bool > &, reco::PFClusterCollection &)
B2DGT &	operator= (const B2DGT &)=delete
virtual	~Basic2DGenericTopoClusterizer ()
Public Member Functions inherited from InitialClusteringStepBase
	_algoName (conf.getParameter< std::string >("algoName"))
	InitialClusteringStepBase (const edm::ParameterSet &conf)
	InitialClusteringStepBase (const ICSB &)=delete
std::ostream &	operator<< (std::ostream &o)
ICSB &	operator= (const ICSB &)=delete
void	reset ()
virtual void	update (const edm::EventSetup &)
virtual	~InitialClusteringStepBase ()

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 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, std::pair< double, double > >	_thresholds

Detailed Description

Definition at line 7 of file Basic2DGenericTopoClusterizer.h.

Member Typedef Documentation

typedef Basic2DGenericTopoClusterizer Basic2DGenericTopoClusterizer::B2DGT

private

Definition at line 8 of file Basic2DGenericTopoClusterizer.h.

Constructor & Destructor Documentation

Basic2DGenericTopoClusterizer::Basic2DGenericTopoClusterizer ( const edm::ParameterSet &  conf )

inline

Definition at line 10 of file Basic2DGenericTopoClusterizer.h.

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

virtual Basic2DGenericTopoClusterizer::~Basic2DGenericTopoClusterizer ( )

inlinevirtual

Definition at line 13 of file Basic2DGenericTopoClusterizer.h.

{}

Basic2DGenericTopoClusterizer::Basic2DGenericTopoClusterizer ( const B2DGT &  )

delete

Member Function Documentation

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

virtual

Implements InitialClusteringStepBase.

Definition at line 17 of file Basic2DGenericTopoClusterizer.cc.

References buildTopoCluster(), i, input, j, fileCollector::seed, 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().size() ) output.push_back(temp);
  }
}

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 45 of file Basic2DGenericTopoClusterizer.cc.

References InitialClusteringStepBase::_thresholds, _useCornerCells, funct::abs(), reco::PFCluster::addRecHitFraction(), PFLayer::HCAL_BARREL2, relval_2017::k, LOGDRESSED, and InitialClusteringStepBase::makeRefhit().

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;
    }    
  const std::pair<double,double>& thresholds =
      _thresholds.find(cell_layer)->second;
  if( cell.energy() < thresholds.first || 
      cell.pt2() < thresholds.second ) {
    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);
  }
}

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

delete

Member Data Documentation

const bool Basic2DGenericTopoClusterizer::_useCornerCells

private

Definition at line 23 of file Basic2DGenericTopoClusterizer.h.

Referenced by buildTopoCluster().