41 inline double a(
int p,
int l,
int m) {
44 return p%2==0 ?
r : -
r ;
47 void checkCoeffs(
int m1,
int l1,
int m2,
int l2) {
48 if (m1 < 0 || m2 < 0) {
49 throw std::invalid_argument(
"RooLegendre: m coefficients need to be >= 0.");
51 if (l1 < m1 || l2 < m2) {
52 throw std::invalid_argument(
"RooLegendre: m coefficients need to be smaller than corresponding l.");
60 _l1(1),_m1(1),_l2(0),_m2(0)
70 , _ctheta(
"ctheta",
"ctheta", this, ctheta)
71 , _l1(
l),_m1(
m),_l2(0),_m2(0)
80 , _ctheta(
"ctheta",
"ctheta", this, ctheta)
81 , _l1(l1),_m1(m1),_l2(l2),_m2(m2)
90 , _ctheta(
"ctheta", this, other._ctheta)
91 , _l1(other._l1), _m1(other._m1)
92 , _l2(other._l2), _m2(other._m2)
102 double ctheta = std::max(-1., std::min((
double)
_ctheta, +1.));
115void compute(
size_t batchSize,
const int l1,
const int m1,
const int l2,
const int m2,
116 double * __restrict
output,
117 double const * __restrict TH)
119 double legendre1=1.0, legendreMinus1=1.0;
129 for (
size_t i=0; i<batchSize; i++) {
131 output[i] = legendreMinus1;
132 }
else if (TH[i] >= 1.0) {
159 return range ==
nullptr || strlen(range) == 0
160 ? std::abs(
x.min() + 1.) < 1.e-8 && std::abs(
x.max() - 1.) < 1.e-8
161 : std::abs(
x.min(range) + 1.) < 1.e-8 && std::abs(
x.max(range) - 1.) < 1.e-8;
192 for (
int p1=0; 2*p1 <=
_l1-
_m1 ;++p1) {
194 for (
int p2=0; 2*p2 <=
_l2-
_m2 ; ++p2) {
208 if (
_m1==0&&
_m2==0)
return 1;
210 if (
_l1<3&&
_l2<3)
return 1;
215 inline double maxSingle(
int i,
int j) {
223 static const double m2[3] = { 3,3 };
size_t size(const MatrixT &matrix)
retrieve the size of a square matrix
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t r
Abstract base class for objects that represent a real value and implements functionality common to al...
bool matchArgs(const RooArgSet &allDeps, RooArgSet &numDeps, const RooArgProxy &a) const
Utility function for use in getAnalyticalIntegral().
RooArgSet is a container object that can hold multiple RooAbsArg objects.
std::span< const double > at(RooAbsArg const *arg, RooAbsArg const *caller=nullptr)
Compute the associated Legendre polynomials using ROOT::Math::assoc_legendre().
double evaluate() const override
Note: P_0^0 = 1, so P_l^m = P_l^m P_0^0.
Int_t getMaxVal(const RooArgSet &vars) const override
Advertise capability to determine maximum value of function for given set of observables.
double analyticalIntegral(Int_t code, const char *rangeName=nullptr) const override
this was verified to match mathematica for l1 in [0,2], m1 in [0,l1], l2 in [l1,4],...
Int_t getAnalyticalIntegral(RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=nullptr) const override
Interface function getAnalyticalIntergral advertises the analytical integrals that are supported.
void computeBatch(double *output, size_t size, RooFit::Detail::DataMap const &) const override
Base function for computing multiple values of a RooAbsReal.
double maxVal(Int_t code) const override
Return maximum value for set of observables identified by code assigned in getMaxVal.
double assoc_legendre(unsigned l, unsigned m, double x)
Computes the associated Legendre polynomials.
RVec< PromoteTypes< T0, T1 > > pow(const T0 &x, const RVec< T1 > &v)
double legendre(unsigned l, unsigned m, double x)
Double_t Factorial(Int_t i)
Computes factorial(n).
Double_t Gamma(Double_t z)
Computation of gamma(z) for all z.