HGCalMulticlusteringImpl Class Reference

#include <HGCalMulticlusteringImpl.h>

Public Member Functions
void	clusterizeDBSCAN (const std::vector< edm::Ptr< l1t::HGCalCluster >> &clustersPtr, l1t::HGCalMulticlusterBxCollection &multiclusters, const HGCalTriggerGeometryBase &triggerGeometry)
void	clusterizeDR (const std::vector< edm::Ptr< l1t::HGCalCluster >> &clustersPtr, l1t::HGCalMulticlusterBxCollection &multiclusters, const HGCalTriggerGeometryBase &triggerGeometry)
	HGCalMulticlusteringImpl (const edm::ParameterSet &conf)
bool	isPertinent (const l1t::HGCalCluster &clu, const l1t::HGCalMulticluster &mclu, double dR) const
void	setGeometry (const HGCalTriggerGeometryBase *const geom)

Private Member Functions
void	finalizeClusters (std::vector< l1t::HGCalMulticluster > &, l1t::HGCalMulticlusterBxCollection &, const HGCalTriggerGeometryBase &)
void	findNeighbor (const std::vector< std::pair< unsigned int, double >> &rankedList, unsigned int searchInd, const std::vector< edm::Ptr< l1t::HGCalCluster >> &clustersPtr, std::vector< unsigned int > &neigbors)

Private Attributes
double	distDbscan_ = 0.005
double	dr_
std::unique_ptr< HGCalTriggerClusterIdentificationBase >	id_
unsigned	minNDbscan_ = 3
std::string	multiclusterAlgoType_
double	ptC3dThreshold_
HGCalShowerShape	shape_
HGCalTriggerTools	triggerTools_

Detailed Description

Definition at line 14 of file HGCalMulticlusteringImpl.h.

Constructor & Destructor Documentation

HGCalMulticlusteringImpl()

HGCalMulticlusteringImpl::HGCalMulticlusteringImpl

(

const edm::ParameterSet &

conf

)

Definition at line 5 of file HGCalMulticlusteringImpl.cc.

References distDbscan_, dr_, get, edm::ParameterSet::getParameter(), id_, minNDbscan_, multiclusterAlgoType_, and ptC3dThreshold_.

    : dr_(conf.getParameter<double>("dR_multicluster")),
      ptC3dThreshold_(conf.getParameter<double>("minPt_multicluster")),
      multiclusterAlgoType_(conf.getParameter<string>("type_multicluster")),
      distDbscan_(conf.getParameter<double>("dist_dbscan_multicluster")),
      minNDbscan_(conf.getParameter<unsigned>("minN_dbscan_multicluster")) {
  edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster dR for Near Neighbour search: " << dr_;
  edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster minimum transverse-momentum: " << ptC3dThreshold_;
  edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster DBSCAN Clustering distance: " << distDbscan_;
  edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster clustering min number of subclusters: " << minNDbscan_;
  edm::LogInfo("HGCalMulticlusterParameters")
      << "Multicluster type of multiclustering algortihm: " << multiclusterAlgoType_;
  id_ = std::unique_ptr<HGCalTriggerClusterIdentificationBase>{
      HGCalTriggerClusterIdentificationFactory::get()->create("HGCalTriggerClusterIdentificationBDT")};
  id_->initialize(conf.getParameter<edm::ParameterSet>("EGIdentification"));
}

Member Function Documentation

clusterizeDBSCAN()

void HGCalMulticlusteringImpl::clusterizeDBSCAN	(	const std::vector< edm::Ptr< l1t::HGCalCluster >> &	clustersPtr,
		l1t::HGCalMulticlusterBxCollection &	multiclusters,
		const HGCalTriggerGeometryBase &	triggerGeometry
	)

Definition at line 103 of file HGCalMulticlusteringImpl.cc.

