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

share: replace online makefiles with proposed versions for ER7

parent 8f5d846c
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
gstlal-inspiral/share/Makefile.BNSERBank deleted 100644 → 0
+ 0
− 88
View file @ 8f5d846c
FSTART = 871147316
FSTOP = 871149864
START = 871147516
STOP = 871149564
# Template bank parameters
MIN_MASS = 0.98
MAX_MASS = 2.05
MIN_TOTAL_MASS = 1.96
MAX_TOTAL_MASS = 4.1
LOW_FREQUENCY_CUTOFF = 40.0
HIGH_PASS_FREQ = 30
HIGH_FREQUENCY_CUTOFF = 2047.0
SAMPLE_RATE = 4096
NUM_SPLIT_TEMPLATES = 200
OVERLAP = 40
MM = 0.975
NUMBANKS = 10
all : bank.dag
FAKE :
gstlal_fake_frames --gps-start-time $(FSTART) --gps-end-time $(FSTOP) --channel-name H1=FAKE-STRAIN --verbose --data-source AdvLIGO --frame-type FAKE --output-path FAKE
frame.cache : FAKE
ls FAKE/H-FAKE-871/*.gwf | lalapps_path2cache > frame.cache
H1-TMPLTBANK-871147516-2048.xml : frame.cache
lalapps_tmpltbank \
--disable-compute-moments \
--grid-spacing Hexagonal \
--dynamic-range-exponent 69.0 \
--enable-high-pass $(HIGH_PASS_FREQ) \
--high-pass-order 8 \
--strain-high-pass-order 8 \
--minimum-mass $(MIN_MASS) \
--maximum-mass $(MAX_MASS) \
--min-total-mass $(MIN_TOTAL_MASS) \
--max-total-mass $(MAX_TOTAL_MASS) \
--gps-start-time $(START) \
--gps-end-time $(STOP) \
--calibrated-data real_8 \
--channel-name H1:FAKE-STRAIN \
--space Tau0Tau3 \
--number-of-segments 15 \
--minimal-match $(MM) \
--candle-snr 8 \
--high-pass-attenuation 0.1 \
--min-high-freq-cutoff SchwarzISCO \
--segment-length 1048576 \
--low-frequency-cutoff $(LOW_FREQUENCY_CUTOFF) \
--pad-data 8 \
--num-freq-cutoffs 1 \
--sample-rate $(SAMPLE_RATE) \
--high-frequency-cutoff $(HIGH_FREQUENCY_CUTOFF) \
--resample-filter ldas \
--strain-high-pass-atten 0.1 \
--strain-high-pass-freq $(HIGH_PASS_FREQ) \
--frame-cache frame.cache \
--max-high-freq-cutoff SchwarzISCO \
--approximant TaylorF2 \
--order twoPN \
--spectrum-type median \
--verbose
%_split_bank.cache: H1-TMPLTBANK-871147516-2048.xml
mkdir -p $*_split_bank
gstlal_bank_splitter --overlap $(OVERLAP) --instrument $* --n $(NUM_SPLIT_TEMPLATES) --sort-by mchirp --add-f-final --max-f-final $(HIGH_FREQUENCY_CUTOFF) H1-TMPLTBANK-871147516-2048.xml --output-cache $@ --output-path $*_split_bank
%psd.xml.gz :
gstlal_psd_xml_from_asd_txt --instrument $* --output $@ /usr/share/gstlal/ZERO_DET_high_P.txt
V1psd.xml.gz :
cp /usr/share/gstlal/virgo_online_psd.xml.gz V1psd.xml.gz
reference_psd.xml: H1psd.xml.gz L1psd.xml.gz V1psd.xml.gz
ligolw_add --output reference_psd.xml H1psd.xml.gz L1psd.xml.gz V1psd.xml.gz
%_bank.dag : %_split_bank.cache reference_psd.xml
gstlal_inspiral_svd_bank_pipe --autocorrelation-length 351 --instrument $* --reference-psd reference_psd.xml --bank-cache $< --overlap $(OVERLAP) --flow $(LOW_FREQUENCY_CUTOFF) --output-name $@ --num-banks $(NUMBANKS)
bank.dag : H1_bank.dag L1_bank.dag V1_bank.dag
cat H1_bank.dag L1_bank.dag V1_bank.dag > bank.dag
rm -f H1_bank.dag L1_bank.dag V1_bank.dag
clean :
rm -rf *.sub* *.dag* *.cache *.sh *.xml *.gz logs gstlal_svd_bank* *_split_bank
gstlal-inspiral/share/Makefile.ERHMBank deleted 100644 → 0
+ 0
− 98
View file @ 8f5d846c
FSTART = 871147316
FSTOP = 871149864
START = 871147516
STOP = 871149564
# Template bank parameters
MIN_MASS = 5.0
MAX_MASS = 95.0
MIN_TOTAL_MASS = 10.0
MAX_TOTAL_MASS = 100.0
MIN_ETA = 0.12
MAX_ETA = 0.25
LOW_FREQUENCY_CUTOFF = 15.0
HIGH_PASS_FREQ = 10
HIGH_FREQUENCY_CUTOFF = 2047.0
SAMPLE_RATE = 4096
NUM_SPLIT_TEMPLATES = 200
OVERLAP = 50
MM = 0.97
all : dag
frame.cache :
gstlal_fake_aligo_frames --gps-start-time $(FSTART) --gps-end-time $(FSTOP) --channel-name H1=FAKE-STRAIN --verbose --data-source AdvLIGO --frame-type FAKE
ls H1/H-H1_FAKE-8711/*.gwf | lalapps_path2cache > frame.cache
H1-TMPLTBANK-871147516-2048.xml : frame.cache
lalapps_tmpltbank \
--disable-compute-moments \
--grid-spacing Hexagonal \
--dynamic-range-exponent 69.0 \
--enable-high-pass $(HIGH_PASS_FREQ) \
--high-pass-order 8 \
--strain-high-pass-order 8 \
--minimum-mass $(MIN_MASS) \
--maximum-mass $(MAX_MASS) \
--min-total-mass $(MIN_TOTAL_MASS) \
--max-total-mass $(MAX_TOTAL_MASS) \
--min-eta $(MIN_ETA) \
--max-eta $(MAX_ETA) \
--gps-start-time $(START) \
--gps-end-time $(STOP) \
--standard-candle \
--calibrated-data real_8 \
--candle-mass1 1 \
--candle-mass2 1 \
--channel-name H1:FAKE-STRAIN \
--space Tau0Tau3 \
--number-of-segments 15 \
--minimal-match $(MM) \
--candle-snr 8 \
--high-pass-attenuation 0.1 \
--min-high-freq-cutoff ERD \
--segment-length 1048576 \
--low-frequency-cutoff $(LOW_FREQUENCY_CUTOFF) \
--pad-data 8 \
--num-freq-cutoffs 1 \
--sample-rate $(SAMPLE_RATE) \
--high-frequency-cutoff $(HIGH_FREQUENCY_CUTOFF) \
--resample-filter ldas \
--strain-high-pass-atten 0.1 \
--strain-high-pass-freq $(HIGH_PASS_FREQ) \
--frame-cache frame.cache \
--max-high-freq-cutoff ERD \
--approximant TaylorF2 \
--order twoPN \
--spectrum-type median \
--verbose \
&& gstlal_add_spins_to_bank --input H1-TMPLTBANK-871147516-2048.xml --output H1-IMR_TMPLTBANK-871147516-2048.xml --approximant IMRPhenomB --chi -0.6 --chi -0.3 --chi 0. --chi 0.6 --chi 0.3 \
&& gstlal_bank_splitter --overlap $(OVERLAP) -i H1 -n $(NUM_SPLIT_TEMPLATES) -s mchirp -F -M $(HIGH_FREQUENCY_CUTOFF) H1-IMR_TMPLTBANK-871147516-2048.xml \
&& gstlal_bank_splitter --overlap $(OVERLAP) -i V1 -n $(NUM_SPLIT_TEMPLATES) -s mchirp -F -M $(HIGH_FREQUENCY_CUTOFF) H1-IMR_TMPLTBANK-871147516-2048.xml \
&& gstlal_bank_splitter --overlap $(OVERLAP) -i L1 -n $(NUM_SPLIT_TEMPLATES) -s mchirp -F -M $(HIGH_FREQUENCY_CUTOFF) H1-IMR_TMPLTBANK-871147516-2048.xml
H1psd.xml.gz :
gstlal_psd_xml_from_asd_txt --instrument H1 --output H1psd.xml.gz /usr/share/gstlal/ZERO_DET_high_P.txt
L1psd.xml.gz :
gstlal_psd_xml_from_asd_txt --instrument L1 --output L1psd.xml.gz /usr/share/gstlal/ZERO_DET_high_P.txt
V1psd.xml.gz :
cp /usr/share/gstlal/virgo_online_psd.xml.gz V1psd.xml.gz
reference_psd.xml: H1psd.xml.gz L1psd.xml.gz V1psd.xml.gz
ligolw_add --output reference_psd.xml H1psd.xml.gz L1psd.xml.gz V1psd.xml.gz
dag : H1-TMPLTBANK-871147516-2048.xml reference_psd.xml
ls *H1_split_bank-H1-IMR_TMPLTBANK-871147516-2048.xml | lalapps_path2cache > H1_split_bank.cache
ls *V1_split_bank-H1-IMR_TMPLTBANK-871147516-2048.xml | lalapps_path2cache > V1_split_bank.cache
ls *L1_split_bank-H1-IMR_TMPLTBANK-871147516-2048.xml | lalapps_path2cache > L1_split_bank.cache
gstlal_inspiral_svd_bank_pipe --autocorrelation-length 351 --instrument H1 --reference-psd reference_psd.xml --bank-cache H1_split_bank.cache --overlap $(OVERLAP) --flow $(LOW_FREQUENCY_CUTOFF) --output-name H1_bank
gstlal_inspiral_svd_bank_pipe --autocorrelation-length 351 --instrument V1 --reference-psd reference_psd.xml --bank-cache V1_split_bank.cache --overlap $(OVERLAP) --flow $(LOW_FREQUENCY_CUTOFF) --output-name V1_bank
gstlal_inspiral_svd_bank_pipe --autocorrelation-length 351 --instrument L1 --reference-psd reference_psd.xml --bank-cache L1_split_bank.cache --overlap $(OVERLAP) --flow $(LOW_FREQUENCY_CUTOFF) --output-name L1_bank
cat H1_bank.dag L1_bank.dag V1_bank.dag > bank.dag
rm -f H1_bank.dag L1_bank.dag V1_bank.dag
realclean :
rm -rf *.sub *.dag* *.cache *.sh H1-TMPLTBANK-871147516-2048.xml *split_bank-H1-IMR_TMPLTBANK-871147516-2048.xml H1-IMR_TMPLTBANK-871147516-2048.xml H1psd.xml.gz L1psd.xml.gz V1psd.xml.gz segments.H.xml logs
gstlal-inspiral/share/Makefile.ERTrigs deleted 100644 → 0
+ 0
− 37
View file @ 8f5d846c
H1_BANK_CACHE = ../bns_bank/H1_bank.cache
L1_BANK_CACHE = ../bns_bank/L1_bank.cache
V1_BANK_CACHE = ../bns_bank/V1_bank.cache
NUMBANKS = 4,5,6
MAXJOBS = 10
FINAL_FAR_THRESH = 1e-4
HT_GATE_THRESH = 15
MIDDLE_FAR_THRESH = 1e-4
INITIAL_FAR_THRESH = -1 # disable uploads
FIRSTRIDE = 2
THINCAINT = 4
PEAKTIME = 2
all : dag
seed : marginalized_likelihood.xml.gz likelihood.xml.gz
gstlal_ll_trigger_pipe --bank-cache H1=$(H1_BANK_CACHE),L1=$(L1_BANK_CACHE),V1=$(V1_BANK_CACHE) --num-banks $(NUMBANKS) --max-jobs $(MAXJOBS) --likelihood-file likelihood.xml.gz --channel=H1=FAKE-STRAIN --channel=V1=FAKE_h_16384Hz_4R --channel=L1=FAKE-STRAIN --gracedb-far-threshold $(INITIAL_FAR_THRESH) --control-peak-time $(PEAKTIME) --fir-stride $(FIRSTRIDE) --marginalized-likelihood-file marginalized_likelihood.xml.gz --gracedb-group CBC --gracedb-type LowMass --thinca-interval $(THINCAINT) --ht-gate-threshold $(HT_GATE_THRESH) --copy-likelihoods --state-vector-on-bits=H1=0x7 --state-vector-on-bits=L1=0x7 --state-vector-on-bits=V1=0x7 --state-vector-off-bits=H1=0x160 --state-vector-off-bits=L1=0x160 --state-vector-off-bits=V1=0x100 --data-source framexmit
dag : marginalized_likelihood.xml.gz likelihood.xml.gz
gstlal_ll_trigger_pipe --bank-cache H1=$(H1_BANK_CACHE),L1=$(L1_BANK_CACHE),V1=$(V1_BANK_CACHE) --num-banks $(NUMBANKS) --max-jobs $(MAXJOBS) --likelihood-file likelihood.xml.gz --channel=H1=FAKE-STRAIN --channel=V1=FAKE_h_16384Hz_4R --channel=L1=FAKE-STRAIN --gracedb-far-threshold $(MIDDLE_FAR_THRESH) --control-peak-time $(PEAKTIME) --fir-stride $(FIRSTRIDE) --marginalized-likelihood-file marginalized_likelihood.xml.gz --gracedb-group CBC --gracedb-type LowMass --thinca-interval $(THINCAINT) --ht-gate-threshold $(HT_GATE_THRESH) --state-vector-on-bits=H1=0x7 --state-vector-on-bits=L1=0x7 --state-vector-on-bits=V1=0x7 --state-vector-off-bits=H1=0x160 --state-vector-off-bits=L1=0x160 --state-vector-off-bits=V1=0x100 --data-source framexmit
# Note that this does not change the marginalized_likelihood
reset-likelihood :
gstlal_inspiral_reset_likelihood --marginalized-likelihood-file marginalized_likelihood.xml.gz --verbose 0*likelihood.xml.gz
set-far-thresh :
gstlal_ll_inspiral_gracedb_threshold --gracedb-far-threshold $(FINAL_FAR_THRESH) *registry.txt
likelihood.xml.gz :
gstlal_inspiral_create_prior_diststats --trials-far-thresh 1e-3 --instrument H1 --instrument L1 --instrument V1 --write-likelihood $@ --verbose
marginalized_likelihood.xml.gz : likelihood.xml.gz
gstlal_inspiral_marginalize_likelihood --output $@ --verbose likelihood.xml.gz
realclean :
rm -r *.sub *.dag* *.cache *.sh logs *.xml *.gz *.sqlite *.txt
gstlal-inspiral/share/Makefile.online_analysis 0 → 100644
+ 63
− 0
View file @ da918757
H1_BANK_CACHE = ../bank_red_region_spin_early_aligo_2048Hz/H1_bank.cache
L1_BANK_CACHE = ../bank_red_region_spin_early_aligo_2048Hz/L1_bank.cache
V1_BANK_CACHE = ../bank_red_region_spin_early_aligo_2048Hz/V1_bank.cache
dag : marginalized_likelihood.xml.gz prior.cache
gstlal_ll_inspiral_pipe \
--bank-cache H1=$(H1_BANK_CACHE),L1=$(L1_BANK_CACHE),V1=$(V1_BANK_CACHE) \
--max-jobs 40 \
--likelihood-cache prior.cache \
--channel=H1=GDS-CALIB_STRAIN \
--channel=V1=FAKE_h_16384Hz_4R \
--channel=L1=GDS-CALIB_STRAIN \
--dq-channel-name=L1=GDS-CALIB_STATE_VECTOR \
--dq-channel-name=H1=GDS-CALIB_STATE_VECTOR \
--dq-channel-name=V1=LLD-DQ_VECTOR --psd-fft-length 4 \
--state-vector-on-bits=H1=0x1 \
--state-vector-on-bits=L1=0x1 \
--state-vector-on-bits=V1=0x7 \
--state-vector-off-bits=H1=0x0 \
--state-vector-off-bits=L1=0x0 \
--state-vector-off-bits=V1=0x0 \
--gracedb-far-threshold -1 \
--control-peak-time 0 \
--fir-stride 1 \
--marginalized-likelihood-file marginalized_likelihood.xml.gz \
--gracedb-group CBC \
--gracedb-search LowMass \
--thinca-interval 1 \
--ht-gate-threshold 15 \
--data-source framexmit \
--likelihood-snapshot-interval 14400 \
--lvalert-listener-program gstlal_inspiral_followups_from_gracedb \
--lvalert-listener-program gstlal_inspiral_lvalert_psd_plotter
set-far-thresh :
gstlal_ll_inspiral_gracedb_threshold \
--gracedb-far-threshold $(FINAL_FAR_THRESH) \
*registry.txt
prior.cache :
gstlal_ll_inspiral_create_prior_diststats \
--write-likelihood-cache $@ \
--segment-and-horizon=H1:1000000000:1000010000:60 \
--segment-and-horizon=L1:1000000000:1000010000:120 \
--segment-and-horizon=V1:1000000000:1000010000:30 \
--segment-and-horizon=H1:1000010000:1000020000:60 \
--segment-and-horizon=L1:1000010000:1000020000:120 \
--segment-and-horizon=H1:1000020000:1000030000:60 \
--segment-and-horizon=V1:1000020000:1000030000:30 \
--segment-and-horizon=L1:1000030000:1000040000:120 \
--segment-and-horizon=V1:1000030000:1000040000:30 \
--num-banks $(shell wc -l $(H1_BANK_CACHE) | awk '{print $1}') \
--verbose
marginalized_likelihood.xml.gz : prior.cache
gstlal_inspiral_marginalize_likelihood \
--output $@ \
--verbose \
--likelihood-cache $<
clean :
rm -rf gstlal_inspiral gstlal_inspiral_marginalize_likelihoods_online gstlal_ll_inspiral_get_urls lvalert_listen
rm -rf *.txt lvalert.ini *.gz trigger_pipe.* *.sub logs lvalert.sh node* *.xml
gstlal-inspiral/share/Makefile.online_bank 0 → 100644
+ 42
− 0
View file @ da918757
FSTART = 871147316
FSTOP = 871149864
START = 871147516
STOP = 871149564
# Template bank parameters
LOW_FREQUENCY_CUTOFF = 30.0
HIGH_PASS_FREQ = 30
HIGH_FREQUENCY_CUTOFF = 1024
SAMPLE_RATE = 2048
NUM_SPLIT_TEMPLATES = 130
OVERLAP = 20
APPROXIMANT = TaylorF2
BANK_PROGRAM = pycbc_geom_aligned_bank
#FIXME
NUMBANKS = 21,22,23,24,25,26,27,28,29,30,31,32
all : bank.dag
LL_5mc_BHspin_red_30early.xml.gz :
cp $(shell echo $(dirname $(dirname $(which gstlal_inspiral))))/share/gstlal/LL_5mc_BHspin_red_30early.xml.gz .
%_split_bank.cache: LL_5mc_BHspin_red_30early.xml.gz
mkdir -p $*_split_bank
gstlal_bank_splitter --group-by-chi --overlap $(OVERLAP) --instrument $* --n $(NUM_SPLIT_TEMPLATES) --sort-by mchirp --add-f-final --max-f-final $(HIGH_FREQUENCY_CUTOFF) $< --output-cache $@ --output-path $*_split_bank --approximant $(APPROXIMANT) --bank-program $(BANK_PROGRAM)
%psd.xml.gz :
gstlal_psd_xml_from_asd_txt --instrument $* --output $@ /usr/share/gstlal/early_aligo_asd.txt
reference_psd.xml: H1psd.xml.gz L1psd.xml.gz V1psd.xml.gz
ligolw_add --output reference_psd.xml H1psd.xml.gz L1psd.xml.gz V1psd.xml.gz
%_bank.dag : %_split_bank.cache reference_psd.xml
cp $< tmp
gstlal_inspiral_svd_bank_pipe --autocorrelation-length 351 --instrument $* --reference-psd reference_psd.xml --bank-cache $< --overlap $(OVERLAP) --flow $(LOW_FREQUENCY_CUTOFF) --output-name $@ --num-banks $(NUMBANKS)
bank.dag : H1_bank.dag L1_bank.dag V1_bank.dag
cat H1_bank.dag L1_bank.dag V1_bank.dag > bank.dag
rm -f H1_bank.dag L1_bank.dag V1_bank.dag
clean :
rm -rf *.sub* *.dag* *.cache *.sh *.xml *.gz logs gstlal_svd_bank* *split_bank
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