diff --git a/eos-toys/plot-sample.py b/eos-toys/plot-sample.py
index 6d4bd0c..844e01b 100644
--- a/eos-toys/plot-sample.py
+++ b/eos-toys/plot-sample.py
@@ -2,7 +2,7 @@
 # @Author: Elena Graverini
 # @Date:   2017-03-16 11:20:45
 # @Last Modified by:   Elena Graverini
-# @Last Modified time: 2017-05-10 16:52:40
+# @Last Modified time: 2017-07-05 11:04:23
 # import matplotlib
 from __future__ import division
 import sys, os, gc
@@ -120,19 +120,23 @@
                 # 'linear/mcmc_gamma_muons_linear_rho0_05.pkl',
                 # 'HQT/mcmc_gamma_muons_HQT_rho-0_35_rho3b0_00.pkl',
                 # 'HQT/mcmc_gamma_muons_HQT_rho-0_15_rho3b0_00.pkl',
-                'HQT/mcmc_gamma_muons_HQT_rho-0_05_rho3b0_00.pkl',
+                # 'MSBAR/mcmc_gamma_muons_MSBAR_rho-0_05_rho3b0_00.pkl',
                 'HQT/mcmc_gamma_muons_HQT_rho0_00_rho3b0_00.pkl',
                 'HQT/mcmc_gamma_muons_HQT_rho0_05_rho3b0_00.pkl',
+                'MSBAR/mcmc_gamma_muons_MSBAR_rho0_00_rho3b0_00.pkl',
+                'MSBAR/mcmc_gamma_muons_MSBAR_rho0_05_rho3b0_00.pkl',
                 # 'HQT/mcmc_gamma_muons_HQT_rho-0_35_rho3b-0_05.pkl',
                 # 'HQT/mcmc_gamma_muons_HQT_rho-0_15_rho3b-0_05.pkl',
-                'HQT/mcmc_gamma_muons_HQT_rho-0_05_rho3b-0_05.pkl',
-                'HQT/mcmc_gamma_muons_HQT_rho0_00_rho3b-0_05.pkl',
-                'HQT/mcmc_gamma_muons_HQT_rho0_05_rho3b-0_05.pkl',
+                # 'MSBAR/mcmc_gamma_muons_MSBAR_rho-0_05_rho3b-0_05.pkl',
+                # 'MSBAR/mcmc_gamma_muons_MSBAR_rho0_00_rho3b-0_05.pkl',
+                # 'MSBAR/mcmc_gamma_muons_MSBAR_rho0_05_rho3b-0_05.pkl',
                 # 'HQT/mcmc_gamma_muons_HQT_rho-0_35_rho3b0_05.pkl',
                 # 'HQT/mcmc_gamma_muons_HQT_rho-0_15_rho3b0_05.pkl',
-                'HQT/mcmc_gamma_muons_HQT_rho-0_05_rho3b0_05.pkl',
+                # 'MSBAR/mcmc_gamma_muons_MSBAR_rho-0_05_rho3b0_05.pkl',
                 'HQT/mcmc_gamma_muons_HQT_rho0_00_rho3b0_05.pkl',
                 'HQT/mcmc_gamma_muons_HQT_rho0_05_rho3b0_05.pkl',
