#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 =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 40 of file EgammaTowerIsolation.h.

Constructor & Destructor Documentation

template<unsigned int NC>

EgammaTowerIsolationNew< NC >::EgammaTowerIsolationNew ( )

inline

Definition at line 57 of file EgammaTowerIsolation.h.

57 : nt(0){}

EgammaTowerIsolationNew::nt

const uint32_t nt

Definition: EgammaTowerIsolation.h:86

template<unsigned int NC>

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

inline

Definition at line 120 of file EgammaTowerIsolation.h.

References prof2calltree::count, alignCSCRings::e, EgammaTowerIsolationNew< NC >::eta, EgammaTowerIsolationNew< NC >::extRadius2_, EgammaTowerIsolationNew< NC >::h2, EgammaTowerIsolationNew< NC >::he, i, cmsHarvester::index, EgammaTowerIsolationNew< NC >::initSoa(), EgammaTowerIsolationNew< NC >::intRadius2_, j, gen::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];
   float e[nt];
   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 63 of file EgammaTowerIsolation.h.

63 { delete[] mem;}

EgammaTowerIsolationNew::mem

uint32_t * mem

Definition: EgammaTowerIsolation.h:93

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 65 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolationNew< 1 >::compute(), and EgammaTowerIsolation::getSum().

                                                                                                                                 {
     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 163 of file EgammaTowerIsolation.h.

References prof2calltree::count, reco::deltaR2(), reco::CaloCluster::eta(), eta(), spr::find(), EgammaTowerIsolationNew< NC >::Sum::h2, EgammaTowerIsolationNew< NC >::Sum::h2BC, EgammaTowerIsolationNew< NC >::Sum::he, EgammaTowerIsolationNew< NC >::Sum::heBC, i, j, prof2calltree::last, maxEta, bookConverter::min, nt, convertSQLiteXML::ok, phi, reco::CaloCluster::phi(), 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
 
   bool ok[iu-il];
   for (std::size_t i=il;i!=iu; ++i)
     ok[i-il] = (std::find(first,last,id[i])==last);
   
 
 
   // 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(ok[i-il]) {
       sum.heBC[j] +=he[i]*tt;
       sum.h2BC[j] +=h2[i]*tt;
     }
       }
     }
   }
 }

template<unsigned int NC>

void EgammaTowerIsolationNew< NC >::initSoa ( )

inline

Definition at line 94 of file EgammaTowerIsolation.h.

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 71 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolation::getSum().

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