HGCalDepthPreClusterer Class Reference

#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
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_
std::unique_ptr< hgcal::ClusterTools >	clusterTools
uint32_t	minClusters
std::vector< float >	radii
bool	realSpaceCone
hgcal::RecHitTools	rhtools_

Detailed Description

Definition at line 13 of file HGCalDepthPreClusterer.h.

Member Typedef Documentation

ClusterCollection

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

Definition at line 43 of file HGCalDepthPreClusterer.h.

Constructor & Destructor Documentation

HGCalDepthPreClusterer() [1/2]

HGCalDepthPreClusterer::HGCalDepthPreClusterer ( )

inline

Definition at line 15 of file HGCalDepthPreClusterer.h.

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

HGCalDepthPreClusterer() [2/2]

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 25 of file HGCalDepthPreClusterer.h.

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

Member Function Documentation

getEvent()

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

inline

Definition at line 36 of file HGCalDepthPreClusterer.h.

References clusterTools, and makeMEIFBenchmarkPlots::ev.

{ clusterTools->getEvent(ev); }

getEventSetup()

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

inline

Definition at line 37 of file HGCalDepthPreClusterer.h.

References caloGeomToken_, clusterTools, relativeConstraints::geom, edm::EventSetup::getHandle(), rhtools_, and hgcal::RecHitTools::setGeometry().

                                              {
    clusterTools->getEventSetup(es);
    edm::ESHandle<CaloGeometry> geom = es.getHandle(caloGeomToken_);
    rhtools_.setGeometry(*geom);
  }

makePreClusters()

std::vector< reco::HGCalMultiCluster > HGCalDepthPreClusterer::makePreClusters

(

const reco::HGCalMultiCluster::ClusterCollection &

thecls

)

const

Definition at line 33 of file HGCalDepthPreClusterer.cc.

References clusterTools, hgcal::RecHitTools::firstLayerBH(), mps_fire::i, dqmiolumiharvest::j, hgcal::RecHitTools::lastLayerEE(), nano_mu_digi_cff::layer, minClusters, radii, CosmicsPD_Skims::radius, realSpaceCone, rhtools_, hgcal_clustering::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);

  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;
      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 radius = radii[2];
          if (layer <= rhtools_.lastLayerEE())
            radius = radii[0];
          else if (layer < rhtools_.firstLayerBH())
            radius = radii[1];
          float radius2 = radius * radius;
          if (distanceCheck<radius2 &&int(thecls[es[j]]->z() * vused[i])> 0) {
            temp.push_back(thecls[es[j]]);
            vused[j] = vused[i];
          }
        }
      }
      if (temp.size() > minClusters) {
        thePreClusters.push_back(temp);
        auto &back = thePreClusters.back();
        back.setPosition(clusterTools->getMultiClusterPosition(back));
        back.setEnergy(clusterTools->getMultiClusterEnergy(back));
      }
    }
  }

  return thePreClusters;
}

Member Data Documentation

caloGeomToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HGCalDepthPreClusterer::caloGeomToken_

private

Definition at line 55 of file HGCalDepthPreClusterer.h.

Referenced by getEventSetup().

clusterTools

std::unique_ptr<hgcal::ClusterTools> HGCalDepthPreClusterer::clusterTools

private

Definition at line 53 of file HGCalDepthPreClusterer.h.

Referenced by getEvent(), getEventSetup(), and makePreClusters().

minClusters

uint32_t HGCalDepthPreClusterer::minClusters

private

Definition at line 50 of file HGCalDepthPreClusterer.h.

Referenced by makePreClusters().

radii

std::vector<float> HGCalDepthPreClusterer::radii

private

Definition at line 49 of file HGCalDepthPreClusterer.h.

Referenced by makePreClusters().

realSpaceCone

bool HGCalDepthPreClusterer::realSpaceCone

private

flag to use cartesian space clustering.

Definition at line 51 of file HGCalDepthPreClusterer.h.

Referenced by makePreClusters().

rhtools_

hgcal::RecHitTools HGCalDepthPreClusterer::rhtools_

private

instance of tools to access RecHit information.

Definition at line 54 of file HGCalDepthPreClusterer.h.

Referenced by getEventSetup(), and makePreClusters().