diff --git a/eos-toys/make_toys.py b/eos-toys/make_toys.py index 085f795..5377c21 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-04-06 15:12:03 +# @Last Modified time: 2017-05-10 12:27:28 import os # import subprocess @@ -11,7 +11,7 @@ leptons = ['mu', 'tau'] commands = [] seed = 12345 -form_factors = "BBGIOvD2017-linear" +form_factors = "BBGIOvD2017" signal_pdf = "Lambda_b->Lambda_c(2625)lnu::dGamma" model = "SM" prerun_min = 500000 @@ -21,27 +21,36 @@ 'tau': (3.1570893124, 8.947815864099999), } +os.system('mkdir -p HQT') + commands = [] for l in leptons: rho = { - -0.35: 'linear/mcmc_samples-rho-0.35-%s.hdf5' % l, - # -0.15: 'linear/mcmc_samples-rho-0.15-%s.hdf5' % l, - # -0.05: 'linear/mcmc_samples-rho-0.05-%s.hdf5' % l, - # 0.00: 'linear/mcmc_samples-rho0.00-%s.hdf5' % l, - # 0.05: 'linear/mcmc_samples-rho0.05-%s.hdf5' % l, + # -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, } for r in rho.keys(): - commands.append(['eos-sample-events-mcmc', - ' --fix \"Lambda_b->Lambda_c(2625)::rho@BBGIOvD2017-linear\" %s' % r, - ' --global-option \"form-factors\" \"%s\"' % form_factors, - ' --signal-pdf \"%s,l=%s,model=%s\"' % (signal_pdf, l, model), - ' --seed %s' % seed, - ' --output %s' % rho[r], - ' --prerun-min %s' % prerun_min, - ' --prerun-update %s' % prerun_update, - ' --chunk-size %s' % chunk_size, - ' --chunks %s' % chunks, - ]) + 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), + } + for r3 in rho3b.keys(): + commands.append(['eos-sample-events-mcmc', + ' --fix \"Lambda_b->Lambda_c(2625)::rho@BBGIOvD2017-HQT\" %s' % r, + ' --fix \"Lambda_b->Lambda_c(2625)::rho3b@BBGIOvD2017-HQT\" %s' % r3, + ' --global-option \"form-factors\" \"%s\"' % form_factors, + ' --signal-pdf \"%s;l=%s,model=%s\"' % (signal_pdf, l, model), + ' --seed %s' % seed, + ' --output %s' % rho3b[r3], + ' --prerun-min %s' % prerun_min, + ' --prerun-update %s' % prerun_update, + ' --chunk-size %s' % chunk_size, + ' --chunks %s' % chunks, + ]) for c in commands: print(' '.join(c)) diff --git a/eos-toys/plot-sample.py b/eos-toys/plot-sample.py index 13a3d4b..bbec625 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-04-06 15:40:13 +# @Last Modified time: 2017-05-10 15:50:43 # import matplotlib from __future__ import division import sys, os, gc @@ -62,6 +62,7 @@ def superimpose(files, datafiles, labels, outname): colors = ['#0078FF', '#FF6600', '#0AAFB6', '#FF3333', '#0000FF', '#00CC00', '#BF8040', '#FF33CC', '#FF7733'] + colors.extend(colors) ds = [] for i, f in enumerate(files): try: @@ -99,54 +100,129 @@ 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, - 'linear/mcmc_gamma_muons_linear_rho-0_35.pkl', - 'linear/mcmc_gamma_muons_linear_rho-0_15.pkl', - 'linear/mcmc_gamma_muons_linear_rho-0_05.pkl', - 'linear/mcmc_gamma_muons_linear_rho0_00.pkl', - 'linear/mcmc_gamma_muons_linear_rho0_05.pkl', + # marzia_f, + # 'linear/mcmc_gamma_muons_linear_rho-0_35.pkl', + # 'linear/mcmc_gamma_muons_linear_rho-0_15.pkl', + # 'linear/mcmc_gamma_muons_linear_rho-0_05.pkl', + # 'linear/mcmc_gamma_muons_linear_rho0_00.pkl', + # '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', + 'HQT/mcmc_gamma_muons_HQT_rho0_00_rho3b0_00.pkl', + 'HQT/mcmc_gamma_muons_HQT_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', + # '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', + 'HQT/mcmc_gamma_muons_HQT_rho0_00_rho3b0_05.pkl', + 'HQT/mcmc_gamma_muons_HQT_rho0_05_rho3b0_05.pkl', ] datafiles_mu = ['constantFF/mcmc_gamma-mu.hdf5', '1pole/mcmc_samples-mu.hdf5', '1pole_x_sqrtlambda/mcmc_samples-mu-sqrtlambda.hdf5', - None, - 'linear/mcmc_samples-rho-0.35-mu.hdf5', - 'linear/mcmc_samples-rho-0.15-mu.hdf5', - 'linear/mcmc_samples-rho-0.05-mu.hdf5', - 'linear/mcmc_samples-rho0.00-mu.hdf5', - 'linear/mcmc_samples-rho0.05-mu.hdf5' + # None, + # 'linear/mcmc_samples-rho-0.35-mu.hdf5', + # 'linear/mcmc_samples-rho-0.15-mu.hdf5', + # 'linear/mcmc_samples-rho-0.05-mu.hdf5', + # 'linear/mcmc_samples-rho0.00-mu.hdf5', + # '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', + 'HQT/mcmc_samples-HQT-rho0.0-rho3b0.00-mu.hdf5', + 'HQT/mcmc_samples-HQT-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', + # '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', + 'HQT/mcmc_samples-HQT-rho0.0-rho3b0.05-mu.hdf5', + 'HQT/mcmc_samples-HQT-rho0.05-rho3b0.05-mu.hdf5', ] 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', - marzia_ft, - 'linear/mcmc_gamma_taus_linear_rho-0_35.pkl', - 'linear/mcmc_gamma_taus_linear_rho-0_15.pkl', - 'linear/mcmc_gamma_taus_linear_rho-0_05.pkl', - 'linear/mcmc_gamma_taus_linear_rho0_00.pkl', - 'linear/mcmc_gamma_taus_linear_rho0_05.pkl', + # marzia_ft, + # 'linear/mcmc_gamma_taus_linear_rho-0_35.pkl', + # 'linear/mcmc_gamma_taus_linear_rho-0_15.pkl', + # 'linear/mcmc_gamma_taus_linear_rho-0_05.pkl', + # 'linear/mcmc_gamma_taus_linear_rho0_00.pkl', + # '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', + 'HQT/mcmc_gamma_taus_HQT_rho0_00_rho3b0_00.pkl', + 'HQT/mcmc_gamma_taus_HQT_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', + # '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', + 'HQT/mcmc_gamma_taus_HQT_rho0_00_rho3b0_05.pkl', + 'HQT/mcmc_gamma_taus_HQT_rho0_05_rho3b0_05.pkl', ] datafiles_tau = ['constantFF/mcmc_gamma-tau.hdf5', '1pole/mcmc_samples-tau.hdf5', '1pole_x_sqrtlambda/mcmc_samples-tau-sqrtlambda.hdf5', - None, - 'linear/mcmc_samples-rho-0.35-tau.hdf5', - 'linear/mcmc_samples-rho-0.15-tau.hdf5', - 'linear/mcmc_samples-rho-0.05-tau.hdf5', - 'linear/mcmc_samples-rho0.00-tau.hdf5', - 'linear/mcmc_samples-rho0.05-tau.hdf5', + # None, + # 'linear/mcmc_samples-rho-0.35-tau.hdf5', + # 'linear/mcmc_samples-rho-0.15-tau.hdf5', + # 'linear/mcmc_samples-rho-0.05-tau.hdf5', + # 'linear/mcmc_samples-rho0.00-tau.hdf5', + # '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', + 'HQT/mcmc_samples-HQT-rho0.0-rho3b0.00-tau.hdf5', + 'HQT/mcmc_samples-HQT-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', + # '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', + 'HQT/mcmc_samples-HQT-rho0.0-rho3b0.05-tau.hdf5', + 'HQT/mcmc_samples-HQT-rho0.05-rho3b0.05-tau.hdf5', ] out_tau = 'comparison_tau.pdf' labels = ['Constant FF', 'One pole', r'One pole $\times\; \sqrt{\lambda}$', - 'MB: constant FF', - r'linear, $\rho=-0.35$', - r'linear, $\rho=-0.15$', - r'linear, $\rho=-0.05$', - r'linear, $\rho=0.0$', - r'linear, $\rho=0.05$', + # 'MB: constant FF', + # r'linear, $\rho=-0.35$', + # r'linear, $\rho=-0.15$', + # r'linear, $\rho=-0.05$', + # r'linear, $\rho=0.0$', + # 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'HQT, $\rho=0.0$, $\rho_{3b}=0.00$', + r'HQT, $\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'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$', ] superimpose(files_mu, datafiles_mu, labels, out_mu) superimpose(files_tau, datafiles_tau, labels, out_tau) diff --git a/eos-toys/send_to_zurich.sh b/eos-toys/send_to_zurich.sh index e117244..568cab2 100755 --- a/eos-toys/send_to_zurich.sh +++ b/eos-toys/send_to_zurich.sh @@ -1,15 +1,11 @@ # @Author: Elena Graverini # @Date: 2017-03-21 13:55:51 # @Last Modified by: Elena Graverini -# @Last Modified time: 2017-03-21 16:41:30 +# @Last Modified time: 2017-05-02 17:06:23 #!/bin/bash - -# Needs the following environment variables: -# - $zurich_cluster_mountpoint pointing to the desired mountpoint -# - $zurich_user with your username on the Zurich cluster -mkdir -p $zurich_cluster_mountpoint -sshfs ${zurich_user}@grid-ui.physik.uzh.ch: $zurich_cluster_mountpoint +mkdir -p /home/elena/Desktop/PhD_Work/zurich-cluster +sshfs egraveri@grid-ui.physik.uzh.ch: /home/elena/Desktop/PhD_Work/zurich-cluster cd .. -rsync -rd eos-toys $zurich_cluster_mountpoint/Public/eos-toys +rsync -rd eos-toys /home/elena/Desktop/PhD_Work/zurich-cluster/Public/eos-toys cd eos-toys -fusermount -u $zurich_cluster_mountpoint || sudo umount -l $zurich_cluster_mountpoint +fusermount -u /home/elena/Desktop/PhD_Work/zurich-cluster || sudo umount -l /home/elena/Documents/University/PhD/UZH/zurich-cluster