pftools::CaloEllipse Class Reference

#include <CaloEllipse.h>

Public Member Functions
void	addPoint (double x, double y)
	CaloEllipse ()
void	clearPoints ()
double	getEccentricity () const
Point	getMajorMinorAxes (double sigma=1.0) const
Point	getPosition () const
double	getTheta () const
void	makeCaches ()
void	reset ()
void	resetCaches ()
virtual	~CaloEllipse ()

Public Attributes
double	cachedEccentricity_
double	cachedMajor_
double	cachedMinor_
double	cachedTheta_

Private Attributes
std::vector< Point >	dataPoints_

Detailed Description

Definition at line 23 of file CaloEllipse.h.

Constructor & Destructor Documentation

CaloEllipse::CaloEllipse

(

)

Definition at line 15 of file CaloEllipse.cc.

                         {


}

CaloEllipse::~CaloEllipse ( )

virtual

Definition at line 20 of file CaloEllipse.cc.

                          {

}

Member Function Documentation

void pftools::CaloEllipse::addPoint ( double  x, double  y  )

inline

Definition at line 29 of file CaloEllipse.h.

References dataPoints_, and point.

                                      {
        std::pair<double, double> point(x, y);
        dataPoints_.push_back(point);
    }

void pftools::CaloEllipse::clearPoints ( )

inline

Definition at line 34 of file CaloEllipse.h.

References dataPoints_.

                       {
        dataPoints_.clear();
    }

double CaloEllipse::getEccentricity

(

)

const

Definition at line 87 of file CaloEllipse.cc.

References a, b, getMajorMinorAxes(), AlCaHLTBitMon_ParallelJobs::p, and mathSSE::sqrt().

Referenced by pftools::operator<<().

                                          {
    Point p = getMajorMinorAxes();
    double a = p.first;
    double b = p.second;
    if(a == 0)
        return 0;
    double ecc = sqrt((a * a - b * b) / (a * a));
    return ecc;
}

Point CaloEllipse::getMajorMinorAxes

(

double

sigma = 1.0

)

const

Definition at line 48 of file CaloEllipse.cc.

References a, b, funct::cos(), dataPoints_, getPosition(), getTheta(), timingPdfMaker::mean, AlCaHLTBitMon_ParallelJobs::p, funct::sin(), mathSSE::sqrt(), makeMuonMisalignmentScenario::sum_xx, makeMuonMisalignmentScenario::sum_yy, theta(), and X.

Referenced by getEccentricity(), makeCaches(), and pftools::operator<<().

                                                       {

    if (dataPoints_.size() < 2) {
        return Point(0, 0);
    }

    Point mean = getPosition();

    PointCit it = dataPoints_.begin();
    double sum_xx(0.0);
    double sum_yy(0.0);
    double theta = getTheta();

    for (; it != dataPoints_.end(); ++it) {
        const Point& p = *it;
        double X = cos(theta) * (p.first - mean.first) - sin(theta) * (p.second
                - mean.second);
        double Y = sin(theta) * (p.first - mean.first) + cos(theta) * (p.second
                - mean.second);
        sum_xx += X * X;
        sum_yy += Y * Y;
    }

    double a = sigma * sqrt(sum_xx / dataPoints_.size());
    double b = sigma * sqrt(sum_yy / dataPoints_.size());

    double major, minor;

    if(a > b) {
        major = a;
        minor = b;
    } else {
        major = b;
        minor = a;
    }

    return Point(major, minor);
}

Point CaloEllipse::getPosition

(

)

const

Definition at line 97 of file CaloEllipse.cc.

References dataPoints_, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by getMajorMinorAxes(), getTheta(), and pftools::operator<<().

                                     {

    if (dataPoints_.size() < 1) {
        return Point(0, 0);
    }

    double x_tot(0.0);
    double y_tot(0.0);
    PointCit it = dataPoints_.begin();
    for (; it != dataPoints_.end(); ++it) {
        const Point& p = *it;
        x_tot += p.first;
        y_tot += p.second;

    }

    return Point(x_tot / dataPoints_.size(), y_tot / dataPoints_.size());
}

double CaloEllipse::getTheta

(

)

const

Definition at line 23 of file CaloEllipse.cc.

References dataPoints_, getPosition(), timingPdfMaker::mean, AlCaHLTBitMon_ParallelJobs::p, lumiQueryAPI::q, and theta().

Referenced by getMajorMinorAxes(), makeCaches(), and pftools::operator<<().

                                   {
    if (dataPoints_.size() < 2) {
        return 0.0;
    }

    PointVector meanAdj;

    Point mean = getPosition();

    double sum_dxdx(0.0);
    double sum_dydy(0.0);
    double sum_dxdy(0.0);
    PointCit it = dataPoints_.begin();
    for (; it != dataPoints_.end(); ++it) {
        const Point& p = *it;
        Point q(p.first - mean.first, p.second - mean.second);
        meanAdj.push_back(q);
        sum_dxdx += q.first * q.first;
        sum_dydy += q.second * q.second;
        sum_dxdy += q.first * q.second;
    }
    double theta = 0.5 * atan(2.0 * sum_dxdy / (sum_dydy - sum_dxdx));
    return theta;
}

void CaloEllipse::makeCaches

(

)

Definition at line 123 of file CaloEllipse.cc.

References cachedMajor_, cachedMinor_, cachedTheta_, getMajorMinorAxes(), and getTheta().

Referenced by pftools::CandidateWrapper::recompute().

                             {
    cachedTheta_ = getTheta();
    Point axes = getMajorMinorAxes();
    cachedMajor_ = axes.first;
    cachedMinor_ = axes.second;
}

void CaloEllipse::reset

(

void

)

Definition at line 130 of file CaloEllipse.cc.

References dataPoints_, and resetCaches().

Referenced by pftools::CandidateWrapper::CandidateWrapper(), and pftools::CandidateWrapper::reset().

                        {
    dataPoints_.clear();
    resetCaches();
}

void CaloEllipse::resetCaches

(

)

Definition at line 116 of file CaloEllipse.cc.

References cachedMajor_, cachedMinor_, and cachedTheta_.

Referenced by reset().

                              {
    cachedTheta_ = 0.0;
    cachedMinor_ = 0.0;
    cachedMajor_ = 0.0;

}

Member Data Documentation

double pftools::CaloEllipse::cachedEccentricity_

Definition at line 49 of file CaloEllipse.h.

double pftools::CaloEllipse::cachedMajor_

Definition at line 47 of file CaloEllipse.h.

Referenced by makeCaches(), and resetCaches().

double pftools::CaloEllipse::cachedMinor_

Definition at line 48 of file CaloEllipse.h.

Referenced by makeCaches(), and resetCaches().

double pftools::CaloEllipse::cachedTheta_

Definition at line 46 of file CaloEllipse.h.

Referenced by makeCaches(), and resetCaches().

std::vector<Point> pftools::CaloEllipse::dataPoints_

private

Definition at line 59 of file CaloEllipse.h.

Referenced by addPoint(), clearPoints(), getMajorMinorAxes(), getPosition(), getTheta(), and reset().