Merge a9e15c1 into 7e0502c

flavio/data/measurements.yml

@@ -72,7 +72,11 @@ LHCb rad 2012:
   experiment: LHCb
   inspire: Aaij:2012ita
   values:
-    BR(Bs->phigamma): 3.5 ± 0.4 e-5
+    # The originally measured value is 1.23 ± 0.06 ± 0.04 ± 0.10
+    # Using a more precise value of fs/fd from https://cds.cern.ch/record/1559262/
+    # one gets 1.19 ± 0.06 ± 0.04 ± 0.07
+    # see https://github.com/flav-io/flavio/issues/145#issuecomment-788085708
+    BR(B0->K*gamma)/BR(Bs->phigamma): 1.19 ± 0.06 ± 0.04 ± 0.07
 
 LHCb phigamma 2019:
   experiment: LHCb

flavio/physics/bdecays/bvgamma.py

@@ -136,6 +136,12 @@ def BR_timeint(wc_obj, par, B, V):
 def BVgamma_function(function, B, V):
     return lambda wc_obj, par: function(wc_obj, par, B, V)
 
+def BVgamma_ratio_function(function_num, B_num, V_num, function_den, B_den, V_den):
+    return lambda wc_obj, par: (
+        function_num(wc_obj, par, B_num, V_num)
+        /
+        function_den(wc_obj, par, B_den, V_den)
+    )
 
 # Observable and Prediction instances
 
@@ -179,6 +185,23 @@ def BVgamma_function(function, B, V):
 _obs.add_taxonomy(_process_taxonomy + _process_tex + r"$")
 Prediction(_obs_name, BVgamma_function(BR_timeint, 'Bs', 'phi'))
 
+_obs_name = "BR(B0->K*gamma)/BR(Bs->phigamma)"
+_obs = Observable(_obs_name)
+_process_num_tex = r"B^0\to K^{*0}\gamma"
+_process_den_tex = r"B_s\to \phi\gamma"
+_obs.set_description(
+    fr"Ratio of branching ratio of ${_process_num_tex}$ and time-integrated "
+    fr"branching ratio of ${_process_den_tex}$"
+)
+_obs.tex = (
+    fr"$\frac{{\text{{BR}}({_process_num_tex})}}"
+    fr"{{\overline{{\text{{BR}}}}({_process_den_tex})}}$"
+)
+_obs.add_taxonomy(_process_taxonomy + _process_den_tex + r"$")
+Prediction(_obs_name, BVgamma_ratio_function(
+    BR, 'B0', 'K*0', BR_timeint, 'Bs', 'phi'
+))
+
 _obs_name = "ADeltaGamma(Bs->phigamma)"
 _obs = Observable(_obs_name)
 _process_tex = r"B_s\to \phi\gamma"