ticl::PatternRecognitionbyCLUE3D< TILES > Class Template Reference

#include <PatternRecognitionbyCLUE3D.h>

Inheritance diagram for ticl::PatternRecognitionbyCLUE3D< TILES >:

Classes
struct	ClustersOnLayer

Public Member Functions
void	energyRegressionAndID (const std::vector< reco::CaloCluster > &layerClusters, std::vector< Trackster > &result)
void	makeTracksters (const typename PatternRecognitionAlgoBaseT< TILES >::Inputs &input, std::vector< Trackster > &result, std::unordered_map< int, std::vector< int >> &seedToTracksterAssociation) override
	PatternRecognitionbyCLUE3D (const edm::ParameterSet &conf, const CacheBase *cache, edm::ConsumesCollector)
	~PatternRecognitionbyCLUE3D () override=default
Public Member Functions inherited from ticl::PatternRecognitionAlgoBaseT< TILES >
virtual void	makeTracksters (const Inputs &input, std::vector< Trackster > &result, std::unordered_map< int, std::vector< int >> &seedToTracksterAssociation)=0
	PatternRecognitionAlgoBaseT (const edm::ParameterSet &conf, const CacheBase *cache, edm::ConsumesCollector)
virtual	~PatternRecognitionAlgoBaseT ()

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &iDesc)

Private Member Functions
void	calculateDistanceToHigher (const TILES &, const unsigned int layerId, const std::vector< std::pair< int, int >> &)
void	calculateLocalDensity (const TILES &, const unsigned int layerId, const std::vector< std::pair< int, int >> &)
void	dumpClusters (const std::vector< std::pair< int, int >> &layerIdx2layerandSoa, const int) const
void	dumpTiles (const TILES &) const
void	dumpTracksters (const std::vector< std::pair< int, int >> &layerIdx2layerandSoa, const int, const std::vector< Trackster > &) const
int	findAndAssignTracksters (const TILES &, const std::vector< std::pair< int, int >> &)
void	reset ()

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_
std::vector< ClustersOnLayer >	clusters_
const double	criticalDensity_
const double	criticalEtaPhiDistance_
const double	densityEtaPhiDistanceSqr_
const double	densityOnSameLayer_
const int	densitySiblingLayers_
const std::string	eidInputName_
const float	eidMinClusterEnergy_
const int	eidNClusters_
const int	eidNLayers_
const std::string	eidOutputNameEnergy_
const std::string	eidOutputNameId_
tensorflow::Session *	eidSession_
const std::vector< int >	filter_on_categories_
const int	minNumLayerCluster_
const double	outlierMultiplier_
hgcal::RecHitTools	rhtools_

Static Private Attributes
static const int	eidNFeatures_ = 3

Additional Inherited Members
Public Types inherited from ticl::PatternRecognitionAlgoBaseT< TILES >
enum	VerbosityLevel {   None = 0, Basic, Advanced, Expert,   Guru }
Protected Attributes inherited from ticl::PatternRecognitionAlgoBaseT< TILES >
int	algo_verbosity_

Detailed Description

template<typename TILES>
class ticl::PatternRecognitionbyCLUE3D< TILES >

Definition at line 12 of file PatternRecognitionbyCLUE3D.h.

Constructor & Destructor Documentation

template<typename TILES >

PatternRecognitionbyCLUE3D::PatternRecognitionbyCLUE3D	(	const edm::ParameterSet &	conf,
		const CacheBase *	cache,
		edm::ConsumesCollector	iC
	)

Definition at line 24 of file PatternRecognitionbyCLUE3D.cc.

References utilities::cache(), tensorflow::createSession(), ticl::TrackstersCache::eidGraphDef, ticl::PatternRecognitionbyCLUE3D< TILES >::eidSession_, and Exception.

    : PatternRecognitionAlgoBaseT<TILES>(conf, cache, iC),
      caloGeomToken_(iC.esConsumes<CaloGeometry, CaloGeometryRecord>()),
      criticalDensity_(conf.getParameter<double>("criticalDensity")),
      densitySiblingLayers_(conf.getParameter<int>("densitySiblingLayers")),
      densityEtaPhiDistanceSqr_(conf.getParameter<double>("densityEtaPhiDistanceSqr")),
      densityOnSameLayer_(conf.getParameter<bool>("densityOnSameLayer")),
      criticalEtaPhiDistance_(conf.getParameter<double>("criticalEtaPhiDistance")),
      outlierMultiplier_(conf.getParameter<double>("outlierMultiplier")),
      minNumLayerCluster_(conf.getParameter<int>("minNumLayerCluster")),
      eidInputName_(conf.getParameter<std::string>("eid_input_name")),
      eidOutputNameEnergy_(conf.getParameter<std::string>("eid_output_name_energy")),
      eidOutputNameId_(conf.getParameter<std::string>("eid_output_name_id")),
      eidMinClusterEnergy_(conf.getParameter<double>("eid_min_cluster_energy")),
      eidNLayers_(conf.getParameter<int>("eid_n_layers")),
      eidNClusters_(conf.getParameter<int>("eid_n_clusters")),
      eidSession_(nullptr) {
  // mount the tensorflow graph onto the session when set
  const TrackstersCache *trackstersCache = dynamic_cast<const TrackstersCache *>(cache);
  if (trackstersCache == nullptr || trackstersCache->eidGraphDef == nullptr) {
    throw cms::Exception("MissingGraphDef")
        << "PatternRecognitionbyCLUE3D received an empty graph definition from the global cache";
  }
  eidSession_ = tensorflow::createSession(trackstersCache->eidGraphDef);
}

template<typename TILES >

ticl::PatternRecognitionbyCLUE3D< TILES >::~PatternRecognitionbyCLUE3D ( )

overridedefault

Member Function Documentation

template<typename TILES >

void PatternRecognitionbyCLUE3D::calculateDistanceToHigher ( const TILES &  tiles, const unsigned int  layerId, const std::vector< std::pair< int, int >> &  layerIdx2layerandSoa  )

