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

Inheritance diagram for IntegrandThetaFunction:

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

	IntegrandThetaFunction ()

	IntegrandThetaFunction (const IntegrandThetaFunction &bluePrint)

IntegrandThetaFunction &	operator= (const IntegrandThetaFunction &bluePrint)

void	SetParameterAlpha (double alpha)

void	SetParameterPhi0 (double phi0)

void	SetParameterTheta0 (double theta0)

	~IntegrandThetaFunction () override

Private 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

Private Attributes
double	alpha_

IntegralOverPhiFunction *	fPhi_

double	phi0_

double	theta0_

Static Private Attributes
static const unsigned int	debugLevel_ = 0

Detailed Description

Description: auxialiary class for fixed area isolation cone computation (this class performs the integration over the polar angle)

Implementation: imported into CMSSW on 05/18/2007

Definition at line 33 of file IntegrandThetaFunction.h.

Constructor & Destructor Documentation

IntegrandThetaFunction::IntegrandThetaFunction ( )

Definition at line 40 of file IntegrandThetaFunction.cc.

References alpha_, fPhi_, phi0_, and theta0_.

Referenced by Clone().

   : ROOT::Math::ParamFunction<ROOT::Math::IParametricGradFunctionOneDim>(3)
 {
   theta0_ = 0.; 
   phi0_ = 0.; 
   alpha_ = 0.; 
 
   fPhi_ = new IntegralOverPhiFunction();
 }

IntegrandThetaFunction::IntegrandThetaFunction ( const IntegrandThetaFunction & bluePrint )

Definition at line 50 of file IntegrandThetaFunction.cc.

References alpha_, fPhi_, phi0_, and theta0_.

 {
   theta0_ = bluePrint.theta0_;
   phi0_ = bluePrint.phi0_;
   alpha_ = bluePrint.alpha_;
 
   fPhi_ = new IntegralOverPhiFunction(*bluePrint.fPhi_);
 }

IntegrandThetaFunction::~IntegrandThetaFunction ( )

override

Definition at line 59 of file IntegrandThetaFunction.cc.

References fPhi_.

 {
   delete fPhi_;
 }

Member Function Documentation

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

inlineoverride

Definition at line 46 of file IntegrandThetaFunction.h.

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

46 { return new IntegrandThetaFunction(*this); }

IntegrandThetaFunction::IntegrandThetaFunction

IntegrandThetaFunction()

Definition: IntegrandThetaFunction.cc:40

double IntegrandThetaFunction::DoDerivative ( double x ) const

private

Definition at line 160 of file IntegrandThetaFunction.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 IntegrandThetaFunction::DoEval ( double x ) const

overrideprivate

Definition at line 114 of file IntegrandThetaFunction.cc.

References Abs(), alpha_, debugLevel_, geometryDiff::epsilon, fPhi_, phi0_, Pi, IntegralOverPhiFunction::SetParameterAlpha(), IntegralOverPhiFunction::SetParameterPhi0(), IntegralOverPhiFunction::SetParameterTheta0(), theta0_, and x.

Referenced by Clone(), and DoEvalPar().

 {
 //--- return zero if theta either close  to zero or close to Pi
 //    (numerical expressions might become "NaN"s)
   const double epsilon = 1.e-3;
   if ( x < epsilon || x > (TMath::Pi() - epsilon) ) return 0.;
 
 //--- calculate trigonometric expressions
 //    (dependend on angle theta;
 //     polar angle of point within cone)
   double sinTheta = TMath::Sin(x);
   double cscTheta = 1./sinTheta;
 
   double detJacobi = -cscTheta; // partial derrivative dEta/dTheta (for constant particle density in tau id. isolation cone)
   //double detJacobi = 1.; // ordinary solid angle (FOR TESTING ONLY)
 
 //--- evaluate integral over angle phi
 //    (azimuth angle of point within cone)
   fPhi_->SetParameterTheta0(theta0_);
   fPhi_->SetParameterPhi0(phi0_);
   fPhi_->SetParameterAlpha(alpha_);
 
   double integralOverPhi = (*fPhi_)(x);
 
   if ( debugLevel_ > 0 ) {
     edm::LogVerbatim("") << "integralOverPhi = " << integralOverPhi << std::endl
              << " theta0 = " << theta0_ << std::endl
              << " phi0 = " << phi0_ << std::endl
              << " alpha = " << alpha_ << std::endl
              << " theta = " << x << std::endl
              << std::endl;
   }
   
 //--- integrand for integration over theta
 //    equals 
 //      |dEta/dTheta| * integral over phi * sin(theta)
 //
 //    (o the factor dEta/dTheta represents the particle density as function of theta, 
 //       assuming that the particle density is flat in eta;
 //     o the factor sin(theta) originates from the solid angle surface element 
 //       expressed in spherical polar coordinates)
 //
   //return TMath::Abs(detJacobi)*integralOverPhi*sinTheta;
   return TMath::Abs(detJacobi)*integralOverPhi;
 }  

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

overrideprivate

Definition at line 105 of file IntegrandThetaFunction.cc.

References alpha_, DoEval(), phi0_, theta0_, and x.

Referenced by Clone().

 {
   theta0_ = param[0];
   phi0_ = param[1];
   alpha_ = param[2];
 
   return DoEval(x);
 }

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

overrideprivate

Definition at line 169 of file IntegrandThetaFunction.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 IntegrandThetaFunction::DoParameterGradient	(	double	x,
		double *	paramGradient
	)		const

private

Definition at line 180 of file IntegrandThetaFunction.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;
 }

IntegrandThetaFunction & IntegrandThetaFunction::operator= ( const IntegrandThetaFunction & bluePrint )

Definition at line 68 of file IntegrandThetaFunction.cc.

References alpha_, fPhi_, phi0_, and theta0_.

 {
   theta0_ = bluePrint.theta0_;
   phi0_ = bluePrint.phi0_;
   alpha_ = bluePrint.alpha_;
   
   (*fPhi_) = (*bluePrint.fPhi_);
 
   return (*this);
 }

void IntegrandThetaFunction::SetParameterAlpha ( double alpha )

Definition at line 93 of file IntegrandThetaFunction.cc.

References alpha, and alpha_.

Referenced by ConeAreaFunction::DoEval().

 {
   alpha_ = alpha;
 }

void IntegrandThetaFunction::SetParameterPhi0 ( double phi0 )

Definition at line 88 of file IntegrandThetaFunction.cc.

References normalizedPhi(), and phi0_.

Referenced by ConeAreaFunction::DoEval().

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

void IntegrandThetaFunction::SetParameters ( double const * param )

overrideprivate

Definition at line 98 of file IntegrandThetaFunction.cc.

References alpha_, phi0_, and theta0_.

Referenced by Clone().

 {
   theta0_ = param[0];
   phi0_ = param[1];
   alpha_ = param[2];
 }

void IntegrandThetaFunction::SetParameterTheta0 ( double theta0 )

Definition at line 83 of file IntegrandThetaFunction.cc.

References theta0_.

Referenced by ConeAreaFunction::DoEval().

 {
   theta0_ = theta0;
 }

Member Data Documentation

double IntegrandThetaFunction::alpha_

mutableprivate

Definition at line 59 of file IntegrandThetaFunction.h.

Referenced by DoEval(), DoEvalPar(), IntegrandThetaFunction(), operator=(), SetParameterAlpha(), and SetParameters().

const unsigned int IntegrandThetaFunction::debugLevel_ = 0

staticprivate

Definition at line 63 of file IntegrandThetaFunction.h.

Referenced by DoEval().

IntegralOverPhiFunction* IntegrandThetaFunction::fPhi_

mutableprivate

Definition at line 61 of file IntegrandThetaFunction.h.

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

double IntegrandThetaFunction::phi0_

mutableprivate

Definition at line 58 of file IntegrandThetaFunction.h.

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

double IntegrandThetaFunction::theta0_

mutableprivate

Definition at line 57 of file IntegrandThetaFunction.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation