flavour
Public Member Functions | Public Attributes | List of all members
BGL2 Class Reference

A second implementation of the BGL model, for testing purposes. More...

Inheritance diagram for BGL2:
Inheritance graph
[legend]
Collaboration diagram for BGL2:
Collaboration graph
[legend]

Public Member Functions

 BGL2 (int genL=2, int genQ=2, int lup=0, int qup=0, int mssm=0)
 
 ~BGL2 ()
 
parameters generateparameters (int max=0) const
 
parameters getlist (const parameters &p) const
 
ex mesonmixing (ex mesonmass, const Fermion &f1, const Fermion &f2) const
 
double bsgammawidth (double tanb_, double McH_, double MR_, double MI_, int option=0)
 
double epsK (double tanb_, double McH_, double MR_, double MI_, int option=0)
 
- Public Member Functions inherited from Model
 Model ()
 
virtual ~Model ()
 
virtual int veto (const parameters &p, int max=0) const
 
double likelihood (const parameters &p, bool check=1, int max=0) const
 calculates the probability of getting all the experimental measures if the model describes the reality More...
 
double loglike (const parameters &p, bool check=1, int max=0) const
 

Public Attributes

const double planck
 
const possymbol GF
 
const possymbol MZ
 
const possymbol MW
 
const possymbol Mh
 
const constant Mpip
 
const constant Mpi0
 
const constant MBp
 
const constant MB0
 
const constant MBs0
 
const constant MKp
 
const constant MK0
 
const constant MDp
 
const constant MD0
 
const constant MDsp
 
const constant MDs0
 
const constant Fpi
 
const constant FB
 
const constant FBs
 
const constant FK
 
const constant FD
 
const constant FDs
 
ex cos2
 
ex g
 
ex alpha
 
const possymbol tanb
 
const possymbol cp
 
const possymbol McH
 
const possymbol MR
 
const possymbol MI
 
const possymbol rho
 
const realsymbol Tparam
 
const realsymbol Sparam
 
const realsymbol QCD1
 
const realsymbol QCD2
 
possymbol Mu [3]
 
possymbol Md [3]
 
vector< Boson > bosons
 
lst replacements
 
ex Btaunu
 
ex BR_Htotaunu
 
ex BR_toptoHq
 
ex BtotaunuR
 
ex BtoDtaunuR
 
ex BtoD2taunuR
 
const Mixes mixes
 
lst conjtoabs
 
realsymbol mu
 
int iBtaunu
 
int iBDtaunu
 
int iBD2taunu
 
vector< int > BGLtype
 
ROOT::Math::Interpolator inter1
 
ROOT::Math::Interpolator inter2
 
ROOT::Math::Interpolator Mu_ [3]
 
ROOT::Math::Interpolator Md_ [3]
 
double mmmax
 
double stepsize
 
calcuex * epsilonK
 
- Public Attributes inherited from Model
TRandom3 * r
 

Detailed Description

A second implementation of the BGL model, for testing purposes.

Definition at line 28 of file draw.cpp.

Constructor & Destructor Documentation

BGL2::BGL2 ( int  genL = 2,
int  genQ = 2,
int  lup = 0,
int  qup = 0,
int  mssm = 0 
)
inline

Definition at line 31 of file draw.cpp.

References BGLmodels::Boson::C, BGLmodels::cParticle, BGLmodels::fElectron, BGLmodels::fMuon, BGLmodels::fTau, BGLmodels::hLeft, BGLmodels::hRight, BGLmodels::iDown, BGLmodels::iUp, BGLmodels::M_MW, BGLmodels::M_MZ, BGLmodels::Boson::mass, BGLmodels::Boson::reset(), BGLmodels::Boson::s, BGLmodels::sScalar, BGLmodels::sVector, BGLmodels::tLepton, and BGLmodels::tQuark.

