RdCheckFootprint.cc

Go to the documentation of this file.

#include "RdCheckFootprint.h"

using namespace std;
using namespace revt;
using namespace evt;
using namespace fwk;
using namespace utl;
//using namespace sevt;

#define OUT(x) if ((x) <= fInfoLevel) cerr << "  "

namespace RdCheckFootprint {

  RdCheckFootprint::RdCheckFootprint() :
    fInfoLevel(eNone),
    fCherenkovRadius(150.),
    fCheckOutliers(0),
    fSigmaMultiplier(1.),
    fTolerance(1),
    fLDFUncertainty(0.15)
  {
  }

  RdCheckFootprint::~RdCheckFootprint()
  {
  }

  VModule::ResultFlag RdCheckFootprint::Init()
  {

    INFO("RdCheckFootprint::Init()");

    Branch topBranch = CentralConfig::GetInstance()->GetTopBranch("RdCheckFootprint");
    topBranch.GetChild("InfoLevel").GetData(fInfoLevel);
    topBranch.GetChild("CherenkovRadius").GetData(fCherenkovRadius);
    topBranch.GetChild("Tolerance").GetData(fTolerance);
    topBranch.GetChild("CheckOutliers").GetData(fCheckOutliers);
    topBranch.GetChild("SigmaMultiplier").GetData(fSigmaMultiplier);
    topBranch.GetChild("LDFUncertainty").GetData(fLDFUncertainty);
    fLDF.clear();

    return eSuccess;
  }

  VModule::ResultFlag RdCheckFootprint::Run(evt::Event& event)
  {
    REvent& rEvent = event.GetREvent();
    evt::ShowerRRecData& showerrrec = event.GetRecShower().GetRRecShower();
    const rdet::RDetector& rDetector = det::Detector::GetInstance().GetRDetector();

    utl::Vector showerAxis = showerrrec.GetAxis();
    utl::Point coordinateOrigin = showerrrec.GetCoordinateOrigin();

    fLDF.clear();

    for(REvent::StationIterator sIt = rEvent.StationsBegin(); sIt != rEvent.StationsEnd(); ++sIt) {
      Station& station = *sIt;
      const utl::Point sPos = rDetector.GetStation(station).GetPosition();
      const utl::Point core = showerrrec.GetCorePosition();
      double distanceToShowerAxis = utl::RadioGeometryUtilities::GetDistanceToAxis(showerAxis, core, sPos);
      fLDF.push_back(ldf(station.GetId(),distanceToShowerAxis));
    }
    std::sort(fLDF.begin(),fLDF.end(),sort_ldf);

    int empty_ant = 0;

    double ldfUpperLimit = 0.0;

    double currentSignal = 0.0;

    for(unsigned int i = 0; i < fLDF.size(); ++i) {
      if(rEvent.GetStation(fLDF[i].id).GetRejectionStatus() & (revt::eCheckAntennaBackgroundLevel | revt::ePolarizationDeselected)) {
        continue;
      }
      if(rEvent.GetStation(fLDF[i].id).GetRecData().GetPulseFound()) {
        if(fLDF[i].distance < fCherenkovRadius) {
          continue;
        }
        if(ldfUpperLimit == 0.0) {
          ldfUpperLimit = (1. + fLDFUncertainty) * (rEvent.GetStation(fLDF[i].id).GetRecData().GetParameter(revt::eSignal)
                          + fSigmaMultiplier * rEvent.GetStation(fLDF[i].id).GetRecData().GetParameterError(revt::eSignal));
        }
        currentSignal = rEvent.GetStation(fLDF[i].id).GetRecData().GetParameter(revt::eSignal);
        if(((currentSignal > ldfUpperLimit) && fCheckOutliers) || (empty_ant > fTolerance)) {
          INFO("Applying footprint cut");
          if(currentSignal > ldfUpperLimit)
            rEvent.GetStation(fLDF[i].id).SetRejected(revt::eCheckFootprintDeselected);
          if(empty_ant > fTolerance)
            rEvent.GetStation(fLDF[i].id).SetRejected(revt::eLonelyStation);
          rEvent.GetStation(fLDF[i].id).GetRecData().SetPulseFound(false);
          showerrrec.SetParameter(revt::eNumberOfStationsWithPulseFound, showerrrec.GetParameter(revt::eNumberOfStationsWithPulseFound)-1);
        } else {
          empty_ant = 0;
          if( (1. + fLDFUncertainty) * rEvent.GetStation(fLDF[i].id).GetRecData().GetParameter(revt::eSignal)
              + fSigmaMultiplier * rEvent.GetStation(fLDF[i].id).GetRecData().GetParameterError(revt::eSignal) < ldfUpperLimit) {
            ldfUpperLimit = (1. + fLDFUncertainty) * (rEvent.GetStation(fLDF[i].id).GetRecData().GetParameter(revt::eSignal)
                            + fSigmaMultiplier * rEvent.GetStation(fLDF[i].id).GetRecData().GetParameterError(revt::eSignal));
          }
        }
      } else {
        ++empty_ant;
      }
    }
    return eSuccess;
  }

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


}