flavour
Public Member Functions | Public Attributes | Static Public Attributes | List of all members
BGLmodels::calcubtosgamma2 Class Reference

calculus of the constraints coming from the b->s gamma decay More...

#include <Formulas.h>

Inheritance diagram for BGLmodels::calcubtosgamma2:
Inheritance graph
[legend]
Collaboration diagram for BGLmodels::calcubtosgamma2:
Collaboration graph
[legend]

Public Member Functions

 calcubtosgamma2 (const Mixes &mixes)
 
double operator() (const parameters &p) const
 
double width (const parameters &p, int option=0) const
 
double A0 (double x) const
 
double A1 (double x) const
 
double A2 (double x) const
 
double A3 (double x) const
 

Public Attributes

ROOT::Math::Interpolator inter1
 
ROOT::Math::Interpolator inter2
 
ROOT::Math::Interpolator inter3
 
ROOT::Math::Interpolator inter4
 
ROOT::Math::Interpolator Mu_ [3]
 
ROOT::Math::Interpolator Md_ [3]
 
const uint ii
 
FUNCP_CUBA fp
 
const gauss2obs g1
 
const gauss2obs g2
 
double ratio
 

Static Public Attributes

static constexpr double calN =2.567e-3
 
static constexpr double a =7.8221
 
static constexpr double aee =0.4384
 
static constexpr double aer =-1.6981
 
static constexpr double a77 =0.8161
 
static constexpr double a7r =4.8802
 
static constexpr double a7er =-0.7827
 
static constexpr double a88 =0.0197
 
static constexpr double a8r =0.5680
 
static constexpr double a8er =-0.0601
 
static constexpr double a87r =0.1923
 
static constexpr double a7i =0.3546
 
static constexpr double a8i =-0.0987
 
static constexpr double aei =2.4997
 
static constexpr double a87i =-0.0487
 
static constexpr double a7ei =-0.9067
 
static constexpr double a8ei =-0.0661
 

Detailed Description

calculus of the constraints coming from the b->s gamma decay

Definition at line 279 of file Formulas.h.

Constructor & Destructor Documentation

BGLmodels::calcubtosgamma2::calcubtosgamma2 ( const Mixes &  mixes)
inline

Definition at line 286 of file Formulas.h.

References BGLmodels::C7SM_Mt, BGLmodels::C8SM_Mt, BGLmodels::iDown, BGLmodels::iUp, BGLmodels::Mixes::N_, BGLmodels::Mixes::replacements, BGLmodels::Mixes::tanb, BGLmodels::tQuark, BGLmodels::Mixes::VN_, and BGLmodels::Vud().