private

Definition at line 561 of file PatternRecognitionbyCLUE3D.cc.

References reco::deltaR2(), mps_fire::i, SiStripPI::max, min(), HLT_FULL_cff::minLayer, L1TMuonDQMOffline_cfi::nEtaBins, ecaldqm::binning::nPhiBins, hltrates_dqm_sourceclient-live_cfg::offset, and mathSSE::sqrt().

                                                                                                              {
  constexpr int nEtaBin = TILES::constants_type_t::nEtaBins;
  constexpr int nPhiBin = TILES::constants_type_t::nPhiBins;
  auto &clustersOnLayer = clusters_[layerId];
  unsigned int numberOfClusters = clustersOnLayer.x.size();

  for (unsigned int i = 0; i < numberOfClusters; i++) {
    if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
      edm::LogVerbatim("PatternRecogntionbyCLUE3D")
          << "Starting searching nearestHigher on " << layerId << " with rho: " << clustersOnLayer.rho[i]
          << " at eta, phi: " << tiles[layerId].etaBin(clustersOnLayer.eta[i]) << ", "
          << tiles[layerId].etaBin(clustersOnLayer.phi[i]);
    }
    // We need to partition the two sides of the HGCAL detector
    auto lastLayerPerSide = static_cast<unsigned int>(rhtools_.lastLayer(false));
    unsigned int minLayer = 0;
    unsigned int maxLayer = 2 * lastLayerPerSide - 1;
    if (layerId < lastLayerPerSide) {
      minLayer = std::max(layerId - densitySiblingLayers_, minLayer);
      maxLayer = std::min(layerId + densitySiblingLayers_, lastLayerPerSide - 1);
    } else {
      minLayer = std::max(layerId - densitySiblingLayers_, lastLayerPerSide + 1);
      maxLayer = std::min(layerId + densitySiblingLayers_, maxLayer);
    }
    float maxDelta = std::numeric_limits<float>::max();
    float i_delta = maxDelta;
    std::pair<int, int> i_nearestHigher(-1, -1);
    for (unsigned int currentLayer = minLayer; currentLayer <= maxLayer; currentLayer++) {
      const auto &tileOnLayer = tiles[currentLayer];
      int etaWindow = 2;
      int phiWindow = 2;
      int etaBinMin = std::max(tileOnLayer.etaBin(clustersOnLayer.eta[i]) - etaWindow, 0);
      int etaBinMax = std::min(tileOnLayer.etaBin(clustersOnLayer.eta[i]) + etaWindow, nEtaBin);
      int phiBinMin = tileOnLayer.phiBin(clustersOnLayer.phi[i]) - phiWindow;
      int phiBinMax = tileOnLayer.phiBin(clustersOnLayer.phi[i]) + phiWindow;
      for (int ieta = etaBinMin; ieta <= etaBinMax; ++ieta) {
        auto offset = ieta * nPhiBin;
        for (int iphi_it = phiBinMin; iphi_it <= phiBinMax; ++iphi_it) {
          int iphi = ((iphi_it % nPhiBin + nPhiBin) % nPhiBin);
          if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
            edm::LogVerbatim("PatternRecogntionbyCLUE3D")
                << "Searching nearestHigher on " << currentLayer << " eta, phi: " << ieta << ", " << iphi_it;
          }
          for (auto otherClusterIdx : tileOnLayer[offset + iphi]) {
            auto const &layerandSoa = layerIdx2layerandSoa[otherClusterIdx];
            // Skip masked layer clusters
            if ((layerandSoa.first == -1) && (layerandSoa.second == -1))
              continue;
            auto const &clustersOnOtherLayer = clusters_[layerandSoa.first];
            float dist = reco::deltaR2(clustersOnLayer.eta[i],
                                       clustersOnLayer.phi[i],
                                       clustersOnOtherLayer.eta[layerandSoa.second],
                                       clustersOnOtherLayer.phi[layerandSoa.second]);
            bool foundHigher = (clustersOnOtherLayer.rho[layerandSoa.second] > clustersOnLayer.rho[i]) ||
                               (clustersOnOtherLayer.rho[layerandSoa.second] == clustersOnLayer.rho[i] &&
                                clustersOnOtherLayer.layerClusterOriginalIdx[layerandSoa.second] >
                                    clustersOnLayer.layerClusterOriginalIdx[i]);
            if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
              edm::LogVerbatim("PatternRecogntionbyCLUE3D")
                  << "Searching nearestHigher on " << currentLayer
                  << " with rho: " << clustersOnOtherLayer.rho[layerandSoa.second]
                  << " on layerIdxInSOA: " << layerandSoa.first << ", " << layerandSoa.second
                  << " with distance: " << sqrt(dist) << " foundHigher: " << foundHigher;
            }
            if (foundHigher && dist <= i_delta) {
              // update i_delta
              i_delta = dist;
              // update i_nearestHigher
              i_nearestHigher = layerandSoa;
            }
          }
        }
      }
    }

    bool foundNearestHigherInEtaPhiCylinder = (i_delta != maxDelta);
    if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
      edm::LogVerbatim("PatternRecogntionbyCLUE3D")
          << "i_delta: " << sqrt(i_delta) << " passed: " << foundNearestHigherInEtaPhiCylinder << " "
          << i_nearestHigher.first << " " << i_nearestHigher.second;
    }
    if (foundNearestHigherInEtaPhiCylinder) {
      clustersOnLayer.delta[i] = sqrt(i_delta);
      clustersOnLayer.nearestHigher[i] = i_nearestHigher;
    } else {
      // otherwise delta is guaranteed to be larger outlierDeltaFactor_*delta_c
      // we can safely maximize delta to be maxDelta
      clustersOnLayer.delta[i] = maxDelta;
      clustersOnLayer.nearestHigher[i] = {-1, -1};
    }
  }
}

