Skip to content
Snippets Groups Projects
Commit 66dfdc9c authored by Madeline Wade's avatar Madeline Wade
Browse files

Changes to a few plotting scripts for testing

parent 8834ed8e
No related branches found
No related tags found
No related merge requests found
...@@ -6,7 +6,6 @@ from gwpy.timeseries import TimeSeriesDict ...@@ -6,7 +6,6 @@ from gwpy.timeseries import TimeSeriesDict
from gwpy.timeseries import TimeSeries from gwpy.timeseries import TimeSeries
import glob import glob
from math import pi from math import pi
from gwpy.plotter import BodePlot
import numpy import numpy
from optparse import OptionParser, Option from optparse import OptionParser, Option
......
...@@ -4,7 +4,7 @@ ...@@ -4,7 +4,7 @@
########################################################## ##########################################################
import matplotlib as mpl; mpl.use('Agg') import matplotlib as mpl; mpl.use('Agg')
from gwpy.plotter import BodePlot from gwpy.plot import BodePlot
from gwpy.timeseries import TimeSeriesDict from gwpy.timeseries import TimeSeriesDict
from gwpy.timeseries import TimeSeries from gwpy.timeseries import TimeSeries
from gwpy.frequencyseries import FrequencySeries from gwpy.frequencyseries import FrequencySeries
...@@ -154,7 +154,7 @@ if options.analyze_calcs_hoft: ...@@ -154,7 +154,7 @@ if options.analyze_calcs_hoft:
plot.add_frequencyseries(CALCS_tf, dB=False, color='#4ba6ff', linewidth=2) plot.add_frequencyseries(CALCS_tf, dB=False, color='#4ba6ff', linewidth=2)
if options.analyze_additional_hoft: if options.analyze_additional_hoft:
plot.add_frequencyseries(additional_hoft_tf*3995.1, dB = False, color="#94ded7", linewidth=2) plot.add_frequencyseries(additional_hoft_tf*3995.1, dB = False, color="#94ded7", linewidth=2)
plot.add_legend([r'Reference model response function', r'h(t) derived response function'], loc='upper right', fontsize='x-small') #plot.add_legend([r'Reference model response function', r'h(t) derived response function'], loc='upper right', fontsize='x-small')
# FIXME: Figure out how to make the legend and title appropriate and flexible # FIXME: Figure out how to make the legend and title appropriate and flexible
plot.maxes.set_yscale('log') plot.maxes.set_yscale('log')
plot.paxes.set_yscale("linear") plot.paxes.set_yscale("linear")
...@@ -172,7 +172,7 @@ if options.analyze_calcs_hoft: ...@@ -172,7 +172,7 @@ if options.analyze_calcs_hoft:
plot.add_frequencyseries(ratio_CALCS, dB = False, color='#4ba6ff',linewidth=2) plot.add_frequencyseries(ratio_CALCS, dB = False, color='#4ba6ff',linewidth=2)
if options.analyze_additional_hoft: if options.analyze_additional_hoft:
plot.add_frequencyseries(ratio_additional_hoft*3995.1, dB = False, color="#94ded7", linewidth=2) plot.add_frequencyseries(ratio_additional_hoft*3995.1, dB = False, color="#94ded7", linewidth=2)
plot.add_legend([r'h(t) derived response / Reference model response'], loc='upper right', fontsize='small') #plot.add_legend([r'h(t) derived response / Reference model response'], loc='upper right', fontsize='small')
plot.maxes.set_yscale('linear') plot.maxes.set_yscale('linear')
plot.paxes.set_yscale('linear') plot.paxes.set_yscale('linear')
plot.maxes.set_ylim(0,5) plot.maxes.set_ylim(0,5)
...@@ -183,7 +183,7 @@ if options.analyze_calcs_hoft: ...@@ -183,7 +183,7 @@ if options.analyze_calcs_hoft:
plot.add_frequencyseries(ratio_CALCS, dB = False, color='#4ba6ff',linewidth=2) plot.add_frequencyseries(ratio_CALCS, dB = False, color='#4ba6ff',linewidth=2)
if options.analyze_additional_hoft: if options.analyze_additional_hoft:
plot.add_frequencyseries(ratio_additional_hoft*3995.1, dB = False, color="#94ded7", linewidth=2) plot.add_frequencyseries(ratio_additional_hoft*3995.1, dB = False, color="#94ded7", linewidth=2)
plot.add_legend([r'h(t) derived response / Reference model response'], loc='upper right', fontsize='small') #plot.add_legend([r'h(t) derived response / Reference model response'], loc='upper right', fontsize='small')
plot.maxes.set_yscale('linear') plot.maxes.set_yscale('linear')
plot.paxes.set_yscale('linear') plot.paxes.set_yscale('linear')
plot.maxes.set_ylim(.9, 1.1) plot.maxes.set_ylim(.9, 1.1)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment