RdStationBandpassFilter.cc

Go to the documentation of this file.

#include "RdStationBandpassFilter.h"

#include <fwk/CentralConfig.h>

#include <utl/config.h>
#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/RectangleFilter.h>

#include <evt/Event.h>
#include <revt/REvent.h>
#include <revt/Station.h>


using namespace revt;
using namespace fwk;
using namespace utl;
using namespace std;


namespace RdStationBandpassFilter {

  VModule::ResultFlag
  RdStationBandpassFilter::Init()
  {
    Branch topBranch =
      CentralConfig::GetInstance()->GetTopBranch("RdStationBandpassFilter");
    topBranch.GetChild("LowerCutoffFrequency").GetData(fLowerCutoffFrequency);
    topBranch.GetChild("UpperCutoffFrequency").GetData(fUpperCutoffFrequency);
    topBranch.GetChild("Delay").GetData(fDelay);

    fFilterTypeString = topBranch.GetChild("FilterType").Get<string>();
    if (fFilterTypeString == "Rectangle") {
      fFilterType = Rectangle;
      fFilter = new RectangleFilter(fLowerCutoffFrequency, fUpperCutoffFrequency, fDelay);
    } /* else if (fFilterTypeString == "Gauss") {
      fFilterType = Gauss;
      fFilter = new GaussFilter(fLowerCutoffFrequency, fUpperCutoffFreqeuncy, fDelay);
    } */

    ostringstream info;
    info << "\n Station Bandpass Filter is set up as " << fFilterTypeString
         << " filter spanning " << fLowerCutoffFrequency/megahertz
         << " MHz to " << fUpperCutoffFrequency/megahertz
         << " MHz with delay of " << fDelay/ns << " ns.";
    INFO(info);

    return eSuccess;
  }


  VModule::ResultFlag
  RdStationBandpassFilter::Run(evt::Event& event)
  {

    // Check if there are events at all
    if(!event.HasREvent()) {
      WARNING("No radio event found!");
      return eContinueLoop;
    }

    REvent& rEvent = event.GetREvent();

    // loop through stations
    for (auto& station : rEvent.StationsRange()) {

      // Work on the frequency spectrum of this station and apply the bandpass filter
      StationFrequencySpectrum& spectrum = station.GetStationFrequencySpectrum();

      complex<double> currentTransmissivity;
      for (StationFrequencySpectrum::SizeType i = 0; i < spectrum.GetSize(); ++i) {
        currentTransmissivity = fFilter->TransmissivityAtFrequency(station.GetFrequencyOfBin(i));
        spectrum[i] *= currentTransmissivity;
      }
    }
    return eSuccess;
  }


  VModule::ResultFlag
  RdStationBandpassFilter::Finish()
  {
    delete fFilter;
    return eSuccess;
  }

}