31  :
32  planck(6.58211928e-25),
33  GF("G_F"),
34  MZ("M_Z"),
35  MW("M_W"),
36  Mpip("Mpip",0.1396,"M_{\\pi^+}",domain::real),
37  Mpi0("Mpi0",0.1349766,"M_{\\pi^0}",domain::real),
38  MBp("MBp",5.279,"M_{B^+}",domain::real),
39  MB0("MB0",5.2795,"M_{B^0}",domain::real),
40  MBs0("MBs0",5.3663,"M_{B_s^0}",domain::real),
41  MKp("MKp",0.493677,"MKp",domain::real),
42  MK0("MK0",0.497614,"MK0",domain::real),
43  MDp("MDp",1.86957,"MDp",domain::real),
44  MD0("MD0",1.86480,"MD0",domain::real),
45  MDsp("MDsp",1.96845,"MDsp",domain::real),
46  MDs0("MDs0",0),
47  Fpi("Fpi",0.132,"Fpi",domain::real),
48  FB("FB",0.189,"FB",domain::real),
49  FBs("FBs",0.225,"FBs",domain::real),
50  FK("FK",0.159,"FK",domain::real),
51  FD("FD",0.208,"FD",domain::real),
52  FDs("FDs",0.248,"FDs",domain::real),
53  //alpha(7.297352e-3*4*M_PI),
54  cos2(pow(MW/MZ,2)),
55  g(sqrt(GF*8/sqrt(ex(2)))*MW),
56  //g(sqrt(4*Pi*alpha/(1-cos2))),
57  tanb("tg\\beta"),
58  cp("cp"),
59  McH("M_{H^+}"),
60  MR("M_{R}"),
61  MI("M_{I}"),
62  Tparam("T_param"),
63  Sparam("S_param"),
64  QCD1("QCD_1"),
65  QCD2("QCD_2"),
66  mixes(tanb,cp, genL,genQ, lup, qup, mssm),
67  mu("\\mu"),
68  BGLtype(4,0),
69  mmmax(1000),
70  stepsize(1e-2)
71  {
72  alpha=pow(g,2)*(1-cos2)/(4*Pi);
73  replacements.append(GF==1.166371e-5);
74  replacements.append(MZ==M_MZ);
75  replacements.append(MW==M_MW);
76 
77  mixes.appendtolst(replacements);
78 
79  replacements.append(Pi==M_PI);
80  replacements.append(sqrt(ex(2))==sqrt(2));
81  replacements.append(Pi==M_PI);
82  replacements.append(sqrt(ex(2))==sqrt(2));
83 
84  Boson boson;
85 
86  realsymbol q3("q3");
87  ex vq3=dirac_slash(q3,4);
88  varidx jmu(mu,4,1);
89 
90  for(uint i=0;i<2;i++)
91  for(uint j=0;j<3;j++)
92  for(uint k=0;k<3;k++){
93  conjtoabs.append(conjugate(mixes.V[i][j][k])==pow(abs(mixes.V[i][j][k]),2)/mixes.V[i][j][k]);
94  }
95 
96  //W+ boson
97  boson.mass=MW;
98  boson.s=sVector;
99 
100  for(uint t=tLepton;t<=tQuark;t++) boson.C[t][iUp][iDown][hLeft]=mixes.V[t]*Matrix(g/sqrt(ex(2)));
101  Boson wboson=boson;
102  bosons.push_back(boson);
103  boson.reset();
104 
105  //H+ boson
106  boson.mass=McH;
107  boson.s=sScalar;
108 
109  for(uint t=tLepton;t<=tQuark;t++)
110  for(uint i=iUp;i<=iDown;i++) boson.C[t][iUp][iDown][i]=mixes.VN[t][i]*Matrix(g/MW/sqrt(ex(2)));
111  Boson chiggs=boson;
112  bosons.push_back(boson);
113  boson.reset();
114 
115  for(int b=bosons.size()-1;b>=0;b--){
116  boson.mass=bosons[b].mass;
117  boson.s=bosons[b].s;
118  if(boson.s==sVector)
119  for(uint t=tLepton;t<=tQuark;t++)
120  for(uint i=iUp;i<=iDown;i++)
121  for(uint j=iUp;j<=iDown;j++)
122  for(uint h=hLeft;h<=hRight;h++){
123  boson.C[t][i][j][h]=bosons[b].C[t][j][i][h].conjugate();
124  }
125  else for(uint t=tLepton;t<=tQuark;t++)
126  for(uint i=iUp;i<=iDown;i++)
127  for(uint j=iUp;j<=iDown;j++)
128  for(uint h=hLeft;h<=hRight;h++){
129  boson.C[t][i][j][hLeft]=bosons[b].C[t][j][i][hRight].conjugate();
130  boson.C[t][i][j][hRight]=bosons[b].C[t][j][i][hLeft].conjugate();
131  }
132  bosons.push_back(boson);
133  boson.reset();
134  }
135 
136  //(R+iI)/sqrt(2) boson
137  boson.mass=MR;
138  boson.s=sScalar;
139 
140  for(uint t=tLepton;t<=tQuark;t++){
141  boson.C[t][iDown][iDown][hRight]=mixes.N[t][iDown]*Matrix(g/MW/ex(2));
142  boson.C[t][iUp][iUp][hLeft]=mixes.N[t][iUp].conjugate()*Matrix(g/MW/ex(2));
143  boson.C[t][iDown][iDown][hLeft]=mixes.N[t][iDown].conjugate()*Matrix(g/MW/ex(2));
144  boson.C[t][iUp][iUp][hRight]=mixes.N[t][iUp]*Matrix(g/MW/ex(2));
145  }
146  bosons.push_back(boson);
147  boson.reset();
148 
149  //(R+iI)/sqrt(2) boson
150  boson.mass=MI;
151  boson.s=sScalar;
152 
153  for(uint t=tLepton;t<=tQuark;t++){
154  boson.C[t][iDown][iDown][hRight]=mixes.N[t][iDown]*Matrix(I*g/MW/ex(2));
155  boson.C[t][iUp][iUp][hLeft]=mixes.N[t][iUp].conjugate()*Matrix(I*g/MW/ex(2));
156  boson.C[t][iDown][iDown][hLeft]=mixes.N[t][iDown].conjugate()*Matrix(-I*g/MW/ex(2));
157  boson.C[t][iUp][iUp][hRight]=mixes.N[t][iUp]*Matrix(-I*g/MW/ex(2));
158  }
159  bosons.push_back(boson);
160  boson.reset();
161 
162  Fermion electron(tLepton,iDown,fElectron);
163  Fermion electronR(tLepton,iDown,fElectron,cParticle,hRight);
164 
165  Fermion muon(tLepton,iDown,fMuon);
166  Fermion muonR(tLepton,iDown,fMuon,cParticle,hRight);
167 
168  Fermion tau(tLepton,iDown,fTau);
169  Fermion tauR(tLepton,iDown,fTau,cParticle,hRight);
170  Fermion neutrino(tLepton,iUp);
171  Fermion neutrinotau(tLepton,iUp,fTau);
172  Fermion neutrinomuon(tLepton,iUp,fMuon);
173  Fermion neutrinoe(tLepton,iUp,fElectron);
174 
175  Fermion up(tQuark,iUp,fElectron);
176  Fermion down(tQuark,iDown,fElectron);
177  Fermion bottom(tQuark,iDown,fTau);
178  Fermion strange(tQuark,iDown,fMuon);
179  Fermion charm(tQuark,iUp,fMuon);
180  Fermion top(tQuark,iUp,fTau);
181 
182  Meson Pi0d(down,down,Mpi0,Fpi);
183  Meson Pi0u(down,down,Mpi0,Fpi);
184  Meson Pip(up,down,Mpip,Fpi);
185  Meson Pim(down,up,Mpip,Fpi);
186 
187  Meson K0(down,strange,MK0,FK);
188  Meson Kp(up,strange,MKp,FK);
189 
190  Meson D0(charm,up,MD0,FD);
191  Meson Dp(charm,down,MDp,FD);
192  Meson Dsp(charm,strange,MDsp,FDs);
193 
194  Meson B0(down,bottom,MB0,FB);
195  Meson Bp(up,bottom,MBp,FB);
196  Meson Bs0(strange,bottom,MBs0,FBs);
197 
198  lst sb;
199  //sb.append(mixes.M[tQuark][iUp][0][0]==0);
200  sb.append(pow(abs(mixes.V[0][2][2]),2)==1-pow(abs(mixes.V[0][1][2]),2)-pow(abs(mixes.V[0][0][2]),2));
201  sb.append(pow(abs(mixes.V[0][2][1]),2)==1-pow(abs(mixes.V[0][1][1]),2)-pow(abs(mixes.V[0][0][1]),2));
202 
203  //cout<<pow(sqrt(2)/8*pow(g/MW,2),2)<<endl;
204  //cout<<pow(1.166,2)<<endl;
205  double fK=0.156;
206  ex KKbar=ex(std::pow(fK,2))*mesonmixing(MK0,strange,down);
207  ex eK=0.94*imag_part(KKbar)/3.5e-15/sqrt(2);
208  cout<<"KKbar "<<KKbar<<endl;
209  KKbar=expand(KKbar.subs(replacements).subs(lst(abs(wild()*pow(MR,-2))==abs(wild())*pow(MR,-2))).subs(lst(log(wild()*pow(MR,-2))==log(wild())-2*log(MR))));
210  KKbar=expand(KKbar.evalf());
211  ex aKKbar=sqrt(KKbar.real_part()*KKbar.real_part()+KKbar.imag_part()*KKbar.imag_part());
212 
213  eK=0.94*imag_part(KKbar)/3.5e-15/sqrt(2);
214  eK=eK.subs(replacements).real_part();
215  eK=collect_common_factors(expand(eK.evalf()));
216  cout<<"eK"<<eK<<endl;
217  //add("a_eK",abs(eK),new limitedobs(2.2e-3));
218  epsilonK=new calcuex(new limitedobs(2*0.011e-3),abs(eK));
219 
220 
221 }
BGL2::mu
realsymbol mu
Definition: draw.cpp:341
BGL2::GF
const possymbol GF
Definition: draw.cpp:322
BGLmodels::Mixes::V
vector< Matrix > V
Definition: Formulas.h:204
BGLmodels::Boson
Gauge boson.
Definition: BGL.h:21
BGL2::mmmax
double mmmax
Definition: draw.cpp:347
BGL2::MB0
const constant MB0
Definition: draw.cpp:323
BGLmodels::Meson
a meson properties
Definition: Formulas.h:47
BGL2::conjtoabs
lst conjtoabs
Definition: draw.cpp:340
BGL2::QCD2
const realsymbol QCD2
Definition: draw.cpp:327
BGL2::Sparam
const realsymbol Sparam
Definition: draw.cpp:327
BGL2::MBp
const constant MBp
Definition: draw.cpp:323
BGLmodels::Boson::reset
void reset()
Definition: BGL.h:63
BGL2::planck
const double planck
Definition: draw.cpp:321
BGL2::FBs
const constant FBs
Definition: draw.cpp:324
BGL2::MW
const possymbol MW
Definition: draw.cpp:322
BGL2::MKp
const constant MKp
Definition: draw.cpp:323
BGL2::MZ
const possymbol MZ
Definition: draw.cpp:322
BGLmodels::M_MZ
constexpr double M_MZ
Definition: Formulas.h:58
BGL2::mixes
const Mixes mixes
Definition: draw.cpp:339
BGLmodels::Mixes::N
multivector< Matrix, 2 > N
Definition: Formulas.h:206
BGL2::FK
const constant FK
Definition: draw.cpp:324
BGL2::McH
const possymbol McH
Definition: draw.cpp:326
BGL2::replacements
lst replacements
Definition: draw.cpp:331
BGL2::Fpi
const constant Fpi
Definition: draw.cpp:324
uint
unsigned int uint
Definition: script.cpp:4
BGL2::MDp
const constant MDp
Definition: draw.cpp:323
limitedobs
An experimental measure which is an upper limit on a parameter with a given Confidence Level...
Definition: model.h:52
BGLmodels::Mixes::VN
multivector< Matrix, 2 > VN
Definition: Formulas.h:207
BGL2::mesonmixing
ex mesonmixing(ex mesonmass, const Fermion &f1, const Fermion &f2) const
Definition: draw.cpp:272
BGLmodels::Boson::C
multivector< Matrix, 4 > C
Definition: BGL.h:68
calcuex
class to do the calculus of a constraint based on a GiNaC symbolic expression
Definition: model.h:282
BGL2::cp
const possymbol cp
Definition: draw.cpp:326
BGL2::Tparam
const realsymbol Tparam
Definition: draw.cpp:327
BGL2::BGLtype
vector< int > BGLtype
Definition: draw.cpp:344
BGL2::MI
const possymbol MI
Definition: draw.cpp:326
BGL2::g
ex g
Definition: draw.cpp:325
BGL2::MDs0
const constant MDs0
Definition: draw.cpp:323
BGL2::FB
const constant FB
Definition: draw.cpp:324
BGLmodels::Mixes::appendtolst
void appendtolst(lst &reps) const
Definition: Formulas.h:181
BGL2::FD
const constant FD
Definition: draw.cpp:324
BGL2::MDsp
const constant MDsp
Definition: draw.cpp:323
BGL2::cos2
ex cos2
Definition: draw.cpp:325
BGL2::alpha
ex alpha
Definition: draw.cpp:325
BGL2::MK0
const constant MK0
Definition: draw.cpp:323
BGL2::QCD1
const realsymbol QCD1
Definition: draw.cpp:327
BGL2::MBs0
const constant MBs0
Definition: draw.cpp:323
BGL2::Mpi0
const constant Mpi0
Definition: draw.cpp:323
BGL2::bosons
vector< Boson > bosons
Definition: draw.cpp:329
BGL2::FDs
const constant FDs
Definition: draw.cpp:324
BGLmodels::M_MW
constexpr double M_MW
Definition: Formulas.h:59
BGL2::epsilonK
calcuex * epsilonK
Definition: draw.cpp:349
BGL2::Mpip
const constant Mpip
Definition: draw.cpp:323
BGLmodels::Fermion
a fermion properties
Definition: Formulas.h:32
BGL2::MR
const possymbol MR
Definition: draw.cpp:326
BGL2::tanb
const possymbol tanb
Definition: draw.cpp:326
BGLmodels::Boson::s
BSpin s
Definition: BGL.h:66
BGL2::stepsize
double stepsize
Definition: draw.cpp:347
BGL2::MD0
const constant MD0
Definition: draw.cpp:323

