Skip to content

Commit

Permalink
Small output update
Browse files Browse the repository at this point in the history
  • Loading branch information
LucaMalavolta committed Apr 24, 2024
1 parent 27ef182 commit f90e066
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions SLOPpy/sky_correction.py
Original file line number Diff line number Diff line change
Expand Up @@ -51,7 +51,6 @@ def compute_sky_correction(config_in):
processed[obs] = {}

" computing the ratio between the lamp flux of fiber A and B"
print()
print(" Computing the ratio between the lamp flux of fiber A and B")

processed['ratioAB'] = calib_data_A['lamp'][map_orders_A, :]/calib_data_B['lamp'][map_orders_B, :]
Expand All @@ -61,7 +60,6 @@ def compute_sky_correction(config_in):
wave_difference = \
input_data_A[first_obs]['wave'][map_orders_A, :] - input_data_B[first_obs]['wave'][map_orders_B, :]

print()
print(" Wavelength difference between fiber A and B: ", \
np.average(wave_difference), " +- ", np.std(wave_difference), " \AA")

Expand All @@ -74,6 +72,8 @@ def compute_sky_correction(config_in):
is minimal """

for obs in lists['observations']:
print(obs)

processed[obs]['sky_fibA'] = np.zeros([input_data_A['n_orders'], input_data_A['n_pixels']])
processed[obs]['sky_fibA'][map_orders_A, :] = \
processed['ratioAB'] * input_data_B[obs]['e2ds'][map_orders_B, :]
Expand All @@ -88,7 +88,7 @@ def compute_sky_correction(config_in):
#processed[obs]['e2ds'][processed[obs]['null']] = replacement

save_to_cpickle('skycorrected_fibA', processed, config_in['output'], night)

print()

def plot_sky_correction(config_in, night_input=''):

Expand Down

0 comments on commit f90e066

Please sign in to comment.