PFClusterWidthCalculator Class Reference

#include <PFClusterWidthCalculator.h>

Collaboration diagram for PFClusterWidthCalculator:

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

Private Attributes
double	m_etaPhiLowerLimit {}

Detailed Description

Definition at line 18 of file PFClusterWidthCalculator.h.

Constructor & Destructor Documentation

PFClusterWidthCalculator()

PFClusterWidthCalculator::PFClusterWidthCalculator

(

)

Definition at line 11 of file PFClusterWidthCalculator.cxx.

: m_etaPhiLowerLimit(0.0025) {}

~PFClusterWidthCalculator()

PFClusterWidthCalculator::~PFClusterWidthCalculator ( )

default

Member Function Documentation

getPFClusterCoordinateWidth()

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

Definition at line 14 of file PFClusterWidthCalculator.cxx.

                                                                                                                                                                                                {
  PFClusterWidth pfClusterWidth;
  if (nCells <= 1) {
    pfClusterWidth.etaMean = clusterEta;
    pfClusterWidth.phiMean = clusterPhi;
    pfClusterWidth.etaVariance = m_etaPhiLowerLimit;
    pfClusterWidth.phiVariance = m_etaPhiLowerLimit;
    return pfClusterWidth;
  }
 
  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;
  }
 
  pfClusterWidth.etaMean = etaSum/static_cast<double>(nCells);
  pfClusterWidth.phiMean = phiSum/static_cast<double>(nCells);
 
  double varianceCorrection = (double)nCells / (double)(nCells-1);
  pfClusterWidth.etaVariance = std::max(m_etaPhiLowerLimit, varianceCorrection * (etaSum2/static_cast<double>(nCells) - pfClusterWidth.etaMean*pfClusterWidth.etaMean));
  pfClusterWidth.phiVariance = std::max(m_etaPhiLowerLimit, varianceCorrection * (phiSum2/static_cast<double>(nCells) - pfClusterWidth.phiMean*pfClusterWidth.phiMean));
 
  return pfClusterWidth;
}

Member Data Documentation

m_etaPhiLowerLimit

double PFClusterWidthCalculator::m_etaPhiLowerLimit {}

private

Definition at line 28 of file PFClusterWidthCalculator.h.

{};

---

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

• PFClusterWidthCalculator.h
• PFClusterWidthCalculator.cxx