Use Resesponse instead of only ArrayFactor in stationresponse simulation

demo/comparison-oskar/stationresponse.cpp

@@ -21,7 +21,7 @@ int main(int argc, char** argv){
     options.element_response_model = response_model;
 
 //     casacore::MeasurementSet ms("/home/vdtol/skalowmini/skalowmini-coef1.MS");
-    casacore::MeasurementSet ms("/home/vdtol//src/EveryBeam/build-release/demo/comparison-oskar/blah.ms");
+    casacore::MeasurementSet ms("blah.ms");
 
     casacore::ScalarMeasColumn<casacore::MDirection> referenceDirColumn(
         ms.field(),
@@ -89,16 +89,19 @@ int main(int argc, char** argv){
             };
             real_t freq0 = 50e6;
 
-//             auto result = station->Response(time, freq, direction, freq0, station0, tile0);
-//             result_arr[i][j][0][0] = result[0][0];
-//             result_arr[i][j][0][1] = result[0][1];
-//             result_arr[i][j][1][0] = result[1][0];
-//             result_arr[i][j][1][1] = result[1][1];
-            auto result = station->ArrayFactor(time, freq, direction, freq0, station0, tile0);
-            result_arr[i][j][0][0] = result[0];
-            result_arr[i][j][0][1] = 0.0;
-            result_arr[i][j][1][0] = 0.0;
-            result_arr[i][j][1][1] = result[1];
+            auto result = station->Response(time, freq, direction, freq0, station0, tile0);
+            result_arr[i][j][0][0] = result[0][0];
+            result_arr[i][j][0][1] = result[0][1];
+            result_arr[i][j][1][0] = result[1][0];
+            result_arr[i][j][1][1] = result[1][1];
+
+//             std::cout << result_arr[i][j][0][0] << ", " << result_arr[i][j][1][0] << std::endl;
+
+//             auto result = station->ArrayFactor(time, freq, direction, freq0, station0, tile0);
+//             result_arr[i][j][0][0] = result[0];
+//             result_arr[i][j][0][1] = 0.0;
+//             result_arr[i][j][1][0] = 0.0;
+//             result_arr[i][j][1][1] = result[1];
 
         }
     }