Skip to content
Snippets Groups Projects
Commit 5e93194c authored by Mattia Mancini's avatar Mattia Mancini
Browse files

SSB-47: Refactor the application of the gains

parent b095c8b3
No related branches found
No related tags found
1 merge request!44Merge back holography to master
...@@ -102,7 +102,7 @@ class CalibrationTable: ...@@ -102,7 +102,7 @@ class CalibrationTable:
new_calibration_table.calibration_date = calibration_date new_calibration_table.calibration_date = calibration_date
new_calibration_table.comment = comment new_calibration_table.comment = comment
new_calibration_table.data = gains new_calibration_table.data *= gains
return new_calibration_table return new_calibration_table
def frequencies(self) -> ndarray: def frequencies(self) -> ndarray:
......
from numpy import ndarray from numpy import ndarray
from lofar.calibration.common.datacontainers.holography_dataset import HolographyDataset from lofar.calibration.common.datacontainers.holography_dataset import HolographyDataset
from numpy import array, exp from numpy import zeros_like, exp
from datetime import datetime, date from datetime import datetime, date
from lofar.calibration.common.coordinates import mjd_to_astropy_time from lofar.calibration.common.coordinates import mjd_to_astropy_time
import logging import logging
...@@ -45,7 +45,7 @@ def compute_calibration_table(dataset: HolographyDataset, comment): ...@@ -45,7 +45,7 @@ def compute_calibration_table(dataset: HolographyDataset, comment):
for mode, calibration_table in zip(modes, calibration_tables): for mode, calibration_table in zip(modes, calibration_tables):
frequencies = calibration_table.frequencies() * _MHZ_TO_HZ frequencies = calibration_table.frequencies() * _MHZ_TO_HZ
data = array(calibration_table.data) data = zeros_like(calibration_table.data)
for polarization in interpolated_gains: for polarization in interpolated_gains:
for antenna in interpolated_gains[polarization]: for antenna in interpolated_gains[polarization]:
gains_per_antenna = interpolated_gains[polarization][antenna] gains_per_antenna = interpolated_gains[polarization][antenna]
...@@ -54,7 +54,7 @@ def compute_calibration_table(dataset: HolographyDataset, comment): ...@@ -54,7 +54,7 @@ def compute_calibration_table(dataset: HolographyDataset, comment):
rcu_id = antenna_polarization_mode_to_rcu(int(antenna), rcu_id = antenna_polarization_mode_to_rcu(int(antenna),
polarization, polarization,
int(mode)) int(mode))
data[:, rcu_id] *= derived_gains data[:, rcu_id] = derived_gains
observation_date = mjds_to_datetime(dataset.start_time) observation_date = mjds_to_datetime(dataset.start_time)
observation_date_formatted = observation_date.strftime(_DATETIME_FORMAT) observation_date_formatted = observation_date.strftime(_DATETIME_FORMAT)
......
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