References finalizeClusters(), findNeighbor(), isotrackNtupler::merged, minNDbscan_, eostools::move(), jetUpdater_cfi::sort, triggerTools_, trackerHitRTTI::vector, class-composition::visited, and HGCalTriggerTools::zside().

                                                                                                 {
  std::vector<l1t::HGCalMulticluster> multiclustersTmp;
  l1t::HGCalMulticluster mcluTmp;
  std::vector<bool> visited(clustersPtrs.size(), false);
  std::vector<bool> merged(clustersPtrs.size(), false);
  std::vector<std::pair<unsigned int, double>> rankedList;
  rankedList.reserve(clustersPtrs.size());
  std::vector<std::vector<unsigned int>> neighborList;
  neighborList.reserve(clustersPtrs.size());

  int iclu = 0, imclu = 0, neighNo;
  double dist = 0.;

  for (std::vector<edm::Ptr<l1t::HGCalCluster>>::const_iterator clu = clustersPtrs.begin(); clu != clustersPtrs.end();
       ++clu, ++iclu) {
    dist = (*clu)->centreProj().mag() * triggerTools_.zside((*clu)->detId());
    rankedList.push_back(std::make_pair(iclu, dist));
  }
  iclu = 0;
  std::sort(rankedList.begin(), rankedList.end(), [](auto& left, auto& right) { return left.second < right.second; });

  for (const auto& cluRanked : rankedList) {
    std::vector<unsigned int> neighbors;

    if (!visited.at(iclu)) {
      visited.at(iclu) = true;
      findNeighbor(rankedList, iclu, clustersPtrs, neighbors);
      neighborList.push_back(std::move(neighbors));

      if (neighborList.at(iclu).size() >= minNDbscan_) {
        multiclustersTmp.emplace_back(clustersPtrs[cluRanked.first]);
        merged.at(iclu) = true;
        /* dynamic range loop: range-based loop syntax cannot be employed */
        for (unsigned int neighInd = 0; neighInd < neighborList.at(iclu).size(); neighInd++) {
          neighNo = neighborList.at(iclu).at(neighInd);
          /* This condition also ensures merging of clusters visited by other clusters but not merged. */
          if (!merged.at(neighNo)) {
            merged.at(neighNo) = true;
            multiclustersTmp.at(imclu).addConstituent(clustersPtrs[rankedList.at(neighNo).first]);

            if (!visited.at(neighNo)) {
              visited.at(neighNo) = true;
              std::vector<unsigned int> secNeighbors;
              findNeighbor(rankedList, neighNo, clustersPtrs, secNeighbors);

              if (secNeighbors.size() >= minNDbscan_) {
                neighborList.at(iclu).insert(neighborList.at(iclu).end(), secNeighbors.begin(), secNeighbors.end());
              }
            }
          }
        }
        imclu++;
      }
    } else
      neighborList.push_back(std::move(neighbors));
    iclu++;
  }
  /* making the collection of multiclusters */
  finalizeClusters(multiclustersTmp, multiclusters, triggerGeometry);
}

clusterizeDR()

void HGCalMulticlusteringImpl::clusterizeDR	(	const std::vector< edm::Ptr< l1t::HGCalCluster >> &	clustersPtr,
		l1t::HGCalMulticlusterBxCollection &	multiclusters,
		const HGCalTriggerGeometryBase &	triggerGeometry
	)

Definition at line 72 of file HGCalMulticlusteringImpl.cc.

References ztail::d, dr_, finalizeClusters(), createfilelist::int, isPertinent(), mag(), and trackerHitRTTI::vector.

                                                                                             {
  std::vector<l1t::HGCalMulticluster> multiclustersTmp;

  int iclu = 0;
  for (std::vector<edm::Ptr<l1t::HGCalCluster>>::const_iterator clu = clustersPtrs.begin(); clu != clustersPtrs.end();
       ++clu, ++iclu) {
    double minDist = dr_;
    int targetMulticlu = -1;

    for (unsigned imclu = 0; imclu < multiclustersTmp.size(); imclu++) {
      if (!this->isPertinent(**clu, multiclustersTmp.at(imclu), dr_))
        continue;

      double d = (multiclustersTmp.at(imclu).centreProj() - (*clu)->centreProj()).mag();
      if (d < minDist) {
        minDist = d;
        targetMulticlu = int(imclu);
      }
    }

    if (targetMulticlu < 0)
      multiclustersTmp.emplace_back(*clu);
    else
      multiclustersTmp.at(targetMulticlu).addConstituent(*clu);
  }

  /* making the collection of multiclusters */
  finalizeClusters(multiclustersTmp, multiclusters, triggerGeometry);
}

finalizeClusters()

void HGCalMulticlusteringImpl::finalizeClusters ( std::vector< l1t::HGCalMulticluster > &  multiclusters_in, l1t::HGCalMulticlusterBxCollection &  multiclusters_out, const HGCalTriggerGeometryBase &  triggerGeometry  )

private

Definition at line 166 of file HGCalMulticlusteringImpl.cc.

References bsc_activity_cfg::clusters, HGCalShowerShape::fillShapes(), id_, ptC3dThreshold_, BXVector< T >::push_back(), shape_, and TtFullHadEvtBuilder_cfi::sumPt.

Referenced by clusterizeDBSCAN(), and clusterizeDR().

                                                                                                 {
  for (auto& multicluster : multiclusters_in) {
    // compute the eta, phi from its barycenter
    // + pT as scalar sum of pT of constituents
    double sumPt = 0.;
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clusters = multicluster.constituents();
    for (const auto& id_cluster : clusters)
      sumPt += id_cluster.second->pt();

    math::PtEtaPhiMLorentzVector multiclusterP4(sumPt, multicluster.centre().eta(), multicluster.centre().phi(), 0.);
    multicluster.setP4(multiclusterP4);

    if (multicluster.pt() > ptC3dThreshold_) {
      //compute shower shapes
      shape_.fillShapes(multicluster, triggerGeometry);
      // fill quality flag
      multicluster.setHwQual(id_->decision(multicluster));
      // fill H/E
      multicluster.saveHOverE();

      multiclusters_out.push_back(0, multicluster);
    }
  }
}

findNeighbor()

void HGCalMulticlusteringImpl::findNeighbor ( const std::vector< std::pair< unsigned int, double >> &  rankedList, unsigned int  searchInd, const std::vector< edm::Ptr< l1t::HGCalCluster >> &  clustersPtr, std::vector< unsigned int > &  neigbors  )

