from spectractor.simulation.atmosphere import plot_transmission_simple
from spectractor import parameters
fig = plt.figure()
ax = plt.gca()
parameters.LAMBDA_MIN = 500
lambdas, transmissions, errors = load_transmission_file(os.path.join(parameters.THROUGHPUT_DIR, "qecurve.txt"))
plot_transmission_simple(ax, lambdas, transmissions, errors, title="CTIO", label="Quantum efficiency")
lambdas, transmissions, errors = load_transmission_file(os.path.join(parameters.THROUGHPUT_DIR, "lsst_mirrorthroughput.txt"))
plot_transmission_simple(ax, lambdas, transmissions, errors, title="CTIO", label="Mirror 1")
lambdas, transmissions, errors = load_transmission_file(os.path.join(parameters.THROUGHPUT_DIR, "FGB37.txt"))
plot_transmission_simple(ax, lambdas, transmissions, errors, title="CTIO", label="FGB37")
lambdas, transmissions, errors = load_transmission_file(os.path.join(parameters.THROUGHPUT_DIR, "RG715.txt"))
plot_transmission_simple(ax, lambdas, transmissions, errors, title="CTIO", label="RG715")
lambdas, transmissions, errors = load_transmission_file(os.path.join(parameters.THROUGHPUT_DIR, parameters.OBS_FULL_INSTRUMENT_TRANSMISSON))
plot_transmission_simple(ax, lambdas, transmissions, errors, title="CTIO", label="Full instrument")
if parameters.DISPLAY: plt.show()