#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

template<unsigned int NC>

EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew ( )

inline

Definition at line 51 of file EgammaTowerIsolation.h.

51 : nt(0) {}

EgammaTowerIsolationNew::nt

const uint32_t nt

Definition: EgammaTowerIsolation.h:78

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.

     : 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.

54 { delete[] mem; }

EgammaTowerIsolationNew::mem

uint32_t * mem

Definition: EgammaTowerIsolation.h:85

Member Function Documentation

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 reco::Candidate::get().

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

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(), PVValHelper::eta, reco::CaloCluster::eta(), spr::find(), EgammaTowerIsolationNew< NC >::Sum::h2, EgammaTowerIsolationNew< NC >::Sum::h2BC, EgammaTowerIsolationNew< NC >::Sum::he, hcalSimParameters_cfi::he, EgammaTowerIsolationNew< NC >::Sum::heBC, mps_fire::i, dqmiolumiharvest::j, cuda_std::lower_bound(), maxEta, min(), nt, PVValHelper::phi, reco::CaloCluster::phi(), groupFilesInBlocks::tt, and cuda_std::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;
         }
       }
     }
   }
 }

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

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 EgammaTowerIsolationNew< NC >::extRadius2_, mps_fire::i, 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);
   }