diff options
| author | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-11-06 14:23:10 -0600 |
|---|---|---|
| committer | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-11-06 14:23:10 -0600 |
| commit | 441e5e8101b03d6c15bc2370700ee14d07375e92 (patch) | |
| tree | 543544161ceb6c5a1732cfa084191b61fedd99d2 /submitter | |
| parent | 84c05583ad3d582d107fe09be700311ea466b1af (diff) | |
| download | GolemFlavor-441e5e8101b03d6c15bc2370700ee14d07375e92.tar.gz GolemFlavor-441e5e8101b03d6c15bc2370700ee14d07375e92.zip | |
Tue 6 Nov 14:23:10 CST 2018
Diffstat (limited to 'submitter')
| -rw-r--r-- | submitter/mcmc_dag.py | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/submitter/mcmc_dag.py b/submitter/mcmc_dag.py index b7e2c84..5a44411 100644 --- a/submitter/mcmc_dag.py +++ b/submitter/mcmc_dag.py @@ -25,7 +25,7 @@ GLOBAL_PARAMS = {} # MCMC GLOBAL_PARAMS.update(dict( - run_mcmc = 'False', + run_mcmc = 'True', burnin = 250, nsteps = 1000, nwalkers = 60, @@ -34,9 +34,9 @@ GLOBAL_PARAMS.update(dict( )) # FR -# dimension = [3, 6] +dimension = [6] # dimension = [4, 5, 7, 8] -dimension = [3, 4, 5, 6, 7, 8] +# dimension = [3, 4, 5, 6, 7, 8] GLOBAL_PARAMS.update(dict( threads = 1, binning = '6e4 1e7 20', @@ -44,7 +44,7 @@ GLOBAL_PARAMS.update(dict( scale_region = "1E10", energy_dependance = 'spectral', spectral_index = -2, - fix_mixing = 'None', + fix_mixing = 'T23', fix_mixing_almost = 'False', fold_index = 'True' )) @@ -67,7 +67,8 @@ GLOBAL_PARAMS.update(dict( plot_elements = 'False', )) -outfile = 'dagman_FR_MCMC_{0}.submit'.format(GLOBAL_PARAMS['likelihood']) +outfile = 'dagman_FR_MCMC_{0}_{1}.submit'.format(GLOBAL_PARAMS['likelihood'], + GLOBAL_PARAMS['fix_mixing']) golemfitsourcepath = os.environ['GOLEMSOURCEPATH'] + '/GolemFit' condor_script = golemfitsourcepath + '/scripts/flavour_ratio/submitter/mcmc_submit.sub' @@ -80,7 +81,7 @@ with open(outfile, 'w') as f: ) for frs in fix_sfr_mfr: print 'frs', frs - outchains = outchain_head + '/fix_ifr/' + outchains = outchain_head + '/fix_ifr/' + '{0}/'.format(GLOBAL_PARAMS['fix_mixing']) if GLOBAL_PARAMS['likelihood'].lower() == 'gaussian': outchains += '{0}/'.format(str(GLOBAL_PARAMS['sigma_ratio']).replace('.', '_')) outchains += '{0}/'.format(GLOBAL_PARAMS['data'].lower()) |