EgammaTowerIsolationNew< NC > Class Template Reference

#include <EgammaTowerIsolation.h>

Classes
struct	Sum

Public Member Functions
template<typename I >
void	compute (bool et, Sum &sum, reco::Candidate const &cand, I first, I last) const
template<typename I >
void	compute (bool et, Sum &sum, reco::SuperCluster const &sc, I first, I last) const
	EgammaTowerIsolationNew ()
	EgammaTowerIsolationNew (float extRadius[NC], float intRadius[NC], CaloTowerCollection const &towers)
void	initSoa ()
void	setRadius (float const extRadius[NC], float const intRadius[NC])
	~EgammaTowerIsolationNew ()

Public Attributes
float *	eta
float	extRadius2_ [NCuts]
float *	h2
float *	he
uint32_t *	id
float	intRadius2_ [NCuts]
float	maxEta
uint32_t *	mem = nullptr
const uint32_t	nt
float *	phi
float *	st

Static Public Attributes
static constexpr unsigned int	NCuts = NC

Detailed Description

template<unsigned int NC>
class EgammaTowerIsolationNew< NC >

Definition at line 36 of file EgammaTowerIsolation.h.

Constructor & Destructor Documentation

EgammaTowerIsolationNew() [1/2]

template<unsigned int NC>

EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew ( )

inline

Definition at line 51 of file EgammaTowerIsolation.h.

: nt(0) {}

EgammaTowerIsolationNew() [2/2]

template<unsigned int NC>

EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew ( float  extRadius[NC], float  intRadius[NC], CaloTowerCollection const &  towers  )

inline

Definition at line 113 of file EgammaTowerIsolation.h.

References submitPVResolutionJobs::count, MillePedeFileConverter_cfg::e, EgammaTowerIsolationNew< NC >::eta, electronEcalRecHitIsolationLcone_cfi::extRadius, EgammaTowerIsolationNew< NC >::extRadius2_, EgammaTowerIsolationNew< NC >::h2, EgammaTowerIsolationNew< NC >::he, mps_fire::i, EgammaTowerIsolationNew< NC >::initSoa(), electronHcalTowerIsolationLcone_cfi::intRadius, EgammaTowerIsolationNew< NC >::intRadius2_, dqmiolumiharvest::j, dqmdumpme::k, EgammaTowerIsolationNew< NC >::NCuts, EgammaTowerIsolationNew< NC >::nt, EgammaTowerIsolationNew< NC >::phi, EgammaTowerIsolationNew< NC >::st, and HLT_2023v12_cff::towers.

    : maxEta(*std::max_element(extRadius, extRadius + NC)), nt(towers.size()) {
  if (nt == 0)
    return;
  initSoa();

#ifdef ETISTATDEBUG
  etiStat::Count::count.create++;
#endif

  for (std::size_t i = 0; i != NCuts; ++i) {
    extRadius2_[i] = extRadius[i] * extRadius[i];
    intRadius2_[i] = intRadius[i] * intRadius[i];
  }

  // sort in eta  (kd-tree anoverkill,does not vectorize...)
  uint32_t index[nt];
#ifdef __clang__
  std::vector<float> e(nt);
#else
  float e[nt];
#endif
  for (std::size_t k = 0; k != nt; ++k) {
    e[k] = towers[k].eta();
    index[k] = k;
    std::push_heap(index, index + k + 1, [&e](uint32_t i, uint32_t j) { return e[i] < e[j]; });
  }
  std::sort_heap(index, index + nt, [&e](uint32_t i, uint32_t j) { return e[i] < e[j]; });

  for (std::size_t i = 0; i != nt; ++i) {
    auto j = index[i];
    eta[i] = towers[j].eta();
    phi[i] = towers[j].phi();
    id[i] = towers[j].id();
    st[i] = 1.f / std::cosh(eta[i]);  // std::sin(towers[j].theta());   //;
    he[i] = towers[j].hadEnergy();
    h2[i] = towers[j].hadEnergyHeOuterLayer();
  }
}

~EgammaTowerIsolationNew()

template<unsigned int NC>

EgammaTowerIsolationNew< NC >::~EgammaTowerIsolationNew ( )

inline

Definition at line 54 of file EgammaTowerIsolation.h.

References EgammaTowerIsolationNew< NC >::mem.

{ delete[] mem; }

Member Function Documentation

compute() [1/2]

template<unsigned int NC>

template<typename I >

void EgammaTowerIsolationNew< NC >::compute ( bool  et, Sum &  sum, reco::Candidate const &  cand, I  first, I  last  ) const

inline

Definition at line 57 of file EgammaTowerIsolation.h.

References EgHLTOffHistBins_cfi::et, dqmdumpme::first, and dqmdumpme::last.

                                                                                    {
    reco::SuperCluster const& sc = *cand.get<reco::SuperClusterRef>().get();
    return compute(et, sum, sc, first, last);
  }

compute() [2/2]

template<unsigned int NC>

template<typename I >

void EgammaTowerIsolationNew< NC >::compute ( bool  et, Sum &  sum, reco::SuperCluster const &  sc, I  first, I  last  ) const

inline

Definition at line 157 of file EgammaTowerIsolation.h.

