RdChannelBandpassFilter.cc

Go to the documentation of this file.

#include "RdChannelBandpassFilter.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/Channel.h>
#include <revt/Station.h>

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


namespace RdChannelBandpassFilter {

  VModule::ResultFlag
  RdChannelBandpassFilter::Init()
  {
    Branch topBranch =
      CentralConfig::GetInstance()->GetTopBranch("RdChannelBandpassFilter");
    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 Channel 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
  RdChannelBandpassFilter::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 and for every station through every channel
    for (auto& station : rEvent.StationsRange()) {
      for (auto& channel : station.ChannelsRange()) {

        if (!channel.IsActive())
          continue;

        // Work on the frequency spectrum of this channel and apply the bandpass filter
        ChannelFrequencySpectrum& spectrum = channel.GetChannelFrequencySpectrum();

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


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

}