template<typename TILES >

void PatternRecognitionbyCLUE3D::calculateLocalDensity ( const TILES &  tiles, const unsigned int  layerId, const std::vector< std::pair< int, int >> &  layerIdx2layerandSoa  )

private

Definition at line 442 of file PatternRecognitionbyCLUE3D.cc.

References CommonMethods::delta(), reco::deltaR2(), mps_fire::i, SiStripPI::max, min(), HLT_FULL_cff::minLayer, L1TMuonDQMOffline_cfi::nEtaBins, ecaldqm::binning::nPhiBins, and hltrates_dqm_sourceclient-live_cfg::offset.

                                                                                                              {
  constexpr int nEtaBin = TILES::constants_type_t::nEtaBins;
  constexpr int nPhiBin = TILES::constants_type_t::nPhiBins;
  auto &clustersOnLayer = clusters_[layerId];
  unsigned int numberOfClusters = clustersOnLayer.x.size();

  auto isReachable = [&](float x1, float x2, float y1, float y2, float delta_sqr) -> bool {
    if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
      edm::LogVerbatim("PatternRecogntionbyCLUE3D")
          << ((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2)) << " vs " << delta_sqr << "["
          << ((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2) < delta_sqr) << "]\n";
    }
    return (x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2) < delta_sqr;
  };

  for (unsigned int i = 0; i < numberOfClusters; i++) {
    // We need to partition the two sides of the HGCAL detector
    auto lastLayerPerSide = static_cast<unsigned int>(rhtools_.lastLayer(false)) - 1;
    unsigned int minLayer = 0;
    unsigned int maxLayer = 2 * lastLayerPerSide - 1;
    if (layerId < lastLayerPerSide) {
      minLayer = std::max(layerId - densitySiblingLayers_, minLayer);
      maxLayer = std::min(layerId + densitySiblingLayers_, lastLayerPerSide - 1);
    } else {
      minLayer = std::max(layerId - densitySiblingLayers_, lastLayerPerSide + 1);
      maxLayer = std::min(layerId + densitySiblingLayers_, maxLayer);
    }
    for (unsigned int currentLayer = minLayer; currentLayer <= maxLayer; currentLayer++) {
      if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "RefLayer: " << layerId << " SoaIDX: " << i;
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "NextLayer: " << currentLayer;
      }
      const auto &tileOnLayer = tiles[currentLayer];
      bool onSameLayer = (currentLayer == layerId);
      if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "onSameLayer: " << onSameLayer;
      }
      int etaWindow = onSameLayer ? 2 : 1;
      int phiWindow = onSameLayer ? 2 : 1;
      int etaBinMin = std::max(tileOnLayer.etaBin(clustersOnLayer.eta[i]) - etaWindow, 0);
      int etaBinMax = std::min(tileOnLayer.etaBin(clustersOnLayer.eta[i]) + etaWindow, nEtaBin);
      int phiBinMin = tileOnLayer.phiBin(clustersOnLayer.phi[i]) - phiWindow;
      int phiBinMax = tileOnLayer.phiBin(clustersOnLayer.phi[i]) + phiWindow;
      if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "eta: " << clustersOnLayer.eta[i];
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "phi: " << clustersOnLayer.phi[i];
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "etaBinMin: " << etaBinMin << ", etaBinMax: " << etaBinMax;
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "phiBinMin: " << phiBinMin << ", phiBinMax: " << phiBinMax;
      }
      for (int ieta = etaBinMin; ieta <= etaBinMax; ++ieta) {
        auto offset = ieta * nPhiBin;
        if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
          edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "offset: " << offset;
        }
        for (int iphi_it = phiBinMin; iphi_it <= phiBinMax; ++iphi_it) {
          int iphi = ((iphi_it % nPhiBin + nPhiBin) % nPhiBin);
          if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
            edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "iphi: " << iphi;
            edm::LogVerbatim("PatternRecogntionbyCLUE3D")
                << "Entries in tileBin: " << tileOnLayer[offset + iphi].size();
          }
          for (auto otherClusterIdx : tileOnLayer[offset + iphi]) {
            auto const &layerandSoa = layerIdx2layerandSoa[otherClusterIdx];
            // Skip masked layer clusters
            if ((layerandSoa.first == -1) && (layerandSoa.second == -1)) {
              if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
                edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Skipping masked layerIdx " << otherClusterIdx;
              }
              continue;
            }
            auto const &clustersLayer = clusters_[layerandSoa.first];
            if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
              edm::LogVerbatim("PatternRecogntionbyCLUE3D")
                  << "OtherLayer: " << layerandSoa.first << " SoaIDX: " << layerandSoa.second;
              edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "OtherEta: " << clustersLayer.eta[layerandSoa.second];
              edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "OtherPhi: " << clustersLayer.phi[layerandSoa.second];
            }
            // extend by 26 mm, roughly 2 cells, more wrt sum of radii
            float delta = clustersOnLayer.radius[i] + clustersLayer.radius[layerandSoa.second] + 2.6f;
            if (densityOnSameLayer_ && onSameLayer) {
              if (isReachable(clustersOnLayer.x[i],
                              clustersLayer.x[layerandSoa.second],
                              clustersOnLayer.y[i],
                              clustersLayer.y[layerandSoa.second],
                              delta * delta)) {
                clustersOnLayer.rho[i] += (clustersOnLayer.layerClusterOriginalIdx[i] == otherClusterIdx ? 1.f : 0.2f) *
                                          clustersLayer.energy[layerandSoa.second];
              }
            } else {
              if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
                edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Distance: "
                                                              << reco::deltaR2(clustersOnLayer.eta[i],
                                                                               clustersOnLayer.phi[i],
                                                                               clustersLayer.eta[layerandSoa.second],
                                                                               clustersLayer.phi[layerandSoa.second]);
              }
              if (reco::deltaR2(clustersOnLayer.eta[i],
                                clustersOnLayer.phi[i],
                                clustersLayer.eta[layerandSoa.second],
                                clustersLayer.phi[layerandSoa.second]) < densityEtaPhiDistanceSqr_) {
                auto energyToAdd = (clustersOnLayer.layerClusterOriginalIdx[i] == otherClusterIdx ? 1.f : 0.5f) *
                                   clustersLayer.energy[layerandSoa.second];
                clustersOnLayer.rho[i] += energyToAdd;
                edm::LogVerbatim("PatternRecogntionbyCLUE3D")
                    << "Adding " << energyToAdd << " partial " << clustersOnLayer.rho[i];
              }
            }
          }  // end of loop on possible compatible clusters
        }    // end of loop over phi-bin region
      }      // end of loop over eta-bin region
    }        // end of loop on the sibling layers
  }          // end of loop over clusters on this layer
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    edm::LogVerbatim("PatternRecogntionbyCLUE3D") << std::endl;
  }
}