286  :
287  ii(2),
288  g1(3.43e-4,sqrt(2)*0.23e-4),
289  g2(9.2e-6,4e-6),
290  ratio(0){
291  //cout<<"C7 "<<C7SM_Mt<<" "<<C7SM_MW<<" "<<C7SM(std::pow(261.8/M_MW,2))<<endl;
292  double res[2];
293  constexpr double C7SM_[2]={C7SM_Mt,C7SM_Mt};
294  constexpr double C8SM_[2]={C8SM_Mt,C8SM_Mt};
295  for(uint j=0; j<2; j++){
296  const uint i=2;
297  const CD epsilon=conj(Vud[0][j])*Vud[0][i]/conj(Vud[2][j])/Vud[2][i];
298  const double upsilon=norm(conj(Vud[2][j])*Vud[2][i]/Vud[1][i]);
299  const CD R7=(C7SM_Mt)/C7SM_MW;
300  const CD R8=(C8SM_Mt)/C8SM_MW;
301  const CD R7_=0;
302  const CD R8_=0;
303 
304  res[j]=a+aee*norm(epsilon)+aer*epsilon.real()+aei*epsilon.imag();
305  res[j]+=a77*(norm(R7)+norm(R7_))+a7r*R7.real()+a7i*R7.imag();
306  res[j]+=a88*(norm(R8)+norm(R8_))+a8r*R8.real()+a8i*R8.imag();
307  res[j]+=a87r*(R8*conj(R7)+R8_*conj(R7_)).real()+a7er*(R7*conj(epsilon)).real()+a8er*(R8*conj(epsilon)).real();
308  res[j]+=a87i*(R8*conj(R7)+R8_*conj(R7_)).imag()+a7ei*(R7*conj(epsilon)).imag()+a8er*(R8*conj(epsilon)).imag();
309  res[j]*=calN/100*upsilon;
310  }
311  //cout<<"Btosgamma "<<res[0]/9.2e-6<<" "<<res[1]/3.15e-4<<endl;
312  ifstream finter("interpolation.dat");
313 
314  if(!finter.is_open()){
315  cout<<"ERROR: interpolation.dat not found"<<endl;
316  exit(1);
317  }
318  vector<double> vinter0, vinter1, vinter2;
319  while(!finter.eof()){
320  double a=0,b=0,c=0;
321  finter>>a>>b>>c;
322  if(a!=0){
323  // cout<<a<<" "<<b<<" "<<c<<endl;
324  vinter0.push_back(a);
325  vinter1.push_back(b);
326  vinter2.push_back(c);
327  }
328  }
329 
330  inter1.SetData(vinter0,vinter1);
331  inter2.SetData(vinter0,vinter2);
332 
333  finter.close();
334 
335  ifstream finter2("masses.dat");
336 
337  if(!finter2.is_open()){
338  cout<<"ERROR: masses.dat not found"<<endl;
339  exit(1);
340  }
341  vector<vector<double> > m_(7);
342  while(!finter2.eof()){
343  for(uint i=0; i<7;i++) {
344  double a=0;
345  finter2>>a;
346 
347  if(a!=0) {
348  if(i==0) a=log(a);
349  else if(i<4) a*=1e-3;
350  m_[i].push_back(a);
351  // cout<<a<<" ";
352  }
353  } //cout<<endl;
354  }
355  for(uint i=0; i<3;i++) {
356  Md_[i].SetData(m_[0],m_[2*i+1]);
357  Mu_[i].SetData(m_[0],m_[2*i+2]);
358  }
359 // cout<<"Eval "<<Mu_[2].Eval(log(100.0))<<endl;
360 // cout<<"Eval "<<Md_[2].Eval(log(100.0))<<endl;
361 
362  finter2.close();
363 
364  ifstream finter3("interpolation2.dat");
365 
366  if(!finter3.is_open()){
367  cout<<"ERROR: interpolation2.dat not found"<<endl;
368  exit(1);
369  }
370  vector<double> vinter20, vinter21, vinter22;
371  while(!finter3.eof()){
372  double a=0,b=0,c=0;
373  finter3>>a>>b>>c;
374  if(a!=0){
375  // cout<<a<<" "<<b<<" "<<c<<endl;
376  vinter20.push_back(a);
377  vinter21.push_back(b);
378  vinter22.push_back(c);
379  }
380  }
381 
382  inter3.SetData(vinter20,vinter21);
383  inter4.SetData(vinter20,vinter22);
384 
385  finter3.close();
386 
387  vector<ex> vex(24);
388 
389  const uint i=ii;
390  for(uint j=0;j<2;j++)
391  for(uint k=0;k<3;k++){
392  vex[j*6+k*2+0]=mixes.VN_[tQuark][iUp][k][j].conjugate()*mixes.VN_[tQuark][iUp][k][i];
393  vex[j*6+k*2+1]=mixes.N_[tQuark][iDown][j][k]*mixes.N_[tQuark][iDown][i][k].conjugate();
394  vex[j+k*2+12]=-mixes.VN_[tQuark][iUp][k][j].conjugate()*mixes.VN_[tQuark][iDown][k][i];
395  vex[j+k*2+18]=mixes.VN_[tQuark][iDown][k][j].conjugate()*mixes.VN_[tQuark][iDown][k][i];
396  }
397  lst l;
398  for(uint k=0;k<vex.size();k++){
399  vex[k]=vex[k].subs(mixes.replacements).evalf();
400  l.append(vex[k].real_part());
401  l.append(vex[k].imag_part());
402  }
403  compile_ex(l, lst(mixes.tanb), fp);
404  }
BGLmodels::calcubtosgamma2::g2
const gauss2obs g2
Definition: Formulas.h:661
BGLmodels::C8SM_Mt
constexpr double C8SM_Mt
Definition: Formulas.h:273
BGLmodels::calcubtosgamma2::inter1
ROOT::Math::Interpolator inter1
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::a7ei
static constexpr double a7ei
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::aer
static constexpr double aer
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::a8i
static constexpr double a8i
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a7r
static constexpr double a7r
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::inter3
ROOT::Math::Interpolator inter3
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::g1
const gauss2obs g1
Definition: Formulas.h:661
BGLmodels::calcubtosgamma2::a
static constexpr double a
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::aee
static constexpr double aee
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::Md_
ROOT::Math::Interpolator Md_[3]
Definition: Formulas.h:657
BGLmodels::calcubtosgamma2::a7i
static constexpr double a7i
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a87r
static constexpr double a87r
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a8er
static constexpr double a8er
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::Mu_
ROOT::Math::Interpolator Mu_[3]
Definition: Formulas.h:657
uint
unsigned int uint
Definition: script.cpp:4
BGLmodels::calcubtosgamma2::inter4
ROOT::Math::Interpolator inter4
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::a88
static constexpr double a88
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::aei
static constexpr double aei
Definition: Formulas.h:284
BGLmodels::C7SM_Mt
constexpr double C7SM_Mt
Definition: Formulas.h:272
BGLmodels::calcubtosgamma2::calN
static constexpr double calN
Definition: Formulas.h:282
BGLmodels::C8SM_MW
constexpr double C8SM_MW
Definition: Formulas.h:273
BGLmodels::calcubtosgamma2::a8r
static constexpr double a8r
Definition: Formulas.h:283
BGLmodels::CD
std::complex< double > CD
Definition: Formulas.h:65
BGLmodels::Vud
const Matrixx Vud(13.04 *M_PI/180, 0.201 *M_PI/180, 2.38 *M_PI/180, 1.2)
BGLmodels::calcubtosgamma2::inter2
ROOT::Math::Interpolator inter2
Definition: Formulas.h:656
BGLmodels::C7SM_MW
constexpr double C7SM_MW
Definition: Formulas.h:272
BGLmodels::calcubtosgamma2::a87i
static constexpr double a87i
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a7er
static constexpr double a7er
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::a77
static constexpr double a77
Definition: Formulas.h:283

