HGCalShowerShape Class Reference

#include <HGCalShowerShape.h>

Classes
struct	DeltaPhi

Public Types
typedef math::XYZTLorentzVector	LorentzVector

Public Member Functions
int	bitmap (const l1t::HGCalMulticluster &c3d, int start=1, int end=14, float threshold=0) const
int	coreShowerLength (const l1t::HGCalMulticluster &c3d, const HGCalTriggerGeometryBase &triggerGeometry) const
float	eMax (const l1t::HGCalMulticluster &c3d) const
void	fillShapes (l1t::HGCalMulticluster &, const HGCalTriggerGeometryBase &) const
int	firstLayer (const l1t::HGCalMulticluster &c3d) const
	HGCalShowerShape ()
	HGCalShowerShape (const edm::ParameterSet &conf)
int	lastLayer (const l1t::HGCalMulticluster &c3d) const
int	maxLayer (const l1t::HGCalMulticluster &c3d) const
float	meanZ (const l1t::HGCalMulticluster &c3d) const
float	percentileLayer (const l1t::HGCalMulticluster &c3d, const HGCalTriggerGeometryBase &triggerGeometry, float quantile=0.5) const
float	percentileTriggerCells (const l1t::HGCalMulticluster &c3d, float quantile=0.5) const
void	setGeometry (const HGCalTriggerGeometryBase *const geom)
int	showerLength (const l1t::HGCalMulticluster &c3d) const
float	sigmaEtaEtaMax (const l1t::HGCalMulticluster &c3d) const
float	sigmaEtaEtaTot (const l1t::HGCalMulticluster &c3d) const
float	sigmaEtaEtaTot (const l1t::HGCalCluster &c2d) const
float	sigmaPhiPhiMax (const l1t::HGCalMulticluster &c3d) const
float	sigmaPhiPhiTot (const l1t::HGCalMulticluster &c3d) const
float	sigmaPhiPhiTot (const l1t::HGCalCluster &c2d) const
float	sigmaRRMax (const l1t::HGCalMulticluster &c3d) const
float	sigmaRRMean (const l1t::HGCalMulticluster &c3d, float radius=5.) const
float	sigmaRRTot (const l1t::HGCalMulticluster &c3d) const
float	sigmaRRTot (const l1t::HGCalCluster &c2d) const
float	sigmaZZ (const l1t::HGCalMulticluster &c3d) const
float	sumLayers (const l1t::HGCalMulticluster &c3d, int start=1, int end=0) const
float	varEtaEta (const l1t::HGCalMulticluster &c3d) const
float	varPhiPhi (const l1t::HGCalMulticluster &c3d) const
float	varRR (const l1t::HGCalMulticluster &c3d) const
float	varZZ (const l1t::HGCalMulticluster &c3d) const
	~HGCalShowerShape ()

Private Member Functions
float	meanX (const std::vector< pair< float, float >> &energy_X_tc) const
template<typename T , typename Tref >
bool	pass (const T &obj, const Tref &ref) const
float	sigmaPhiPhi (const std::vector< pair< float, float >> &energy_phi_tc, const float phi_cluster) const
float	sigmaXX (const std::vector< pair< float, float >> &energy_X_tc, const float X_cluster) const
float	varPhiPhi (const std::vector< pair< float, float >> &energy_phi_tc, const float phi_cluster) const
template<typename Delta = std::minus<float>>
float	varXX (const std::vector< pair< float, float >> &energy_X_tc, const float X_cluster) const

Private Attributes
double	distance_ = 1.
double	threshold_ = 0.
HGCalTriggerTools	triggerTools_

Detailed Description

Definition at line 12 of file HGCalShowerShape.h.

Member Typedef Documentation

LorentzVector

typedef math::XYZTLorentzVector HGCalShowerShape::LorentzVector

Definition at line 14 of file HGCalShowerShape.h.

Constructor & Destructor Documentation

HGCalShowerShape() [1/2]

HGCalShowerShape::HGCalShowerShape ( )

inline

Definition at line 16 of file HGCalShowerShape.h.

: threshold_(0.) {}

HGCalShowerShape() [2/2]

HGCalShowerShape::HGCalShowerShape

(

const edm::ParameterSet &

conf

)

Definition at line 9 of file HGCalShowerShape.cc.

    : threshold_(conf.getParameter<double>("shape_threshold")),
      distance_(conf.getParameter<double>("shape_distance")) {}

~HGCalShowerShape()

HGCalShowerShape::~HGCalShowerShape ( )

inline

Definition at line 19 of file HGCalShowerShape.h.

{}

Member Function Documentation

bitmap()

int HGCalShowerShape::bitmap	(	const l1t::HGCalMulticluster &	c3d,
		int	start = 1,
		int	end = 14,
		float	threshold = 0
	)		const

Definition at line 552 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), mps_fire::end, HGCalTriggerTools::getTriggerGeometry(), mps_fire::i, HGCalTriggerGeometryBase::lastTriggerLayer(), nano_mu_digi_cff::layer, hgcalTBTopologyTester_cfi::layers, pass(), DiMuonV_cfg::threshold, caloTruthCellsProducer_cfi::triggerCells, HGCalTriggerTools::triggerLayer(), and triggerTools_.

