-
Notifications
You must be signed in to change notification settings - Fork 60
/
test_lambdab.py
58 lines (43 loc) · 2.84 KB
/
test_lambdab.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
import unittest
import flavio
from flavio.classes import Implementation
from flavio.physics.bdecays.test_lambdablambdall import ass_sm
wc_sm = flavio.WilsonCoefficients()
def pred_sm(s, name, q2val, target, delta, scalef=1):
obs = flavio.classes.Observable[name]
c = flavio.sm_prediction(name, q2=q2val)*scalef
s.assertAlmostEqual(c, target, delta=delta)
class TestLambdabLambda1520_FF(unittest.TestCase):
def test_lambdab_lambda1520_ff(self):
flavio.config['implementation']['Lambdab->Lambda(1520) form factor'] = 'Lambdab->Lambda(1520) LatticeQCD'
par = flavio.default_parameters.get_central_all()
mLst = par['m_Lambda(1520)']
mLb = par['m_Lambdab']
BR = par['BR(Lambda(1520)->NKbar)_exp']
BRinv = 2/BR
# Comparison to figure on slide 12 S.Meinel b-baryon FEST 2020
pred_sm(self, 'dBR/dq2(Lambdab->Lambda(1520)mumu)', 16.0, pow(2.5, -9), pow(0.1, -9), BRinv)
pred_sm(self, 'dBR/dq2(Lambdab->Lambda(1520)mumu)', 16.4, pow(1.2, -9), pow(0.1, -9), BRinv)
pred_sm(self, 'dBR/dq2(Lambdab->Lambda(1520)mumu)', 16.8, 0.0, pow(0.01, -9), BRinv)
# Comparison to figure on slide 14 S.Meinel b-baryon FEST 2020
pred_sm(self, 'AFBl(Lambdab->Lambda(1520)mumu)', 16.0, -0.09, 0.01)
pred_sm(self, 'AFBl(Lambdab->Lambda(1520)mumu)', 16.4, 0.05, 0.01)
pred_sm(self, 'AFBl(Lambdab->Lambda(1520)mumu)', 16.8, 0.25, 0.02)
# Comparison to figure on slide 13 S.Meinel b-baryon FEST 2020
pred_sm(self, 'S_1cc(Lambdab->Lambda(1520)mumu)', 16.0, 0.56, 0.01)
pred_sm(self, 'S_1cc(Lambdab->Lambda(1520)mumu)', 16.4, 0.56, 0.02)
pred_sm(self, 'S_1cc(Lambdab->Lambda(1520)mumu)', 16.8, 0.36, 0.02)
flavio.config['implementation']['Lambdab->Lambda(1520) form factor'] = 'Lambdab->Lambda(1520) MCN'
# Comparison to figure on slide 12 S.Meinel b-baryon FEST 2020
pred_sm(self, 'dBR/dq2(Lambdab->Lambda(1520)mumu)', 16.0, pow(5.4, -9), pow(0.2, -9), BRinv)
pred_sm(self, 'dBR/dq2(Lambdab->Lambda(1520)mumu)', 16.4, pow(2.2, -9), pow(0.2, -9), BRinv)
pred_sm(self, 'dBR/dq2(Lambdab->Lambda(1520)mumu)', 16.6, 0.0, pow(0.01, -9), BRinv)
# Comparison to figure on slide 14 S.Meinel b-baryon FEST 2020
pred_sm(self, 'AFBl(Lambdab->Lambda(1520)mumu)', 16.0, -0.16, 0.02)
pred_sm(self, 'AFBl(Lambdab->Lambda(1520)mumu)', 16.4, -0.08, 0.02)
pred_sm(self, 'AFBl(Lambdab->Lambda(1520)mumu)', 16.6, 0.0, 0.02)
pred_sm(self, 'AFBl(Lambdab->Lambda(1520)mumu)', 16.8, 0.22, 0.02)
# Comparison to figure on slide 13 S.Meinel b-baryon FEST 2020
pred_sm(self, 'S_1cc(Lambdab->Lambda(1520)mumu)', 16.0, 0.6, 0.02)
pred_sm(self, 'S_1cc(Lambdab->Lambda(1520)mumu)', 16.4, 0.6, 0.02)
pred_sm(self, 'S_1cc(Lambdab->Lambda(1520)mumu)', 16.8, 0.38, 0.02)