lalsuite 7.26.2.dev20251206__cp312-cp312-manylinux_2_28_x86_64.whl

lalsuite-7.26.2.dev20251206.data/scripts/lalburst_power_plot_binjtf

@@ -0,0 +1,302 @@
+#!python
+#
+# Copyright (C) 2007  Kipp Cannon
+#
+# This program is free software; you can redistribute it and/or modify it
+# under the terms of the GNU General Public License as published by the
+# Free Software Foundation; either version 2 of the License, or (at your
+# option) any later version.
+#
+# This program is distributed in the hope that it will be useful, but
+# WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General
+# Public License for more details.
+#
+# You should have received a copy of the GNU General Public License along
+# with this program; if not, write to the Free Software Foundation, Inc.,
+# 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
+
+
+#
+# =============================================================================
+#
+#                                   Preamble
+#
+# =============================================================================
+#
+
+
+import math
+from matplotlib import cm, colors, collections
+import numpy
+from optparse import OptionParser
+import sqlite3
+import sys
+
+
+from lalburst import git_version
+from lalburst import SnglBurstUtils
+
+
+__author__ = "Kipp Cannon <kipp.cannon@ligo.org>"
+__version__ = "git id %s" % git_version.id
+__date__ = git_version.date
+
+
+#
+# =============================================================================
+#
+#                                 Command Line
+#
+# =============================================================================
+#
+
+
+def parse_command_line():
+	parser = OptionParser(
+		version = "Name: %%prog\n%s" % git_version.verbose_msg
+	)
+	parser.add_option("-b", "--base", metavar = "base", default = "plotbinjtf_", help = "set the prefix for output filenames (default = plotbinj_)")
+	parser.add_option("-f", "--format", metavar = "format", default = "png", help = "set the output image format (default = png)")
+	parser.add_option("-v", "--verbose", action = "store_true", help = "be verbose")
+	options, filenames = parser.parse_args()
+
+	return options, (filenames or [None])
+
+options, filenames = parse_command_line()
+
+
+#
+# =============================================================================
+#
+#                             Time-Frequency Plane
+#
+# =============================================================================
+#
+
+
+def time_freq_plot(database, instrument, sim):
+	fig = SnglBurstUtils.figure.Figure()
+	SnglBurstUtils.FigureCanvas(fig)
+	# 6.5" wide, golden ratio high
+	fig.set_size_inches(6.5, 6.5 / ((1 + math.sqrt(5)) / 2))
+
+	#
+	# top plot --- triggers that matched injection
+	#
+
+	t_sim = sim.time_at_instrument(instrument)
+
+	axes = fig.add_subplot(211)
+	axes.grid(True)
+	#axes.set_xlabel("$t - t_{\\mathrm{injection}}$ (s)")
+	axes.set_ylabel("$f - f_{\\mathrm{injection}}$ (Hz)")
+	axes.set_title("%s Triggers Matching %g Hz Injection at GPS %s" % (instrument, sim.frequency, t_sim))
+
+	xmin = xmax = 0.0
+	ymin = ymax = 0.0
+	verts = []
+	colours = []
+	peakx = []
+	peaky = []
+	match_ids = []
+	# find triggers from the desired instrument that are coincident
+	# with the injection, and iterate over them in order from least to
+	# most confident
+	for burst in map(database.sngl_burst_table.row_from_cols, database.connection.cursor().execute("""
+SELECT sngl_burst.* FROM
+	sngl_burst
+	JOIN coinc_event_map AS a ON (
+		sngl_burst.event_id == a.event_id
+		AND a.table_name == 'sngl_burst'
+	)
+	JOIN coinc_event_map AS b ON (
+		a.coinc_event_id == b.coinc_event_id
+		AND b.table_name == 'sim_burst'
+	)
+WHERE
+	sngl_burst.ifo == ?
+	AND b.event_id == ?
+ORDER BY
+	sngl_burst.confidence ASC
+	""", (instrument, sim.simulation_id))):
+		match_ids.append(burst.event_id)
+
+		# Add time-frequency tile to collection
+		tmin = float(burst.start - t_sim)
+		tmax = float(burst.start + burst.duration - t_sim)
+		fmin = burst.central_freq - burst.bandwidth / 2 - sim.frequency
+		fmax = burst.central_freq + burst.bandwidth / 2 - sim.frequency
+		verts.append(((tmin, fmin), (tmax, fmin), (tmax, fmax), (tmin, fmax)))
+		colours.append(burst.confidence)
+
+		try:
+			# draw most significant tile if there is one
+			tmin = float(burst.ms_start - t_sim)
+			tmax = float(burst.ms_start + burst.ms_duration - t_sim)
+			fmin = burst.ms_flow - sim.frequency
+			fmax = burst.ms_flow + burst.ms_bandwidth - sim.frequency
+			verts.append(((tmin, fmin), (tmax, fmin), (tmax, fmax), (tmin, fmax)))
+			colours.append(burst.ms_confidence)
+		except AttributeError:
+			pass
+
+		peakx.append(float(burst.peak - t_sim))
+		try:
+			# use peak_frequency col if it exists
+			peaky.append(burst.peak_frequency - sim.frequency)
+		except AttributeError:
+			peaky.append(burst.central_freq - sim.frequency)
+
+		# update bounding box
+		tmin = float(burst.start - t_sim)
+		tmax = float(burst.start + burst.duration - t_sim)
+		fmin = burst.central_freq - burst.bandwidth / 2 - sim.frequency
+		fmax = burst.central_freq + burst.bandwidth / 2 - sim.frequency
+		xmin = min(xmin, tmin)
+		xmax = max(xmax, tmax)
+		ymin = min(ymin, fmin)
+		ymax = max(ymax, fmax)
+
+	polys = collections.PolyCollection(verts)
+	polys.set_array(numpy.array(colours))
+	polys.set_alpha(0.3)
+	polys.set_cmap(cm.get_cmap())
+	polys.set_norm(colors.normalize())
+	axes.add_collection(polys)
+
+	axes.plot(peakx, peaky, "k+")
+
+	axes.axvline(0, color = "k")
+	axes.axhline(0, color = "k")
+
+	# set the bounding box
+	axes.set_xlim([1.4 * xmin, 1.4 * xmax])
+	axes.set_ylim([1.4 * ymin, 1.4 * ymax])
+
+	#
+	# bottom plot --- triggers near injection
+	#
+
+	axes = fig.add_subplot(212)
+	axes.grid(True)
+	axes.set_xlabel("$t - t_{\\mathrm{injection}}$ (s)")
+	axes.set_ylabel("$f - f_{\\mathrm{injection}}$ (Hz)")
+
+	xmin = xmax = 0.0
+	ymin = ymax = 0.0
+	verts = []
+	colours = []
+	edgecolours = []
+	peakx = []
+	peaky = []
+	# find triggers from the desired instrument that are near the
+	# injection, and iterate over them in order from least to most
+	# confident
+	for burst in map(database.sngl_burst_table.row_from_cols, database.connection.cursor().execute("""
+SELECT * FROM
+	sngl_burst
+WHERE
+	ifo == ?
+	AND start_time BETWEEN ? AND ?
+	AND central_freq BETWEEN ? AND ?
+ORDER BY
+	sngl_burst.confidence ASC
+	""", (instrument, int(t_sim - 2), int(t_sim + 2), sim.frequency - 300, sim.frequency + 300))):
+		# Add time-frequency tile to collection
+		tmin = float(burst.start - t_sim)
+		tmax = float(burst.start + burst.duration - t_sim)
+		fmin = burst.central_freq - burst.bandwidth / 2 - sim.frequency
+		fmax = burst.central_freq + burst.bandwidth / 2 - sim.frequency
+		verts.append(((tmin, fmin), (tmax, fmin), (tmax, fmax), (tmin, fmax)))
+		colours.append(burst.confidence)
+		if burst.event_id in match_ids:
+			edgecolours.append("g")
+		else:
+			edgecolours.append("k")
+
+		peakx.append(float(burst.peak - t_sim))
+		try:
+			# use peak_frequency col if it exists
+			peaky.append(burst.peak_frequency - sim.frequency)
+		except:
+			peaky.append(burst.central_freq - sim.frequency)
+
+		# update bounding box
+		xmin = min(xmin, tmin)
+		xmax = max(xmax, tmax)
+		ymin = min(ymin, fmin)
+		ymax = max(ymax, fmax)
+
+	polys = collections.PolyCollection(verts, edgecolors = edgecolours)
+	polys.set_array(numpy.array(colours))
+	polys.set_alpha(0.3)
+	polys.set_cmap(cm.get_cmap())
+	polys.set_norm(colors.normalize())
+	axes.add_collection(polys)
+
+	axes.plot(peakx, peaky, "k+")
+
+	axes.axvline(0, color = "k")
+	axes.axhline(0, color = "k")
+
+	# set the bounding box
+	axes.set_xlim([1.4 * xmin, 1.4 * xmax])
+	axes.set_ylim([1.4 * ymin, 1.4 * ymax])
+
+
+	return fig
+
+
+#
+# =============================================================================
+#
+#                                     Plot
+#
+# =============================================================================
+#
+
+
+def found_injections(contents, instrument):
+	for values in contents.connection.cursor().execute("""
+SELECT
+	*
+FROM
+	sim_burst
+WHERE
+	EXISTS (
+		-- Find a link through the coinc_event_map table to a row
+		-- in the sngl_burst table with the correct ifo value.
+		SELECT
+			*
+		FROM
+			coinc_event_map AS a
+			JOIN coinc_event_map AS b ON (
+				a.coinc_event_id == b.coinc_event_id
+			)
+			JOIN sngl_burst ON (
+				b.table_name == 'sngl_burst'
+				AND b.event_id == sngl_burst.event_id
+			)
+		WHERE
+			a.table_name == 'sim_burst'
+			AND a.event_id == sim_burst.simulation_id
+			AND sngl_burst.ifo == ?
+	)
+	""", (instrument,)):
+		yield contents.sim_burst_table.row_from_cols(values)
+
+
+for n, filename in enumerate(filenames):
+	if options.verbose:
+		print("%d/%d: %s" % (n + 1, len(filenames), filename), file=sys.stderr)
+	database = SnglBurstUtils.CoincDatabase(sqlite3.connect(filename), "lalapps_power")
+	if options.verbose:
+		SnglBurstUtils.summarize_coinc_database(database)
+	for instrument in database.instruments:
+		for sim in found_injections(database, instrument):
+			plotname = "%s%d_%s.%s" % (options.base, sim.time_at_instrument(instrument).seconds, instrument, options.format)
+			if options.verbose:
+				print("--> %s" % plotname, file=sys.stderr)
+			time_freq_plot(database, instrument, sim).savefig(plotname)
+	database.connection.close()

lalsuite-7.26.2.dev20251206.data/scripts/lalburst_version

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalburst_version') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-cat

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-cat') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-cksum

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-cksum') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-cut

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-cut') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-dump

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-dump') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-fmt

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-fmt') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-paste

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-paste') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-print

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-print') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-split

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-split') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-stat

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-stat') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-stream

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-stream') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalfr-vis

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalfr-vis') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalframe_version

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalframe_version') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalinference_bench

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalinference_bench') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalinference_burst

@@ -0,0 +1,6 @@
+#!python
+import os
+import sys
+from importlib.resources import files, as_file
+with as_file(files('lalapps') / 'bin' / 'lalinference_burst') as path:
+    os.execv(path, sys.argv)

lalsuite-7.26.2.dev20251206.data/scripts/lalinference_burst_pp_pipe

@@ -0,0 +1,220 @@
+#!python
+# End-to-end LALInference test pipeline
+# (C) 2014 John Veitch
+
+from __future__ import print_function
+
+import os
+import sys
+import uuid
+from optparse import OptionParser
+
+from configparser import ConfigParser
+
+from lal import pipeline
+
+from lalinference import lalinference_pipe_utils as pipe_utils
+from lalinference.lalinference_pipe_utils import mkdirs
+
+usage=""" %prog [options] config.ini
+Setup a DAG to run an end-to-end lalinference test:
+ 1) Generate samples from prior
+ 2) Analyse a set of injections drawn from the prior
+ 3) Run P vs P test on results
+"""
+
+parser=OptionParser(usage)
+parser.add_option("-r","--run-path",default='./',action="store",type="string",help="Directory to run pipeline in (default: $PWD)",metavar="RUNDIR")
+parser.add_option("-p","--daglog-path",default=None,action="store",type="string",help="Path to directory to contain DAG log file. SHOULD BE LOCAL TO SUBMIT NODE",metavar="LOGDIR")
+parser.add_option("-I","--injections",default=None,action="store",type="string",help="Path to injection file, will bypass the prior-sampling stage",metavar="injections.xml")
+parser.add_option('-N','--trials',action='store',type='int',metavar='NUM',help='Number of prior samples to analyse')
+
+
+(opts,args)=parser.parse_args()
+
+if len(args)==0:
+        parser.print_help()
+        sys.exit(1)
+
+inifile=args[0]
+
+# Set up the configuration for the sub-dags
+
+prior_cp=ConfigParser()
+prior_cp.optionxform = str
+try:
+  prior_cp.read_file(open(inifile))
+except AttributeError:
+  prior_cp.readfp(open(inifile))
+
+main_cp=ConfigParser()
+main_cp.optionxform = str
+try:
+  main_cp.read_file(open(inifile))
+except AttributeError:
+  main_cp.readfp(open(inifile))
+
+
+rundir=os.path.abspath(opts.run_path)
+
+if opts.daglog_path is not None:
+  prior_cp.set('paths','daglogdir',os.path.join(os.path.abspath(opts.daglog_path),'prior'))
+  main_cp.set('paths','daglogdir',os.path.join(os.path.abspath(opts.daglog_path),'main'))
+  daglogdir=os.path.abspath(opts.daglog_path)
+else:
+  prior_cp.set('paths','daglogdir',os.path.join(os.path.abspath(opts.run_path),'prior'))
+  main_cp.set('paths','daglogdir',os.path.join(os.path.abspath(opts.run_path),'main'))
+  daglogdir=os.path.abspath(opts.run_path)
+
+webdir=main_cp.get('ppanalysis','webdir')
+priordir=os.path.join(rundir,'prior')
+maindir=os.path.join(rundir,'main')
+priorwebdir=os.path.join(webdir,'prior')
+mainwebdir=os.path.join(webdir,'injections')
+prior_cp.set('paths','basedir',priordir)
+main_cp.set('paths','basedir',maindir)
+prior_cp.set('paths','webdir',priorwebdir)
+main_cp.set('paths','webdir',mainwebdir)
+
+outerlogdir=os.path.join(daglogdir,'log')
+mkdirs(outerlogdir)
+mkdirs(priordir)
+mkdirs(maindir)
+mkdirs(priorwebdir)
+mkdirs(mainwebdir)
+
+# Add the prior options to the sub dag
+if prior_cp.get('analysis','engine')=='lalinferencenest':
+  prior_cp.set('engine','sampleprior',str(20*opts.trials)) # more samples drawn since not all will end up in posterior
+  prior_cp.set('engine','zeroLogLike','')
+  prior_cp.set('engine','nlive',str(20*opts.trials))
+if prior_cp.get('analysis','engine')=='lalinferenceburst':
+  prior_cp.set('engine','sampleprior',str(20*opts.trials)) # more samples drawn since not all will end up in posterior
+  prior_cp.set('engine','zeroLogLike','')
+  prior_cp.set('engine','nlive',str(20*opts.trials))
+elif prior_cp.get('analysis','engine')=='lalinferencemcmc':
+  prior_cp.set('engine','Neff',str(opts.trials))
+  prior_cp.set('engine','zeroLogLike','')
+elif prior_cp.get('analysis','engine')=='lalinferencebambi':
+  prior_cp.set('engine','zeroLogLike','')
+  prior_cp.set('engine','nlive',str(opts.trials))
+elif prior_cp.get('analysis','engine')=='lalinferencebambimpi':
+  prior_cp.set('engine','zeroLogLike','')
+  prior_cp.set('engine','nlive',str(opts.trials))
+
+# Create a DAG to contain the other scripts
+outerdaglog=os.path.join(daglogdir,'lalinference_injection_test_'+str(uuid.uuid1())+'.log')
+outerdag=pipeline.CondorDAG(outerdaglog)
+outerdag.set_dag_file(os.path.join(rundir,'priortest'))
+
+# Run code with prior sampling
+trig_time=1085855789
+fake_event=pipe_utils.Event(trig_time=trig_time)
+tfpath=os.path.join(rundir,'time.txt')
+tfile=open(tfpath,'w')
+print('%i\n'%(trig_time), file=tfile)
+tfile.close()
+prior_cp.set('input','gps-time-file',tfpath)
+
+priordag=pipe_utils.LALInferencePipelineDAG(prior_cp)
+priordag.set_dag_file(os.path.join(priordir,'lalinference_priorsample'))
+priordagjob=pipeline.CondorDAGManJob(priordag.get_dag_file(),dir=priordir)
+priordagnode=pipeline.CondorDAGManNode(priordagjob)
+# Find the output file
+pagenode=filter(lambda n:isinstance(n,pipe_utils.ResultsPageNode), priordag.get_nodes())[0]
+priorfile=pagenode.get_pos_file()
+
+# Convert prior samples to injections
+convertsub=os.path.join(rundir,'samples2injections.sub')
+converterr=os.path.join(outerlogdir,'samples2injection-$(cluster)-$(process)-$(node).err')
+convertout=os.path.join(outerlogdir,'samples2injection-$(cluster)-$(process)-$(node).out')
+
+if opts.injections:
+        injfile=os.path.abspath(opts.injections)
+else:
+        injfile=os.path.join(rundir,'priorsamples.xml')
+approx=prior_cp.get('engine','approx')
+prior2injexe=prior_cp.get('condor','pos_to_sim_burst')
+prior2injjob=pipeline.CondorDAGJob('vanilla',prior2injexe)
+if main_cp.has_option('analysis','accounting_group'):
+  prior2injjob.add_condor_cmd('accounting_group',main_cp.get('analysis','accounting_group'))
+prior2injjob.set_sub_file(convertsub)
+prior2injjob.set_stderr_file(converterr)
+prior2injjob.set_stdout_file(convertout)
+prior2injjob.add_condor_cmd('getenv','True')
+prior2injnode=pipeline.CondorDAGNode(prior2injjob)
+prior2injnode.add_var_opt('output',injfile)
+prior2injnode.add_var_opt('num-of-injs',str(opts.trials))
+prior2injnode.add_var_opt('approx',approx)
+prior2injnode.add_var_arg(priorfile)
+prior2injnode.add_parent(priordagnode)
+
+# Create the pipeline based on the injections
+#main_cp.set('input','injection-file',injfile)
+main_cp.set('input','gps-start-time',str(trig_time-1000))
+main_cp.set('input','gps-end-time',str(trig_time+1000))
+maindag=pipe_utils.LALInferencePipelineDAG(main_cp)
+maindag.set_dag_file(os.path.join(maindir,'lalinference_pipeline'))
+maindagjob=pipeline.CondorDAGManJob(maindag.get_dag_file(),dir=maindir)
+maindagnode=pipeline.CondorDAGManNode(maindagjob)
+maindag.config.set('input','burst-injection-file',injfile)
+for i in range(int(opts.trials)):
+  ev=pipe_utils.Event(trig_time=trig_time,event_id=i)
+  e=maindag.add_full_analysis(ev)
+
+outerdag.add_node(maindagnode)
+
+if not opts.injections:
+    outerdag.add_node(priordagnode)
+    outerdag.add_node(prior2injnode)
+    maindagnode.add_parent(prior2injnode)
+
+# Get a list of posterior samples files
+resultspagenodes=filter(lambda n: isinstance(n, pipe_utils.ResultsPageNode), maindag.get_nodes())
+posteriorfiles=[n.get_pos_file() for n in resultspagenodes]
+
+# Analyse results of injection runs to generate PP plot
+ppsub=os.path.join(rundir,'ppanalysis.sub')
+pperr=os.path.join(outerlogdir,'ppanalysis-$(cluster)-$(process)-$(node).err')
+ppout=os.path.join(outerlogdir,'ppanalysis-$(cluster)-$(process)-$(node).out')
+ppexe=prior_cp.get('condor','ppanalysis')
+ppjob=pipeline.CondorDAGJob('vanilla',ppexe)
+ppjob.set_sub_file(ppsub)
+ppjob.set_stderr_file(pperr)
+ppjob.set_stdout_file(ppout)
+ppjob.add_condor_cmd('getenv','True')
+if main_cp.has_option('analysis','accounting_group'):
+  ppjob.add_condor_cmd('accounting_group',main_cp.get('analysis','accounting_group'))
+
+ppnode=pipeline.CondorDAGNode(ppjob)
+ppnode.add_var_opt('injXML',injfile)
+if main_cp.has_option('ppanalysis','webdir'):
+  outdir=main_cp.get('ppanalysis','webdir')
+else:
+  outdir=os.path.join(rundir,'ppanalysis')
+
+mkdirs(outdir)
+ppnode.add_var_opt('outdir',outdir)
+for f in posteriorfiles:
+  ppnode.add_var_arg(f)
+
+ppnode.add_parent(maindagnode)
+outerdag.add_node(ppnode)
+
+
+outerdag.write_sub_files()
+outerdag.write_dag()
+outerdag.write_script()
+
+if not opts.injections:
+    priordag.write_sub_files()
+    priordag.write_dag()
+    priordag.write_script()
+
+maindag.write_sub_files()
+maindag.write_dag()
+maindag.write_script()
+
+# End of program
+print('Successfully created DAG file.')
+print('Now run condor_submit_dag %s\n'%(outerdag.get_dag_file()))