#!/home/conda/feedstock_root/build_artifacts/bld/rattler-build_gstlal-ugly_1767612080/host_env_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_placehold_pl/bin/python

import numpy
import sys
import glob
import lal
from ligo.lw import utils as ligolw_utils
from gstlal.psd import read_psd, write_psd
from optparse import OptionParser

parser = OptionParser(description = __doc__)
parser.add_option("--instrument", help = "The instrument name (required).")
parser.add_option("--output-name", metavar = "filename", help = "The output xml file (required)")
parser.add_option("--verbose", action = "store_true", help = "Be verbose.")
options, filenames = parser.parse_args()


# Identify the various files

psd = read_psd(filenames[0], verbose=options.verbose)[options.instrument]

# FIXME assumes all psds have the same frequency vectors !!!
psd.data = numpy.median(numpy.array([read_psd(f, verbose=options.verbose)[options.instrument].data for f in filenames]), axis=0)

# Write it to disk
write_psd(options.output_name, {options.instrument: psd}, verbose=options.verbose)
