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
262BR_MVll_tot::BR_MVll_tot(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
263: GammaPrime(SM_i, meson_i, vector_i, lep_i)
264{
265 lep = lep_i;
266 meson = meson_i;
267 vectorM = vector_i;
268
269 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
270}
271
272double BR_MVll_tot::computeThValue()
273{
274 double q_min = getBinMin();
275 double q_max = getBinMax();
276
277 double ys = SM.getMesons(QCD::B_S).getDgamma_gamma()/2.;
278
279 switch(vectorM){
280 case StandardModel::K_star:
281 case StandardModel::K_star_P:
282 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));
283 break;
284 case StandardModel::PHI:
285 return computeGammaPrime(q_min, q_max, lep)/SM.getFlavour().getMVll(meson, vectorM, lep).getwidth()/(1. - ys*ys) ;
286 break;
287 default:
288 std::stringstream out;
289 out << vectorM;
290 throw std::runtime_error("BR_MVll_tot: vector " + out.str() + " not implemented");
291 }
292}
293
294
295R_MVll::R_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
296: GammaPrime(SM_i, meson_i, vector_i, lep_1)
297{
298 lep1 = lep_1;
299 lep2 = lep_2;
300 meson = meson_i;
301 vectorM = vector_i;
302
303 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
304 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
305}
306
307double R_MVll::computeThValue()
308{
309 double q_min = getBinMin();
310 double q_max = getBinMax();
311
312 if (vectorM == StandardModel::PHI)
313 return computeGammaPrime(q_min, q_max, lep2) / computeGammaPrime(q_min, q_max, lep1);
314 else
315 return computeGammaPrime(q_min, q_max, lep1) / computeGammaPrime(q_min, q_max, lep2);
316}
317
318
319RL_MVll::RL_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
320: F_L(SM_i, meson_i, vector_i, lep_1)
321{
322 lep1 = lep_1;
323 lep2 = lep_2;
324 meson = meson_i;
325 vectorM = vector_i;
326
327 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
328 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
329}
330
331double RL_MVll::computeThValue()
332{
333 double q_min = getBinMin();
334 double q_max = getBinMax();
335
336 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));
337}
338
339
340RT_MVll::RT_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
341: F_L(SM_i, meson_i, vector_i, lep_1)
342{
343 lep1 = lep_1;
344 lep2 = lep_2;
345 meson = meson_i;
346 vectorM = vector_i;
347
348 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
349 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
350}
351
352double RT_MVll::computeThValue()
353{
354 double q_min = getBinMin();
355 double q_max = getBinMax();
356
357 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)));
358}
359
360
361R_6::R_6(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_1, QCD::lepton lep_2)
362: ThObservable(SM_i)
363{
364 lep1 = lep_1;
365 lep2 = lep_2;
366 meson = meson_i;
367 vectorM = vector_i;
368
369 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep1).initializeMVllParameters());
370 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep2).initializeMVllParameters());
371}
372
373double R_6::computeThValue()
374{
375 double q_min = getBinMin();
376 double q_max = getBinMax();
377
378 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);
379}
380
381
382ACP_MVll::ACP_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
383: GammaPrime(SM_i, meson_i, vector_i, lep_i)
384{
385 lep = lep_i;
386 meson = meson_i;
387 vectorM = vector_i;
388
389 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
390}
391
392double ACP_MVll::computeThValue()
393{
394 double q_min = getBinMin();
395 double q_max = getBinMax();
396
397 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));
398}
399
400
401P3CP::P3CP(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
402: ThObservable(SM_i)
403{
404 lep = lep_i;
405 meson = meson_i;
406 vectorM = vector_i;
407
408 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
409}
410
411double P3CP::computeThValue()
412{
413 double q_min = getBinMin();
414 double q_max = getBinMax();
415
416 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));
417}
418
419
420F_L::F_L(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
421: GammaPrime(SM_i, meson_i, vector_i, lep_i)
422{
423 lep = lep_i;
424 meson = meson_i;
425 vectorM = vector_i;
426
427 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
428}
429
430double F_L::computeFL(double qmin, double qmax, QCD::lepton lep)
431{
432 double q_min = qmin;
433 double q_max = qmax;
434 QCD::lepton lep_i = lep;
435
436 double sigma0 = SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(0,q_min,q_max);
437 double sigma2 = SM.getFlavour().getMVll(meson, vectorM, lep_i).integrateSigma(2,q_min,q_max);
438
439 return (3.*sigma0 - sigma2) / (4. * computeGammaPrime(q_min, q_max, lep_i)) ;
440}
441
442double F_L::computeThValue()
443{
444 double q_min = getBinMin();
445 double q_max = getBinMax();
446
447 return computeFL(q_min, q_max, lep);
448}
449
450
451M_1Prime::M_1Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
452: ThObservable(SM_i)
453{
454 lep = lep_i;
455 meson = meson_i;
456 vectorM = vector_i;
457
458 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
459}
460
461double M_1Prime::computeThValue()
462{
463 double q_min = getBinMin();
464
465 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() )/
466 ( 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() ) );
467}
468
469
470M_2Prime::M_2Prime(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
471: ThObservable(SM_i)
472{
473 lep = lep_i;
474 meson = meson_i;
475 vectorM = vector_i;
476
477 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
478}
479
480double M_2Prime::computeThValue()
481{
482 double q_min = getBinMin();
483
484 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() )/
485 ( 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() );
486}
487
488
489S_1c::S_1c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
490: GammaPrime(SM_i, meson_i, vector_i, lep_i)
491{
492 lep = lep_i;
493 meson = meson_i;
494 vectorM = vector_i;
495
496 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
497}
498
499
500double S_1c::computeThValue()
501{
502 double q_min = getBinMin();
503 double q_max = getBinMax();
504
505 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(0,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
506}
507
508
509S_1s::S_1s(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
510: GammaPrime(SM_i, meson_i, vector_i, lep_i)
511{
512 lep = lep_i;
513 meson = meson_i;
514 vectorM = vector_i;
515
516 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
517}
518
519
520double S_1s::computeThValue()
521{
522 double q_min = getBinMin();
523 double q_max = getBinMax();
524
525 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(1,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
526}
527
528
529S_2c::S_2c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
530: GammaPrime(SM_i, meson_i, vector_i, lep_i)
531{
532 lep = lep_i;
533 meson = meson_i;
534 vectorM = vector_i;
535
536 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
537}
538
539
540double S_2c::computeThValue()
541{
542 double q_min = getBinMin();
543 double q_max = getBinMax();
544
545 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(2,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
546}
547
548
549S_2s::S_2s(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
550: GammaPrime(SM_i, meson_i, vector_i, lep_i)
551{
552 lep = lep_i;
553 meson = meson_i;
554 vectorM = vector_i;
555
556 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
557}
558
559
560double S_2s::computeThValue()
561{
562 double q_min = getBinMin();
563 double q_max = getBinMax();
564
565 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(3,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
566}
567
568
569S_3::S_3(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
570: GammaPrime(SM_i, meson_i, vector_i, lep_i)
571{
572 lep = lep_i;
573 meson = meson_i;
574 vectorM = vector_i;
575
576 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
577}
578
579
580double S_3::computeThValue()
581{
582 double q_min = getBinMin();
583 double q_max = getBinMax();
584
585 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(4,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
586}
587
588
589S_4::S_4(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
590: GammaPrime(SM_i, meson_i, vector_i, lep_i)
591{
592 lep = lep_i;
593 meson = meson_i;
594 vectorM = vector_i;
595
596 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
597}
598
599/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
600double S_4::computeThValue()
601{
602 double q_min = getBinMin();
603 double q_max = getBinMax();
604
605 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(5,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
606}
607
608
609S_5::S_5(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
610: GammaPrime(SM_i, meson_i, vector_i, lep_i)
611{
612 lep = lep_i;
613 meson = meson_i;
614 vectorM = vector_i;
615
616 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
617}
618
619
620double S_5::computeThValue()
621{
622 double q_min = getBinMin();
623 double q_max = getBinMax();
624
625 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(6,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
626}
627
628
629S_6c::S_6c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
630: GammaPrime(SM_i, meson_i, vector_i, lep_i)
631{
632 lep = lep_i;
633 meson = meson_i;
634 vectorM = vector_i;
635
636 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
637}
638
639
640double S_6c::computeThValue()
641{
642 double q_min = getBinMin();
643 double q_max = getBinMax();
644
645 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(8,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
646}
647
648
649S_7::S_7(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
650: GammaPrime(SM_i, meson_i, vector_i, lep_i)
651{
652 lep = lep_i;
653 meson = meson_i;
654 vectorM = vector_i;
655
656 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
657}
658
659/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
660double S_7::computeThValue()
661{
662 double q_min = getBinMin();
663 double q_max = getBinMax();
664
665 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(9,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
666}
667
668
669S_8::S_8(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
670: GammaPrime(SM_i, meson_i, vector_i, lep_i)
671{
672 lep = lep_i;
673 meson = meson_i;
674 vectorM = vector_i;
675
676 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
677}
678
679
680double S_8::computeThValue()
681{
682 double q_min = getBinMin();
683 double q_max = getBinMax();
684
685 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(10,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
686}
687
688
689S_9::S_9(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
690: GammaPrime(SM_i, meson_i, vector_i, lep_i)
691{
692 lep = lep_i;
693 meson = meson_i;
694 vectorM = vector_i;
695
696 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
697}
698
699/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
700double S_9::computeThValue()
701{
702 double q_min = getBinMin();
703 double q_max = getBinMax();
704
705 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateSigma(11,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
706}
707
708
709A_5::A_5(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
710: GammaPrime(SM_i, meson_i, vector_i, lep_i)
711{
712 lep = lep_i;
713 meson = meson_i;
714 vectorM = vector_i;
715
716 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
717}
718
719
720double A_5::computeThValue()
721{
722 double q_min = getBinMin();
723 double q_max = getBinMax();
724
725 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(6,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
726}
727
728
729A_6::A_6(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
730: GammaPrime(SM_i, meson_i, vector_i, lep_i)
731{
732 lep = lep_i;
733 meson = meson_i;
734 vectorM = vector_i;
735
736 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
737}
738
739/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
740double A_6::computeThValue()
741{
742 double q_min = getBinMin();
743 double q_max = getBinMax();
744
745 return -SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(7,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
746}
747
748
749A_6c::A_6c(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
750: GammaPrime(SM_i, meson_i, vector_i, lep_i)
751{
752 lep = lep_i;
753 meson = meson_i;
754 vectorM = vector_i;
755
756 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
757}
758
759
760double A_6c::computeThValue()
761{
762 double q_min = getBinMin();
763 double q_max = getBinMax();
764
765 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(8,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
766}
767
768
769A_8::A_8(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
770: GammaPrime(SM_i, meson_i, vector_i, lep_i)
771{
772 lep = lep_i;
773 meson = meson_i;
774 vectorM = vector_i;
775
776 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
777}
778
779
780double A_8::computeThValue()
781{
782 double q_min = getBinMin();
783 double q_max = getBinMax();
784
785 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(10,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
786}
787
788
789A_9::A_9(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
790: GammaPrime(SM_i, meson_i, vector_i, lep_i)
791{
792 lep = lep_i;
793 meson = meson_i;
794 vectorM = vector_i;
795
796 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
797}
798
799
800double A_9::computeThValue()
801{
802 double q_min = getBinMin();
803 double q_max = getBinMax();
804
805 return SM.getFlavour().getMVll(meson, vectorM, lep).integrateDelta(11,q_min,q_max) / computeGammaPrime(q_min, q_max, lep);
806}
807
808V0::V0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
809: ThObservable(SM_i)
810{
811 lep = lep_i;
812 meson = meson_i;
813 vectorM = vector_i;
814
815 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
816}
817
818double V0::computeThValue()
819{
820 return SM.getFlavour().getMVll(meson, vectorM, lep).getV0(getBinMin());
821}
822
823Vp::Vp(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
824: ThObservable(SM_i)
825{
826 lep = lep_i;
827 meson = meson_i;
828 vectorM = vector_i;
829
830 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
831}
832
833double Vp::computeThValue()
834{
835 return SM.getFlavour().getMVll(meson, vectorM, lep).getVp(getBinMin());
836}
837
838Vm::Vm(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
839: ThObservable(SM_i)
840{
841 lep = lep_i;
842 meson = meson_i;
843 vectorM = vector_i;
844
845 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
846}
847
848double Vm::computeThValue()
849{
850 return SM.getFlavour().getMVll(meson, vectorM, lep).getVm(getBinMin());
851}
852
853T0::T0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
854: ThObservable(SM_i)
855{
856 lep = lep_i;
857 meson = meson_i;
858 vectorM = vector_i;
859
860 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
861}
862
863double T0::computeThValue()
864{
865 return SM.getFlavour().getMVll(meson, vectorM, lep).getT0(getBinMin());
866}
867
868Tp::Tp(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
869: ThObservable(SM_i)
870{
871 lep = lep_i;
872 meson = meson_i;
873 vectorM = vector_i;
874
875 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
876}
877
878double Tp::computeThValue()
879{
880 return SM.getFlavour().getMVll(meson, vectorM, lep).getTp(getBinMin());
881}
882
883Tm::Tm(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
884: ThObservable(SM_i)
885{
886 lep = lep_i;
887 meson = meson_i;
888 vectorM = vector_i;
889
890 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
891}
892
893double Tm::computeThValue()
894{
895 return SM.getFlavour().getMVll(meson, vectorM, lep).getTm(getBinMin());
896}
897
898S::S(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
899: ThObservable(SM_i)
900{
901 lep = lep_i;
902 meson = meson_i;
903 vectorM = vector_i;
904
905 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
906}
907
908double S::computeThValue()
909{
910 return SM.getFlavour().getMVll(meson, vectorM, lep).getS(getBinMin());
911}
912
913unitarity_bound_f_F1::unitarity_bound_f_F1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
914: ThObservable(SM_i)
915{
916 lep = lep_i;
917 meson = meson_i;
918 vectorM = vector_i;
919
920 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
921}
922
923double unitarity_bound_f_F1::computeThValue()
924{
925 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_f_F1();
926}
927
928unitarity_bound_g::unitarity_bound_g(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
929: ThObservable(SM_i)
930{
931 lep = lep_i;
932 meson = meson_i;
933 vectorM = vector_i;
934
935 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
936}
937
938double unitarity_bound_g::computeThValue()
939{
940 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_g();
941}
942
943unitarity_bound_F2::unitarity_bound_F2(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
944: ThObservable(SM_i)
945{
946 lep = lep_i;
947 meson = meson_i;
948 vectorM = vector_i;
949
950 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
951}
952
953double unitarity_bound_F2::computeThValue()
954{
955 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_F2();
956}
957
958unitarity_bound_T1::unitarity_bound_T1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
959: ThObservable(SM_i)
960{
961 lep = lep_i;
962 meson = meson_i;
963 vectorM = vector_i;
964
965 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
966}
967
968double unitarity_bound_T1::computeThValue()
969{
970 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_T1();
971}
972
973unitarity_bound_T2_T0::unitarity_bound_T2_T0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
974: ThObservable(SM_i)
975{
976 lep = lep_i;
977 meson = meson_i;
978 vectorM = vector_i;
979
980 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
981}
982
983double unitarity_bound_T2_T0::computeThValue()
984{
985 return SM.getFlavour().getMVll(meson, vectorM, lep).get_unitarity_bound_T2_T0();
986}
987
988
989gtilde_1::gtilde_1(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
990: ThObservable(SM_i)
991{
992 lep = lep_i;
993 meson = meson_i;
994 vectorM = vector_i;
995 typ = typ_i;
996
997 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
998}
999
1000double gtilde_1::computeThValue()
1001{
1002 double q_min = getBinMin();
1003
1004 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_re(q_min));
1005 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_1_im(q_min));
1006 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());
1007 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());
1008 else throw std::runtime_error("MVllObservables::gtilde_1: incorrect type");
1009}
1010
1011
1012gtilde_2::gtilde_2(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1013: ThObservable(SM_i)
1014{
1015 lep = lep_i;
1016 meson = meson_i;
1017 vectorM = vector_i;
1018 typ = typ_i;
1019
1020 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1021}
1022
1023double gtilde_2::computeThValue()
1024{
1025 double q_min = getBinMin();
1026
1027 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_re(q_min));
1028 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_2_im(q_min));
1029 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());
1030 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());
1031 else throw std::runtime_error("MVllObservables::gtilde_2: incorrect type");
1032}
1033
1034
1035gtilde_3::gtilde_3(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1036: ThObservable(SM_i)
1037{
1038 lep = lep_i;
1039 meson = meson_i;
1040 vectorM = vector_i;
1041 typ = typ_i;
1042
1043 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1044}
1045
1046double gtilde_3::computeThValue()
1047{
1048 double q_min = getBinMin();
1049
1050 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_re(q_min));
1051 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).getgtilde_3_im(q_min));
1052 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());
1053 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());
1054 else throw std::runtime_error("MVllObservables::gtilde_3: incorrect type");
1055}
1056
1057h_0::h_0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1058: ThObservable(SM_i)
1059{
1060 lep = lep_i;
1061 meson = meson_i;
1062 vectorM = vector_i;
1063 typ = typ_i;
1064
1065 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1066}
1067
1068double h_0::computeThValue()
1069{
1070 double q_min = getBinMin();
1071
1072 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_re(q_min));
1073 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_0_im(q_min));
1074 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());
1075 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());
1076 else throw std::runtime_error("MVllObservables::h_0: type can only be 1:real, 2:imaginary, 3:absolute and 4:argument");
1077}
1078
1079
1080h_p::h_p(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1081: ThObservable(SM_i)
1082{
1083 lep = lep_i;
1084 meson = meson_i;
1085 vectorM = vector_i;
1086 typ = typ_i;
1087
1088 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1089}
1090
1091double h_p::computeThValue()
1092{
1093 double q_min = getBinMin();
1094
1095 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_re(q_min));
1096 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_p_im(q_min));
1097 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());
1098 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());
1099 else throw std::runtime_error("MVllObservables::h_p: type can only be 1:real, 2:imaginary, 3:absolute and 4:argument");
1100}
1101
1102
1103h_m::h_m(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i, unsigned int typ_i)
1104: ThObservable(SM_i)
1105{
1106 lep = lep_i;
1107 meson = meson_i;
1108 vectorM = vector_i;
1109 typ = typ_i;
1110
1111 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1112}
1113
1114double h_m::computeThValue()
1115{
1116 double q_min = getBinMin();
1117
1118 if (typ == 1) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_re(q_min));
1119 else if (typ == 2) return (SM.getFlavour().getMVll(meson, vectorM, lep).geth_m_im(q_min));
1120 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());
1121 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());
1122 else throw std::runtime_error("MVllObservables::h_m: type can only be 1:real, 2:imaginary, 3:absolute and 4:argument");
1123}
1124
1125/***********************************************************************************************************************************
1126FUNCTIONAL
1127***********************************************************************************************************************************/
1128
1129P_1f::P_1f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1130: ThObservable(SM_i)
1131{
1132 lep = lep_i;
1133 meson = meson_i;
1134 vectorM = vector_i;
1135
1136 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1137}
1138
1139double P_1f::computeThValue()
1140{
1141 double q_min = getBinMin();
1142
1143 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min)/(2.* SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1144}
1145
1146
1147P_2f::P_2f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1148: ThObservable(SM_i)
1149{
1150 lep = lep_i;
1151 meson = meson_i;
1152 vectorM = vector_i;
1153
1154 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1155}
1156
1157/*Returns experimental value, defined according to 1510.04239.*/
1158double P_2f::computeThValue()
1159{
1160 double q_min = getBinMin();
1161
1162 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));
1163
1164 return -P_2f_theory;
1165}
1166
1167
1168P_3f::P_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1169: ThObservable(SM_i)
1170{
1171 lep = lep_i;
1172 meson = meson_i;
1173 vectorM = vector_i;
1174
1175 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1176}
1177
1178/*Returns experimental value, defined according to 1510.04239.*/
1179double P_3f::computeThValue()
1180{
1181 double q_min = getBinMin();
1182
1183 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));
1184
1185 return -P_3f_theory;
1186}
1187
1188
1189P_4Primef::P_4Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1190: ThObservable(SM_i)
1191{
1192 lep = lep_i;
1193 meson = meson_i;
1194 vectorM = vector_i;
1195
1196 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1197}
1198
1199/*Returns experimental value, defined according to 1510.04239.*/
1200double P_4Primef::computeThValue()
1201{
1202 double q_min = getBinMin();
1203
1204 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));
1205
1206 return -1./2.*P_4pf_theory;
1207}
1208
1209
1210P_5Primef::P_5Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1211: ThObservable(SM_i)
1212{
1213 lep = lep_i;
1214 meson = meson_i;
1215 vectorM = vector_i;
1216
1217 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1218}
1219
1220double P_5Primef::computeThValue()
1221{
1222 double q_min = getBinMin();
1223
1224 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)));
1225}
1226
1227
1228P_6Primef::P_6Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1229: ThObservable(SM_i)
1230{
1231 lep = lep_i;
1232 meson = meson_i;
1233 vectorM = vector_i;
1234
1235 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1236}
1237
1238double P_6Primef::computeThValue()
1239{
1240 double q_min = getBinMin();
1241
1242 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)));
1243}
1244
1245/*Returns experimental value, defined according to 1510.04239.*/
1246P_8Primef::P_8Primef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1247: ThObservable(SM_i)
1248{
1249 lep = lep_i;
1250 meson = meson_i;
1251 vectorM = vector_i;
1252
1253 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1254}
1255
1256double P_8Primef::computeThValue()
1257{
1258 double q_min = getBinMin();
1259
1260 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)));
1261
1262 return -1./2.*P_8pf_theory;
1263}
1264
1265
1266GammaPrimef::GammaPrimef(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1267: ThObservable(SM_i)
1268{
1269 lep = lep_i;
1270 meson = meson_i;
1271 vectorM = vector_i;
1272
1273 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1274}
1275
1276double GammaPrimef::computeGammaPrimef(double qmin, QCD::lepton lep)
1277{
1278 double q_min = qmin;
1279 QCD::lepton lep_i = lep;
1280
1281 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.;
1282}
1283
1284double GammaPrimef::computeThValue()
1285{
1286 double q_min = getBinMin();
1287
1288 return computeGammaPrimef(q_min, lep);
1289}
1290
1291
1292A_FBf::A_FBf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1293: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1294{
1295 lep = lep_i;
1296 meson = meson_i;
1297 vectorM = vector_i;
1298
1299 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1300}
1301
1302double A_FBf::computeThValue()
1303{
1304 double q_min = getBinMin();
1305
1306 return -3. * SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(7,q_min) / 4. / computeGammaPrimef(q_min, lep);
1307}
1308
1309F_Lf::F_Lf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1310: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1311{
1312 lep = lep_i;
1313 meson = meson_i;
1314 vectorM = vector_i;
1315
1316 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1317}
1318
1319double F_Lf::computeFLf(double qmin, QCD::lepton lep)
1320{
1321 double q_min = qmin;
1322 QCD::lepton lep_i = lep;
1323
1324 double sigma0 = SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(0,q_min);
1325 double sigma2 = SM.getFlavour().getMVll(meson, vectorM, lep_i).getSigma(2,q_min);
1326
1327 return (3.*sigma0 - sigma2) / (4. * computeGammaPrimef(q_min, lep_i)) ;
1328}
1329
1330
1331double F_Lf::computeThValue()
1332{
1333 double q_min = getBinMin();
1334
1335 return computeFLf(q_min, lep);
1336}
1337
1338S_3f::S_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1339: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1340{
1341 lep = lep_i;
1342 meson = meson_i;
1343 vectorM = vector_i;
1344
1345 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1346}
1347
1348
1349double S_3f::computeThValue()
1350{
1351 double q_min = getBinMin();
1352
1353 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min) / computeGammaPrimef(q_min, lep);
1354}
1355
1356
1357S_4f::S_4f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1358: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1359{
1360 lep = lep_i;
1361 meson = meson_i;
1362 vectorM = vector_i;
1363
1364 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1365}
1366
1367/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
1368double S_4f::computeThValue()
1369{
1370 double q_min = getBinMin();
1371
1372 return -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(5,q_min) / computeGammaPrimef(q_min, lep);
1373}
1374
1375
1376S_5f::S_5f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1377: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1378{
1379 lep = lep_i;
1380 meson = meson_i;
1381 vectorM = vector_i;
1382
1383 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1384}
1385
1386
1387double S_5f::computeThValue()
1388{
1389 double q_min = getBinMin();
1390
1391 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(6,q_min) / computeGammaPrimef(q_min, lep);
1392}
1393
1394
1395S_7f::S_7f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1396: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1397{
1398 lep = lep_i;
1399 meson = meson_i;
1400 vectorM = vector_i;
1401
1402 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1403}
1404
1405/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
1406double S_7f::computeThValue()
1407{
1408 double q_min = getBinMin();
1409
1410 return -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(9,q_min) / computeGammaPrimef(q_min, lep);
1411}
1412
1413
1414S_8f::S_8f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1415: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1416{
1417 lep = lep_i;
1418 meson = meson_i;
1419 vectorM = vector_i;
1420
1421 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1422}
1423
1424double S_8f::computeThValue()
1425{
1426 double q_min = getBinMin();
1427
1428 return SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(10,q_min) / computeGammaPrimef(q_min, lep);
1429}
1430
1431
1432S_9f::S_9f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1433: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1434{
1435 lep = lep_i;
1436 meson = meson_i;
1437 vectorM = vector_i;
1438
1439 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1440}
1441
1442/*Returns experimental value, due to difference in \theta_l and \phi definitions*/
1443double S_9f::computeThValue()
1444{
1445 double q_min = getBinMin();
1446
1447 return -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(11,q_min) / computeGammaPrimef(q_min, lep);
1448}
1449
1450BRf_MVll::BRf_MVll(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1451: GammaPrimef(SM_i, meson_i, vector_i, lep_i)
1452{
1453 lep = lep_i;
1454 meson = meson_i;
1455 vectorM = vector_i;
1456
1457 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1458}
1459
1460double BRf_MVll::computeThValue()
1461{
1462 double q_min = getBinMin();
1463
1464 return computeGammaPrimef(q_min, lep)/SM.getFlavour().getMVll(meson, vectorM, lep).getwidth();
1465}
1466
1467P_relationf::P_relationf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1468: ThObservable(SM_i)
1469{
1470 lep = lep_i;
1471 meson = meson_i;
1472 vectorM = vector_i;
1473
1474 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1475}
1476
1477double P_relationf::computeThValue()
1478{
1479 double q_min = getBinMin();
1480 double P1 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min)/(2.* SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1481 double P2 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(7,q_min)/(8.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1482 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));
1483 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)));
1484 double beta = sqrt(1. - 4.*SM.getLeptons(lep).getMass()*SM.getLeptons(lep).getMass()/q_min/q_min);
1485
1486 return 1./2.*(P4p*P5p + 1./beta *sqrt(std::abs((-1. + P1 + P4p*P4p)*(-1. - P1 + beta*beta*P5p*P5p)))) - P2;
1487}
1488
1489P_relation_exactf::P_relation_exactf(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1490: ThObservable(SM_i)
1491{
1492 lep = lep_i;
1493 meson = meson_i;
1494 vectorM = vector_i;
1495
1496 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1497}
1498
1499double P_relation_exactf::computeThValue()
1500{
1501 double q_min = getBinMin();
1502 double P1 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(4,q_min)/(2.* SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1503 double P2 = SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(7,q_min)/(8.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1504 double P3 = -SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(11,q_min)/(4.*SM.getFlavour().getMVll(meson, vectorM, lep).getSigma(3,q_min));
1505 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));
1506 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)));
1507 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)));
1508 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)));
1509
1510 double beta = sqrt(1. - 4.*SM.getLeptons(lep).getMass()*SM.getLeptons(lep).getMass()/q_min/q_min);
1511
1512 /* k1 and k2 not fully implemented. arXiv:1402.6855*/
1513 double k1 = 1.;
1514 double k2 = 1.;
1515
1516 double delta_1 = P6p*P8p;
1517 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));
1518 double delta_3 = (1. - k1)*P4p*P4p + beta*beta*((-1. + k1)*P5p*P5p - k1*P6p*P6p) + k1*P8p*P8p;
1519 double delta_4 = 1. - k1*k1;
1520
1521 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;
1522}
1523
1524QCDfC9_1f::QCDfC9_1f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1525: ThObservable(SM_i)
1526{
1527 lep = lep_i;
1528 meson = meson_i;
1529 vectorM = vector_i;
1530
1531 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1532}
1533
1534double QCDfC9_1f::computeThValue()
1535{
1536 double q2 = getBinMin();
1537 double cutoff = getBinMax();
1538
1539 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9_1(q2, cutoff);
1540}
1541
1542QCDfC9_2f::QCDfC9_2f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1543: ThObservable(SM_i)
1544{
1545 lep = lep_i;
1546 meson = meson_i;
1547 vectorM = vector_i;
1548
1549 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1550}
1551
1552double QCDfC9_2f::computeThValue()
1553{
1554 double q2 = getBinMin();
1555 double cutoff = getBinMax();
1556
1557 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9_2(q2, cutoff);
1558}
1559
1560QCDfC9_3f::QCDfC9_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1561: ThObservable(SM_i)
1562{
1563 lep = lep_i;
1564 meson = meson_i;
1565 vectorM = vector_i;
1566
1567 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1568}
1569
1570double QCDfC9_3f::computeThValue()
1571{
1572 double q2 = getBinMin();
1573 double cutoff = getBinMax();
1574
1575 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9_3(q2, cutoff);
1576}
1577
1578QCDfC9p_1f::QCDfC9p_1f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1579: ThObservable(SM_i)
1580{
1581 lep = lep_i;
1582 meson = meson_i;
1583 vectorM = vector_i;
1584
1585 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1586}
1587
1588double QCDfC9p_1f::computeThValue()
1589{
1590 double cutoff = getBinMin();
1591
1592 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9p_1(cutoff);
1593}
1594
1595QCDfC9p_2f::QCDfC9p_2f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1596: ThObservable(SM_i)
1597{
1598 lep = lep_i;
1599 meson = meson_i;
1600 vectorM = vector_i;
1601
1602 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1603}
1604
1605double QCDfC9p_2f::computeThValue()
1606{
1607 double cutoff = getBinMin();
1608
1609 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9p_2(cutoff);
1610}
1611
1612QCDfC9p_3f::QCDfC9p_3f(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1613: ThObservable(SM_i)
1614{
1615 lep = lep_i;
1616 meson = meson_i;
1617 vectorM = vector_i;
1618
1619 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1620}
1621
1622double QCDfC9p_3f::computeThValue()
1623{
1624 double cutoff = getBinMin();
1625
1626 return SM.getFlavour().getMVll(meson, vectorM, lep).getQCDfC9p_3(cutoff);
1627}
1628
1629BR_MVpsi::BR_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1630: ThObservable(SM_i)
1631{
1632 lep = lep_i;
1633 meson = meson_i;
1634 vectorM = vector_i;
1635
1636 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1637}
1638
1639double BR_MVpsi::computeThValue()
1640{
1641 double q2 = getBinMin();
1642 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1643 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1644 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1645 return 2.*(A2_0+A2_par+A2_perp)/SM.getFlavour().getMVll(meson, vectorM, lep).getwidth();
1646}
1647
1648BR_MVpsi_ratio::BR_MVpsi_ratio(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1649: ThObservable(SM_i)
1650{
1651 lep = lep_i;
1652 meson = meson_i;
1653 vectorM = vector_i;
1654
1655 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1656}
1657
1658double BR_MVpsi_ratio::computeThValue()
1659{
1660 double A2_0_JPSI = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.0969,0).abs2();
1661 double A2_par_JPSI = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.0969,1).abs2();
1662 double A2_perp_JPSI = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.0969,2).abs2();
1663
1664 double A2_0_PSI2S = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.6861,0).abs2();
1665 double A2_par_PSI2S = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.6861,1).abs2();
1666 double A2_perp_PSI2S = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(3.6861,2).abs2();
1667
1668 return (A2_0_PSI2S+A2_par_PSI2S+A2_perp_PSI2S)/(A2_0_JPSI+A2_par_JPSI+A2_perp_JPSI);
1669}
1670
1671Abs2Ampar_MVpsi::Abs2Ampar_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1672: ThObservable(SM_i)
1673{
1674 lep = lep_i;
1675 meson = meson_i;
1676 vectorM = vector_i;
1677
1678 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1679}
1680
1681double Abs2Ampar_MVpsi::computeThValue()
1682{
1683 double q2 = getBinMin();
1684 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1685 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1686 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1687 return A2_par/(A2_0+A2_par+A2_perp);
1688}
1689
1690ArgAmpar_MVpsi::ArgAmpar_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1691: ThObservable(SM_i)
1692{
1693 lep = lep_i;
1694 meson = meson_i;
1695 vectorM = vector_i;
1696
1697 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1698}
1699
1700double ArgAmpar_MVpsi::computeThValue()
1701{
1702 double q2 = getBinMin();
1703 return (SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1)/SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0)).arg();
1704}
1705
1706Abs2Amperp_MVpsi::Abs2Amperp_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1707: ThObservable(SM_i)
1708{
1709 lep = lep_i;
1710 meson = meson_i;
1711 vectorM = vector_i;
1712
1713 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1714}
1715
1716double Abs2Amperp_MVpsi::computeThValue()
1717{
1718 double q2 = getBinMin();
1719 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1720 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1721 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1722 return A2_perp/(A2_0+A2_par+A2_perp);
1723}
1724
1725ArgAmperp_MVpsi::ArgAmperp_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1726: ThObservable(SM_i)
1727{
1728 lep = lep_i;
1729 meson = meson_i;
1730 vectorM = vector_i;
1731
1732 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1733}
1734
1735double ArgAmperp_MVpsi::computeThValue()
1736{
1737 double q2 = getBinMin();
1738 return (SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2)/SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0)).arg();
1739}
1740
1741Abs2Ampzero_MVpsi::Abs2Ampzero_MVpsi(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1742: ThObservable(SM_i)
1743{
1744 lep = lep_i;
1745 meson = meson_i;
1746 vectorM = vector_i;
1747
1748 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1749}
1750
1751double Abs2Ampzero_MVpsi::computeThValue()
1752{
1753 double q2 = getBinMin();
1754 double A2_0 = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,0).abs2();
1755 double A2_par = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,1).abs2();
1756 double A2_perp = SM.getFlavour().getMVll(meson, vectorM, lep).AmpMVpsi_zExpansion(q2,2).abs2();
1757 return A2_0/(A2_0+A2_par+A2_perp);
1758}
1759
1760Delta_C9_zExp_0::Delta_C9_zExp_0(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1761: ThObservable(SM_i)
1762{
1763 lep = lep_i;
1764 meson = meson_i;
1765 vectorM = vector_i;
1766
1767 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1768}
1769
1770double Delta_C9_zExp_0::computeThValue()
1771{
1772 return SM.getFlavour().getMVll(meson, vectorM, lep).getDelta_C9_zExp_0();
1773}
1774
1775Delta_C9_zExp_p::Delta_C9_zExp_p(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1776: ThObservable(SM_i)
1777{
1778 lep = lep_i;
1779 meson = meson_i;
1780 vectorM = vector_i;
1781
1782 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1783}
1784
1785double Delta_C9_zExp_p::computeThValue()
1786{
1787 return SM.getFlavour().getMVll(meson, vectorM, lep).getDelta_C9_zExp_p();
1788}
1789
1790Delta_C9_zExp_m::Delta_C9_zExp_m(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
1791: ThObservable(SM_i)
1792{
1793 lep = lep_i;
1794 meson = meson_i;
1795 vectorM = vector_i;
1796
1797 setParametersForObservable(SM.getFlavour().getMVll(meson, vectorM, lep).initializeMVllParameters());
1798}
1799
1800double Delta_C9_zExp_m::computeThValue()
1801{
1802 return SM.getFlavour().getMVll(meson, vectorM, lep).getDelta_C9_zExp_m();
1803}
A_5::meson
QCD::meson meson
Definition: MVllObservables.h:1328
A_5::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1329
A_5::lep
QCD::lepton lep
Definition: MVllObservables.h:1327
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:709
A_5::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:720
A_6::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:740
A_6::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1368
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:729
A_6::lep
QCD::lepton lep
Definition: MVllObservables.h:1366
A_6::meson
QCD::meson meson
Definition: MVllObservables.h:1367
A_6c::lep
QCD::lepton lep
Definition: MVllObservables.h:1405
A_6c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1407
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:749
A_6c::meson
QCD::meson meson
Definition: MVllObservables.h:1406
A_6c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:760
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:769
A_8::meson
QCD::meson meson
Definition: MVllObservables.h:1445
A_8::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1446
A_8::lep
QCD::lepton lep
Definition: MVllObservables.h:1444
A_8::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:780
A_9::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1485
A_9::meson
QCD::meson meson
Definition: MVllObservables.h:1484
A_9::lep
QCD::lepton lep
Definition: MVllObservables.h:1483
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:789
A_9::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:800
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:2496
A_FBf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1302
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:1292
A_FBf::lep
QCD::lepton lep
Definition: MVllObservables.h:2495
A_FBf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2497
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:743
ACP_MVll::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:392
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:382
ACP_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:742
ACP_MVll::lep
QCD::lepton lep
Definition: MVllObservables.h:741
Abs2Ampar_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3171
Abs2Ampar_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3172
Abs2Ampar_MVpsi::computeThValue
double computeThValue()
The observable Abs(A//) of .
Definition: MVllObservables.cpp:1681
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:1671
Abs2Ampar_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3170
Abs2Amperp_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3239
Abs2Amperp_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3238
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:1706
Abs2Amperp_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3240
Abs2Amperp_MVpsi::computeThValue
double computeThValue()
The observable Abs(A_|_) of .
Definition: MVllObservables.cpp:1716
Abs2Ampzero_MVpsi::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1751
Abs2Ampzero_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3306
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:1741
Abs2Ampzero_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3307
Abs2Ampzero_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3308
ArgAmpar_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3206
ArgAmpar_MVpsi::computeThValue
double computeThValue()
The observable Arg(A//) of .
Definition: MVllObservables.cpp:1700
ArgAmpar_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3205
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:1690
ArgAmpar_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3204
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:1725
ArgAmperp_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3274
ArgAmperp_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3273
ArgAmperp_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3272
ArgAmperp_MVpsi::computeThValue
double computeThValue()
The observable Arg(A_|_) of .
Definition: MVllObservables.cpp:1735
BR_MVll_tot::lep
QCD::lepton lep
Definition: MVllObservables.h:491
BR_MVll_tot::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:272
BR_MVll_tot::BR_MVll_tot
BR_MVll_tot(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:262
BR_MVll_tot::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:493
BR_MVll_tot::meson
QCD::meson meson
Definition: MVllObservables.h:492
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:3137
BR_MVpsi_ratio::lep
QCD::lepton lep
Definition: MVllObservables.h:3136
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:1648
BR_MVpsi_ratio::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3138
BR_MVpsi_ratio::computeThValue
double computeThValue()
The observable BR of .
Definition: MVllObservables.cpp:1658
BR_MVpsi::meson
QCD::meson meson
Definition: MVllObservables.h:3103
BR_MVpsi::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3104
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:1629
BR_MVpsi::computeThValue
double computeThValue()
The observable BR of .
Definition: MVllObservables.cpp:1639
BR_MVpsi::lep
QCD::lepton lep
Definition: MVllObservables.h:3102
BRf_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2536
BRf_MVll::lep
QCD::lepton lep
Definition: MVllObservables.h:2534
BRf_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:2535
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:1450
BRf_MVll::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1460
Delta_C9_zExp_0::lep
QCD::lepton lep
Definition: MVllObservables.h:3340
Delta_C9_zExp_0::meson
QCD::meson meson
Definition: MVllObservables.h:3341
Delta_C9_zExp_0::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1770
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:1760
Delta_C9_zExp_0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3342
Delta_C9_zExp_m::lep
QCD::lepton lep
Definition: MVllObservables.h:3408
Delta_C9_zExp_m::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3410
Delta_C9_zExp_m::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1800
Delta_C9_zExp_m::meson
QCD::meson meson
Definition: MVllObservables.h:3409
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:1790
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:1775
Delta_C9_zExp_p::lep
QCD::lepton lep
Definition: MVllObservables.h:3374
Delta_C9_zExp_p::meson
QCD::meson meson
Definition: MVllObservables.h:3375
Delta_C9_zExp_p::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3376
Delta_C9_zExp_p::computeThValue
double computeThValue()
The observable Abs(A0) of .
Definition: MVllObservables.cpp:1785
F_L
A class for the binned observable in .
Definition: MVllObservables.h:511
F_L::lep
QCD::lepton lep
Definition: MVllObservables.h:539
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:430
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:420
F_L::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:442
F_L::meson
QCD::meson meson
Definition: MVllObservables.h:540
F_L::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:541
F_Lf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2583
F_Lf::lep
QCD::lepton lep
Definition: MVllObservables.h:2581
F_Lf::meson
QCD::meson meson
Definition: MVllObservables.h:2582
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:1309
F_Lf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1331
F_Lf::computeFLf
double computeFLf(double qmin, QCD::lepton lep)
A method to compute the observable in at a generic .
Definition: MVllObservables.cpp:1319
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:2429
GammaPrimef::meson
QCD::meson meson
Definition: MVllObservables.h:2457
GammaPrimef::computeGammaPrimef
double computeGammaPrimef(double qmin, QCD::lepton lep)
A method to compute the observable in at a generic .
Definition: MVllObservables.cpp:1276
GammaPrimef::lep
QCD::lepton lep
Definition: MVllObservables.h:2456
GammaPrimef::GammaPrimef
GammaPrimef(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:1266
GammaPrimef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1284
GammaPrimef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2458
M_1Prime::meson
QCD::meson meson
Definition: MVllObservables.h:821
M_1Prime::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:461
M_1Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:822
M_1Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:820
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:451
M_2Prime::meson
QCD::meson meson
Definition: MVllObservables.h:860
M_2Prime::lep
QCD::lepton lep
Definition: MVllObservables.h:859
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:470
M_2Prime::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:861
M_2Prime::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:480
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:2590
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:2602
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:3040
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:2566
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:2992
MVll::H_V_p
gslpp::complex H_V_p(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2572
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:3139
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:2578
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:2660
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:2584
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:2596
MVll::H_P
gslpp::complex H_P(double q2, bool bar)
The helicity amplitude .
Definition: MVll.cpp:2610
MVll::integrateSigma
double integrateSigma(int i, double q_min, double q_max)
The integral of from to .
Definition: MVll.cpp:2865
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:2692
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:781
P3CP::P3CP
P3CP(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:401
P3CP::computeThValue
double computeThValue()
The binned observable in .
Definition: MVllObservables.cpp:411
P3CP::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:783
P3CP::meson
QCD::meson meson
Definition: MVllObservables.h:782
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:2172
P_1f::lep
QCD::lepton lep
Definition: MVllObservables.h:2170
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:1129
P_1f::meson
QCD::meson meson
Definition: MVllObservables.h:2171
P_1f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1139
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:1147
P_2f::meson
QCD::meson meson
Definition: MVllObservables.h:2210
P_2f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2211
P_2f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1158
P_2f::lep
QCD::lepton lep
Definition: MVllObservables.h:2209
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:2249
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:1168
P_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2250
P_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:2248
P_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1179
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:1200
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:1189
P_4Primef::lep
QCD::lepton lep
Definition: MVllObservables.h:2290
P_4Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2292
P_4Primef::meson
QCD::meson meson
Definition: MVllObservables.h:2291
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:2332
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:1210
P_5Primef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1220
P_5Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2333
P_5Primef::lep
QCD::lepton lep
Definition: MVllObservables.h:2331
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:2370
P_6Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2372
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:1228
P_6Primef::meson
QCD::meson meson
Definition: MVllObservables.h:2371
P_6Primef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1238
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:1246
P_8Primef::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1256
P_8Primef::meson
QCD::meson meson
Definition: MVllObservables.h:2410
P_8Primef::lep
QCD::lepton lep
Definition: MVllObservables.h:2409
P_8Primef::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2411
P_relation_exactf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1499
P_relation_exactf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2866
P_relation_exactf::meson
QCD::meson meson
Definition: MVllObservables.h:2865
P_relation_exactf::lep
QCD::lepton lep
Definition: MVllObservables.h:2864
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:1489
P_relationf::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2841
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:1467
P_relationf::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1477
P_relationf::meson
QCD::meson meson
Definition: MVllObservables.h:2840
P_relationf::lep
QCD::lepton lep
Definition: MVllObservables.h:2839
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:2900
QCDfC9_1f::lep
QCD::lepton lep
Definition: MVllObservables.h:2898
QCDfC9_1f::meson
QCD::meson meson
Definition: MVllObservables.h:2899
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:1524
QCDfC9_1f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1534
QCDfC9_2f::lep
QCD::lepton lep
Definition: MVllObservables.h:2932
QCDfC9_2f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2934
QCDfC9_2f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1552
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:1542
QCDfC9_2f::meson
QCD::meson meson
Definition: MVllObservables.h:2933
QCDfC9_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1570
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:1560
QCDfC9_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:2966
QCDfC9_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2968
QCDfC9_3f::meson
QCD::meson meson
Definition: MVllObservables.h:2967
QCDfC9p_1f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3002
QCDfC9p_1f::lep
QCD::lepton lep
Definition: MVllObservables.h:3000
QCDfC9p_1f::meson
QCD::meson meson
Definition: MVllObservables.h:3001
QCDfC9p_1f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1588
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:1578
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:1595
QCDfC9p_2f::meson
QCD::meson meson
Definition: MVllObservables.h:3035
QCDfC9p_2f::lep
QCD::lepton lep
Definition: MVllObservables.h:3034
QCDfC9p_2f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1605
QCDfC9p_2f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3036
QCDfC9p_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:3068
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:1612
QCDfC9p_3f::meson
QCD::meson meson
Definition: MVllObservables.h:3069
QCDfC9p_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:3070
QCDfC9p_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1622
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:361
R_6::meson
QCD::meson meson
Definition: MVllObservables.h:704
R_6::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:705
R_6::lep1
QCD::lepton lep1
Definition: MVllObservables.h:702
R_6::lep2
QCD::lepton lep2
Definition: MVllObservables.h:703
R_6::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:373
R_MVll::lep1
QCD::lepton lep1
Definition: MVllObservables.h:579
R_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:581
R_MVll::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:307
R_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:582
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:295
R_MVll::lep2
QCD::lepton lep2
Definition: MVllObservables.h:580
RL_MVll::lep2
QCD::lepton lep2
Definition: MVllObservables.h:621
RL_MVll::lep1
QCD::lepton lep1
Definition: MVllObservables.h:620
RL_MVll::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:331
RL_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:623
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:319
RL_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:622
RT_MVll::lep1
QCD::lepton lep1
Definition: MVllObservables.h:661
RT_MVll::lep2
QCD::lepton lep2
Definition: MVllObservables.h:662
RT_MVll::computeThValue
double computeThValue()
The binned observable ratio in .
Definition: MVllObservables.cpp:352
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:340
RT_MVll::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:664
RT_MVll::meson
QCD::meson meson
Definition: MVllObservables.h:663
S_1c::meson
QCD::meson meson
Definition: MVllObservables.h:899
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:489
S_1c::lep
QCD::lepton lep
Definition: MVllObservables.h:898
S_1c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:900
S_1c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:500
S_1s::lep
QCD::lepton lep
Definition: MVllObservables.h:937
S_1s::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:939
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:509
S_1s::meson
QCD::meson meson
Definition: MVllObservables.h:938
S_1s::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:520
S_2c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:978
S_2c::lep
QCD::lepton lep
Definition: MVllObservables.h:976
S_2c::meson
QCD::meson meson
Definition: MVllObservables.h:977
S_2c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:540
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:529
S_2s::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1017
S_2s::meson
QCD::meson meson
Definition: MVllObservables.h:1016
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:549
S_2s::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:560
S_2s::lep
QCD::lepton lep
Definition: MVllObservables.h:1015
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:569
S_3::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1056
S_3::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:580
S_3::meson
QCD::meson meson
Definition: MVllObservables.h:1055
S_3::lep
QCD::lepton lep
Definition: MVllObservables.h:1054
S_3f::meson
QCD::meson meson
Definition: MVllObservables.h:2620
S_3f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1349
S_3f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2621
S_3f::lep
QCD::lepton lep
Definition: MVllObservables.h:2619
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:1338
S_4::lep
QCD::lepton lep
Definition: MVllObservables.h:1093
S_4::meson
QCD::meson meson
Definition: MVllObservables.h:1094
S_4::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1095
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:589
S_4::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:600
S_4f::lep
QCD::lepton lep
Definition: MVllObservables.h:2658
S_4f::meson
QCD::meson meson
Definition: MVllObservables.h:2659
S_4f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2660
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:1357
S_4f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1368
S_5::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:620
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:609
S_5::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1134
S_5::lep
QCD::lepton lep
Definition: MVllObservables.h:1132
S_5::meson
QCD::meson meson
Definition: MVllObservables.h:1133
S_5f::lep
QCD::lepton lep
Definition: MVllObservables.h:2697
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:1376
S_5f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1387
S_5f::meson
QCD::meson meson
Definition: MVllObservables.h:2698
S_5f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2699
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:629
S_6c::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1173
S_6c::meson
QCD::meson meson
Definition: MVllObservables.h:1172
S_6c::lep
QCD::lepton lep
Definition: MVllObservables.h:1171
S_6c::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:640
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:649
S_7::lep
QCD::lepton lep
Definition: MVllObservables.h:1210
S_7::meson
QCD::meson meson
Definition: MVllObservables.h:1211
S_7::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1212
S_7::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:660
S_7f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2738
S_7f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1406
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:1395
S_7f::meson
QCD::meson meson
Definition: MVllObservables.h:2737
S_7f::lep
QCD::lepton lep
Definition: MVllObservables.h:2736
S_8::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1251
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:669
S_8::meson
QCD::meson meson
Definition: MVllObservables.h:1250
S_8::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:680
S_8::lep
QCD::lepton lep
Definition: MVllObservables.h:1249
S_8f::meson
QCD::meson meson
Definition: MVllObservables.h:2776
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:1414
S_8f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2777
S_8f::lep
QCD::lepton lep
Definition: MVllObservables.h:2775
S_8f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1424
S_9::lep
QCD::lepton lep
Definition: MVllObservables.h:1288
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:689
S_9::meson
QCD::meson meson
Definition: MVllObservables.h:1289
S_9::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1290
S_9::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:700
S_9f::computeThValue
double computeThValue()
The observable in .
Definition: MVllObservables.cpp:1443
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:1432
S_9f::lep
QCD::lepton lep
Definition: MVllObservables.h:2814
S_9f::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2816
S_9f::meson
QCD::meson meson
Definition: MVllObservables.h:2815
S::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1730
S::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:908
S::S
S(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:898
S::meson
QCD::meson meson
Definition: MVllObservables.h:1729
S::lep
QCD::lepton lep
Definition: MVllObservables.h:1728
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:1623
T0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1625
T0::T0
T0(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:853
T0::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:863
T0::meson
QCD::meson meson
Definition: MVllObservables.h:1624
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:1693
Tm::meson
QCD::meson meson
Definition: MVllObservables.h:1694
Tm::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:893
Tm::Tm
Tm(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:883
Tm::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1695
Tp::meson
QCD::meson meson
Definition: MVllObservables.h:1659
Tp::Tp
Tp(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:868
Tp::lep
QCD::lepton lep
Definition: MVllObservables.h:1658
Tp::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1660
Tp::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:878
V0::V0
V0(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:808
V0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1520
V0::lep
QCD::lepton lep
Definition: MVllObservables.h:1518
V0::meson
QCD::meson meson
Definition: MVllObservables.h:1519
V0::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:818
Vm::Vm
Vm(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:838
Vm::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1590
Vm::meson
QCD::meson meson
Definition: MVllObservables.h:1589
Vm::lep
QCD::lepton lep
Definition: MVllObservables.h:1588
Vm::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:848
Vp::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1555
Vp::Vp
Vp(const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
Constructor.
Definition: MVllObservables.cpp:823
Vp::meson
QCD::meson meson
Definition: MVllObservables.h:1554
Vp::computeThValue
double computeThValue()
The averaged form factor in .
Definition: MVllObservables.cpp:833
Vp::lep
QCD::lepton lep
Definition: MVllObservables.h:1553
gtilde_1::typ
unsigned int typ
Definition: MVllObservables.h:1937
gtilde_1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1936
gtilde_1::computeThValue
double computeThValue()
The averaged correction in .
Definition: MVllObservables.cpp:1000
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:989
gtilde_1::meson
QCD::meson meson
Definition: MVllObservables.h:1935
gtilde_1::lep
QCD::lepton lep
Definition: MVllObservables.h:1934
gtilde_2::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1975
gtilde_2::lep
QCD::lepton lep
Definition: MVllObservables.h:1973
gtilde_2::typ
unsigned int typ
Definition: MVllObservables.h:1976
gtilde_2::meson
QCD::meson meson
Definition: MVllObservables.h:1974
gtilde_2::computeThValue
double computeThValue()
The averaged correction in .
Definition: MVllObservables.cpp:1023
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:1012
gtilde_3::computeThValue
double computeThValue()
The averaged correction in .
Definition: MVllObservables.cpp:1046
gtilde_3::meson
QCD::meson meson
Definition: MVllObservables.h:2013
gtilde_3::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2014
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:1035
gtilde_3::typ
unsigned int typ
Definition: MVllObservables.h:2015
gtilde_3::lep
QCD::lepton lep
Definition: MVllObservables.h:2012
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:1057
h_0::meson
QCD::meson meson
Definition: MVllObservables.h:2050
h_0::typ
unsigned int typ
Definition: MVllObservables.h:2052
h_0::lep
QCD::lepton lep
Definition: MVllObservables.h:2049
h_0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2051
h_0::computeThValue
double computeThValue()
The correction in .
Definition: MVllObservables.cpp:1068
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:1103
h_m::typ
unsigned int typ
Definition: MVllObservables.h:2128
h_m::computeThValue
double computeThValue()
The correction in .
Definition: MVllObservables.cpp:1114
h_m::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2127
h_m::lep
QCD::lepton lep
Definition: MVllObservables.h:2125
h_m::meson
QCD::meson meson
Definition: MVllObservables.h:2126
h_p::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:2089
h_p::typ
unsigned int typ
Definition: MVllObservables.h:2090
h_p::meson
QCD::meson meson
Definition: MVllObservables.h:2088
h_p::lep
QCD::lepton lep
Definition: MVllObservables.h:2087
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:1080
h_p::computeThValue
double computeThValue()
The correction in .
Definition: MVllObservables.cpp:1091
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:943
unitarity_bound_F2::computeThValue
double computeThValue()
The unitarity bound on form factors in decays.
Definition: MVllObservables.cpp:953
unitarity_bound_F2::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1832
unitarity_bound_F2::meson
QCD::meson meson
Definition: MVllObservables.h:1831
unitarity_bound_F2::lep
QCD::lepton lep
Definition: MVllObservables.h:1830
unitarity_bound_T1::computeThValue
double computeThValue()
The unitarity bound on form factors in decays.
Definition: MVllObservables.cpp:968
unitarity_bound_T1::meson
QCD::meson meson
Definition: MVllObservables.h:1865
unitarity_bound_T1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1866
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:958
unitarity_bound_T1::lep
QCD::lepton lep
Definition: MVllObservables.h:1864
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:973
unitarity_bound_T2_T0::computeThValue
double computeThValue()
The unitarity bound on form factors and in decays.
Definition: MVllObservables.cpp:983
unitarity_bound_T2_T0::lep
QCD::lepton lep
Definition: MVllObservables.h:1898
unitarity_bound_T2_T0::meson
QCD::meson meson
Definition: MVllObservables.h:1899
unitarity_bound_T2_T0::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1900
unitarity_bound_f_F1::lep
QCD::lepton lep
Definition: MVllObservables.h:1762
unitarity_bound_f_F1::computeThValue
double computeThValue()
The unitarity bound on form factors and in decays.
Definition: MVllObservables.cpp:923
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:913
unitarity_bound_f_F1::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1764
unitarity_bound_f_F1::meson
QCD::meson meson
Definition: MVllObservables.h:1763
unitarity_bound_g::vectorM
QCD::meson vectorM
Definition: MVllObservables.h:1798
unitarity_bound_g::lep
QCD::lepton lep
Definition: MVllObservables.h:1796
unitarity_bound_g::computeThValue
double computeThValue()
The unitarity bound on form factors in decays.
Definition: MVllObservables.cpp:938
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:928
unitarity_bound_g::meson
QCD::meson meson
Definition: MVllObservables.h:1797