Maaijke Mevius
--- a/scintillation/averaging.py
+++ b/scintillation/averaging.py
-                                                                                 data_array)
-        make_plot(averaged_data_array, time_axis, frequency_axis, station_name, full_path + '.png')
-        create_fits_from_dataset(sample_info, averaged_data_array, full_path + '.fits')
+            'sample_end_frequency': end_frequency}
+        sample_info.update(metadata)
+        
+        
+        frequency_axis = numpy.linspace(start_frequency, end_frequency, data_array.shape[1])
+        filtered_data,flags,flux = apply_bandpass(data_array[start_index:end_index], frequency_axis, 
+                                                  freqaxis=1, timeaxis=0, target = metadata["TARGET"],
+                                                  sample_rate= int(averaging_window_in_samples*3./averaging_window_in_seconds),  
+                                                  #sample every 3 seconds
+                                                  filter_length = 600,    # window size 30 minutes
+                                                  flag_value=5, replace_value="interpolate")
+        S4_data_array = getS4(filtered_data, 
+                              window_size = int(60*averaging_window_in_samples/averaging_window_in_seconds),