14namespace Belle2::InvariantMassMuMuCalib {
23 void init(
double Mean,
38 double eval(
double t);
The integrator aims to evaluate convolution of PDFgenLevel and resolution function.
double integralKronrod(double a)
Integration of the PDF which avoids steps and uses variable transformation (Gauss-Konrod rule as back...
double m_tauR
1/slope of the right exponential tail
double m_eps
cut-off term for the power-spectrum caused by the ISR (in GeV)
double m_m0
invariant mass of the collisions
double m_sigmaK
sigma of the gaus in the convolution
double m_C
the coefficient between part below eps and above eps cut-off
double m_sigmaE
sigma of the external gaus
double m_x
the resulting PDF is function of this variable the actual rec-level mass
double m_tauL
1/slope of the left exponential tail
void init(double Mean, double Sigma, double SigmaK, double BMean, double BDelta, double Tau, double SigmaE, double Frac, double M0, double Eps, double CC, double Slope, double X)
Init the parameters of the PDF integrator.
double m_frac
fraction of events in the external gaus
double m_sigma
sigma of the resolution function
double m_bMean
(bRight + bLeft)/2 where bLeft, bRight are the transition points between gaus and exp (in sigma)
double m_slope
power in the power-like spectrum from the ISR
double eval(double t)
evaluate the PDF integrand for given t - the integration variable
double integralTrap(double a, double b)
Simple integration of the PDF for a to b based on the Trapezoidal rule (for validation)
double m_mean
mean position of the resolution function, i.e. (Gaus+Exp tails) conv Gaus
double m_bDelta
(bRight - bLeft)/2 where bLeft, bRight are the transition points between gaus and exp (in sigma)
double integralTrapImp(double Eps, double a)
Integration of the PDF which avoids steps and uses variable transformation (Trapezoidal rule as back-...