RdChannelTimeSeriesTaperer.cc

Go to the documentation of this file.

#include "RdChannelTimeSeriesTaperer.h"

#include <fwk/CentralConfig.h>

#include <utl/config.h>
#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/HannWindow.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 RdChannelTimeSeriesTaperer {

  VModule::ResultFlag
  RdChannelTimeSeriesTaperer::Init()
  {
    Branch topBranch =
      CentralConfig::GetInstance()->GetTopBranch("RdChannelTimeSeriesTaperer");
    topBranch.GetChild("RelativeWindowWidthOnEachSide").GetData(fWindowSizeOnEachSide);
    topBranch.GetChild("Renormalize").GetData(fRenormalize);

    string tmpstring = topBranch.GetChild("WindowType").Get<string>();
    if (tmpstring == "Hann") {
      fWindowType = Hann;
      fWindow = new HannWindow(fWindowSizeOnEachSide);
    }

    ostringstream info;
    info << "Channel TimeSeries Windower is set up as " << tmpstring
         << " window with a relative width of " << fWindowSizeOnEachSide
         << " on each side.\n";
    INFO(info);

    return eSuccess;
  }


  VModule::ResultFlag
  RdChannelTimeSeriesTaperer::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 time series this channel
        ChannelTimeSeries& timeSeries = channel.GetChannelTimeSeries();

        // set the trace length of the window for the specific channel
        fWindow->SetTraceLength(timeSeries.GetSize());

        // apply the window
        for (ChannelTimeSeries::SizeType i = 0; i < timeSeries.GetSize(); ++i) {
          timeSeries[i] *= fWindow->GetWeightAtBin(i);
          if (fRenormalize)
            timeSeries[i] *= fWindow->GetRenormalizationFactor();
        }
      }
    }

    return eSuccess;
  }


  VModule::ResultFlag
  RdChannelTimeSeriesTaperer::Finish()
  {
    delete fWindow;

    return eSuccess;
  }

}