Here is the call graph for this function:

BGL2::~BGL2 ( )
inline

Definition at line 223 of file draw.cpp.

223 {epsilonK->~calcuex();}
calcuex::~calcuex
~calcuex()
Definition: model.h:285
BGL2::epsilonK
calcuex * epsilonK
Definition: draw.cpp:349

Member Function Documentation

double BGL2::bsgammawidth ( double  tanb_,
double  McH_,
double  MR_,
double  MI_,
int  option = 0 
)
inline

Definition at line 298 of file draw.cpp.

References BGLmodels::calcubtosgamma2::width().

298  {
299  parameters p=generateparameters();
300  p[0].value=pow(10.0,tanb_);
301  p[1].value=McH_;
302  p[2].value=MR_;
303  p[3].value=MI_;
304  calcubtosgamma2 cal(mixes);
305 
306  return cal.width(p,option);
307 }
BGLmodels::calcubtosgamma2
calculus of the constraints coming from the b->s gamma decay
Definition: Formulas.h:279
BGL2::generateparameters
parameters generateparameters(int max=0) const
Definition: draw.cpp:225
BGL2::mixes
const Mixes mixes
Definition: draw.cpp:339
parameters
vector of parameters
Definition: model.h:177

Here is the call graph for this function:

double BGL2::epsK ( double  tanb_,
double  McH_,
double  MR_,
double  MI_,
int  option = 0 
)
inline