Referenced by fillShapes().

                                                                                                       {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();
  unsigned nlayers = triggerTools_.getTriggerGeometry()->lastTriggerLayer();
  std::vector<double> layers(nlayers, 0);
  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      unsigned layer = triggerTools_.triggerLayer(id_tc.second->detId());
      if (layer == 0 || layer > nlayers)
        continue;
      layers[layer - 1] += id_tc.second->pt();  //shift by -1 because layer 0 doesn't exist
    }
  }
  uint32_t bitmap = 0;
  const int bitmap_size = 32;
  if (start == 0) {
    edm::LogWarning("DataNotFound") << "Trying to read layer 0 that doesn't exist, defaulted start to layer 1";
    start = 1;
  }
  if ((end - start) + 1 > bitmap_size) {
    edm::LogWarning("TooMuchData") << " Specified bounds cannot fit into bitmap size, defaulting to 0.";
  } else {
    for (int i = start; i <= end; i++) {
      bitmap += (layers[i - 1] > threshold) << (end - (i));
    }
  }
  return bitmap;
}

coreShowerLength()

int HGCalShowerShape::coreShowerLength	(	const l1t::HGCalMulticluster &	c3d,
		const HGCalTriggerGeometryBase &	triggerGeometry
	)		const

Definition at line 81 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), ForwardEmpty, HFNose, HGCalTriggerTools::isNose(), nano_mu_digi_cff::layer, hgcalTBTopologyTester_cfi::layers, HGCalTriggerTools::layers(), pass(), HGCalTriggerGeometryBase::triggerLayer(), and triggerTools_.

Referenced by fillShapes().

                                                                                              {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();
  unsigned nlayers = triggerTools_.layers(ForwardSubdetector::ForwardEmpty);
  std::vector<bool> layers(nlayers);
  for (const auto& id_cluster : clustersPtrs) {
    if (!pass(*id_cluster.second, c3d))
      continue;
    unsigned layer = triggerGeometry.triggerLayer(id_cluster.second->detId());
    if (triggerTools_.isNose(id_cluster.second->detId()))
      nlayers = triggerTools_.layers(ForwardSubdetector::HFNose);
    else {
      nlayers = triggerTools_.layers(ForwardSubdetector::ForwardEmpty);
    }
    if (layer == 0 || layer > nlayers)
      continue;
    layers[layer - 1] = true;  //layer 0 doesn't exist, so shift by -1
  }
  int length = 0;
  int maxlength = 0;
  for (bool layer : layers) {
    if (layer)
      length++;
    else
      length = 0;
    if (length > maxlength)
      maxlength = length;
  }
  return maxlength;
}

eMax()

float HGCalShowerShape::eMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 407 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), pass(), and triggerTools_.

Referenced by fillShapes().

                                                                  {
  std::unordered_map<int, float> layer_energy;

  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  for (const auto& id_clu : clustersPtrs) {
    if (!pass(*id_clu.second, c3d))
      continue;
    unsigned layer = triggerTools_.layerWithOffset(id_clu.second->detId());
    layer_energy[layer] += id_clu.second->energy();
  }

  float EMax = 0;

  for (const auto& layer : layer_energy) {
    if (layer.second > EMax)
      EMax = layer.second;
  }

  return EMax;
}

fillShapes()

void HGCalShowerShape::fillShapes	(	l1t::HGCalMulticluster &	c3d,
		const HGCalTriggerGeometryBase &	triggerGeometry
	)		const

Definition at line 584 of file HGCalShowerShape.cc.

References bitmap(), coreShowerLength(), eMax(), firstLayer(), HGCEE, HGCalTriggerGeometryBase::lastTriggerLayer(), HGCalTriggerTools::layers(), maxLayer(), meanZ(), percentileLayer(), percentileTriggerCells(), l1t::HGCalClusterT< C >::setCoreShowerLength(), l1t::HGCalClusterT< C >::setEbm0(), l1t::HGCalClusterT< C >::setEbm1(), l1t::HGCalClusterT< C >::setEMax(), l1t::HGCalClusterT< C >::setEmax1layers(), l1t::HGCalClusterT< C >::setEmax3layers(), l1t::HGCalClusterT< C >::setEmax5layers(), l1t::HGCalClusterT< C >::setEot(), l1t::HGCalClusterT< C >::setFirst1layers(), l1t::HGCalClusterT< C >::setFirst3layers(), l1t::HGCalClusterT< C >::setFirst5layers(), l1t::HGCalClusterT< C >::setFirstHcal1layers(), l1t::HGCalClusterT< C >::setFirstHcal3layers(), l1t::HGCalClusterT< C >::setFirstHcal5layers(), l1t::HGCalClusterT< C >::setFirstLayer(), l1t::HGCalClusterT< C >::setHbm(), l1t::HGCalClusterT< C >::setLast1layers(), l1t::HGCalClusterT< C >::setLast3layers(), l1t::HGCalClusterT< C >::setLast5layers(), l1t::HGCalClusterT< C >::setLayer10percent(), l1t::HGCalClusterT< C >::setLayer50percent(), l1t::HGCalClusterT< C >::setLayer90percent(), l1t::HGCalClusterT< C >::setMaxLayer(), l1t::HGCalClusterT< C >::setShowerLength(), l1t::HGCalClusterT< C >::setSigmaEtaEtaMax(), l1t::HGCalClusterT< C >::setSigmaEtaEtaTot(), l1t::HGCalClusterT< C >::setSigmaPhiPhiMax(), l1t::HGCalClusterT< C >::setSigmaPhiPhiTot(), l1t::HGCalClusterT< C >::setSigmaRRMax(), l1t::HGCalClusterT< C >::setSigmaRRMean(), l1t::HGCalClusterT< C >::setSigmaRRTot(), l1t::HGCalClusterT< C >::setSigmaZZ(), l1t::HGCalClusterT< C >::setTriggerCells67percent(), l1t::HGCalClusterT< C >::setTriggerCells90percent(), l1t::HGCalClusterT< C >::setVarEtaEta(), l1t::HGCalClusterT< C >::setVarPhiPhi(), l1t::HGCalClusterT< C >::setVarRR(), l1t::HGCalClusterT< C >::setVarZZ(), l1t::HGCalClusterT< C >::setZBarycenter(), showerLength(), sigmaEtaEtaMax(), sigmaEtaEtaTot(), sigmaPhiPhiMax(), sigmaPhiPhiTot(), sigmaRRMax(), sigmaRRMean(), sigmaRRTot(), sigmaZZ(), sumLayers(), triggerTools_, varEtaEta(), varPhiPhi(), varRR(), and varZZ().

