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

Inheritance diagram for ConeAreaFunction:

Public Member Functions
virtual ROOT::Math::IGenFunction *	Clone () const

	ConeAreaFunction ()

	ConeAreaFunction (const ConeAreaFunction &bluePrint)

ConeAreaFunction &	operator= (const ConeAreaFunction &bluePrint)

void	SetAcceptanceLimit (double etaMax)

void	SetParameterPhi0 (double phi0)

void	SetParameterTheta0 (double theta0)

	~ConeAreaFunction ()

Protected Member Functions
double	DoDerivative (double x) const

double	DoEval (double x) const

virtual double	DoEvalPar (double, const double *) const

virtual double	DoParameterDerivative (double, const double *, unsigned int) const

void	DoParameterGradient (double x, double *paramGradient) const

void	SetParameters (double *param)

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 52 of file ConeAreaFunction.cc.

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

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

ConeAreaFunction::~ConeAreaFunction ( )

Definition at line 63 of file ConeAreaFunction.cc.

References fTheta_, and integrator_.

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

Member Function Documentation

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

inlinevirtual

Reimplemented in ConeAreaRootFunction.

Definition at line 47 of file ConeAreaFunction.h.

References ConeAreaFunction().

47 { return new ConeAreaFunction(*this); }

ConeAreaFunction::ConeAreaFunction

ConeAreaFunction()

Definition: ConeAreaFunction.cc:40

double ConeAreaFunction::DoDerivative ( double x ) const

protected

Definition at line 160 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

protected

Definition at line 139 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

protectedvirtual

Definition at line 124 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

protectedvirtual

Definition at line 169 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 180 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 73 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 112 of file ConeAreaFunction.cc.

References 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 95 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 * param )

protected

Definition at line 100 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];
 }