Definition at line 309 of file draw.cpp.

References parameters::p.

309  {
310  parameters p=generateparameters();
311  p[0].value=pow(10.0,tanb_);
312  p[1].value=McH_;
313  p[2].value=MR_;
314  p[3].value=MI_;
315  p.p=lst(tanb==p[0].value,McH==p[1].value,MR==p[2].value,MI==p[3].value);
316 
317 
318  return epsilonK->error(p);
319 }
BGL2::generateparameters
parameters generateparameters(int max=0) const
Definition: draw.cpp:225
calcuex::error
double error(const parameters &p) const
Definition: model.h:310
BGL2::McH
const possymbol McH
Definition: draw.cpp:326
parameters
vector of parameters
Definition: model.h:177
parameters::p
lst p
Definition: model.h:231
BGL2::MI
const possymbol MI
Definition: draw.cpp:326
BGL2::epsilonK
calcuex * epsilonK
Definition: draw.cpp:349
BGL2::MR
const possymbol MR
Definition: draw.cpp:326
BGL2::tanb
const possymbol tanb
Definition: draw.cpp:326
parameters BGL2::generateparameters ( int  max = 0) const
inlinevirtual

Implements Model.

Definition at line 225 of file draw.cpp.

225  {
226  parameters p;
227  //x=log_10(tanb)
228  p.push_back(freeparameter(-3,3,r,stepsize));
229  //y=log_10(McH)
230  if(max==1) p.push_back(freeparameter(10,10000,r,stepsize));
231  else p.push_back(freeparameter(10,mmmax,r,stepsize));
232  //log_10(massR)
233  p.push_back(freeparameter(-200,200,r,stepsize));
234  //log_10(massI)
235  p.push_back(freeparameter(-50,50,r,stepsize));
236 
237  return p;
238 }
BGL2::mmmax
double mmmax
Definition: draw.cpp:347
Model::r
TRandom3 * r
Definition: model.h:405
parameters
vector of parameters
Definition: model.h:177
freeparameter
A parameter which will be fitted in the simulation.
Definition: model.h:124
BGL2::stepsize
double stepsize
Definition: draw.cpp:347
parameters BGL2::getlist ( const parameters &  p) const
inlinevirtual

