QualityCut.h

/* -*- C++ -*-.*********************************************************************************************
 Author: Ben Strutt
 Email: strutt@physics.ucla.edu

 Description:
             My ANITA-3 data quality cuts, to remove bad events
***********************************************************************************************************/

#ifndef QUALITYCUTS_H
#define QUALITYCUTS_H

#include "TObject.h"
#include "TString.h"

#include "FilteredAnitaEvent.h"
#include "AnalysisWaveform.h"
#include "RootTools.h"

#include "AnitaEventSummary.h"

namespace Acclaim{


  class QualityCut : public TObject {

  public:
    static Bool_t applyAll(const UsefulAnitaEvent* usefulEvent, AnitaEventSummary* sum=NULL); // static utility function, applies all cuts defined here
    static Bool_t passedAll(const AnitaEventSummary* sum, bool describe = false); // static utility function, reads flags in AnitaEventSummary, returns true if event passed
    

    
    QualityCut(){;}
    virtual ~QualityCut(){;}
    // virtual void apply(FilteredAnitaEvent* fEv) = 0;
    virtual void apply(const UsefulAnitaEvent* useful, AnitaEventSummary* sum = NULL) = 0;    
    TString description;
    Bool_t eventPassesCut;
    ClassDef(QualityCut, 1)
  };



  class SurfSaturationCut : public QualityCut {
    ClassDef(Acclaim::SurfSaturationCut, 1);

  protected:
    double maxLimit;
    double minLimit;
    double asymLimit;

    double maxVolts;
    int maxVoltsAnt;
    AnitaPol::AnitaPol_t maxVoltsPol;

    double minVolts;
    int minVoltsAnt;
    AnitaPol::AnitaPol_t minVoltsPol;
    
    double asymVolts;
    int asymVoltsAnt;
    AnitaPol::AnitaPol_t asymVoltsPol;

  public:
    SurfSaturationCut();
    virtual ~SurfSaturationCut(){;}    
    // virtual void apply(FilteredAnitaEvent* fEv);
    virtual void apply(const UsefulAnitaEvent* useful, AnitaEventSummary* sum = NULL);    
  };




  
  class PayloadBlastCut : public QualityCut {
    ClassDef(Acclaim::PayloadBlastCut, 2);
    
  protected:
    double ratioCutHigh;
    double ratioCutLow;


    double maxRatio;
    int maxRatioPhi;
    AnitaPol::AnitaPol_t maxRatioPol;

  public:
    PayloadBlastCut();   
    virtual ~PayloadBlastCut(){;}
    // virtual void apply(FilteredAnitaEvent* fEv);
    virtual void apply(const UsefulAnitaEvent* useful, AnitaEventSummary* sum = NULL);
  };
  


  class NumPointsCut : public QualityCut {
    ClassDef(Acclaim::NumPointsCut, 1);
    
  protected:
    int numPointsCutLow; // presumably don't need a cut on the other side of this...
    int minNumPoints;

  public:
    NumPointsCut();
    virtual ~NumPointsCut(){;}
    // virtual void apply(FilteredAnitaEvent* fEv);
    virtual void apply(const UsefulAnitaEvent* useful, AnitaEventSummary* sum = NULL);
  };
  
}

#endif //ANALYSISCUTS