private

Definition at line 37 of file HGCalMulticlusteringImpl.cc.

References distDbscan_, Exception, and mag().

Referenced by clusterizeDBSCAN().

                                                                                {
  if (clustersPtrs.size() <= searchInd || clustersPtrs.size() < rankedList.size()) {
    throw cms::Exception("IndexOutOfBound: clustersPtrs in 'findNeighbor'");
  }

  for (unsigned int ind = searchInd + 1;
       ind < rankedList.size() && fabs(rankedList.at(ind).second - rankedList.at(searchInd).second) < distDbscan_;
       ind++) {
    if (clustersPtrs.size() <= rankedList.at(ind).first) {
      throw cms::Exception("IndexOutOfBound: clustersPtrs in 'findNeighbor'");

    } else if (((*(clustersPtrs[rankedList.at(ind).first])).centreProj() -
                (*(clustersPtrs[rankedList.at(searchInd).first])).centreProj())
                   .mag() < distDbscan_) {
      neighbors.push_back(ind);
    }
  }

  for (unsigned int ind = 0;
       ind < searchInd && fabs(rankedList.at(searchInd).second - rankedList.at(ind).second) < distDbscan_;
       ind++) {
    if (clustersPtrs.size() <= rankedList.at(ind).first) {
      throw cms::Exception("IndexOutOfBound: clustersPtrs in 'findNeighbor'");

    } else if (((*(clustersPtrs[rankedList.at(ind).first])).centreProj() -
                (*(clustersPtrs[rankedList.at(searchInd).first])).centreProj())
                   .mag() < distDbscan_) {
      neighbors.push_back(ind);
    }
  }
}

isPertinent()

bool HGCalMulticlusteringImpl::isPertinent	(	const l1t::HGCalCluster &	clu,
		const l1t::HGCalMulticluster &	mclu,
		double	dR
	)		const

Definition at line 22 of file HGCalMulticlusteringImpl.cc.

References l1t::HGCalClusterT< C >::centreProj(), l1t::HGCalClusterT< C >::detId(), HGC3DClusterGenMatchSelector_cfi::dR, mag(), triggerTools_, and HGCalTriggerTools::zside().

Referenced by clusterizeDR().

                                                            {
  DetId cluDetId(clu.detId());
  DetId firstClusterDetId(mclu.detId());

  if (triggerTools_.zside(cluDetId) != triggerTools_.zside(firstClusterDetId)) {
    return false;
  }
  if ((mclu.centreProj() - clu.centreProj()).mag() < dR) {
    return true;
  }
  return false;
}

setGeometry()

void HGCalMulticlusteringImpl::setGeometry ( const HGCalTriggerGeometryBase *const  geom )

inline

Definition at line 18 of file HGCalMulticlusteringImpl.h.

References relativeConstraints::geom, HGCalShowerShape::setGeometry(), HGCalTriggerTools::setGeometry(), shape_, and triggerTools_.

                                                               {
    triggerTools_.setGeometry(geom);
    shape_.setGeometry(geom);
  }

Member Data Documentation

distDbscan_

double HGCalMulticlusteringImpl::distDbscan_ = 0.005

private

Definition at line 45 of file HGCalMulticlusteringImpl.h.

Referenced by findNeighbor(), and HGCalMulticlusteringImpl().

dr_

double HGCalMulticlusteringImpl::dr_

private

Definition at line 42 of file HGCalMulticlusteringImpl.h.

Referenced by clusterizeDR(), and HGCalMulticlusteringImpl().

id_

std::unique_ptr<HGCalTriggerClusterIdentificationBase> HGCalMulticlusteringImpl::id_

private

Definition at line 50 of file HGCalMulticlusteringImpl.h.

Referenced by finalizeClusters(), and HGCalMulticlusteringImpl().

minNDbscan_

unsigned HGCalMulticlusteringImpl::minNDbscan_ = 3

private

Definition at line 46 of file HGCalMulticlusteringImpl.h.

Referenced by clusterizeDBSCAN(), and HGCalMulticlusteringImpl().

multiclusterAlgoType_

std::string HGCalMulticlusteringImpl::multiclusterAlgoType_

private

Definition at line 44 of file HGCalMulticlusteringImpl.h.

Referenced by HGCalMulticlusteringImpl().

ptC3dThreshold_

double HGCalMulticlusteringImpl::ptC3dThreshold_

private

Definition at line 43 of file HGCalMulticlusteringImpl.h.

Referenced by finalizeClusters(), and HGCalMulticlusteringImpl().

shape_

HGCalShowerShape HGCalMulticlusteringImpl::shape_

private

Definition at line 48 of file HGCalMulticlusteringImpl.h.

Referenced by finalizeClusters(), and setGeometry().

triggerTools_

HGCalTriggerTools HGCalMulticlusteringImpl::triggerTools_

private

Definition at line 49 of file HGCalMulticlusteringImpl.h.

Referenced by clusterizeDBSCAN(), isPertinent(), and setGeometry().