Modified plotting script for testing and makefile

gstlal-calibration/examples/L1_GDS_CALCS_comparison_plots

@@ -41,19 +41,25 @@ calcs_asd = calcs_asd.filter([30]*6, [0.3]*6, 1e-12 / 4e3)
 
 #plot spectrum
 plot=calcs_asd.plot(label='CALCS h(t) ASD')
-plot.gca().plot(gds_asd,label='GDS h(t) ASD')
+plot.gca().plot(gds_asd,label='DCS h(t) ASD')
 ax = plot.gca()
-ax.set_ylabel = 'Strain [Hz$^{-1/2}$]'
-ax.set_xlabel = 'Frequency [Hz]'
+#ax.set_ylabel = 'Strain [Hz$^{-1/2}$]'
+#ax.set_xlabel = 'Frequency [Hz]'
+plot.set_ylabel('Strain [Hz $^{-1/2}$]', fontname = 'Times', fontsize = 18)
+plot.set_xlabel('Frequency [Hz]', fontname = 'Times', fontsize = 18)
 ax.set_xlim(10,8192)
 ax.set_ylim(1e-24,1e-16)
+ax.legend()
 plot.save('spectrum_comparison.png')
 
 diff = calcs_asd / gds_asd
-plot = diff.plot(label="ASD ratio CALCS / GDS", logy = False)
+plot = diff.plot(label="ASD ratio CALCS / DCS", logy = False)
 ax = plot.gca()
-ax.set_ylabel = 'Strain [Hz$^{-1/2}$]'
-ax.set_xlabel = 'Frequency [Hz]'
+#ax.set_ylabel = 'Strain [Hz$^{-1/2}$]'
+#ax.set_xlabel = 'Frequency [Hz]'
+plot.set_ylabel('Strain [Hz $^{-1/2}$]', fontname = 'Times', fontsize = 18)
+plot.set_xlabel('Frequency [Hz]', fontname = 'Times', fontsize = 18)
 ax.set_xlim(10,5000)
 ax.set_ylim(0.7, 1.3)
-plot.save('CALCS_GDS_residual.png')
+ax.legend()
+plot.save('CALCS_DCS_residual.png')

gstlal-calibration/tests/Makefile.DCS_calibration_test

+# time of GW170817 1187008882
+START = 1187007882 # Thu Aug 17 12:24:24 GMT 2017
+END = 1187008582 # Thu Aug 17 12:36:04 GMT 2017
+FILTERS = 'L1DCS_1175961600.npz'
+# Filters file from aligocalibration/trunk/Runs/O2/GDSFilters/L1DCS_1175961600.npz
+
+all: L1_hoft_frames.cache
+
+L1_raw_frames.cache: 
+	gw_data_find -o L -t L1_R -s $(START) -e $(END) -l --url-type file > $@
+
+L1_hoft_frames.cache: L1_raw_frames.cache
+	gstlal_compute_strain --data-source frames --frame-cache L1_raw_frames.cache --gps-start-time $(START) --gps-end-time $(END) --frame-duration=4 --frames-per-file=1 --filters-file $(FILTERS) --ifo L1 --frame-type L1_TEST --compression-scheme=6 --compression-level=3 --full-calibration --control-sample-rate 16384 --factors-from-filters-file --expected-fcc 376.0 --no-fs --no-srcQ --coherence-uncertainty-threshold 0.02 --kappas-default-to-median --apply-kappatst --apply-kappapu --apply-kappac --verbose
+	ls *.gwf | lalapps_path2cache > $@
+
+clean: 
+	rm *.cache *.gwf