Referenced by HGCalMulticlusteringImpl::finalizeClusters(), HGCalHistoClusteringImpl::finalizeClusters(), and CaloTruthCellsProducer::produce().

                                                                                                                  {
  unsigned hcal_offset = triggerTools_.layers(ForwardSubdetector::HGCEE) / 2;
  unsigned lastlayer = triggerGeometry.lastTriggerLayer();
  const unsigned showermaxlayer = 7;
  c3d.setShowerLength(showerLength(c3d));
  c3d.setCoreShowerLength(coreShowerLength(c3d, triggerGeometry));
  c3d.setFirstLayer(firstLayer(c3d));
  c3d.setMaxLayer(maxLayer(c3d));
  c3d.setSigmaEtaEtaTot(sigmaEtaEtaTot(c3d));
  c3d.setSigmaEtaEtaMax(sigmaEtaEtaMax(c3d));
  c3d.setVarEtaEta(varEtaEta(c3d));
  c3d.setSigmaPhiPhiTot(sigmaPhiPhiTot(c3d));
  c3d.setSigmaPhiPhiMax(sigmaPhiPhiMax(c3d));
  c3d.setVarPhiPhi(varPhiPhi(c3d));
  c3d.setSigmaZZ(sigmaZZ(c3d));
  c3d.setVarZZ(varZZ(c3d));
  c3d.setSigmaRRTot(sigmaRRTot(c3d));
  c3d.setVarRR(varRR(c3d));
  c3d.setSigmaRRMax(sigmaRRMax(c3d));
  c3d.setSigmaRRMean(sigmaRRMean(c3d));
  c3d.setEMax(eMax(c3d));
  c3d.setZBarycenter(meanZ(c3d));
  c3d.setLayer10percent(percentileLayer(c3d, triggerGeometry, 0.10));
  c3d.setLayer50percent(percentileLayer(c3d, triggerGeometry, 0.50));
  c3d.setLayer90percent(percentileLayer(c3d, triggerGeometry, 0.90));
  c3d.setTriggerCells67percent(percentileTriggerCells(c3d, 0.67));
  c3d.setTriggerCells90percent(percentileTriggerCells(c3d, 0.90));
  c3d.setFirst1layers(sumLayers(c3d, 1, 1));
  c3d.setFirst3layers(sumLayers(c3d, 1, 3));
  c3d.setFirst5layers(sumLayers(c3d, 1, 5));
  c3d.setFirstHcal1layers(sumLayers(c3d, hcal_offset, hcal_offset));
  c3d.setFirstHcal3layers(sumLayers(c3d, hcal_offset, hcal_offset + 2));
  c3d.setFirstHcal5layers(sumLayers(c3d, hcal_offset, hcal_offset + 4));
  c3d.setLast1layers(sumLayers(c3d, lastlayer, lastlayer));
  c3d.setLast3layers(sumLayers(c3d, lastlayer - 2, lastlayer));
  c3d.setLast5layers(sumLayers(c3d, lastlayer - 4, lastlayer));
  c3d.setEmax1layers(sumLayers(c3d, showermaxlayer, showermaxlayer));
  c3d.setEmax3layers(sumLayers(c3d, showermaxlayer - 1, showermaxlayer + 1));
  c3d.setEmax5layers(sumLayers(c3d, showermaxlayer - 2, showermaxlayer + 2));
  c3d.setEot(sumLayers(c3d, 1, hcal_offset));
  c3d.setEbm0(bitmap(c3d, 1, hcal_offset, 0));
  c3d.setEbm1(bitmap(c3d, 1, hcal_offset, 1));
  c3d.setHbm(bitmap(c3d, hcal_offset, lastlayer, 0));
}

firstLayer()

int HGCalShowerShape::firstLayer

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 30 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), pass(), and triggerTools_.