Implements Model.

Definition at line 241 of file draw.cpp.

References parameters::p, and parameters::values.

241  {
242  //cout<<aux<<endl;
243  //double c2=(1+sqrt(1-4*sqrt(ex_to<numeric>(mudecay.subs(lst(tanb==exp(p[0].value),McH==p[1].value))).to_double())))/2;
244 
245  double x=pow(10.0,p[0].value);
246  //double y=pow(10.0,p[1].value);
247  //double z=pow(10.0,p[2].value);
248  //double w=pow(10.0,p[3].value);
249 
250  double y=p[1].value;
251  double z=y+p[2].value;
252  double w=z+p[3].value;
253 
254  parameters pp(p);
255  pp[0].value=x;
256  pp[2].value+=pp[1].value;
257  pp[3].value+=pp[2].value;
258  pp.values=vector<double>();
259  for(uint i=0; i<4; i++) pp.values.push_back(pp[i].value);
260  lst &l=pp.p;
261  l=lst(tanb==x,McH==y,MR==z,MI==w);
262  l.append(QCD1==inter1.Eval(y));
263  l.append(QCD2==inter2.Eval(y));
264 
265  for(uint i=0;i<3;i++){
266  l.append(Mu[i]==Mu_[i].Eval(log(y)));
267  l.append(Md[i]==Md_[i].Eval(log(y)));
268  }
269  return pp;
270 }
BGL2::inter1
ROOT::Math::Interpolator inter1
Definition: draw.cpp:345
BGL2::QCD2
const realsymbol QCD2
Definition: draw.cpp:327
BGL2::Md
possymbol Md[3]
Definition: draw.cpp:328
BGL2::inter2
ROOT::Math::Interpolator inter2
Definition: draw.cpp:345
BGL2::McH
const possymbol McH
Definition: draw.cpp:326
uint
unsigned int uint
Definition: script.cpp:4
parameters
vector of parameters
Definition: model.h:177
BGL2::Md_
ROOT::Math::Interpolator Md_[3]
Definition: draw.cpp:346
BGL2::MI
const possymbol MI
Definition: draw.cpp:326
BGL2::QCD1
const realsymbol QCD1
Definition: draw.cpp:327
BGL2::Mu_
ROOT::Math::Interpolator Mu_[3]
Definition: draw.cpp:346
BGL2::MR
const possymbol MR
Definition: draw.cpp:326
BGL2::tanb
const possymbol tanb
Definition: draw.cpp:326
BGL2::Mu
possymbol Mu[3]
Definition: draw.cpp:328
ex BGL2::mesonmixing ( ex  mesonmass,
const Fermion &  f1,
const Fermion &  f2 
) const
inline

