HGCGraph Class Reference

#include <HGCGraph.h>

Public Types
enum	VerbosityLevel {   None = 0, Basic, Advanced, Expert,   Guru }

Public Member Functions
bool	areTimeCompatible (int innerIdx, int outerIdx, const edm::ValueMap< float > &layerClustersTime, float maxDeltaTime)
void	clear ()
void	findNtuplets (std::vector< HGCDoublet::HGCntuplet > &foundNtuplets, std::vector< int > &seedIndices, const unsigned int minClustersPerNtuplet, const bool outInDFS, const unsigned int maxOutInHops)
std::vector< HGCDoublet > &	getAllDoublets ()
void	makeAndConnectDoublets (const TICLLayerTiles &h, const std::vector< TICLSeedingRegion > &regions, int nEtaBins, int nPhiBins, const std::vector< reco::CaloCluster > &layerClusters, const std::vector< float > &mask, const edm::ValueMap< float > &layerClustersTime, int deltaIEta, int deltaIPhi, float minCosThetai, float maxCosPointing, int missing_layers, int maxNumberOfLayers, float maxDeltaTime)
void	setVerbosity (int level)

Private Attributes
std::vector< HGCDoublet >	allDoublets_
std::vector< std::vector< int > >	isOuterClusterOfDoublets_
std::vector< std::vector< HGCDoublet * > >	theNtuplets_
std::vector< unsigned int >	theRootDoublets_
int	verbosity_

Detailed Description

Definition at line 14 of file HGCGraph.h.

Member Enumeration Documentation

enum HGCGraph::VerbosityLevel

Enumerator
None
Basic
Advanced
Expert
Guru

Definition at line 46 of file HGCGraph.h.

{ None = 0, Basic, Advanced, Expert, Guru };

Member Function Documentation

bool HGCGraph::areTimeCompatible	(	int	innerIdx,
		int	outerIdx,
		const edm::ValueMap< float > &	layerClustersTime,
		float	maxDeltaTime
	)

Definition at line 126 of file HGCGraph.cc.

References funct::abs(), and edm::ValueMap< T >::get().

Referenced by makeAndConnectDoublets().

                                                     {
  float timeIn = layerClustersTime.get(innerIdx);
  float timeOut = layerClustersTime.get(outerIdx);

  return (timeIn == -99 || timeOut == -99 || std::abs(timeIn - timeOut) < maxDeltaTime);
}

void HGCGraph::clear ( void  )

inline

Definition at line 39 of file HGCGraph.h.

References allDoublets_, isOuterClusterOfDoublets_, theNtuplets_, and theRootDoublets_.

               {
    allDoublets_.clear();
    theRootDoublets_.clear();
    theNtuplets_.clear();
    isOuterClusterOfDoublets_.clear();
  }

void HGCGraph::findNtuplets	(	std::vector< HGCDoublet::HGCntuplet > &	foundNtuplets,
		std::vector< int > &	seedIndices,
		const unsigned int	minClustersPerNtuplet,
		const bool	outInDFS,
		const unsigned int	maxOutInHops
	)

Definition at line 137 of file HGCGraph.cc.

References allDoublets_, and theRootDoublets_.

Referenced by getAllDoublets().

                                                       {
  HGCDoublet::HGCntuplet tmpNtuplet;
  tmpNtuplet.reserve(minClustersPerNtuplet);
  std::vector<std::pair<unsigned int, unsigned int> > outInToVisit;
  for (auto rootDoublet : theRootDoublets_) {
    tmpNtuplet.clear();
    outInToVisit.clear();
    int seedIndex = allDoublets_[rootDoublet].seedIndex();
    int outInHops = 0;
    allDoublets_[rootDoublet].findNtuplets(
        allDoublets_, tmpNtuplet, seedIndex, outInDFS, outInHops, maxOutInHops, outInToVisit);
    while (!outInToVisit.empty()) {
      allDoublets_[outInToVisit.back().first].findNtuplets(
          allDoublets_, tmpNtuplet, seedIndex, outInDFS, outInToVisit.back().second, maxOutInHops, outInToVisit);
      outInToVisit.pop_back();
    }

    if (tmpNtuplet.size() > minClustersPerNtuplet) {
      foundNtuplets.push_back(tmpNtuplet);
      seedIndices.push_back(seedIndex);
    }
  }
}