Referenced by fillShapes(), and showerLength().

                                                                      {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  int firstLayer = 999;

  for (const auto& id_clu : clustersPtrs) {
    if (!pass(*id_clu.second, c3d))
      continue;
    int layer = triggerTools_.layerWithOffset(id_clu.second->detId());
    if (layer < firstLayer)
      firstLayer = layer;
  }

  return firstLayer;
}

lastLayer()

int HGCalShowerShape::lastLayer

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 65 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), pass(), and triggerTools_.

Referenced by showerLength().

                                                                     {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  int lastLayer = -999;

  for (const auto& id_clu : clustersPtrs) {
    if (!pass(*id_clu.second, c3d))
      continue;
    int layer = triggerTools_.layerWithOffset(id_clu.second->detId());
    if (layer > lastLayer)
      lastLayer = layer;
  }

  return lastLayer;
}

maxLayer()

int HGCalShowerShape::maxLayer

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 46 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), pass(), and triggerTools_.

Referenced by fillShapes().

                                                                    {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();
  std::unordered_map<int, float> layers_pt;
  float max_pt = 0.;
  int max_layer = 0;
  for (const auto& id_cluster : clustersPtrs) {
    if (!pass(*id_cluster.second, c3d))
      continue;
    unsigned layer = triggerTools_.layerWithOffset(id_cluster.second->detId());
    auto itr_insert = layers_pt.emplace(layer, 0.);
    itr_insert.first->second += id_cluster.second->pt();
    if (itr_insert.first->second > max_pt) {
      max_pt = itr_insert.first->second;
      max_layer = layer;
    }
  }
  return max_layer;
}

meanX()

float HGCalShowerShape::meanX ( const std::vector< pair< float, float >> &  energy_X_tc ) const

private

Definition at line 15 of file HGCalShowerShape.cc.

Referenced by meanZ(), sigmaRRMax(), sigmaRRMean(), sigmaRRTot(), varRR(), and varZZ().

                                                                                    {
  float Etot = 0;
  float X_sum = 0;

  for (const auto& energy_X : energy_X_tc) {
    X_sum += energy_X.first * energy_X.second;
    Etot += energy_X.first;
  }

  float X_mean = 0;
  if (Etot > 0)
    X_mean = X_sum / Etot;
  return X_mean;
}

meanZ()

float HGCalShowerShape::meanZ

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 429 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), meanX(), pass(), and caloTruthCellsProducer_cfi::triggerCells.

Referenced by fillShapes().

                                                                   {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  std::vector<std::pair<float, float>> tc_energy_z;

  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      tc_energy_z.emplace_back(id_tc.second->energy(), id_tc.second->position().z());
    }
  }

  return meanX(tc_energy_z);
}

pass()

template<typename T , typename Tref >

bool HGCalShowerShape::pass ( const T &  obj, const Tref &  ref  ) const

inlineprivate

Definition at line 99 of file HGCalShowerShape.h.

References funct::abs(), distance_, mag(), getGTfromDQMFile::obj, amptDefault_cfi::proj, and threshold_.

Referenced by bitmap(), coreShowerLength(), eMax(), firstLayer(), lastLayer(), maxLayer(), meanZ(), percentileLayer(), percentileTriggerCells(), sigmaEtaEtaMax(), sigmaEtaEtaTot(), sigmaPhiPhiMax(), sigmaPhiPhiTot(), sigmaRRMax(), sigmaRRMean(), sigmaRRTot(), sumLayers(), varEtaEta(), varPhiPhi(), varRR(), and varZZ().

                                                 {
    bool pass_threshold = (obj.mipPt() > threshold_);
    GlobalPoint proj(Basic3DVector<float>(obj.position()) / std::abs(obj.position().z()));
    bool pass_distance = ((proj - ref.centreProj()).mag() < distance_);
    return pass_threshold && pass_distance;
  }

percentileLayer()

float HGCalShowerShape::percentileLayer	(	const l1t::HGCalMulticluster &	c3d,
		const HGCalTriggerGeometryBase &	triggerGeometry,
		float	quantile = 0.5
	)		const

Definition at line 112 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), ForwardEmpty, HFNose, HGCalTriggerTools::isNose(), nano_mu_digi_cff::layer, hgcalTBTopologyTester_cfi::layers, HGCalTriggerTools::layers(), pass(), DiDispStaMuonMonitor_cfi::pt, HLT_2024v14_cff::pt1, SusyPostProcessor_cff::quantile, caloTruthCellsProducer_cfi::triggerCells, HGCalTriggerGeometryBase::triggerLayer(), and triggerTools_.

Referenced by fillShapes().

                                                              {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();
  unsigned nlayers = triggerTools_.layers(ForwardSubdetector::ForwardEmpty);
  std::vector<double> layers(nlayers, 0);
  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      unsigned layer = triggerGeometry.triggerLayer(id_tc.second->detId());
      if (triggerTools_.isNose(id_tc.second->detId()))
        nlayers = triggerTools_.layers(ForwardSubdetector::HFNose);
      else {
        nlayers = triggerTools_.layers(ForwardSubdetector::ForwardEmpty);
      }
      if (layer == 0 || layer > nlayers)
        continue;
      layers[layer - 1] += id_tc.second->pt();  //layer 0 doesn't exist, so shift by -1
    }
  }
  std::partial_sum(layers.begin(), layers.end(), layers.begin());
  double pt_threshold = layers.back() * quantile;
  unsigned percentile = 0;
  for (double pt : layers) {
    if (pt > pt_threshold) {
      break;
    }
    percentile++;
  }
  // Linear interpolation of percentile value
  double pt0 = (percentile > 0 ? layers[percentile - 1] : 0.);
  double pt1 = (percentile < layers.size() ? layers[percentile] : layers.back());
  return percentile + (pt1 - pt0 > 0. ? (pt_threshold - pt0) / (pt1 - pt0) : 0.);
}

