FlasherFinder.cxx

Go to the documentation of this file.

// C++ includes
#include <vector>
#include <algorithm>
#include <iostream>

// NOvASoft includes
#include "DataQuality/FlasherFinder/FlasherFinder.h"
#include "RawData/RawDigit.h"
#include "RecoBase/CellHit.h"
#include "RecoBase/Cluster.h"

using namespace dq;


//......................................................................
static bool gsFEBTimeSort(const FEBHit& h1, const FEBHit& h2)
{
  if (h1.fFEB == h2.fFEB) return (h1.fTDC < h2.fTDC);
  return h1.fFEB < h2.fFEB;
}


//......................................................................
FlasherFinder::FlasherFinder():
  fInstigatorQ(10000), //12-10,000 for FarDet, ~30,000 for NearDet, based on DocDB-13008, slides 21,22
  fInstigatorT(64),
  fVetoT(960), //640 (10 us), 960 (15 us), 1280 (20 us)
  fVetoQ(500)
{ }


//......................................................................
unsigned int FlasherFinder::FindInstigators(std::vector<art::Ptr<rawdata::RawDigit> >& d)
{
  //
  // Extract the information I need from the digit list then sort by
  // FEB and time
  //

  unsigned int i, j;
  std::vector<FEBHit> febhit(d.size());
  for (i=0; i<d.size(); ++i) {
    febhit[i].fFEB = d[i]->DaqChannel() & 0xFFFFFF00; ///< Pixel is low 8 bits
    febhit[i].fTDC = d[i]->TDC();
    febhit[i].fADC = d[i]->ADC();
    febhit[i].fDAQ = d[i]->DaqChannel();
  }

  //
  // Search for time windows where summed ADC crosses threshold
  //
  sort(febhit.begin(), febhit.end(), gsFEBTimeSort);
  for (i=0; i<febhit.size(); ++i) {
    unsigned int sumADC = 0;
    for (j=i+1; j<febhit.size(); ++j) {
      bool kInSameFEB = (febhit[j].fFEB == febhit[i].fFEB);
      bool kInTime    = ((febhit[j].fTDC - febhit[i].fTDC) < fInstigatorT);

      if (!kInSameFEB || !kInTime) break;
      sumADC += febhit[j].fADC;
    }

    //
    // Check the ADC to see if it crossed threshold
    //
    if (sumADC>fInstigatorQ) {
      fInstigators.push_back(febhit[i]);
    }
  }
  fInstigatorsOK = true;
  return fInstigators.size();
}


//......................................................................
bool FlasherFinder::IsInstigator(const art::Ptr<rawdata::RawDigit>& d) const
{
  assert(fInstigatorsOK);
  unsigned int j;

  int          DAQ = d->DaqChannel();
  unsigned int TDC = d->TDC();

  for (j=0; j<fInstigators.size(); j++) {
    if (fInstigators[j].fDAQ == DAQ && fInstigators[j].fTDC == TDC)
      return true;
  }

  return false;
}


//......................................................................
bool FlasherFinder::IsInstigator(const art::Ptr<rb::CellHit>& ch) const
{
  assert(fInstigatorsOK);
  unsigned int j;

  int          DAQ = ch->DaqChannel();
  unsigned int TDC = ch->TDC();

  for (j=0; j<fInstigators.size(); j++) {
    if (fInstigators[j].fDAQ == DAQ && fInstigators[j].fTDC == TDC)
      return true;
  }

  return false;
}


//......................................................................
bool FlasherFinder::Veto(const art::Ptr<rawdata::RawDigit>& d) const
{
  if ((unsigned int)d->ADC() > fVetoQ) return false;

  assert(fInstigatorsOK);
  unsigned int j;

  unsigned int feb = d->DaqChannel() & 0xFFFFFF00;
  unsigned int tdc = d->TDC();

  for (j=0; j<fInstigators.size(); ++j) {
    bool kInSameFEB = (fInstigators[j].fFEB == feb);
    bool kInTime    = (tdc > fInstigators[j].fTDC && tdc < fInstigators[j].fTDC + fVetoT);

    if (kInSameFEB && kInTime) return true;
  }

  return false;
}


//......................................................................
bool FlasherFinder::Veto(const art::Ptr<rb::CellHit>& ch) const
{
  if ((unsigned int)ch->ADC() > fVetoQ) return false;

  assert(fInstigatorsOK);
  unsigned int j;

  unsigned int feb = ch->DaqChannel()&0xFFFFFF00;
  unsigned int tdc = ch->TDC();

  for (j=0; j<fInstigators.size(); ++j) {
    bool kInSameFEB = (fInstigators[j].fFEB == feb);
    bool kInTime    = (tdc > fInstigators[j].fTDC && tdc < fInstigators[j].fTDC + fVetoT);

    if (kInSameFEB && kInTime) return true;
  }

  return false;
}


//......................................................................
unsigned int FlasherFinder::Filter(std::vector<art::Ptr<rawdata::RawDigit>>& d)
{
  std::vector<art::Ptr<rawdata::RawDigit>> filtered;
  unsigned int i;
  unsigned int nfilt = 0;

  for (i=0; i<d.size(); ++i) {
    if (this->Veto(d[i])) {
      ++nfilt;
    }
    else {
      filtered.push_back(d[i]);
    }
  }

  swap(d, filtered);
  return nfilt;
}


//......................................................................
std::vector<art::Ptr<rawdata::RawDigit>> FlasherFinder::FEBFlashHits(std::vector<art::Ptr<rawdata::RawDigit>>& d)
{
  std::vector<art::Ptr<rawdata::RawDigit>> filtered;
  for (unsigned int i=0; i<d.size(); ++i) {
    if (this->Veto(d[i])) {
      filtered.push_back(d[i]);
    }
  }

  return filtered;
}


//......................................................................
unsigned int FlasherFinder::Nflash(const art::Ptr<rb::Cluster>& c) const
{
  unsigned int i;
  unsigned int nhit = 0;

  art::PtrVector<rb::CellHit> slicehits = c->AllCells();
  for (i=0; i<slicehits.size(); ++i) {
    if (this->Veto(slicehits[i])) ++nhit;
  }

  return nhit;
}