IsolationUtils/interface/IntegrandThetaFunction.h

0001 #ifndef PhysicsTools_IsolationUtils_IntegrandThetaFunction_h
0002 #define PhysicsTools_IsolationUtils_IntegrandThetaFunction_h
0003
0004 // -*- C++ -*-
0005 //
0006 // Package:    IntegrandThetaFunction
0007 // Class:      IntegrandThetaFunction
0008 //
0009 /**\class IntegrandThetaFunction IntegrandThetaFunction.cc PhysicsTools/IsolationUtils/src/IntegrandThetaFunction.cc
0010
0011  Description: auxialiary class for fixed area isolation cone computation
0012               (this class performs the integration over the polar angle)
0013
0014  Implementation:
0015      imported into CMSSW on 05/18/2007
0016 */
0017 //
0018 // Original Author:  Christian Veelken, UC Davis
0019 //         Created:  Thu Nov  2 13:47:40 CST 2006
0020 //
0021 //
0022
0023 // ROOT include files
0024 #include <Math/ParamFunction.h>
0025 #include <Math/Integrator.h>
0026
0027 class IntegralOverPhiFunction;
0028
0029 //
0030 // class declaration
0031 //
0032
0033 class IntegrandThetaFunction : public ROOT::Math::ParamFunction<ROOT::Math::IParametricGradFunctionOneDim> {
0034 public:
0035   IntegrandThetaFunction();
0036   IntegrandThetaFunction(const IntegrandThetaFunction& bluePrint);
0037   ~IntegrandThetaFunction() override;
0038
0039   IntegrandThetaFunction& operator=(const IntegrandThetaFunction& bluePrint);
0040
0041   void SetParameterTheta0(double theta0);
0042   void SetParameterPhi0(double phi0);
0043   void SetParameterAlpha(double alpha);
0044
0045   ROOT::Math::IGenFunction* Clone() const override { return new IntegrandThetaFunction(*this); }
0046
0047 private:
0048   void SetParameters(double const* param) override;
0049
0050   double DoEval(double x) const override;
0051   double DoEvalPar(double, const double*) const override;
0052   double DoDerivative(double x) const;
0053   double DoParameterDerivative(double, const double*, unsigned int) const override;
0054   void DoParameterGradient(double x, double* paramGradient) const;
0055
0056   mutable double theta0_;  // polar angle of cone axis
0057   mutable double phi0_;    // azimuth angle of cone axis
0058   mutable double alpha_;   // opening angle of cone (measured from cone axis)
0059
0060   mutable IntegralOverPhiFunction* fPhi_;
0061
0062   static const unsigned int debugLevel_ = 0;
0063 };
0064
0065 #endif