percentileTriggerCells()

float HGCalShowerShape::percentileTriggerCells	(	const l1t::HGCalMulticluster &	c3d,
		float	quantile = 0.5
	)		const

Definition at line 150 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), pass(), SusyPostProcessor_cff::quantile, and caloTruthCellsProducer_cfi::triggerCells.

Referenced by fillShapes().

                                                                                                    {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();
  std::set<double> ordered_tcs;
  double pt_sum = 0.;
  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      ordered_tcs.emplace(id_tc.second->pt());
      pt_sum += id_tc.second->pt();
    }
  }
  double pt_threshold = pt_sum * quantile;
  double partial_sum = 0.;
  double partial_sum_prev = 0.;
  int ntc = 0;
  for (auto itr = ordered_tcs.rbegin(); itr != ordered_tcs.rend(); ++itr) {
    partial_sum_prev = partial_sum;
    partial_sum += *itr;
    ntc++;
    if (partial_sum > pt_threshold) {
      break;
    }
  }
  // Linear interpolation of ntc
  return ntc - 1 +
         (partial_sum - partial_sum_prev > 0. ? (pt_threshold - partial_sum_prev) / (partial_sum - partial_sum_prev)
                                              : 0.);
}

setGeometry()

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

inline

Definition at line 21 of file HGCalShowerShape.h.

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

Referenced by CaloTruthCellsProducer::produce(), HGCalHistoClusteringImpl::setGeometry(), and HGCalMulticlusteringImpl::setGeometry().

{ triggerTools_.setGeometry(geom); }

showerLength()

int HGCalShowerShape::showerLength ( const l1t::HGCalMulticluster &  c3d ) const

inline

Definition at line 26 of file HGCalShowerShape.h.

References firstLayer(), and lastLayer().

Referenced by fillShapes().

                                                          {
    return lastLayer(c3d) - firstLayer(c3d) + 1;
  }  //in number of layers

sigmaEtaEtaMax()

float HGCalShowerShape::sigmaEtaEtaMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 253 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), pass(), sigmaXX(), caloTruthCellsProducer_cfi::triggerCells, and triggerTools_.

Referenced by fillShapes().

                                                                            {
  std::unordered_map<int, std::vector<std::pair<float, float>>> tc_layer_energy_eta;
  std::unordered_map<int, LorentzVector> layer_LV;

  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  for (const auto& id_clu : clustersPtrs) {
    unsigned layer = triggerTools_.layerWithOffset(id_clu.second->detId());

    layer_LV[layer] += id_clu.second->p4();

    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      tc_layer_energy_eta[layer].emplace_back(std::make_pair(id_tc.second->energy(), id_tc.second->eta()));
    }
  }

  float SigmaEtaEtaMax = 0;

  for (auto& tc_iter : tc_layer_energy_eta) {
    const std::vector<std::pair<float, float>>& energy_eta_layer = tc_iter.second;
    const LorentzVector& LV_layer = layer_LV[tc_iter.first];
    float SigmaEtaEtaLayer = sigmaXX(energy_eta_layer, LV_layer.eta());  //RMS wrt layer eta, not wrt c3d eta
    if (SigmaEtaEtaLayer > SigmaEtaEtaMax)
      SigmaEtaEtaMax = SigmaEtaEtaLayer;
  }

  return SigmaEtaEtaMax;
}

sigmaEtaEtaTot() [1/2]

float HGCalShowerShape::sigmaEtaEtaTot

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 182 of file HGCalShowerShape.cc.

References mathSSE::sqrt(), and varEtaEta().

Referenced by fillShapes().

{ return sqrt(varEtaEta(c3d)); }

sigmaEtaEtaTot() [2/2]

float HGCalShowerShape::sigmaEtaEtaTot

(

const l1t::HGCalCluster &

c2d

)

const

Definition at line 468 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), reco::LeafCandidate::eta(), pass(), and sigmaXX().

                                                                       {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& cellsPtrs = c2d.constituents();

  std::vector<std::pair<float, float>> tc_energy_eta;

  for (const auto& id_cell : cellsPtrs) {
    if (!pass(*id_cell.second, c2d))
      continue;
    tc_energy_eta.emplace_back(std::make_pair(id_cell.second->energy(), id_cell.second->eta()));
  }

  float See = sigmaXX(tc_energy_eta, c2d.eta());

  return See;
}

sigmaPhiPhi()

float HGCalShowerShape::sigmaPhiPhi ( const std::vector< pair< float, float >> &  energy_phi_tc, const float  phi_cluster  ) const

inlineprivate

Definition at line 87 of file HGCalShowerShape.h.

References mathSSE::sqrt(), and varPhiPhi().

Referenced by sigmaPhiPhiMax(), and sigmaPhiPhiTot().

                                                                                                       {
    return sqrt(varPhiPhi(energy_phi_tc, phi_cluster));
  }