template<typename TILES >

void PatternRecognitionbyCLUE3D::dumpClusters ( const std::vector< std::pair< int, int >> &  layerIdx2layerandSoa, const int  eventNumber  ) const

private

Definition at line 107 of file PatternRecognitionbyCLUE3D.cc.

References phase1PixelTopology::layer, pileupDistInMC::num, and findQualityFiles::v.

                                                                                  {
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "[evt, layer, x, y, eta, phi, cells, energy, radius, rho, delta, "
                                                     "isSeed, clusterIdx, layerClusterOriginalIdx";
  }

  for (unsigned int layer = 0; layer < clusters_.size(); layer++) {
    if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
      edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "On Layer " << layer;
    }
    auto const &thisLayer = clusters_[layer];
    int num = 0;
    for (auto v : thisLayer.x) {
      if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
        edm::LogVerbatim("PatternRecogntionbyCLUE3D")
            << "ClusterInfo: " << eventNumber << ", " << layer << ", " << v << ", " << thisLayer.y[num] << ", "
            << thisLayer.eta[num] << ", " << thisLayer.phi[num] << ", " << thisLayer.cells[num] << ", "
            << thisLayer.energy[num] << ", " << thisLayer.radius[num] << ", " << thisLayer.rho[num] << ", "
            << thisLayer.delta[num] << ", " << thisLayer.isSeed[num] << ", " << thisLayer.clusterIndex[num] << ", "
            << thisLayer.layerClusterOriginalIdx[num];
      }
      ++num;
    }
  }
  for (unsigned int lcIdx = 0; lcIdx < layerIdx2layerandSoa.size(); lcIdx++) {
    auto const &layerandSoa = layerIdx2layerandSoa[lcIdx];
    // Skip masked layer clusters
    if ((layerandSoa.first == -1) && (layerandSoa.second == -1))
      continue;
    if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
      edm::LogVerbatim("PatternRecogntionbyCLUE3D")
          << "lcIdx: " << lcIdx << " on Layer: " << layerandSoa.first << " SOA: " << layerandSoa.second;
    }
  }
}

template<typename TILES >

void PatternRecognitionbyCLUE3D::dumpTiles ( const TILES &  tiles ) const

private

Definition at line 53 of file PatternRecognitionbyCLUE3D.cc.

References relativeConstraints::empty, phase1PixelTopology::layer, L1TMuonDQMOffline_cfi::nEtaBins, ecaldqm::binning::nPhiBins, and hltrates_dqm_sourceclient-live_cfg::offset.

                                                                          {
  constexpr int nEtaBin = TILES::constants_type_t::nEtaBins;
  constexpr int nPhiBin = TILES::constants_type_t::nPhiBins;
  auto lastLayerPerSide = (unsigned int)(rhtools_.lastLayer(false));
  unsigned int maxLayer = 2 * lastLayerPerSide - 1;
  for (unsigned int layer = 0; layer <= maxLayer; layer++) {
    for (int ieta = 0; ieta < nEtaBin; ieta++) {
      auto offset = ieta * nPhiBin;
      for (int phi = 0; phi < nPhiBin; phi++) {
        int iphi = ((phi % nPhiBin + nPhiBin) % nPhiBin);
        if (!tiles[layer][offset + iphi].empty()) {
          if (this->algo_verbosity_ > this->Advanced) {
            edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Layer: " << layer << " ieta: " << ieta << " phi: " << phi
                                                          << " " << tiles[layer][offset + iphi].size();
          }
        }
      }
    }
  }
}

template<typename TILES >

void PatternRecognitionbyCLUE3D::dumpTracksters ( const std::vector< std::pair< int, int >> &  layerIdx2layerandSoa, const int  eventNumber, const std::vector< Trackster > &  tracksters  ) const

private

Definition at line 75 of file PatternRecognitionbyCLUE3D.cc.

References ticl::Trackster::charged_hadron, ticl::Trackster::electron, ticl::Trackster::neutral_hadron, pileupDistInMC::num, ticl::Trackster::photon, AlCaHLTBitMon_QueryRunRegistry::string, submitPVValidationJobs::t, and findQualityFiles::v.

                                                                                                     {
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    edm::LogVerbatim("PatternRecogntionbyCLUE3D")
        << "[evt, tracksterId, cells, prob_photon, prob_ele, prob_chad, prob_nhad, layer_i, x_i, y_i, eta_i, phi_i, "
           "energy_i, radius_i, rho_i, delta_i, isSeed_i";
  }

  int num = 0;
  const std::string sep(", ");
  for (auto const &t : tracksters) {
    for (auto v : t.vertices()) {
      auto [lyrIdx, soaIdx] = layerIdx2layerandSoa[v];
      auto const &thisLayer = clusters_[lyrIdx];
      if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
        edm::LogVerbatim("PatternRecogntionbyCLUE3D")
            << "TracksterInfo: " << eventNumber << sep << num << sep << t.vertices().size() << sep
            << t.id_probability(ticl::Trackster::ParticleType::photon) << sep
            << t.id_probability(ticl::Trackster::ParticleType::electron) << sep
            << t.id_probability(ticl::Trackster::ParticleType::charged_hadron) << sep
            << t.id_probability(ticl::Trackster::ParticleType::neutral_hadron) << sep << lyrIdx << sep
            << thisLayer.x[soaIdx] << sep << thisLayer.y[soaIdx] << sep << thisLayer.eta[soaIdx] << sep
            << thisLayer.phi[soaIdx] << sep << thisLayer.energy[soaIdx] << sep << thisLayer.radius[soaIdx] << sep
            << thisLayer.rho[soaIdx] << sep << thisLayer.delta[soaIdx] << sep << thisLayer.isSeed[soaIdx] << '\n';
      }
    }
    num++;
  }
}

