#include <RAT/NhitCut.hh>
#include <RAT/DU/Utility.hh>
#include <RAT/DU/PMTInfo.hh>

using namespace std;

namespace RAT {
  void NhitCut::BeginOfRun(DS::Run&)
  {
    fLClean = DB::Get()->GetLink("DATACLEANING",fName);
    fNhitValue = fLClean->GetI("nhit_value");
  }

  Processor::Result NhitCut::DSEvent(DS::Run&, DS::Entry& ds)
  {
    bool pass = true; // If any triggered event fails, they all fail
    for( size_t iEV = 0; iEV < ds.GetEVCount(); iEV++ )
      if( Event(ds, ds.GetEV(iEV)) != OKTRUE )
        pass = false;
    return pass ? OKTRUE : OKFALSE;
  }

  Processor::Result NhitCut::Event(DS::Entry&, DS::EV& ev)
  {

    fPassFlag = true;

    // fail any event that has multiple entries for the same pmt
    int nhit = ev.GetNhits();
    if (nhit<fNhitValue)
       fPassFlag=false;

    // now update the mask
    UpdateMask(ev);

    return fPassFlag ? OKFALSE : OKTRUE;
  }

} // namespace RAT