Here is the call graph for this function:

Member Function Documentation

double BGLmodels::calcubtosgamma2::A0 ( double  x) const
inline

Definition at line 640 of file Formulas.h.

640  {
641  return x*(2+3*x-6*x*x+ x*x*x+6*x*std::log(x))/(24*std::pow(1-x,4));
642  }
double BGLmodels::calcubtosgamma2::A1 ( double  x) const
inline

Definition at line 644 of file Formulas.h.

644  {
645  return x*(-3+4*x-x*x-2*std::log(x))/(4*std::pow(1-x,3));
646  }
double BGLmodels::calcubtosgamma2::A2 ( double  x) const
inline

Definition at line 648 of file Formulas.h.

648  {
649  return x/(6*std::pow(1-x,3))*((-7+5*x+8*x*x)/6.0+x*std::log(x)/(1-x)*(-2+3*x));
650  }
double BGLmodels::calcubtosgamma2::A3 ( double  x) const
inline

Definition at line 652 of file Formulas.h.

652  {
653  return (-3+8*x-5*x*x+(6*x-4)*std::log(x))*x/(6*std::pow(1-x,3));
654  }
double BGLmodels::calcubtosgamma2::operator() ( const parameters &  p) const
inlinevirtual
Parameters
hipothesisthe theoretical hypothesis
Returns
the logarithm of the probability of measuring what was measured, assuming that the hypothesis is true

Implements calcu.

Definition at line 406 of file Formulas.h.

References BGLmodels::C7SM_Mt, BGLmodels::C8SM_Mt, BGLmodels::mt_mt, and BGLmodels::Vud().

