Skip to content
Snippets Groups Projects
Commit d94b23a9 authored by Chad Hanna's avatar Chad Hanna
Browse files

treebank: update makefile for simulations

parent d6751498
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
gstlal-ugly/share/Makefile.treebank
+ 70
− 60
View file @ d94b23a9
# SMALL TEST CASE
#M1bi=1.02
M1bi=1
M1bf=2
M2bi=1
M2bf=2
S1bi=-0.0
S1bf=0.0
S2bi=-0.0
S2bf=0.0
M1bi=1.4
M1bf=2.1
M2bi=1.4
M2bf=2.1
S1bi=0
S1bf=0
S2bi=0
S2bf=0
FLOW = 30
MM=0.97
SBMM=0.97
M1i=5.0
M1f=10
M2i=5.0
M2f=10
M1i=1.6
M1f=1.9
M2i=1.6
M2f=1.9
#S1i = 0
#S1f = .85
S1i = 0
......@@ -24,87 +24,97 @@ S1f = 0
#all : banksim_match_vs_injm1.png
all: BBH1_treebank.xml.gz BBH2_treebank.xml.gz BNS_treebank.xml.gz NSBH1_treebank.xml.gz NSBH2_treebank.xml.gz
all: BNS_treebank.xml.gz NSBH1_treebank.xml.gz NSBH2_treebank.xml.gz BBH1_treebank.xml.gz BBH2_treebank.xml.gz
ligolw_add --output spin_treebank.xml.gz $^ --verbose
###
### FIRST PART OF BBH REGION
###
BBH1_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 30.00 --max-mass2 30.00 --min-mass2 2.00 --min-mass1 2.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.985 --max-spin2z 0.985 --output-name $@ --max-mtotal 400 --max-q 20
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 30.00 --max-mass2 30.00 --min-mass2 2.00 --min-mass1 2.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.985 --max-spin2z 0.985 --output-name $@ --max-mtotal 100 --max-q 20
###
### SECOND PART OF BBH REGION
###
BBH2_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 399.0 --max-mass2 399.0 --min-mass2 2.0 --min-mass1 30.0 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.985 --max-spin2z 0.985 --output-name $@ --max-mtotal 400 --max-q 20
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 99.0 --max-mass2 99.0 --min-mass2 2.0 --min-mass1 30.0 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.985 --max-spin2z 0.985 --output-name $@ --max-mtotal 100 --max-q 20
###
### BNS REGION
###
BNS_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant TaylorF2 --max-mass1 2.00 --max-mass2 2.00 --min-mass2 1.00 --min-mass1 1.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.4 --max-spin1z 0.4 --min-spin2z -0.4 --max-spin2z 0.4 --output-name $@ --max-mtotal 400 --max-q 20
gstlal_inspiral_treebank --verbose --approximant TaylorF2 --max-mass1 2.00 --max-mass2 2.00 --min-mass2 1.00 --min-mass1 1.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.0 --max-spin1z 0.0 --min-spin2z -0.0 --max-spin2z 0.0 --output-name $@ --max-mtotal 100 --max-q 20
NSBH1_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 15.0 --max-mass2 2.00 --min-mass2 1.00 --min-mass1 2.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.0 --max-spin2z 0.0 --output-name $@ --max-mtotal 400 --max-q 20
BNS_sbank.xml.gz:
lalapps_cbc_sbank --verbose --approximant IMRPhenomD --mass1-min 1 --mass1-max 2 --mass2-min 1 --mass2-max 2 --flow $(FLOW) --reference-psd psd_for_treebank.xml.gz --instrument H1 --spin1-min 0 --spin1-max 0 --checkpoint 50 --iterative-match-df-max 2.0 --coarse-match-df 2.0 --aligned-spin --output-filename $@ --match-min $(SBMM) --fhigh-max 1024 --spin2-max 0 --spin2-min 0
NSBH2_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 40.0 --max-mass2 2.00 --min-mass2 1.00 --min-mass1 15.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.0 --max-spin2z 0.0 --output-name $@ --max-mtotal 400 --max-q 20
nonspin_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 399.0 --max-mass2 399.0 --min-mass2 1.0 --min-mass1 1.0 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.0 --max-spin1z 0.0 --min-spin2z -0.0 --max-spin2z 0.0 --output-name $@ --max-mtotal 400 --max-q 20
###
### FIRST PART OF NSBH REGION
###
test_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 $(M1bf) --max-mass2 $(M2bf) --min-mass2 $(M2bi) --min-mass1 $(M1bi) --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z $(S1bi) --max-spin1z $(S1bf) --min-spin2z $(S2bi) --max-spin2z $(S2bf) --output-name $@ --max-mtotal 40 --max-q 4
NSBH1_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 15.0 --max-mass2 2.00 --min-mass2 1.00 --min-mass1 2.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.0 --max-spin2z 0.0 --output-name $@ --max-mtotal 100 --max-q 20
###
### SECOND PART OF NSBH REGION
###
NSBH2_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 40.0 --max-mass2 2.00 --min-mass2 1.00 --min-mass1 15.00 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.985 --max-spin1z 0.985 --min-spin2z -0.0 --max-spin2z 0.0 --output-name $@ --max-mtotal 100 --max-q 20
###
### NONSPIN FULL REGION
###
#Number of leaf nodes: 72495
#Number of templates (tiles): 72395
#Expected number 45549.9045249
#
#real 6m38.043s
#user 6m36.020s
#sys 0m1.308s
nonspin_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 99.0 --max-mass2 99.0 --min-mass2 1.0 --min-mass1 1.0 --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z -0.0 --max-spin1z 0.0 --min-spin2z -0.0 --max-spin2z 0.0 --output-name $@ --max-mtotal 100 --max-q 20
# > 80 minutes
# > 31,000 templates
# Jon, can you finish this one?
nonspin_sbank.xml.gz:
lalapps_cbc_sbank --verbose --approximant IMRPhenomD --mass1-min 1 --mass1-max 99 --mass2-min 1 --mass2-max 99 --flow $(FLOW) --reference-psd psd_for_treebank.xml.gz --instrument H1 --spin1-min 0 --spin1-max 0 --checkpoint 50 --iterative-match-df-max 2.0 --coarse-match-df 2.0 --aligned-spin --output-filename $@ --match-min $(SBMM) --fhigh-max 1024 --spin2-max 0 --spin2-min 0
###
### Used for doing quick tests
###
test_treebank.xml.gz:
gstlal_inspiral_treebank --verbose --approximant IMRPhenomD --max-mass1 $(M1bf) --max-mass2 $(M2bf) --min-mass2 $(M2bi) --min-mass1 $(M1bi) --min-match $(MM) --flow $(FLOW) --psd-file psd_for_treebank.xml.gz --min-spin1z $(S1bi) --max-spin1z $(S1bf) --min-spin2z $(S2bi) --max-spin2z $(S2bf) --output-name $@
H1-SBANK_SMALL-0-999999999.xml.gz:
# This produces 1269 templates
lalapps_cbc_sbank --verbose --approximant IMRPhenomD --mass1-min $(M1bi) --mass1-max $(M1bf) --mass2-min $(M2bi) --mass2-max $(M2bf) --flow $(FLOW) --reference-psd psd_for_treebank.xml.gz --instrument H1 --spin1-min $(S1bi) --spin1-max $(S1bf) --checkpoint 50 --iterative-match-df-max 2.0 --coarse-match-df 2.0 --aligned-spin --output-filename $@ --match-min $(SBMM) --fhigh-max 1024 --spin2-max $(S2bf) --spin2-min $(S2bi)
H1-TMPLTBANK_H1_BNS_NSBH_BBH-941365351-2048.xml.gz:
# I've written here PhenomB, but I don't think lalapps_tmpltbank actually uses the approximant for anything.
# I'm using an aLIGO HPZD noise model here
# Is the PN order correct?
# This produces 57465 templates!
lalapps_tmpltbank --grid-spacing Hexagonal --space Tau0Tau3 \
--order twoPN --minimal-match 0.97 --spectrum-type aLIGOZDHiP \
--minimum-mass 1.0 --maximum-mass 99.0 --max-total-mass 100.0 \
--approximant IMRPhenomB --low-frequency-cutoff $(FLOW) \
--dynamic-range-exponent 69.0 --disable-high-pass --user-tag BNS_NSBH_BBH \
--gps-end-time 941367399 --gps-start-time 941365351 \
--min-high-freq-cutoff LRD --segment-length 65536 --num-freq-cutoffs 1 \
--sample-rate 16384 --high-frequency-cutoff 8192.0 \
--max-high-freq-cutoff LRD --write-compress --ifo-tag H1 --verbose
H1-TMPLTBANK_H1_SMALL-941365351-2048.xml.gz:
# This produces 3248 templates!
lalapps_tmpltbank --grid-spacing Hexagonal --space Tau0Tau3 \
--order twoPN --minimal-match 0.97 --spectrum-type aLIGONoSRMLoP \
--minimum-mass $(M1bi) --maximum-mass $(M1bf) \
--approximant IMRPhenomB --low-frequency-cutoff $(FLOW) \
--dynamic-range-exponent 69.0 --disable-high-pass --user-tag SMALL \
--gps-end-time 941367399 --gps-start-time 941365351 \
--min-high-freq-cutoff LRD --segment-length 65536 --num-freq-cutoffs 1 \
--sample-rate 16384 --high-frequency-cutoff 8192.0 \
--max-high-freq-cutoff LRD --write-compress --ifo-tag H1 --verbose
# For some crazy reason, if I add --max-total-mass 100, I get 2x more templates!?
HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml:
lalapps_inspinj --aligned --max-mass1 $(M1f) --max-mass2 $(M2f) --min-mass2 $(M2i) --min-mass1 $(M1i) --m-distr componentMass --max-spin1 $(S1f) --max-spin2 0.0 --time-step 10 --taper-injection startend --gps-end-time 1000010000 --min-distance 10000 --max-distance 1000000 --waveform IMRPhenomD --polarization uniform --gps-start-time 1000000000 --l-distr random --d-distr volume --i-distr uniform --max-inc 179.99 --min-spin2 0.0 --min-spin1 $(S1i) --enable-spin --f-lower $(FLOW) --seed 1234 --output HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml
treebank.h5: test_treebank.xml.gz HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml
lalapps_cbc_sbank_sim --template-bank $< --injection-file HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml --injection-approx IMRPhenomD --verbose --flow 30.0 --template-approx IMRPhenomD --cache-waveforms --reference-psd psd.xml.gz --instrument=H1 --user-tag treebank
lalapps_cbc_sbank_sim --template-bank $< --injection-file HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml --injection-approx IMRPhenomD --verbose --flow $(FLOW) --template-approx IMRPhenomD --cache-waveforms --reference-psd psd.xml.gz --instrument=H1 --user-tag treebank
sbank.h5: H1-SBANK_SMALL-0-999999999.xml.gz HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml
lalapps_cbc_sbank_sim --template-bank $< --injection-file HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml --injection-approx IMRPhenomD --verbose --flow 30.0 --template-approx IMRPhenomD --cache-waveforms --reference-psd psd.xml.gz --instrument=H1 --user-tag sbank
tmpltbank.h5: H1-TMPLTBANK_H1_SMALL-941365351-2048.xml.gz HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml
lalapps_cbc_sbank_sim --template-bank $< --injection-file HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml --injection-approx IMRPhenomD --verbose --flow 30.0 --template-approx IMRPhenomD --cache-waveforms --reference-psd psd.xml.gz --instrument=H1 --user-tag tmpltbank
lalapps_cbc_sbank_sim --template-bank $< --injection-file HL-INJECTIONS_1234_IMRPD-1000000000-10000.xml --injection-approx IMRPhenomD --verbose --flow $(FLOW) --template-approx IMRPhenomD --cache-waveforms --reference-psd psd.xml.gz --instrument=H1 --user-tag sbank
banksim_match_vs_injm1.png: treebank.h5
lalapps_cbc_sbank_plot_sim $<
clean :
rm BBH*_treebank.xml.gz BNS_treebank.xml.gz NSBH*_treebank.xml.gz spin_treebank.xml.gz
rm *BBH*_treebank.xml.gz *BNS_treebank.xml.gz *NSBH*_treebank.xml.gz *spin_treebank.xml.gz
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment