#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, std::vector< float > radii_in, uint32_t min_clusters, bool real_space_cone)

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

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

uint32_t	minClusters

std::vector< float >	radii

bool	realSpaceCone

Static Private Attributes
static const unsigned int	lastLayerBH = 52

static const unsigned int	lastLayerEE = 28

static const unsigned int	lastLayerFH = 40

Detailed Description

Definition at line 13 of file HGCalDepthPreClusterer.h.

Member Typedef Documentation

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

Definition at line 31 of file HGCalDepthPreClusterer.h.

Constructor & Destructor Documentation

HGCalDepthPreClusterer::HGCalDepthPreClusterer ( )

inline

Definition at line 17 of file HGCalDepthPreClusterer.h.

References clusterTools, minClusters, and realSpaceCone.

                           : radii({0.,0.,0.,}), minClusters(0), realSpaceCone(false), clusterTools(nullptr)
     {
   }

HGCalDepthPreClusterer::HGCalDepthPreClusterer	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	sumes,
		std::vector< float >	radii_in,
		uint32_t	min_clusters,
		bool	real_space_cone
	)

inline

Definition at line 21 of file HGCalDepthPreClusterer.h.

                                                                                                                                                           :
   radii(radii_in),
   minClusters(min_clusters),
   realSpaceCone(real_space_cone),
   clusterTools(std::make_unique<hgcal::ClusterTools>(conf,sumes)) {
   }

Member Function Documentation

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

inline

Definition at line 28 of file HGCalDepthPreClusterer.h.

References clusterTools.

28 { clusterTools->getEvent(ev); }

HGCalDepthPreClusterer::clusterTools

std::unique_ptr< hgcal::ClusterTools > clusterTools

Definition: HGCalDepthPreClusterer.h:41

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

inline

Definition at line 29 of file HGCalDepthPreClusterer.h.

References clusterTools.

29 { clusterTools->getEventSetup(es); }

HGCalDepthPreClusterer::clusterTools

std::unique_ptr< hgcal::ClusterTools > clusterTools

Definition: HGCalDepthPreClusterer.h:41

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

Definition at line 36 of file HGCalDepthPreClusterer.cc.

References clusterTools, mps_fire::i, lastLayerBH, lastLayerEE, lastLayerFH, minClusters, reco::HGCalMultiCluster::push_back(), radii, realSpaceCone, 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;
 
   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) {
           float distanceCheck = 9999.;
           if( realSpaceCone ) distanceCheck = distAxisCluster2(thecls[es[i]],thecls[es[j]]);
           else distanceCheck = dist2(thecls[es[i]],thecls[es[j]]);
           DetId detid = thecls[es[j]]->hitsAndFractions()[0].first();
           unsigned int layer = clusterTools->getLayer(detid);
           float radius2 = 9999.;
           if(layer <= lastLayerEE) radius2 = radii[0]*radii[0];
           else if(layer <= lastLayerFH) radius2 = radii[1]*radii[1];
           else if(layer <= lastLayerBH) radius2 = radii[2]*radii[2];
           else assert(radius2<100. && "nonsense layer value - cannot assign multicluster radius");
       if( distanceCheck<radius2 && int(thecls[es[j]]->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;
 }