Definition at line 272 of file draw.cpp.

272  {
273 
274  ex ret=0;
275 
276  ex v1=0, v2=0;
277  ex mq1=mixes.mass(f1),mq2=mixes.mass(f2);
278  ex m2q1=mq1*mq1, m2q2=mq2*mq2;
279 
280  for(uint i=0;i<bosons.size();i++)
281  if(bosons[i].s==0){
282  ex a=(bosons[i].couplingdaggerR(f2,f1)-bosons[i].couplingdaggerL(f2,f1));
283  v1=v1+pow(a/bosons[i].mass,2);
284 
285  ex b=(bosons[i].couplingdaggerR(f2,f1)+bosons[i].couplingdaggerL(f2,f1));
286  v2=v2+pow(b/bosons[i].mass,2);
287  }
288 
289  ex fc=mesonmass/(mixes.massnum(f1)+mixes.massnum(f2));
290  fc=pow(fc,2);
291 
292  ret=2*(-v1*(1+11*fc)+v2*(1+fc))*mesonmass/96;
293 
294  return collect_common_factors(ret.subs(conjtoabs));
295  //return expand(ret.subs(lst(exp(-I*wild())==1/exp(I*wild()),sin(wild())==sqrt(1-pow(cos(wild()),2)))));
296 }
BGL2::conjtoabs
lst conjtoabs
Definition: draw.cpp:340
BGL2::mixes
const Mixes mixes
Definition: draw.cpp:339
uint
unsigned int uint
Definition: script.cpp:4
BGLmodels::Mixes::massnum
double massnum(const Fermion &f) const
Definition: Formulas.h:179
BGLmodels::Mixes::mass
ex mass(const Fermion &f) const
Definition: Formulas.h:178
BGL2::bosons
vector< Boson > bosons
Definition: draw.cpp:329

