IntegrandThetaFunction.h

Go to the documentation of this file.

#ifndef PhysicsTools_IsolationUtils_IntegrandThetaFunction_h
#define PhysicsTools_IsolationUtils_IntegrandThetaFunction_h

// -*- C++ -*-
//
// Package:    IntegrandThetaFunction
// Class:      IntegrandThetaFunction
// 
//
// Original Author:  Christian Veelken, UC Davis
//         Created:  Thu Nov  2 13:47:40 CST 2006
// $Id: IntegrandThetaFunction.h,v 1.1 2007/05/23 20:23:08 veelken Exp $
//
//

// ROOT include files
#include <Math/ParamFunction.h>
#include <Math/Integrator.h>

class IntegralOverPhiFunction;

//
// class declaration
//

class IntegrandThetaFunction : public ROOT::Math::ParamFunction<ROOT::Math::IParametricGradFunctionOneDim>
{
 public:
  IntegrandThetaFunction();
  IntegrandThetaFunction(const IntegrandThetaFunction& bluePrint);
  ~IntegrandThetaFunction();

  IntegrandThetaFunction& operator=(const IntegrandThetaFunction& bluePrint);
  
  void SetParameterTheta0(double theta0);
  void SetParameterPhi0(double phi0);
  void SetParameterAlpha(double alpha);

  virtual ROOT::Math::IGenFunction* Clone () const { return new IntegrandThetaFunction(*this); }

 private:
  void SetParameters(double* param);

  double DoEval(double x) const;
  virtual double DoEvalPar(double, const double*) const;
  double DoDerivative(double x) const;
  virtual double DoParameterDerivative(double, const double*, unsigned int) const;
  void DoParameterGradient(double x, double* paramGradient) const;

  mutable double theta0_; // polar angle of cone axis
  mutable double phi0_; // azimuth angle of cone axis
  mutable double alpha_; // opening angle of cone (measured from cone axis)

  mutable IntegralOverPhiFunction* fPhi_;

  static const unsigned int debugLevel_ = 0;
};

#endif