testHierarchicalSearch.py

Go to the documentation of this file.

# Copyright (C) 2022 Rodrigo Tenorio, David Keitel
#
# 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.
 
import subprocess
 
# run this test with "make check TESTS=testHierarchicalSearch"
# everything will be automatically put into testHierarchicalSearch.testdir
testdir = "."
 
mfd_cml = " ".join(
    [
        "lalpulsar_Makefakedata_v5",
        "--outSingleSFT=TRUE",
        f"--outSFTdir={testdir}",
        '--outLabel="simulatedsignal"',
        '--IFOs="H1","L1"',
        '--sqrtSX="1e-22"',
        "--startTime=1000000000",
        f"--duration={10 * 1800}",
        "--fmin=100",
        "--Band=1",
        "--Tsft=1800",
    ]
)
 
hs_cml = " ".join(
    [
        "lalpulsar_HierarchicalSearch",
        "--ephemSun=sun00-40-DE405.dat.gz",
        "--ephemEarth=earth00-40-DE405.dat.gz",
        "--method=0",
        "--Freq=100.5",
        "--dFreq=0.001",
        "--FreqBand=0.01",
        '--skyRegion="allsky"',
        f'--DataFiles1="{testdir}/*-10*.sft"',
        f"--tStack={1800 * 5}",
        "--printCand1",
        f"--fnameout=HS_test",
    ]
)
 
print(f"Running MFD to generate SFTs: {mfd_cml}")
subprocess.check_call(mfd_cml, shell=True)
print(f"Running test: {hs_cml}")
subprocess.check_call(hs_cml, shell=True)