Member Data Documentation

ex BGL2::alpha

Definition at line 325 of file draw.cpp.

vector<int> BGL2::BGLtype

Definition at line 344 of file draw.cpp.

vector< Boson > BGL2::bosons

Definition at line 329 of file draw.cpp.

ex BGL2::BR_Htotaunu

Definition at line 333 of file draw.cpp.

ex BGL2::BR_toptoHq

Definition at line 334 of file draw.cpp.

ex BGL2::Btaunu

Definition at line 332 of file draw.cpp.

ex BGL2::BtoD2taunuR

Definition at line 337 of file draw.cpp.

ex BGL2::BtoDtaunuR

Definition at line 336 of file draw.cpp.

ex BGL2::BtotaunuR

Definition at line 335 of file draw.cpp.

lst BGL2::conjtoabs

Definition at line 340 of file draw.cpp.

ex BGL2::cos2

Definition at line 325 of file draw.cpp.

const possymbol BGL2::cp

Definition at line 326 of file draw.cpp.

calcuex* BGL2::epsilonK

Definition at line 349 of file draw.cpp.

const constant BGL2::FB

Definition at line 324 of file draw.cpp.

const constant BGL2::FBs

Definition at line 324 of file draw.cpp.

const constant BGL2::FD

Definition at line 324 of file draw.cpp.