sigmaPhiPhiMax()

float HGCalShowerShape::sigmaPhiPhiMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 286 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), pass(), sigmaPhiPhi(), caloTruthCellsProducer_cfi::triggerCells, and triggerTools_.

Referenced by fillShapes().

                                                                            {
  std::unordered_map<int, std::vector<std::pair<float, float>>> tc_layer_energy_phi;
  std::unordered_map<int, LorentzVector> layer_LV;

  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  for (const auto& id_clu : clustersPtrs) {
    unsigned layer = triggerTools_.layerWithOffset(id_clu.second->detId());

    layer_LV[layer] += id_clu.second->p4();

    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      tc_layer_energy_phi[layer].emplace_back(std::make_pair(id_tc.second->energy(), id_tc.second->phi()));
    }
  }

  float SigmaPhiPhiMax = 0;

  for (auto& tc_iter : tc_layer_energy_phi) {
    const std::vector<std::pair<float, float>>& energy_phi_layer = tc_iter.second;
    const LorentzVector& LV_layer = layer_LV[tc_iter.first];
    float SigmaPhiPhiLayer = sigmaPhiPhi(energy_phi_layer, LV_layer.phi());  //RMS wrt layer phi, not wrt c3d phi
    if (SigmaPhiPhiLayer > SigmaPhiPhiMax)
      SigmaPhiPhiMax = SigmaPhiPhiLayer;
  }

  return SigmaPhiPhiMax;
}

sigmaPhiPhiTot() [1/2]

float HGCalShowerShape::sigmaPhiPhiTot

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 204 of file HGCalShowerShape.cc.

References mathSSE::sqrt(), and varPhiPhi().

Referenced by fillShapes().

{ return sqrt(varPhiPhi(c3d)); }

sigmaPhiPhiTot() [2/2]

float HGCalShowerShape::sigmaPhiPhiTot

(

const l1t::HGCalCluster &

c2d

)

const

Definition at line 484 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), pass(), reco::LeafCandidate::phi(), and sigmaPhiPhi().

                                                                       {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& cellsPtrs = c2d.constituents();

  std::vector<std::pair<float, float>> tc_energy_phi;

  for (const auto& id_cell : cellsPtrs) {
    if (!pass(*id_cell.second, c2d))
      continue;
    tc_energy_phi.emplace_back(std::make_pair(id_cell.second->energy(), id_cell.second->phi()));
  }

  float Spp = sigmaPhiPhi(tc_energy_phi, c2d.phi());

  return Spp;
}

sigmaRRMax()

float HGCalShowerShape::sigmaRRMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 319 of file HGCalShowerShape.cc.

References funct::abs(), l1t::HGCalClusterT< C >::constituents(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), meanX(), pass(), funct::pow(), alignCSCRings::r, sigmaXX(), mathSSE::sqrt(), caloTruthCellsProducer_cfi::triggerCells, and triggerTools_.

Referenced by fillShapes().

                                                                        {
  std::unordered_map<int, std::vector<std::pair<float, float>>> tc_layer_energy_r;

  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  for (const auto& id_clu : clustersPtrs) {
    unsigned layer = triggerTools_.layerWithOffset(id_clu.second->detId());

    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      float r = (id_tc.second->position().z() != 0.
                     ? std::sqrt(pow(id_tc.second->position().x(), 2) + pow(id_tc.second->position().y(), 2)) /
                           std::abs(id_tc.second->position().z())
                     : 0.);
      tc_layer_energy_r[layer].emplace_back(std::make_pair(id_tc.second->energy(), r));
    }
  }

  float SigmaRRMax = 0;

  for (auto& tc_iter : tc_layer_energy_r) {
    const std::vector<std::pair<float, float>>& energy_r_layer = tc_iter.second;
    float r_mean_layer = meanX(energy_r_layer);
    float SigmaRRLayer = sigmaXX(energy_r_layer, r_mean_layer);
    if (SigmaRRLayer > SigmaRRMax)
      SigmaRRMax = SigmaRRLayer;
  }

  return SigmaRRMax;
}

sigmaRRMean()

float HGCalShowerShape::sigmaRRMean	(	const l1t::HGCalMulticluster &	c3d,
		float	radius = 5.
	)		const

Definition at line 353 of file HGCalShowerShape.cc.

