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