diff options
| author | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-11-06 17:49:20 -0600 |
|---|---|---|
| committer | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-11-06 17:49:20 -0600 |
| commit | 4f358ee507c06c1ccbce472817f4eeac8da5cfa0 (patch) | |
| tree | 0009ba46d75b9c573aea8f106ccfef4ee49ff50a /submitter/mcmc_dag.py | |
| parent | eb98a187b1bc4ee8e8c962d0e5d1a3d0ae77bc35 (diff) | |
| parent | c56c15762ae69e92b805367d0a84091cecaad059 (diff) | |
| download | GolemFlavor-4f358ee507c06c1ccbce472817f4eeac8da5cfa0.tar.gz GolemFlavor-4f358ee507c06c1ccbce472817f4eeac8da5cfa0.zip | |
Merge branch 'master' of github.com:ShiveshM/flavour_ratio
Diffstat (limited to 'submitter/mcmc_dag.py')
| -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()) |