ConeAreaFunction Class Reference

#include <PhysicsTools/IsolationUtils//src/ConeAreaFunction.cc>

Inheritance diagram for ConeAreaFunction:

Public Member Functions
ROOT::Math::IGenFunction *	Clone () const override
	ConeAreaFunction ()
	ConeAreaFunction (const ConeAreaFunction &bluePrint)
ConeAreaFunction &	operator= (const ConeAreaFunction &bluePrint)
void	SetAcceptanceLimit (double etaMax)
void	SetParameterPhi0 (double phi0)
void	SetParameterTheta0 (double theta0)
	~ConeAreaFunction () override

Protected Member Functions
double	DoDerivative (double x) const
double	DoEval (double x) const override
double	DoEvalPar (double, const double *) const override
double	DoParameterDerivative (double, const double *, unsigned int) const override
void	DoParameterGradient (double x, double *paramGradient) const
void	SetParameters (double const *param) override

Protected Attributes
double	etaMax_
double	phi0_
double	theta0_

Private Attributes
IntegrandThetaFunction *	fTheta_
ROOT::Math::Integrator *	integrator_

Static Private Attributes
static const unsigned int	debugLevel_ = 0

Detailed Description

Description: low level class to compute area of signal cone corresponding to three-dimensional opening angle alpha given as function argument

Implementation: imported into CMSSW on 05/18/2007

Definition at line 33 of file ConeAreaFunction.h.

Constructor & Destructor Documentation

ConeAreaFunction::ConeAreaFunction

(

)

Definition at line 40 of file ConeAreaFunction.cc.

References etaMax_, fTheta_, integrator_, phi0_, and theta0_.

Referenced by Clone().

                                   : ROOT::Math::ParamFunction<ROOT::Math::IParametricGradFunctionOneDim>(2) {
  theta0_ = 0.;
  phi0_ = 0.;

  etaMax_ = -1;

  fTheta_ = new IntegrandThetaFunction();
  integrator_ = new ROOT::Math::Integrator(*fTheta_);
}

ConeAreaFunction::ConeAreaFunction

(

const ConeAreaFunction &

bluePrint

)

Definition at line 50 of file ConeAreaFunction.cc.

References etaMax_, fTheta_, integrator_, phi0_, and theta0_.

    : ROOT::Math::ParamFunction<ROOT::Math::IParametricGradFunctionOneDim>(bluePrint) {
  theta0_ = bluePrint.theta0_;
  phi0_ = bluePrint.phi0_;

  etaMax_ = bluePrint.etaMax_;

  fTheta_ = new IntegrandThetaFunction(*bluePrint.fTheta_);
  integrator_ = new ROOT::Math::Integrator(*fTheta_);
}

ConeAreaFunction::~ConeAreaFunction ( )

override

Definition at line 61 of file ConeAreaFunction.cc.

References fTheta_, and integrator_.

                                    {
  delete fTheta_;  // function gets deleted automatically by Integrator ?
  delete integrator_;
}

Member Function Documentation

ROOT::Math::IGenFunction* ConeAreaFunction::Clone ( ) const

inlineoverride

Definition at line 46 of file ConeAreaFunction.h.

References ConeAreaFunction(), DoDerivative(), DoEval(), DoEvalPar(), DoParameterDerivative(), DoParameterGradient(), SetParameters(), and x.

{ return new ConeAreaFunction(*this); }

double ConeAreaFunction::DoDerivative ( double  x ) const

protected

Definition at line 147 of file ConeAreaFunction.cc.

Referenced by Clone().

                                                    {
  //--- virtual function inherited from ROOT::Math::ParamFunction base class;
  //    not implemented, because not neccessary, but needs to be defined to make code compile...
  edm::LogWarning("") << "Function not implemented yet !" << std::endl;

  return 0.;
}

double ConeAreaFunction::DoEval ( double  x ) const

overrideprotected

Definition at line 127 of file ConeAreaFunction.cc.

References etaMax_, fTheta_, integrator_, phi0_, Pi, IntegrandThetaFunction::SetParameterAlpha(), IntegrandThetaFunction::SetParameterPhi0(), IntegrandThetaFunction::SetParameterTheta0(), theta0_, muonRecoAnalyzer_cfi::thetaMax, and muonRecoAnalyzer_cfi::thetaMin.

Referenced by Clone(), ConeAreaRootFunction::DoEval(), and DoEvalPar().

                                              {
  //--- calculate area covered by cone of opening angle alpha
  //    (measured from cone axis);
  //    evaluate integral over angle theta
  //    (polar angle of point within cone)

  fTheta_->SetParameterTheta0(theta0_);
  fTheta_->SetParameterPhi0(phi0_);
  fTheta_->SetParameterAlpha(x);

  integrator_->SetFunction(*fTheta_);  // set updated parameter values in Integrator

  double thetaMin = (etaMax_ > 0) ? 2 * TMath::ATan(TMath::Exp(-etaMax_)) : 0.;
  double thetaMax = TMath::Pi() - thetaMin;

  double integralOverTheta = integrator_->Integral(thetaMin, thetaMax);

  return integralOverTheta;
}

double ConeAreaFunction::DoEvalPar ( double  x, const double *  param  ) const

overrideprotected