406  {
407  double tanb=p[0].value;
408  double y=p[1].value;
409  double z=p[2].value;
410  double w=p[3].value;
411  double McH=y, MR=z, MI=w;
412 
413  double y0=y;
414  if(y<mt_mt) y0=mt_mt;
415  double QCD1[2]={inter3.Eval(y0),inter1.Eval(y)};
416  double QCD2[2]={inter4.Eval(y0),inter2.Eval(y)};
417 
418  double Mu[3],Md[3];
419 
420  for(uint i=0;i<3;i++){
421  Mu[i]=Mu_[i].Eval(log(y));
422  Md[i]=Md_[i].Eval(log(z));
423  }
424  const uint i=ii;
425  CD CC7[2],DD7[2],CC8[2],DD8[2];
426  double res[2];
427  // constexpr double C7SM_[2]={C7SM_Mt,C7SM_Mb};
428  // constexpr double C8SM_[2]={C8SM_Mt,C8SM_Mb};
429 
430  std::array<double,48> ret;
431  const int n=1,m=48;
432  fp(&n,&(tanb),&m,&(ret[0]));
433  for(uint j=0;j<2;j++){
434  const double mbottom=Md[i];
435  const double mstrange=Md[j];
436  //ex mbottom=mixes.M[tQuark][iDown][i][i];
437  //ex mstrange=mixes.M[tQuark][iDown][j][j];
438 
439  CD C7,D7,C8,D8;
440  for(uint k=0;k<3;k++){
441  double mup=Mu[k];
442  double mdown=Md[k];
443  //ex mup=mixes.M[tQuark][iUp][k][k];
444  //ex mdown=mixes.M[tQuark][iDown][k][k];
445  //f1+=
446  double mmu=std::pow(mup/McH,2);
447  double mmdR=std::pow(mdown/MR,2);
448  double mmdI=std::pow(mdown/MI,2);
449 
450  double A0u=A0(mmu);
451  double A1u=A1(mmu);
452  double A2u=A2(mmu);
453  double A3u=A3(mmu);
454  double A0d=(A0(mmdR)+A0(mmdI));
455  double A1d=(A1(mmdR)-A1(mmdI));
456 
457  CD f1(ret[j*12+4*k+0],ret[j*12+4*k+1]);
458  C7+=f1*A2u;
459  C8+=-2.0*f1*A0u;
460 
461  CD f2=CD(ret[36+j*2+4*k+0],ret[36+j*2+4*k+1])*mstrange*mbottom/mup/mup;
462  //CD f2=f1*mstrange*mbottom/mup/mup;
463  D7+=f2*A2u;
464  D8+=-2.0*f2*A0u;
465 
466  CD f12(ret[24+j*2+4*k+0],ret[24+j*2+4*k+1]);
467  C7+=f12*A3u;
468  C8+=2.0*f12*A1u;
469 
470  CD f4(ret[j*12+4*k+2],ret[j*12+4*k+3]);
471  C7+=f4*A0d/3.0;
472  C8+=-f4*A0d;
473 
474  C7+=f4*A1d/3.0;
475  C8+=-f4*A1d;
476 
477  CD f6=f4*mstrange*mbottom/mdown/mdown;
478  D7+=f6*A0d/3.0;
479  D8+=-f6*A0d;
480  }
481  uint j0=j;
482  CC7[j]=(QCD1[j]*C7+QCD2[j]*C8)/2.0/conj(Vud[2][j])/Vud[2][i];
483  DD7[j]=(QCD1[j]*D7+QCD2[j]*D8)/2.0/conj(Vud[2][j])/Vud[2][i];
484  const double QCD3=(3*QCD2[j]/8+QCD1[j]);
485  CC8[j]=QCD3*C8/2.0/conj(Vud[2][j])/Vud[2][i];
486  DD8[j]=QCD3*D8/2.0/conj(Vud[2][j])/Vud[2][i];
487  const CD epsilon=conj(Vud[0][j])*Vud[0][i]/conj(Vud[2][j])/Vud[2][i];
488  const double upsilon=norm(conj(Vud[2][j])*Vud[2][i]/Vud[1][i]);
489  const CD R7=(C7SM_Mt+CC7[j])/C7SM_MW;
490  const CD R8=(C8SM_Mt+CD(0)*CC8[j])/C8SM_MW;
491  const CD R7_=(DD7[j])/C7SM_MW;
492  const CD R8_=CD(0)*(DD8[j])/C8SM_MW;
493 
494 
495  res[j]=a+aee*norm(epsilon)+aer*epsilon.real()+aei*epsilon.imag();
496  res[j]+=a77*(norm(R7)+norm(R7_))+a7r*R7.real()+a7i*R7.imag();
497  res[j]+=a88*(norm(R8)+norm(R8_))+a8r*R8.real()+a8i*R8.imag();
498  res[j]+=a87r*(R8*conj(R7)+R8_*conj(R7_)).real()+a7er*(R7*conj(epsilon)).real()+a8er*(R8*conj(epsilon)).real();
499  res[j]+=a87i*(R8*conj(R7)+R8_*conj(R7_)).imag()+a7ei*(R7*conj(epsilon)).imag()+a8er*(R8*conj(epsilon)).imag();
500  res[j]*=calN/100*upsilon;
501 
502  /*res[j]=a+aee*norm(epsilon)+aer*epsilon.real()+aei*epsilon.imag();
503  res[j]+=a77*(norm(R7)+norm(R7_))+a7r*1+a7i*0;
504  res[j]+=a88*(norm(R8)+norm(R8_))+a8r*R8.real()+a8i*R8.imag();
505  res[j]+=a87r*1+a7er*(conj(epsilon)).real()+a8er*(R8*conj(epsilon)).real();
506  res[j]+=a87i*0+a7ei*(conj(epsilon)).imag()+a8er*(R8*conj(epsilon)).imag();
507  res[j]*=calN/100*upsilon;
508  */
509  }
510  double r1=3.15e-4+0.00247*(norm(CC7[1])+norm(DD7[1])-0.706*CC7[1].real());
511 
512  //ratio=res[0]/9.2e-6;
513  //cout<<"RATIO "<<ratio<<endl;
514  return g1.loglikelihood(r1)+0*g2.loglikelihood(res[0]);
515  }
BGLmodels::calcubtosgamma2::g2
const gauss2obs g2
Definition: Formulas.h:661
BGLmodels::C8SM_Mt
constexpr double C8SM_Mt
Definition: Formulas.h:273
BGLmodels::calcubtosgamma2::inter1
ROOT::Math::Interpolator inter1
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::a7ei
static constexpr double a7ei
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::aer
static constexpr double aer
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::A1
double A1(double x) const
Definition: Formulas.h:644
BGLmodels::calcubtosgamma2::a8i
static constexpr double a8i
Definition: Formulas.h:284
BGLmodels::mt_mt
constexpr double mt_mt
Definition: Formulas.h:262
BGLmodels::calcubtosgamma2::A0
double A0(double x) const
Definition: Formulas.h:640
BGLmodels::calcubtosgamma2::a7r
static constexpr double a7r
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::inter3
ROOT::Math::Interpolator inter3
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::g1
const gauss2obs g1
Definition: Formulas.h:661
BGLmodels::calcubtosgamma2::a
static constexpr double a
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::aee
static constexpr double aee
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::Md_
ROOT::Math::Interpolator Md_[3]
Definition: Formulas.h:657
BGLmodels::calcubtosgamma2::a7i
static constexpr double a7i
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a87r
static constexpr double a87r
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a8er
static constexpr double a8er
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::Mu_
ROOT::Math::Interpolator Mu_[3]
Definition: Formulas.h:657
uint
unsigned int uint
Definition: script.cpp:4
BGLmodels::calcubtosgamma2::inter4
ROOT::Math::Interpolator inter4
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::a88
static constexpr double a88
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::aei
static constexpr double aei
Definition: Formulas.h:284
gauss2obs::loglikelihood
double loglikelihood(double hipothesis) const
Definition: model.h:109
BGLmodels::C7SM_Mt
constexpr double C7SM_Mt
Definition: Formulas.h:272
BGLmodels::calcubtosgamma2::calN
static constexpr double calN
Definition: Formulas.h:282
BGLmodels::C8SM_MW
constexpr double C8SM_MW
Definition: Formulas.h:273
BGLmodels::calcubtosgamma2::a8r
static constexpr double a8r
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::A3
double A3(double x) const
Definition: Formulas.h:652
BGLmodels::CD
std::complex< double > CD
Definition: Formulas.h:65
BGLmodels::Vud
const Matrixx Vud(13.04 *M_PI/180, 0.201 *M_PI/180, 2.38 *M_PI/180, 1.2)
BGLmodels::calcubtosgamma2::inter2
ROOT::Math::Interpolator inter2
Definition: Formulas.h:656
BGLmodels::C7SM_MW
constexpr double C7SM_MW
Definition: Formulas.h:272
BGLmodels::calcubtosgamma2::a87i
static constexpr double a87i
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a7er
static constexpr double a7er
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::a77
static constexpr double a77
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::A2
double A2(double x) const
Definition: Formulas.h:648

