RdChannelResampler.cc

Go to the documentation of this file.

#include "RdChannelResampler.h"

#include <fwk/CentralConfig.h>
#include <fwk/RunController.h>

#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/config.h>
#include <utl/FFTDataContainerAlgorithm.h>
#include <utl/AugerUnits.h>

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

#include <det/Detector.h>
#include <rdet/RDetector.h>
#include <rdet/Channel.h>


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


namespace RdChannelResampler {

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

    fOverrideChannelSampleRate =
      topBranch.GetChild("OverrideChannelSampleRate").Get<string>() == "yes";

    topBranch.GetChild("OverrideResampleBinning").GetData(fOverrideResampleBinning);

    ostringstream info;
    if (fOverrideChannelSampleRate) {
      info << "\n\tRdChannelResampler uses an OverrideChannelSampleRate with an OverrideResampleBinning of "
           << fOverrideResampleBinning / nanosecond << " ns.\n";
    } else {
      info << "\n\tRdChannelResampler uses sampling rate provided by the detector.\n";
    }
    INFODebug(info);

    return eSuccess;
  }


  VModule::ResultFlag
  RdChannelResampler::Run(evt::Event& event)
  {
    // Check if there is event data at all
    if(!event.HasREvent()) {
      WARNING("No radio event found!");
      return eContinueLoop;
    }

    det::Detector& thedet = det::Detector::GetInstance();
    revt::REvent& rEvent = event.GetREvent();

    ostringstream info;

    // loop through stations and for every station through every channel
    for (auto& station : rEvent.StationsRange()) {
      for (auto& channel : station.ChannelsRange()) {

        const rdet::Channel& channeldet =
          thedet.GetRDetector().GetStation(station.GetId()).GetChannel(channel.GetId());

        double channelBinning = 1. / channeldet.GetSamplingFrequency();
        if (fOverrideChannelSampleRate) {
          channelBinning = fOverrideResampleBinning;
        }

        info << "Resampling channel " << channel.GetId() << " of station "
             << channel.GetStationId() << " from a binning of "
             << channel.GetChannelTimeSeries().GetBinning()<< " to "
             <<  channelBinning << ".";
        INFODebug(info);

        FFTDataContainerAlgorithm::ResampleTimeSeries(channel.GetFFTDataContainer(), channelBinning);
      }
    }

    return eSuccess;
  }


  VModule::ResultFlag
  RdChannelResampler::Finish()
  {
    return eSuccess;
  }

}