Makefile.nonspinning_BNS_gaussian

# Template bank parameters
LOW_FREQUENCY_CUTOFF = 40.0
HIGH_FREQUENCY_CUTOFF = 2048.0
NUM_SPLIT_TEMPLATES = 100
OVERLAP = 10
BANK_PROGRAM = pycbc_geom_nonspinbank
APPROXIMANT = TaylorF2

# controls triggers
IFOS = H1 L1 V1
START = 966384015
STOP = 966484015
TAG = T1200307_LV_gaussian_4s_w_zerolag_injections_40Hz_5ms_test
WEBDIR = ~/public_html/mdc/non_spinning_BNS_gaussian/
NUMBANKS = 4
PEAK = 0
AC_LENGTH = 351
# additional options, e.g.,
#ADDITIONAL_DAG_OPTIONS = "--blind-injections BNS-MDC1-WIDE.xml"

# Injections
# The seed is the string before the suffix _injections.xml
# Change as appropriate, whitespace is important
INJECTIONS := BNS-SpinMDC-ISOTROPIC.xml BNS-SpinMDC-ALIGNED.xml
comma:=,
INJECTION_REGEX = $(subst $(space),$(comma),$(INJECTIONS))

# Segment and frame type info
SEG_SERVER=https://segdb.ligo.caltech.edu
LIGO_FRAME_TYPE='$*_T1200307_V4_EARLY_GAUSSIAN'
LIGO_SEGMENTS="$*:CBC-MDC1_SCIENCE_EARLY_GAUSSIAN"
VIRGO_FRAME_TYPE='V1_T1300121_V1_EARLY_GAUSSIAN'
CHANNEL_NAMES:=--channel-name=H1=GAUSSIAN --channel-name=L1=GAUSSIAN --channel-name=V1=GAUSSIAN

# Get some basic definitions
include Makefile.offline_analysis_rules

#
# Workflow
#

all : dag

BNS_NonSpin_30Hz_earlyaLIGO.xml:
    gsiscp sugar-dev1.phy.syr.edu:/home/jveitch/public_html/bns/mdc/spin/tmpltbank/BNS_NonSpin_30Hz_earlyaLIGO.xml .

$(INJECTIONS):
    gsiscp sugar-dev1.phy.syr.edu:/home/jveitch/public_html/bns/mdc/spin/"{$(INJECTION_REGEX)}" .

%_split_bank.cache : BNS_NonSpin_30Hz_earlyaLIGO.xml
    mkdir -p $*_split_bank
    gstlal_bank_splitter --output-path $*_split_bank --approximant $(APPROXIMANT) --bank-program $(BANK_PROGRAM) --output-cache $@ --overlap $(OVERLAP) --instrument $* --n $(NUM_SPLIT_TEMPLATES) --sort-by mchirp --add-f-final --max-f-final $(HIGH_FREQUENCY_CUTOFF) $<

plots :
    mkdir plots

$(WEBDIR) : 
    mkdir -p $(WEBDIR)

tisi.xml :
    ligolw_tisi --instrument=H1=0:0:0 --instrument=H2=0:0:0 --instrument=L1=0:0:0 --instrument=V1=0:0:0 tisi0.xml
    ligolw_tisi --instrument=H1=0:0:0 --instrument=H2=0:0:0 --instrument=L1=3.14159:3.14159:3.14159 --instrument=V1=7.892:7.892:7.892 tisi1.xml
    ligolw_add --output $@ tisi0.xml tisi1.xml

dag : segments.xml frame.cache tisi.xml plots $(WEBDIR) $(INJECTIONS) $(BANK_CACHE_FILES)
    gstlal_inspiral_pipe --data-source frames --gps-start-time $(START) --gps-end-time $(STOP) --frame-cache frame.cache --frame-segments-file segments.xml --frame-segments-name datasegments  --control-peak-time $(PEAK) --num-banks $(NUMBANKS) --fir-stride 4 --web-dir $(WEBDIR) --time-slide-file tisi.xml $(INJECTION_LIST) --bank-cache $(BANK_CACHE_STRING) --tolerance 0.9999 --overlap $(OVERLAP) --flow $(LOW_FREQUENCY_CUTOFF) $(CHANNEL_NAMES) --autocorrelation-length $(AC_LENGTH) $(ADDITIONAL_DAG_OPTIONS)

%_segmentspadded.xml:
    ligolw_segment_query --segment-url=${SEG_SERVER} -q --gps-start-time ${START} --gps-end-time ${STOP} --include-segments=$(LIGO_SEGMENTS) --result-name=datasegments > $@

V1_frame.cache:
    ligo_data_find -o V -t $(VIRGO_FRAME_TYPE) -l  -s $(START) -e $(STOP) --url-type file > $@

%_frame.cache:
    #FIXME horrible hack to get the observatory, not guaranteed to work
    $(eval OBS:=$*)
    $(eval OBS:=$(subst 1,$(empty),$(OBS)))
    $(eval OBS:=$(subst 2,$(empty),$(OBS)))
    ligo_data_find -o $(OBS) -t $(LIGO_FRAME_TYPE) -l  -s $(START) -e $(STOP) --url-type file > $@

frame.cache: $(FRAME_CACHE_FILES)
    cat $(FRAME_CACHE_FILES) > frame.cache

segments.xml: $(SEGMENTS_FILES) frame.cache
    ligolw_add --output segdb.xml $(SEGMENTS_FILES)
    ligolw_cut --delete-column segment:segment_def_cdb --delete-column segment:creator_db --delete-column segment_definer:insertion_time segdb.xml
    gstlal_cache_to_segments frame.cache nogaps.xml
    gstlal_segments_operations --segment-file1 segdb.xml --segment-file2 nogaps.xml --intersection --output-file $@
    -rm -vf nogaps.xml segdb.xml
    gstlal_segments_trim --trim 8 --min-length 512 --output $@ $@


clean:
    -rm -rvf *.sub *.dag* *.cache *.sh logs *.sqlite plots *.html Images *.css *.js
    -rm -rvf lalapps_run_sqlite/ ligolw_* gstlal_*
    -rm -vf segments.xml tisi.xml H*.xml L*.xml V*.xml ?_injections.xml ????-*_split_bank-*.xml
    -rm -vf *marginalized*.xml.gz *-ALL_LLOID*.xml.gz
    -rm -vf tisi0.xml tisi1.xml
    -rm -rf *_split_bank
    -rm -vf BNS*.xml