diff options
| author | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-05-08 13:17:01 -0500 |
|---|---|---|
| committer | shivesh <s.p.mandalia@qmul.ac.uk> | 2018-05-08 13:17:01 -0500 |
| commit | e9b62e357c95eec6bc3cb28e5eae9d99a6556e37 (patch) | |
| tree | 9c04971e23bcadbfede577845f12bf5ec701d71a /plot_llh/mcmc_scan.py | |
| parent | 5f76bd5bb0bbe4e1ebc85a60722498f3747fd722 (diff) | |
| download | GolemFlavor-e9b62e357c95eec6bc3cb28e5eae9d99a6556e37.tar.gz GolemFlavor-e9b62e357c95eec6bc3cb28e5eae9d99a6556e37.zip | |
naming in mcmc_scan for plot_llh
Diffstat (limited to 'plot_llh/mcmc_scan.py')
| -rw-r--r-- | plot_llh/mcmc_scan.py | 10 |
1 files changed, 8 insertions, 2 deletions
diff --git a/plot_llh/mcmc_scan.py b/plot_llh/mcmc_scan.py index 455d5dd..d1f9e1b 100644 --- a/plot_llh/mcmc_scan.py +++ b/plot_llh/mcmc_scan.py @@ -20,6 +20,11 @@ MEASURED_FR = [1, 1, 1] SIGMA = 0.001 +def solve_ratio(fr): + denominator = reduce(gcd, fr) + return [int(x/denominator) for x in fr] + + def angles_to_fr(angles): sphi4, c2psi = angles @@ -145,8 +150,9 @@ def main(): pass print "Finished" - outfile = args.outfile+'_{0:03d}_{1:03d}_{2:03d}_{3:.1E}'.format( - int(MEASURED_FR[0]*100), int(MEASURED_FR[1]*100), int(MEASURED_FR[2]*100), SIGMA + mr = solve_ratio(MEASURED_FR) + outfile = args.outfile+'_{0}_{1}_{2}_{3:.1E}'.format( + mr[0], mr[1], mr[2], SIGMA ) samples = sampler.chain[0, :, :, :].reshape((-1, ndim)) |