DataQualityMonitor.cxx

#include "DataQualityMonitor.h"

DataQualityMonitor::DataQualityMonitor(){

  Double_t defaultThresholds[NUM_POL][NUM_SEAVEYS] = {
    {109, 104, 139, 102, 173, 104, 116, 107, 85, 186, 224, 157, 245, 118, 137, 108,
     95, 117, 108, 111, 109, 111, 103, 103, 97, 199, 208, 200, 184, 205, 175, 109,
     106, 147, 117, 115, 117, 116, 114, 106, 95, 179, 208, 189, 184, 198, 110, 99},
    {186, 205, 235, 269, 180, 239, 247, 175, 197, 310, 288, 327, 291, 270, 374, 232,
     247, 257, 264, 249, 218, 283, 216, 250, 220, 360, 293, 347, 375, 423, 373, 247,
     114, 177, 201, 266, 207, 239, 206, 184, 190, 327, 326, 271, 397, 415, 247, 255}};
  
  fillThresholds(defaultThresholds);  
}


DataQualityMonitor::~DataQualityMonitor(){
  
}


Int_t DataQualityMonitor::ProcessEvent(RawAnitaEvent* rawEvent){

  UsefulAnitaEvent* pedEvent = new UsefulAnitaEvent(rawEvent, WaveCalType::kAddPeds);
  pedEvent->setAlfaFilterFlag(false); // Turn off the ALFA filtering
  
  fillMaxAbsSecondDerivArray(pedEvent);

  Int_t numBadChannels = 0;

  for(Int_t pol=0; pol<NUM_POL; pol++){
    for(Int_t ant=0; ant<NUM_SEAVEYS; ant++){

      if(lastMaxAbsSecondDeriv[pol][ant] > threshMaxAbsSecondDeriv[pol][ant]){
    numBadChannels++;
      }
    }
  }
  
  delete pedEvent;

  return numBadChannels;
}


void DataQualityMonitor::fillMaxAbsSecondDerivArray(UsefulAnitaEvent* pedEvent){

  for(Int_t pol=0; pol<NUM_POL; pol++){
    for(Int_t ant=0; ant<NUM_SEAVEYS; ant++){
      lastMaxAbsSecondDeriv[pol][ant] = -9999;

      TGraph* gr = pedEvent->getGraph(ant, (AnitaPol::AnitaPol_t) pol);

      for(int samp = 0; samp < gr->GetN(); samp++){

    Double_t V = gr->GetY()[samp];

    if(samp < gr->GetN() - 2){
      Double_t V2 = gr->GetY()[samp+1];
      Double_t V3 = gr->GetY()[samp+2];
      double thisAbsSecondDeriv = TMath::Abs(2*V2 - V3 - V);
      if(thisAbsSecondDeriv > lastMaxAbsSecondDeriv[pol][ant]){
        lastMaxAbsSecondDeriv[pol][ant] = thisAbsSecondDeriv;
      }
    }

      }

      delete gr;

    }
  }
}


void DataQualityMonitor::fillThresholds(Double_t vals[NUM_POL][NUM_SEAVEYS]){
  
  for(Int_t pol=0; pol<NUM_POL; pol++){
    for(Int_t ant=0; ant<NUM_SEAVEYS; ant++){
      threshMaxAbsSecondDeriv[pol][ant] = vals[pol][ant];
    }
  }
  
}