Belle II Software development
InvariantMassMuMuIntegrator.h
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9
10
11#pragma once
12
13
14namespace Belle2::InvariantMassMuMuCalib {
15
16
19
20 public:
21
23 void init(double Mean,
24 double Sigma,
25 double SigmaK,
26 double BMean,
27 double BDelta,
28 double Tau,
29 double SigmaE,
30 double Frac,
31 double M0,
32 double Eps,
33 double CC,
34 double Slope,
35 double X);
36
38 double eval(double t);
39
41 double integralTrap(double a, double b);
42
44 double integralTrapImp(double Eps, double a);
45
47 double integralKronrod(double a);
48
49 private:
50
51 double m_mean = 4;
52 double m_sigma = 30;
53 double m_sigmaK = 30;
54 double m_bMean = 0;
55 double m_bDelta = 2.6;
56 double m_tauL = 60;
57 double m_tauR = 60;
58
59 double m_sigmaE = 30;
60 double m_frac = 0.1;
61
62 double m_m0 = 10500;
63 double m_eps = 0.01;
64 double m_slope = 0.95;
65
66 double m_x = 10400;
67
68 double m_C = 16;
69
70 };
71
72}
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_eps
cut-off term for the power-spectrum caused by the ISR (in GeV)
double m_C
the coefficient between part below eps and above eps cut-off
double m_x
the resulting PDF is function of this variable the actual rec-level mass
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_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-...