diff --git a/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs_dag b/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs_dag index eee4e17c519f431f004da3202dcc608dab60c603..3dcf5ef1c84393b9f5edacdefeceb7bfe109ce37 100755 --- a/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs_dag +++ b/gstlal-inspiral/bin/gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs_dag @@ -25,8 +25,13 @@ parser.add_argument('--psd-xml', help = 'XML containing HLV psd') parser.add_argument('--H-snr', type = float, help = 'H characteristic SNR') parser.add_argument('--L-snr', type = float, help = 'L characteristic SNR') parser.add_argument('--V-snr', type = float, help = 'V characteristic SNR') +parser.add_argument('--flow', type = float, default = 10., help = 'Low frequency cut-off. Default 10 Hz') +parser.add_argument('--fhigh', type = float, default = 1024., help = 'High frequency cut-off. Default 1024 Hz') args = parser.parse_args() +if args.flow >= args.fhigh: + raise ValueError("flow cannot be greater than fhigh") + try: os.mkdir("logs") except: @@ -37,7 +42,7 @@ covJob = dagparts.DAGJob("gstlal_inspiral_compute_dtdphideff_cov_matrix", condor margJob = dagparts.DAGJob("gstlal_inspiral_create_dt_dphi_snr_ratio_pdfs", condor_commands = {"request_memory":"7GB", "want_graceful_removal":"True", "kill_sig":"15", "accounting_group":"ligo.prod.o3.cbc.uber.gstlaloffline"}) addJob = dagparts.DAGJob("gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs", condor_commands = {"request_memory":"4GB", "want_graceful_removal":"True", "kill_sig":"15", "accounting_group":"ligo.prod.o3.cbc.uber.gstlaloffline"}) -covnode = dagparts.DAGNode(covJob, dag, parent_nodes = [], opts = {"H-snr": args.H_snr, "L-snr": args.L_snr, "V-snr": args.V_snr}, input_files = {"psd-xml": args.psd_xml}, output_files = {"output":"covmat.h5"}) +covnode = dagparts.DAGNode(covJob, dag, parent_nodes = [], opts = {"H-snr": args.H_snr, "L-snr": args.L_snr, "V-snr": args.V_snr, "flow": args.flow, "fhigh": args.fhigh}, input_files = {"psd-xml": args.psd_xml}, output_files = {"output":"covmat.h5"}) num = 1000 margnodes = []