a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
MPll.h
Go to the documentation of this file.
1/*
2 * Copyright (C) 2014 HEPfit Collaboration
3 *
4 *
5 * For the licensing terms see doc/COPYING.
6 */
7
8#ifndef MPLL_H
9#define MPLL_H
10
11class StandardModel;
12class F_1;
13class F_2;
14#include <gsl/gsl_integration.h>
15#include <TF1.h>
16#include <TGraph.h>
17#include <TFitResultPtr.h>
18#include <gsl/gsl_spline.h>
19#include <memory>
20
21#define MPllSWITCH 8.2
22#define LATTICE true
23#define BSZ false
24#define GSL_INTERP_DIM_DC 10
25#define SPLINE true
26
27#define NFPOLARBASIS_MPLL false
28
173class MPll{
174public:
182 MPll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i);
183
187 virtual ~MPll();
188
189 gslpp::complex funct_g(double q2);
190
191 gslpp::complex DeltaC9_KD(double q2);
192
199 gslpp::complex h_lambda(double q2);
200
206 gslpp::complex V_L(double q2);
207
208
214 gslpp::complex T_L(double q2);
215
216
222 double S_L(double q2);
223
229 gslpp::complex H_V(double q2);
230
231
237 gslpp::complex H_A(double q2);
238
239
245 gslpp::complex H_S(double q2);
246
247
253 gslpp::complex H_P(double q2);
254
261 gslpp::complex H_nunu(double q2, QCD::lepton lep);
262
263
271 double integrateSigma(int i, double q_min, double q_max);
272
279 double getSigma(int i, double q_2);
280
288 double integrateDelta(int i, double q_min, double q_max);
289
296 double integrateSigmaTree(double q_min, double q_max);
297
302 double getwidth(){
303 updateParameters();
304 return width;
305 }
306
311 std::vector<std::string> initializeMPllParameters();
312
313
314private:
315 const StandardModel& mySM;
316 QCD::lepton lep;
317 QCD::meson meson;
318 QCD::meson pseudoscalar;
319 std::vector<std::string> mpllParameters;
320 std::unique_ptr<F_1> myF_1;
321 std::unique_ptr<F_2> myF_2;
322 bool dispersion = false;
323 bool FixedWCbtos = false;
324 bool MPll_Lattice_flag;
325 bool MPll_GRvDV_flag;
326 bool MPll_DM_flag;
327 bool NeutrinoTree_flag;
328 double mJ2;
329
330 double GF;
331 double ale;
332 double Mlep;
333 double MM;
334 double MP;
335 double Mb;
336 double mu_b;
337 double mu_h;
338 double Mc;
339 double mb_pole;
340 double mc_pole;
341 double Ms;
342 double spectator_charge;
343 double width;
344 double alpha_s_mub;
345 double angmomP;
346 int etaP;
348 double MW;
349 gslpp::complex lambda_t;
350 gslpp::complex h_0;
351 gslpp::complex h_1;
352 gslpp::complex h_2;
353 gslpp::complex r_1;
354 gslpp::complex r_2;
355 gslpp::complex Delta_C9;
356 gslpp::complex exp_Phase;
357
358 double Delta_C7_U;
359 double Delta_C9_U;
360
361 /*LATTICE fit parameters*/
362 double b_0_fplus;
363 double b_1_fplus;
364 double b_2_fplus;
365 double m_fit2_fplus_lat;
366 double b_0_fT;
367 double b_1_fT;
368 double b_2_fT;
369 double m_fit2_fT_lat;
370 double b_0_f0;
371 double b_1_f0;
372 double b_2_f0;
373 double m_fit2_f0_lat;
374 /*LCSR fit parameters*/
375 double r_1_fplus;
376 double r_2_fplus;
377 double m_fit2_fplus;
378 double r_1_fT;
379 double r_2_fT;
380 double m_fit2_fT;
381 double r_2_f0;
382 double m_fit2_f0;
385 gslpp::vector<gslpp::complex> ** allcoeff;
386 gslpp::vector<gslpp::complex> ** allcoeffh;
387 gslpp::vector<gslpp::complex> ** allcoeffprime;
389 gslpp::vector<gslpp::complex> ** allcoeff_noSM;
391 gslpp::vector<gslpp::complex> ** allcoeff_nu;
393 gslpp::vector<gslpp::complex> ** allcoeff_noSM_nu;
395 gslpp::complex C_1;
396 gslpp::complex C_1L_bar;
397 gslpp::complex C_1Lh_bar;
398 gslpp::complex C_2;
399 gslpp::complex C_2L_bar;
400 gslpp::complex C_2Lh_bar;
401 gslpp::complex C_3;
402 gslpp::complex C_4;
403 gslpp::complex C_5;
404 gslpp::complex C_6;
405 gslpp::complex C_7;
406 gslpp::complex C_8;
407 gslpp::complex C_8L;
408 gslpp::complex C_8Lh;
409 gslpp::complex C_9;
410 gslpp::complex C_10;
411 gslpp::complex C_S;
412 gslpp::complex C_P;
414 gslpp::complex C_7p;
415 gslpp::complex C_9p;
416 gslpp::complex C_10p;
417 gslpp::complex C_Sp;
418 gslpp::complex C_Pp;
420 gslpp::complex C_L_nunu_e;
421 gslpp::complex C_R_nunu_e;
422 gslpp::complex C_L_nunu_mu;
423 gslpp::complex C_R_nunu_mu;
424 gslpp::complex C_L_nunu_tau;
425 gslpp::complex C_R_nunu_tau;
427 gsl_interp_accel *acc_Re_deltaC7_QCDF;
428 gsl_interp_accel *acc_Im_deltaC7_QCDF;
429 gsl_interp_accel *acc_Re_deltaC9_QCDF;
430 gsl_interp_accel *acc_Im_deltaC9_QCDF;
431
432 gsl_spline *spline_Re_deltaC7_QCDF;
433 gsl_spline *spline_Im_deltaC7_QCDF;
434 gsl_spline *spline_Re_deltaC9_QCDF;
435 gsl_spline *spline_Im_deltaC9_QCDF;
436
437
438 std::vector<double> ReDeltaC9;
439 std::vector<double> ImDeltaC9;
440 std::vector<double> myq2;
441 TFitResultPtr refres;
442 TFitResultPtr imfres;
444 TGraph gr1;
445 TGraph gr2;
447 TF1 reffit;
448 TF1 imffit;
450 double tmpq2;
452 gslpp::complex H_0_pre;
453 gslpp::complex H_0_WC;
454 gslpp::complex H_c_WC;
455 gslpp::complex H_b_WC;
457 gslpp::complex ihalfMPI;
458 double fournineth;
459 double half;
460 double twothird;
461 double Mc2;
462 double Mb2;
463 double logMc;
464 double logMb;
465 double mu_b2;
466 double fourMc2;
467 double fourMb2;
468 double Mlep2;
469 double NN;
470 double MM2;
471 double MM4;
472 double MP2;
473 double MP4;
474 double MM2mMP2;
475 double twoMP2;
476 double twoMM;
477 double twoMM2;
478 double twoMM2_MMpMP;
479 double twoMM_MbpMs;
480 double S_L_pre;
481 double fourMM2;
482 double twoMboMM;
483 double sixteenM_PI2;
484 double ninetysixM_PI3MM3;
485 double MboMW;
486 double MboMM;
487 double MsoMb;
488 double twoMlepMb;
489 double threeGegen0;
490 double threeGegen1otwo;
491 double M_PI2osix;
492 double twoMc2;
493 double sixMMoMb;
494 double CF;
495 double deltaT_0;
496 double deltaT_1par;
498 gslpp::complex ubar;
499 gslpp::complex arg1;
500 gslpp::complex B01;
501 gslpp::complex B00;
502 gslpp::complex xp;
503 gslpp::complex xm;
504 gslpp::complex yp;
505 gslpp::complex ym;
506 gslpp::complex L1xp;
507 gslpp::complex L1xm;
508 gslpp::complex L1yp;
509 gslpp::complex L1ym;
510 gslpp::complex F87_0;
511 gslpp::complex F87_1;
512 gslpp::complex F87_2;
513 gslpp::complex F87_3;
514 gslpp::complex F29_0;
515 gslpp::complex F29_L1;
516 gslpp::complex F29_1;
517 gslpp::complex F29_2;
518 gslpp::complex F29_3;
519 gslpp::complex F29_L1_1;
520 gslpp::complex F29_L1_2;
521 gslpp::complex F29_L1_3;
522 gslpp::complex F27_0;
523 gslpp::complex F27_1;
524 gslpp::complex F27_2;
525 gslpp::complex F27_3;
526 gslpp::complex F27_L1_1;
527 gslpp::complex F27_L1_2;
528 gslpp::complex F27_L1_3;
529 double F89_0;
530 double F89_1;
531 double F89_2;
532 double F89_3;
533 double Ee;
535 //additional variables for B to K nu nu
536 const double M_PI2 = M_PI * M_PI;
537 double GF4;
538 double MM3;
539 double fM2;
540 double fP2;
541
542 double mtau;
543 double mtau2;
544 double Gammatau;
545 double VusVub_abs2;
546
547 unsigned int fplus_lat_updated;
548 gslpp::vector<double> fplus_lat_cache;
550 unsigned int fT_lat_updated;
551 gslpp::vector<double> fT_lat_cache;
553 unsigned int f0_lat_updated;
554 gslpp::vector<double> f0_lat_cache;
556 unsigned int fplus_updated;
557 gslpp::vector<double> fplus_cache;
559 unsigned int fT_updated;
560 gslpp::vector<double> fT_cache;
562 unsigned int f0_updated;
563 double f0_cache;
565 unsigned int k2_updated;
566 gslpp::vector<double> k2_cache;
568 unsigned int beta_updated;
569 double beta_cache;
571 unsigned int lambda_updated;
573 unsigned int F_updated;
575 unsigned int VL_updated;
577 unsigned int TL_updated;
579 unsigned int SL_updated;
580 gslpp::vector<double> SL_cache;
582 unsigned int N_updated;
583 gslpp::vector<double> N_cache;
584 gslpp::complex Nc_cache;
586 unsigned int C_1_updated;
587 gslpp::complex C_1_cache;
589 unsigned int C_2_updated;
590 gslpp::complex C_2_cache;
592 unsigned int C_3_updated;
593 gslpp::complex C_3_cache;
595 unsigned int C_4_updated;
596 gslpp::complex C_4_cache;
598 unsigned int C_5_updated;
599 gslpp::complex C_5_cache;
601 unsigned int C_6_updated;
602 gslpp::complex C_6_cache;
604 unsigned int C_7_updated;
605 gslpp::complex C_7_cache;
607 unsigned int C_9_updated;
608 gslpp::complex C_9_cache;
610 unsigned int C_10_updated;
611 gslpp::complex C_10_cache;
613 unsigned int C_7p_updated;
614 gslpp::complex C_7p_cache;
616 unsigned int C_9p_updated;
617 gslpp::complex C_9p_cache;
619 unsigned int C_10p_updated;
620 gslpp::complex C_10p_cache;
622 unsigned int C_S_updated;
623 gslpp::complex C_S_cache;
625 unsigned int C_P_updated;
626 gslpp::complex C_P_cache;
628 unsigned int C_Sp_updated;
629 gslpp::complex C_Sp_cache;
631 unsigned int C_Pp_updated;
632 gslpp::complex C_Pp_cache;
634 unsigned int C_2Lh_updated;
635 gslpp::complex C_2Lh_cache;
637 unsigned int C_8Lh_updated;
638 gslpp::complex C_8Lh_cache;
640 unsigned int C_L_nunu_e_updated;
641 unsigned int C_L_nunu_mu_updated;
642 unsigned int C_L_nunu_tau_updated;
643 gslpp::complex C_L_nunu_e_cache;
644 gslpp::complex C_L_nunu_mu_cache;
645 gslpp::complex C_L_nunu_tau_cache;
647 unsigned int C_R_nunu_e_updated;
648 unsigned int C_R_nunu_mu_updated;
649 unsigned int C_R_nunu_tau_updated;
650 gslpp::complex C_R_nunu_e_cache;
651 gslpp::complex C_R_nunu_mu_cache;
652 gslpp::complex C_R_nunu_tau_cache;
654 unsigned int Yupdated;
655 gslpp::vector<double> Ycache;
657 unsigned int H_V0updated;
658 gslpp::vector<double> H_V0cache;
659 gslpp::complex H_V0Ccache[3];
660 gslpp::complex H_V0Ccache_dispersion[4];
662 unsigned int H_A0updated;
664 unsigned int H_Supdated;
665 gslpp::vector<double> H_Scache;
667 unsigned int H_P_updated;
668 gslpp::vector<double> H_P_cache;
670 unsigned int I0_updated;
671 unsigned int I2_updated;
672 unsigned int I8_updated;
674 unsigned int Itree_updated;
675 gslpp::vector<double> Itree_cache;
677 std::map<std::pair<double, double>, unsigned int > I1Cached;
679 std::map<std::pair<double, double>, unsigned int > sigma0Cached;
680 std::map<std::pair<double, double>, unsigned int > sigma2Cached;
681 std::map<std::pair<double, double>, unsigned int > sigma8Cached;
683 std::map<std::pair<double, double>, unsigned int > delta0Cached;
684 std::map<std::pair<double, double>, unsigned int > delta2Cached;
686 std::map<std::pair<double, double>, unsigned int > sigmaTreeCached;
688 double avaSigma;
689 double avaDelta;
690 double avaSigmaTree;
691 double avaDTPPR;
693 double errSigma;
694 double errDelta;
695 double errSigmaTree;
696 double errDTPPR;
698 gsl_function FS;
699 gsl_function FD;
700 gsl_function DTPPR;
702 gsl_integration_cquad_workspace * w_sigma;
703 gsl_integration_cquad_workspace * w_delta;
704 gsl_integration_cquad_workspace * w_sigmaTree;
705 gsl_integration_cquad_workspace * w_DTPPR;
707 gsl_error_handler_t * old_handler;
709 std::map<std::pair<double, double>, gslpp::complex > cacheI1;
711 std::map<std::pair<double, double>, double > cacheSigma0;
712 std::map<std::pair<double, double>, double > cacheSigma2;
713 std::map<std::pair<double, double>, double > cacheSigma8;
715 std::map<std::pair<double, double>, double > cacheDelta0;
716 std::map<std::pair<double, double>, double > cacheDelta2;
718 std::map<std::pair<double, double>, double > cacheSigmaTree;
720 std::map<double, unsigned int> deltaTparpCached;
721 std::map<double, unsigned int> deltaTparmCached;
723 std::map<double, gslpp::complex> cacheDeltaTparp;
724 std::map<double, gslpp::complex> cacheDeltaTparm;
726 unsigned int deltaTparpupdated;
727 unsigned int deltaTparmupdated;
729 unsigned int T_updated;
730 gslpp::vector<double> T_cache;
737 void updateParameters();
738
742 void checkCache();
743
752 double LCSR_fit1(double q2, double r_1, double r_2, double m_fit2);
753
754
762 double LCSR_fit2(double q2, double r_2, double m_fit2);
763
764
774 double LCSR_fit3(double q2, double b_0, double b_1, double b_2, double m_fit2);
775
776
782 double zeta(double q2);
783
784
790 double zeta_DM(double q2);
791
798 double phiplus_DM(double q2, double m_fit2);
799
806 double phi0_DM(double q2);
807
814 double phiT_DM(double q2, double m_fit2);
815
825 double fplus_DM(double q2, double b_0, double b_1, double b_2, double m_fit2);
826
836 double f0_DM(double q2, double b_0, double b_1, double b_2);
837
847 double fT_DM(double q2, double b_0, double b_1, double b_2, double m_fit2);
848
849
859 double LATTICE_fit1(double q2, double b_0, double b_1, double b_2, double m_fit2);
860
861
871 double LATTICE_fit2(double q2, double b_0, double b_1, double b_2, double m_fit2);
872
873
879 double f_plus(double q2);
880
881
887 double f_T(double q2);
888
889
895 double f_0(double q2);
896
897
903 gslpp::complex H_0(double q2);
904
911 gslpp::complex H_c(double q2, double mu2);
912
919 gslpp::complex H_b(double q2, double mu2);
920
921
927 gslpp::complex Y(double q2);
928
929
935 double k2 (double q2);
936
937
943 double beta (double q2);
944
950 double beta2 (double q2);
951
952
958 double lambda(double q2);
959
960
966 double F(double q2);
967
973 double I_1c(double q2);
974
980 double I_2c(double q2);
981
987 double I_6c(double q2);
988
989
996 double Delta(int i, double q2);
997
1003 double SigmaTree(double q2);
1004
1009 double getintegratedSigmaTree();
1010
1016 double getSigma1c(double q2)
1017 {
1018 return I_1c(q2);
1019 };
1020
1026 double getSigma2c(double q2)
1027 {
1028 return I_2c(q2);
1029 };
1030
1036 double getSigma6c(double q2)
1037 {
1038 return I_6c(q2);
1039 };
1040
1046 double getDelta1c(double q2)
1047 {
1048 return Delta(0, q2);
1049 };
1050
1056 double getDelta2c(double q2)
1057 {
1058 return Delta(2, q2);
1059 };
1060
1067 gslpp::complex I1(double u, double q2);
1068
1075 gslpp::complex Tparplus(double u, double q2);
1076
1083 gslpp::complex Tparminus(double u, double q2);
1084
1090 double Integrand_ReTparplus(double up);
1091
1097 double Integrand_ImTparplus(double up);
1098
1104 double Integrand_ReTparminus(double up);
1105
1111 double Integrand_ImTparminus(double up);
1112
1118 double Integrand_ReTpar_pm(double up);
1119
1125 double Integrand_ImTpar_pm(double up);
1126
1132 gslpp::complex F19(double q2);
1133
1139 gslpp::complex F27(double q2);
1140
1146 gslpp::complex F29(double q2);
1147
1153 gslpp::complex F87(double q2);
1154
1160 double F89(double q2);
1161
1167 gslpp::complex Cpar(double q2);
1168
1174 gslpp::complex deltaTpar(double q2);
1175
1182 double reDC9fit(double* x, double* p);
1183
1190 double imDC9fit(double* x, double* p);
1191
1195 void fit_DeltaC9_mumu();
1196
1202 gslpp::complex fDeltaC9(double q2);
1203
1209 gslpp::complex DeltaC9(double q2);
1210
1216 gslpp::complex deltaC7_QCDF(double q2, bool spline = false);
1217
1223 gslpp::complex deltaC9_QCDF(double q2, bool spline = false);
1224
1225 void spline_QCDF_func();
1226
1227};
1228
1229#endif /* MPLL_H */
F_1
Definition: F_1.h:15
F_2
Definition: F_2.h:15
MPll
A class for the decay.
Definition: MPll.h:173
MPll::initializeMPllParameters
std::vector< std::string > initializeMPllParameters()
A method for initializing the parameters necessary for MPll.
Definition: MPll.cpp:80
MPll::Integrand_ReTpar_pm
double Integrand_ReTpar_pm(double up)
The sum of Integrand_ReTparplus() and Integrand_ReTparminus().
Definition: MPll.cpp:1058
MPll::getSigma1c
double getSigma1c(double q2)
The CP average .
Definition: MPll.h:1016
MPll::FixedWCbtos
bool FixedWCbtos
Definition: MPll.h:323
MPll::SigmaTree
double SigmaTree(double q2)
Definition: MPll.cpp:1559
MPll::MPll_GRvDV_flag
bool MPll_GRvDV_flag
Definition: MPll.h:325
MPll::H_A
gslpp::complex H_A(double q2)
The helicity amplitude .
Definition: MPll.cpp:1345
MPll::Mb
double Mb
Definition: MPll.h:335
MPll::h_lambda
gslpp::complex h_lambda(double q2)
The non-pertubative ccbar contributions to the helicity amplitudes.
Definition: MPll.cpp:1321
MPll::F89
double F89(double q2)
The correction from .
Definition: MPll.cpp:1102
MPll::mc_pole
double mc_pole
Definition: MPll.h:340
MPll::deltaC9_QCDF
gslpp::complex deltaC9_QCDF(double q2, bool spline=false)
QCDF Correction from various BFS papers (hep-ph/0403185, hep-ph/0412400) and Greub et....
Definition: MPll.cpp:1223
MPll::DeltaC9_KD
gslpp::complex DeltaC9_KD(double q2)
Definition: MPll.cpp:1315
MPll::H_P
gslpp::complex H_P(double q2)
The helicity amplitude .
Definition: MPll.cpp:1355
MPll::pseudoscalar
QCD::meson pseudoscalar
Definition: MPll.h:318
MPll::H_S
gslpp::complex H_S(double q2)
The helicity amplitude .
Definition: MPll.cpp:1350
MPll::F87
gslpp::complex F87(double q2)
The correction from .
Definition: MPll.cpp:1095
MPll::mpllParameters
std::vector< std::string > mpllParameters
Definition: MPll.h:319
MPll::Integrand_ImTparminus
double Integrand_ImTparminus(double up)
The imaginary part of the integral involving at fixed , according to .
Definition: MPll.cpp:1042
MPll::myF_2
std::unique_ptr< F_2 > myF_2
Definition: MPll.h:321
MPll::MPll_Lattice_flag
bool MPll_Lattice_flag
Definition: MPll.h:324
MPll::Mlep
double Mlep
Definition: MPll.h:332
MPll::mu_h
double mu_h
Definition: MPll.h:337
MPll::I1
gslpp::complex I1(double u, double q2)
The function from .
Definition: MPll.cpp:973
MPll::Integrand_ReTparminus
double Integrand_ReTparminus(double up)
The real part of the integral involving at fixed , according to .
Definition: MPll.cpp:1031
MPll::MPll_DM_flag
bool MPll_DM_flag
Definition: MPll.h:326
MPll::fit_DeltaC9_mumu
void fit_DeltaC9_mumu()
The fitting routine for the QCDF correction in the muon channel.
Definition: MPll.cpp:1152
MPll::Tparminus
gslpp::complex Tparminus(double u, double q2)
The function from .
Definition: MPll.cpp:1008
MPll::imDC9fit
double imDC9fit(double *x, double *p)
The fit function for the imaginary part of the QCDF correction .
Definition: MPll.cpp:1145
MPll::getDelta2c
double getDelta2c(double q2)
The CP asymmetry .
Definition: MPll.h:1056
MPll::meson
QCD::meson meson
Definition: MPll.h:317
MPll::width
double width
Definition: MPll.h:343
MPll::F19
gslpp::complex F19(double q2)
The correction from .
Definition: MPll.cpp:1063
MPll::integrateSigma
double integrateSigma(int i, double q_min, double q_max)
The integral of from to .
Definition: MPll.cpp:1425
MPll::Cpar
gslpp::complex Cpar(double q2)
The correction from .
Definition: MPll.cpp:1109
MPll::reDC9fit
double reDC9fit(double *x, double *p)
The fit function for the real part of the QCDF correction .
Definition: MPll.cpp:1140
MPll::Mc
double Mc
Definition: MPll.h:338
MPll::mySM
const StandardModel & mySM
Definition: MPll.h:315
MPll::V_L
gslpp::complex V_L(double q2)
The helicity form factor .
Definition: MPll.cpp:954
MPll::H_nunu
gslpp::complex H_nunu(double q2, QCD::lepton lep)
The helicity amplitude for the invisible decay.
Definition: MPll.cpp:1360
MPll::getSigma
double getSigma(int i, double q_2)
The value of from to .
Definition: MPll.cpp:1471
MPll::lep
QCD::lepton lep
Definition: MPll.h:316
MPll::~MPll
virtual ~MPll()
Destructor.
Definition: MPll.cpp:76
MPll::funct_g
gslpp::complex funct_g(double q2)
Definition: MPll.cpp:1307
MPll::mJ2
double mJ2
Definition: MPll.h:328
MPll::GF
double GF
Definition: MPll.h:330
MPll::Integrand_ImTpar_pm
double Integrand_ImTpar_pm(double up)
The sum of Integrand_ImTparplus() and Integrand_ImTparminus().
Definition: MPll.cpp:1053
MPll::Integrand_ReTparplus
double Integrand_ReTparplus(double up)
The real part of the integral involving at fixed , according to .
Definition: MPll.cpp:1015
MPll::T_L
gslpp::complex T_L(double q2)
The helicity form factor .
Definition: MPll.cpp:959
MPll::alpha_s_mub
double alpha_s_mub
Definition: MPll.h:344
MPll::getSigma2c
double getSigma2c(double q2)
The CP average .
Definition: MPll.h:1026
MPll::Tparplus
gslpp::complex Tparplus(double u, double q2)
The function from .
Definition: MPll.cpp:995
MPll::MM
double MM
Definition: MPll.h:333
MPll::myF_1
std::unique_ptr< F_1 > myF_1
Definition: MPll.h:320
MPll::NeutrinoTree_flag
bool NeutrinoTree_flag
Definition: MPll.h:327
MPll::integrateSigmaTree
double integrateSigmaTree(double q_min, double q_max)
The integral of from to (arxiv/2301.06990)
Definition: MPll.cpp:1529
MPll::getintegratedSigmaTree
double getintegratedSigmaTree()
The integral of from 0 to .
Definition: MPll.cpp:1564
MPll::Ms
double Ms
Definition: MPll.h:341
MPll::getDelta1c
double getDelta1c(double q2)
The CP asymmetry .
Definition: MPll.h:1046
MPll::MP
double MP
Definition: MPll.h:334
MPll::H_V
gslpp::complex H_V(double q2)
The helicity amplitude .
Definition: MPll.cpp:1338
MPll::DeltaC9
gslpp::complex DeltaC9(double q2)
The total QCDF correction computed integrating over .
Definition: MPll.cpp:1193
MPll::deltaC7_QCDF
gslpp::complex deltaC7_QCDF(double q2, bool spline=false)
QCDF Correction from various BFS papers (hep-ph/0403185, hep-ph/0412400) and Greub et....
Definition: MPll.cpp:1198
MPll::dispersion
bool dispersion
Definition: MPll.h:322
MPll::getSigma6c
double getSigma6c(double q2)
The CP average .
Definition: MPll.h:1036
MPll::mb_pole
double mb_pole
Definition: MPll.h:339
MPll::F27
gslpp::complex F27(double q2)
The correction from .
Definition: MPll.cpp:1077
MPll::ale
double ale
Definition: MPll.h:331
MPll::getwidth
double getwidth()
The width of the meson M.
Definition: MPll.h:302
MPll::mu_b
double mu_b
Definition: MPll.h:336
MPll::fDeltaC9
gslpp::complex fDeltaC9(double q2)
The total QCDF correction computed fitting over .
Definition: MPll.cpp:1184
MPll::S_L
double S_L(double q2)
The helicity form factor .
Definition: MPll.cpp:964
MPll::spectator_charge
double spectator_charge
Definition: MPll.h:342
MPll::MPll
MPll(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i)
Constructor.
Definition: MPll.cpp:20
MPll::deltaTpar
gslpp::complex deltaTpar(double q2)
The total correction from .
Definition: MPll.cpp:1115
MPll::integrateDelta
double integrateDelta(int i, double q_min, double q_max)
The integral of from to .
Definition: MPll.cpp:1492
MPll::F29
gslpp::complex F29(double q2)
The correction from .
Definition: MPll.cpp:1086
MPll::spline_QCDF_func
void spline_QCDF_func()
Definition: MPll.cpp:1248
MPll::Integrand_ImTparplus
double Integrand_ImTparplus(double up)
The imaginary part of the integral involving at fixed , according to .
Definition: MPll.cpp:1023
QCD::meson
meson
An enum type for mesons.
Definition: QCD.h:336
QCD::lepton
lepton
An enum type for leptons.
Definition: QCD.h:310
StandardModel
A model class for the Standard Model.
Definition: StandardModel.h:521
h_0
A class for the correction in .
Definition: MVllObservables.h:1991