Definition at line 114 of file ConeAreaFunction.cc.

References DoEval(), phi0_, and theta0_.

Referenced by Clone().

                                                                      {
  //--- calculate area covered by cone of opening angle alpha
  //    (measured from cone axis);
  //    evaluate integral over angle theta
  //    (polar angle of point within cone)
  // FIXME: the const above is actually not true as it is implemented now.

  theta0_ = param[0];
  phi0_ = param[1];

  return DoEval(x);
}

double ConeAreaFunction::DoParameterDerivative ( double  , const double *  , unsigned  int  ) const

overrideprotected

Definition at line 155 of file ConeAreaFunction.cc.

Referenced by Clone().

                                                                                        {
  //--- virtual function inherited from ROOT::Math::ParamFunction base class;
  //    not implemented, because not neccessary, but needs to be defined to make code compile...
  edm::LogWarning("") << "Function not implemented yet !" << std::endl;

  return 0.;
}

void ConeAreaFunction::DoParameterGradient ( double  x, double *  paramGradient  ) const

protected

Definition at line 163 of file ConeAreaFunction.cc.

Referenced by Clone().

                                                                                {
  //--- virtual function inherited from ROOT::Math::ParamFunction base class;
  //    not implemented, because not neccessary, but needs to be defined to make code compile...
  edm::LogWarning("") << "Function not implemented yet !" << std::endl;
}

ConeAreaFunction & ConeAreaFunction::operator=

(

const ConeAreaFunction &

bluePrint

)

Definition at line 70 of file ConeAreaFunction.cc.

References etaMax_, fTheta_, integrator_, phi0_, and theta0_.

Referenced by ConeAreaRootFunction::operator=().

                                                                               {
  theta0_ = bluePrint.theta0_;
  phi0_ = bluePrint.phi0_;

  etaMax_ = bluePrint.etaMax_;

  (*fTheta_) = (*bluePrint.fTheta_);
  integrator_->SetFunction(*fTheta_);

  return (*this);
}

void ConeAreaFunction::SetAcceptanceLimit

(

double

etaMax

)

Definition at line 103 of file ConeAreaFunction.cc.

References ALCARECOTkAlBeamHalo_cff::etaMax, and etaMax_.

Referenced by FixedAreaIsolationCone::setAcceptanceLimit().

                                                       {
  //--- check that pseudo-rapidity given as function argument is positive
  //    (assume equal acceptance for positive and negative pseudo-rapidities)

  if (etaMax > 0) {
    etaMax_ = etaMax;
  } else {
    edm::LogError("") << "etaMax cannot be negative !" << std::endl;
  }
}

void ConeAreaFunction::SetParameterPhi0

(

double

phi0

)

Definition at line 88 of file ConeAreaFunction.cc.

References normalizedPhi(), and phi0_.

Referenced by FixedAreaIsolationCone::operator()().

                                                   {
  phi0_ = normalizedPhi(phi0);  // map azimuth angle into interval [-pi,+pi]
}

void ConeAreaFunction::SetParameters ( double const *  param )

overrideprotected

Definition at line 92 of file ConeAreaFunction.cc.

References debugLevel_, phi0_, and theta0_.

Referenced by Clone(), and ConeAreaRootFunction::SetParameters().

                                                        {
  if (debugLevel_ > 0) {
    edm::LogVerbatim("") << "<ConeAreaFunction::SetParameters>:" << std::endl
                         << " theta0 = " << param[0] << std::endl
                         << " phi0 = " << param[1] << std::endl;
  }

  theta0_ = param[0];
  phi0_ = param[1];
}

void ConeAreaFunction::SetParameterTheta0

(

double

theta0

)

Definition at line 86 of file ConeAreaFunction.cc.

References theta0_.

Referenced by FixedAreaIsolationCone::operator()().

{ theta0_ = theta0; }

Member Data Documentation

const unsigned int ConeAreaFunction::debugLevel_ = 0

staticprivate

Definition at line 74 of file ConeAreaFunction.h.

Referenced by SetParameters().

double ConeAreaFunction::etaMax_

mutableprotected

Definition at line 61 of file ConeAreaFunction.h.

Referenced by ConeAreaFunction(), DoEval(), operator=(), and SetAcceptanceLimit().

IntegrandThetaFunction* ConeAreaFunction::fTheta_

mutableprivate

Definition at line 71 of file ConeAreaFunction.h.

Referenced by ConeAreaFunction(), DoEval(), operator=(), and ~ConeAreaFunction().

ROOT::Math::Integrator* ConeAreaFunction::integrator_

mutableprivate

Definition at line 72 of file ConeAreaFunction.h.

Referenced by ConeAreaFunction(), DoEval(), operator=(), and ~ConeAreaFunction().

double ConeAreaFunction::phi0_

mutableprotected

Definition at line 58 of file ConeAreaFunction.h.

Referenced by ConeAreaFunction(), DoEval(), DoEvalPar(), operator=(), SetParameterPhi0(), and SetParameters().

double ConeAreaFunction::theta0_

mutableprotected

Definition at line 57 of file ConeAreaFunction.h.

Referenced by ConeAreaFunction(), DoEval(), DoEvalPar(), operator=(), SetParameters(), and SetParameterTheta0().