EgammaTowerIsolationNew< NC > Class Template Reference

#include <EgammaTowerIsolation.h>

Classes
struct	Sum

Public Member Functions
void	compute (bool et, Sum &sum, reco::Candidate const &cand, CaloTowerDetId const *first, CaloTowerDetId const *last) const
void	compute (bool et, Sum &sum, reco::SuperCluster const &sc, CaloTowerDetId const *first, CaloTowerDetId const *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 = 0
const uint32_t	nt
float *	phi
float *	st

Static Public Attributes
static unsigned int	NCuts = NC

Detailed Description

template<unsigned int NC>
class EgammaTowerIsolationNew< NC >

Definition at line 36 of file EgammaTowerIsolation.h.

Constructor & Destructor Documentation

template<unsigned int NC>

EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew ( )

inline

Definition at line 51 of file EgammaTowerIsolation.h.

References electronEcalRecHitIsolationLcone_cfi::extRadius, electronHcalTowerIsolationLcone_cfi::intRadius, and HLT_2018_cff::towers.

: nt(0) {}

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 KineDebug3::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, and EgammaTowerIsolationNew< NC >::st.

    : 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();
  }
}

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

template<unsigned int NC>

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

inline

Definition at line 56 of file EgammaTowerIsolation.h.

References EgHLTOffHistBins_cfi::et, dqmdumpme::first, reco::Candidate::get(), dqmdumpme::last, and SimDataFormats::CaloAnalysis::sc.

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

template<unsigned int NC>

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

inline

Definition at line 156 of file EgammaTowerIsolation.h.

References KineDebug3::count(), reco::deltaR2(), EgammaTowerIsolationNew< NC >::eta, reco::CaloCluster::eta(), EgammaTowerIsolationNew< NC >::extRadius2_, spr::find(), EgammaTowerIsolationNew< NC >::Sum::h2, EgammaTowerIsolationNew< NC >::h2, EgammaTowerIsolationNew< NC >::Sum::h2BC, EgammaTowerIsolationNew< NC >::Sum::he, EgammaTowerIsolationNew< NC >::he, EgammaTowerIsolationNew< NC >::Sum::heBC, mps_fire::i, EgammaTowerIsolationNew< NC >::intRadius2_, dqmiolumiharvest::j, EgammaTowerIsolationNew< NC >::maxEta, min(), EgammaTowerIsolationNew< NC >::NCuts, EgammaTowerIsolationNew< NC >::nt, EgammaTowerIsolationNew< NC >::phi, reco::CaloCluster::phi(), EgammaTowerIsolationNew< NC >::st, and groupFilesInBlocks::tt.

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

template<unsigned int NC>

void EgammaTowerIsolationNew< NC >::initSoa ( )

inline

Definition at line 86 of file EgammaTowerIsolation.h.

References AlCaHLTBitMon_QueryRunRegistry::comp, KineDebug3::count(), beamerCreator::create(), and EgammaTowerIsolationNew< NC >::nt.

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

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

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::eta

Definition at line 79 of file EgammaTowerIsolation.h.

Referenced by Particle.Particle::__str__(), EgammaTowerIsolationNew< NC >::compute(), EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

template<unsigned int NC>

float EgammaTowerIsolationNew< NC >::extRadius2_[NCuts]

Definition at line 73 of file EgammaTowerIsolation.h.

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

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(), EgammaTowerIsolationNew< NC >::compute(), utils.KS::do_test(), utils.Chi2::do_test(), EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), utils.StatisticalTest::get_rank(), and utils.StatisticalTest::set_operands().

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::he

Definition at line 81 of file EgammaTowerIsolation.h.

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

template<unsigned int NC>

uint32_t* EgammaTowerIsolationNew< NC >::id

Definition at line 84 of file EgammaTowerIsolation.h.

template<unsigned int NC>

float EgammaTowerIsolationNew< NC >::intRadius2_[NCuts]

Definition at line 74 of file EgammaTowerIsolation.h.

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

template<unsigned int NC>

float EgammaTowerIsolationNew< NC >::maxEta

Definition at line 76 of file EgammaTowerIsolation.h.

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

template<unsigned int NC>

uint32_t* EgammaTowerIsolationNew< NC >::mem = 0

Definition at line 85 of file EgammaTowerIsolation.h.

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

template<unsigned int NC>

unsigned int EgammaTowerIsolationNew< NC >::NCuts = NC

static

Definition at line 47 of file EgammaTowerIsolation.h.

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

template<unsigned int NC>

const uint32_t EgammaTowerIsolationNew< NC >::nt

Definition at line 78 of file EgammaTowerIsolation.h.

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

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::phi

Definition at line 80 of file EgammaTowerIsolation.h.

Referenced by Particle.Particle::__str__(), EgammaTowerIsolationNew< NC >::compute(), EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew(), and ntupleDataFormat.Track::phiPull().

template<unsigned int NC>

float* EgammaTowerIsolationNew< NC >::st

Definition at line 83 of file EgammaTowerIsolation.h.

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