References funct::abs(), l1t::HGCalClusterT< C >::constituents(), PVValHelper::dx, PVValHelper::dy, reco::LeafCandidate::energy(), nano_mu_digi_cff::layer, HGCalTriggerTools::layerWithOffset(), meanX(), pass(), l1t::HGCalTriggerCell::position(), alignCSCRings::r, CosmicsPD_Skims::radius, sigmaXX(), mathSSE::sqrt(), caloTruthCellsProducer_cfi::triggerCells, triggerTools_, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by fillShapes().

                                                                                       {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();
  // group trigger cells by layer
  std::unordered_map<int, std::vector<edm::Ptr<l1t::HGCalTriggerCell>>> layers_tcs;
  for (const auto& id_clu : clustersPtrs) {
    unsigned layer = triggerTools_.layerWithOffset(id_clu.second->detId());
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();
    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      layers_tcs[layer].emplace_back(id_tc.second);
    }
  }

  // Select trigger cells within X cm of the max TC in the layer
  std::unordered_map<int, std::vector<std::pair<float, float>>> tc_layers_energy_r;
  for (const auto& layer_tcs : layers_tcs) {
    int layer = layer_tcs.first;
    edm::Ptr<l1t::HGCalTriggerCell> max_tc = layer_tcs.second.front();
    for (const auto& tc : layer_tcs.second) {
      if (tc->energy() > max_tc->energy())
        max_tc = tc;
    }
    for (const auto& tc : layer_tcs.second) {
      double dx = tc->position().x() - max_tc->position().x();
      double dy = tc->position().y() - max_tc->position().y();
      double distance_to_max = std::sqrt(dx * dx + dy * dy);
      if (distance_to_max < radius) {
        float r = (tc->position().z() != 0.
                       ? std::sqrt(tc->position().x() * tc->position().x() + tc->position().y() * tc->position().y()) /
                             std::abs(tc->position().z())
                       : 0.);
        tc_layers_energy_r[layer].emplace_back(std::make_pair(tc->energy(), r));
      }
    }
  }

  // Compute srr layer by layer
  std::vector<std::pair<float, float>> layers_energy_srr2;
  for (const auto& layer_energy_r : tc_layers_energy_r) {
    const auto& energies_r = layer_energy_r.second;
    float r_mean_layer = meanX(energies_r);
    float srr = sigmaXX(energies_r, r_mean_layer);
    double energy_sum = 0.;
    for (const auto& energy_r : energies_r) {
      energy_sum += energy_r.first;
    }
    layers_energy_srr2.emplace_back(std::make_pair(energy_sum, srr * srr));
  }
  // Combine all layer srr
  float srr2_mean = meanX(layers_energy_srr2);
  return std::sqrt(srr2_mean);
}

sigmaRRTot() [1/2]

float HGCalShowerShape::sigmaRRTot

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 226 of file HGCalShowerShape.cc.

References mathSSE::sqrt(), and varRR().

Referenced by fillShapes().

{ return sqrt(varRR(c3d)); }

sigmaRRTot() [2/2]

float HGCalShowerShape::sigmaRRTot

(

const l1t::HGCalCluster &

c2d

)

const

Definition at line 500 of file HGCalShowerShape.cc.

References funct::abs(), l1t::HGCalClusterT< C >::constituents(), meanX(), pass(), funct::pow(), alignCSCRings::r, sigmaXX(), and mathSSE::sqrt().

                                                                   {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& cellsPtrs = c2d.constituents();

  std::vector<std::pair<float, float>> tc_energy_r;

  for (const auto& id_cell : cellsPtrs) {
    if (!pass(*id_cell.second, c2d))
      continue;
    float r = (id_cell.second->position().z() != 0.
                   ? std::sqrt(pow(id_cell.second->position().x(), 2) + pow(id_cell.second->position().y(), 2)) /
                         std::abs(id_cell.second->position().z())
                   : 0.);
    tc_energy_r.emplace_back(std::make_pair(id_cell.second->energy(), r));
  }

  float r_mean = meanX(tc_energy_r);
  float Srr = sigmaXX(tc_energy_r, r_mean);

  return Srr;
}

sigmaXX()

float HGCalShowerShape::sigmaXX ( const std::vector< pair< float, float >> &  energy_X_tc, const float  X_cluster  ) const

inlineprivate

Definition at line 81 of file HGCalShowerShape.h.

References mathSSE::sqrt(), and varXX().

Referenced by sigmaEtaEtaMax(), sigmaEtaEtaTot(), sigmaRRMax(), sigmaRRMean(), and sigmaRRTot().

                                                                                               {
    float X_RMS = sqrt(varXX(energy_X_tc, X_cluster));
    return X_RMS;
  }

sigmaZZ()

float HGCalShowerShape::sigmaZZ

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 447 of file HGCalShowerShape.cc.

References mathSSE::sqrt(), and varZZ().

Referenced by fillShapes().

{ return sqrt(varZZ(c3d)); }

sumLayers()

float HGCalShowerShape::sumLayers	(	const l1t::HGCalMulticluster &	c3d,
		int	start = 1,
		int	end = 0
	)		const

Definition at line 521 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), mps_fire::end, DivergingColor::frac, HGCalTriggerTools::getTriggerGeometry(), mps_fire::i, HGCalTriggerGeometryBase::lastTriggerLayer(), nano_mu_digi_cff::layer, hgcalTBTopologyTester_cfi::layers, pass(), compareTotals::tot, caloTruthCellsProducer_cfi::triggerCells, HGCalTriggerTools::triggerLayer(), and triggerTools_.

Referenced by fillShapes().

                                                                                           {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();
  unsigned nlayers = triggerTools_.getTriggerGeometry()->lastTriggerLayer();
  std::vector<double> layers(nlayers, 0);
  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      unsigned layer = triggerTools_.triggerLayer(id_tc.second->detId());
      if (layer == 0 || layer > nlayers)
        continue;
      layers[layer - 1] += id_tc.second->pt();  //shift by -1 because layer 0 doesn't exist
    }
  }
  double sum_pt = 0;
  for (int i = start - 1; i <= end - 1; i++) {  //shift by -1 because layer 1 is layers[0]
    sum_pt += layers[i];
  }
  double tot = 0;
  for (unsigned i = 0; i < layers.size(); ++i) {
    tot += layers[i];
  }
  float frac = 0;
  if (tot > 0) {
    frac = sum_pt / tot;
  }
  return frac;
}