+                'MSBAR/mcmc_gamma_muons_MSBAR_rho0_00_rho3b0_05.pkl',
+                'MSBAR/mcmc_gamma_muons_MSBAR_rho0_05_rho3b0_05.pkl',
                 ]
     datafiles_mu = ['constantFF/mcmc_gamma-mu.hdf5',
                     '1pole/mcmc_samples-mu.hdf5',
@@ -145,19 +149,23 @@
                     # 'linear/mcmc_samples-rho0.05-mu.hdf5',
                     # 'HQT/mcmc_samples-HQT-rho-0.35-rho3b0.00-mu.hdf5',
                     # 'HQT/mcmc_samples-HQT-rho-0.15-rho3b0.00-mu.hdf5',
-                    'HQT/mcmc_samples-HQT-rho-0.05-rho3b0.00-mu.hdf5',
+                    # 'MSBAR/mcmc_samples-MSBAR-rho-0.05-rho3b0.00-mu.hdf5',
                     'HQT/mcmc_samples-HQT-rho0.0-rho3b0.00-mu.hdf5',
                     'HQT/mcmc_samples-HQT-rho0.05-rho3b0.00-mu.hdf5',
+                    'MSBAR/mcmc_samples-MSBAR-rho0.0-rho3b0.00-mu.hdf5',
+                    'MSBAR/mcmc_samples-MSBAR-rho0.05-rho3b0.00-mu.hdf5',
                     # 'HQT/mcmc_samples-HQT-rho-0.35-rho3b-0.05-mu.hdf5',
                     # 'HQT/mcmc_samples-HQT-rho-0.15-rho3b-0.05-mu.hdf5',
-                    'HQT/mcmc_samples-HQT-rho-0.05-rho3b-0.05-mu.hdf5',
-                    'HQT/mcmc_samples-HQT-rho0.0-rho3b-0.05-mu.hdf5',
-                    'HQT/mcmc_samples-HQT-rho0.05-rho3b-0.05-mu.hdf5',
+                    # 'MSBAR/mcmc_samples-MSBAR-rho-0.05-rho3b-0.05-mu.hdf5',
+                    # 'MSBAR/mcmc_samples-MSBAR-rho0.0-rho3b-0.05-mu.hdf5',
+                    # 'MSBAR/mcmc_samples-MSBAR-rho0.05-rho3b-0.05-mu.hdf5',
                     # 'HQT/mcmc_samples-HQT-rho-0.35-rho3b0.05-mu.hdf5',
                     # 'HQT/mcmc_samples-HQT-rho-0.15-rho3b0.05-mu.hdf5',
-                    'HQT/mcmc_samples-HQT-rho-0.05-rho3b0.05-mu.hdf5',
+                    # 'MSBAR/mcmc_samples-MSBAR-rho-0.05-rho3b0.05-mu.hdf5',
                     'HQT/mcmc_samples-HQT-rho0.0-rho3b0.05-mu.hdf5',
                     'HQT/mcmc_samples-HQT-rho0.05-rho3b0.05-mu.hdf5',
+                    'MSBAR/mcmc_samples-MSBAR-rho0.0-rho3b0.05-mu.hdf5',
+                    'MSBAR/mcmc_samples-MSBAR-rho0.05-rho3b0.05-mu.hdf5',
                     ]
     out_mu = 'comparison_mu.pdf'
     files_tau = ['constantFF/mcmc_gamma_taus.pkl',
@@ -171,19 +179,23 @@
                  # 'linear/mcmc_gamma_taus_linear_rho0_05.pkl',
                  # 'HQT/mcmc_gamma_taus_HQT_rho-0_35_rho3b0_00.pkl',
                  # 'HQT/mcmc_gamma_taus_HQT_rho-0_15_rho3b0_00.pkl',
-                 'HQT/mcmc_gamma_taus_HQT_rho-0_05_rho3b0_00.pkl',
+                 # 'MSBAR/mcmc_gamma_taus_MSBAR_rho-0_05_rho3b0_00.pkl',
                  'HQT/mcmc_gamma_taus_HQT_rho0_00_rho3b0_00.pkl',
                  'HQT/mcmc_gamma_taus_HQT_rho0_05_rho3b0_00.pkl',
+                 'MSBAR/mcmc_gamma_taus_MSBAR_rho0_00_rho3b0_00.pkl',
+                 'MSBAR/mcmc_gamma_taus_MSBAR_rho0_05_rho3b0_00.pkl',
                  # 'HQT/mcmc_gamma_taus_HQT_rho-0_35_rho3b-0_05.pkl',
                  # 'HQT/mcmc_gamma_taus_HQT_rho-0_15_rho3b-0_05.pkl',
-                 'HQT/mcmc_gamma_taus_HQT_rho-0_05_rho3b-0_05.pkl',
-                 'HQT/mcmc_gamma_taus_HQT_rho0_00_rho3b-0_05.pkl',
-                 'HQT/mcmc_gamma_taus_HQT_rho0_05_rho3b-0_05.pkl',
+                 # 'MSBAR/mcmc_gamma_taus_MSBAR_rho-0_05_rho3b-0_05.pkl',
+                 # 'MSBAR/mcmc_gamma_taus_MSBAR_rho0_00_rho3b-0_05.pkl',
+                 # 'MSBAR/mcmc_gamma_taus_MSBAR_rho0_05_rho3b-0_05.pkl',
                  # 'HQT/mcmc_gamma_taus_HQT_rho-0_35_rho3b0_05.pkl',
                  # 'HQT/mcmc_gamma_taus_HQT_rho-0_15_rho3b0_05.pkl',
-                 'HQT/mcmc_gamma_taus_HQT_rho-0_05_rho3b0_05.pkl',
+                 # 'MSBAR/mcmc_gamma_taus_MSBAR_rho-0_05_rho3b0_05.pkl',
                  'HQT/mcmc_gamma_taus_HQT_rho0_00_rho3b0_05.pkl',
                  'HQT/mcmc_gamma_taus_HQT_rho0_05_rho3b0_05.pkl',
+                 'MSBAR/mcmc_gamma_taus_MSBAR_rho0_00_rho3b0_05.pkl',
+                 'MSBAR/mcmc_gamma_taus_MSBAR_rho0_05_rho3b0_05.pkl',
                  ]
     datafiles_tau = ['constantFF/mcmc_gamma-tau.hdf5',
                      '1pole/mcmc_samples-tau.hdf5',
@@ -196,19 +208,23 @@
                      # 'linear/mcmc_samples-rho0.05-tau.hdf5',
                      # 'HQT/mcmc_samples-HQT-rho-0.35-rho3b0.00-tau.hdf5',
                      # 'HQT/mcmc_samples-HQT-rho-0.15-rho3b0.00-tau.hdf5',
-                     'HQT/mcmc_samples-HQT-rho-0.05-rho3b0.00-tau.hdf5',
+                     # 'MSBAR/mcmc_samples-MSBAR-rho-0.05-rho3b0.00-tau.hdf5',
                      'HQT/mcmc_samples-HQT-rho0.0-rho3b0.00-tau.hdf5',
                      'HQT/mcmc_samples-HQT-rho0.05-rho3b0.00-tau.hdf5',
+                     'MSBAR/mcmc_samples-MSBAR-rho0.0-rho3b0.00-tau.hdf5',
+                     'MSBAR/mcmc_samples-MSBAR-rho0.05-rho3b0.00-tau.hdf5',
                      # 'HQT/mcmc_samples-HQT-rho-0.35-rho3b-0.05-tau.hdf5',
                      # 'HQT/mcmc_samples-HQT-rho-0.15-rho3b-0.05-tau.hdf5',
-                     'HQT/mcmc_samples-HQT-rho-0.05-rho3b-0.05-tau.hdf5',
-                     'HQT/mcmc_samples-HQT-rho0.0-rho3b-0.05-tau.hdf5',
-                     'HQT/mcmc_samples-HQT-rho0.05-rho3b-0.05-tau.hdf5',
+                     # 'MSBAR/mcmc_samples-MSBAR-rho-0.05-rho3b-0.05-tau.hdf5',
+                     # 'MSBAR/mcmc_samples-MSBAR-rho0.0-rho3b-0.05-tau.hdf5',
+                     # 'MSBAR/mcmc_samples-MSBAR-rho0.05-rho3b-0.05-tau.hdf5',
                      # 'HQT/mcmc_samples-HQT-rho-0.35-rho3b0.05-tau.hdf5',
                      # 'HQT/mcmc_samples-HQT-rho-0.15-rho3b0.05-tau.hdf5',
-                     'HQT/mcmc_samples-HQT-rho-0.05-rho3b0.05-tau.hdf5',
+                     # 'MSBAR/mcmc_samples-MSBAR-rho-0.05-rho3b0.05-tau.hdf5',
                      'HQT/mcmc_samples-HQT-rho0.0-rho3b0.05-tau.hdf5',
                      'HQT/mcmc_samples-HQT-rho0.05-rho3b0.05-tau.hdf5',
+                     'MSBAR/mcmc_samples-MSBAR-rho0.0-rho3b0.05-tau.hdf5',
+                     'MSBAR/mcmc_samples-MSBAR-rho0.05-rho3b0.05-tau.hdf5',
                      ]
     out_tau = 'comparison_tau.pdf'
     labels = ['Constant FF',
@@ -222,19 +238,23 @@
               # r'linear, $\rho=0.05$',
               # r'HQT, $\rho=-0.35$, $\rho_{3b}=0.00$',
               # r'HQT, $\rho=-0.15$, $\rho_{3b}=0.00$',
-              r'HQT, $\rho=-0.05$, $\rho_{3b}=0.00$',
+              # r'MSBAR, $\rho=-0.05$, $\rho_{3b}=0.00$',
               r'HQT, $\rho=0.0$, $\rho_{3b}=0.00$',
               r'HQT, $\rho=0.05$, $\rho_{3b}=0.00$',
+              r'MSBAR, $\rho=0.0$, $\rho_{3b}=0.00$',
+              r'MSBAR, $\rho=0.05$, $\rho_{3b}=0.00$',
               # r'HQT, $\rho=-0.35$, $\rho_{3b}=-0.05$',
               # r'HQT, $\rho=-0.15$, $\rho_{3b}=-0.05$',
-              r'HQT, $\rho=-0.05$, $\rho_{3b}=-0.05$',
-              r'HQT, $\rho=0.0$, $\rho_{3b}=-0.05$',
-              r'HQT, $\rho=0.05$, $\rho_{3b}=-0.05$',
+              # r'MSBAR, $\rho=-0.05$, $\rho_{3b}=-0.05$',
+              # r'MSBAR, $\rho=0.0$, $\rho_{3b}=-0.05$',
+              # r'MSBAR, $\rho=0.05$, $\rho_{3b}=-0.05$',
               # r'HQT, $\rho=-0.35$, $\rho_{3b}=0.05$',
               # r'HQT, $\rho=-0.15$, $\rho_{3b}=0.05$',
-              r'HQT, $\rho=-0.05$, $\rho_{3b}=0.05$',
+              # r'MSBAR, $\rho=-0.05$, $\rho_{3b}=0.05$',
               r'HQT, $\rho=0.0$, $\rho_{3b}=0.05$',
               r'HQT, $\rho=0.05$, $\rho_{3b}=0.05$',
+              r'MSBAR, $\rho=0.0$, $\rho_{3b}=0.05$',
+              r'MSBAR, $\rho=0.05$, $\rho_{3b}=0.05$',
               ]
     superimpose(files_mu, datafiles_mu, labels, out_mu, 'mu')
     superimpose(files_tau, datafiles_tau, labels, out_tau, 'tau')