RdStationQualityAssessor.cc

Go to the documentation of this file.

#include "RdStationQualityAssessor.h"

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

#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/config.h>
#include <utl/Trace.h>
#include <utl/TraceAlgorithm.h>
#include <utl/TimeStamp.h>
#include <utl/TimeInterval.h>
#include <utl/AugerUnits.h>
#include <utl/AugerException.h>
#include <utl/FFTDataContainerAlgorithm.h>
#include <utl/Math.h>
#include <utl/Vector.h>
#include <utl/AxialVector.h>
#include <utl/AugerCoordinateSystem.h>
#include <utl/CoordinateSystemPtr.h>
#include <utl/UTMPoint.h>

#include <evt/Event.h>
#include <revt/REvent.h>
#include <revt/Header.h>
#include <evt/ShowerRecData.h>
#include <evt/ShowerSimData.h>
#include <evt/ShowerRRecData.h>
#include <evt/ShowerFRecData.h>
#include <evt/ShowerSRecData.h>
#include <revt/Channel.h>

#include <rdet/RDetector.h>
#include <rdet/Station.h>

#include <fwk/LocalCoordinateSystem.h>


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

namespace RdStationQualityAssessor {
    
  RdStationQualityAssessor::RdStationQualityAssessor()
  : fInfoLevel(0)
  {}

  RdStationQualityAssessor::~RdStationQualityAssessor() {}
  
  VModule::ResultFlag RdStationQualityAssessor::Init()
  {
    INFO("RdStationQualityAssessor::Init()");

    // Read in the configurations of the xml file
    Branch topBranch =
      CentralConfig::GetInstance()->GetTopBranch("RdStationQualityAssessor");
    topBranch.GetChild("InfoLevel").GetData(fInfoLevel);
    topBranch.GetChild("VectorialComponent").GetData(fVectorialComponent);
    return eSuccess;
  }


  VModule::ResultFlag RdStationQualityAssessor::Run(evt::Event& event)
  {
    INFO("RdStationQualityAssessor::Run()");
    
    // Check wether there are events at all
    if(!event.HasREvent()) {
      WARNING("No radio event found!");
      return eContinueLoop;
    }

    REvent& rEvent = event.GetREvent();      
    // loop through stations
    for (REvent::SignalStationIterator sIt = rEvent.SignalStationsBegin();
         sIt != rEvent.SignalStationsEnd(); ++sIt) {
      
      StationRecData &recData(sIt->GetRecData());
      double signalSearchWindowStart = recData.GetParameter(eSignalSearchWindowStart);
      double signalSearchWindowStop = recData.GetParameter(eSignalSearchWindowStop);
      
      StationTimeSeries &timeSeries = sIt->GetFFTDataContainer().GetTimeSeries();
      double binning = timeSeries.GetBinning();

      double maxOutsideSignalSearchWindow = 0;
      double amplitude = 0;
      for (unsigned int i = 0; i < timeSeries.GetSize(); ++i) {
        if (fVectorialComponent == 0) // absolute value of vectorial sum
          amplitude = timeSeries[i].GetMag();
        else if (fVectorialComponent == 1) // x
          amplitude = fabs(timeSeries[i][0]);
        else if (fVectorialComponent == 2) // y
          amplitude = fabs(timeSeries[i][1]);
        else if (fVectorialComponent == 3) // z
          amplitude = fabs(timeSeries[i][2]);
        else if (fVectorialComponent == 4) // xy-plane
          amplitude = sqrt(pow(timeSeries[i][0], 2) + pow(timeSeries[i][1], 2)); 
        bool isOutside = i * binning < signalSearchWindowStart or i * binning > signalSearchWindowStop;
        if ( amplitude > maxOutsideSignalSearchWindow and isOutside) {
          maxOutsideSignalSearchWindow = amplitude;
        }
      }
      recData.SetParameter(eMaxAmplitudeOutsideSignalSearchWindow, maxOutsideSignalSearchWindow);
    }
    return eSuccess;
  }

  VModule::ResultFlag
  RdStationQualityAssessor::Finish()
  {

    INFO("RdStationQualityAssessor::Finish()");

    return eSuccess;
  }

}