IntegralOverPhiFunction Class Reference

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

Inheritance diagram for IntegralOverPhiFunction:

Public Member Functions
ROOT::Math::IGenFunction *	Clone () const override
	IntegralOverPhiFunction ()
void	SetParameterAlpha (double alpha)
void	SetParameterPhi0 (double phi0)
void	SetParameterTheta0 (double theta0)
	~IntegralOverPhiFunction () override

Private Member Functions
double	DoDerivative (double x) const
double	DoEval (double x) const override
double	DoEvalPar (double x, const double *param) 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_
unsigned int	numSolutionMax1_
unsigned int	numSolutionMax2_
unsigned int	numSolutionMax3_
unsigned int	numSolutionMax4_
unsigned int	numSolutionMin1_
unsigned int	numSolutionMin2_
unsigned int	numSolutionMin3_
unsigned int	numSolutionMin4_
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 azimuthal angle)

Implementation: imported into CMSSW on 05/18/2007

Definition at line 30 of file IntegralOverPhiFunction.h.

Constructor & Destructor Documentation

IntegralOverPhiFunction()

IntegralOverPhiFunction::IntegralOverPhiFunction

(

)

Definition at line 49 of file IntegralOverPhiFunction.cc.

References alpha_, numSolutionMax1_, numSolutionMax2_, numSolutionMax3_, numSolutionMax4_, numSolutionMin1_, numSolutionMin2_, numSolutionMin3_, numSolutionMin4_, phi0_, and theta0_.

Referenced by Clone().

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

  // !!! ONLY FOR TESTING
  numSolutionMin1_ = 0;
  numSolutionMax1_ = 0;
  numSolutionMin2_ = 0;
  numSolutionMax2_ = 0;
  numSolutionMin3_ = 0;
  numSolutionMax3_ = 0;
  numSolutionMin4_ = 0;
  numSolutionMax4_ = 0;
  //     FOR TESTING ONLY !!!
}

~IntegralOverPhiFunction()

IntegralOverPhiFunction::~IntegralOverPhiFunction ( )

override

Definition at line 67 of file IntegralOverPhiFunction.cc.

References debugLevel_, numSolutionMax1_, numSolutionMax2_, numSolutionMax3_, numSolutionMax4_, numSolutionMin1_, numSolutionMin2_, numSolutionMin3_, and numSolutionMin4_.

                                                  {
  // !!! ONLY FOR TESTING
  if (debugLevel_ > 0) {
    edm::LogVerbatim("") << "<IntegralOverPhiFunction::~IntegralOverPhiFunction>:" << std::endl
                         << " numSolutionMin1 = " << numSolutionMin1_ << std::endl
                         << " numSolutionMax1 = " << numSolutionMax1_ << std::endl
                         << " numSolutionMin2 = " << numSolutionMin2_ << std::endl
                         << " numSolutionMax2 = " << numSolutionMax2_ << std::endl
                         << " numSolutionMin3 = " << numSolutionMin3_ << std::endl
                         << " numSolutionMax3 = " << numSolutionMax3_ << std::endl
                         << " numSolutionMin4 = " << numSolutionMin4_ << std::endl
                         << " numSolutionMax4 = " << numSolutionMax4_ << std::endl;
  }
  //     FOR TESTING ONLY !!!
}

Member Function Documentation

Clone()

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

inlineoverride

Definition at line 39 of file IntegralOverPhiFunction.h.

References IntegralOverPhiFunction().

{ return new IntegralOverPhiFunction(*this); }

DoDerivative()

double IntegralOverPhiFunction::DoDerivative ( double  x ) const

private

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

overrideprivate

Definition at line 111 of file IntegralOverPhiFunction.cc.

