RdStationHilbertEnveloper.cc

Go to the documentation of this file.

#include "RdStationHilbertEnveloper.h"

#include <cmath>

#include <fwk/CentralConfig.h>

#include <utl/config.h>
#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/RectangleFilter.h>
#include <utl/FFTDataContainerAlgorithm.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 RdStationHilbertEnveloper {

  RdStationHilbertEnveloper::RdStationHilbertEnveloper()
  {
  }

  RdStationHilbertEnveloper::~RdStationHilbertEnveloper()
  {
  }

  VModule::ResultFlag
  RdStationHilbertEnveloper::Init()
  {

    INFO("RdStationHilbertEnveloper::Init()");


    // Read in the configurations of the xml file
    Branch topBranch =
        CentralConfig::GetInstance()->GetTopBranch("RdStationHilbertEnveloper");
    topBranch.GetChild("InfoLevel").GetData(fInfoLevel);

    return eSuccess;
  }


  VModule::ResultFlag
  RdStationHilbertEnveloper::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();
    vector<double> testtrace;

    // loop through stations and for every station through every channel
    for (REvent::StationIterator sIt = rEvent.StationsBegin();
         sIt != rEvent.StationsEnd(); ++sIt) {
      Station& station = *sIt;

      ostringstream info;
      if (fInfoLevel >= eIntermediate)
        info <<  "Apply Hilbert envelope on station " << station.GetId();

      FFTDataContainerAlgorithm::HilbertEnvelope(station.GetFFTDataContainer());

    }
    return eSuccess;

  }

  VModule::ResultFlag
  RdStationHilbertEnveloper::Finish()
  {

    INFO("RdStationHilbertEnveloper::Finish()");

    return eSuccess;
  }


 }