rescale bvll subleading a_i parameters

Author: peterstangl

flavio/data/parameters_correlated.yml

@@ -106,9 +106,9 @@
 # Subleading hadronic effects in B->Vll
 -
   values:
-    - B0->K*0 deltaC7 a_0 Re: 0 ± 0.12
-    - B+->K*+ deltaC7 a_0 Re: 0 ± 0.12
-    - Bs->phi deltaC7 a_0 Re: 0 ± 0.12
+    - B0->K*0 deltaC7 a_0 Re: 0 ± 0.186 # ± 0.12 * (2*mb/mB)
+    - B+->K*+ deltaC7 a_0 Re: 0 ± 0.186
+    - Bs->phi deltaC7 a_0 Re: 0 ± 0.186
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
@@ -120,9 +120,9 @@
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
   values:
-    - B0->K*0 deltaC7p a_+ Re: 0 ± 0.004
-    - B+->K*+ deltaC7p a_+ Re: 0 ± 0.004
-    - Bs->phi deltaC7p a_+ Re: 0 ± 0.004
+    - B0->K*0 deltaC7p a_+ Re: 0 ± 0.006 # ± 0.004 * (2*mb/mB)
+    - B+->K*+ deltaC7p a_+ Re: 0 ± 0.006
+    - Bs->phi deltaC7p a_+ Re: 0 ± 0.006
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
@@ -133,10 +133,10 @@
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
-  values:
-    - B0->K*0 deltaC7 a_- Re: 0 ± 0.015 # see eq. (26) of arXiv:1503.05534v1
-    - B+->K*+ deltaC7 a_- Re: 0 ± 0.015
-    - Bs->phi deltaC7 a_- Re: 0 ± 0.015
+  values: # see eq. (26) of arXiv:1503.05534v1
+    - B0->K*0 deltaC7 a_- Re: 0 ± 0.023 # ± 0.015 * (2*mb/mB)
+    - B+->K*+ deltaC7 a_- Re: 0 ± 0.023
+    - Bs->phi deltaC7 a_- Re: 0 ± 0.023
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
@@ -147,10 +147,10 @@
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
-  values:
-    - B0->K*0 deltaC7 a_0 Im: 0 ± 0.12 # see eq. (26) of arXiv:1503.05534v1
-    - B+->K*+ deltaC7 a_0 Im: 0 ± 0.12
-    - Bs->phi deltaC7 a_0 Im: 0 ± 0.12
+  values: # see eq. (26) of arXiv:1503.05534v1
+    - B0->K*0 deltaC7 a_0 Im: 0 ± 0.186 # ± 0.12 * (2*mb/mB)
+    - B+->K*+ deltaC7 a_0 Im: 0 ± 0.186
+    - Bs->phi deltaC7 a_0 Im: 0 ± 0.186
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
@@ -162,9 +162,9 @@
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
   values:
-    - B0->K*0 deltaC7p a_+ Im: 0 ± 0.004
-    - B+->K*+ deltaC7p a_+ Im: 0 ± 0.004
-    - Bs->phi deltaC7p a_+ Im: 0 ± 0.004
+    - B0->K*0 deltaC7p a_+ Im: 0 ± 0.006 # ± 0.004 * (2*mb/mB)
+    - B+->K*+ deltaC7p a_+ Im: 0 ± 0.006
+    - Bs->phi deltaC7p a_+ Im: 0 ± 0.006
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
@@ -175,10 +175,10 @@
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -
-  values:
-    - B0->K*0 deltaC7 a_- Im: 0 ± 0.015 # see eq. (26) of arXiv:1503.05534v1
-    - B+->K*+ deltaC7 a_- Im: 0 ± 0.015
-    - Bs->phi deltaC7 a_- Im: 0 ± 0.015
+  values: # see eq. (26) of arXiv:1503.05534v1
+    - B0->K*0 deltaC7 a_- Im: 0 ± 0.023 # ± 0.015 * (2*mb/mB)
+    - B+->K*+ deltaC7 a_- Im: 0 ± 0.023
+    - Bs->phi deltaC7 a_- Im: 0 ± 0.023
   correlation:
     [[1, 0.99, 0.9], [1, 0.9], [1]] # B+/B0 correlation is 99%, B0/Bs is 90%
 -

flavio/physics/bdecays/bvll/subleading.py

@@ -52,13 +52,12 @@ def __call__(self):
         q2 = self.q2
         mB = self.mB
         mb = self.mb
-        q2_factor = q2/(2*mb*mB)
-        deltaC7_0   =( par[B+'->'+V+' deltaC7 a_0 Re']  + par[B+'->'+V+' deltaC7 b_0 Re'] *q2_factor
-                 +1j*( par[B+'->'+V+' deltaC7 a_0 Im']  + par[B+'->'+V+' deltaC7 b_0 Im'] *q2_factor ))
-        deltaC7p_pl  =( par[B+'->'+V+' deltaC7p a_+ Re']  + par[B+'->'+V+' deltaC7p b_+ Re'] *q2_factor
-                 +1j*( par[B+'->'+V+' deltaC7p a_+ Im']  + par[B+'->'+V+' deltaC7p b_+ Im'] *q2_factor ))
-        deltaC7_mi  =( par[B+'->'+V+' deltaC7 a_- Re']  + par[B+'->'+V+' deltaC7 b_- Re'] *q2_factor
-                 +1j*( par[B+'->'+V+' deltaC7 a_- Im']  + par[B+'->'+V+' deltaC7 b_- Im'] *q2_factor ))
+        deltaC7_0   =( par[B+'->'+V+' deltaC7 a_0 Re']  + par[B+'->'+V+' deltaC7 b_0 Re']  *q2/mB**2
+                 +1j*( par[B+'->'+V+' deltaC7 a_0 Im']  + par[B+'->'+V+' deltaC7 b_0 Im']  *q2/mB**2 )) *mB/(2*mb)
+        deltaC7p_pl =( par[B+'->'+V+' deltaC7p a_+ Re'] + par[B+'->'+V+' deltaC7p b_+ Re'] *q2/mB**2
+                 +1j*( par[B+'->'+V+' deltaC7p a_+ Im'] + par[B+'->'+V+' deltaC7p b_+ Im'] *q2/mB**2 )) *mB/(2*mb)
+        deltaC7_mi  =( par[B+'->'+V+' deltaC7 a_- Re']  + par[B+'->'+V+' deltaC7 b_- Re']  *q2/mB**2
+                 +1j*( par[B+'->'+V+' deltaC7 a_- Im']  + par[B+'->'+V+' deltaC7 b_- Im']  *q2/mB**2 )) *mB/(2*mb)
         ha = {}
         ha['0', 'V'] = self.ha_deltaC(deltaC7_0, '7')['0', 'V']
         ha['pl', 'V'] = self.ha_deltaC(deltaC7p_pl, '7p')['pl', 'V']