References alpha_, DMR_cfg::cerr, checkSolutions(), debugLevel_, geometryDiff::epsilon, mps_fire::i, dqmiolumiharvest::j, METSkim_cff::Max, METSkim_cff::Min, normalizedPhi(), numSolutionMax1_, numSolutionMax2_, numSolutionMax3_, numSolutionMax4_, numSolutionMin1_, numSolutionMin2_, numSolutionMin3_, numSolutionMin4_, phi, phi0_, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, Pi, alignCSCRings::r, alignCSCRings::s, submitPVValidationJobs::t, theta0_, and x.

                                                     {
  //--- return zero if theta either close to zero or close to Pi
  //    (phi not well-defined in that case;
  //     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 theta0;
  //     polar angle of cone axis);
  //    avoid theta0 exactly zero or exactly Pi
  //    (numerical expressions become "NaN"s)
  double theta0 = theta0_;
  if (theta0 < epsilon)
    theta0 = epsilon;
  if (theta0 > (TMath::Pi() - epsilon))
    theta0 = (TMath::Pi() - epsilon);
  double cosTheta0 = TMath::Cos(theta0);
  double cos2Theta0 = TMath::Cos(2 * theta0);
  double sinTheta0 = TMath::Sin(theta0);
  double cotTheta0 = 1. / TMath::Tan(theta0);
  double cscTheta0 = 1. / TMath::Sin(theta0);
  //    (dependend on angle phi0;
  //     azimuth angle of cone axis)
  double cosPhi0 = TMath::Cos(phi0_);
  double tanPhi0 = TMath::Tan(phi0_);
  //    (dependend on angle theta;
  //     polar angle of point within cone)
  double cosTheta = TMath::Cos(x);
  double cos2Theta = TMath::Cos(2 * x);
  double sinTheta = TMath::Sin(x);
  double cotTheta = 1. / TMath::Tan(x);
  double cscTheta = 1. / TMath::Sin(x);
  //    (dependent on angle alpha;
  //     opening angle of cone, measured from cone axis)
  double cosAlpha = TMath::Cos(alpha_);

  double s = -cosPhi0 * cosPhi0 *
             (2 * cosAlpha * cosAlpha + cos2Theta - 4 * cosAlpha * cosTheta * cosTheta0 + cos2Theta0) * sinTheta *
             sinTheta * sinTheta0 * sinTheta0;

  //--- return either zero or 2 Pi
  //    in case argument of square-root is zero or negative
  //    (negative values solely arise from rounding errors);
  //    check whether to return zero or 2 Pi:
  //     o return zero
  //       if |theta- theta0| > alpha,
  //       (theta outside cone of opening angle alpha, so vanishing integral over phi)
  //     o return 2 Pi
  //       if |theta- theta0| < alpha
  //       (theta within cone of opening angle alpha;
  //        actually theta0 < alpha in forward/backward direction,
  //        so that integral includes all phi values, hence yielding 2 Pi)
  if (s <= 0) {
    if (TMath::Abs(x - theta0) > alpha_)
      return 0;
    if (TMath::Abs(x - theta0) <= alpha_)
      return 2 * TMath::Pi();
    std::cerr << "Error in <IntegralOverPhiFunction::operator()>: failed to compute return value !" << std::endl;
  }

  double r = (1. / TMath::Sqrt(2.)) * (cscTheta * cscTheta * cscTheta0 * cscTheta0 * TMath::Sqrt(s) * tanPhi0);
  double t = cosPhi0 * (-cotTheta * cotTheta0 + cosAlpha * cscTheta * cscTheta0);

  double phi[4];
  phi[0] = -TMath::ACos(t - r);
  phi[1] = TMath::ACos(t - r);
  phi[2] = -TMath::ACos(t + r);
  phi[3] = TMath::ACos(t + r);

  if (debugLevel_ > 0) {
    edm::LogVerbatim("") << "phi0 = " << phi0_ << std::endl
                         << "phi[0] = " << phi[0] << " (phi[0] - phi0 = " << (phi[0] - phi0_) * 180 / TMath::Pi() << ")"
                         << std::endl
                         << "phi[1] = " << phi[1] << " (phi[1] - phi0 = " << (phi[1] - phi0_) * 180 / TMath::Pi() << ")"
                         << std::endl
                         << "phi[2] = " << phi[2] << " (phi[2] - phi0 = " << (phi[2] - phi0_) * 180 / TMath::Pi() << ")"
                         << std::endl
                         << "phi[3] = " << phi[3] << " (phi[3] - phi0 = " << (phi[3] - phi0_) * 180 / TMath::Pi() << ")"
                         << std::endl;
  }

  double phiMin = 0.;
  double phiMax = 0.;
  for (unsigned int i = 0; i < 4; ++i) {
    for (unsigned int j = (i + 1); j < 4; ++j) {
      //--- search for the two solutions for phi
      //    that have an equal distance to phi0
      //    and are on either side
      double dPhi_i = phi[i] - phi0_;
      double dPhi_j = phi0_ - phi[j];
      if (TMath::Abs(normalizedPhi(dPhi_i - dPhi_j)) <
          epsilon) {  // map difference in azimuth angle into interval [-pi,+pi] and require it to be negligible
        //if ( TMath::Abs((phi[i] - phi0_) - (phi0_ - phi[j])) < epsilon ) { // map difference in azimuth angle into interval [-pi,+pi] and require it to be negligible
        phiMin = TMath::Min(phi[i], phi[j]);
        phiMax = TMath::Max(phi[i], phi[j]);

        // !!! ONLY FOR TESTING
        if (phi[i] == phiMin)
          checkSolutions(i, numSolutionMin1_, numSolutionMin2_, numSolutionMin3_, numSolutionMin4_);
        if (phi[i] == phiMax)
          checkSolutions(i, numSolutionMax1_, numSolutionMax2_, numSolutionMax3_, numSolutionMax4_);
        if (phi[j] == phiMin)
          checkSolutions(j, numSolutionMin1_, numSolutionMin2_, numSolutionMin3_, numSolutionMin4_);
        if (phi[j] == phiMax)
          checkSolutions(j, numSolutionMax1_, numSolutionMax2_, numSolutionMax3_, numSolutionMax4_);
        //     FOR TESTING ONLY !!!
      }
    }
  }

  if (phiMin == 0 && phiMax == 0) {
    edm::LogError("") << "failed to compute Return Value !" << std::endl;
  }

  return TMath::Abs(normalizedPhi(phi0_ - phiMin)) + TMath::Abs(normalizedPhi(phiMax - phi0_));
}

