From 23899b5a0832fd92984cfa0a56bd2fdae03b9b0d Mon Sep 17 00:00:00 2001 From: Richard O'Shaughnessy Date: Sat, 12 Jan 2019 19:20:06 -0500 Subject: [PATCH] CIP: Modify to allow distance fitting --- ...util_ConstructIntrinsicPosterior_GenericCoordinates.py | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/MonteCarloMarginalizeCode/Code/util_ConstructIntrinsicPosterior_GenericCoordinates.py b/MonteCarloMarginalizeCode/Code/util_ConstructIntrinsicPosterior_GenericCoordinates.py index 3fad2fd6..bc569ca5 100755 --- a/MonteCarloMarginalizeCode/Code/util_ConstructIntrinsicPosterior_GenericCoordinates.py +++ b/MonteCarloMarginalizeCode/Code/util_ConstructIntrinsicPosterior_GenericCoordinates.py @@ -175,7 +175,8 @@ def add_field(a, descr): parser = argparse.ArgumentParser() parser.add_argument("--fname",help="filename of *.dat file [standard ILE output]") -parser.add_argument("--input-tides",action='store_true',help="Use input format with tidal fields included.") +parser.add_argument("--input-tides",action='store_true',help="Use input format with tidal fields included.") +parser.add_argument("--input-distance",action='store_true',help="Use input format with distance fields (but not tidal fields?) enabled.") parser.add_argument("--fname-lalinference",help="filename of posterior_samples.dat file [standard LI output], to overlay on corner plots") parser.add_argument("--fname-output-samples",default="output-ILE-samples",help="output posterior samples (default output-ILE-samples -> output-ILE)") parser.add_argument("--fname-output-integral",default="integral_result",help="output filename for integral result. Postfixes appended") @@ -924,6 +925,9 @@ col_lnL = 9 if opts.input_tides: print " Tides input" col_lnL +=2 +if opts.input_distance: + print " Distance input" + col_lnL +=1 dat_orig = dat = np.loadtxt(opts.fname) dat_orig = dat[dat[:,col_lnL].argsort()] # sort http://stackoverflow.com/questions/2828059/sorting-arrays-in-numpy-by-column print " Original data size = ", len(dat), dat.shape @@ -994,6 +998,8 @@ for line in dat: if opts.input_tides: P.lambda1 = line[9] P.lambda2 = line[10] + if opts.input_distance: + P.dist = lal.PC_SI*1e6*line[9] # Incompatible with tides, note! # INPUT GRID: Evaluate binary parameters on fitting coordinates line_out = np.zeros(len(coord_names)+2) -- GitLab