#include <HGCalDepthPreClusterer.h>

Public Types
typedef std::vector < reco::BasicCluster >	ClusterCollection

Public Member Functions
void	getEvent (const edm::Event &ev)

void	getEventSetup (const edm::EventSetup &es)

	HGCalDepthPreClusterer ()

	HGCalDepthPreClusterer (const edm::ParameterSet &conf, edm::ConsumesCollector &sumes, double radius_in, uint32_t min_clusters)

std::vector < reco::HGCalMultiCluster >	makePreClusters (const reco::HGCalMultiCluster::ClusterCollection &) const

Private Attributes
std::unique_ptr < hgcal::ClusterTools >	clusterTools

uint32_t	minClusters

float	radius

Detailed Description

Definition at line 13 of file HGCalDepthPreClusterer.h.

Member Typedef Documentation

typedef std::vector<reco::BasicCluster> HGCalDepthPreClusterer::ClusterCollection

Definition at line 30 of file HGCalDepthPreClusterer.h.

Constructor & Destructor Documentation

HGCalDepthPreClusterer::HGCalDepthPreClusterer ( )

inline

Definition at line 17 of file HGCalDepthPreClusterer.h.

                           : radius(0.), minClusters(0.), clusterTools(nullptr)
     {
   }

HGCalDepthPreClusterer::HGCalDepthPreClusterer	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	sumes,
		double	radius_in,
		uint32_t	min_clusters
	)

inline

Definition at line 21 of file HGCalDepthPreClusterer.h.

                                                                                                                            : 
   radius(radius_in),
   minClusters(min_clusters),
   clusterTools(std::make_unique<hgcal::ClusterTools>(conf,sumes)) {
   }

Member Function Documentation

void HGCalDepthPreClusterer::getEvent ( const edm::Event & ev )

inline

Definition at line 27 of file HGCalDepthPreClusterer.h.

References clusterTools.

27 { clusterTools->getEvent(ev); }

HGCalDepthPreClusterer::clusterTools

std::unique_ptr< hgcal::ClusterTools > clusterTools

Definition: HGCalDepthPreClusterer.h:39

void HGCalDepthPreClusterer::getEventSetup ( const edm::EventSetup & es )

inline

Definition at line 28 of file HGCalDepthPreClusterer.h.

References clusterTools.

28 { clusterTools->getEventSetup(es); }

HGCalDepthPreClusterer::clusterTools

std::unique_ptr< hgcal::ClusterTools > clusterTools

Definition: HGCalDepthPreClusterer.h:39

std::vector< reco::HGCalMultiCluster > HGCalDepthPreClusterer::makePreClusters ( const reco::HGCalMultiCluster::ClusterCollection & thecls ) const

Definition at line 26 of file HGCalDepthPreClusterer.cc.

References clusterTools, i, j, minClusters, reco::HGCalMultiCluster::push_back(), radius, reco::HGCalMultiCluster::size(), sorted_indices(), groupFilesInBlocks::temp, and z.

                                                                                                                                    {
 
   std::vector<reco::HGCalMultiCluster> thePreClusters;
   std::vector<size_t> es = sorted_indices(thecls);
   std::vector<int> vused(es.size(),0);
   unsigned int used = 0;
   const float radius2 = radius*radius;
 
   for(unsigned int i = 0; i < es.size(); ++i) {
     if(vused[i]==0) {
       reco::HGCalMultiCluster temp;      
       temp.push_back(thecls[es[i]]);
       vused[i]=(thecls[es[i]]->z()>0)? 1 : -1;
       ++used;
       for(unsigned int j = i+1; j < es.size(); ++j) {
     if(vused[j]==0) {
       if( dist2(thecls[es[i]],thecls[es[j]]) < radius2 && int(thecls[es[i]]->z()*vused[i])>0 ) {
         temp.push_back(thecls[es[j]]);
         vused[j]=vused[i];
         ++used;
       } 
     }
       }
       if( temp.size() > minClusters ) {
         thePreClusters.push_back(temp);
         auto& back = thePreClusters.back();
         back.setPosition(clusterTools->getMultiClusterPosition(back));
         back.setEnergy(clusterTools->getMultiClusterEnergy(back));
       }
     }
   }
   
   
 
   return thePreClusters;
 }