Here is the call graph for this function:

double BGLmodels::calcubtosgamma2::width ( const parameters &  p,
int  option = 0 
) const
inline

Definition at line 517 of file Formulas.h.

References BGLmodels::C7SM_Mt, BGLmodels::C8SM_Mt, BGLmodels::mt_mt, and BGLmodels::Vud().

Referenced by BGL2::bsgammawidth(), and BGLmodels::BGL::bsgammawidth().

517  {
518  double tanb=p[0].value;
519  double y=p[1].value;
520  double z=p[2].value;
521  double w=p[3].value;
522  double McH=y, MR=z, MI=w;
523 
524  double y0=y;
525  if(y<mt_mt) y0=mt_mt;
526  double QCD1[2]={inter3.Eval(y0),inter1.Eval(y)};
527  double QCD2[2]={inter4.Eval(y0),inter2.Eval(y)};
528 
529  double Mu[3],Md[3];
530 
531  for(uint i=0;i<3;i++){
532  Mu[i]=Mu_[i].Eval(log(y));
533  Md[i]=Md_[i].Eval(log(z));
534  }
535  const uint i=ii;
536  CD CC7[2],DD7[2],CC8[2],DD8[2];
537  double res[2];
538  // constexpr double C7SM_[2]={C7SM_Mt,C7SM_Mb};
539  // constexpr double C8SM_[2]={C8SM_Mt,C8SM_Mb};
540 
541  std::array<double,24> ret;
542  const int n=1,m=24;
543  fp(&n,&(tanb),&m,&(ret[0]));
544  for(uint j=0;j<2;j++){
545  const double mbottom=Md[i];
546  const double mstrange=Md[j];
547  //ex mbottom=mixes.M[tQuark][iDown][i][i];
548  //ex mstrange=mixes.M[tQuark][iDown][j][j];
549 
550  CD C7,D7,C8,D8;
551  for(uint k=0;k<3;k++){
552  double mup=Mu[k];
553  double mdown=Md[k];
554  //ex mup=mixes.M[tQuark][iUp][k][k];
555  //ex mdown=mixes.M[tQuark][iDown][k][k];
556  //f1+=
557  double mmu=std::pow(mup/McH,2);
558  double mmdR=std::pow(mdown/MR,2);
559  double mmdI=std::pow(mdown/MI,2);
560  double A0u=0,A1u=0, A2u=0, A3u=0, A0d=0, A1d=0;
561 
562  if(option==0 || option==1){
563  A0u=A0(mmu);
564  A1u=A1(mmu);
565  A2u=A2(mmu);
566  A3u=A3(mmu);
567  }
568  if(option==0 || option==2){
569  A0d=(A0(mmdR)+A0(mmdI));
570  A1d=(A1(mmdR)-A1(mmdI));
571  }
572  if(option==3){
573  A0d=(A0(mmdR));
574  A1d=(A1(mmdR));
575  }
576  if(option==4){
577  A0d=(A0(mmdI));
578  A1d=(-A1(mmdI));
579  }
580 
581  CD f1(ret[j*12+4*k+0],ret[j*12+4*k+1]);
582  C7+=f1*A2u;
583  C8+=-2.0*f1*A0u;
584 
585  CD f2=f1*mstrange*mbottom/mup/mup;
586  D7+=f2*A2u;
587  D8+=-2.0*f2*A0u;
588 
589  C7+=-f1*A3u;
590  C8+=-2.0*f1*A1u;
591 
592  CD f4(ret[j*12+4*k+2],ret[j*12+4*k+3]);
593  C7+=f4*A0d/3.0;
594  C8+=-f4*A0d;
595 
596  C7+=f4*A1d/3.0;
597  C8+=-f4*A1d;
598 
599  CD f6=f4*mstrange*mbottom/mdown/mdown;
600  D7+=f6*A0d/3.0;
601  D8+=-f6*A0d;
602 
603  }
604  uint j0=j;
605  CC7[j]=(QCD1[j]*C7+QCD2[j]*C8)/2.0/conj(Vud[2][j])/Vud[2][i];
606  DD7[j]=(QCD1[j]*D7+QCD2[j]*D8)/2.0/conj(Vud[2][j])/Vud[2][i];
607  const double QCD3=(3*QCD2[j]/8+QCD1[j]);
608  CC8[j]=QCD3*C8/2.0/conj(Vud[2][j])/Vud[2][i];
609  DD8[j]=QCD3*D8/2.0/conj(Vud[2][j])/Vud[2][i];
610  const CD epsilon=conj(Vud[0][j])*Vud[0][i]/conj(Vud[2][j])/Vud[2][i];
611  const double upsilon=norm(conj(Vud[2][j])*Vud[2][i]/Vud[1][i]);
612  const CD R7=(C7SM_Mt+CC7[j])/C7SM_MW;
613  const CD R8=(C8SM_Mt+CD(0)*CC8[j])/C8SM_MW;
614  const CD R7_=(DD7[j])/C7SM_MW;
615  const CD R8_=CD(0)*(DD8[j])/C8SM_MW;
616 
617 
618  res[j]=a+aee*norm(epsilon)+aer*epsilon.real()+aei*epsilon.imag();
619  res[j]+=a77*(norm(R7)+norm(R7_))+a7r*R7.real()+a7i*R7.imag();
620  res[j]+=a88*(norm(R8)+norm(R8_))+a8r*R8.real()+a8i*R8.imag();
621  res[j]+=a87r*(R8*conj(R7)+R8_*conj(R7_)).real()+a7er*(R7*conj(epsilon)).real()+a8er*(R8*conj(epsilon)).real();
622  res[j]+=a87i*(R8*conj(R7)+R8_*conj(R7_)).imag()+a7ei*(R7*conj(epsilon)).imag()+a8er*(R8*conj(epsilon)).imag();
623  res[j]*=calN/100*upsilon;
624 
625  /*res[j]=a+aee*norm(epsilon)+aer*epsilon.real()+aei*epsilon.imag();
626  res[j]+=a77*(norm(R7)+norm(R7_))+a7r*1+a7i*0;
627  res[j]+=a88*(norm(R8)+norm(R8_))+a8r*R8.real()+a8i*R8.imag();
628  res[j]+=a87r*1+a7er*(conj(epsilon)).real()+a8er*(R8*conj(epsilon)).real();
629  res[j]+=a87i*0+a7ei*(conj(epsilon)).imag()+a8er*(R8*conj(epsilon)).imag();
630  res[j]*=calN/100*upsilon;
631  */
632  }
633  double r1=3.15e-4+0.00247*(norm(CC7[1])+norm(DD7[1])-0.706*CC7[1].real());
634 
635  //ratio=res[0]/9.2e-6;
636  //cout<<"RATIO "<<ratio<<endl;
637  return g1.error(r1);
638  }
BGLmodels::C8SM_Mt
constexpr double C8SM_Mt
Definition: Formulas.h:273
BGLmodels::calcubtosgamma2::inter1
ROOT::Math::Interpolator inter1
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::a7ei
static constexpr double a7ei
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::aer
static constexpr double aer
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::A1
double A1(double x) const
Definition: Formulas.h:644
BGLmodels::calcubtosgamma2::a8i
static constexpr double a8i
Definition: Formulas.h:284
BGLmodels::mt_mt
constexpr double mt_mt
Definition: Formulas.h:262
BGLmodels::calcubtosgamma2::A0
double A0(double x) const
Definition: Formulas.h:640
BGLmodels::calcubtosgamma2::a7r
static constexpr double a7r
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::inter3
ROOT::Math::Interpolator inter3
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::g1
const gauss2obs g1
Definition: Formulas.h:661
BGLmodels::calcubtosgamma2::a
static constexpr double a
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::aee
static constexpr double aee
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::Md_
ROOT::Math::Interpolator Md_[3]
Definition: Formulas.h:657
BGLmodels::calcubtosgamma2::a7i
static constexpr double a7i
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a87r
static constexpr double a87r
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a8er
static constexpr double a8er
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::Mu_
ROOT::Math::Interpolator Mu_[3]
Definition: Formulas.h:657
uint
unsigned int uint
Definition: script.cpp:4
BGLmodels::calcubtosgamma2::inter4
ROOT::Math::Interpolator inter4
Definition: Formulas.h:656
BGLmodels::calcubtosgamma2::a88
static constexpr double a88
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::aei
static constexpr double aei
Definition: Formulas.h:284
BGLmodels::C7SM_Mt
constexpr double C7SM_Mt
Definition: Formulas.h:272
BGLmodels::calcubtosgamma2::calN
static constexpr double calN
Definition: Formulas.h:282
gauss2obs::error
double error(double hipothesis) const
Definition: model.h:113
BGLmodels::C8SM_MW
constexpr double C8SM_MW
Definition: Formulas.h:273
BGLmodels::calcubtosgamma2::a8r
static constexpr double a8r
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::A3
double A3(double x) const
Definition: Formulas.h:652
BGLmodels::CD
std::complex< double > CD
Definition: Formulas.h:65
BGLmodels::Vud
const Matrixx Vud(13.04 *M_PI/180, 0.201 *M_PI/180, 2.38 *M_PI/180, 1.2)
BGLmodels::calcubtosgamma2::inter2
ROOT::Math::Interpolator inter2
Definition: Formulas.h:656
BGLmodels::C7SM_MW
constexpr double C7SM_MW
Definition: Formulas.h:272
BGLmodels::calcubtosgamma2::a87i
static constexpr double a87i
Definition: Formulas.h:284
BGLmodels::calcubtosgamma2::a7er
static constexpr double a7er
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::a77
static constexpr double a77
Definition: Formulas.h:283
BGLmodels::calcubtosgamma2::A2
double A2(double x) const
Definition: Formulas.h:648

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

