#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().

46 { return new ConeAreaFunction(*this); }

ConeAreaFunction::ConeAreaFunction

ConeAreaFunction()

Definition: ConeAreaFunction.cc:40

double ConeAreaFunction::DoDerivative ( double x ) const

protected

Definition at line 147 of file ConeAreaFunction.cc.

                                                     {
   //--- 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(), and theta0_.

Referenced by 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_.

                                                                       {
   //--- 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.

                                                                                         {
   //--- 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.

                                                                                 {
   //--- 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 Puppi_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 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];
 }