diff --git a/eos-toys/plot-sample.py b/eos-toys/plot-sample.py index bf979bc..e3dc8b8 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-03-21 11:38:40 +# @Last Modified time: 2017-03-22 08:25:24 # import matplotlib from __future__ import division import sys, os, gc @@ -69,13 +69,17 @@ if __name__ == '__main__': # main() - marzia_f = 'Marzia/Marzia_constantFF.pkl' + marzia_f = 'Marzia/Marzia_mu_constantFF.pkl' + marzia_ft = 'Marzia/Marzia_tau_constantFF.pkl' mx = [0.2, 0.4, 0.80, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5] - my = [0.0102553, 0.0110076, 0.0123026, 0.0128953, 0.0142323, 0.015366, 0.0162975, 0.0170258, 0.0175493, 0.0178644, 0.0179659, 0.0178455, 0.017491, 0.016884, 0.0159959, 0.0147803, 0.0131556, 0.0109562, 0.00773188] + my = [0.147214, 0.14805, 0.147707, 0.147292, 0.145776, 0.143623, 0.140846, 0.137443, 0.133406, 0.12872, 0.123359, 0.117288, 0.110451, 0.102761, 0.0940879, 0.0842141, 0.0727531, 0.0589096, 0.0404804] + mxt = [3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5] + myt = [0.00641833, 0.0194386, 0.0789258, 0.140122, 0.189971, 0.225482, 0.24678, 0.254695, 0.249817, 0.231917, 0.198965, 0.143346] makeFileFromMarzia(mx, my, marzia_f) + makeFileFromMarzia(mxt, myt, marzia_ft) files_mu = ['constantFF/mcmc_gamma_muons.pkl', '1pole/mcmc_gamma_muons_1pole.pkl', '1pole_x_sqrtlambda/mcmc_gamma_muonss_1pole_x_sqrtlambda.pkl', marzia_f] out_mu = 'comparison_mu.pdf' - files_tau = ['constantFF/mcmc_gamma_taus.pkl', '1pole/mcmc_gamma_taus_1pole.pkl', '1pole_x_sqrtlambda/mcmc_gamma_taus_1pole_x_sqrtlambda.pkl'] + files_tau = ['constantFF/mcmc_gamma_taus.pkl', '1pole/mcmc_gamma_taus_1pole.pkl', '1pole_x_sqrtlambda/mcmc_gamma_taus_1pole_x_sqrtlambda.pkl', marzia_ft] out_tau = 'comparison_tau.pdf' labels = ['Constant FF', 'One pole', r'One pole $\times\; \sqrt{\lambda}$', 'MB: constant FF'] superimpose(files_mu, labels, out_mu)