constexpr double BGLmodels::calcubtosgamma2::a =7.8221
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a77 =0.8161
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a7ei =-0.9067
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a7er =-0.7827
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a7i =0.3546
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a7r =4.8802
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a87i =-0.0487
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a87r =0.1923
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a88 =0.0197
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a8ei =-0.0661
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a8er =-0.0601
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a8i =-0.0987
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::a8r =0.5680
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::aee =0.4384
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::aei =2.4997
static

Definition at line 284 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::aer =-1.6981
static

Definition at line 283 of file Formulas.h.

constexpr double BGLmodels::calcubtosgamma2::calN =2.567e-3
static

Definition at line 282 of file Formulas.h.

FUNCP_CUBA BGLmodels::calcubtosgamma2::fp

Definition at line 660 of file Formulas.h.

const gauss2obs BGLmodels::calcubtosgamma2::g1

Definition at line 661 of file Formulas.h.

const gauss2obs BGLmodels::calcubtosgamma2::g2

Definition at line 661 of file Formulas.h.

const uint BGLmodels::calcubtosgamma2::ii

Definition at line 659 of file Formulas.h.

ROOT::Math::Interpolator BGLmodels::calcubtosgamma2::inter1

Definition at line 656 of file Formulas.h.

ROOT::Math::Interpolator BGLmodels::calcubtosgamma2::inter2

Definition at line 656 of file Formulas.h.

ROOT::Math::Interpolator BGLmodels::calcubtosgamma2::inter3

Definition at line 656 of file Formulas.h.

ROOT::Math::Interpolator BGLmodels::calcubtosgamma2::inter4

Definition at line 656 of file Formulas.h.

ROOT::Math::Interpolator BGLmodels::calcubtosgamma2::Md_[3]

Definition at line 657 of file Formulas.h.

ROOT::Math::Interpolator BGLmodels::calcubtosgamma2::Mu_[3]

Definition at line 657 of file Formulas.h.

double BGLmodels::calcubtosgamma2::ratio
mutable

Definition at line 662 of file Formulas.h.

The documentation for this class was generated from the following file: