IntegrandThetaFunction Class Reference

#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() [1/2]

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() [2/2]

IntegrandThetaFunction::IntegrandThetaFunction

(

const IntegrandThetaFunction &

bluePrint

)

Definition at line 49 of file IntegrandThetaFunction.cc.

References alpha_, fPhi_, phi0_, and theta0_.

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

  fPhi_ = new IntegralOverPhiFunction(*bluePrint.fPhi_);
}

~IntegrandThetaFunction()

IntegrandThetaFunction::~IntegrandThetaFunction ( )

override

Definition at line 58 of file IntegrandThetaFunction.cc.

References fPhi_.

{ delete fPhi_; }

Member Function Documentation

Clone()

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

inlineoverride

Definition at line 45 of file IntegrandThetaFunction.h.

References IntegrandThetaFunction().

{ return new IntegrandThetaFunction(*this); }

DoDerivative()

double IntegrandThetaFunction::DoDerivative ( double  x ) const

private

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

DoEval()

double IntegrandThetaFunction::DoEval ( double  x ) const

overrideprivate

Definition at line 101 of file IntegrandThetaFunction.cc.

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

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

DoEvalPar()

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

overrideprivate

Definition at line 92 of file IntegrandThetaFunction.cc.

References x.

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

  return DoEval(x);
}

DoParameterDerivative()

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

overrideprivate

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

DoParameterGradient()

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

private

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

operator=()

IntegrandThetaFunction & IntegrandThetaFunction::operator=

(

const IntegrandThetaFunction &

bluePrint

)

Definition at line 64 of file IntegrandThetaFunction.cc.

References alpha_, fPhi_, phi0_, and theta0_.

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

  (*fPhi_) = (*bluePrint.fPhi_);

  return (*this);
}

SetParameterAlpha()

void IntegrandThetaFunction::SetParameterAlpha

(

double

alpha

)

Definition at line 84 of file IntegrandThetaFunction.cc.

References alpha, and alpha_.

Referenced by ConeAreaFunction::DoEval().

{ alpha_ = alpha; }

SetParameterPhi0()

void IntegrandThetaFunction::SetParameterPhi0

(

double

phi0

)

Definition at line 80 of file IntegrandThetaFunction.cc.

References normalizedPhi(), and phi0_.

Referenced by ConeAreaFunction::DoEval().

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

SetParameters()

void IntegrandThetaFunction::SetParameters ( double const *  param )

overrideprivate

Definition at line 86 of file IntegrandThetaFunction.cc.

References alpha_, phi0_, and theta0_.

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

SetParameterTheta0()

void IntegrandThetaFunction::SetParameterTheta0

(

double

theta0

)

Definition at line 78 of file IntegrandThetaFunction.cc.

References theta0_.

Referenced by ConeAreaFunction::DoEval().

{ theta0_ = theta0; }

Member Data Documentation

alpha_

double IntegrandThetaFunction::alpha_

mutableprivate

Definition at line 58 of file IntegrandThetaFunction.h.

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

debugLevel_

const unsigned int IntegrandThetaFunction::debugLevel_ = 0

staticprivate

Definition at line 62 of file IntegrandThetaFunction.h.

Referenced by DoEval().

fPhi_

IntegralOverPhiFunction* IntegrandThetaFunction::fPhi_

mutableprivate

Definition at line 60 of file IntegrandThetaFunction.h.

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

phi0_

double IntegrandThetaFunction::phi0_

mutableprivate

Definition at line 57 of file IntegrandThetaFunction.h.

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

theta0_

double IntegrandThetaFunction::theta0_

mutableprivate

Definition at line 56 of file IntegrandThetaFunction.h.

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