/// Class to get a decent gas analysis measurements
/// at a particular time.  Uses averaging over
/// several measurements to help reduce any
/// noise effects.
/// Blair Jamieson (C) 2009, 2010

#ifndef TTPCGasSmooth_hxx__
#define TTPCGasSmooth_hxx__

#include "ITPCGas.hxx"
#include "IAvgRMSPlot.hxx"
#include <vector>
#include <TFile.h>

/// Class to get a decent gas analysis measurement
/// at a particular time.  Uses averaging over
/// several measurements to help reduce any
/// noise effects.
/// Blair Jamieson (C) 2009, 2010
class ITPCGasSmooth {
  
public:
  ITPCGasSmooth( long autmin, long autmax, int dtaverage = 3600 );
  
  ITPCGasSmooth( int ayearmin, int amonthmin, int adaymin, int ahourmin, int aminmin, int asecmin,
     int ayearmax, int amonthmax, int adaymax, int ahourmax, int aminmax, int asecmax, 
     int dtaverage = 3600);
  
  /// Get the averaged temperature at atime for itpc
  float T( long atime ){ 
    if ( atime > fTimeMin && atime < fTimeMax  ){
      for (unsigned int i=0; i<fTimes.size(); i++){
  if ( fTimes[i] > atime ) return fTAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::TTPC> time or tpc outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the pressure at atime 
  float P( long atime ){ 
    if ( atime > fTimeMin && atime < fTimeMax ){
      for (unsigned int i=0; i<fPAvg.size(); i++){
  if ( fTimes[i] > atime ) {
    return fPAvg[i].fmean;
  }
      }
    }
    std::cout<<"ITPCGasSmooth::P> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged density at atime
  /// in units of mbar/K
  float Density( long atime ){ 
    if ( atime > fTimeMin && atime < fTimeMax ){
      for (unsigned int i=0; i<fDensAvg.size(); i++){
  if ( fTimes[i] > atime ) {
    return fDensAvg[i].fmean;
  }
      }
    }
    std::cout<<"ITPCGasSmooth::Density> time outside of range"<<std::endl;
    return -1.0;

  }
  
  /// Get the averaged H2O reading (ppm)
  float H2O( long atime ){ 
    for (unsigned int i=0; i<fH2OAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fH2OAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::H2O> time outside of range"<<std::endl;
    return -1.0;
  }
  
  
  /// Get the averaged O2 reading (ppm)
  float O2( long atime ){ 
    for (unsigned int i=0; i<fO2Avg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fO2Avg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::O2> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged IbuMon reading (%)
  float IbuMon( long atime ){ 
    for (unsigned int i=0; i<fIbuMonAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fIbuMonAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::IbuMon> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged CO2 reading (ppm)
  float CO2( long atime ){ 
    for (unsigned int i=0; i<fCO2Avg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fCO2Avg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::CO2> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged CF4 reading (ppm)
  float CF4( long atime ){ 
    for (unsigned int i=0; i<fCF4Avg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fCF4Avg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::CF4> time outside of range"<<std::endl;
    return -1.0;
  }

  /// Get the averaged VdIn reading (cm/us)
  float VdIn( long atime ){ 
    for (unsigned int i=0; i<fVdInAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fVdInAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::VdIn> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged VdOut reading (cm/us)
  float VdOut( long atime ){ 
    for (unsigned int i=0; i<fVdOutAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fVdOutAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::VdOut> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged GainIn reading ()
  float GainIn( long atime ){ 
    for (unsigned int i=0; i<fGainInAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fGainInAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::GainIn> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged GainOut reading ()
  float GainOut( long atime ){ 
    for (unsigned int i=0; i<fGainOutAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fGainOutAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::GainOut> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged Argon Fraction reading ()
  float ArFrac( long atime ){ 
    for (unsigned int i=0; i<fArFracAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fArFracAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::ArFrac> time outside of range"<<std::endl;
    return -1.0;
  }

  /// Get the averaged CF4 fraction reading ()
  float CF4Frac( long atime ){ 
    for (unsigned int i=0; i<fCF4FracAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fCF4FracAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::CF4Frac> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// Get the averaged Isobutane Fraction reading ()
  float IbuFrac( long atime ){ 
    for (unsigned int i=0; i<fIbuFracAvg.size(); i++){
      if ( fTimes[i] > atime ) {
  return fIbuFracAvg[i].fmean;
      }
    }
    std::cout<<"ITPCGasSmooth::IbuFrac> time outside of range"<<std::endl;
    return -1.0;
  }
  
  /// plot the averaged quantities
  void Plot();
  
  
  /// Add to average and rms objects
  void AddUpdate( ITPCGas &agas, bool lastflag );
  
  
private:
  
  
  long fTimeMin;  ///<minimum time reading
  long fTimeMax;  ///<maximum time reading
  int  fDtAvg;    ///<length of time in seconds to average over
  unsigned int  fCurIdx;   ///<index to maps of first unfinished avgrmsplot 
  
  // map of unix time to averaged gas quantities
  std::vector< long > fTimes;
  std::vector< IAvgRMSPlot > fTAvg       ;  ///< O2 Analyzer temperature
  std::vector< IAvgRMSPlot > fPAvg       ;  ///< Atmospheric pressure
  std::vector< IAvgRMSPlot > fDensAvg    ;  ///< Density (mbar/K)
  std::vector< IAvgRMSPlot > fH2OAvg     ;  ///< H2O in ~ppm
  std::vector< IAvgRMSPlot > fO2Avg      ;  ///< Inner volume O2 in ppm
  std::vector< IAvgRMSPlot > fIbuMonAvg  ;  ///< Isobutane % from safety monitor
  std::vector< IAvgRMSPlot > fCO2Avg     ;  ///< CO2 ppm from MGA3000 (crap)
  std::vector< IAvgRMSPlot > fCF4Avg     ;  ///< CF4 % from MGA3000 (fine but offset)
  std::vector< IAvgRMSPlot > fVdInAvg    ;  ///< Drift speed of input gas ( cm / microseconds )
  std::vector< IAvgRMSPlot > fVdOutAvg   ;  ///< Drift speed of return gas ( cm / microseconds )
  std::vector< IAvgRMSPlot > fGainInAvg  ;  ///< Gain of input gas 
  std::vector< IAvgRMSPlot > fGainOutAvg ;  ///< Drift speed of return gas
  std::vector< IAvgRMSPlot > fArFracAvg  ;  ///< Argon flow ()
  std::vector< IAvgRMSPlot > fCF4FracAvg ;  ///< CF4 flow ()
  std::vector< IAvgRMSPlot > fIbuFracAvg ;  ///< isobutane flow ()
  
  // Initialize values for later lookup
  void Init();
  
  // Make some plots of the averaged quantities during initialization
  IAvgRMSPlot fT   ;     ///< Tpc temperature (degC)
  IAvgRMSPlot fP   ;     ///< Pressure of inner volume (mbar)
  IAvgRMSPlot fDens;     ///< Density  (mbar/K)
  IAvgRMSPlot fH2O    ;  ///< H2O in ~ppm
  IAvgRMSPlot fO2     ;  ///< Inner volume O2 in ppm
  IAvgRMSPlot fIbuMon ;  ///< Isobutane % from safety monitor
  IAvgRMSPlot fCO2    ;  ///< CO2 ppm from MGA3000 (crap)
  IAvgRMSPlot fCF4    ;  ///< CF4 % from MGA3000 (fine but offset)
  IAvgRMSPlot fVdIn   ;  ///< Drift speed of input gas ( cm / microseconds )
  IAvgRMSPlot fVdOut  ;  ///< Drift speed of return gas ( cm / microseconds )
  IAvgRMSPlot fGainIn ;  ///< Gain of input gas 
  IAvgRMSPlot fGainOut;  ///< Drift speed of return gas
  IAvgRMSPlot fArFrac ;  ///< Argon flow ()
  IAvgRMSPlot fCF4Frac;  ///< CF4 flow ()
  IAvgRMSPlot fIbuFrac;  ///< isobutane flow ()
  
};

#endif