PFClusterWidthCalculator Class Reference

#include <PFClusterWidthCalculator.h>

Collaboration diagram for PFClusterWidthCalculator:

Public Member Functions
	PFClusterWidthCalculator ()
	~PFClusterWidthCalculator ()=default
std::pair< double, double >	getPFClusterCoordinateWidth (const std::vector< double > &eta, const std::vector< double > &phi, const double &clusterEta, const double &clusterPhi, unsigned int nCells)
double	getEtaMean () const
double	getPhiMean () const

Private Attributes
double	m_etaPhiLowerLimit
double	m_etaMean
double	m_phiMean

Detailed Description

Definition at line 6 of file PFClusterWidthCalculator.h.

Constructor & Destructor Documentation

PFClusterWidthCalculator()

PFClusterWidthCalculator::PFClusterWidthCalculator

(

)

Definition at line 7 of file PFClusterWidthCalculator.cxx.

: m_etaPhiLowerLimit(0.0025),m_etaMean(0.0),m_phiMean(0.0) {}

~PFClusterWidthCalculator()

PFClusterWidthCalculator::~PFClusterWidthCalculator ( )

default

Member Function Documentation

getEtaMean()

double PFClusterWidthCalculator::getEtaMean ( ) const

inline

Definition at line 14 of file PFClusterWidthCalculator.h.

{return m_etaMean;}

getPFClusterCoordinateWidth()

std::pair< double, double > PFClusterWidthCalculator::getPFClusterCoordinateWidth	(	const std::vector< double > &	eta,
		const std::vector< double > &	phi,
		const double &	clusterEta,
		const double &	clusterPhi,
		unsigned int	nCells
	)

Definition at line 11 of file PFClusterWidthCalculator.cxx.

                                                                                                                                                                                                               {
 
    if (nCells <= 1) {
        m_etaMean = clusterEta;
        m_phiMean = clusterPhi;
        return std::make_pair(m_etaPhiLowerLimit,m_etaPhiLowerLimit);
    }
 
    double etaSum(0.0);
    double etaSum2(0.0);
    double phiSum(0.0);
    double phiSum2(0.0);
 
    for(unsigned int iCell=0; iCell<nCells; ++iCell){
        etaSum += eta[iCell];
        etaSum2 += eta[iCell]*eta[iCell];
        double thisCellPhi = eflowAzimuth(phi[iCell]).cycle(clusterPhi);
        phiSum += thisCellPhi;
        phiSum2 += thisCellPhi*thisCellPhi;
 
    }
 
    m_etaMean = etaSum/static_cast<double>(nCells);
    m_phiMean = phiSum/static_cast<double>(nCells);
 
    double varianceCorrection = (double)nCells / (double)(nCells-1);
    double etaVariance = std::max(m_etaPhiLowerLimit,varianceCorrection * (etaSum2/static_cast<double>(nCells) - m_etaMean*m_etaMean));
    double phiVariance = std::max(m_etaPhiLowerLimit,varianceCorrection * (phiSum2/static_cast<double>(nCells) - m_phiMean*m_phiMean));
 
    return std::make_pair(etaVariance,phiVariance);
}

getPhiMean()

double PFClusterWidthCalculator::getPhiMean ( ) const

inline

Definition at line 15 of file PFClusterWidthCalculator.h.

{return m_phiMean;}

Member Data Documentation

m_etaMean

double PFClusterWidthCalculator::m_etaMean

private

Definition at line 20 of file PFClusterWidthCalculator.h.

m_etaPhiLowerLimit

double PFClusterWidthCalculator::m_etaPhiLowerLimit

private

Definition at line 19 of file PFClusterWidthCalculator.h.

m_phiMean

double PFClusterWidthCalculator::m_phiMean

private

Definition at line 21 of file PFClusterWidthCalculator.h.

---

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

• PFClusterWidthCalculator.h
• PFClusterWidthCalculator.cxx