References submitPVResolutionJobs::count, reco::deltaR2(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, reco::CaloCluster::eta(), spr::find(), dqmdumpme::first, EgammaTowerIsolationNew< NC >::Sum::h2, EgammaTowerIsolationNew< NC >::Sum::h2BC, EgammaTowerIsolationNew< NC >::Sum::he, hcalSimParameters_cfi::he, EgammaTowerIsolationNew< NC >::Sum::heBC, mps_fire::i, dqmiolumiharvest::j, dqmdumpme::last, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, razorScouting_cff::maxEta, SiStripPI::min, nt, PVValHelper::phi, reco::CaloCluster::phi(), groupFilesInBlocks::tt, and pfDeepBoostedJetPreprocessParams_cfi::upper_bound.

                                                                          {
  if (nt == 0)
    return;

#ifdef ETISTATDEBUG
  etiStat::Count::count.comp++;
#endif

  float candEta = sc.eta();
  float candPhi = sc.phi();

  auto lb = std::lower_bound(eta, eta + nt, candEta - maxEta);
  auto ub = std::upper_bound(lb, eta + nt, candEta + maxEta);
  uint32_t il = lb - eta;
  uint32_t iu = std::min(nt, uint32_t(ub - eta + 1));

#ifdef ETISTATDEBUG
  etiStat::Count::count.span += (iu - il);
#endif

  // should be restricted in eta....
  for (std::size_t i = il; i != iu; ++i) {
    float dr2 = reco::deltaR2(candEta, candPhi, eta[i], phi[i]);
    float tt = et ? st[i] : 1.f;
    for (std::size_t j = 0; j != NCuts; ++j) {
      if (dr2 < extRadius2_[j]) {
        if (dr2 >= intRadius2_[j]) {
          sum.he[j] += he[i] * tt;
          sum.h2[j] += h2[i] * tt;
        }
        if (std::find(first, last, id[i]) == last) {
          sum.heBC[j] += he[i] * tt;
          sum.h2BC[j] += h2[i] * tt;
        }
      }
    }
  }
}

initSoa()

template<unsigned int NC>

void EgammaTowerIsolationNew< NC >::initSoa ( )

inline

Definition at line 86 of file EgammaTowerIsolation.h.

References EgammaTowerIsolationNew< NC >::eta, EgammaTowerIsolationNew< NC >::h2, EgammaTowerIsolationNew< NC >::he, EgammaTowerIsolationNew< NC >::mem, EgammaTowerIsolationNew< NC >::nt, EgammaTowerIsolationNew< NC >::phi, and EgammaTowerIsolationNew< NC >::st.

Referenced by EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew().

                 {
    mem = new uint32_t[nt * 6];
    eta = (float*)(mem);
    phi = eta + nt;
    he = phi + nt;
    h2 = he + nt;
    st = h2 + nt;
    id = (uint32_t*)(st) + nt;
  }

setRadius()

template<unsigned int NC>

void EgammaTowerIsolationNew< NC >::setRadius ( float const  extRadius[NC], float const  intRadius[NC]  )

inline

Definition at line 64 of file EgammaTowerIsolation.h.

References electronEcalRecHitIsolationLcone_cfi::extRadius, EgammaTowerIsolationNew< NC >::extRadius2_, mps_fire::i, electronHcalTowerIsolationLcone_cfi::intRadius, EgammaTowerIsolationNew< NC >::intRadius2_, EgammaTowerIsolationNew< NC >::maxEta, and EgammaTowerIsolationNew< NC >::NCuts.

                                                                       {
    for (std::size_t i = 0; i != NCuts; ++i) {
      extRadius2_[i] = extRadius[i] * extRadius[i];
      intRadius2_[i] = intRadius[i] * intRadius[i];
    }
    maxEta = *std::max_element(extRadius, extRadius + NC);
  }

Member Data Documentation

eta

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::eta

Definition at line 79 of file EgammaTowerIsolation.h.

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), EgammaTowerIsolationNew< NC >::initSoa(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

extRadius2_

template<unsigned int NC>

float EgammaTowerIsolationNew< NC >::extRadius2_[NCuts]

Definition at line 73 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), and EgammaTowerIsolationNew< NC >::setRadius().

h2

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::h2

Definition at line 82 of file EgammaTowerIsolation.h.

Referenced by utils.Chi2::absval(), utils.Chi2::check_filled_bins(), utils.BinToBin::checkBinningMatches(), utils.BinToBin1percent::checkBinningMatches(), utils.KS::do_test(), utils.Chi2::do_test(), utils.BinToBin::do_test(), utils.BinToBin1percent::do_test(), EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), utils.StatisticalTest::get_rank(), EgammaTowerIsolationNew< NC >::initSoa(), and utils.StatisticalTest::set_operands().

he

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::he

Definition at line 81 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), and EgammaTowerIsolationNew< NC >::initSoa().

id

template<unsigned int NC>

uint32_t* EgammaTowerIsolationNew< NC >::id

Definition at line 84 of file EgammaTowerIsolation.h.

intRadius2_

template<unsigned int NC>

float EgammaTowerIsolationNew< NC >::intRadius2_[NCuts]

Definition at line 74 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), and EgammaTowerIsolationNew< NC >::setRadius().

maxEta

template<unsigned int NC>

float EgammaTowerIsolationNew< NC >::maxEta

Definition at line 76 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::setRadius().

mem

template<unsigned int NC>

uint32_t* EgammaTowerIsolationNew< NC >::mem = nullptr

Definition at line 85 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::initSoa(), and EgammaTowerIsolationNew< NC >::~EgammaTowerIsolationNew().

NCuts

template<unsigned int NC>

constexpr unsigned int EgammaTowerIsolationNew< NC >::NCuts = NC

static

Definition at line 47 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), and EgammaTowerIsolationNew< NC >::setRadius().

nt

template<unsigned int NC>

const uint32_t EgammaTowerIsolationNew< NC >::nt

Definition at line 78 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), and EgammaTowerIsolationNew< NC >::initSoa().

phi

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::phi

Definition at line 80 of file EgammaTowerIsolation.h.

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), EgammaTowerIsolationNew< NC >::initSoa(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

st

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::st

Definition at line 83 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), and EgammaTowerIsolationNew< NC >::initSoa().