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

#ifndef TTPCGasDensity_hxx__
#define TTPCGasDensity_hxx__

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

/// Class to get a decent density measurement
/// at a particular time.  Uses averaging over
/// several measurements to help reduce any
/// noise effects.
/// Blair Jamieson (C) 2009, 2010
class ITPCGasDensity {
  
public:
  ITPCGasDensity( long autmin, long autmax, int dtaverage = 3600 );
  
  ITPCGasDensity( 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,  int itpc ){ 
    if ( atime > fTimeMin && atime < fTimeMax &&
	 itpc>=0 && itpc<3 ){
      for (unsigned int i=0; i<fTimes.size(); i++){
	if ( fTimes[i] > atime ) {
	  if (itpc==0) return fT1Avg[i].fmean;
	  if (itpc==1) return fT2Avg[i].fmean;
	  if (itpc==2) return fT3Avg[i].fmean;
	}
      }
    }
    std::cout<<"ITPCGasDensity::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<<"ITPCGasDensity::P> time outside of range"<<std::endl;
    return -1.0;
  }
 
  /// Get the averaged density at atime for itpc
  /// in units of mbar/K
  float Density( long atime,  int itpc ){ 
    float aP = P(atime);
    float aT = T(atime,itpc);
    if (aP>0.0&&aT>0.0){
      return aP/(aT+273.15);
    } else {
      std::cout<<"ITPCGasDensity::Density> 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 > fT1Avg;
  std::vector< IAvgRMSPlot > fT2Avg;
  std::vector< IAvgRMSPlot > fT3Avg;
  std::vector< IAvgRMSPlot > fDens1Avg;
  std::vector< IAvgRMSPlot > fDens2Avg;
  std::vector< IAvgRMSPlot > fDens3Avg;
  std::vector< IAvgRMSPlot > fPAvg;
  
  // Initialize values for later lookup
  void Init();

  // Make some plots of the averaged quantities during initialization
  IAvgRMSPlot fT[3];       ///< Tpc temperature (degC)
  IAvgRMSPlot fP  ;        ///< Pressure of inner volume (mbar)
  IAvgRMSPlot fDensity[3]; ///< Density  (mbar/K)

};

#endif