const constant BGL2::FDs

Definition at line 324 of file draw.cpp.

const constant BGL2::FK

Definition at line 324 of file draw.cpp.

const constant BGL2::Fpi

Definition at line 324 of file draw.cpp.

ex BGL2::g

Definition at line 325 of file draw.cpp.

const possymbol BGL2::GF

Definition at line 322 of file draw.cpp.

int BGL2::iBD2taunu

Definition at line 343 of file draw.cpp.

int BGL2::iBDtaunu

Definition at line 343 of file draw.cpp.

int BGL2::iBtaunu

Definition at line 343 of file draw.cpp.

ROOT::Math::Interpolator BGL2::inter1

Definition at line 345 of file draw.cpp.

ROOT::Math::Interpolator BGL2::inter2

Definition at line 345 of file draw.cpp.

const constant BGL2::MB0

Definition at line 323 of file draw.cpp.

const constant BGL2::MBp

Definition at line 323 of file draw.cpp.

const constant BGL2::MBs0

Definition at line 323 of file draw.cpp.

const possymbol BGL2::McH

Definition at line 326 of file draw.cpp.

possymbol BGL2::Md[3]

Definition at line 328 of file draw.cpp.

const constant BGL2::MD0

Definition at line 323 of file draw.cpp.

ROOT::Math::Interpolator BGL2::Md_[3]

Definition at line 346 of file draw.cpp.

const constant BGL2::MDp

Definition at line 323 of file draw.cpp.

const constant BGL2::MDs0

Definition at line 323 of file draw.cpp.

const constant BGL2::MDsp

Definition at line 323 of file draw.cpp.

const possymbol BGL2::Mh

Definition at line 322 of file draw.cpp.

const possymbol BGL2::MI

Definition at line 326 of file draw.cpp.

const Mixes BGL2::mixes

Definition at line 339 of file draw.cpp.

const constant BGL2::MK0

Definition at line 323 of file draw.cpp.

const constant BGL2::MKp

Definition at line 323 of file draw.cpp.

double BGL2::mmmax

Definition at line 347 of file draw.cpp.

const constant BGL2::Mpi0

Definition at line 323 of file draw.cpp.

const constant BGL2::Mpip

Definition at line 323 of file draw.cpp.

const possymbol BGL2::MR

Definition at line 326 of file draw.cpp.

possymbol BGL2::Mu[3]

Definition at line 328 of file draw.cpp.

realsymbol BGL2::mu

Definition at line 341 of file draw.cpp.

ROOT::Math::Interpolator BGL2::Mu_[3]

Definition at line 346 of file draw.cpp.

const possymbol BGL2::MW

Definition at line 322 of file draw.cpp.

const possymbol BGL2::MZ

Definition at line 322 of file draw.cpp.

const double BGL2::planck

Definition at line 321 of file draw.cpp.

const realsymbol BGL2::QCD1

Definition at line 327 of file draw.cpp.

const realsymbol BGL2::QCD2

Definition at line 327 of file draw.cpp.

lst BGL2::replacements

Definition at line 331 of file draw.cpp.

const possymbol BGL2::rho

Definition at line 326 of file draw.cpp.

const realsymbol BGL2::Sparam

Definition at line 327 of file draw.cpp.

double BGL2::stepsize

Definition at line 347 of file draw.cpp.

const possymbol BGL2::tanb

Definition at line 326 of file draw.cpp.

const realsymbol BGL2::Tparam

Definition at line 327 of file draw.cpp.

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