template<typename TILES >

void PatternRecognitionbyCLUE3D::energyRegressionAndID	(	const std::vector< reco::CaloCluster > &	layerClusters,
		std::vector< Trackster > &	result
	)

Definition at line 303 of file PatternRecognitionbyCLUE3D.cc.

References a, funct::abs(), b, data, relval_parameters_module::energy, reco::CaloCluster::energy(), reco::CaloCluster::eta(), validate-o2o-wbm::f, lowptgsfeleseed::features(), reco::CaloCluster::hitsAndFractions(), mps_fire::i, input, dqmiolumiharvest::j, isotrackApplyRegressor::k, cmsLHEtoEOSManager::l, reco::CaloCluster::phi(), tensorflow::run(), and beam_dqm_sourceclient-live_cfg::vertices.

                                                                                                {
  // Energy regression and particle identification strategy:
  //
  // 1. Set default values for regressed energy and particle id for each trackster.
  // 2. Store indices of tracksters whose total sum of cluster energies is above the
  //    eidMinClusterEnergy_ (GeV) treshold. Inference is not applied for soft tracksters.
  // 3. When no trackster passes the selection, return.
  // 4. Create input and output tensors. The batch dimension is determined by the number of
  //    selected tracksters.
  // 5. Fill input tensors with layer cluster features. Per layer, clusters are ordered descending
  //    by energy. Given that tensor data is contiguous in memory, we can use pointer arithmetic to
  //    fill values, even with batching.
  // 6. Zero-fill features for empty clusters in each layer.
  // 7. Batched inference.
  // 8. Assign the regressed energy and id probabilities to each trackster.
  //
  // Indices used throughout this method:
  // i -> batch element / trackster
  // j -> layer
  // k -> cluster
  // l -> feature

  // set default values per trackster, determine if the cluster energy threshold is passed,
  // and store indices of hard tracksters
  std::vector<int> tracksterIndices;
  for (int i = 0; i < static_cast<int>(tracksters.size()); i++) {
    // calculate the cluster energy sum (2)
    // note: after the loop, sumClusterEnergy might be just above the threshold which is enough to
    // decide whether to run inference for the trackster or not
    float sumClusterEnergy = 0.;
    for (const unsigned int &vertex : tracksters[i].vertices()) {
      sumClusterEnergy += static_cast<float>(layerClusters[vertex].energy());
      // there might be many clusters, so try to stop early
      if (sumClusterEnergy >= eidMinClusterEnergy_) {
        // set default values (1)
        tracksters[i].setRegressedEnergy(0.f);
        tracksters[i].zeroProbabilities();
        tracksterIndices.push_back(i);
        break;
      }
    }
  }

  // do nothing when no trackster passes the selection (3)
  int batchSize = static_cast<int>(tracksterIndices.size());
  if (batchSize == 0) {
    return;
  }

  // create input and output tensors (4)
  tensorflow::TensorShape shape({batchSize, eidNLayers_, eidNClusters_, eidNFeatures_});
  tensorflow::Tensor input(tensorflow::DT_FLOAT, shape);
  tensorflow::NamedTensorList inputList = {{eidInputName_, input}};

  std::vector<tensorflow::Tensor> outputs;
  std::vector<std::string> outputNames;
  if (!eidOutputNameEnergy_.empty()) {
    outputNames.push_back(eidOutputNameEnergy_);
  }
  if (!eidOutputNameId_.empty()) {
    outputNames.push_back(eidOutputNameId_);
  }

  // fill input tensor (5)
  for (int i = 0; i < batchSize; i++) {
    const Trackster &trackster = tracksters[tracksterIndices[i]];

    // per layer, we only consider the first eidNClusters_ clusters in terms of energy, so in order
    // to avoid creating large / nested structures to do the sorting for an unknown number of total
    // clusters, create a sorted list of layer cluster indices to keep track of the filled clusters
    std::vector<int> clusterIndices(trackster.vertices().size());
    for (int k = 0; k < (int)trackster.vertices().size(); k++) {
      clusterIndices[k] = k;
    }
    sort(clusterIndices.begin(), clusterIndices.end(), [&layerClusters, &trackster](const int &a, const int &b) {
      return layerClusters[trackster.vertices(a)].energy() > layerClusters[trackster.vertices(b)].energy();
    });

    // keep track of the number of seen clusters per layer
    std::vector<int> seenClusters(eidNLayers_);

    // loop through clusters by descending energy
    for (const int &k : clusterIndices) {
      // get features per layer and cluster and store the values directly in the input tensor
      const reco::CaloCluster &cluster = layerClusters[trackster.vertices(k)];
      int j = rhtools_.getLayerWithOffset(cluster.hitsAndFractions()[0].first) - 1;
      if (j < eidNLayers_ && seenClusters[j] < eidNClusters_) {
        // get the pointer to the first feature value for the current batch, layer and cluster
        float *features = &input.tensor<float, 4>()(i, j, seenClusters[j], 0);

        // fill features
        *(features++) = float(cluster.energy() / float(trackster.vertex_multiplicity(k)));
        *(features++) = float(std::abs(cluster.eta()));
        *(features) = float(cluster.phi());

        // increment seen clusters
        seenClusters[j]++;
      }
    }

    // zero-fill features of empty clusters in each layer (6)
    for (int j = 0; j < eidNLayers_; j++) {
      for (int k = seenClusters[j]; k < eidNClusters_; k++) {
        float *features = &input.tensor<float, 4>()(i, j, k, 0);
        for (int l = 0; l < eidNFeatures_; l++) {
          *(features++) = 0.f;
        }
      }
    }
  }

  // run the inference (7)
  tensorflow::run(eidSession_, inputList, outputNames, &outputs);

  // store regressed energy per trackster (8)
  if (!eidOutputNameEnergy_.empty()) {
    // get the pointer to the energy tensor, dimension is batch x 1
    float *energy = outputs[0].flat<float>().data();

    for (const int &i : tracksterIndices) {
      tracksters[i].setRegressedEnergy(*(energy++));
    }
  }

  // store id probabilities per trackster (8)
  if (!eidOutputNameId_.empty()) {
    // get the pointer to the id probability tensor, dimension is batch x id_probabilities.size()
    int probsIdx = eidOutputNameEnergy_.empty() ? 0 : 1;
    float *probs = outputs[probsIdx].flat<float>().data();

    for (const int &i : tracksterIndices) {
      tracksters[i].setProbabilities(probs);
      probs += tracksters[i].id_probabilities().size();
    }
  }
}

template<typename TILES >

void PatternRecognitionbyCLUE3D::fillPSetDescription ( edm::ParameterSetDescription &  iDesc )

static

Definition at line 716 of file PatternRecognitionbyCLUE3D.cc.

References edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                             {
  iDesc.add<int>("algo_verbosity", 0);
  iDesc.add<double>("criticalDensity", 4)->setComment("in GeV");
  iDesc.add<int>("densitySiblingLayers", 3);
  iDesc.add<double>("densityEtaPhiDistanceSqr", 0.0008);
  iDesc.add<bool>("densityOnSameLayer", false);
  iDesc.add<double>("criticalEtaPhiDistance", 0.035);
  iDesc.add<double>("outlierMultiplier", 2);
  iDesc.add<int>("minNumLayerCluster", 5)->setComment("Not Inclusive");
  iDesc.add<std::string>("eid_input_name", "input");
  iDesc.add<std::string>("eid_output_name_energy", "output/regressed_energy");
  iDesc.add<std::string>("eid_output_name_id", "output/id_probabilities");
  iDesc.add<double>("eid_min_cluster_energy", 1.);
  iDesc.add<int>("eid_n_layers", 50);
  iDesc.add<int>("eid_n_clusters", 10);
}

template<typename TILES >

int PatternRecognitionbyCLUE3D::findAndAssignTracksters ( const TILES &  tiles, const std::vector< std::pair< int, int >> &  layerIdx2layerandSoa  )

private

Definition at line 656 of file PatternRecognitionbyCLUE3D.cc.

References mps_fire::i, and phase1PixelTopology::layer.

                                                                                  {
  unsigned int nTracksters = 0;

  std::vector<std::pair<int, int>> localStack;
  // find cluster seeds and outlier
  for (unsigned int layer = 0; layer < 2 * rhtools_.lastLayer(); layer++) {
    auto &clustersOnLayer = clusters_[layer];
    unsigned int numberOfClusters = clustersOnLayer.x.size();
    for (unsigned int i = 0; i < numberOfClusters; i++) {
      // initialize clusterIndex
      clustersOnLayer.clusterIndex[i] = -1;
      bool isSeed =
          (clustersOnLayer.delta[i] > criticalEtaPhiDistance_) && (clustersOnLayer.rho[i] >= criticalDensity_);
      bool isOutlier = (clustersOnLayer.delta[i] > outlierMultiplier_ * criticalEtaPhiDistance_) &&
                       (clustersOnLayer.rho[i] < criticalDensity_);
      if (isSeed) {
        if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
          edm::LogVerbatim("PatternRecogntionbyCLUE3D")
              << "Found seed on Layer " << layer << " SOAidx: " << i << " assigned ClusterIdx: " << nTracksters;
        }
        clustersOnLayer.clusterIndex[i] = nTracksters++;
        clustersOnLayer.isSeed[i] = true;
        localStack.emplace_back(layer, i);
      } else if (!isOutlier) {
        auto [lyrIdx, soaIdx] = clustersOnLayer.nearestHigher[i];
        if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
          edm::LogVerbatim("PatternRecogntionbyCLUE3D")
              << "Found follower on Layer " << layer << " SOAidx: " << i << " attached to cluster on layer: " << lyrIdx
              << " SOAidx: " << soaIdx;
        }
        clusters_[lyrIdx].followers[soaIdx].emplace_back(layer, i);
      } else {
        if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
          edm::LogVerbatim("PatternRecogntionbyCLUE3D")
              << "Found Outlier on Layer " << layer << " SOAidx: " << i << " with rho: " << clustersOnLayer.rho[i]
              << " and delta: " << clustersOnLayer.delta[i];
        }
      }
    }
  }

  // Propagate cluster index
  while (!localStack.empty()) {
    auto [lyrIdx, soaIdx] = localStack.back();
    auto &thisSeed = clusters_[lyrIdx].followers[soaIdx];
    localStack.pop_back();

    // loop over followers
    for (auto [follower_lyrIdx, follower_soaIdx] : thisSeed) {
      // pass id to a follower
      clusters_[follower_lyrIdx].clusterIndex[follower_soaIdx] = clusters_[lyrIdx].clusterIndex[soaIdx];
      // push this follower to localStack
      localStack.emplace_back(follower_lyrIdx, follower_soaIdx);
    }
  }
  return nTracksters;
}

template<typename TILES >

void PatternRecognitionbyCLUE3D::makeTracksters ( const typename PatternRecognitionAlgoBaseT< TILES >::Inputs &  input, std::vector< Trackster > &  result, std::unordered_map< int, std::vector< int >> &  seedToTracksterAssociation  )

override

Definition at line 145 of file PatternRecognitionbyCLUE3D.cc.

References cms::cuda::assert(), ticl::assignPCAtoTracksters(), SplitLinear::begin, caloGeomToken_, dataset::end, ticl::PatternRecognitionAlgoBaseT< TILES >::Inputs::es, ticl::PatternRecognitionAlgoBaseT< TILES >::Inputs::ev, edm::EventAuxiliary::event(), edm::Event::eventAuxiliary(), validate-o2o-wbm::f, relativeConstraints::geom, edm::EventSetup::getData(), mps_fire::i, phase1PixelTopology::layer, ticl::PatternRecognitionAlgoBaseT< TILES >::Inputs::layerClusters, ticl::PatternRecognitionAlgoBaseT< TILES >::Inputs::layerClustersTime, ticl::PatternRecognitionAlgoBaseT< TILES >::Inputs::mask, SiStripPI::max, point, ticl::PatternRecognitionAlgoBaseT< TILES >::Inputs::regions, reset(), mathSSE::sqrt(), makeMuonMisalignmentScenario::sum_x, makeMuonMisalignmentScenario::sum_y, submitPVValidationJobs::t, ticl::PatternRecognitionAlgoBaseT< TILES >::Inputs::tiles, and findQualityFiles::v.

                                                                       {
  // Protect from events with no seeding regions
  if (input.regions.empty())
    return;

  const int eventNumber = input.ev.eventAuxiliary().event();
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "New Event";
  }

  edm::EventSetup const &es = input.es;
  const CaloGeometry &geom = es.getData(caloGeomToken_);
  rhtools_.setGeometry(geom);

  clusters_.clear();
  clusters_.resize(2 * rhtools_.lastLayer(false));
  std::vector<std::pair<int, int>> layerIdx2layerandSoa;  //used everywhere also to propagate cluster masking

  layerIdx2layerandSoa.reserve(input.layerClusters.size());
  unsigned int layerIdx = 0;
  for (auto const &lc : input.layerClusters) {
    if (input.mask[layerIdx] == 0.) {
      if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Skipping masked clustrer: " << layerIdx;
      }
      layerIdx2layerandSoa.emplace_back(-1, -1);
      layerIdx++;
      continue;
    }
    const auto firstHitDetId = lc.hitsAndFractions()[0].first;
    int layer = rhtools_.getLayerWithOffset(firstHitDetId) - 1 +
                rhtools_.lastLayer(false) * ((rhtools_.zside(firstHitDetId) + 1) >> 1);
    assert(layer >= 0);

    layerIdx2layerandSoa.emplace_back(layer, clusters_[layer].x.size());
    float sum_x = 0.;
    float sum_y = 0.;
    float sum_sqr_x = 0.;
    float sum_sqr_y = 0.;
    float ref_x = lc.x();
    float ref_y = lc.y();
    float invClsize = 1. / lc.hitsAndFractions().size();
    for (auto const &hitsAndFractions : lc.hitsAndFractions()) {
      auto const &point = rhtools_.getPosition(hitsAndFractions.first);
      sum_x += point.x() - ref_x;
      sum_sqr_x += (point.x() - ref_x) * (point.x() - ref_x);
      sum_y += point.y() - ref_y;
      sum_sqr_y += (point.y() - ref_y) * (point.y() - ref_y);
    }
    // The variance of X for X uniform in circle of radius R is R^2/4,
    // therefore we multiply the sqrt(var) by 2 to have a rough estimate of the
    // radius. On the other hand, while averaging the x and y radius, we would
    // end up dividing by 2. Hence we omit the value here and in the average
    // below, too.
    float radius_x = sqrt((sum_sqr_x - (sum_x * sum_x) * invClsize) * invClsize);
    float radius_y = sqrt((sum_sqr_y - (sum_y * sum_y) * invClsize) * invClsize);
    clusters_[layer].x.emplace_back(lc.x());
    clusters_[layer].y.emplace_back(lc.y());
    clusters_[layer].radius.emplace_back(radius_x + radius_y);
    clusters_[layer].eta.emplace_back(lc.eta());
    clusters_[layer].phi.emplace_back(lc.phi());
    clusters_[layer].cells.push_back(lc.hitsAndFractions().size());
    clusters_[layer].energy.emplace_back(lc.energy());
    clusters_[layer].isSeed.push_back(false);
    clusters_[layer].clusterIndex.emplace_back(-1);
    clusters_[layer].layerClusterOriginalIdx.emplace_back(layerIdx++);
    clusters_[layer].nearestHigher.emplace_back(-1, -1);
    clusters_[layer].rho.emplace_back(0.f);
    clusters_[layer].delta.emplace_back(std::numeric_limits<float>::max());
  }
  for (unsigned int layer = 0; layer < clusters_.size(); layer++) {
    clusters_[layer].followers.resize(clusters_[layer].x.size());
  }

  auto lastLayerPerSide = (unsigned int)(rhtools_.lastLayer(false));
  unsigned int maxLayer = 2 * lastLayerPerSide - 1;
  std::vector<int> numberOfClustersPerLayer(maxLayer, 0);
  for (unsigned int i = 0; i <= maxLayer; i++) {
    calculateLocalDensity(input.tiles, i, layerIdx2layerandSoa);
  }
  for (unsigned int i = 0; i <= maxLayer; i++) {
    calculateDistanceToHigher(input.tiles, i, layerIdx2layerandSoa);
  }

  auto nTracksters = findAndAssignTracksters(input.tiles, layerIdx2layerandSoa);
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Reconstructed " << nTracksters << " tracksters" << std::endl;
    dumpClusters(layerIdx2layerandSoa, eventNumber);
  }

  // Build Trackster
  result.resize(nTracksters);

  for (unsigned int layer = 0; layer < clusters_.size(); ++layer) {
    const auto &thisLayer = clusters_[layer];
    if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
      edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Examining Layer: " << layer;
    }
    for (unsigned int lc = 0; lc < thisLayer.x.size(); ++lc) {
      if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
        edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Trackster " << thisLayer.clusterIndex[lc];
      }
      if (thisLayer.clusterIndex[lc] >= 0) {
        if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
          edm::LogVerbatim("PatternRecogntionbyCLUE3D") << " adding lcIdx: " << thisLayer.layerClusterOriginalIdx[lc];
        }
        result[thisLayer.clusterIndex[lc]].vertices().push_back(thisLayer.layerClusterOriginalIdx[lc]);
        result[thisLayer.clusterIndex[lc]].vertex_multiplicity().push_back(1);
        // loop over followers
        for (auto [follower_lyrIdx, follower_soaIdx] : thisLayer.followers[lc]) {
          std::array<unsigned int, 2> edge = {
              {(unsigned int)thisLayer.layerClusterOriginalIdx[lc],
               (unsigned int)clusters_[follower_lyrIdx].layerClusterOriginalIdx[follower_soaIdx]}};
          result[thisLayer.clusterIndex[lc]].edges().push_back(edge);
        }
      }
    }
  }

  result.erase(
      std::remove_if(std::begin(result),
                     std::end(result),
                     [&](auto const &v) { return static_cast<int>(v.vertices().size()) < minNumLayerCluster_; }),
      result.end());
  result.shrink_to_fit();

  ticl::assignPCAtoTracksters(result,
                              input.layerClusters,
                              input.layerClustersTime,
                              rhtools_.getPositionLayer(rhtools_.lastLayerEE(false), false).z());

  // run energy regression and ID
  energyRegressionAndID(input.layerClusters, result);
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    for (auto const &t : result) {
      edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Barycenter: " << t.barycenter();
      edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "LCs: " << t.vertices().size();
      edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Energy: " << t.raw_energy();
      edm::LogVerbatim("PatternRecogntionbyCLUE3D") << "Regressed: " << t.regressed_energy();
    }
  }

  // Dump Tracksters information
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    dumpTracksters(layerIdx2layerandSoa, eventNumber, result);
  }

  // Reset internal clusters_ structure of array for next event
  reset();
  if (PatternRecognitionAlgoBaseT<TILES>::algo_verbosity_ > PatternRecognitionAlgoBaseT<TILES>::Advanced) {
    edm::LogVerbatim("PatternRecogntionbyCLUE3D") << std::endl;
  }
}