varEtaEta()

float HGCalShowerShape::varEtaEta

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 184 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), reco::LeafCandidate::eta(), pass(), caloTruthCellsProducer_cfi::triggerCells, and varXX().

Referenced by fillShapes(), and sigmaEtaEtaTot().

                                                                       {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  std::vector<std::pair<float, float>> tc_energy_eta;

  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      tc_energy_eta.emplace_back(std::make_pair(id_tc.second->energy(), id_tc.second->eta()));
    }
  }

  float varee = varXX(tc_energy_eta, c3d.eta());

  return varee;
}

varPhiPhi() [1/2]

float HGCalShowerShape::varPhiPhi

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 206 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), pass(), reco::LeafCandidate::phi(), and caloTruthCellsProducer_cfi::triggerCells.

Referenced by fillShapes(), sigmaPhiPhi(), and sigmaPhiPhiTot().

                                                                       {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  std::vector<std::pair<float, float>> tc_energy_phi;

  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      tc_energy_phi.emplace_back(std::make_pair(id_tc.second->energy(), id_tc.second->phi()));
    }
  }

  float varphiphi = varPhiPhi(tc_energy_phi, c3d.phi());

  return varphiphi;
}

varPhiPhi() [2/2]

float HGCalShowerShape::varPhiPhi ( const std::vector< pair< float, float >> &  energy_phi_tc, const float  phi_cluster  ) const

inlineprivate

Definition at line 94 of file HGCalShowerShape.h.

                                                                                                     {
    return varXX<DeltaPhi<float>>(energy_phi_tc, phi_cluster);
  }

varRR()

float HGCalShowerShape::varRR

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 228 of file HGCalShowerShape.cc.

References funct::abs(), l1t::HGCalClusterT< C >::constituents(), meanX(), pass(), funct::pow(), alignCSCRings::r, mathSSE::sqrt(), caloTruthCellsProducer_cfi::triggerCells, and varXX().

Referenced by fillShapes(), and sigmaRRTot().

                                                                   {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  std::vector<std::pair<float, float>> tc_energy_r;

  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      float r = (id_tc.second->position().z() != 0.
                     ? std::sqrt(pow(id_tc.second->position().x(), 2) + pow(id_tc.second->position().y(), 2)) /
                           std::abs(id_tc.second->position().z())
                     : 0.);
      tc_energy_r.emplace_back(std::make_pair(id_tc.second->energy(), r));
    }
  }

  float r_mean = meanX(tc_energy_r);
  float vrr = varXX(tc_energy_r, r_mean);

  return vrr;
}

varXX()

template<typename Delta = std::minus<float>>

float HGCalShowerShape::varXX ( const std::vector< pair< float, float >> &  energy_X_tc, const float  X_cluster  ) const

inlineprivate

Definition at line 67 of file HGCalShowerShape.h.

References dumpMFGeometry_cfg::delta, and funct::pow().

Referenced by sigmaXX(), varEtaEta(), varRR(), and varZZ().

                                                                                             {
    Delta delta;
    float Etot = 0;
    float deltaX2_sum = 0;
    for (const auto& energy_X : energy_X_tc) {
      deltaX2_sum += energy_X.first * pow(delta(energy_X.second, X_cluster), 2);
      Etot += energy_X.first;
    }
    float X_MSE = 0;
    if (Etot > 0)
      X_MSE = deltaX2_sum / Etot;
    return X_MSE;
  }

varZZ()

float HGCalShowerShape::varZZ

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 448 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), meanX(), pass(), caloTruthCellsProducer_cfi::triggerCells, and varXX().

Referenced by fillShapes(), and sigmaZZ().

                                                                   {
  const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>>& clustersPtrs = c3d.constituents();

  std::vector<std::pair<float, float>> tc_energy_z;

  for (const auto& id_clu : clustersPtrs) {
    const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalTriggerCell>>& triggerCells = id_clu.second->constituents();

    for (const auto& id_tc : triggerCells) {
      if (!pass(*id_tc.second, c3d))
        continue;
      tc_energy_z.emplace_back(std::make_pair(id_tc.second->energy(), id_tc.second->position().z()));
    }
  }

  float z_mean = meanX(tc_energy_z);
  float varzz = varXX(tc_energy_z, z_mean);

  return varzz;
}

Member Data Documentation

distance_

double HGCalShowerShape::distance_ = 1.

private

Definition at line 108 of file HGCalShowerShape.h.

Referenced by pass().

threshold_

double HGCalShowerShape::threshold_ = 0.

private

Definition at line 107 of file HGCalShowerShape.h.

Referenced by pass().

triggerTools_

HGCalTriggerTools HGCalShowerShape::triggerTools_

private

Definition at line 106 of file HGCalShowerShape.h.

Referenced by bitmap(), coreShowerLength(), eMax(), fillShapes(), firstLayer(), lastLayer(), maxLayer(), percentileLayer(), setGeometry(), sigmaEtaEtaMax(), sigmaPhiPhiMax(), sigmaRRMax(), sigmaRRMean(), and sumLayers().