a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
MVllObservables.cpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2014 HEPfit Collaboration
3 *
4 *
5 * For the licensing terms see doc/COPYING.
6 */
7
8#include "MVllObservables.h"
9#include "MVll.h"
10#include "StandardModel.h"
11//#include "gslpp.h"
12
13/*******************************************************************************
14 * Observables *
15 * ****************************************************************************/
16
17
18
19P_1::P_1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
20: ThObservable(SM_i)
21{
22 lep = lep_i;
23 meson = meson_i;
24 vectorM = vector_i;
25
26 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
27}
28
29double P_1::computeThValue()
30{
31 double q_min = getBinMin();
32 double q_max = getBinMax();
33
34 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(4,q_min,q_max)/(2.* SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max));
35}
36
37/*Returns experimental value, defined according to 1510.04239.*/
38P_2::P_2(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
39: ThObservable(SM_i)
40{
41 lep = lep_i;
42 meson = meson_i;
43 vectorM = vector_i;
44
45 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
46}
47
48double P_2::computeThValue()
49{
50 double q_min = getBinMin();
51 double q_max = getBinMax();
52
53 double P_2_theory = SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(7,q_min,q_max)/(8.*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max));
54
55 return -P_2_theory;
56}
57
58
59P_3::P_3(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
60: ThObservable(SM_i)
61{
62 lep = lep_i;
63 meson = meson_i;
64 vectorM = vector_i;
65
66 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
67}
68
69/*Returns experimental value, defined according to 1510.04239.*/
70double P_3::computeThValue()
71{
72 double q_min = getBinMin();
73 double q_max = getBinMax();
74
75 double P_3_theory = -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(11,q_min,q_max)/(4.*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max));
76
77 return -P_3_theory;
78}
79
80
81P_4Prime::P_4Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
82: ThObservable(SM_i)
83{
84 lep = lep_i;
85 meson = meson_i;
86 vectorM = vector_i;
87
88 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
89}
90
91/*Returns experimental value, defined according to 1510.04239.*/
92double P_4Prime::computeThValue()
93{
94 double q_min = getBinMin();
95 double q_max = getBinMax();
96
97 double P_4p_theory = SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(5,q_min,q_max)/sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(2,q_min,q_max)*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max));
98
99 return -1./2.*P_4p_theory;
100}
101
102
103P_5Prime::P_5Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
104: ThObservable(SM_i)
105{
106 lep = lep_i;
107 meson = meson_i;
108 vectorM = vector_i;
109
110 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
111}
112
113double P_5Prime::computeThValue()
114{
115 double q_min = getBinMin();
116 double q_max = getBinMax();
117
118 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(6,q_min,q_max)/(2.*sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(2,q_min,q_max)*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max)));
119}
120
121
122P_6Prime::P_6Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
123: ThObservable(SM_i)
124{
125 lep = lep_i;
126 meson = meson_i;
127 vectorM = vector_i;
128
129 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
130}
131
132double P_6Prime::computeThValue()
133{
134 double q_min = getBinMin();
135 double q_max = getBinMax();
136
137 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(9,q_min,q_max)/(2.*sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(2,q_min,q_max)*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max)));
138}
139
140
141P_8Prime::P_8Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
142: ThObservable(SM_i)
143{
144 lep = lep_i;
145 meson = meson_i;
146 vectorM = vector_i;
147
148 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
149}
150
151/*Returns experimental value, defined according to 1510.04239.*/
152double P_8Prime::computeThValue()
153{
154 double q_min = getBinMin();
155 double q_max = getBinMax();
156
157 double P_8p_theory = -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(10,q_min,q_max)/(sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(2,q_min,q_max)*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max)));
158
159 return -1./2.*P_8p_theory;
160}
161
162
163GammaPrime::GammaPrime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
164: ThObservable(SM_i)
165{
166 lep = lep_i;
167 meson = meson_i;
168 vectorM = vector_i;
169
170 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
171}
172
173double GammaPrime::computeGammaPrime(double qmin, double qmax, QCD::lepton lep)
174{
175 double q_min = qmin;
176 double q_max = qmax;
177 QCD::lepton lep_i = lep;
178
179 return ((3.*SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(0,q_min,q_max) - SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(2,q_min,q_max)) + 2.*(3.*SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(1,q_min,q_max) - SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(3,q_min,q_max)))/4.;
180}
181
182double GammaPrime::computeThValue()
183{
184 double q_min = getBinMin();
185 double q_max = getBinMax();
186
187 return computeGammaPrime(q_min, q_max, lep);
188}
189
190
191A_FB::A_FB(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
192: GammaPrime(SM_i, meson_i, vector_i, lep_i)
193{
194 lep = lep_i;
195 meson = meson_i;
196 vectorM = vector_i;
197
198 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
199}
200
201double A_FB::computeThValue()
202{
203 double q_min = getBinMin();
204 double q_max = getBinMax();
205
206 return -3. * SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(7,q_min,q_max) / 4. / computeGammaPrime(q_min, q_max, lep);
207}
208
209
210A_T1::A_T1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
211: GammaPrime(SM_i, meson_i, vector_i, lep_i)
212{
213 lep = lep_i;
214 meson = meson_i;
215 vectorM = vector_i;
216
217 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
218}
219
220double A_T1::computeThValue()
221{
222 double q_min = getBinMin();
223 double q_max = getBinMax();
224
225 return 2. * SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(7,q_min,q_max) / (3.*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(1,q_min,q_max) - SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max));
226}
227
228
229BR_MVll::BR_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
230: GammaPrime(SM_i, meson_i, vector_i, lep_i)
231{
232 lep = lep_i;
233 meson = meson_i;
234 vectorM = vector_i;
235
236 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
237}
238
239double BR_MVll::computeThValue()
240{
241 double q_min = getBinMin();
242 double q_max = getBinMax();
243
244 double ys = SM.getMesons(QCD::B_S).getDgamma_gamma()/2.;
245
246 switch(vectorM){
247 case StandardModel::K_star:
248 case StandardModel::K_star_P:
249 return (computeGammaPrime(q_min, q_max, lep)/SM.getFlavour().getMVll(meson, vectorM, lep).getwidth() + SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigmaTree(q_min, q_max)) / ( q_max - q_min );
250 break;
251 case StandardModel::PHI:
252 return computeGammaPrime(q_min, q_max, lep)/SM.getFlavour().getMVll(meson, vectorM, lep).getwidth()/(1. - ys*ys) / ( q_max - q_min );
253 break;
254 default:
255 std::stringstream out;
256 out << vectorM;
257 throw std::runtime_error("BR_MVll: vector " + out.str() + " not implemented");
258 }
259}
260
261
262R_MVll::R_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
263: GammaPrime(SM_i, meson_i, vector_i, lep_1)
264{
265 lep1 = lep_1;
266 lep2 = lep_2;
267 meson = meson_i;
268 vectorM = vector_i;
269
270 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
271 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
272}
273
274double R_MVll::computeThValue()
275{
276 double q_min = getBinMin();
277 double q_max = getBinMax();
278
279 if (vectorM == StandardModel::PHI)
280 return computeGammaPrime(q_min, q_max, lep2) / computeGammaPrime(q_min, q_max, lep1);
281 else
282 return computeGammaPrime(q_min, q_max, lep1) / computeGammaPrime(q_min, q_max, lep2);
283}
284
285
286RL_MVll::RL_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
287: F_L(SM_i, meson_i, vector_i, lep_1)
288{
289 lep1 = lep_1;
290 lep2 = lep_2;
291 meson = meson_i;
292 vectorM = vector_i;
293
294 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
295 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
296}
297
298double RL_MVll::computeThValue()
299{
300 double q_min = getBinMin();
301 double q_max = getBinMax();
302
303 return (computeGammaPrime(q_min, q_max, lep1) * computeFL(q_min, q_max, lep1)) / (computeGammaPrime(q_min, q_max, lep2) * computeFL(q_min, q_max, lep2));
304}
305
306
307RT_MVll::RT_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
308: F_L(SM_i, meson_i, vector_i, lep_1)
309{
310 lep1 = lep_1;
311 lep2 = lep_2;
312 meson = meson_i;
313 vectorM = vector_i;
314
315 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
316 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
317}
318
319double RT_MVll::computeThValue()
320{
321 double q_min = getBinMin();
322 double q_max = getBinMax();
323
324 return (computeGammaPrime(q_min, q_max, lep1) * (1 - computeFL(q_min, q_max, lep1))) / (computeGammaPrime(q_min, q_max, lep2) * ( 1 - computeFL(q_min, q_max, lep2)));
325}
326
327
328R_6::R_6(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
329: ThObservable(SM_i)
330{
331 lep1 = lep_1;
332 lep2 = lep_2;
333 meson = meson_i;
334 vectorM = vector_i;
335
336 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
337 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
338}
339
340double R_6::computeThValue()
341{
342 double q_min = getBinMin();
343 double q_max = getBinMax();
344
345 return SM.getFlavour().getMVll(meson, vectorM, lep1).integrateSigma(7,q_min,q_max) / SM.getFlavour().getMVll(meson, vectorM, lep2).integrateSigma(7,q_min,q_max);
346}
347
348
349ACP_MVll::ACP_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
350: GammaPrime(SM_i, meson_i, vector_i, lep_i)
351{
352 lep = lep_i;
353 meson = meson_i;
354 vectorM = vector_i;
355
356 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
357}
358
359double ACP_MVll::computeThValue()
360{
361 double q_min = getBinMin();
362 double q_max = getBinMax();
363
364 return (3.*SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(0,q_min,q_max) - SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(2,q_min,q_max) + 2. * ( 3.*SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(1,q_min,q_max) - SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(3,q_min,q_max) ) )/(4.*computeGammaPrime(q_min, q_max, lep));
365}
366
367
368P3CP::P3CP(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
369: ThObservable(SM_i)
370{
371 lep = lep_i;
372 meson = meson_i;
373 vectorM = vector_i;
374
375 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
376}
377
378double P3CP::computeThValue()
379{
380 double q_min = getBinMin();
381 double q_max = getBinMax();
382
383 return - SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(11,q_min,q_max)/(4.*SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max));
384}
385
386
387F_L::F_L(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
388: GammaPrime(SM_i, meson_i, vector_i, lep_i)
389{
390 lep = lep_i;
391 meson = meson_i;
392 vectorM = vector_i;
393
394 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
395}
396
397double F_L::computeFL(double qmin, double qmax, QCD::lepton lep)
398{
399 double q_min = qmin;
400 double q_max = qmax;
401 QCD::lepton lep_i = lep;
402
403 double sigma0 = SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(0,q_min,q_max);
404 double sigma2 = SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(2,q_min,q_max);
405
406 return (3.*sigma0 - sigma2) / (4. * computeGammaPrime(q_min, q_max, lep_i)) ;
407}
408
409double F_L::computeThValue()
410{
411 double q_min = getBinMin();
412 double q_max = getBinMax();
413
414 return computeFL(q_min, q_max, lep);
415}
416
417
418M_1Prime::M_1Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
419: ThObservable(SM_i)
420{
421 lep = lep_i;
422 meson = meson_i;
423 vectorM = vector_i;
424
425 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
426}
427
428double M_1Prime::computeThValue()
429{
430 double q_min = getBinMin();
431
432 return ( SM.getFlavour().getMVll(meson, vectorM, lep).H_V_p(q_min, 0).abs2() + SM.getFlavour().getMVll(meson, vectorM, lep).H_V_m(q_min, 0).abs2() - SM.getFlavour().getMVll(meson, vectorM, lep).H_A_p(q_min, 0).abs2() - SM.getFlavour().getMVll(meson, vectorM, lep).H_A_m(q_min, 0).abs2() )/
433 ( 2.*( SM.getFlavour().getMVll(meson, vectorM, lep).H_V_p(q_min, 0).abs2() + SM.getFlavour().getMVll(meson, vectorM, lep).H_V_m(q_min, 0).abs2() + SM.getFlavour().getMVll(meson, vectorM, lep).H_A_p(q_min, 0).abs2() + SM.getFlavour().getMVll(meson, vectorM, lep).H_A_m(q_min, 0).abs2() ) );
434}
435
436
437M_2Prime::M_2Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
438: ThObservable(SM_i)
439{
440 lep = lep_i;
441 meson = meson_i;
442 vectorM = vector_i;
443
444 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
445}
446
447double M_2Prime::computeThValue()
448{
449 double q_min = getBinMin();
450
451 return ( q_min/(2.*SM.getFlavour().getMVll(meson, vectorM, lep).getMlep()*SM.getFlavour().getMVll(meson, vectorM, lep).getMlep())*( SM.getFlavour().getMVll(meson, vectorM, lep).H_P(q_min, 0).abs2() + SM.getFlavour().getMVll(meson, vectorM, lep).beta(q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).beta(q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).H_S(q_min, 0).abs2() ) + SM.getFlavour().getMVll(meson, vectorM, lep).H_V_0(q_min, 0).abs2() - SM.getFlavour().getMVll(meson, vectorM, lep).H_A_0(q_min, 0).abs2() )/
452 ( SM.getFlavour().getMVll(meson, vectorM, lep).H_V_0(q_min, 0).abs2() + SM.getFlavour().getMVll(meson, vectorM, lep).H_A_0(q_min, 0).abs2() );
453}
454
455
456S_1c::S_1c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
457: GammaPrime(SM_i, meson_i, vector_i, lep_i)
458{
459 lep = lep_i;
460 meson = meson_i;
461 vectorM = vector_i;
462
463 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
464}
465
466
467double S_1c::computeThValue()
468{
469 double q_min = getBinMin();
470 double q_max = getBinMax();
471
472 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(0,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
473}
474
475
476S_1s::S_1s(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
477: GammaPrime(SM_i, meson_i, vector_i, lep_i)
478{
479 lep = lep_i;
480 meson = meson_i;
481 vectorM = vector_i;
482
483 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
484}
485
486
487double S_1s::computeThValue()
488{
489 double q_min = getBinMin();
490 double q_max = getBinMax();
491
492 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(1,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
493}
494
495
496S_2c::S_2c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
497: GammaPrime(SM_i, meson_i, vector_i, lep_i)
498{
499 lep = lep_i;
500 meson = meson_i;
501 vectorM = vector_i;
502
503 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
504}
505
506
507double S_2c::computeThValue()
508{
509 double q_min = getBinMin();
510 double q_max = getBinMax();
511
512 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(2,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
513}
514
515
516S_2s::S_2s(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
517: GammaPrime(SM_i, meson_i, vector_i, lep_i)
518{
519 lep = lep_i;
520 meson = meson_i;
521 vectorM = vector_i;
522
523 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
524}
525
526
527double S_2s::computeThValue()
528{
529 double q_min = getBinMin();
530 double q_max = getBinMax();
531
532 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
533}
534
535
536S_3::S_3(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
537: GammaPrime(SM_i, meson_i, vector_i, lep_i)
538{
539 lep = lep_i;
540 meson = meson_i;
541 vectorM = vector_i;
542
543 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
544}
545
546
547double S_3::computeThValue()
548{
549 double q_min = getBinMin();
550 double q_max = getBinMax();
551
552 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(4,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
553}
554
555
556S_4::S_4(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
557: GammaPrime(SM_i, meson_i, vector_i, lep_i)
558{
559 lep = lep_i;
560 meson = meson_i;
561 vectorM = vector_i;
562
563 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
564}
565
566/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
567double S_4::computeThValue()
568{
569 double q_min = getBinMin();
570 double q_max = getBinMax();
571
572 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(5,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
573}
574
575
576S_5::S_5(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
577: GammaPrime(SM_i, meson_i, vector_i, lep_i)
578{
579 lep = lep_i;
580 meson = meson_i;
581 vectorM = vector_i;
582
583 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
584}
585
586
587double S_5::computeThValue()
588{
589 double q_min = getBinMin();
590 double q_max = getBinMax();
591
592 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(6,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
593}
594
595
596S_6c::S_6c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
597: GammaPrime(SM_i, meson_i, vector_i, lep_i)
598{
599 lep = lep_i;
600 meson = meson_i;
601 vectorM = vector_i;
602
603 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
604}
605
606
607double S_6c::computeThValue()
608{
609 double q_min = getBinMin();
610 double q_max = getBinMax();
611
612 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(8,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
613}
614
615
616S_7::S_7(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
617: GammaPrime(SM_i, meson_i, vector_i, lep_i)
618{
619 lep = lep_i;
620 meson = meson_i;
621 vectorM = vector_i;
622
623 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
624}
625
626/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
627double S_7::computeThValue()
628{
629 double q_min = getBinMin();
630 double q_max = getBinMax();
631
632 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(9,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
633}
634
635
636S_8::S_8(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
637: GammaPrime(SM_i, meson_i, vector_i, lep_i)
638{
639 lep = lep_i;
640 meson = meson_i;
641 vectorM = vector_i;
642
643 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
644}
645
646
647double S_8::computeThValue()
648{
649 double q_min = getBinMin();
650 double q_max = getBinMax();
651
652 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(10,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
653}
654
655
656S_9::S_9(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
657: GammaPrime(SM_i, meson_i, vector_i, lep_i)
658{
659 lep = lep_i;
660 meson = meson_i;
661 vectorM = vector_i;
662
663 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
664}
665
666/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
667double S_9::computeThValue()
668{
669 double q_min = getBinMin();
670 double q_max = getBinMax();
671
672 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(11,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
673}
674
675
676A_5::A_5(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
677: GammaPrime(SM_i, meson_i, vector_i, lep_i)
678{
679 lep = lep_i;
680 meson = meson_i;
681 vectorM = vector_i;
682
683 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
684}
685
686
687double A_5::computeThValue()
688{
689 double q_min = getBinMin();
690 double q_max = getBinMax();
691
692 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(6,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
693}
694
695
696A_6::A_6(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
697: GammaPrime(SM_i, meson_i, vector_i, lep_i)
698{
699 lep = lep_i;
700 meson = meson_i;
701 vectorM = vector_i;
702
703 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
704}
705
706/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
707double A_6::computeThValue()
708{
709 double q_min = getBinMin();
710 double q_max = getBinMax();
711
712 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(7,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
713}
714
715
716A_6c::A_6c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
717: GammaPrime(SM_i, meson_i, vector_i, lep_i)
718{
719 lep = lep_i;
720 meson = meson_i;
721 vectorM = vector_i;
722
723 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
724}
725
726
727double A_6c::computeThValue()
728{
729 double q_min = getBinMin();
730 double q_max = getBinMax();
731
732 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(8,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
733}
734
735
736A_8::A_8(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
737: GammaPrime(SM_i, meson_i, vector_i, lep_i)
738{
739 lep = lep_i;
740 meson = meson_i;
741 vectorM = vector_i;
742
743 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
744}
745
746
747double A_8::computeThValue()
748{
749 double q_min = getBinMin();
750 double q_max = getBinMax();
751
752 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(10,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
753}
754
755
756A_9::A_9(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
757: GammaPrime(SM_i, meson_i, vector_i, lep_i)
758{
759 lep = lep_i;
760 meson = meson_i;
761 vectorM = vector_i;
762
763 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
764}
765
766
767double A_9::computeThValue()
768{
769 double q_min = getBinMin();
770 double q_max = getBinMax();
771
772 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(11,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
773}
774
775V0::V0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
776: ThObservable(SM_i)
777{
778 lep = lep_i;
779 meson = meson_i;
780 vectorM = vector_i;
781
782 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
783}
784
785double V0::computeThValue()
786{
787 return SM.getFlavour().getMVll(meson, vectorM, lep).getV0(getBinMin());
788}
789
790Vp::Vp(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
791: ThObservable(SM_i)
792{
793 lep = lep_i;
794 meson = meson_i;
795 vectorM = vector_i;
796
797 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
798}
799
800double Vp::computeThValue()
801{
802 return SM.getFlavour().getMVll(meson, vectorM, lep).getVp(getBinMin());
803}
804
805Vm::Vm(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
806: ThObservable(SM_i)
807{
808 lep = lep_i;
809 meson = meson_i;
810 vectorM = vector_i;
811
812 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
813}
814
815double Vm::computeThValue()
816{
817 return SM.getFlavour().getMVll(meson, vectorM, lep).getVm(getBinMin());
818}
819
820T0::T0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
821: ThObservable(SM_i)
822{
823 lep = lep_i;
824 meson = meson_i;
825 vectorM = vector_i;
826
827 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
828}
829
830double T0::computeThValue()
831{
832 return SM.getFlavour().getMVll(meson, vectorM, lep).getT0(getBinMin());
833}
834
835Tp::Tp(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
836: ThObservable(SM_i)
837{
838 lep = lep_i;
839 meson = meson_i;
840 vectorM = vector_i;
841
842 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
843}
844
845double Tp::computeThValue()
846{
847 return SM.getFlavour().getMVll(meson, vectorM, lep).getTp(getBinMin());
848}
849
850Tm::Tm(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
851: ThObservable(SM_i)
852{
853 lep = lep_i;
854 meson = meson_i;
855 vectorM = vector_i;
856
857 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
858}
859
860double Tm::computeThValue()
861{
862 return SM.getFlavour().getMVll(meson, vectorM, lep).getTm(getBinMin());
863}
864
865S::S(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
866: ThObservable(SM_i)
867{
868 lep = lep_i;
869 meson = meson_i;
870 vectorM = vector_i;
871
872 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
873}
874
875double S::computeThValue()
876{
877 return SM.getFlavour().getMVll(meson, vectorM, lep).getS(getBinMin());
878}
879
880unitarity_bound_f_F1::unitarity_bound_f_F1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
881: ThObservable(SM_i)
882{
883 lep = lep_i;
884 meson = meson_i;
885 vectorM = vector_i;
886
887 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
888}
889
890double unitarity_bound_f_F1::computeThValue()
891{
892 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_f_F1();
893}
894
895unitarity_bound_g::unitarity_bound_g(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
896: ThObservable(SM_i)
897{
898 lep = lep_i;
899 meson = meson_i;
900 vectorM = vector_i;
901
902 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
903}
904
905double unitarity_bound_g::computeThValue()
906{
907 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_g();
908}
909
910unitarity_bound_F2::unitarity_bound_F2(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
911: ThObservable(SM_i)
912{
913 lep = lep_i;
914 meson = meson_i;
915 vectorM = vector_i;
916
917 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
918}
919
920double unitarity_bound_F2::computeThValue()
921{
922 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_F2();
923}
924
925unitarity_bound_T1::unitarity_bound_T1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
926: ThObservable(SM_i)
927{
928 lep = lep_i;
929 meson = meson_i;
930 vectorM = vector_i;
931
932 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
933}
934
935double unitarity_bound_T1::computeThValue()
936{
937 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_T1();
938}
939
940unitarity_bound_T2_T0::unitarity_bound_T2_T0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
941: ThObservable(SM_i)
942{
943 lep = lep_i;
944 meson = meson_i;
945 vectorM = vector_i;
946
947 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
948}
949
950double unitarity_bound_T2_T0::computeThValue()
951{
952 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_T2_T0();
953}
954
955
956gtilde_1::gtilde_1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
957: ThObservable(SM_i)
958{
959 lep = lep_i;
960 meson = meson_i;
961 vectorM = vector_i;
962 typ = typ_i;
963
964 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
965}
966
967double gtilde_1::computeThValue()
968{
969 double q_min = getBinMin();
970
971 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_re(q_min));
972 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_im(q_min));
973 else if (typ == 3) return ((SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_im(q_min)).abs());
974 else if (typ == 4) return ((SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_im(q_min)).arg());
975 else throw std::runtime_error("MVllObservables::gtilde_1: incorrect type");
976}
977
978
979gtilde_2::gtilde_2(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
980: ThObservable(SM_i)
981{
982 lep = lep_i;
983 meson = meson_i;
984 vectorM = vector_i;
985 typ = typ_i;
986
987 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
988}
989
990double gtilde_2::computeThValue()
991{
992 double q_min = getBinMin();
993
994 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_re(q_min));
995 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_im(q_min));
996 else if (typ == 3) return ((SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_im(q_min)).abs());
997 else if (typ == 4) return ((SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_im(q_min)).arg());
998 else throw std::runtime_error("MVllObservables::gtilde_2: incorrect type");
999}
1000
1001
1002gtilde_3::gtilde_3(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1003: ThObservable(SM_i)
1004{
1005 lep = lep_i;
1006 meson = meson_i;
1007 vectorM = vector_i;
1008 typ = typ_i;
1009
1010 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1011}
1012
1013double gtilde_3::computeThValue()
1014{
1015 double q_min = getBinMin();
1016
1017 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_re(q_min));
1018 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_im(q_min));
1019 else if (typ == 3) return ((SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_im(q_min)).abs());
1020 else if (typ == 4) return ((SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_im(q_min)).arg());
1021 else throw std::runtime_error("MVllObservables::gtilde_3: incorrect type");
1022}
1023
1024h_0::h_0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1025: ThObservable(SM_i)
1026{
1027 lep = lep_i;
1028 meson = meson_i;
1029 vectorM = vector_i;
1030 typ = typ_i;
1031
1032 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1033}
1034
1035double h_0::computeThValue()
1036{
1037 double q_min = getBinMin();
1038
1039 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_re(q_min));
1040 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_im(q_min));
1041 else if (typ == 3) return ((SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_im(q_min)).abs());
1042 else if (typ == 4) return ((SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_im(q_min)).arg());
1043 else throw std::runtime_error("MVllObservables::h_0: type can only be 1:real, 2:imaginary, 3:absolute and 4:argument");
1044}
1045
1046
1047h_p::h_p(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1048: ThObservable(SM_i)
1049{
1050 lep = lep_i;
1051 meson = meson_i;
1052 vectorM = vector_i;
1053 typ = typ_i;
1054
1055 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1056}
1057
1058double h_p::computeThValue()
1059{
1060 double q_min = getBinMin();
1061
1062 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_re(q_min));
1063 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_im(q_min));
1064 else if (typ == 3) return ((SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_im(q_min)).abs());
1065 else if (typ == 4) return ((SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_im(q_min)).arg());
1066 else throw std::runtime_error("MVllObservables::h_p: type can only be 1:real, 2:imaginary, 3:absolute and 4:argument");
1067}
1068
1069
1070h_m::h_m(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1071: ThObservable(SM_i)
1072{
1073 lep = lep_i;
1074 meson = meson_i;
1075 vectorM = vector_i;
1076 typ = typ_i;
1077
1078 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1079}
1080
1081double h_m::computeThValue()
1082{
1083 double q_min = getBinMin();
1084
1085 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_re(q_min));
1086 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_im(q_min));
1087 else if (typ == 3) return ((SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_im(q_min)).abs());
1088 else if (typ == 4) return ((SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_re(q_min) + gslpp::complex::i() * SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_im(q_min)).arg());
1089 else throw std::runtime_error("MVllObservables::h_m: type can only be 1:real, 2:imaginary, 3:absolute and 4:argument");
1090}
1091
1092/***********************************************************************************************************************************
1093FUNCTIONAL
1094***********************************************************************************************************************************/
1095
1096P_1f::P_1f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1097: ThObservable(SM_i)
1098{
1099 lep = lep_i;
1100 meson = meson_i;
1101 vectorM = vector_i;
1102
1103 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1104}
1105
1106double P_1f::computeThValue()
1107{
1108 double q_min = getBinMin();
1109
1110 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min)/(2.* SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1111}
1112
1113
1114P_2f::P_2f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1115: ThObservable(SM_i)
1116{
1117 lep = lep_i;
1118 meson = meson_i;
1119 vectorM = vector_i;
1120
1121 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1122}
1123
1124/*Returns experimental value, defined according to 1510.04239.*/
1125double P_2f::computeThValue()
1126{
1127 double q_min = getBinMin();
1128
1129 double P_2f_theory = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(7,q_min)/(8.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1130
1131 return -P_2f_theory;
1132}
1133
1134
1135P_3f::P_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1136: ThObservable(SM_i)
1137{
1138 lep = lep_i;
1139 meson = meson_i;
1140 vectorM = vector_i;
1141
1142 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1143}
1144
1145/*Returns experimental value, defined according to 1510.04239.*/
1146double P_3f::computeThValue()
1147{
1148 double q_min = getBinMin();
1149
1150 double P_3f_theory = -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(11,q_min)/(4.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1151
1152 return -P_3f_theory;
1153}
1154
1155
1156P_4Primef::P_4Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1157: ThObservable(SM_i)
1158{
1159 lep = lep_i;
1160 meson = meson_i;
1161 vectorM = vector_i;
1162
1163 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1164}
1165
1166/*Returns experimental value, defined according to 1510.04239.*/
1167double P_4Primef::computeThValue()
1168{
1169 double q_min = getBinMin();
1170
1171 double P_4pf_theory = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(5,q_min)/sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1172
1173 return -1./2.*P_4pf_theory;
1174}
1175
1176
1177P_5Primef::P_5Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1178: ThObservable(SM_i)
1179{
1180 lep = lep_i;
1181 meson = meson_i;
1182 vectorM = vector_i;
1183
1184 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1185}
1186
1187double P_5Primef::computeThValue()
1188{
1189 double q_min = getBinMin();
1190
1191 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(6,q_min)/(2.*sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min)));
1192}
1193
1194
1195P_6Primef::P_6Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1196: ThObservable(SM_i)
1197{
1198 lep = lep_i;
1199 meson = meson_i;
1200 vectorM = vector_i;
1201
1202 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1203}
1204
1205double P_6Primef::computeThValue()
1206{
1207 double q_min = getBinMin();
1208
1209 return -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(9,q_min)/(2.*sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min)));
1210}
1211
1212/*Returns experimental value, defined according to 1510.04239.*/
1213P_8Primef::P_8Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1214: ThObservable(SM_i)
1215{
1216 lep = lep_i;
1217 meson = meson_i;
1218 vectorM = vector_i;
1219
1220 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1221}
1222
1223double P_8Primef::computeThValue()
1224{
1225 double q_min = getBinMin();
1226
1227 double P_8pf_theory = -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(10,q_min)/(sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min)));
1228
1229 return -1./2.*P_8pf_theory;
1230}
1231
1232
1233GammaPrimef::GammaPrimef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1234: ThObservable(SM_i)
1235{
1236 lep = lep_i;
1237 meson = meson_i;
1238 vectorM = vector_i;
1239
1240 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1241}
1242
1243double GammaPrimef::computeGammaPrimef(double qmin, QCD::lepton lep)
1244{
1245 double q_min = qmin;
1246 QCD::lepton lep_i = lep;
1247
1248 return ((3.*SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(0,q_min) - SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(2,q_min)) + 2.*(3.*SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(1,q_min) - SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(3,q_min)))/4.;
1249}
1250
1251double GammaPrimef::computeThValue()
1252{
1253 double q_min = getBinMin();
1254
1255 return computeGammaPrimef(q_min, lep);
1256}
1257
1258
1259A_FBf::A_FBf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1260: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1261{
1262 lep = lep_i;
1263 meson = meson_i;
1264 vectorM = vector_i;
1265
1266 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1267}
1268
1269double A_FBf::computeThValue()
1270{
1271 double q_min = getBinMin();
1272
1273 return -3. * SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(7,q_min) / 4. / computeGammaPrimef(q_min, lep);
1274}
1275
1276F_Lf::F_Lf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1277: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1278{
1279 lep = lep_i;
1280 meson = meson_i;
1281 vectorM = vector_i;
1282
1283 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1284}
1285
1286double F_Lf::computeFLf(double qmin, QCD::lepton lep)
1287{
1288 double q_min = qmin;
1289 QCD::lepton lep_i = lep;
1290
1291 double sigma0 = SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(0,q_min);
1292 double sigma2 = SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(2,q_min);
1293
1294 return (3.*sigma0 - sigma2) / (4. * computeGammaPrimef(q_min, lep_i)) ;
1295}
1296
1297
1298double F_Lf::computeThValue()
1299{
1300 double q_min = getBinMin();
1301
1302 return computeFLf(q_min, lep);
1303}
1304
1305S_3f::S_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1306: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1307{
1308 lep = lep_i;
1309 meson = meson_i;
1310 vectorM = vector_i;
1311
1312 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1313}
1314
1315
1316double S_3f::computeThValue()
1317{
1318 double q_min = getBinMin();
1319
1320 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min) / computeGammaPrimef(q_min, lep);
1321}
1322
1323
1324S_4f::S_4f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1325: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1326{
1327 lep = lep_i;
1328 meson = meson_i;
1329 vectorM = vector_i;
1330
1331 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1332}
1333
1334/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
1335double S_4f::computeThValue()
1336{
1337 double q_min = getBinMin();
1338
1339 return -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(5,q_min) / computeGammaPrimef(q_min, lep);
1340}
1341
1342
1343S_5f::S_5f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1344: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1345{
1346 lep = lep_i;
1347 meson = meson_i;
1348 vectorM = vector_i;
1349
1350 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1351}
1352
1353
1354double S_5f::computeThValue()
1355{
1356 double q_min = getBinMin();
1357
1358 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(6,q_min) / computeGammaPrimef(q_min, lep);
1359}
1360
1361
1362S_7f::S_7f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1363: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1364{
1365 lep = lep_i;
1366 meson = meson_i;
1367 vectorM = vector_i;
1368
1369 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1370}
1371
1372/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
1373double S_7f::computeThValue()
1374{
1375 double q_min = getBinMin();
1376
1377 return -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(9,q_min) / computeGammaPrimef(q_min, lep);
1378}
1379
1380
1381S_8f::S_8f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1382: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1383{
1384 lep = lep_i;
1385 meson = meson_i;
1386 vectorM = vector_i;
1387
1388 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1389}
1390
1391double S_8f::computeThValue()
1392{
1393 double q_min = getBinMin();
1394
1395 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(10,q_min) / computeGammaPrimef(q_min, lep);
1396}
1397
1398
1399S_9f::S_9f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1400: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1401{
1402 lep = lep_i;
1403 meson = meson_i;
1404 vectorM = vector_i;
1405
1406 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1407}
1408
1409/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
1410double S_9f::computeThValue()
1411{
1412 double q_min = getBinMin();
1413
1414 return -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(11,q_min) / computeGammaPrimef(q_min, lep);
1415}
1416
1417BRf_MVll::BRf_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1418: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1419{
1420 lep = lep_i;
1421 meson = meson_i;
1422 vectorM = vector_i;
1423
1424 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1425}
1426
1427double BRf_MVll::computeThValue()
1428{
1429 double q_min = getBinMin();
1430
1431 return computeGammaPrimef(q_min, lep)/SM.getFlavour().getMVll(meson, vectorM, lep).getwidth();
1432}
1433
1434P_relationf::P_relationf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1435: ThObservable(SM_i)
1436{
1437 lep = lep_i;
1438 meson = meson_i;
1439 vectorM = vector_i;
1440
1441 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1442}
1443
1444double P_relationf::computeThValue()
1445{
1446 double q_min = getBinMin();
1447 double P1 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min)/(2.* SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1448 double P2 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(7,q_min)/(8.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1449 double P4p = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(5,q_min)/sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1450 double P5p = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(6,q_min)/(2.*sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min)));
1451 double beta = sqrt(1. - 4.*SM.getLeptons(lep).getMass()*SM.getLeptons(lep).getMass()/q_min/q_min);
1452
1453 return 1./2.*(P4p*P5p + 1./beta *sqrt(std::abs((-1. + P1 + P4p*P4p)*(-1. - P1 + beta*beta*P5p*P5p)))) - P2;
1454}
1455
1456P_relation_exactf::P_relation_exactf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1457: ThObservable(SM_i)
1458{
1459 lep = lep_i;
1460 meson = meson_i;
1461 vectorM = vector_i;
1462
1463 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1464}
1465
1466double P_relation_exactf::computeThValue()
1467{
1468 double q_min = getBinMin();
1469 double P1 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min)/(2.* SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1470 double P2 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(7,q_min)/(8.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1471 double P3 = -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(11,q_min)/(4.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1472 double P4p = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(5,q_min)/sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1473 double P5p = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(6,q_min)/(2.*sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min)));
1474 double P6p = -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(9,q_min)/(2.*sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min)));
1475 double P8p = -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(10,q_min)/(sqrt(-SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(2,q_min)*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min)));
1476
1477 double beta = sqrt(1. - 4.*SM.getLeptons(lep).getMass()*SM.getLeptons(lep).getMass()/q_min/q_min);
1478
1479 /* k1 and k2 not fully implemented. arXiv:1402.6855*/
1480 double k1 = 1.;
1481 double k2 = 1.;
1482
1483 double delta_1 = P6p*P8p;
1484 double delta_2 = -1. + k1*k1*k2*k2 + (1. - k1*k2)*(P4p*P4p + beta*beta*P5p*P5p) - 4.*k1*k1*P3*P3 + beta*beta*P6p*P8p*(2.*P4p*P5p + P6p*P8p) + k1*(beta*beta*P6p*(4.*P3*P5p - k2*P6p) - P8p*(4.*P3*P4p + k2*P8p));
1485 double delta_3 = (1. - k1)*P4p*P4p + beta*beta*((-1. + k1)*P5p*P5p - k1*P6p*P6p) + k1*P8p*P8p;
1486 double delta_4 = 1. - k1*k1;
1487
1488 return 1./2./k1*((P4p*P5p + delta_1) + 1./beta *sqrt(std::abs((-1. + P1 + P4p*P4p)*(-1. - P1 + beta*beta*P5p*P5p) + delta_2 + delta_3*P1 + delta_4*P1*P1))) - P2;
1489}
1490
1491QCDfC9_1f::QCDfC9_1f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1492: ThObservable(SM_i)
1493{
1494 lep = lep_i;
1495 meson = meson_i;
1496 vectorM = vector_i;
1497
1498 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1499}
1500
1501double QCDfC9_1f::computeThValue()
1502{
1503 double q2 = getBinMin();
1504 double cutoff = getBinMax();
1505
1506 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9_1(q2, cutoff);
1507}
1508
1509QCDfC9_2f::QCDfC9_2f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1510: ThObservable(SM_i)
1511{
1512 lep = lep_i;
1513 meson = meson_i;
1514 vectorM = vector_i;
1515
1516 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1517}
1518
1519double QCDfC9_2f::computeThValue()
1520{
1521 double q2 = getBinMin();
1522 double cutoff = getBinMax();
1523
1524 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9_2(q2, cutoff);
1525}
1526
1527QCDfC9_3f::QCDfC9_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1528: ThObservable(SM_i)
1529{
1530 lep = lep_i;
1531 meson = meson_i;
1532 vectorM = vector_i;
1533
1534 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1535}
1536
1537double QCDfC9_3f::computeThValue()
1538{
1539 double q2 = getBinMin();
1540 double cutoff = getBinMax();
1541
1542 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9_3(q2, cutoff);
1543}
1544
1545QCDfC9p_1f::QCDfC9p_1f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1546: ThObservable(SM_i)
1547{
1548 lep = lep_i;
1549 meson = meson_i;
1550 vectorM = vector_i;
1551
1552 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1553}
1554
1555double QCDfC9p_1f::computeThValue()
1556{
1557 double cutoff = getBinMin();
1558
1559 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9p_1(cutoff);
1560}
1561
1562QCDfC9p_2f::QCDfC9p_2f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1563: ThObservable(SM_i)
1564{
1565 lep = lep_i;
1566 meson = meson_i;
1567 vectorM = vector_i;
1568
1569 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1570}
1571
1572double QCDfC9p_2f::computeThValue()
1573{
1574 double cutoff = getBinMin();
1575
1576 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9p_2(cutoff);
1577}
1578
1579QCDfC9p_3f::QCDfC9p_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1580: ThObservable(SM_i)
1581{
1582 lep = lep_i;
1583 meson = meson_i;
1584 vectorM = vector_i;
1585
1586 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1587}
1588
1589double QCDfC9p_3f::computeThValue()
1590{
1591 double cutoff = getBinMin();
1592
1593 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9p_3(cutoff);
1594}
1595
1596BR_MVpsi::BR_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1597: ThObservable(SM_i)
1598{
1599 lep = lep_i;
1600 meson = meson_i;
1601 vectorM = vector_i;
1602
1603 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1604}
1605
1606double BR_MVpsi::computeThValue()
1607{
1608 double q2 = getBinMin();
1609 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1610 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1611 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1612 return 2.*(A2_0+A2_par+A2_perp)/SM.getFlavour().getMVll(meson, vectorM, lep).getwidth();
1613}
1614
1615BR_MVpsi_ratio::BR_MVpsi_ratio(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1616: ThObservable(SM_i)
1617{
1618 lep = lep_i;
1619 meson = meson_i;
1620 vectorM = vector_i;
1621
1622 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1623}
1624
1625double BR_MVpsi_ratio::computeThValue()
1626{
1627 double A2_0_JPSI = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.0969,0).abs2();
1628 double A2_par_JPSI = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.0969,1).abs2();
1629 double A2_perp_JPSI = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.0969,2).abs2();
1630
1631 double A2_0_PSI2S = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.6861,0).abs2();
1632 double A2_par_PSI2S = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.6861,1).abs2();
1633 double A2_perp_PSI2S = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.6861,2).abs2();
1634
1635 return (A2_0_PSI2S+A2_par_PSI2S+A2_perp_PSI2S)/(A2_0_JPSI+A2_par_JPSI+A2_perp_JPSI);
1636}
1637
1638Abs2Ampar_MVpsi::Abs2Ampar_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1639: ThObservable(SM_i)
1640{
1641 lep = lep_i;
1642 meson = meson_i;
1643 vectorM = vector_i;
1644
1645 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1646}
1647
1648double Abs2Ampar_MVpsi::computeThValue()
1649{
1650 double q2 = getBinMin();
1651 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1652 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1653 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1654 return A2_par/(A2_0+A2_par+A2_perp);
1655}
1656
1657ArgAmpar_MVpsi::ArgAmpar_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1658: ThObservable(SM_i)
1659{
1660 lep = lep_i;
1661 meson = meson_i;
1662 vectorM = vector_i;
1663
1664 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1665}
1666
1667double ArgAmpar_MVpsi::computeThValue()
1668{
1669 double q2 = getBinMin();
1670 return (SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1)/SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0)).arg();
1671}
1672
1673Abs2Amperp_MVpsi::Abs2Amperp_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1674: ThObservable(SM_i)
1675{
1676 lep = lep_i;
1677 meson = meson_i;
1678 vectorM = vector_i;
1679
1680 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1681}
1682
1683double Abs2Amperp_MVpsi::computeThValue()
1684{
1685 double q2 = getBinMin();
1686 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1687 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1688 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1689 return A2_perp/(A2_0+A2_par+A2_perp);
1690}
1691
1692ArgAmperp_MVpsi::ArgAmperp_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1693: ThObservable(SM_i)
1694{
1695 lep = lep_i;
1696 meson = meson_i;
1697 vectorM = vector_i;
1698
1699 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1700}
1701
1702double ArgAmperp_MVpsi::computeThValue()
1703{
1704 double q2 = getBinMin();
1705 return (SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2)/SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0)).arg();
1706}
1707
1708Abs2Ampzero_MVpsi::Abs2Ampzero_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1709: ThObservable(SM_i)
1710{
1711 lep = lep_i;
1712 meson = meson_i;
1713 vectorM = vector_i;
1714
1715 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1716}
1717
1718double Abs2Ampzero_MVpsi::computeThValue()
1719{
1720 double q2 = getBinMin();
1721 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1722 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1723 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1724 return A2_0/(A2_0+A2_par+A2_perp);
1725}
1726
1727Delta_C9_zExp_0::Delta_C9_zExp_0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1728: ThObservable(SM_i)
1729{
1730 lep = lep_i;
1731 meson = meson_i;
1732 vectorM = vector_i;
1733
1734 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1735}
1736
1737double Delta_C9_zExp_0::computeThValue()
1738{
1739 return SM.getFlavour().getMVll(meson, vectorM, lep).getDelta_C9_zExp_0();
1740}
1741
1742Delta_C9_zExp_p::Delta_C9_zExp_p(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1743: ThObservable(SM_i)
1744{
1745 lep = lep_i;
1746 meson = meson_i;
1747 vectorM = vector_i;
1748
1749 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1750}
1751
1752double Delta_C9_zExp_p::computeThValue()
1753{
1754 return SM.getFlavour().getMVll(meson, vectorM, lep).getDelta_C9_zExp_p();
1755}
1756
1757Delta_C9_zExp_m::Delta_C9_zExp_m(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1758: ThObservable(SM_i)
1759{
1760 lep = lep_i;
1761 meson = meson_i;
1762 vectorM = vector_i;
1763
1764 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1765}
1766
1767double Delta_C9_zExp_m::computeThValue()
1768{
1769 return SM.getFlavour().getMVll(meson, vectorM, lep).getDelta_C9_zExp_m();
1770}
A_5::meson
QCD::meson meson
Definition: MVllObservables.h:1289
A_5::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1290
A_5::lep
QCD::lepton lep
Definition: MVllObservables.h:1288
A_5::A_5
A_5(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:676
A_5::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:687
A_6::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:707
A_6::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1329
A_6::A_6
A_6(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:696
A_6::lep
QCD::lepton lep
Definition: MVllObservables.h:1327
A_6::meson
QCD::meson meson
Definition: MVllObservables.h:1328
A_6c::lep
QCD::lepton lep
Definition: MVllObservables.h:1366
A_6c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1368
A_6c::A_6c
A_6c(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:716
A_6c::meson
QCD::meson meson
Definition: MVllObservables.h:1367
A_6c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:727
A_8::A_8
A_8(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:736
A_8::meson
QCD::meson meson
Definition: MVllObservables.h:1406
A_8::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1407
A_8::lep
QCD::lepton lep
Definition: MVllObservables.h:1405
A_8::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:747
A_9::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1446
A_9::meson
QCD::meson meson
Definition: MVllObservables.h:1445
A_9::lep
QCD::lepton lep
Definition: MVllObservables.h:1444
A_9::A_9
A_9(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:756
A_9::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:767
A_FB::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:376
A_FB::A_FB
A_FB(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:191
A_FB::lep
QCD::lepton lep
Definition: MVllObservables.h:374
A_FB::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:201
A_FB::meson
QCD::meson meson
Definition: MVllObservables.h:375
A_FBf::meson
QCD::meson meson
Definition: MVllObservables.h:2457
A_FBf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1269
A_FBf::A_FBf
A_FBf(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1259
A_FBf::lep
QCD::lepton lep
Definition: MVllObservables.h:2456
A_FBf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2458
A_T1::meson
QCD::meson meson
Definition: MVllObservables.h:414
A_T1::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:220
A_T1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:415
A_T1::A_T1
A_T1(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:210
A_T1::lep
QCD::lepton lep
Definition: MVllObservables.h:413
ACP_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:704
ACP_MVll::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:359
ACP_MVll::ACP_MVll
ACP_MVll(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:349
ACP_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:703
ACP_MVll::lep
QCD::lepton lep
Definition: MVllObservables.h:702
Abs2Ampar_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3132
Abs2Ampar_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3133
Abs2Ampar_MVpsi::computeThValue
double computeThValue()
The observable Abs(A//) of .
Definition: MVllObservables.cpp:1648
Abs2Ampar_MVpsi::Abs2Ampar_MVpsi
Abs2Ampar_MVpsi(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1638
Abs2Ampar_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3131
Abs2Amperp_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3200
Abs2Amperp_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3199
Abs2Amperp_MVpsi::Abs2Amperp_MVpsi
Abs2Amperp_MVpsi(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1673
Abs2Amperp_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3201
Abs2Amperp_MVpsi::computeThValue
double computeThValue()
The observable Abs(A_|_) of .
Definition: MVllObservables.cpp:1683
Abs2Ampzero_MVpsi::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1718
Abs2Ampzero_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3267
Abs2Ampzero_MVpsi::Abs2Ampzero_MVpsi
Abs2Ampzero_MVpsi(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1708
Abs2Ampzero_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3268
Abs2Ampzero_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3269
ArgAmpar_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3167
ArgAmpar_MVpsi::computeThValue
double computeThValue()
The observable Arg(A//) of .
Definition: MVllObservables.cpp:1667
ArgAmpar_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3166
ArgAmpar_MVpsi::ArgAmpar_MVpsi
ArgAmpar_MVpsi(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1657
ArgAmpar_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3165
ArgAmperp_MVpsi::ArgAmperp_MVpsi
ArgAmperp_MVpsi(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1692
ArgAmperp_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3235
ArgAmperp_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3234
ArgAmperp_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3233
ArgAmperp_MVpsi::computeThValue
double computeThValue()
The observable Arg(A_|_) of .
Definition: MVllObservables.cpp:1702
BR_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:454
BR_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:453
BR_MVll::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:239
BR_MVll::lep
QCD::lepton lep
Definition: MVllObservables.h:452
BR_MVll::BR_MVll
BR_MVll(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:229
BR_MVpsi_ratio::meson
QCD::meson meson
Definition: MVllObservables.h:3098
BR_MVpsi_ratio::lep
QCD::lepton lep
Definition: MVllObservables.h:3097
BR_MVpsi_ratio::BR_MVpsi_ratio
BR_MVpsi_ratio(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1615
BR_MVpsi_ratio::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3099
BR_MVpsi_ratio::computeThValue
double computeThValue()
The observable BR of .
Definition: MVllObservables.cpp:1625
BR_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3064
BR_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3065
BR_MVpsi::BR_MVpsi
BR_MVpsi(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1596
BR_MVpsi::computeThValue
double computeThValue()
The observable BR of .
Definition: MVllObservables.cpp:1606
BR_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3063
BRf_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2497
BRf_MVll::lep
QCD::lepton lep
Definition: MVllObservables.h:2495
BRf_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:2496
BRf_MVll::BRf_MVll
BRf_MVll(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1417
BRf_MVll::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1427
Delta_C9_zExp_0::lep
QCD::lepton lep
Definition: MVllObservables.h:3301
Delta_C9_zExp_0::meson
QCD::meson meson
Definition: MVllObservables.h:3302
Delta_C9_zExp_0::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1737
Delta_C9_zExp_0::Delta_C9_zExp_0
Delta_C9_zExp_0(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1727
Delta_C9_zExp_0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3303
Delta_C9_zExp_m::lep
QCD::lepton lep
Definition: MVllObservables.h:3369
Delta_C9_zExp_m::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3371
Delta_C9_zExp_m::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1767
Delta_C9_zExp_m::meson
QCD::meson meson
Definition: MVllObservables.h:3370
Delta_C9_zExp_m::Delta_C9_zExp_m
Delta_C9_zExp_m(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1757
Delta_C9_zExp_p::Delta_C9_zExp_p
Delta_C9_zExp_p(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1742
Delta_C9_zExp_p::lep
QCD::lepton lep
Definition: MVllObservables.h:3335
Delta_C9_zExp_p::meson
QCD::meson meson
Definition: MVllObservables.h:3336
Delta_C9_zExp_p::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3337
Delta_C9_zExp_p::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1752
F_L
A class for the binned observable in .
Definition: MVllObservables.h:472
F_L::lep
QCD::lepton lep
Definition: MVllObservables.h:500
F_L::computeFL
double computeFL(double qmin, double qmax, QCD::lepton lep)
A method to compute the binned observable in in a generic bin.
Definition: MVllObservables.cpp:397
F_L::F_L
F_L(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:387
F_L::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:409
F_L::meson
QCD::meson meson
Definition: MVllObservables.h:501
F_L::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:502
F_Lf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2544
F_Lf::lep
QCD::lepton lep
Definition: MVllObservables.h:2542
F_Lf::meson
QCD::meson meson
Definition: MVllObservables.h:2543
F_Lf::F_Lf
F_Lf(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1276
F_Lf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1298
F_Lf::computeFLf
double computeFLf(double qmin, QCD::lepton lep)
A method to compute the observable in at a generic .
Definition: MVllObservables.cpp:1286
Flavour::getMVll
MVll & getMVll(QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i) const
Returns the initial and final state dependent object for .
Definition: Flavour.cpp:210
GammaPrime
A class for the binned observable in .
Definition: MVllObservables.h:307
GammaPrime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:337
GammaPrime::GammaPrime
GammaPrime(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:163
GammaPrime::meson
QCD::meson meson
Definition: MVllObservables.h:336
GammaPrime::lep
QCD::lepton lep
Definition: MVllObservables.h:335
GammaPrime::computeGammaPrime
double computeGammaPrime(double qmin, double qmax, QCD::lepton lep)
A method to compute the binned observable in in a generic bin.
Definition: MVllObservables.cpp:173
GammaPrime::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:182
GammaPrimef
A class for the observable in .
Definition: MVllObservables.h:2390
GammaPrimef::meson
QCD::meson meson
Definition: MVllObservables.h:2418
GammaPrimef::computeGammaPrimef
double computeGammaPrimef(double qmin, QCD::lepton lep)
A method to compute the observable in at a generic .
Definition: MVllObservables.cpp:1243
GammaPrimef::lep
QCD::lepton lep
Definition: MVllObservables.h:2417
GammaPrimef::GammaPrimef
GammaPrimef(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1233
GammaPrimef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1251
GammaPrimef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2419
M_1Prime::meson
QCD::meson meson
Definition: MVllObservables.h:782
M_1Prime::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:428
M_1Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:783
M_1Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:781
M_1Prime::M_1Prime
M_1Prime(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:418
M_2Prime::meson
QCD::meson meson
Definition: MVllObservables.h:821
M_2Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:820
M_2Prime::M_2Prime
M_2Prime(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:437
M_2Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:822
M_2Prime::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:447
MVll::getDelta_C9_zExp_p
double getDelta_C9_zExp_p()
The non-pertubative ccbar contributions to the helicity amplitudes.
Definition: MVll.h:539
MVll::getgtilde_2_re
double getgtilde_2_re(double q2)
The real part of .
Definition: MVll.h:727
MVll::getS
double getS(double q2)
The form factor .
Definition: MVll.h:454
MVll::getQCDfC9_2
double getQCDfC9_2(double q2, double cutoff)
Definition: MVll.h:670
MVll::getQCDfC9p_1
double getQCDfC9p_1(double cutoff)
Definition: MVll.h:682
MVll::getDelta_C9_zExp_m
double getDelta_C9_zExp_m()
The non-pertubative ccbar contributions to the helicity amplitudes.
Definition: MVll.h:549
MVll::geth_0_re
double geth_0_re(double q2)
The real part of .
Definition: MVll.h:773
MVll::geth_p_im
double geth_p_im(double q2)
The imaginary part of .
Definition: MVll.h:812
MVll::H_A_p
gslpp::complex H_A_p(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2569
MVll::getMlep
double getMlep()
The mass of the lepton l.
Definition: MVll.h:372
MVll::getgtilde_2_im
double getgtilde_2_im(double q2)
The immaginary part of .
Definition: MVll.h:738
MVll::geth_m_im
double geth_m_im(double q2)
The imaginary part of .
Definition: MVll.h:841
MVll::getVm
double getVm(double q2)
The form factor .
Definition: MVll.h:411
MVll::get_unitarity_bound_g
double get_unitarity_bound_g()
The unitarity constraints on form factors .
Definition: MVll.h:474
MVll::getVp
double getVp(double q2)
The form factor .
Definition: MVll.h:400
MVll::H_S
gslpp::complex H_S(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2581
MVll::initializeMVllParameters
std::vector< std::string > initializeMVllParameters()
A method for initializing the parameters necessary for MVll.
Definition: MVll.cpp:160
MVll::integrateDelta
double integrateDelta(int i, double q_min, double q_max)
The integral of from to .
Definition: MVll.cpp:3019
MVll::getgtilde_3_im
double getgtilde_3_im(double q2)
The imaginary part of .
Definition: MVll.h:761
MVll::getTm
double getTm(double q2)
The form factor .
Definition: MVll.h:443
MVll::H_V_0
gslpp::complex H_V_0(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2545
MVll::getT0
double getT0(double q2)
The form factor .
Definition: MVll.h:421
MVll::getSigma
double getSigma(int i, double q_2)
The value of from to .
Definition: MVll.cpp:2971
MVll::H_V_p
gslpp::complex H_V_p(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2551
MVll::getgtilde_3_re
double getgtilde_3_re(double q2)
The real part of .
Definition: MVll.h:749
MVll::getwidth
double getwidth()
The width of the meson M.
Definition: MVll.h:363
MVll::get_unitarity_bound_T1
double get_unitarity_bound_T1()
The unitarity constraints on form factors .
Definition: MVll.h:494
MVll::getV0
double getV0(double q2)
The form factor .
Definition: MVll.h:389
MVll::get_unitarity_bound_T2_T0
double get_unitarity_bound_T2_T0()
The unitarity constraints on form factors and .
Definition: MVll.h:504
MVll::getDelta_C9_zExp_0
double getDelta_C9_zExp_0()
The non-pertubative ccbar contributions to the helicity amplitudes.
Definition: MVll.h:529
MVll::geth_p_re
double geth_p_re(double q2)
The real part of .
Definition: MVll.h:802
MVll::integrateSigmaTree
double integrateSigmaTree(double q_min, double q_max)
The integral of from to (arxiv/2301.06990)
Definition: MVll.cpp:3118
MVll::geth_0_im
double geth_0_im(double q2)
The imaginary part of .
Definition: MVll.h:783
MVll::getTp
double getTp(double q2)
The form factor .
Definition: MVll.h:432
MVll::H_V_m
gslpp::complex H_V_m(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2557
MVll::AmpMVpsi_zExpansion
gslpp::complex AmpMVpsi_zExpansion(double mpsi, int tran)
Polarization amplitudes for M to V psi, Eq. B.16 of arXiv:2206.03797.
Definition: MVll.cpp:2639
MVll::get_unitarity_bound_F2
double get_unitarity_bound_F2()
The unitarity constraints on form factors .
Definition: MVll.h:484
MVll::geth_m_re
double geth_m_re(double q2)
The real part of .
Definition: MVll.h:831
MVll::H_A_0
gslpp::complex H_A_0(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2563
MVll::getQCDfC9p_2
double getQCDfC9p_2(double cutoff)
Definition: MVll.h:688
MVll::getQCDfC9p_3
double getQCDfC9p_3(double cutoff)
Definition: MVll.h:694
MVll::getgtilde_1_im
double getgtilde_1_im(double q2)
The immaginary part of .
Definition: MVll.h:716
MVll::getQCDfC9_1
double getQCDfC9_1(double q2, double cutoff)
Definition: MVll.h:664
MVll::H_A_m
gslpp::complex H_A_m(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2575
MVll::H_P
gslpp::complex H_P(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2589
MVll::integrateSigma
double integrateSigma(int i, double q_min, double q_max)
The integral of from to .
Definition: MVll.cpp:2844
MVll::getQCDfC9_3
double getQCDfC9_3(double q2, double cutoff)
Definition: MVll.h:676
MVll::beta
double beta(double q2)
The factor used in the angular coefficients .
Definition: MVll.cpp:2671
MVll::getgtilde_1_re
double getgtilde_1_re(double q2)
The real part of .
Definition: MVll.h:705
MVll::get_unitarity_bound_f_F1
double get_unitarity_bound_f_F1()
The unitarity constraints on form factors and .
Definition: MVll.h:464
Meson::getDgamma_gamma
const double & getDgamma_gamma() const
Definition: Meson.h:411
P3CP::lep
QCD::lepton lep
Definition: MVllObservables.h:742
P3CP::P3CP
P3CP(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:368
P3CP::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:378
P3CP::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:744
P3CP::meson
QCD::meson meson
Definition: MVllObservables.h:743
P_1::meson
QCD::meson meson
Definition: MVllObservables.h:49
P_1::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:29
P_1::P_1
P_1(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:19
P_1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:50
P_1::lep
QCD::lepton lep
Definition: MVllObservables.h:48
P_1f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2133
P_1f::lep
QCD::lepton lep
Definition: MVllObservables.h:2131
P_1f::P_1f
P_1f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1096
P_1f::meson
QCD::meson meson
Definition: MVllObservables.h:2132
P_1f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1106
P_2::P_2
P_2(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:38
P_2::lep
QCD::lepton lep
Definition: MVllObservables.h:87
P_2::meson
QCD::meson meson
Definition: MVllObservables.h:88
P_2::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:89
P_2::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:48
P_2f::P_2f
P_2f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1114
P_2f::meson
QCD::meson meson
Definition: MVllObservables.h:2171
P_2f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2172
P_2f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1125
P_2f::lep
QCD::lepton lep
Definition: MVllObservables.h:2170
P_3::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:128
P_3::meson
QCD::meson meson
Definition: MVllObservables.h:127
P_3::P_3
P_3(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:59
P_3::lep
QCD::lepton lep
Definition: MVllObservables.h:126
P_3::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:70
P_3f::meson
QCD::meson meson
Definition: MVllObservables.h:2210
P_3f::P_3f
P_3f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1135
P_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2211
P_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:2209
P_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1146
P_4Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:170
P_4Prime::meson
QCD::meson meson
Definition: MVllObservables.h:169
P_4Prime::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:92
P_4Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:168
P_4Prime::P_4Prime
P_4Prime(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:81
P_4Primef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1167
P_4Primef::P_4Primef
P_4Primef(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1156
P_4Primef::lep
QCD::lepton lep
Definition: MVllObservables.h:2251
P_4Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2253
P_4Primef::meson
QCD::meson meson
Definition: MVllObservables.h:2252
P_5Prime::meson
QCD::meson meson
Definition: MVllObservables.h:210
P_5Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:211
P_5Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:209
P_5Prime::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:113
P_5Prime::P_5Prime
P_5Prime(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:103
P_5Primef::meson
QCD::meson meson
Definition: MVllObservables.h:2293
P_5Primef::P_5Primef
P_5Primef(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1177
P_5Primef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1187
P_5Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2294
P_5Primef::lep
QCD::lepton lep
Definition: MVllObservables.h:2292
P_6Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:248
P_6Prime::meson
QCD::meson meson
Definition: MVllObservables.h:249
P_6Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:250
P_6Prime::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:132
P_6Prime::P_6Prime
P_6Prime(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:122
P_6Primef::lep
QCD::lepton lep
Definition: MVllObservables.h:2331
P_6Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2333
P_6Primef::P_6Primef
P_6Primef(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1195
P_6Primef::meson
QCD::meson meson
Definition: MVllObservables.h:2332
P_6Primef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1205
P_8Prime::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:152
P_8Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:287
P_8Prime::meson
QCD::meson meson
Definition: MVllObservables.h:288
P_8Prime::P_8Prime
P_8Prime(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:141
P_8Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:289
P_8Primef::P_8Primef
P_8Primef(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1213
P_8Primef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1223
P_8Primef::meson
QCD::meson meson
Definition: MVllObservables.h:2371
P_8Primef::lep
QCD::lepton lep
Definition: MVllObservables.h:2370
P_8Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2372
P_relation_exactf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1466
P_relation_exactf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2827
P_relation_exactf::meson
QCD::meson meson
Definition: MVllObservables.h:2826
P_relation_exactf::lep
QCD::lepton lep
Definition: MVllObservables.h:2825
P_relation_exactf::P_relation_exactf
P_relation_exactf(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1456
P_relationf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2802
P_relationf::P_relationf
P_relationf(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1434
P_relationf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1444
P_relationf::meson
QCD::meson meson
Definition: MVllObservables.h:2801
P_relationf::lep
QCD::lepton lep
Definition: MVllObservables.h:2800
Particle::getMass
const double & getMass() const
A get method to access the particle mass.
Definition: Particle.h:61
QCD::meson
meson
An enum type for mesons.
Definition: QCD.h:336
QCD::PHI
@ PHI
Definition: QCD.h:348
QCD::K_star
@ K_star
Definition: QCD.h:349
QCD::K_star_P
@ K_star_P
Definition: QCD.h:350
QCD::B_S
@ B_S
Definition: QCD.h:346
QCD::getMesons
const Meson & getMesons(const QCD::meson m) const
A get method to access a meson as an object of the type Meson.
Definition: QCD.h:526
QCD::lepton
lepton
An enum type for leptons.
Definition: QCD.h:310
QCDfC9_1f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2861
QCDfC9_1f::lep
QCD::lepton lep
Definition: MVllObservables.h:2859
QCDfC9_1f::meson
QCD::meson meson
Definition: MVllObservables.h:2860
QCDfC9_1f::QCDfC9_1f
QCDfC9_1f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1491
QCDfC9_1f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1501
QCDfC9_2f::lep
QCD::lepton lep
Definition: MVllObservables.h:2893
QCDfC9_2f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2895
QCDfC9_2f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1519
QCDfC9_2f::QCDfC9_2f
QCDfC9_2f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1509
QCDfC9_2f::meson
QCD::meson meson
Definition: MVllObservables.h:2894
QCDfC9_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1537
QCDfC9_3f::QCDfC9_3f
QCDfC9_3f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1527
QCDfC9_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:2927
QCDfC9_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2929
QCDfC9_3f::meson
QCD::meson meson
Definition: MVllObservables.h:2928
QCDfC9p_1f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2963
QCDfC9p_1f::lep
QCD::lepton lep
Definition: MVllObservables.h:2961
QCDfC9p_1f::meson
QCD::meson meson
Definition: MVllObservables.h:2962
QCDfC9p_1f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1555
QCDfC9p_1f::QCDfC9p_1f
QCDfC9p_1f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1545
QCDfC9p_2f::QCDfC9p_2f
QCDfC9p_2f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1562
QCDfC9p_2f::meson
QCD::meson meson
Definition: MVllObservables.h:2996
QCDfC9p_2f::lep
QCD::lepton lep
Definition: MVllObservables.h:2995
QCDfC9p_2f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1572
QCDfC9p_2f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2997
QCDfC9p_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:3029
QCDfC9p_3f::QCDfC9p_3f
QCDfC9p_3f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1579
QCDfC9p_3f::meson
QCD::meson meson
Definition: MVllObservables.h:3030
QCDfC9p_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3031
QCDfC9p_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1589
R_6::R_6
R_6(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
Constructor.
Definition: MVllObservables.cpp:328
R_6::meson
QCD::meson meson
Definition: MVllObservables.h:665
R_6::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:666
R_6::lep1
QCD::lepton lep1
Definition: MVllObservables.h:663
R_6::lep2
QCD::lepton lep2
Definition: MVllObservables.h:664
R_6::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:340
R_MVll::lep1
QCD::lepton lep1
Definition: MVllObservables.h:540
R_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:542
R_MVll::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:274
R_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:543
R_MVll::R_MVll
R_MVll(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
Constructor.
Definition: MVllObservables.cpp:262
R_MVll::lep2
QCD::lepton lep2
Definition: MVllObservables.h:541
RL_MVll::lep2
QCD::lepton lep2
Definition: MVllObservables.h:582
RL_MVll::lep1
QCD::lepton lep1
Definition: MVllObservables.h:581
RL_MVll::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:298
RL_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:584
RL_MVll::RL_MVll
RL_MVll(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
Constructor.
Definition: MVllObservables.cpp:286
RL_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:583
RT_MVll::lep1
QCD::lepton lep1
Definition: MVllObservables.h:622
RT_MVll::lep2
QCD::lepton lep2
Definition: MVllObservables.h:623
RT_MVll::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:319
RT_MVll::RT_MVll
RT_MVll(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
Constructor.
Definition: MVllObservables.cpp:307
RT_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:625
RT_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:624
S_1c::meson
QCD::meson meson
Definition: MVllObservables.h:860
S_1c::S_1c
S_1c(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:456
S_1c::lep
QCD::lepton lep
Definition: MVllObservables.h:859
S_1c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:861
S_1c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:467
S_1s::lep
QCD::lepton lep
Definition: MVllObservables.h:898
S_1s::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:900
S_1s::S_1s
S_1s(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:476
S_1s::meson
QCD::meson meson
Definition: MVllObservables.h:899
S_1s::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:487
S_2c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:939
S_2c::lep
QCD::lepton lep
Definition: MVllObservables.h:937
S_2c::meson
QCD::meson meson
Definition: MVllObservables.h:938
S_2c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:507
S_2c::S_2c
S_2c(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:496
S_2s::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:978
S_2s::meson
QCD::meson meson
Definition: MVllObservables.h:977
S_2s::S_2s
S_2s(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:516
S_2s::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:527
S_2s::lep
QCD::lepton lep
Definition: MVllObservables.h:976
S_3::S_3
S_3(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:536
S_3::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1017
S_3::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:547
S_3::meson
QCD::meson meson
Definition: MVllObservables.h:1016
S_3::lep
QCD::lepton lep
Definition: MVllObservables.h:1015
S_3f::meson
QCD::meson meson
Definition: MVllObservables.h:2581
S_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1316
S_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2582
S_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:2580
S_3f::S_3f
S_3f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1305
S_4::lep
QCD::lepton lep
Definition: MVllObservables.h:1054
S_4::meson
QCD::meson meson
Definition: MVllObservables.h:1055
S_4::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1056
S_4::S_4
S_4(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:556
S_4::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:567
S_4f::lep
QCD::lepton lep
Definition: MVllObservables.h:2619
S_4f::meson
QCD::meson meson
Definition: MVllObservables.h:2620
S_4f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2621
S_4f::S_4f
S_4f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1324
S_4f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1335
S_5::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:587
S_5::S_5
S_5(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:576
S_5::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1095
S_5::lep
QCD::lepton lep
Definition: MVllObservables.h:1093
S_5::meson
QCD::meson meson
Definition: MVllObservables.h:1094
S_5f::lep
QCD::lepton lep
Definition: MVllObservables.h:2658
S_5f::S_5f
S_5f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1343
S_5f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1354
S_5f::meson
QCD::meson meson
Definition: MVllObservables.h:2659
S_5f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2660
S_6c::S_6c
S_6c(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:596
S_6c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1134
S_6c::meson
QCD::meson meson
Definition: MVllObservables.h:1133
S_6c::lep
QCD::lepton lep
Definition: MVllObservables.h:1132
S_6c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:607
S_7::S_7
S_7(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:616
S_7::lep
QCD::lepton lep
Definition: MVllObservables.h:1171
S_7::meson
QCD::meson meson
Definition: MVllObservables.h:1172
S_7::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1173
S_7::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:627
S_7f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2699
S_7f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1373
S_7f::S_7f
S_7f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1362
S_7f::meson
QCD::meson meson
Definition: MVllObservables.h:2698
S_7f::lep
QCD::lepton lep
Definition: MVllObservables.h:2697
S_8::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1212
S_8::S_8
S_8(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:636
S_8::meson
QCD::meson meson
Definition: MVllObservables.h:1211
S_8::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:647
S_8::lep
QCD::lepton lep
Definition: MVllObservables.h:1210
S_8f::meson
QCD::meson meson
Definition: MVllObservables.h:2737
S_8f::S_8f
S_8f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1381
S_8f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2738
S_8f::lep
QCD::lepton lep
Definition: MVllObservables.h:2736
S_8f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1391
S_9::lep
QCD::lepton lep
Definition: MVllObservables.h:1249
S_9::S_9
S_9(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:656
S_9::meson
QCD::meson meson
Definition: MVllObservables.h:1250
S_9::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1251
S_9::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:667
S_9f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1410
S_9f::S_9f
S_9f(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1399
S_9f::lep
QCD::lepton lep
Definition: MVllObservables.h:2775
S_9f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2777
S_9f::meson
QCD::meson meson
Definition: MVllObservables.h:2776
S::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1691
S::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:875
S::S
S(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:865
S::meson
QCD::meson meson
Definition: MVllObservables.h:1690
S::lep
QCD::lepton lep
Definition: MVllObservables.h:1689
StandardModel
A model class for the Standard Model.
Definition: StandardModel.h:521
StandardModel::getLeptons
const Particle & getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
Definition: StandardModel.h:756
StandardModel::getFlavour
const Flavour & getFlavour() const
Definition: StandardModel.h:1065
T0::lep
QCD::lepton lep
Definition: MVllObservables.h:1584
T0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1586
T0::T0
T0(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:820
T0::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:830
T0::meson
QCD::meson meson
Definition: MVllObservables.h:1585
ThObservable
A class for a model prediction of an observable.
Definition: ThObservable.h:25
ThObservable::setParametersForObservable
void setParametersForObservable(std::vector< std::string > parametersForObservable_i)
A set method to get the parameters for the specific observable.
Definition: ThObservable.h:109
ThObservable::getBinMin
double getBinMin()
A get method to get the minimum value of the bin.
Definition: ThObservable.h:82
ThObservable::getBinMax
double getBinMax()
A get method to get the maximum value of the bin.
Definition: ThObservable.h:91
ThObservable::SM
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121
Tm::lep
QCD::lepton lep
Definition: MVllObservables.h:1654
Tm::meson
QCD::meson meson
Definition: MVllObservables.h:1655
Tm::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:860
Tm::Tm
Tm(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:850
Tm::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1656
Tp::meson
QCD::meson meson
Definition: MVllObservables.h:1620
Tp::Tp
Tp(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:835
Tp::lep
QCD::lepton lep
Definition: MVllObservables.h:1619
Tp::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1621
Tp::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:845
V0::V0
V0(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:775
V0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1481
V0::lep
QCD::lepton lep
Definition: MVllObservables.h:1479
V0::meson
QCD::meson meson
Definition: MVllObservables.h:1480
V0::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:785
Vm::Vm
Vm(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:805
Vm::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1551
Vm::meson
QCD::meson meson
Definition: MVllObservables.h:1550
Vm::lep
QCD::lepton lep
Definition: MVllObservables.h:1549
Vm::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:815
Vp::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1516
Vp::Vp
Vp(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:790
Vp::meson
QCD::meson meson
Definition: MVllObservables.h:1515
Vp::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:800
Vp::lep
QCD::lepton lep
Definition: MVllObservables.h:1514
gtilde_1::typ
unsigned int typ
Definition: MVllObservables.h:1898
gtilde_1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1897
gtilde_1::computeThValue
double computeThValue()
The averaged correction in .
Definition: MVllObservables.cpp:967
gtilde_1::gtilde_1
gtilde_1(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
Constructor.
Definition: MVllObservables.cpp:956
gtilde_1::meson
QCD::meson meson
Definition: MVllObservables.h:1896
gtilde_1::lep
QCD::lepton lep
Definition: MVllObservables.h:1895
gtilde_2::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1936
gtilde_2::lep
QCD::lepton lep
Definition: MVllObservables.h:1934
gtilde_2::typ
unsigned int typ
Definition: MVllObservables.h:1937
gtilde_2::meson
QCD::meson meson
Definition: MVllObservables.h:1935
gtilde_2::computeThValue
double computeThValue()
The averaged correction in .
Definition: MVllObservables.cpp:990
gtilde_2::gtilde_2
gtilde_2(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
Constructor.
Definition: MVllObservables.cpp:979
gtilde_3::computeThValue
double computeThValue()
The averaged correction in .
Definition: MVllObservables.cpp:1013
gtilde_3::meson
QCD::meson meson
Definition: MVllObservables.h:1974
gtilde_3::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1975
gtilde_3::gtilde_3
gtilde_3(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
Constructor.
Definition: MVllObservables.cpp:1002
gtilde_3::typ
unsigned int typ
Definition: MVllObservables.h:1976
gtilde_3::lep
QCD::lepton lep
Definition: MVllObservables.h:1973
h_0::h_0
h_0(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
Constructor.
Definition: MVllObservables.cpp:1024
h_0::meson
QCD::meson meson
Definition: MVllObservables.h:2011
h_0::typ
unsigned int typ
Definition: MVllObservables.h:2013
h_0::lep
QCD::lepton lep
Definition: MVllObservables.h:2010
h_0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2012
h_0::computeThValue
double computeThValue()
The correction in .
Definition: MVllObservables.cpp:1035
h_m::h_m
h_m(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
Constructor.
Definition: MVllObservables.cpp:1070
h_m::typ
unsigned int typ
Definition: MVllObservables.h:2089
h_m::computeThValue
double computeThValue()
The correction in .
Definition: MVllObservables.cpp:1081
h_m::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2088
h_m::lep
QCD::lepton lep
Definition: MVllObservables.h:2086
h_m::meson
QCD::meson meson
Definition: MVllObservables.h:2087
h_p::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2050
h_p::typ
unsigned int typ
Definition: MVllObservables.h:2051
h_p::meson
QCD::meson meson
Definition: MVllObservables.h:2049
h_p::lep
QCD::lepton lep
Definition: MVllObservables.h:2048
h_p::h_p
h_p(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
Constructor.
Definition: MVllObservables.cpp:1047
h_p::computeThValue
double computeThValue()
The correction in .
Definition: MVllObservables.cpp:1058
unitarity_bound_F2::unitarity_bound_F2
unitarity_bound_F2(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:910
unitarity_bound_F2::computeThValue
double computeThValue()
The unitarity bound on form factors in decays.
Definition: MVllObservables.cpp:920
unitarity_bound_F2::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1793
unitarity_bound_F2::meson
QCD::meson meson
Definition: MVllObservables.h:1792
unitarity_bound_F2::lep
QCD::lepton lep
Definition: MVllObservables.h:1791
unitarity_bound_T1::computeThValue
double computeThValue()
The unitarity bound on form factors in decays.
Definition: MVllObservables.cpp:935
unitarity_bound_T1::meson
QCD::meson meson
Definition: MVllObservables.h:1826
unitarity_bound_T1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1827
unitarity_bound_T1::unitarity_bound_T1
unitarity_bound_T1(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:925
unitarity_bound_T1::lep
QCD::lepton lep
Definition: MVllObservables.h:1825
unitarity_bound_T2_T0::unitarity_bound_T2_T0
unitarity_bound_T2_T0(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:940
unitarity_bound_T2_T0::computeThValue
double computeThValue()
The unitarity bound on form factors and in decays.
Definition: MVllObservables.cpp:950
unitarity_bound_T2_T0::lep
QCD::lepton lep
Definition: MVllObservables.h:1859
unitarity_bound_T2_T0::meson
QCD::meson meson
Definition: MVllObservables.h:1860
unitarity_bound_T2_T0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1861
unitarity_bound_f_F1::lep
QCD::lepton lep
Definition: MVllObservables.h:1723
unitarity_bound_f_F1::computeThValue
double computeThValue()
The unitarity bound on form factors and in decays.
Definition: MVllObservables.cpp:890
unitarity_bound_f_F1::unitarity_bound_f_F1
unitarity_bound_f_F1(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:880
unitarity_bound_f_F1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1725
unitarity_bound_f_F1::meson
QCD::meson meson
Definition: MVllObservables.h:1724
unitarity_bound_g::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1759
unitarity_bound_g::lep
QCD::lepton lep
Definition: MVllObservables.h:1757
unitarity_bound_g::computeThValue
double computeThValue()
The unitarity bound on form factors in decays.
Definition: MVllObservables.cpp:905
unitarity_bound_g::unitarity_bound_g
unitarity_bound_g(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:895
unitarity_bound_g::meson
QCD::meson meson
Definition: MVllObservables.h:1758