PixelSelector.cc

Go to the documentation of this file.

#include "PixelSelector.h"
#include <fwk/CentralConfig.h>
#include <evt/Event.h>
#include <fevt/FEvent.h>
#include <fevt/Eye.h>
#include <fevt/EyeRecData.h>
#include <fevt/Pixel.h>
#include <utl/Vector.h>
#include <utl/Reader.h>
#include <det/Detector.h>
#include <fdet/FDetector.h>
#include <fdet/Pixel.h>
#include <utl/MathConstants.h>

#include <sstream>

using namespace utl;


namespace PixelSelectorOG {

  fwk::VModule::ResultFlag
  PixelSelector::Init()
  {
    auto topB = fwk::CentralConfig::GetInstance()->GetTopBranch("PixelSelector");

    if (!topB) {
      ERROR("Configuration for module PixelSelector not found.");
      return eFailure;
    }

    topB.GetChild("minAngle").GetData(fMinAngle);

    // info output
    std::ostringstream info;
    info << "\n"
            " Version: " << GetVersionInfo(VModule::eRevisionNumber) << "\n"
            " Parameters:\n"
            "   minimum angle: " << fMinAngle/deg << "deg";
    INFO(info);

    return eSuccess;
  }


  fwk::VModule::ResultFlag
  PixelSelector::CheckIsolated(evt::Event& event)
  {
    if (!event.HasFEvent())
      return eSuccess;

    auto& fevent = event.GetFEvent();

    //unsigned int neyes = 0;  // number of eyes with enough pixels

    for (auto eyeiter = fevent.EyesBegin(fevt::ComponentSelector::eHasData);
         eyeiter != fevent.EyesEnd(fevt::ComponentSelector::eHasData); ++eyeiter) {

      // no rec data here, continue ...
      if (!eyeiter->HasRecData())
        continue;

      auto& eyerecdata = eyeiter->GetRecData();

      for (auto pix1 = eyerecdata.PulsedPixelsBegin();
           pix1 != eyerecdata.PulsedPixelsEnd(); ) {

        const auto& dir1 = det::Detector::GetInstance().GetFDetector().GetPixel(*pix1).GetDirection();
        bool isisolated = true;
        for (auto pix2 = eyerecdata.PulsedPixelsBegin();
             pix2 != eyerecdata.PulsedPixelsEnd(); ++pix2) {

          if (pix1 == pix2)
            continue;

          const auto& dir2 = det::Detector::GetInstance().GetFDetector().GetPixel(*pix2).GetDirection();

          if (Angle(dir1, dir2) < fMinAngle) {
            isisolated = false;
            break;
          }

        }
        if (isisolated)
          pix1 = eyerecdata.RemovePulsedPixel(pix1);
        else
          ++pix1;

      }

    }

    return eSuccess;
  }

}