a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Loading...
Searching...
No Matches
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 double Chi1minus;
375 double Chi0plus;
376 double ChiTT;
377 /*LCSR fit parameters*/
378 double r_1_fplus;
379 double r_2_fplus;
380 double m_fit2_fplus;
381 double r_1_fT;
382 double r_2_fT;
383 double m_fit2_fT;
384 double r_2_f0;
385 double m_fit2_f0;
388 gslpp::vector<gslpp::complex> ** allcoeff;
389 gslpp::vector<gslpp::complex> ** allcoeffh;
390 gslpp::vector<gslpp::complex> ** allcoeffprime;
392 gslpp::vector<gslpp::complex> ** allcoeff_noSM;
394 gslpp::vector<gslpp::complex> ** allcoeff_nu;
396 gslpp::vector<gslpp::complex> ** allcoeff_noSM_nu;
398 gslpp::complex C_1;
399 gslpp::complex C_1L_bar;
400 gslpp::complex C_1Lh_bar;
401 gslpp::complex C_2;
402 gslpp::complex C_2L_bar;
403 gslpp::complex C_2Lh_bar;
404 gslpp::complex C_3;
405 gslpp::complex C_4;
406 gslpp::complex C_5;
407 gslpp::complex C_6;
408 gslpp::complex C_7;
409 gslpp::complex C_8;
410 gslpp::complex C_8L;
411 gslpp::complex C_8Lh;
412 gslpp::complex C_9;
413 gslpp::complex C_10;
414 gslpp::complex C_S;
415 gslpp::complex C_P;
417 gslpp::complex C_7p;
418 gslpp::complex C_9p;
419 gslpp::complex C_10p;
420 gslpp::complex C_Sp;
421 gslpp::complex C_Pp;
423 gslpp::complex C_L_nunu_e;
424 gslpp::complex C_R_nunu_e;
425 gslpp::complex C_L_nunu_mu;
426 gslpp::complex C_R_nunu_mu;
427 gslpp::complex C_L_nunu_tau;
428 gslpp::complex C_R_nunu_tau;
430 gsl_interp_accel *acc_Re_deltaC7_QCDF;
431 gsl_interp_accel *acc_Im_deltaC7_QCDF;
432 gsl_interp_accel *acc_Re_deltaC9_QCDF;
433 gsl_interp_accel *acc_Im_deltaC9_QCDF;
434
435 gsl_spline *spline_Re_deltaC7_QCDF;
436 gsl_spline *spline_Im_deltaC7_QCDF;
437 gsl_spline *spline_Re_deltaC9_QCDF;
438 gsl_spline *spline_Im_deltaC9_QCDF;
439
440
441 std::vector<double> ReDeltaC9;
442 std::vector<double> ImDeltaC9;
443 std::vector<double> myq2;
444 TFitResultPtr refres;
445 TFitResultPtr imfres;
447 TGraph gr1;
448 TGraph gr2;
450 TF1 reffit;
451 TF1 imffit;
453 double tmpq2;
455 gslpp::complex H_0_pre;
456 gslpp::complex H_0_WC;
457 gslpp::complex H_c_WC;
458 gslpp::complex H_b_WC;
460 gslpp::complex ihalfMPI;
461 double fournineth;
462 double half;
463 double twothird;
464 double Mc2;
465 double Mb2;
466 double logMc;
467 double logMb;
468 double mu_b2;
469 double fourMc2;
470 double fourMb2;
471 double Mlep2;
472 double NN;
473 double MM2;
474 double MM4;
475 double MP2;
476 double MP4;
477 double MM2mMP2;
478 double twoMP2;
479 double twoMM;
480 double twoMM2;
481 double twoMM2_MMpMP;
482 double twoMM_MbpMs;
483 double S_L_pre;
484 double fourMM2;
485 double twoMboMM;
486 double sixteenM_PI2;
487 double ninetysixM_PI3MM3;
488 double MboMW;
489 double MboMM;
490 double MsoMb;
491 double twoMlepMb;
492 double threeGegen0;
493 double threeGegen1otwo;
494 double M_PI2osix;
495 double twoMc2;
496 double sixMMoMb;
497 double CF;
498 double deltaT_0;
499 double deltaT_1par;
501 gslpp::complex ubar;
502 gslpp::complex arg1;
503 gslpp::complex B01;
504 gslpp::complex B00;
505 gslpp::complex xp;
506 gslpp::complex xm;
507 gslpp::complex yp;
508 gslpp::complex ym;
509 gslpp::complex L1xp;
510 gslpp::complex L1xm;
511 gslpp::complex L1yp;
512 gslpp::complex L1ym;
513 gslpp::complex F87_0;
514 gslpp::complex F87_1;
515 gslpp::complex F87_2;
516 gslpp::complex F87_3;
517 gslpp::complex F29_0;
518 gslpp::complex F29_L1;
519 gslpp::complex F29_1;
520 gslpp::complex F29_2;
521 gslpp::complex F29_3;
522 gslpp::complex F29_L1_1;
523 gslpp::complex F29_L1_2;
524 gslpp::complex F29_L1_3;
525 gslpp::complex F27_0;
526 gslpp::complex F27_1;
527 gslpp::complex F27_2;
528 gslpp::complex F27_3;
529 gslpp::complex F27_L1_1;
530 gslpp::complex F27_L1_2;
531 gslpp::complex F27_L1_3;
532 double F89_0;
533 double F89_1;
534 double F89_2;
535 double F89_3;
536 double Ee;
538 //additional variables for B to K nu nu
539 const double M_PI2 = M_PI * M_PI;
540 double GF4;
541 double MM3;
542 double fM2;
543 double fP2;
544
545 double mtau;
546 double mtau2;
547 double Gammatau;
548 double VusVub_abs2;
549
550 unsigned int fplus_lat_updated;
551 gslpp::vector<double> fplus_lat_cache;
553 unsigned int fT_lat_updated;
554 gslpp::vector<double> fT_lat_cache;
556 unsigned int f0_lat_updated;
557 gslpp::vector<double> f0_lat_cache;
559 unsigned int fplus_updated;
560 gslpp::vector<double> fplus_cache;
562 unsigned int fT_updated;
563 gslpp::vector<double> fT_cache;
565 unsigned int f0_updated;
566 double f0_cache;
568 unsigned int k2_updated;
569 gslpp::vector<double> k2_cache;
571 unsigned int beta_updated;
572 double beta_cache;
574 unsigned int lambda_updated;
576 unsigned int F_updated;
578 unsigned int VL_updated;
580 unsigned int TL_updated;
582 unsigned int SL_updated;
583 gslpp::vector<double> SL_cache;
585 unsigned int N_updated;
586 gslpp::vector<double> N_cache;
587 gslpp::complex Nc_cache;
589 unsigned int C_1_updated;
590 gslpp::complex C_1_cache;
592 unsigned int C_2_updated;
593 gslpp::complex C_2_cache;
595 unsigned int C_3_updated;
596 gslpp::complex C_3_cache;
598 unsigned int C_4_updated;
599 gslpp::complex C_4_cache;
601 unsigned int C_5_updated;
602 gslpp::complex C_5_cache;
604 unsigned int C_6_updated;
605 gslpp::complex C_6_cache;
607 unsigned int C_7_updated;
608 gslpp::complex C_7_cache;
610 unsigned int C_9_updated;
611 gslpp::complex C_9_cache;
613 unsigned int C_10_updated;
614 gslpp::complex C_10_cache;
616 unsigned int C_7p_updated;
617 gslpp::complex C_7p_cache;
619 unsigned int C_9p_updated;
620 gslpp::complex C_9p_cache;
622 unsigned int C_10p_updated;
623 gslpp::complex C_10p_cache;
625 unsigned int C_S_updated;
626 gslpp::complex C_S_cache;
628 unsigned int C_P_updated;
629 gslpp::complex C_P_cache;
631 unsigned int C_Sp_updated;
632 gslpp::complex C_Sp_cache;
634 unsigned int C_Pp_updated;
635 gslpp::complex C_Pp_cache;
637 unsigned int C_2Lh_updated;
638 gslpp::complex C_2Lh_cache;
640 unsigned int C_8Lh_updated;
641 gslpp::complex C_8Lh_cache;
643 unsigned int C_L_nunu_e_updated;
644 unsigned int C_L_nunu_mu_updated;
645 unsigned int C_L_nunu_tau_updated;
646 gslpp::complex C_L_nunu_e_cache;
647 gslpp::complex C_L_nunu_mu_cache;
648 gslpp::complex C_L_nunu_tau_cache;
650 unsigned int C_R_nunu_e_updated;
651 unsigned int C_R_nunu_mu_updated;
652 unsigned int C_R_nunu_tau_updated;
653 gslpp::complex C_R_nunu_e_cache;
654 gslpp::complex C_R_nunu_mu_cache;
655 gslpp::complex C_R_nunu_tau_cache;
657 unsigned int Yupdated;
658 gslpp::vector<double> Ycache;
660 unsigned int H_V0updated;
661 gslpp::vector<double> H_V0cache;
662 gslpp::complex H_V0Ccache[3];
663 gslpp::complex H_V0Ccache_dispersion[4];
665 unsigned int H_A0updated;
667 unsigned int H_Supdated;
668 gslpp::vector<double> H_Scache;
670 unsigned int H_P_updated;
671 gslpp::vector<double> H_P_cache;
673 unsigned int I0_updated;
674 unsigned int I2_updated;
675 unsigned int I8_updated;
677 unsigned int Itree_updated;
678 gslpp::vector<double> Itree_cache;
680 std::map<std::pair<double, double>, unsigned int > I1Cached;
682 std::map<std::pair<double, double>, unsigned int > sigma0Cached;
683 std::map<std::pair<double, double>, unsigned int > sigma2Cached;
684 std::map<std::pair<double, double>, unsigned int > sigma8Cached;
686 std::map<std::pair<double, double>, unsigned int > delta0Cached;
687 std::map<std::pair<double, double>, unsigned int > delta2Cached;
689 std::map<std::pair<double, double>, unsigned int > sigmaTreeCached;
691 double avaSigma;
692 double avaDelta;
693 double avaSigmaTree;
694 double avaDTPPR;
696 double errSigma;
697 double errDelta;
698 double errSigmaTree;
699 double errDTPPR;
701 gsl_function FS;
702 gsl_function FD;
703 gsl_function DTPPR;
705 gsl_integration_cquad_workspace * w_sigma;
706 gsl_integration_cquad_workspace * w_delta;
707 gsl_integration_cquad_workspace * w_sigmaTree;
708 gsl_integration_cquad_workspace * w_DTPPR;
710 gsl_error_handler_t * old_handler;
712 std::map<std::pair<double, double>, gslpp::complex > cacheI1;
714 std::map<std::pair<double, double>, double > cacheSigma0;
715 std::map<std::pair<double, double>, double > cacheSigma2;
716 std::map<std::pair<double, double>, double > cacheSigma8;
718 std::map<std::pair<double, double>, double > cacheDelta0;
719 std::map<std::pair<double, double>, double > cacheDelta2;
721 std::map<std::pair<double, double>, double > cacheSigmaTree;
723 std::map<double, unsigned int> deltaTparpCached;
724 std::map<double, unsigned int> deltaTparmCached;
726 std::map<double, gslpp::complex> cacheDeltaTparp;
727 std::map<double, gslpp::complex> cacheDeltaTparm;
729 unsigned int deltaTparpupdated;
730 unsigned int deltaTparmupdated;
732 unsigned int T_updated;
733 gslpp::vector<double> T_cache;
740 void updateParameters();
741
745 void checkCache();
746
755 double LCSR_fit1(double q2, double r_1, double r_2, double m_fit2);
756
757
765 double LCSR_fit2(double q2, double r_2, double m_fit2);
766
767
777 double LCSR_fit3(double q2, double b_0, double b_1, double b_2, double m_fit2);
778
779
785 double zeta(double q2);
786
787
793 double zeta_DM(double q2);
794
801 double phiplus_DM(double q2, double m_fit2);
802
809 double phi0_DM(double q2);
810
817 double phiT_DM(double q2, double m_fit2);
818
828 double fplus_DM(double q2, double b_0, double b_1, double b_2, double m_fit2);
829
839 double f0_DM(double q2, double b_0, double b_1, double b_2);
840
850 double fT_DM(double q2, double b_0, double b_1, double b_2, double m_fit2);
851
852
862 double LATTICE_fit1(double q2, double b_0, double b_1, double b_2, double m_fit2);
863
864
874 double LATTICE_fit2(double q2, double b_0, double b_1, double b_2, double m_fit2);
875
876
882 double f_plus(double q2);
883
884
890 double f_T(double q2);
891
892
898 double f_0(double q2);
899
900
906 gslpp::complex H_0(double q2);
907
914 gslpp::complex H_c(double q2, double mu2);
915
922 gslpp::complex H_b(double q2, double mu2);
923
924
930 gslpp::complex Y(double q2);
931
932
938 double k2 (double q2);
939
940
946 double beta (double q2);
947
953 double beta2 (double q2);
954
955
961 double lambda(double q2);
962
963
969 double F(double q2);
970
976 double I_1c(double q2);
977
983 double I_2c(double q2);
984
990 double I_6c(double q2);
991
992
999 double Delta(int i, double q2);
1000
1006 double SigmaTree(double q2);
1007
1012 double getintegratedSigmaTree();
1013
1019 double getSigma1c(double q2)
1020 {
1021 return I_1c(q2);
1022 };
1023
1029 double getSigma2c(double q2)
1030 {
1031 return I_2c(q2);
1032 };
1033
1039 double getSigma6c(double q2)
1040 {
1041 return I_6c(q2);
1042 };
1043
1049 double getDelta1c(double q2)
1050 {
1051 return Delta(0, q2);
1052 };
1053
1059 double getDelta2c(double q2)
1060 {
1061 return Delta(2, q2);
1062 };
1063
1070 gslpp::complex I1(double u, double q2);
1071
1078 gslpp::complex Tparplus(double u, double q2);
1079
1086 gslpp::complex Tparminus(double u, double q2);
1087
1093 double Integrand_ReTparplus(double up);
1094
1100 double Integrand_ImTparplus(double up);
1101
1107 double Integrand_ReTparminus(double up);
1108
1114 double Integrand_ImTparminus(double up);
1115
1121 double Integrand_ReTpar_pm(double up);
1122
1128 double Integrand_ImTpar_pm(double up);
1129
1135 gslpp::complex F19(double q2);
1136
1142 gslpp::complex F27(double q2);
1143
1149 gslpp::complex F29(double q2);
1150
1156 gslpp::complex F87(double q2);
1157
1163 double F89(double q2);
1164
1170 gslpp::complex Cpar(double q2);
1171
1177 gslpp::complex deltaTpar(double q2);
1178
1185 double reDC9fit(double* x, double* p);
1186
1193 double imDC9fit(double* x, double* p);
1194
1198 void fit_DeltaC9_mumu();
1199
1205 gslpp::complex fDeltaC9(double q2);
1206
1212 gslpp::complex DeltaC9(double q2);
1213
1219 gslpp::complex deltaC7_QCDF(double q2, bool spline = false);
1220
1226 gslpp::complex deltaC9_QCDF(double q2, bool spline = false);
1227
1228 void spline_QCDF_func();
1229
1230};
1231
1232#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:1060
MPll::getSigma1c
double getSigma1c(double q2)
The CP average .
Definition MPll.h:1019
MPll::FixedWCbtos
bool FixedWCbtos
Definition MPll.h:323
MPll::SigmaTree
double SigmaTree(double q2)
Definition MPll.cpp:1561
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:1347
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:1323
MPll::F89
double F89(double q2)
The correction from .
Definition MPll.cpp:1104
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:1225
MPll::DeltaC9_KD
gslpp::complex DeltaC9_KD(double q2)
Definition MPll.cpp:1317
MPll::H_P
gslpp::complex H_P(double q2)
The helicity amplitude .
Definition MPll.cpp:1357
MPll::pseudoscalar
QCD::meson pseudoscalar
Definition MPll.h:318
MPll::H_S
gslpp::complex H_S(double q2)
The helicity amplitude .
Definition MPll.cpp:1352
MPll::F87
gslpp::complex F87(double q2)
The correction from .
Definition MPll.cpp:1097
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:1044
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:975
MPll::Integrand_ReTparminus
double Integrand_ReTparminus(double up)
The real part of the integral involving at fixed , according to .
Definition MPll.cpp:1033
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:1154
MPll::Tparminus
gslpp::complex Tparminus(double u, double q2)
The function from .
Definition MPll.cpp:1010
MPll::imDC9fit
double imDC9fit(double *x, double *p)
The fit function for the imaginary part of the QCDF correction .
Definition MPll.cpp:1147
MPll::getDelta2c
double getDelta2c(double q2)
The CP asymmetry .
Definition MPll.h:1059
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:1065
MPll::integrateSigma
double integrateSigma(int i, double q_min, double q_max)
The integral of from to .
Definition MPll.cpp:1427
MPll::Cpar
gslpp::complex Cpar(double q2)
The correction from .
Definition MPll.cpp:1111
MPll::reDC9fit
double reDC9fit(double *x, double *p)
The fit function for the real part of the QCDF correction .
Definition MPll.cpp:1142
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:956
MPll::H_nunu
gslpp::complex H_nunu(double q2, QCD::lepton lep)
The helicity amplitude for the invisible decay.
Definition MPll.cpp:1362
MPll::getSigma
double getSigma(int i, double q_2)
The value of from to .
Definition MPll.cpp:1473
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:1309
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:1055
MPll::Integrand_ReTparplus
double Integrand_ReTparplus(double up)
The real part of the integral involving at fixed , according to .
Definition MPll.cpp:1017
MPll::T_L
gslpp::complex T_L(double q2)
The helicity form factor .
Definition MPll.cpp:961
MPll::alpha_s_mub
double alpha_s_mub
Definition MPll.h:344
MPll::getSigma2c
double getSigma2c(double q2)
The CP average .
Definition MPll.h:1029
MPll::Tparplus
gslpp::complex Tparplus(double u, double q2)
The function from .
Definition MPll.cpp:997
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:1531
MPll::getintegratedSigmaTree
double getintegratedSigmaTree()
The integral of from 0 to .
Definition MPll.cpp:1566
MPll::Ms
double Ms
Definition MPll.h:341
MPll::getDelta1c
double getDelta1c(double q2)
The CP asymmetry .
Definition MPll.h:1049
MPll::MP
double MP
Definition MPll.h:334
MPll::H_V
gslpp::complex H_V(double q2)
The helicity amplitude .
Definition MPll.cpp:1340
MPll::DeltaC9
gslpp::complex DeltaC9(double q2)
The total QCDF correction computed integrating over .
Definition MPll.cpp:1195
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:1200
MPll::dispersion
bool dispersion
Definition MPll.h:322
MPll::getSigma6c
double getSigma6c(double q2)
The CP average .
Definition MPll.h:1039
MPll::mb_pole
double mb_pole
Definition MPll.h:339
MPll::F27
gslpp::complex F27(double q2)
The correction from .
Definition MPll.cpp:1079
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:1186
MPll::S_L
double S_L(double q2)
The helicity form factor .
Definition MPll.cpp:966
MPll::spectator_charge
double spectator_charge
Definition MPll.h:342
MPll::deltaTpar
gslpp::complex deltaTpar(double q2)
The total correction from .
Definition MPll.cpp:1117
MPll::integrateDelta
double integrateDelta(int i, double q_min, double q_max)
The integral of from to .
Definition MPll.cpp:1494
MPll::F29
gslpp::complex F29(double q2)
The correction from .
Definition MPll.cpp:1088
MPll::spline_QCDF_func
void spline_QCDF_func()
Definition MPll.cpp:1250
MPll::Integrand_ImTparplus
double Integrand_ImTparplus(double up)
The imaginary part of the integral involving at fixed , according to .
Definition MPll.cpp:1025
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:2030