template<typename TILES >

void ticl::PatternRecognitionbyCLUE3D< TILES >::reset ( void  )

inlineprivate

Definition at line 79 of file PatternRecognitionbyCLUE3D.h.

References c, and ticl::PatternRecognitionbyCLUE3D< TILES >::clusters_.

                 {
      for (auto& c : clusters_) {
        c.clear();
        c.shrink_to_fit();
      }
    }

Member Data Documentation

template<typename TILES >

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> ticl::PatternRecognitionbyCLUE3D< TILES >::caloGeomToken_

private

Definition at line 96 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

std::vector<ClustersOnLayer> ticl::PatternRecognitionbyCLUE3D< TILES >::clusters_

private

Definition at line 94 of file PatternRecognitionbyCLUE3D.h.

Referenced by ticl::PatternRecognitionbyCLUE3D< TILES >::reset().

template<typename TILES >

const double ticl::PatternRecognitionbyCLUE3D< TILES >::criticalDensity_

private

Definition at line 97 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const double ticl::PatternRecognitionbyCLUE3D< TILES >::criticalEtaPhiDistance_

private

Definition at line 101 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const double ticl::PatternRecognitionbyCLUE3D< TILES >::densityEtaPhiDistanceSqr_

private

Definition at line 99 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const double ticl::PatternRecognitionbyCLUE3D< TILES >::densityOnSameLayer_

private

Definition at line 100 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const int ticl::PatternRecognitionbyCLUE3D< TILES >::densitySiblingLayers_

private

Definition at line 98 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const std::string ticl::PatternRecognitionbyCLUE3D< TILES >::eidInputName_

private

Definition at line 105 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const float ticl::PatternRecognitionbyCLUE3D< TILES >::eidMinClusterEnergy_

private

Definition at line 108 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const int ticl::PatternRecognitionbyCLUE3D< TILES >::eidNClusters_

private

Definition at line 110 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const int ticl::PatternRecognitionbyCLUE3D< TILES >::eidNFeatures_ = 3

staticprivate

Definition at line 115 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const int ticl::PatternRecognitionbyCLUE3D< TILES >::eidNLayers_

private

Definition at line 109 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const std::string ticl::PatternRecognitionbyCLUE3D< TILES >::eidOutputNameEnergy_

private

Definition at line 106 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const std::string ticl::PatternRecognitionbyCLUE3D< TILES >::eidOutputNameId_

private

Definition at line 107 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

tensorflow::Session* ticl::PatternRecognitionbyCLUE3D< TILES >::eidSession_

private

Definition at line 113 of file PatternRecognitionbyCLUE3D.h.

Referenced by ticl::PatternRecognitionbyCLUE3D< TILES >::PatternRecognitionbyCLUE3D().

template<typename TILES >

const std::vector<int> ticl::PatternRecognitionbyCLUE3D< TILES >::filter_on_categories_

private

Definition at line 104 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const int ticl::PatternRecognitionbyCLUE3D< TILES >::minNumLayerCluster_

private

Definition at line 103 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

const double ticl::PatternRecognitionbyCLUE3D< TILES >::outlierMultiplier_

private

Definition at line 102 of file PatternRecognitionbyCLUE3D.h.

template<typename TILES >

hgcal::RecHitTools ticl::PatternRecognitionbyCLUE3D< TILES >::rhtools_

private

Definition at line 112 of file PatternRecognitionbyCLUE3D.h.