Skip to content
Snippets Groups Projects

Plot horizon distance from ranking statistics

Merged ChiWai Chan requested to merge plot_psd_horizon into master
1 unresolved thread
3 files
+ 17
15
Compare changes
  • Side-by-side
  • Inline
Files
3
  • gstlal_inspiral_pipe and inspiral_pipe.py: replace gstlal_insprial_plot_psd_horizon with gstlal_inspiral_plot_rankingstats_horizon, and temporarily remove horizon distance from summary page since there will be too many horizon distance plots.
@@ -44,6 +44,7 @@ __author__ = 'Chad Hanna <chad.hanna@ligo.org>, Patrick Godwin <patrick.godwin@l
from collections import OrderedDict
from optparse import OptionParser
from pathlib import Path
import os
import warnings
@@ -301,7 +302,7 @@ def set_up_jobs(options):
jobs['svd'] = dagparts.DAGJob("gstlal_svd_bank", condor_commands = svd_condor_commands)
jobs['model'] = dagparts.DAGJob("gstlal_inspiral_mass_model", condor_commands = base_condor_commands)
jobs['modelAdd'] = dagparts.DAGJob("gstlal_inspiral_add_mass_models", condor_commands = base_condor_commands)
jobs['horizon'] = dagparts.DAGJob("gstlal_plot_psd_horizon", condor_commands = base_condor_commands)
jobs['horizon'] = dagparts.DAGJob("gstlal_inspiral_plot_rankingstats_horizon", condor_commands = base_condor_commands)
jobs['gstlalInspiral1IFO'] = dagparts.DAGJob("gstlal_inspiral", tag_base="gstlal_inspiral_1ifo", condor_commands = inspiral_1ifo_condor_commands)
jobs['gstlalInspiral2IFO'] = dagparts.DAGJob("gstlal_inspiral", tag_base="gstlal_inspiral_2ifo", condor_commands = inspiral_2ifo_condor_commands)
jobs['gstlalInspiral3IFO'] = dagparts.DAGJob("gstlal_inspiral", tag_base="gstlal_inspiral_3ifo", condor_commands = inspiral_3ifo_condor_commands)
@@ -362,6 +363,7 @@ if __name__ == '__main__':
# output directories
output_dir = "plots"
try:
os.mkdir(Path(output_dir) / "horizon_distance")
os.mkdir("logs")
except:
pass
@@ -406,9 +408,6 @@ if __name__ == '__main__':
# Compute the PSDs for each segment
psd_nodes = inspiral_pipe.ref_psd_layer(dag, jobs, [], segsdict, channel_dict, options)
# plot the horizon distance
inspiral_pipe.horizon_dist_layer(dag, jobs, psd_nodes, options, boundary_seg, output_dir, instruments)
# compute the median PSD
median_psd_node = inspiral_pipe.median_psd_layer(dag, jobs, psd_nodes, options, boundary_seg, instruments)
ref_psd = median_psd_node.output_files["output-name"]
@@ -437,6 +436,9 @@ if __name__ == '__main__':
# marginalize jobs
marg_nodes = inspiral_pipe.marginalize_layer(dag, jobs, svd_nodes, lloid_output, lloid_diststats, options, boundary_seg, instrument_set, model_node, model_file, ref_psd, svd_dtdphi_map, options.idq_file)
# plot the horizon distance from ranking stats
inspiral_pipe.horizon_dist_layer(dag, jobs, marg_nodes, Path(output_dir) / "horizon_history", instruments)
# calc rank PDF jobs
rankpdf_nodes, rankpdf_zerolag_nodes = inspiral_pipe.calc_rank_pdf_layer(dag, jobs, marg_nodes, options, boundary_seg, instrument_set)
Loading