RdChannelPedestalRemover.cc

Go to the documentation of this file.

#include "RdChannelPedestalRemover.h"

#include <fwk/CentralConfig.h>

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

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


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


namespace RdChannelPedestalRemover {

  VModule::ResultFlag
  RdChannelPedestalRemover::Init()
  {
    return eSuccess;
  }


  VModule::ResultFlag
  RdChannelPedestalRemover::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 trace of this channel
        revt::ChannelTimeSeries& trace = channel.GetChannelTimeSeries();
       
        if (trace.GetSize() == 0)
          continue; 

        RemoveDCOffset(trace);
      }
    }

    return eSuccess;
  }


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


  void
  RdChannelPedestalRemover::RemoveDCOffset(revt::ChannelTimeSeries& trace) {
    double mean = TraceAlgorithm::Mean(trace, 0, trace.GetSize()-1 );
    for (auto& binValue : trace) {
      binValue -= mean;
    }
  }

}