BandwidthMeasure.h

#ifndef BANDWIDTH_MEASURE_H
#define BANDWIDTH_MEASURE_H

#include <vector>


class AnalysisWaveform; 
class TGraph;
class TGraphAligned; 

namespace bandwidth 
{
  double bandwidthMeasure(const AnalysisWaveform * wf, int timeCheck = 0, TGraph* gTest = 0);
  double differenceFromImpulse(const AnalysisWaveform * wf, int timeCheck = 0, TGraph* gTest = 0);
  double maxDifferenceFromImpulse(const AnalysisWaveform * wf, int timeCheck = 0, TGraph* gTest = 0);
  double hooverIndex(const AnalysisWaveform * wf, int timeCheck = 0);
  double theilIndex(const AnalysisWaveform * wf, int timeCheck = 0);
  double powerInBand(const AnalysisWaveform * wf, double minFreq=.17, double maxFreq=.26);


  double lowness(const AnalysisWaveform *wf, double minf = 0.17, double maxf = 0.6) ; 


  void checkNotches(int timeCheck, double& notch0, double& notch1, double& notch2);
  double fillPowers(const TGraphAligned* powd, std::vector<double> &powers, double notch0, double notch1, double notch2);
  TGraph* loadImpulsePower(int timeCheck);
  void normalizePower(TGraph* g);
  TGraph* downsampleImpulse(TGraph* imp, const TGraphAligned* examp);

}




#endif