HGCalShowerShape Class Reference

#include <HGCalShowerShape.h>

Classes
struct	DeltaPhi

Public Types
typedef math::XYZTLorentzVector	LorentzVector

Public Member Functions
int	coreShowerLength (const l1t::HGCalMulticluster &c3d, const HGCalTriggerGeometryBase &triggerGeometry) const
float	eMax (const l1t::HGCalMulticluster &c3d) const
void	eventSetup (const edm::EventSetup &es)
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
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
	~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
template<typename Delta = std::minus<float>>
float	sigmaXX (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

typedef math::XYZTLorentzVector HGCalShowerShape::LorentzVector

Definition at line 14 of file HGCalShowerShape.h.

Constructor & Destructor Documentation

HGCalShowerShape::HGCalShowerShape ( )

inline

Definition at line 16 of file HGCalShowerShape.h.

: threshold_(0.) {}

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 ( )

inline

Definition at line 19 of file HGCalShowerShape.h.

{}

Member Function Documentation

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, hgcalTopologyTester_cfi::layers, HGCalTriggerTools::layers(), nlayers, pass(), HGCalTriggerGeometryBase::triggerLayer(), and triggerTools_.

Referenced by fillShapes(), and showerLength().

                                                                                              {
  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 (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;
}

float HGCalShowerShape::eMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 391 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                  {
  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;
}

void HGCalShowerShape::eventSetup ( const edm::EventSetup &  es )

inline

Definition at line 21 of file HGCalShowerShape.h.

References HGCalTriggerTools::eventSetup(), firstLayer(), lastLayer(), maxLayer(), and triggerTools_.

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

{ triggerTools_.eventSetup(es); }

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

Definition at line 505 of file HGCalShowerShape.cc.

References coreShowerLength(), l1t::HGCalClusterT< C >::coreShowerLength(), eMax(), l1t::HGCalClusterT< C >::eMax(), firstLayer(), l1t::HGCalClusterT< C >::firstLayer(), l1t::HGCalClusterT< C >::layer10percent(), l1t::HGCalClusterT< C >::layer50percent(), l1t::HGCalClusterT< C >::layer90percent(), maxLayer(), l1t::HGCalClusterT< C >::maxLayer(), meanZ(), percentileLayer(), percentileTriggerCells(), showerLength(), l1t::HGCalClusterT< C >::showerLength(), sigmaEtaEtaMax(), l1t::HGCalClusterT< C >::sigmaEtaEtaMax(), sigmaEtaEtaTot(), l1t::HGCalClusterT< C >::sigmaEtaEtaTot(), sigmaPhiPhiMax(), l1t::HGCalClusterT< C >::sigmaPhiPhiMax(), sigmaPhiPhiTot(), l1t::HGCalClusterT< C >::sigmaPhiPhiTot(), sigmaRRMax(), l1t::HGCalClusterT< C >::sigmaRRMax(), sigmaRRMean(), l1t::HGCalClusterT< C >::sigmaRRMean(), sigmaRRTot(), l1t::HGCalClusterT< C >::sigmaRRTot(), sigmaZZ(), l1t::HGCalClusterT< C >::sigmaZZ(), l1t::HGCalClusterT< C >::triggerCells67percent(), l1t::HGCalClusterT< C >::triggerCells90percent(), and l1t::HGCalClusterT< C >::zBarycenter().

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

                                                                                                                  {
  c3d.showerLength(showerLength(c3d));
  c3d.coreShowerLength(coreShowerLength(c3d, triggerGeometry));
  c3d.firstLayer(firstLayer(c3d));
  c3d.maxLayer(maxLayer(c3d));
  c3d.sigmaEtaEtaTot(sigmaEtaEtaTot(c3d));
  c3d.sigmaEtaEtaMax(sigmaEtaEtaMax(c3d));
  c3d.sigmaPhiPhiTot(sigmaPhiPhiTot(c3d));
  c3d.sigmaPhiPhiMax(sigmaPhiPhiMax(c3d));
  c3d.sigmaZZ(sigmaZZ(c3d));
  c3d.sigmaRRTot(sigmaRRTot(c3d));
  c3d.sigmaRRMax(sigmaRRMax(c3d));
  c3d.sigmaRRMean(sigmaRRMean(c3d));
  c3d.eMax(eMax(c3d));
  c3d.zBarycenter(meanZ(c3d));
  c3d.layer10percent(percentileLayer(c3d, triggerGeometry, 0.10));
  c3d.layer50percent(percentileLayer(c3d, triggerGeometry, 0.50));
  c3d.layer90percent(percentileLayer(c3d, triggerGeometry, 0.90));
  c3d.triggerCells67percent(percentileTriggerCells(c3d, 0.67));
  c3d.triggerCells90percent(percentileTriggerCells(c3d, 0.90));
}

int HGCalShowerShape::firstLayer

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 30 of file HGCalShowerShape.cc.

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

Referenced by eventSetup(), 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;
}

int HGCalShowerShape::lastLayer

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 65 of file HGCalShowerShape.cc.

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

Referenced by eventSetup(), and 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;
}

