eflowRec::EtaPhiLUT Class Reference

2D look up table for eflowRecClusters More...

#include <EtaPhiLUT.h>

Collaboration diagram for eflowRec::EtaPhiLUT:

Public Member Functions
	EtaPhiLUT (unsigned int nbins=50)
	constructor taking the desired binsize
void	fill (eflowRecClusterContainer &clustersin)
std::vector< eflowRecCluster * >	clustersInCone (float eta, float phi, float dr) const
	collect eflowRecClusters in a given cone

Private Attributes
unsigned int	m_nphiBins
float	m_phiBinSize
	number of bins
std::vector< std::vector< eflowRecCluster * > >	m_phiBinnedLookUpTable
	bin size

Detailed Description

2D look up table for eflowRecClusters

Definition at line 26 of file EtaPhiLUT.h.

Constructor & Destructor Documentation

EtaPhiLUT()

eflowRec::EtaPhiLUT::EtaPhiLUT

(

unsigned int

nbins = 50

)

constructor taking the desired binsize

Definition at line 25 of file EtaPhiLUT.cxx.

                                           :
    m_nphiBins(nbins),
    m_phiBinSize(TWOPI/m_nphiBins),
    m_phiBinnedLookUpTable(m_nphiBins)
    {}

Member Function Documentation

clustersInCone()

std::vector< eflowRecCluster * > eflowRec::EtaPhiLUT::clustersInCone	(	float	eta,
		float	phi,
		float	dr ) const

collect eflowRecClusters in a given cone

Definition at line 43 of file EtaPhiLUT.cxx.

                                                                                              {
    std::vector<eflowRecCluster*> result;
 
    // Indices corresponding to phi range
    unsigned int iPhiMin = phiIndex( phiInRange(phi-dr), m_phiBinSize );
    unsigned int iPhiMax = phiIndex( phiInRange(phi+dr), m_phiBinSize );
 
    // Extract index ranges to iterate over
    std::vector< std::pair<int,int> > iPhiRanges;
    if( iPhiMin < iPhiMax ) {
      iPhiRanges.emplace_back(iPhiMin,iPhiMax );
    } else { // special treatment for phi-wrapping
      iPhiRanges.emplace_back(0,iPhiMax );
      iPhiRanges.emplace_back(iPhiMin,m_nphiBins-1 );
    }
 
    float dr2Cut = dr*dr;
 
    // loop over ranges
    for( auto& range : iPhiRanges ){
      unsigned int indexMin = range.first;
      unsigned int indexMax = range.second;
      for( ; indexMin <= indexMax; ++indexMin ){
        const std::vector<eflowRecCluster*>& phiClusters = m_phiBinnedLookUpTable[indexMin];
 
        // Get iterators for clusters in relevant eta range
        auto it_min = std::lower_bound (phiClusters.begin(), phiClusters.end(), eta-dr, [] (const eflowRecCluster* cl, float etaval) {return cl->getCluster()->eta()<etaval;} );
        auto it_max = std::upper_bound (it_min, phiClusters.end(), eta+dr, [] (float etaval, const eflowRecCluster* cl) {return etaval<cl->getCluster()->eta();} );
 
        // Apply deltaR cut
        for( ;it_min!=it_max;++it_min ){
          const xAOD::CaloCluster& xcluster = *(*it_min)->getCluster();
          float deta = eta - xcluster.eta();
          float dphi = phiInRange(phi - xcluster.phi());
          float dr2 = deta*deta + dphi*dphi;
          if( dr2 < dr2Cut ) result.push_back(*it_min);
        }
      }
    }
    return result;
  }

fill()

void eflowRec::EtaPhiLUT::fill

(

eflowRecClusterContainer &

clustersin

)

Definition at line 31 of file EtaPhiLUT.cxx.

  {
    for(eflowRecCluster* cluster : clustersin) {
      int index = phiIndex(cluster->getCluster()->phi(), m_phiBinSize);
      m_phiBinnedLookUpTable[index].push_back(cluster);
    }
 
    for( auto& vec : m_phiBinnedLookUpTable ) {
      std::stable_sort(vec.begin(),vec.end(),[](const eflowRecCluster* c1, const eflowRecCluster* c2) { return c1->getCluster()->eta() < c2->getCluster()->eta(); });
    }
  }

Member Data Documentation

m_nphiBins

unsigned int eflowRec::EtaPhiLUT::m_nphiBins

private

Definition at line 38 of file EtaPhiLUT.h.

m_phiBinnedLookUpTable

std::vector< std::vector<eflowRecCluster*> > eflowRec::EtaPhiLUT::m_phiBinnedLookUpTable

private

bin size

Definition at line 40 of file EtaPhiLUT.h.

m_phiBinSize

float eflowRec::EtaPhiLUT::m_phiBinSize

private

number of bins

Definition at line 39 of file EtaPhiLUT.h.

---

The documentation for this class was generated from the following files:

• EtaPhiLUT.h
• EtaPhiLUT.cxx