diff --git a/eos-toys/make_toys.py b/eos-toys/make_toys.py
index 5377c21..58bd8fc 100644
--- a/eos-toys/make_toys.py
+++ b/eos-toys/make_toys.py
@@ -2,7 +2,7 @@
 # @Author: Elena Graverini
 # @Date:   2017-03-16 15:22:00
 # @Last Modified by:   Elena Graverini
-# @Last Modified time: 2017-05-10 12:27:28
+# @Last Modified time: 2017-07-04 09:23:44
 import os
 # import subprocess
 
@@ -28,15 +28,15 @@
     rho = {
         # -0.35: 'HQT/mcmc_samples-HQT-rho-0.35-%s.hdf5' % l,
         # -0.15: 'HQT/mcmc_samples-HQT-rho-0.15-%s.hdf5' % l,
-        -0.05: 'HQT/mcmc_samples-HQT-rho-0.05-%s.hdf5' % l,
-        0.00: 'HQT/mcmc_samples-HQT-rho0.00-%s.hdf5' % l,
-        0.05: 'HQT/mcmc_samples-HQT-rho0.05-%s.hdf5' % l,
+        -0.05: 'MSBAR/mcmc_samples-MSBAR-rho-0.05-%s.hdf5' % l,
+        0.00: 'MSBAR/mcmc_samples-MSBAR-rho0.00-%s.hdf5' % l,
+        0.05: 'MSBAR/mcmc_samples-MSBAR-rho0.05-%s.hdf5' % l,
     }
     for r in rho.keys():
         rho3b = {
-            -0.05: 'HQT/mcmc_samples-HQT-rho%s-rho3b-0.05-%s.hdf5' % (r, l),
-            0.00: 'HQT/mcmc_samples-HQT-rho%s-rho3b0.00-%s.hdf5' % (r, l),
-            0.05: 'HQT/mcmc_samples-HQT-rho%s-rho3b0.05-%s.hdf5' % (r, l),
+            -0.05: 'MSBAR/mcmc_samples-MSBAR-rho%s-rho3b-0.05-%s.hdf5' % (r, l),
+            0.00: 'MSBAR/mcmc_samples-MSBAR-rho%s-rho3b0.00-%s.hdf5' % (r, l),
+            0.05: 'MSBAR/mcmc_samples-MSBAR-rho%s-rho3b0.05-%s.hdf5' % (r, l),
         }
         for r3 in rho3b.keys():
             commands.append(['eos-sample-events-mcmc',
diff --git a/eos-toys/scan-dBRds.py b/eos-toys/scan-dBRds.py
index 775f1ac..243d419 100644
--- a/eos-toys/scan-dBRds.py
+++ b/eos-toys/scan-dBRds.py
@@ -2,7 +2,7 @@
 # @Author: Elena Graverini
 # @Date:   2017-03-22 11:15:24
 # @Last Modified by:   Elena Graverini
-# @Last Modified time: 2017-03-22 12:01:00
+# @Last Modified time: 2017-07-04 11:39:29
 import os
 tempfile = 'out.temp'
 step = 0.5
@@ -18,7 +18,7 @@
 for l in bounds.keys():
     os.system("for i in $(LANG=en_US seq {min} {step} {max});"
               "do out=$(eos-evaluate  --kinematics s $i"
-              " --observable 'Lambda_b->Lambda_c(2625)lnu::dBR/ds;l={l}');"
+              " --observable 'Lambda_b->Lambda_c(2625)lnu::dBR/ds;l={l},model=SM');"
               " echo $out; done > {tempf}".format(min=bounds[l][0],
                     max=bounds[l][1], l=l, tempf=tempfile, step=step))
     with open(tempfile, 'rb') as f:
@@ -35,3 +35,6 @@
           "\t{dbr:>8} = {br}\n".format(l=l,
             q2='q^2', dbr='dBR/dq^2',
             s=scan[l]['s'], br=scan[l]['br']))
+
+# Total BR:
+# eos-evaluate --kinematics s_min 0.012 --kinematics s_max 8.948 --observable 'Lambda_b->Lambda_c(2625)lnu::BR;l=mu,model=SM'