(Generated on Fri Feb 19 2016 13:25:07 by doxygen 1.8.9.1)
MVll.h
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2014 HEPfit Collaboration
3  * All rights reserved.
4  *
5  * For the licensing terms see doc/COPYING.
6  */
7 
8 #ifndef MVLL_H
9 #define MVLL_H
10 
11 #include "StandardModel.h"
12 #include "ThObservable.h"
13 #include <math.h>
14 #include <gsl/gsl_math.h>
15 #include <gsl/gsl_integration.h>
16 #include <assert.h>
17 #include <gsl/gsl_monte_plain.h>
18 #include <TF1.h>
19 #include <TGraph.h>
20 #include <TFitResultPtr.h>
21 
22 #define SWITCH 8.2
23 
24 /*******************************************************************************
25  * GSL Function Conversion BEGIN *
26  * ****************************************************************************/
27 
28 template<class F>
29 static double gslFunctionAdapter( double x, void* p)
30 {
31  // Here I do recover the "right" pointer, safer to use static_cast
32  // than reinterpret_cast.
33  F* function = static_cast<F*>( p );
34  return (*function)( x );
35 }
36 
37 template<class F>
38 gsl_function convertToGslFunction( const F& f )
39 {
40  gsl_function gslFunction;
41 
42  const void* p = &f;
43  assert (p != 0);
44 
45  gslFunction.function = &gslFunctionAdapter<F>;
46  // Just to eliminate the const.
47  gslFunction.params = const_cast<void*>( p );
48 
49  return gslFunction;
50 }
51 
52 /*******************************************************************************
53  * GSL Function conversion END *
54  * ****************************************************************************/
55 
87 class MVll {
88 public:
89 
97  MVll(const StandardModel& SM_i, StandardModel::meson meson_i, StandardModel::meson vector_i, StandardModel::lepton lep_i);
98 
102  virtual ~MVll();
103 
111  double integrateSigma(int i, double q_min, double q_max);
112 
119  double getSigma(int i, double q_2);
120 
128  double integrateDelta(int i, double q_min, double q_max);
129 
134  double getwidth(){
135  updateParameters();
136  return width;
137  }
138 
143  double getMlep(){
144  updateParameters();
145  return Mlep;
146  }
147 
153  double beta (double q2);
154 
160  double getV0(double q2)
161  {
162  updateParameters();
163  return (2. * MM * sqrt(q2))/sqrt(lambda(q2)) * V_0t(q2);
164  };
165 
171  double getVp(double q2)
172  {
173  updateParameters();
174  return V_p(q2);
175  };
176 
182  double getVm(double q2)
183  {
184  return V_m(q2);
185  };
186 
192  double getT0(double q2)
193  {
194  updateParameters();
195  return twoMM3/sqrt(q2 * lambda(q2)) * T_0t(q2);
196  };
197 
203  double getTp(double q2)
204  {
205  updateParameters();
206  return T_p(q2);
207  };
208 
214  double getTm(double q2)
215  {
216  updateParameters();
217  return T_m(q2);
218  };
219 
225  double getS(double q2)
226  {
227  updateParameters();
228  return S_L_pre/sqrt(lambda(q2)) * S_L(q2);
229  };
230 
236  gslpp::complex H_V_0(double q2);
237 
243  gslpp::complex H_V_p(double q2);
244 
250  gslpp::complex H_V_m(double q2);
251 
257  gslpp::complex H_A_0(double q2);
258 
264  gslpp::complex H_A_p(double q2);
265 
271  gslpp::complex H_A_m(double q2);
272 
278  gslpp::complex H_S(double q2);
279 
285  gslpp::complex H_P(double q2);
286 
292  double getgtilde_1_re(double q2)
293  {
294  updateParameters();
295  return C2_inv * (gtilde_1_pre/(sqrt(lambda(q2)) * V(q2)) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 - (h_0[1]/q2 + h_1[1] + h_2[1] * q2))).real();
296  }
297 
303  double getgtilde_1_im(double q2)
304  {
305  updateParameters();
306  return C2_inv * (gtilde_1_pre/(sqrt(lambda(q2)) * V(q2)) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 - (h_0[1]/q2 + h_1[1] + h_2[1] * q2))).imag();
307  }
308 
314  double getgtilde_2_re(double q2)
315  {
316  updateParameters();
317  return C2_inv * (gtilde_2_pre/A_1(q2) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2)).real();
318  }
319 
325  double getgtilde_2_im(double q2)
326  {
327  updateParameters();
328  return C2_inv * (gtilde_2_pre/A_1(q2) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2)).imag();
329  }
330 
336  double getgtilde_3_re(double q2)
337  {
338  updateParameters();
339  return C2_inv * (gtilde_3_pre/(lambda(q2) * A_2(q2)) * (sqrt(q2)*(h_0[0]/q2 +
340  h_1[0] + h_2[0] * q2)-(MM2mMV2 - q2)/(4.*MV) * (h_0[2]/q2 + h_1[2] +
341  h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2))).real();
342  }
343 
349  double getgtilde_3_im(double q2)
350  {
351  updateParameters();
352  return C2_inv * (gtilde_3_pre/(lambda(q2) * A_2(q2)) * (sqrt(q2)*(h_0[0]/q2 +
353  h_1[0] + h_2[0] * q2)-(MM2mMV2 - q2)/(4.*MV) * (h_0[2]/q2 + h_1[2] +
354  h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2))).imag();
355  }
356 
362  double geth_0_re(double q2)
363  {
364  return (sixteenM_PI2MM2 * (h_0[0]/q2 + h_1[0] + h_2[0] * q2)).real();
365  }
366 
372  double geth_0_im(double q2)
373  {
374  return (sixteenM_PI2MM2 * (h_0[0]/q2 + h_1[0] + h_2[0] * q2)).imag();
375  }
376 
382  double geth_p_re(double q2)
383  {
384  return (sixteenM_PI2MM2 * (h_0[1]/q2 + h_1[1] + h_2[1] * q2)).real();
385  }
386 
392  double geth_p_im(double q2)
393  {
394  return (sixteenM_PI2MM2 * (h_0[1]/q2 + h_1[1] + h_2[1] * q2)).imag();
395  }
396 
402  double geth_m_re(double q2)
403  {
404  return (sixteenM_PI2MM2 * (h_0[2]/q2 + h_1[2] + h_2[2] * q2)).real();
405  }
406 
412  double geth_m_im(double q2)
413  {
414  return (sixteenM_PI2MM2 * (h_0[2]/q2 + h_1[2] + h_2[2] * q2)).imag();
415  }
416 
421  double gethp0_hm0_abs()
422  {
423  updateParameters();
424  return (h_0[1]/h_0[2]).abs();
425  }
426 
431  double gethm0_h00_abs()
432  {
433  updateParameters();
434  return (h_0[2]/h_0[0]).abs();
435  }
436 
437 private:
438  const StandardModel& mySM;
439  StandardModel::lepton lep;
440  StandardModel::meson meson;
441  StandardModel::meson vectorM;
443  double GF;
444  double ale;
445  double Mlep;
446  double MM;
447  double MV;
448  double Mb;
449  double mu_b;
450  double mu_h;
451  double Mc;
452  double Ms;
453  double width;
454  double MW;
455  gslpp::complex lambda_t;
456  double b;
457  gslpp::complex h_0[3];
458  gslpp::complex h_1[3];
459  gslpp::complex h_2[3];
461  double t_p;
462  double t_m;
463  double t_0;
464  double z_0;
465  double MMpMV;
466  double MMpMV2;
467  double MMmMV;
468  double MMmMV2;
469  double MM2;
470  double MM4;
471  double MV2;
472  double MV4;
473  double MMMV;
474  double MM2mMV2;
475  double fourMV;
476  double onepMMoMV;
477  double MM_MMpMV;
478  double twoMM2;
479  double twoMV2;
480  double twoMM_mbpms;
481  double fourMM2;
482  double Mlep2;
483  double twoMlepMb;
484  double MboMW;
485  double MsoMb;
486  double M_PI2osix;
487  double twoMM;
488  double m4downcharge;
489  double threeGegen0;
490  double threeGegen1otwo;
491  double twoMc2;
492  double ninetysixM_PI3MM3;
493  double sixteenM_PI2;
494  double sixteenM_PI2MM2;
495  double twoMboMM;
496  gslpp::complex H_0_pre;
497  double mu_b2;
498  double Mc2;
499  double Mb2;
500  double fourMc2;
501  double fourMb2;
502  double logMc;
503  double logMb;
504  gslpp::complex H_0_WC;
505  gslpp::complex H_c_WC;
506  gslpp::complex H_b_WC;
507  double fournineth;
508  double half;
509  double twothird;
510  gslpp::complex ihalfMPI;
511  double twoMM3;
512  double gtilde_1_pre;
513  double gtilde_2_pre;
514  double gtilde_3_pre;
515  double C2_inv;
516  double S_L_pre;
517  double NN;
518  double sixMMoMb;
519  double CF;
520  double deltaT_0;
521  double deltaT_1par;
522  double deltaT_1perp;
523  double Ee;
525  gslpp::complex ubar;
526  gslpp::complex arg1;
527  gslpp::complex B01;
528  gslpp::complex B00;
529  gslpp::complex xp;
530  gslpp::complex xm;
531  gslpp::complex yp;
532  gslpp::complex ym;
533  gslpp::complex L1xp;
534  gslpp::complex L1xm;
535  gslpp::complex L1yp;
536  gslpp::complex L1ym;
537  gslpp::complex F87_0;
538  gslpp::complex F87_1;
539  gslpp::complex F87_2;
540  gslpp::complex F87_3;
541  gslpp::complex F29_0;
542  gslpp::complex F29_L1;
543  gslpp::complex F29_1;
544  gslpp::complex F29_2;
545  gslpp::complex F29_3;
546  gslpp::complex F29_L1_1;
547  gslpp::complex F29_L1_2;
548  gslpp::complex F29_L1_3;
549  gslpp::complex F27_0;
550  gslpp::complex F27_1;
551  gslpp::complex F27_2;
552  gslpp::complex F27_3;
553  gslpp::complex F27_L1_1;
554  gslpp::complex F27_L1_2;
555  gslpp::complex F27_L1_3;
556  double F89_0;
557  double F89_1;
558  double F89_2;
559  double F89_3;
561  double a_0A12_LCSR;
562  double a_0T2_LCSR;
564  double a_0V;
565  double a_1V;
566  double a_2V;
567  double MRV_2;
568  double a_0A0;
569  double a_1A0;
570  double a_2A0;
571  double MRA0_2;
572  double a_0A1;
573  double a_1A1;
574  double a_2A1;
575  double MRA1_2;
576  double a_0A12;
577  double a_1A12;
578  double a_2A12;
579  double MRA12_2;
580  double a_0T1;
581  double a_1T1;
582  double a_2T1;
583  double MRT1_2;
584  double a_0T2;
585  double a_1T2;
586  double a_2T2;
587  double MRT2_2;
588  double a_0T23;
589  double a_1T23;
590  double a_2T23;
591  double MRT23_2;
593  gslpp::vector<gslpp::complex> ** allcoeff;
594  gslpp::vector<gslpp::complex> ** allcoeffh;
595  gslpp::vector<gslpp::complex> ** allcoeffprime;
597  gslpp::complex C_1;
598  gslpp::complex C_1L_bar;
599  gslpp::complex C_1Lh_bar;
600  gslpp::complex C_2;
601  gslpp::complex C_2L_bar;
602  gslpp::complex C_2Lh_bar;
603  gslpp::complex C_3;
604  gslpp::complex C_4;
605  gslpp::complex C_5;
606  gslpp::complex C_6;
607  gslpp::complex C_7;
608  gslpp::complex C_8L;
609  gslpp::complex C_8Lh;
610  gslpp::complex C_9;
611  gslpp::complex C_10;
612  gslpp::complex C_S;
613  gslpp::complex C_P;
615  gslpp::complex C_7p;
616  gslpp::complex C_9p;
617  gslpp::complex C_10p;
618  gslpp::complex C_Sp;
619  gslpp::complex C_Pp;
621  std::vector<double> ReDeltaC9_p_mumu;
622  std::vector<double> ImDeltaC9_p_mumu;
623  std::vector<double> ReDeltaC9_m_mumu;
624  std::vector<double> ImDeltaC9_m_mumu;
625  std::vector<double> ReDeltaC9_0_mumu;
626  std::vector<double> ImDeltaC9_0_mumu;
627  std::vector<double> ReDeltaC9_p_ee;
628  std::vector<double> ImDeltaC9_p_ee;
629  std::vector<double> ReDeltaC9_m_ee;
630  std::vector<double> ImDeltaC9_m_ee;
631  std::vector<double> ReDeltaC9_0_ee;
632  std::vector<double> ImDeltaC9_0_ee;
633  std::vector<double> myq2;
635  TFitResultPtr refres_p_mumu;
636  TFitResultPtr imfres_p_mumu;
637  TFitResultPtr refres_m_mumu;
638  TFitResultPtr imfres_m_mumu;
639  TFitResultPtr refres_0_mumu;
640  TFitResultPtr imfres_0_mumu;
642  TFitResultPtr refres_p_ee;
643  TFitResultPtr imfres_p_ee;
644  TFitResultPtr refres_m_ee;
645  TFitResultPtr imfres_m_ee;
646  TFitResultPtr refres_0_ee;
647  TFitResultPtr imfres_0_ee;
649  TGraph gr1;
650  TGraph gr2;
652  TF1 reffit;
653  TF1 imffit;
655  double tmpq2;
657  double avaSigma;
658  double avaDelta;
659  double avaDTPPR;
661  double errSigma;
662  double errDelta;
663  double errDTPPR;
665  gsl_function FS;
666  gsl_function FD;
667  gsl_function DTPPR;
669  gsl_integration_cquad_workspace * w_DTPPR;
670  gsl_integration_cquad_workspace * w_sigma;
671  gsl_integration_cquad_workspace * w_delta;
673  gsl_error_handler_t * old_handler;
675  std::map<std::pair<double, double>, gslpp::complex > cacheI1;
677  std::map<std::pair<double, double>, double > cacheSigma0;
678  std::map<std::pair<double, double>, double > cacheSigma1;
679  std::map<std::pair<double, double>, double > cacheSigma2;
680  std::map<std::pair<double, double>, double > cacheSigma3;
681  std::map<std::pair<double, double>, double > cacheSigma4;
682  std::map<std::pair<double, double>, double > cacheSigma5;
683  std::map<std::pair<double, double>, double > cacheSigma6;
684  std::map<std::pair<double, double>, double > cacheSigma7;
685  std::map<std::pair<double, double>, double > cacheSigma9;
686  std::map<std::pair<double, double>, double > cacheSigma10;
687  std::map<std::pair<double, double>, double > cacheSigma11;
689  std::map<std::pair<double, double>, double > cacheDelta0;
690  std::map<std::pair<double, double>, double > cacheDelta1;
691  std::map<std::pair<double, double>, double > cacheDelta2;
692  std::map<std::pair<double, double>, double > cacheDelta3;
693  std::map<std::pair<double, double>, double > cacheDelta7;
694  std::map<std::pair<double, double>, double > cacheDelta11;
696  unsigned int N_updated;
697  gslpp::vector<double> N_cache;
698  gslpp::complex Nc_cache;
700  unsigned int V_updated;
701  gslpp::vector<double> V_cache;
703  unsigned int A0_updated;
704  gslpp::vector<double> A0_cache;
706  unsigned int A1_updated;
707  gslpp::vector<double> A1_cache;
709  unsigned int T1_updated;
710  gslpp::vector<double> T1_cache;
712  unsigned int T2_updated;
713  gslpp::vector<double> T2_cache;
715  unsigned int k2_updated;
716  gslpp::vector<double> k2_cache;
718  unsigned int z_updated;
720  unsigned int lambda_updated;
722  unsigned int beta_updated;
723  double beta_cache;
725  unsigned int F_updated;
727  unsigned int VL1_updated;
728  unsigned int VL2_updated;
730  unsigned int TL1_updated;
731  unsigned int TL2_updated;
733  unsigned int VR1_updated;
734  unsigned int VR2_updated;
736  unsigned int TR1_updated;
737  unsigned int TR2_updated;
739  unsigned int VL0_updated;
740  gslpp::vector<double> VL0_cache;
742  unsigned int TL0_updated;
743  gslpp::vector<double> TL0_cache;
745  unsigned int VR0_updated;
747  unsigned int TR0_updated;
749  unsigned int Mb_Ms_updated;
751  unsigned int SL_updated;
752  gslpp::vector<double> SL_cache;
754  unsigned int SR_updated;
756  unsigned int C_1_updated;
757  gslpp::complex C_1_cache;
759  unsigned int C_2_updated;
760  gslpp::complex C_2_cache;
762  unsigned int C_3_updated;
763  gslpp::complex C_3_cache;
765  unsigned int C_4_updated;
766  gslpp::complex C_4_cache;
768  unsigned int C_5_updated;
769  gslpp::complex C_5_cache;
771  unsigned int C_6_updated;
772  gslpp::complex C_6_cache;
774  unsigned int C_7_updated;
775  gslpp::complex C_7_cache;
777  unsigned int C_9_updated;
778  gslpp::complex C_9_cache;
780  unsigned int C_10_updated;
781  gslpp::complex C_10_cache;
783  unsigned int C_7p_updated;
784  gslpp::complex C_7p_cache;
786  unsigned int C_9p_updated;
787  gslpp::complex C_9p_cache;
789  unsigned int C_10p_updated;
790  gslpp::complex C_10p_cache;
792  unsigned int C_S_updated;
793  gslpp::complex C_S_cache;
795  unsigned int C_P_updated;
796  gslpp::complex C_P_cache;
798  unsigned int C_Sp_updated;
799  gslpp::complex C_Sp_cache;
801  unsigned int C_Pp_updated;
802  gslpp::complex C_Pp_cache;
804  unsigned int C_2Lh_updated;
805  gslpp::complex C_2Lh_cache;
807  unsigned int C_8Lh_updated;
808  gslpp::complex C_8Lh_cache;
810  unsigned int Yupdated;
811  gslpp::vector<double> Ycache;
813  gslpp::complex h0Ccache[3];
814  gslpp::complex h1Ccache[3];
815  gslpp::complex h2Ccache[3];
817  unsigned int h0_updated;
818  unsigned int h1_updated;
819  unsigned int h2_updated;
821  unsigned int H_V0updated;
822  gslpp::vector<double> H_V0cache;
824  unsigned int H_V1updated;
825  gslpp::vector<double> H_V1cache;
827  unsigned int H_V2updated;
828  gslpp::vector<double> H_V2cache;
830  unsigned int H_A0updated;
831  unsigned int H_A1updated;
832  unsigned int H_A2updated;
834  unsigned int H_Supdated;
835  gslpp::vector<double> H_Scache;
837  unsigned int H_Pupdated;
838  gslpp::vector<double> H_Pcache;
840  unsigned int I0_updated;
841  unsigned int I1_updated;
842  unsigned int I2_updated;
843  unsigned int I3_updated;
844  unsigned int I4_updated;
845  unsigned int I5_updated;
846  unsigned int I6_updated;
847  unsigned int I7_updated;
848  unsigned int I8_updated;
849  unsigned int I9_updated;
850  unsigned int I10_updated;
851  unsigned int I11_updated;
853  std::map<std::pair<double, double>, unsigned int > I1Cached;
855  std::map<std::pair<double, double>, unsigned int > sigma0Cached;
856  std::map<std::pair<double, double>, unsigned int > sigma1Cached;
857  std::map<std::pair<double, double>, unsigned int > sigma2Cached;
858  std::map<std::pair<double, double>, unsigned int > sigma3Cached;
859  std::map<std::pair<double, double>, unsigned int > sigma4Cached;
860  std::map<std::pair<double, double>, unsigned int > sigma5Cached;
861  std::map<std::pair<double, double>, unsigned int > sigma6Cached;
862  std::map<std::pair<double, double>, unsigned int > sigma7Cached;
863  std::map<std::pair<double, double>, unsigned int > sigma9Cached;
864  std::map<std::pair<double, double>, unsigned int > sigma10Cached;
865  std::map<std::pair<double, double>, unsigned int > sigma11Cached;
867  std::map<std::pair<double, double>, unsigned int > delta0Cached;
868  std::map<std::pair<double, double>, unsigned int > delta1Cached;
869  std::map<std::pair<double, double>, unsigned int > delta2Cached;
870  std::map<std::pair<double, double>, unsigned int > delta3Cached;
871  std::map<std::pair<double, double>, unsigned int > delta7Cached;
872  std::map<std::pair<double, double>, unsigned int > delta11Cached;
874  std::map<double, unsigned int> deltaTparpCached;
875  std::map<double, unsigned int> deltaTparmCached;
876  std::map<double, unsigned int> deltaTperpCached;
878  std::map<double, gslpp::complex> cacheDeltaTparp;
879  std::map<double, gslpp::complex> cacheDeltaTparm;
880  std::map<double, gslpp::complex> cacheDeltaTperp;
882  unsigned int deltaTparpupdated;
883  unsigned int deltaTparmupdated;
884  unsigned int deltaTperpupdated;
886  unsigned int T_updated;
887  gslpp::vector<double> T_cache;
892  void updateParameters();
893 
897  void checkCache();
898 
904  double z(double q2);
905 
911  double V(double q2);
912 
913 
919  double A_0(double q2);
920 
921 
927  double A_1(double q2);
928 
934  double A_2(double q2);
935 
941  double T_1(double q2);
942 
943 
949  double T_2(double q2);
950 
956  double V_0t(double q2);
957 
963  double V_p(double q2);
964 
970  double V_m(double q2);
971 
977  double T_0t(double q2);
978 
984  double T_p(double q2);
985 
991  double T_m(double q2);
992 
998  double S_L(double q2);
999 
1005  gslpp::complex H_0(double q2);
1006 
1013  gslpp::complex H_c(double q2, double mu);
1014 
1020  gslpp::complex H_b(double q2);
1021 
1027  gslpp::complex Y(double q2);
1028 
1034  double k2 (double q2);
1035 
1041  double beta2 (double q2);
1042 
1048  double lambda(double q2);
1049 
1056  double F(double q2, double b_i);
1057 
1063  double I_1c(double q2);
1064 
1070  double I_1s(double q2);
1071 
1077  double I_2c(double q2);
1078 
1084  double I_2s(double q2);
1085 
1091  double I_3(double q2);
1092 
1098  double I_4(double q2);
1099 
1105  double I_5(double q2);
1106 
1112  double I_6c(double q2);
1113 
1119  double I_6s(double q2);
1120 
1126  double I_7(double q2);
1127 
1133  double I_8(double q2);
1134 
1140  double I_9(double q2);
1141 
1148  double Delta(int i, double q2);
1149 
1155  double getSigma1c(double q2)
1156  {
1157  return I_1c(q2);
1158  };
1159 
1165  double getSigma1s(double q2)
1166  {
1167  return I_1s(q2);
1168  };
1169 
1175  double getSigma2c(double q2)
1176  {
1177  return I_2c(q2);
1178  };
1179 
1185  double getSigma2s(double q2)
1186  {
1187  return I_2s(q2);
1188  };
1189 
1195  double getSigma3(double q2)
1196  {
1197  return I_3(q2);
1198  };
1199 
1205  double getSigma4(double q2)
1206  {
1207  return I_4(q2);
1208  };
1209 
1215  double getSigma5(double q2)
1216  {
1217  return I_5(q2);
1218  };
1219 
1225  double getSigma6s(double q2)
1226  {
1227  return I_6s(q2);
1228  };
1229 
1235  double getSigma6c(double q2)
1236  {
1237  return I_6c(q2);
1238  };
1239 
1245  double getSigma7(double q2)
1246  {
1247  return I_7(q2);
1248  };
1249 
1255  double getSigma8(double q2)
1256  {
1257  return I_8(q2);
1258  };
1259 
1265  double getSigma9(double q2)
1266  {
1267  return I_9(q2);
1268  };
1269 
1275  double getDelta0(double q2)
1276  {
1277  return Delta(0, q2);
1278  };
1279 
1285  double getDelta1(double q2)
1286  {
1287  return Delta(1, q2);
1288  };
1289 
1295  double getDelta2(double q2)
1296  {
1297  return Delta(2, q2);
1298  };
1299 
1305  double getDelta3(double q2)
1306  {
1307  return Delta(3, q2);
1308  };
1309 
1315  double getDelta7(double q2)
1316  {
1317  return Delta(7, q2);
1318  };
1319 
1325  double getDelta11(double q2)
1326  {
1327  return Delta(11, q2);
1328  };
1329 
1330 
1340  double FF_fit(double q2, double a_0, double a_1, double a_2, double MR2);
1341 
1348  gslpp::complex I1(double u, double q2);
1349 
1356  gslpp::complex Tperpplus(double u, double q2);
1357 
1364  gslpp::complex Tparplus(double u, double q2);
1365 
1372  gslpp::complex Tparminus(double u, double q2);
1373 
1379  double Integrand_ReTperpplus(double up);
1380 
1386  double Integrand_ImTperpplus(double up);
1387 
1393  double Integrand_ReTparplus(double up);
1394 
1400  double Integrand_ImTparplus(double up);
1401 
1407  double Integrand_ReTparminus(double up);
1408 
1414  double Integrand_ImTparminus(double up);
1415 
1421  double Integrand_ReTpar_pm(double up);
1422 
1428  double Integrand_ImTpar_pm(double up);
1429 
1435  gslpp::complex F19(double q2);
1436 
1442  gslpp::complex F27(double q2);
1443 
1449  gslpp::complex F29(double q2);
1450 
1456  gslpp::complex F87(double q2);
1457 
1463  double F89(double q2);
1464 
1470  gslpp::complex Cperp(double q2);
1471 
1477  gslpp::complex Cpar(double q2);
1478 
1484  gslpp::complex deltaTperp(double q2);
1485 
1491  gslpp::complex deltaTpar(double q2);
1492 
1498  gslpp::complex DeltaC9_0(double q2);
1499 
1505  gslpp::complex DeltaC9_p(double q2);
1506 
1512  gslpp::complex DeltaC9_m(double q2);
1513 
1520  double reDC9fit(double* x, double* p);
1521 
1528  double imDC9fit(double* x, double* p);
1529 
1533  void fit_DeltaC9_p_mumu();
1534 
1538  void fit_DeltaC9_m_mumu();
1539 
1543  void fit_DeltaC9_0_mumu();
1544 
1548  void fit_DeltaC9_p_ee();
1549 
1553  void fit_DeltaC9_m_ee();
1554 
1558  void fit_DeltaC9_0_ee();
1559 
1565  gslpp::complex fDeltaC9_p(double q2);
1566 
1572  gslpp::complex fDeltaC9_m(double q2);
1573 
1579  gslpp::complex fDeltaC9_0(double q2);
1580 
1581 };
1582 
1583 #endif /* MVLL_H */
1584 
MVll::z
double z(double q2)
The lattice parameter from arXiv:1310.3722v3.
Definition: MVll.cpp:754
MVll::a_0A12_LCSR
double a_0A12_LCSR
Definition: MVll.h:561
MVll::cacheSigma10
std::map< std::pair< double, double >, double > cacheSigma10
Definition: MVll.h:686
MVll::Integrand_ReTperpplus
double Integrand_ReTperpplus(double up)
The real part of the integral involving at fixed , according to .
Definition: MVll.cpp:884
MVll::MboMW
double MboMW
Definition: MVll.h:484
MVll::MM4
double MM4
Definition: MVll.h:470
MVll::MRA1_2
double MRA1_2
Definition: MVll.h:575
MVll::A1_updated
unsigned int A1_updated
Definition: MVll.h:706
MVll::h2Ccache
gslpp::complex h2Ccache[3]
Definition: MVll.h:815
MVll::Mlep
double Mlep
Definition: MVll.h:445
MVll::L1yp
gslpp::complex L1yp
Definition: MVll.h:535
MVll::I6_updated
unsigned int I6_updated
Definition: MVll.h:846
MVll::MRT23_2
double MRT23_2
Definition: MVll.h:591
MVll::A_2
double A_2(double q2)
The transverse form factor .
Definition: MVll.cpp:774
MVll::imfres_m_mumu
TFitResultPtr imfres_m_mumu
Definition: MVll.h:638
MVll::twoMc2
double twoMc2
Definition: MVll.h:491
MVll::lambda
double lambda(double q2)
The factor used in the angular coefficients .
Definition: MVll.cpp:1431
MVll::imfres_p_mumu
TFitResultPtr imfres_p_mumu
Definition: MVll.h:636
MVll::F87_3
gslpp::complex F87_3
Definition: MVll.h:540
MVll::I_3
double I_3(double q2)
The angular coefficient .
Definition: MVll.cpp:1463
MVll::imDC9fit
double imDC9fit(double *x, double *p)
The fit function for the imaginary part of the QCDF correction .
Definition: MVll.cpp:1079
MVll::xm
gslpp::complex xm
Definition: MVll.h:530
MVll::integrateSigma
double integrateSigma(int i, double q_min, double q_max)
The integral of from to .
Definition: MVll.cpp:1512
MVll::yp
gslpp::complex yp
Definition: MVll.h:531
MVll::gr1
TGraph gr1
Definition: MVll.h:649
MVll::F89_1
double F89_1
Definition: MVll.h:557
MVll::Mb2
double Mb2
Definition: MVll.h:499
MVll::C_2L_bar
gslpp::complex C_2L_bar
Definition: MVll.h:601
MVll::L1ym
gslpp::complex L1ym
Definition: MVll.h:536
MVll::I10_updated
unsigned int I10_updated
Definition: MVll.h:850
MVll::C_3_updated
unsigned int C_3_updated
Definition: MVll.h:762
MVll::F87_0
gslpp::complex F87_0
Definition: MVll.h:537
MVll::h1_updated
unsigned int h1_updated
Definition: MVll.h:818
MVll::old_handler
gsl_error_handler_t * old_handler
Definition: MVll.h:673
MVll::Ms
double Ms
Definition: MVll.h:452
MVll::getSigma1c
double getSigma1c(double q2)
The CP average .
Definition: MVll.h:1155
MVll::C_6
gslpp::complex C_6
Definition: MVll.h:606
MVll::cacheSigma0
std::map< std::pair< double, double >, double > cacheSigma0
Definition: MVll.h:677
MVll::F29_L1_1
gslpp::complex F29_L1_1
Definition: MVll.h:546
MVll::C_Pp
gslpp::complex C_Pp
Definition: MVll.h:619
MVll::deltaTparmCached
std::map< double, unsigned int > deltaTparmCached
Definition: MVll.h:875
MVll::errDelta
double errDelta
Definition: MVll.h:662
MVll::I1Cached
std::map< std::pair< double, double >, unsigned int > I1Cached
Definition: MVll.h:853
MVll::SL_cache
gslpp::vector< double > SL_cache
Definition: MVll.h:752
MVll::reffit
TF1 reffit
Definition: MVll.h:652
MVll::lep
StandardModel::lepton lep
Definition: MVll.h:439
MVll::h_1
gslpp::complex h_1[3]
Definition: MVll.h:458
MVll::getVp
double getVp(double q2)
The form factor .
Definition: MVll.h:171
MVll::getDelta1
double getDelta1(double q2)
The CP asymmetry .
Definition: MVll.h:1285
MVll::C_1_updated
unsigned int C_1_updated
Definition: MVll.h:756
MVll::twoMM
double twoMM
Definition: MVll.h:487
MVll::getSigma3
double getSigma3(double q2)
The CP average .
Definition: MVll.h:1195
MVll::cacheDeltaTperp
std::map< double, gslpp::complex > cacheDeltaTperp
Definition: MVll.h:880
MVll::deltaTparpCached
std::map< double, unsigned int > deltaTparpCached
Definition: MVll.h:874
MVll::fourMM2
double fourMM2
Definition: MVll.h:481
MVll::fit_DeltaC9_m_mumu
void fit_DeltaC9_m_mumu()
The fitting routine for the QCDF correction in the muon channel.
Definition: MVll.cpp:1148
MVll::C_3_cache
gslpp::complex C_3_cache
Definition: MVll.h:763
MVll::Ee
double Ee
Definition: MVll.h:523
MVll::C_9p_cache
gslpp::complex C_9p_cache
Definition: MVll.h:787
MVll::deltaT_0
double deltaT_0
Definition: MVll.h:520
MVll::Cperp
gslpp::complex Cperp(double q2)
The correction from .
Definition: MVll.cpp:992
MVll::deltaTpar
gslpp::complex deltaTpar(double q2)
The total correction from .
Definition: MVll.cpp:1024
MVll::V_p
double V_p(double q2)
The helicity form factor .
Definition: MVll.cpp:794
MVll::MMmMV
double MMmMV
Definition: MVll.h:467
MVll::getSigma7
double getSigma7(double q2)
The CP average .
Definition: MVll.h:1245
MVll::F27_L1_1
gslpp::complex F27_L1_1
Definition: MVll.h:553
MVll::F29_L1
gslpp::complex F29_L1
Definition: MVll.h:542
MVll::C_1
gslpp::complex C_1
Definition: MVll.h:597
MVll::VL2_updated
unsigned int VL2_updated
Definition: MVll.h:728
MVll::CF
double CF
Definition: MVll.h:519
MVll::MMmMV2
double MMmMV2
Definition: MVll.h:468
MVll::sigma6Cached
std::map< std::pair< double, double >, unsigned int > sigma6Cached
Definition: MVll.h:861
MVll::a_2T1
double a_2T1
Definition: MVll.h:582
MVll::H_V2cache
gslpp::vector< double > H_V2cache
Definition: MVll.h:828
MVll::MRV_2
double MRV_2
Definition: MVll.h:567
MVll::z_updated
unsigned int z_updated
Definition: MVll.h:718
MVll::FF_fit
double FF_fit(double q2, double a_0, double a_1, double a_2, double MR2)
The fit function from , .
Definition: MVll.cpp:749
MVll::C_5_cache
gslpp::complex C_5_cache
Definition: MVll.h:769
MVll::C_P_updated
unsigned int C_P_updated
Definition: MVll.h:795
MVll::H_A_m
gslpp::complex H_A_m(double q2)
The helicity amplitude .
Definition: MVll.cpp:1398
MVll::getT0
double getT0(double q2)
The form factor .
Definition: MVll.h:192
MVll::delta1Cached
std::map< std::pair< double, double >, unsigned int > delta1Cached
Definition: MVll.h:868
StandardModel.h
MVll::ReDeltaC9_0_mumu
std::vector< double > ReDeltaC9_0_mumu
Definition: MVll.h:625
MVll::I1_updated
unsigned int I1_updated
Definition: MVll.h:841
MVll::w_DTPPR
gsl_integration_cquad_workspace * w_DTPPR
Definition: MVll.h:669
MVll::ImDeltaC9_m_mumu
std::vector< double > ImDeltaC9_m_mumu
Definition: MVll.h:624
MVll::getDelta0
double getDelta0(double q2)
The CP asymmetry .
Definition: MVll.h:1275
MVll::B01
gslpp::complex B01
Definition: MVll.h:527
MVll::MV2
double MV2
Definition: MVll.h:471
MVll::integrateDelta
double integrateDelta(int i, double q_min, double q_max)
The integral of from to .
Definition: MVll.cpp:1675
MVll::I_5
double I_5(double q2)
The angular coefficient .
Definition: MVll.cpp:1473
MVll::FD
gsl_function FD
Definition: MVll.h:666
MVll::F_updated
unsigned int F_updated
Definition: MVll.h:725
MVll::F29_0
gslpp::complex F29_0
Definition: MVll.h:541
MVll::H_Pupdated
unsigned int H_Pupdated
Definition: MVll.h:837
MVll::ImDeltaC9_0_mumu
std::vector< double > ImDeltaC9_0_mumu
Definition: MVll.h:626
MVll::cacheSigma7
std::map< std::pair< double, double >, double > cacheSigma7
Definition: MVll.h:684
MVll::F27_0
gslpp::complex F27_0
Definition: MVll.h:549
MVll::C_9_cache
gslpp::complex C_9_cache
Definition: MVll.h:778
MVll::SR_updated
unsigned int SR_updated
Definition: MVll.h:754
MVll::C_P
gslpp::complex C_P
Definition: MVll.h:613
MVll::F89_0
double F89_0
Definition: MVll.h:556
MVll::H_A1updated
unsigned int H_A1updated
Definition: MVll.h:831
MVll::S_L_pre
double S_L_pre
Definition: MVll.h:516
MVll::A_1
double A_1(double q2)
The transverse form factor .
Definition: MVll.cpp:769
MVll::h2_updated
unsigned int h2_updated
Definition: MVll.h:819
MVll::ubar
gslpp::complex ubar
Definition: MVll.h:525
MVll::gethm0_h00_abs
double gethm0_h00_abs()
The absolute value of the ratio .
Definition: MVll.h:431
MVll::sigma9Cached
std::map< std::pair< double, double >, unsigned int > sigma9Cached
Definition: MVll.h:863
MVll::H_0_pre
gslpp::complex H_0_pre
Definition: MVll.h:496
MVll::cacheDelta0
std::map< std::pair< double, double >, double > cacheDelta0
Definition: MVll.h:689
MVll::fourMc2
double fourMc2
Definition: MVll.h:500
MVll::a_0A0
double a_0A0
Definition: MVll.h:568
MVll::C_1_cache
gslpp::complex C_1_cache
Definition: MVll.h:757
MVll::Mb_Ms_updated
unsigned int Mb_Ms_updated
Definition: MVll.h:749
MVll::MMpMV2
double MMpMV2
Definition: MVll.h:466
MVll::twoMboMM
double twoMboMM
Definition: MVll.h:495
MVll::I3_updated
unsigned int I3_updated
Definition: MVll.h:843
MVll::twoMM3
double twoMM3
Definition: MVll.h:511
MVll::Tparminus
gslpp::complex Tparminus(double u, double q2)
The function from .
Definition: MVll.cpp:877
MVll::refres_p_mumu
TFitResultPtr refres_p_mumu
Definition: MVll.h:635
gslFunctionAdapter
static double gslFunctionAdapter(double x, void *p)
Definition: MVll.h:29
MVll::MM2
double MM2
Definition: MVll.h:469
convertToGslFunction
gsl_function convertToGslFunction(const F &f)
Definition: MVll.h:38
MVll::F
double F(double q2, double b_i)
The factor used in the angular coefficients I_i.
Definition: MVll.cpp:1436
MVll::I11_updated
unsigned int I11_updated
Definition: MVll.h:851
MVll::DeltaC9_0
gslpp::complex DeltaC9_0(double q2)
The total QCDF correction computed integrating over .
Definition: MVll.cpp:1064
MVll::fDeltaC9_m
gslpp::complex fDeltaC9_m(double q2)
The total QCDF correction computed fitting over .
Definition: MVll.cpp:1301
MVll::cacheSigma6
std::map< std::pair< double, double >, double > cacheSigma6
Definition: MVll.h:683
MVll::getSigma2c
double getSigma2c(double q2)
The CP average .
Definition: MVll.h:1175
MVll::C_8Lh
gslpp::complex C_8Lh
Definition: MVll.h:609
QCD::meson
meson
An enum type for mesons.
Definition: QCD.h:713
MVll::getgtilde_2_im
double getgtilde_2_im(double q2)
The immaginary part of .
Definition: MVll.h:325
MVll::delta2Cached
std::map< std::pair< double, double >, unsigned int > delta2Cached
Definition: MVll.h:869
MVll::getgtilde_3_re
double getgtilde_3_re(double q2)
The real part of .
Definition: MVll.h:336
MVll::I2_updated
unsigned int I2_updated
Definition: MVll.h:842
MVll::F27_3
gslpp::complex F27_3
Definition: MVll.h:552
MVll::B00
gslpp::complex B00
Definition: MVll.h:528
MVll::h0_updated
unsigned int h0_updated
Definition: MVll.h:817
MVll::sigma2Cached
std::map< std::pair< double, double >, unsigned int > sigma2Cached
Definition: MVll.h:857
MVll::getgtilde_2_re
double getgtilde_2_re(double q2)
The real part of .
Definition: MVll.h:314
MVll::geth_m_re
double geth_m_re(double q2)
The real part of .
Definition: MVll.h:402
MVll::imfres_0_mumu
TFitResultPtr imfres_0_mumu
Definition: MVll.h:640
MVll::Mc
double Mc
Definition: MVll.h:451
MVll::V_updated
unsigned int V_updated
Definition: MVll.h:700
MVll::C_5
gslpp::complex C_5
Definition: MVll.h:605
MVll::H_Supdated
unsigned int H_Supdated
Definition: MVll.h:834
MVll::mu_b
double mu_b
Definition: MVll.h:449
MVll::cacheDelta7
std::map< std::pair< double, double >, double > cacheDelta7
Definition: MVll.h:693
MVll::F27_2
gslpp::complex F27_2
Definition: MVll.h:551
MVll::C_S_updated
unsigned int C_S_updated
Definition: MVll.h:792
MVll::ImDeltaC9_p_ee
std::vector< double > ImDeltaC9_p_ee
Definition: MVll.h:628
MVll::getDelta11
double getDelta11(double q2)
The CP asymmetry .
Definition: MVll.h:1325
MVll::getSigma9
double getSigma9(double q2)
The CP average .
Definition: MVll.h:1265
MVll::F29_2
gslpp::complex F29_2
Definition: MVll.h:544
MVll::a_2T23
double a_2T23
Definition: MVll.h:590
MVll::F87_1
gslpp::complex F87_1
Definition: MVll.h:538
MVll::F19
gslpp::complex F19(double q2)
The correction from .
Definition: MVll.cpp:946
MVll::C_Sp_updated
unsigned int C_Sp_updated
Definition: MVll.h:798
MVll::a_1T1
double a_1T1
Definition: MVll.h:581
MVll::MM2mMV2
double MM2mMV2
Definition: MVll.h:474
MVll::H_0
gslpp::complex H_0(double q2)
The function involved into .
Definition: MVll.cpp:1363
MVll::sigma7Cached
std::map< std::pair< double, double >, unsigned int > sigma7Cached
Definition: MVll.h:862
MVll::a_0V
double a_0V
Definition: MVll.h:564
MVll::twoMV2
double twoMV2
Definition: MVll.h:479
MVll::V_m
double V_m(double q2)
The helicity form factor .
Definition: MVll.cpp:799
MVll::refres_p_ee
TFitResultPtr refres_p_ee
Definition: MVll.h:642
MVll::xp
gslpp::complex xp
Definition: MVll.h:529
MVll::A0_cache
gslpp::vector< double > A0_cache
Definition: MVll.h:704
MVll::t_m
double t_m
Definition: MVll.h:462
MVll::myq2
std::vector< double > myq2
Definition: MVll.h:633
MVll::m4downcharge
double m4downcharge
Definition: MVll.h:488
StandardModel
A model class for the Standard Model.
Definition: StandardModel.h:443
MVll::deltaTperpCached
std::map< double, unsigned int > deltaTperpCached
Definition: MVll.h:876
MVll::cacheI1
std::map< std::pair< double, double >, gslpp::complex > cacheI1
Definition: MVll.h:675
MVll::H_V1cache
gslpp::vector< double > H_V1cache
Definition: MVll.h:825
MVll::H_b_WC
gslpp::complex H_b_WC
Definition: MVll.h:506
MVll::ReDeltaC9_p_ee
std::vector< double > ReDeltaC9_p_ee
Definition: MVll.h:627
MVll::C_2_updated
unsigned int C_2_updated
Definition: MVll.h:759
MVll::tmpq2
double tmpq2
Definition: MVll.h:655
MVll::F29
gslpp::complex F29(double q2)
The correction from .
Definition: MVll.cpp:969
MVll::T_p
double T_p(double q2)
The helicity form factor .
Definition: MVll.cpp:809
MVll::meson
StandardModel::meson meson
Definition: MVll.h:440
MVll::fit_DeltaC9_m_ee
void fit_DeltaC9_m_ee()
The fitting routine for the QCDF correction in the electron channel.
Definition: MVll.cpp:1184
MVll::H_b
gslpp::complex H_b(double q2)
The function involved into .
Definition: MVll.cpp:1352
MVll::C_8Lh_cache
gslpp::complex C_8Lh_cache
Definition: MVll.h:808
MVll::Integrand_ImTpar_pm
double Integrand_ImTpar_pm(double up)
The sum of Integrand_ImTparplus() and Integrand_ImTparminus().
Definition: MVll.cpp:938
MVll
A class for the decay.
Definition: MVll.h:87
MVll::MW
double MW
Definition: MVll.h:454
MVll::cacheSigma1
std::map< std::pair< double, double >, double > cacheSigma1
Definition: MVll.h:678
MVll::ImDeltaC9_p_mumu
std::vector< double > ImDeltaC9_p_mumu
Definition: MVll.h:622
MVll::Integrand_ImTperpplus
double Integrand_ImTperpplus(double up)
The imaginary part of the integral involving at fixed , according to .
Definition: MVll.cpp:892
MVll::Mc2
double Mc2
Definition: MVll.h:498
MVll::getVm
double getVm(double q2)
The form factor .
Definition: MVll.h:182
MVll::V_cache
gslpp::vector< double > V_cache
Definition: MVll.h:701
MVll::NN
double NN
Definition: MVll.h:517
MVll::Ycache
gslpp::vector< double > Ycache
Definition: MVll.h:811
MVll::MsoMb
double MsoMb
Definition: MVll.h:485
MVll::N_updated
unsigned int N_updated
Definition: MVll.h:696
MVll::C_Pp_cache
gslpp::complex C_Pp_cache
Definition: MVll.h:802
MVll::TL2_updated
unsigned int TL2_updated
Definition: MVll.h:731
MVll::C_Pp_updated
unsigned int C_Pp_updated
Definition: MVll.h:801
MVll::mu_b2
double mu_b2
Definition: MVll.h:497
MVll::getTm
double getTm(double q2)
The form factor .
Definition: MVll.h:214
MVll::ReDeltaC9_p_mumu
std::vector< double > ReDeltaC9_p_mumu
Definition: MVll.h:621
MVll::A0_updated
unsigned int A0_updated
Definition: MVll.h:703
MVll::deltaT_1perp
double deltaT_1perp
Definition: MVll.h:522
MVll::ImDeltaC9_0_ee
std::vector< double > ImDeltaC9_0_ee
Definition: MVll.h:632
MVll::beta_cache
double beta_cache
Definition: MVll.h:723
MVll::MM_MMpMV
double MM_MMpMV
Definition: MVll.h:477
MVll::F29_L1_2
gslpp::complex F29_L1_2
Definition: MVll.h:547
MVll::T2_updated
unsigned int T2_updated
Definition: MVll.h:712
MVll::geth_m_im
double geth_m_im(double q2)
The imaginary part of .
Definition: MVll.h:412
MVll::I_7
double I_7(double q2)
The angular coefficient .
Definition: MVll.cpp:1489
MVll::H_V0cache
gslpp::vector< double > H_V0cache
Definition: MVll.h:822
MVll::a_2T2
double a_2T2
Definition: MVll.h:586
MVll::k2
double k2(double q2)
The square of the 3-momentum of the recoiling meson in the M rest frame, .
Definition: MVll.cpp:1416
MVll::fourMV
double fourMV
Definition: MVll.h:475
MVll::SL_updated
unsigned int SL_updated
Definition: MVll.h:751
MVll::MV
double MV
Definition: MVll.h:447
MVll::H_c_WC
gslpp::complex H_c_WC
Definition: MVll.h:505
MVll::T_1
double T_1(double q2)
The transverse form factor .
Definition: MVll.cpp:779
MVll::allcoeffprime
gslpp::vector< gslpp::complex > ** allcoeffprime
Definition: MVll.h:595
MVll::beta2
double beta2(double q2)
The factor used in the angular coefficients .
Definition: MVll.cpp:1426
MVll::C_Sp_cache
gslpp::complex C_Sp_cache
Definition: MVll.h:799
MVll::cacheSigma5
std::map< std::pair< double, double >, double > cacheSigma5
Definition: MVll.h:682
MVll::TR2_updated
unsigned int TR2_updated
Definition: MVll.h:737
MVll::C_P_cache
gslpp::complex C_P_cache
Definition: MVll.h:796
MVll::DeltaC9_m
gslpp::complex DeltaC9_m(double q2)
The total QCDF correction computed integrating over .
Definition: MVll.cpp:1057
MVll::sigma10Cached
std::map< std::pair< double, double >, unsigned int > sigma10Cached
Definition: MVll.h:864
MVll::getgtilde_1_re
double getgtilde_1_re(double q2)
The real part of .
Definition: MVll.h:292
MVll::fournineth
double fournineth
Definition: MVll.h:507
MVll::MRA0_2
double MRA0_2
Definition: MVll.h:571
MVll::TR0_updated
unsigned int TR0_updated
Definition: MVll.h:747
MVll::sigma3Cached
std::map< std::pair< double, double >, unsigned int > sigma3Cached
Definition: MVll.h:858
MVll::a_1V
double a_1V
Definition: MVll.h:565
MVll::C_4_cache
gslpp::complex C_4_cache
Definition: MVll.h:766
MVll::getS
double getS(double q2)
The form factor .
Definition: MVll.h:225
MVll::a_0A12
double a_0A12
Definition: MVll.h:576
MVll::z_0
double z_0
Definition: MVll.h:464
MVll::ReDeltaC9_m_mumu
std::vector< double > ReDeltaC9_m_mumu
Definition: MVll.h:623
MVll::Yupdated
unsigned int Yupdated
Definition: MVll.h:810
MVll::delta0Cached
std::map< std::pair< double, double >, unsigned int > delta0Cached
Definition: MVll.h:867
MVll::S_L
double S_L(double q2)
The helicity form factor .
Definition: MVll.cpp:819
MVll::threeGegen0
double threeGegen0
Definition: MVll.h:489
MVll::F27_L1_3
gslpp::complex F27_L1_3
Definition: MVll.h:555
MVll::H_V2updated
unsigned int H_V2updated
Definition: MVll.h:827
MVll::C_3
gslpp::complex C_3
Definition: MVll.h:603
MVll::a_2V
double a_2V
Definition: MVll.h:566
MVll::gr2
TGraph gr2
Definition: MVll.h:650
MVll::getSigma6c
double getSigma6c(double q2)
The CP average .
Definition: MVll.h:1235
MVll::I0_updated
unsigned int I0_updated
Definition: MVll.h:840
MVll::A_0
double A_0(double q2)
The transverse form factor .
Definition: MVll.cpp:764
MVll::a_1A0
double a_1A0
Definition: MVll.h:569
MVll::I_1c
double I_1c(double q2)
The angular coefficient .
Definition: MVll.cpp:1441
MVll::VL1_updated
unsigned int VL1_updated
Definition: MVll.h:727
MVll::a_0T23
double a_0T23
Definition: MVll.h:588
MVll::cacheSigma11
std::map< std::pair< double, double >, double > cacheSigma11
Definition: MVll.h:687
MVll::cacheDeltaTparp
std::map< double, gslpp::complex > cacheDeltaTparp
Definition: MVll.h:878
MVll::Cpar
gslpp::complex Cpar(double q2)
The correction from .
Definition: MVll.cpp:998
MVll::TL0_cache
gslpp::vector< double > TL0_cache
Definition: MVll.h:743
MVll::b
double b
Definition: MVll.h:456
MVll::ImDeltaC9_m_ee
std::vector< double > ImDeltaC9_m_ee
Definition: MVll.h:630
MVll::twoMM_mbpms
double twoMM_mbpms
Definition: MVll.h:480
MVll::gtilde_1_pre
double gtilde_1_pre
Definition: MVll.h:512
MVll::H_c
gslpp::complex H_c(double q2, double mu)
The function involved into .
Definition: MVll.cpp:1341
MVll::H_A_0
gslpp::complex H_A_0(double q2)
The helicity amplitude .
Definition: MVll.cpp:1388
MVll::a_1T23
double a_1T23
Definition: MVll.h:589
MVll::onepMMoMV
double onepMMoMV
Definition: MVll.h:476
MVll::width
double width
Definition: MVll.h:453
MVll::A1_cache
gslpp::vector< double > A1_cache
Definition: MVll.h:707
MVll::fit_DeltaC9_0_ee
void fit_DeltaC9_0_ee()
The fitting routine for the QCDF correction in the electron channel.
Definition: MVll.cpp:1247
MVll::GF
double GF
Definition: MVll.h:443
MVll::deltaT_1par
double deltaT_1par
Definition: MVll.h:521
MVll::getTp
double getTp(double q2)
The form factor .
Definition: MVll.h:203
MVll::logMc
double logMc
Definition: MVll.h:502
MVll::sixMMoMb
double sixMMoMb
Definition: MVll.h:518
MVll::getV0
double getV0(double q2)
The form factor .
Definition: MVll.h:160
MVll::H_V_p
gslpp::complex H_V_p(double q2)
The helicity amplitude .
Definition: MVll.cpp:1378
MVll::C_S_cache
gslpp::complex C_S_cache
Definition: MVll.h:793
MVll::Integrand_ReTparplus
double Integrand_ReTparplus(double up)
The real part of the integral involving at fixed , according to .
Definition: MVll.cpp:900
MVll::beta_updated
unsigned int beta_updated
Definition: MVll.h:722
MVll::H_V_m
gslpp::complex H_V_m(double q2)
The helicity amplitude .
Definition: MVll.cpp:1383
MVll::allcoeffh
gslpp::vector< gslpp::complex > ** allcoeffh
Definition: MVll.h:594
MVll::C_4
gslpp::complex C_4
Definition: MVll.h:604
MVll::H_S
gslpp::complex H_S(double q2)
The helicity amplitude .
Definition: MVll.cpp:1403
MVll::ihalfMPI
gslpp::complex ihalfMPI
Definition: MVll.h:510
MVll::beta
double beta(double q2)
The factor used in the angular coefficients .
Definition: MVll.cpp:1421
MVll::C_2Lh_cache
gslpp::complex C_2Lh_cache
Definition: MVll.h:805
MVll::T_0t
double T_0t(double q2)
The helicity form factor .
Definition: MVll.cpp:804
MVll::C_4_updated
unsigned int C_4_updated
Definition: MVll.h:765
MVll::VL0_cache
gslpp::vector< double > VL0_cache
Definition: MVll.h:740
MVll::F87
gslpp::complex F87(double q2)
The correction from .
Definition: MVll.cpp:978
MVll::a_0T2_LCSR
double a_0T2_LCSR
Definition: MVll.h:562
MVll::I_2c
double I_2c(double q2)
The angular coefficient .
Definition: MVll.cpp:1453
MVll::allcoeff
gslpp::vector< gslpp::complex > ** allcoeff
Definition: MVll.h:593
MVll::C_7p_updated
unsigned int C_7p_updated
Definition: MVll.h:783
MVll::T1_cache
gslpp::vector< double > T1_cache
Definition: MVll.h:710
MVll::refres_m_mumu
TFitResultPtr refres_m_mumu
Definition: MVll.h:637
MVll::N_cache
gslpp::vector< double > N_cache
Definition: MVll.h:697
MVll::F27
gslpp::complex F27(double q2)
The correction from .
Definition: MVll.cpp:960
MVll::TL1_updated
unsigned int TL1_updated
Definition: MVll.h:730
MVll::geth_p_im
double geth_p_im(double q2)
The imaginary part of .
Definition: MVll.h:392
MVll::C_Sp
gslpp::complex C_Sp
Definition: MVll.h:618
MVll::F29_3
gslpp::complex F29_3
Definition: MVll.h:545
MVll::getDelta2
double getDelta2(double q2)
The CP asymmetry .
Definition: MVll.h:1295
MVll::H_A2updated
unsigned int H_A2updated
Definition: MVll.h:832
MVll::reDC9fit
double reDC9fit(double *x, double *p)
The fit function for the real part of the QCDF correction .
Definition: MVll.cpp:1074
MVll::C_2_cache
gslpp::complex C_2_cache
Definition: MVll.h:760
MVll::t_0
double t_0
Definition: MVll.h:463
MVll::C_9p_updated
unsigned int C_9p_updated
Definition: MVll.h:786
MVll::C_9p
gslpp::complex C_9p
Definition: MVll.h:616
MVll::a_1A12
double a_1A12
Definition: MVll.h:577
MVll::C_6_updated
unsigned int C_6_updated
Definition: MVll.h:771
MVll::MMMV
double MMMV
Definition: MVll.h:473
MVll::sixteenM_PI2MM2
double sixteenM_PI2MM2
Definition: MVll.h:494
MVll::sixteenM_PI2
double sixteenM_PI2
Definition: MVll.h:493
MVll::gethp0_hm0_abs
double gethp0_hm0_abs()
The absolute value of the ratio .
Definition: MVll.h:421
MVll::avaDelta
double avaDelta
Definition: MVll.h:658
MVll::I_2s
double I_2s(double q2)
The angular coefficient .
Definition: MVll.cpp:1458
MVll::I7_updated
unsigned int I7_updated
Definition: MVll.h:847
MVll::T_2
double T_2(double q2)
The transverse form factor .
Definition: MVll.cpp:784
MVll::geth_0_re
double geth_0_re(double q2)
The real part of .
Definition: MVll.h:362
MVll::cacheDelta1
std::map< std::pair< double, double >, double > cacheDelta1
Definition: MVll.h:690
MVll::gtilde_2_pre
double gtilde_2_pre
Definition: MVll.h:513
MVll::fourMb2
double fourMb2
Definition: MVll.h:501
MVll::C_7_cache
gslpp::complex C_7_cache
Definition: MVll.h:775
MVll::TR1_updated
unsigned int TR1_updated
Definition: MVll.h:736
MVll::F89
double F89(double q2)
The correction from .
Definition: MVll.cpp:985
MVll::I_4
double I_4(double q2)
The angular coefficient .
Definition: MVll.cpp:1468
MVll::cacheSigma4
std::map< std::pair< double, double >, double > cacheSigma4
Definition: MVll.h:681
MVll::F29_1
gslpp::complex F29_1
Definition: MVll.h:543
MVll::imfres_0_ee
TFitResultPtr imfres_0_ee
Definition: MVll.h:647
MVll::deltaTparpupdated
unsigned int deltaTparpupdated
Definition: MVll.h:882
MVll::Integrand_ImTparminus
double Integrand_ImTparminus(double up)
The imaginary part of the integral involving at fixed , according to .
Definition: MVll.cpp:927
MVll::L1xm
gslpp::complex L1xm
Definition: MVll.h:534
MVll::Integrand_ImTparplus
double Integrand_ImTparplus(double up)
The imaginary part of the integral involving at fixed , according to .
Definition: MVll.cpp:908
MVll::fit_DeltaC9_0_mumu
void fit_DeltaC9_0_mumu()
The fitting routine for the QCDF correction in the muon channel.
Definition: MVll.cpp:1211
MVll::fit_DeltaC9_p_mumu
void fit_DeltaC9_p_mumu()
The fitting routine for the QCDF correction in the muon channel.
Definition: MVll.cpp:1086
MVll::C_2
gslpp::complex C_2
Definition: MVll.h:600
MVll::cacheDelta3
std::map< std::pair< double, double >, double > cacheDelta3
Definition: MVll.h:692
MVll::h_2
gslpp::complex h_2[3]
Definition: MVll.h:459
MVll::I_6s
double I_6s(double q2)
The angular coefficient .
Definition: MVll.cpp:1479
MVll::fit_DeltaC9_p_ee
void fit_DeltaC9_p_ee()
The fitting routine for the QCDF correction in the electron channel.
Definition: MVll.cpp:1121
MVll::C_7
gslpp::complex C_7
Definition: MVll.h:607
MVll::w_delta
gsl_integration_cquad_workspace * w_delta
Definition: MVll.h:671
MVll::C_7p
gslpp::complex C_7p
Definition: MVll.h:615
MVll::sigma0Cached
std::map< std::pair< double, double >, unsigned int > sigma0Cached
Definition: MVll.h:855
MVll::refres_0_mumu
TFitResultPtr refres_0_mumu
Definition: MVll.h:639
MVll::cacheSigma3
std::map< std::pair< double, double >, double > cacheSigma3
Definition: MVll.h:680
MVll::getgtilde_1_im
double getgtilde_1_im(double q2)
The immaginary part of .
Definition: MVll.h:303
MVll::C_5_updated
unsigned int C_5_updated
Definition: MVll.h:768
MVll::avaDTPPR
double avaDTPPR
Definition: MVll.h:659
MVll::geth_p_re
double geth_p_re(double q2)
The real part of .
Definition: MVll.h:382
MVll::C_9_updated
unsigned int C_9_updated
Definition: MVll.h:777
MVll::T_updated
unsigned int T_updated
Definition: MVll.h:886
MVll::twoMlepMb
double twoMlepMb
Definition: MVll.h:483
MVll::getSigma2s
double getSigma2s(double q2)
The CP average .
Definition: MVll.h:1185
MVll::a_2A12
double a_2A12
Definition: MVll.h:578
MVll::F27_L1_2
gslpp::complex F27_L1_2
Definition: MVll.h:554
MVll::errDTPPR
double errDTPPR
Definition: MVll.h:663
MVll::H_V0updated
unsigned int H_V0updated
Definition: MVll.h:821
MVll::getSigma6s
double getSigma6s(double q2)
The CP average .
Definition: MVll.h:1225
MVll::C_7p_cache
gslpp::complex C_7p_cache
Definition: MVll.h:784
MVll::ReDeltaC9_m_ee
std::vector< double > ReDeltaC9_m_ee
Definition: MVll.h:629
MVll::C_10p_cache
gslpp::complex C_10p_cache
Definition: MVll.h:790
MVll::F89_3
double F89_3
Definition: MVll.h:559
MVll::V
double V(double q2)
The transverse form factor .
Definition: MVll.cpp:759
MVll::t_p
double t_p
Definition: MVll.h:461
MVll::C_10
gslpp::complex C_10
Definition: MVll.h:611
MVll::C2_inv
double C2_inv
Definition: MVll.h:515
MVll::sigma5Cached
std::map< std::pair< double, double >, unsigned int > sigma5Cached
Definition: MVll.h:860
MVll::h1Ccache
gslpp::complex h1Ccache[3]
Definition: MVll.h:814
MVll::a_0T1
double a_0T1
Definition: MVll.h:580
MVll::ym
gslpp::complex ym
Definition: MVll.h:532
MVll::~MVll
virtual ~MVll()
Destructor.
Definition: MVll.cpp:79
MVll::cacheSigma2
std::map< std::pair< double, double >, double > cacheSigma2
Definition: MVll.h:679
MVll::Y
gslpp::complex Y(double q2)
The function involved into .
Definition: MVll.cpp:1368
MVll::gtilde_3_pre
double gtilde_3_pre
Definition: MVll.h:514
MVll::lambda_t
gslpp::complex lambda_t
Definition: MVll.h:455
MVll::twothird
double twothird
Definition: MVll.h:509
MVll::F87_2
gslpp::complex F87_2
Definition: MVll.h:539
MVll::imffit
TF1 imffit
Definition: MVll.h:653
MVll::C_S
gslpp::complex C_S
Definition: MVll.h:612
MVll::I5_updated
unsigned int I5_updated
Definition: MVll.h:845
MVll::I9_updated
unsigned int I9_updated
Definition: MVll.h:849
MVll::H_0_WC
gslpp::complex H_0_WC
Definition: MVll.h:504
MVll::Integrand_ReTparminus
double Integrand_ReTparminus(double up)
The real part of the integral involving at fixed , according to .
Definition: MVll.cpp:916
MVll::cacheDelta11
std::map< std::pair< double, double >, double > cacheDelta11
Definition: MVll.h:694
MVll::H_V1updated
unsigned int H_V1updated
Definition: MVll.h:824
MVll::H_P
gslpp::complex H_P(double q2)
The helicity amplitude .
Definition: MVll.cpp:1408
MVll::checkCache
void checkCache()
The caching method for MVll.
Definition: MVll.cpp:381
MVll::cacheDelta2
std::map< std::pair< double, double >, double > cacheDelta2
Definition: MVll.h:691
MVll::a_1A1
double a_1A1
Definition: MVll.h:573
MVll::avaSigma
double avaSigma
Definition: MVll.h:657
MVll::C_8L
gslpp::complex C_8L
Definition: MVll.h:608
MVll::C_10p
gslpp::complex C_10p
Definition: MVll.h:617
MVll::V_0t
double V_0t(double q2)
The helicity form factor .
Definition: MVll.cpp:789
MVll::T_cache
gslpp::vector< double > T_cache
Definition: MVll.h:887
MVll::mu_h
double mu_h
Definition: MVll.h:450
MVll::H_V_0
gslpp::complex H_V_0(double q2)
The helicity amplitude .
Definition: MVll.cpp:1373
MVll::C_6_cache
gslpp::complex C_6_cache
Definition: MVll.h:772
MVll::I8_updated
unsigned int I8_updated
Definition: MVll.h:848
MVll::I_6c
double I_6c(double q2)
The angular coefficient .
Definition: MVll.cpp:1484
MVll::fDeltaC9_0
gslpp::complex fDeltaC9_0(double q2)
The total QCDF correction computed fitting over .
Definition: MVll.cpp:1320
MVll::MMpMV
double MMpMV
Definition: MVll.h:465
MVll::a_0A1
double a_0A1
Definition: MVll.h:572
MVll::logMb
double logMb
Definition: MVll.h:503
MVll::ReDeltaC9_0_ee
std::vector< double > ReDeltaC9_0_ee
Definition: MVll.h:631
MVll::refres_0_ee
TFitResultPtr refres_0_ee
Definition: MVll.h:646
MVll::getSigma5
double getSigma5(double q2)
The CP average .
Definition: MVll.h:1215
MVll::getDelta3
double getDelta3(double q2)
The CP asymmetry .
Definition: MVll.h:1305
MVll::cacheDeltaTparm
std::map< double, gslpp::complex > cacheDeltaTparm
Definition: MVll.h:879
MVll::MRT2_2
double MRT2_2
Definition: MVll.h:587
MVll::I_8
double I_8(double q2)
The angular coefficient .
Definition: MVll.cpp:1495
MVll::VR1_updated
unsigned int VR1_updated
Definition: MVll.h:733
MVll::getwidth
double getwidth()
The width of the meson M.
Definition: MVll.h:134
MVll::delta11Cached
std::map< std::pair< double, double >, unsigned int > delta11Cached
Definition: MVll.h:872
MVll::sigma4Cached
std::map< std::pair< double, double >, unsigned int > sigma4Cached
Definition: MVll.h:859
gslpp::vector< double >
A class for constructing and defining operations on real vectors.
Definition: gslpp_vector_double.h:33
MVll::C_7_updated
unsigned int C_7_updated
Definition: MVll.h:774
MVll::vectorM
StandardModel::meson vectorM
Definition: MVll.h:441
MVll::twoMM2
double twoMM2
Definition: MVll.h:478
MVll::imfres_m_ee
TFitResultPtr imfres_m_ee
Definition: MVll.h:645
MVll::MV4
double MV4
Definition: MVll.h:472
MVll::MVll
MVll(const StandardModel &SM_i, StandardModel::meson meson_i, StandardModel::meson vector_i, StandardModel::lepton lep_i)
Constructor.
Definition: MVll.cpp:16
MVll::deltaTperp
gslpp::complex deltaTperp(double q2)
The total correction from .
Definition: MVll.cpp:1004
MVll::F89_2
double F89_2
Definition: MVll.h:558
MVll::a_1T2
double a_1T2
Definition: MVll.h:585
MVll::lambda_updated
unsigned int lambda_updated
Definition: MVll.h:720
MVll::C_2Lh_bar
gslpp::complex C_2Lh_bar
Definition: MVll.h:602
MVll::getSigma8
double getSigma8(double q2)
The CP average .
Definition: MVll.h:1255
MVll::refres_m_ee
TFitResultPtr refres_m_ee
Definition: MVll.h:644
MVll::VR2_updated
unsigned int VR2_updated
Definition: MVll.h:734
MVll::H_Pcache
gslpp::vector< double > H_Pcache
Definition: MVll.h:838
MVll::h0Ccache
gslpp::complex h0Ccache[3]
Definition: MVll.h:813
MVll::T_m
double T_m(double q2)
The helicity form factor .
Definition: MVll.cpp:814
MVll::C_1Lh_bar
gslpp::complex C_1Lh_bar
Definition: MVll.h:599
MVll::getMlep
double getMlep()
The mass of the lepton l.
Definition: MVll.h:143
MVll::H_A0updated
unsigned int H_A0updated
Definition: MVll.h:830
MVll::getgtilde_3_im
double getgtilde_3_im(double q2)
The imaginary part of .
Definition: MVll.h:349
MVll::C_10_updated
unsigned int C_10_updated
Definition: MVll.h:780
MVll::updateParameters
void updateParameters()
The update parameter method for MVll.
Definition: MVll.cpp:83
MVll::F27_1
gslpp::complex F27_1
Definition: MVll.h:550
MVll::delta3Cached
std::map< std::pair< double, double >, unsigned int > delta3Cached
Definition: MVll.h:870
MVll::Delta
double Delta(int i, double q2)
The CP asymmetry .
Definition: MVll.cpp:1505
MVll::VL0_updated
unsigned int VL0_updated
Definition: MVll.h:739
MVll::deltaTperpupdated
unsigned int deltaTperpupdated
Definition: MVll.h:884
MVll::TL0_updated
unsigned int TL0_updated
Definition: MVll.h:742
MVll::I1
gslpp::complex I1(double u, double q2)
The function from .
Definition: MVll.cpp:828
MVll::H_Scache
gslpp::vector< double > H_Scache
Definition: MVll.h:835
MVll::C_2Lh_updated
unsigned int C_2Lh_updated
Definition: MVll.h:804
h_0
A class for the correction in .
Definition: MVllObservables.h:1470
MVll::Mb
double Mb
Definition: MVll.h:448
MVll::deltaTparmupdated
unsigned int deltaTparmupdated
Definition: MVll.h:883
MVll::ninetysixM_PI3MM3
double ninetysixM_PI3MM3
Definition: MVll.h:492
MVll::k2_updated
unsigned int k2_updated
Definition: MVll.h:715
MVll::MRT1_2
double MRT1_2
Definition: MVll.h:583
MVll::Tparplus
gslpp::complex Tparplus(double u, double q2)
The function from .
Definition: MVll.cpp:864
MVll::Nc_cache
gslpp::complex Nc_cache
Definition: MVll.h:698
MVll::getDelta7
double getDelta7(double q2)
The CP asymmetry .
Definition: MVll.h:1315
MVll::fDeltaC9_p
gslpp::complex fDeltaC9_p(double q2)
The total QCDF correction computed fitting over .
Definition: MVll.cpp:1283
MVll::Mlep2
double Mlep2
Definition: MVll.h:482
gslpp::complex
A class for defining operations on and functions of complex numbers.
Definition: gslpp_complex.h:35
MVll::C_9
gslpp::complex C_9
Definition: MVll.h:610
MVll::C_8Lh_updated
unsigned int C_8Lh_updated
Definition: MVll.h:807
MVll::geth_0_im
double geth_0_im(double q2)
The imaginary part of .
Definition: MVll.h:372
MVll::w_sigma
gsl_integration_cquad_workspace * w_sigma
Definition: MVll.h:670
MVll::k2_cache
gslpp::vector< double > k2_cache
Definition: MVll.h:716
MVll::mySM
const StandardModel & mySM
Definition: MVll.h:438
MVll::cacheSigma9
std::map< std::pair< double, double >, double > cacheSigma9
Definition: MVll.h:685
ThObservable.h
MVll::a_0T2
double a_0T2
Definition: MVll.h:584
MVll::MM
double MM
Definition: MVll.h:446
MVll::a_2A1
double a_2A1
Definition: MVll.h:574
MVll::L1xp
gslpp::complex L1xp
Definition: MVll.h:533
MVll::imfres_p_ee
TFitResultPtr imfres_p_ee
Definition: MVll.h:643
MVll::getSigma
double getSigma(int i, double q_2)
The value of from to .
Definition: MVll.cpp:1630
MVll::H_A_p
gslpp::complex H_A_p(double q2)
The helicity amplitude .
Definition: MVll.cpp:1393
MVll::M_PI2osix
double M_PI2osix
Definition: MVll.h:486
MVll::T2_cache
gslpp::vector< double > T2_cache
Definition: MVll.h:713
MVll::MRA12_2
double MRA12_2
Definition: MVll.h:579
MVll::Tperpplus
gslpp::complex Tperpplus(double u, double q2)
The function from .
Definition: MVll.cpp:850
MVll::F29_L1_3
gslpp::complex F29_L1_3
Definition: MVll.h:548
MVll::FS
gsl_function FS
Definition: MVll.h:665
MVll::I_1s
double I_1s(double q2)
The angular coefficient .
Definition: MVll.cpp:1447
StandardModel::lepton
lepton
An enum type for leptons.
Definition: StandardModel.h:449
MVll::Integrand_ReTpar_pm
double Integrand_ReTpar_pm(double up)
The sum of Integrand_ReTparplus() and Integrand_ReTparminus().
Definition: MVll.cpp:942
MVll::ale
double ale
Definition: MVll.h:444
MVll::getSigma1s
double getSigma1s(double q2)
The CP average .
Definition: MVll.h:1165
MVll::arg1
gslpp::complex arg1
Definition: MVll.h:526
MVll::threeGegen1otwo
double threeGegen1otwo
Definition: MVll.h:490
MVll::getSigma4
double getSigma4(double q2)
The CP average .
Definition: MVll.h:1205
MVll::C_10p_updated
unsigned int C_10p_updated
Definition: MVll.h:789
MVll::DTPPR
gsl_function DTPPR
Definition: MVll.h:667
MVll::errSigma
double errSigma
Definition: MVll.h:661
MVll::sigma1Cached
std::map< std::pair< double, double >, unsigned int > sigma1Cached
Definition: MVll.h:856
MVll::a_2A0
double a_2A0
Definition: MVll.h:570
MVll::delta7Cached
std::map< std::pair< double, double >, unsigned int > delta7Cached
Definition: MVll.h:871
MVll::C_1L_bar
gslpp::complex C_1L_bar
Definition: MVll.h:598
MVll::DeltaC9_p
gslpp::complex DeltaC9_p(double q2)
The total QCDF correction computed integrating over .
Definition: MVll.cpp:1051
MVll::I4_updated
unsigned int I4_updated
Definition: MVll.h:844
MVll::I_9
double I_9(double q2)
The angular coefficient .
Definition: MVll.cpp:1500
MVll::sigma11Cached
std::map< std::pair< double, double >, unsigned int > sigma11Cached
Definition: MVll.h:865
MVll::T1_updated
unsigned int T1_updated
Definition: MVll.h:709
MVll::VR0_updated
unsigned int VR0_updated
Definition: MVll.h:745
MVll::half
double half
Definition: MVll.h:508
gslpp::sqrt
complex sqrt(const complex &z)
Definition: gslpp_complex.cpp:385
MVll::C_10_cache
gslpp::complex C_10_cache
Definition: MVll.h:781