DoEvalPar()

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

overrideprivate

Definition at line 101 of file IntegralOverPhiFunction.cc.

References x.

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

  return DoEval(x);
}

DoParameterDerivative()

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

overrideprivate

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

private

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

SetParameterAlpha()

void IntegralOverPhiFunction::SetParameterAlpha

(

double

alpha

)

Definition at line 93 of file IntegralOverPhiFunction.cc.

References isotrackApplyRegressor::alpha, and alpha_.

Referenced by IntegrandThetaFunction::DoEval().

{ alpha_ = alpha; }

SetParameterPhi0()

void IntegralOverPhiFunction::SetParameterPhi0

(

double

phi0

)

Definition at line 89 of file IntegralOverPhiFunction.cc.

References normalizedPhi(), and phi0_.

Referenced by IntegrandThetaFunction::DoEval().

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

SetParameters()

void IntegralOverPhiFunction::SetParameters ( double const *  param )

overrideprivate

Definition at line 95 of file IntegralOverPhiFunction.cc.

References alpha_, phi0_, and theta0_.

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

SetParameterTheta0()

void IntegralOverPhiFunction::SetParameterTheta0

(

double

theta0

)

Definition at line 87 of file IntegralOverPhiFunction.cc.

References theta0_.

Referenced by IntegrandThetaFunction::DoEval().

{ theta0_ = theta0; }

Member Data Documentation

alpha_

double IntegralOverPhiFunction::alpha_

mutableprivate

Definition at line 52 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), SetParameterAlpha(), and SetParameters().

debugLevel_

const unsigned int IntegralOverPhiFunction::debugLevel_ = 0

staticprivate

Definition at line 65 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), and ~IntegralOverPhiFunction().

numSolutionMax1_

unsigned int IntegralOverPhiFunction::numSolutionMax1_

mutableprivate

Definition at line 56 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

numSolutionMax2_

unsigned int IntegralOverPhiFunction::numSolutionMax2_

mutableprivate

Definition at line 58 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

numSolutionMax3_

unsigned int IntegralOverPhiFunction::numSolutionMax3_

mutableprivate

Definition at line 60 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

numSolutionMax4_

unsigned int IntegralOverPhiFunction::numSolutionMax4_

mutableprivate

Definition at line 62 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

numSolutionMin1_

unsigned int IntegralOverPhiFunction::numSolutionMin1_

mutableprivate

Definition at line 55 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

numSolutionMin2_

unsigned int IntegralOverPhiFunction::numSolutionMin2_

mutableprivate

Definition at line 57 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

numSolutionMin3_

unsigned int IntegralOverPhiFunction::numSolutionMin3_

mutableprivate

Definition at line 59 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

numSolutionMin4_

unsigned int IntegralOverPhiFunction::numSolutionMin4_

mutableprivate

Definition at line 61 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), and ~IntegralOverPhiFunction().

phi0_

double IntegralOverPhiFunction::phi0_

mutableprivate

Definition at line 51 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), SetParameterPhi0(), and SetParameters().

theta0_

double IntegralOverPhiFunction::theta0_

mutableprivate

Definition at line 50 of file IntegralOverPhiFunction.h.

Referenced by DoEval(), IntegralOverPhiFunction(), SetParameters(), and SetParameterTheta0().