2 # Template bank parameters
5 # The filtering start frequency
6 LOW_FREQUENCY_CUTOFF = 30.0
7 # The maximum frequency to filter to
8 HIGH_FREQUENCY_CUTOFF = 1024.0
9 # Controls the number of templates in each SVD sub bank
10 NUM_SPLIT_TEMPLATES = 100
11 # Controls the overlap from sub bank to sub bank - helps mitigate edge effects
12 # in the SVD. Redundant templates will be removed
14 # The program used to make the template bank. This will be searched for in the
15 # process param table in order to extract some metadata
16 BANK_PROGRAM = pycbc_geom_nonspinbank
17 # The approximant that you wish to filter with
18 APPROXIMANT = TaylorF2
21 # Triggering parameters
24 # The detectors to analyze
30 # A user tag for the run
32 # A web directory for output
33 WEBDIR = ~/public_html/MDC/NSBH/Summer2014/recolored/nonspin/$(START)-$(STOP)-$(TAG)
34 # The number of sub banks to process in parallel
for each gstlal_inspiral job
36 # The control peak time
for the composite detection statistic. If
set to 0 the
37 # statistic is disabled
39 # The length of autocorrelation chi-squared in sample points
41 # The minimum number of samples to include in a given time slice
43 # The maximum number of samples to include in the 256 Hz or above time slices
47 # additional options, e.g.,
50 #ADDITIONAL_DAG_OPTIONS =
"--blind-injections BNS-MDC1-WIDE.xml"
56 # The seed is the
string before the suffix _injections.xml
57 # Change as appropriate, whitespace is important
58 INJECTIONS := HL-INJECTIONS_976731_NSBHMDCINJ_SET1_T4-966384015-5184000.xml HL-INJECTIONS_976765_NSBHMDCINJ_SET1_T2-966384015-5184000.xml
59 FAR_INJECTIONS := HL-FAR_INJECTIONS_976731_NSBHMDCINJ_SET1_T4-966384015-5184000.xml HL-FAR_INJECTIONS_976765_NSBHMDCINJ_SET1_T2-966384015-5184000.xml
60 # NOTE you shouldn
't need to change these next three lines
62 INJECTION_REGEX = $(subst $(space),$(comma),$(INJECTIONS))
63 FAR_INJECTION_REGEX = $(subst $(space),$(comma),$(FAR_INJECTIONS))
66 # Segment and frame type info
69 # The LIGO and Virgo frame types
70 LIGO_FRAME_TYPE='T1200307_V4_EARLY_RECOLORED_V2
'
71 VIRGO_FRAME_TYPE='T1300121_V1_EARLY_RECOLORED_V2
'
72 # The Channel names. FIXME sadly you have to change the CHANNEL_NAMES string if
73 # you want to analyze a different set of IFOS
74 H1_CHANNEL=LDAS-STRAIN
75 L1_CHANNEL=LDAS-STRAIN
77 CHANNEL_NAMES:=--channel-name=H1=$(H1_CHANNEL) --channel-name=L1=$(L1_CHANNEL)
80 # Get some basic definitions. NOTE this comes from the share directory probably.
83 include Makefile.offline_analysis_rules
92 gsiscp h2.atlas.aei.uni-hannover.de:/home/tito/projects/spin_search/nsbh_mdc1/bank_no_spin/bank_no_spin.xml.gz .
95 gsiscp ldas-pcdev1.ligo-wa.caltech.edu:/home/spxiwh/aLIGO/MDC/MDC1/inj_sets/inj_set_1/"{$(INJECTION_REGEX)}" .
98 gsiscp ldas-pcdev1.ligo-wa.caltech.edu:/home/spxiwh/aLIGO/MDC/MDC1/inj_sets/inj_set_1/"{$(FAR_INJECTION_REGEX)}" .
100 %_split_bank.cache : bank_no_spin.xml.gz
101 mkdir -p $*_split_bank
102 gstlal_bank_splitter --f-low $(LOW_FREQUENCY_CUTOFF) --group-by-chi --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) $<
111 ligolw_tisi --instrument=H1=0:0:0 --instrument=H2=0:0:0 --instrument=L1=0:0:0 --instrument=V1=0:0:0 tisi0.xml
112 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
113 ligolw_add --output $@ tisi0.xml tisi1.xml
115 dag : segments.xml.gz vetoes.xml.gz frame.cache tisi.xml plots $(WEBDIR) $(INJECTIONS) $(BANK_CACHE_FILES) $(FAR_INJECTIONS)
116 gstlal_inspiral_pipe --data-source frames --gps-start-time $(START) --gps-end-time $(STOP) --frame-cache frame.cache --frame-segments-file segments.xml.gz --vetoes vetoes.xml.gz --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) --samples-min $(SAMPLES_MIN) --samples-max-256 $(SAMPLES_MAX_256) $(ADDITIONAL_DAG_OPTIONS) $(FAR_INJECTION_LIST)
119 # FIXME force the observatory column to actually be instrument
120 ligo_data_find -o V -t $(VIRGO_FRAME_TYPE) -l -s $(START) -e $(STOP) --url-type file | awk '{ print $$1
" V1_"$$2
" "$$3
" "$$4
" "$$5}
' > $@
123 # FIXME horrible hack to get the observatory, not guaranteed to work
125 $(eval OBS:=$(subst 1,$(empty),$(OBS)))
126 $(eval OBS:=$(subst 2,$(empty),$(OBS)))
127 # FIXME force the observatory column to actually be instrument
128 ligo_data_find -o $(OBS) -t $(LIGO_FRAME_TYPE) -l -s $(START) -e $(STOP) --url-type file | awk '{ print $$1
" $*_"$$2
" "$$3
" "$$4
" "$$5}
' > $@
130 frame.cache: $(FRAME_CACHE_FILES)
131 cat $(FRAME_CACHE_FILES) > frame.cache
133 segments.xml.gz: frame.cache
134 # These segments come from the MDC set
135 gsiscp pcdev3.cgca.uwm.edu:/home/channa/public_html/SELECTED_SEGS.xml.gz $@
136 gstlal_cache_to_segments frame.cache nogaps.xml
137 gstlal_segments_operations --segment-file1 $@ --segment-file2 nogaps.xml --intersection --output-file $@
139 gstlal_segments_trim --trim 8 --gps-start-time $(START) --gps-end-time $(STOP) --min-length 2048 --output $@ $@
142 gsiscp pcdev3.cgca.uwm.edu:/home/channa/public_html/COMBINED_CAT_4_VETO_SEGS.xml.gz $@
143 gstlal_segments_trim --gps-start-time $(START) --gps-end-time $(STOP) --segment-name vetoes --output $@ $@
146 -rm -rvf *.sub *.dag* *.cache *.sh logs *.sqlite plots *.html Images *.css *.js
147 -rm -rvf lalapps_run_sqlite/ ligolw_* gstlal_*
148 -rm -vf segments.xml.gz tisi.xml H*.xml L*.xml V*.xml ?_injections.xml ????-*_split_bank-*.xml vetoes.xml.gz
149 -rm -vf *marginalized*.xml.gz *-ALL_LLOID*.xml.gz
150 -rm -vf tisi0.xml tisi1.xml
152 -rm -rf bank_no_spin.xml.gz
153 -rm -rf $(INJECTIONS) $(FAR_INJECTIONS)