int HGCalShowerShape::maxLayer

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 46 of file HGCalShowerShape.cc.

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

Referenced by eventSetup(), and 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;
}

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(), showerLength(), sigmaRRMax(), sigmaRRMean(), sigmaRRTot(), and sigmaZZ().

                                                                                    {
  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;
}

float HGCalShowerShape::meanZ

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 413 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                   {
  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);
}

template<typename T , typename Tref >

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

inlineprivate

Definition at line 85 of file HGCalShowerShape.h.

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

Referenced by coreShowerLength(), eMax(), firstLayer(), lastLayer(), maxLayer(), meanZ(), percentileLayer(), percentileTriggerCells(), sigmaEtaEtaMax(), sigmaEtaEtaTot(), sigmaPhiPhiMax(), sigmaPhiPhiTot(), sigmaRRMax(), sigmaRRMean(), sigmaRRTot(), and sigmaZZ().

                                                 {
    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;
  }

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

Definition at line 107 of file HGCalShowerShape.cc.

References l1t::HGCalClusterT< C >::constituents(), ForwardEmpty, hgcalTopologyTester_cfi::layers, HGCalTriggerTools::layers(), nlayers, pass(), DiDispStaMuonMonitor_cfi::pt, HLT_2018_cff::pt1, SusyPostProcessor_cff::quantile, caloTruthCellsProducer_cfi::triggerCells, HGCalTriggerGeometryBase::triggerLayer(), and triggerTools_.

Referenced by fillShapes(), and showerLength().

                                                              {
  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 (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;
  int 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 = layers[percentile];
  return percentile + (pt1 - pt0 > 0. ? (pt_threshold - pt0) / (pt1 - pt0) : 0.);
}

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

Definition at line 140 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                                                    {
  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.);
}

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

inline

Definition at line 26 of file HGCalShowerShape.h.

References coreShowerLength(), eMax(), fillShapes(), firstLayer(), lastLayer(), meanX(), meanZ(), percentileLayer(), percentileTriggerCells(), SusyPostProcessor_cff::quantile, CosmicsPD_Skims::radius, sigmaEtaEtaMax(), sigmaEtaEtaTot(), sigmaPhiPhiMax(), sigmaPhiPhiTot(), sigmaRRMax(), sigmaRRMean(), sigmaRRTot(), and sigmaZZ().

Referenced by fillShapes().

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

float HGCalShowerShape::sigmaEtaEtaMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 237 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                            {
  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;
}

float HGCalShowerShape::sigmaEtaEtaTot

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 172 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                            {
  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 SeeTot = sigmaXX(tc_energy_eta, c3d.eta());

  return SeeTot;
}

float HGCalShowerShape::sigmaEtaEtaTot

(

const l1t::HGCalCluster &

c2d

)

const

Definition at line 452 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;
}

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

inlineprivate

Definition at line 81 of file HGCalShowerShape.h.

Referenced by sigmaPhiPhiMax(), and sigmaPhiPhiTot().

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

float HGCalShowerShape::sigmaPhiPhiMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 270 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                            {
  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;
}

float HGCalShowerShape::sigmaPhiPhiTot

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 192 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                            {
  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 SppTot = sigmaPhiPhi(tc_energy_phi, c3d.phi());

  return SppTot;
}

float HGCalShowerShape::sigmaPhiPhiTot

(

const l1t::HGCalCluster &

c2d

)

const

Definition at line 468 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;
}

float HGCalShowerShape::sigmaRRMax

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 303 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                        {
  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;
}

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

Definition at line 337 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                                       {
  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);
}

float HGCalShowerShape::sigmaRRTot

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 212 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                        {
  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 Szz = sigmaXX(tc_energy_r, r_mean);

  return Szz;
}

float HGCalShowerShape::sigmaRRTot

(

const l1t::HGCalCluster &

c2d

)

const

Definition at line 484 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;
}

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

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

inlineprivate

Definition at line 62 of file HGCalShowerShape.h.

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

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

                                                                                               {
    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;
    float X_RMS = sqrt(X_MSE);
    return X_RMS;
  }

float HGCalShowerShape::sigmaZZ

(

const l1t::HGCalMulticluster &

c3d

)

const

Definition at line 431 of file HGCalShowerShape.cc.

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

Referenced by fillShapes(), and showerLength().

                                                                     {
  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 Szz = sigmaXX(tc_energy_z, z_mean);

  return Szz;
}

Member Data Documentation

double HGCalShowerShape::distance_ = 1.

private

Definition at line 94 of file HGCalShowerShape.h.

Referenced by pass().

double HGCalShowerShape::threshold_ = 0.

private

Definition at line 93 of file HGCalShowerShape.h.

Referenced by pass().

HGCalTriggerTools HGCalShowerShape::triggerTools_

private

Definition at line 92 of file HGCalShowerShape.h.

Referenced by coreShowerLength(), eMax(), eventSetup(), firstLayer(), lastLayer(), maxLayer(), percentileLayer(), sigmaEtaEtaMax(), sigmaPhiPhiMax(), sigmaRRMax(), and sigmaRRMean().