std::vector<HGCDoublet>& HGCGraph::getAllDoublets ( )

inline

Definition at line 33 of file HGCGraph.h.

References allDoublets_, and findNtuplets().

{ return allDoublets_; }

void HGCGraph::makeAndConnectDoublets	(	const TICLLayerTiles &	h,
		const std::vector< TICLSeedingRegion > &	regions,
		int	nEtaBins,
		int	nPhiBins,
		const std::vector< reco::CaloCluster > &	layerClusters,
		const std::vector< float > &	mask,
		const edm::ValueMap< float > &	layerClustersTime,
		int	deltaIEta,
		int	deltaIPhi,
		float	minCosThetai,
		float	maxCosPointing,
		int	missing_layers,
		int	maxNumberOfLayers,
		float	maxDeltaTime
	)

Definition at line 10 of file HGCGraph.cc.

References Advanced, allDoublets_, areTimeCompatible(), Expert, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, muons_cff::isGlobal, isOuterClusterOfDoublets_, LogDebug, SiStripPI::max, min(), ticl::constants::nEtaBins, None, ticl::constants::nPhiBins, hltrates_dqm_sourceclient-live_cfg::offset, alignCSCRings::r, theRootDoublets_, and verbosity_.

                                                          {
  isOuterClusterOfDoublets_.clear();
  isOuterClusterOfDoublets_.resize(layerClusters.size());
  allDoublets_.clear();
  theRootDoublets_.clear();
  for (const auto &r : regions) {
    bool isGlobal = (r.index == -1);
    auto zSide = r.zSide;
    int startEtaBin, endEtaBin, startPhiBin, endPhiBin;

    if (isGlobal) {
      startEtaBin = 0;
      startPhiBin = 0;
      endEtaBin = nEtaBins;
      endPhiBin = nPhiBins;
    } else {
      auto firstLayerOnZSide = maxNumberOfLayers * zSide;
      const auto &firstLayerHisto = histo[firstLayerOnZSide];

      int entryEtaBin = firstLayerHisto.etaBin(r.origin.eta());
      int entryPhiBin = firstLayerHisto.phiBin(r.origin.phi());
      startEtaBin = std::max(entryEtaBin - deltaIEta, 0);
      endEtaBin = std::min(entryEtaBin + deltaIEta + 1, nEtaBins);
      startPhiBin = entryPhiBin - deltaIPhi;
      endPhiBin = entryPhiBin + deltaIPhi;
    }

    for (int il = 0; il < maxNumberOfLayers - 1; ++il) {
      for (int outer_layer = 0; outer_layer < std::min(1 + missing_layers, maxNumberOfLayers - 1 - il); ++outer_layer) {
        int currentInnerLayerId = il + maxNumberOfLayers * zSide;
        int currentOuterLayerId = currentInnerLayerId + 1 + outer_layer;
        auto const &outerLayerHisto = histo[currentOuterLayerId];
        auto const &innerLayerHisto = histo[currentInnerLayerId];

        for (int oeta = startEtaBin; oeta < endEtaBin; ++oeta) {
          auto offset = oeta * nPhiBins;
          for (int ophi_it = startPhiBin; ophi_it < endPhiBin; ++ophi_it) {
            int ophi = ((ophi_it % nPhiBins + nPhiBins) % nPhiBins);
            for (auto outerClusterId : outerLayerHisto[offset + ophi]) {
              // Skip masked clusters
              if (mask[outerClusterId] == 0.)
                continue;
              const auto etaRangeMin = std::max(0, oeta - deltaIEta);
              const auto etaRangeMax = std::min(oeta + deltaIEta + 1, nEtaBins);

              for (int ieta = etaRangeMin; ieta < etaRangeMax; ++ieta) {
                // wrap phi bin
                for (int phiRange = 0; phiRange < 2 * deltaIPhi + 1; ++phiRange) {
                  // The first wrapping is to take into account the
                  // cases in which we would have to seach in
                  // negative bins. The second wrap is mandatory to
                  // account for all other cases, since we add in
                  // between a full nPhiBins slot.
                  auto iphi = ((ophi + phiRange - deltaIPhi) % nPhiBins + nPhiBins) % nPhiBins;
                  for (auto innerClusterId : innerLayerHisto[ieta * nPhiBins + iphi]) {
                    // Skip masked clusters
                    if (mask[innerClusterId] == 0.)
                      continue;
                    auto doubletId = allDoublets_.size();
                    if (maxDeltaTime != -1 &&
                        !areTimeCompatible(innerClusterId, outerClusterId, layerClustersTime, maxDeltaTime))
                      continue;
                    allDoublets_.emplace_back(innerClusterId, outerClusterId, doubletId, &layerClusters, r.index);
                    if (verbosity_ > Advanced) {
                      LogDebug("HGCGraph")
                          << "Creating doubletsId: " << doubletId << " layerLink in-out: [" << currentInnerLayerId
                          << ", " << currentOuterLayerId << "] clusterLink in-out: [" << innerClusterId << ", "
                          << outerClusterId << "]" << std::endl;
                    }
                    isOuterClusterOfDoublets_[outerClusterId].push_back(doubletId);
                    auto &neigDoublets = isOuterClusterOfDoublets_[innerClusterId];
                    auto &thisDoublet = allDoublets_[doubletId];
                    if (verbosity_ > Expert) {
                      LogDebug("HGCGraph")
                          << "Checking compatibility of doubletId: " << doubletId
                          << " with all possible inners doublets link by the innerClusterId: " << innerClusterId
                          << std::endl;
                    }
                    bool isRootDoublet = thisDoublet.checkCompatibilityAndTag(allDoublets_,
                                                                              neigDoublets,
                                                                              r.directionAtOrigin,
                                                                              minCosTheta,
                                                                              minCosPointing,
                                                                              verbosity_ > Advanced);
                    if (isRootDoublet)
                      theRootDoublets_.push_back(doubletId);
                  }
                }
              }
            }
          }
        }
      }
    }
  }
  // #ifdef FP_DEBUG
  if (verbosity_ > None) {
    LogDebug("HGCGraph") << "number of Root doublets " << theRootDoublets_.size() << " over a total number of doublets "
                         << allDoublets_.size() << std::endl;
  }
  // #endif
}

void HGCGraph::setVerbosity ( int  level )

inline

Definition at line 45 of file HGCGraph.h.

References personalPlayback::level, and verbosity_.

{ verbosity_ = level; }

Member Data Documentation

std::vector<HGCDoublet> HGCGraph::allDoublets_

private

Definition at line 49 of file HGCGraph.h.

Referenced by clear(), findNtuplets(), getAllDoublets(), and makeAndConnectDoublets().

std::vector<std::vector<int> > HGCGraph::isOuterClusterOfDoublets_

private

Definition at line 52 of file HGCGraph.h.

Referenced by clear(), and makeAndConnectDoublets().

std::vector<std::vector<HGCDoublet *> > HGCGraph::theNtuplets_

private

Definition at line 51 of file HGCGraph.h.

Referenced by clear().

std::vector<unsigned int> HGCGraph::theRootDoublets_

private

Definition at line 50 of file HGCGraph.h.

Referenced by clear(), findNtuplets(), and makeAndConnectDoublets().

int HGCGraph::verbosity_

private

Definition at line 53 of file HGCGraph.h.

Referenced by makeAndConnectDoublets(), and setVerbosity().