/// Class to build average gas properties inside TPC
/// over a particular time period (or run).  Also
/// has facilities to make plots of TPC gas properties
/// over a time period.
/// Blair Jamieson (C) 2009, 2010
#include "ITPCGasAvg.hxx"

/// Initialize TPC gas properties averaging object using unix times
ITPCGasAvg::ITPCGasAvg( int maxsamples, long autmin,  long autmax){

  fTimeMin = autmin;
  fTimeMax = autmax;
  fMaxSamples = maxsamples;

  Init();
}
  
/// Initialize TPC gas properties averaging object using year,
/// month(0-11), day (1-31), ...
ITPCGasAvg:: ITPCGasAvg(int maxsamples,
				int ayearmin, int amonthmin, int adaymin, int ahourmin, int aminmin, int asecmin,
				int ayearmax, int amonthmax, int adaymax, int ahourmax, int aminmax, int asecmax ){
  fMaxSamples = maxsamples;

  tm mt;  
  // use time structure from time.h
  // to build unix time from year,month, etc.
  mt.tm_year = ayearmin - 1900;
  mt.tm_mon  = amonthmin-1;
  mt.tm_mday = adaymin;
  mt.tm_hour = ahourmin;
  mt.tm_min  = aminmin;
  mt.tm_sec  = asecmin;
  fTimeMin = (long)mktime(&mt);

  mt.tm_year = ayearmax - 1900;
  mt.tm_mon  = amonthmax-1;
  mt.tm_mday = adaymax;
  mt.tm_hour = ahourmax;
  mt.tm_min  = aminmax;
  mt.tm_sec  = asecmax;
  fTimeMax = (long)mktime(&mt);

  Init();
}


/// Initialize gas properties / pressure averaging object
void ITPCGasAvg::Init(){
  // first calculate delta-time for maximum number of samples:
  long adt;
  if ( fMaxSamples >0 ){
    adt = (fTimeMax-fTimeMin)/fMaxSamples;
    if (adt<30) adt = 30; // readings are only every so often anyway so set minimum to every 30 seconds
  } else {
    // use a default of every 30 sec?
    adt = 30; 
  }
  
  std::cout<<"<ITPCGasAvg::Init> Look for gas properties readings from time "
	   <<fTimeMin<<" to "<<fTimeMax<<std::endl;
  std::cout<<"<ITPCGasAvg::Init> MaxSamples = "<<fMaxSamples<<" (every "<<adt<<" seconds)"<<std::endl;

  // first make gas reading object
  ITPCGas aGas( fTimeMin, fTimeMax );
  long astart_time;
  long aend_time;
  aGas.GetTimeInterval( astart_time, aend_time );
  int nquery = 0;
  do {
    if (nquery%1 == 0 ){
      std::cout<<"DB Query"<<nquery
	       <<" time = "<<aGas.GetTimeStamp();
    }
    AddUpdate( aGas, false );
    aGas.ReadNext( adt );
    aGas.GetTimeInterval( astart_time, aend_time );
    nquery++;
  } while ( aend_time < fTimeMax );

  // add the final point
  AddUpdate( aGas, true );

  return;
}


/// Update the Average and RMS calculations for each of the gas system conditions
void ITPCGasAvg::AddUpdate( ITPCGas &agas, bool lastflag ){
  if ( lastflag==false && agas.TTPC1() < 0.1 ) return; // dont add bad reading.
  fArFlow .Add( agas.GetTime(),  agas.ArFlow()  ,lastflag );
  fCF4Flow.Add( agas.GetTime(), agas.CF4Flow() ,lastflag ); 
  fIbuFlow.Add( agas.GetTime(), agas.IbuFlow() ,lastflag ); 

  double losum = agas.ArFlow() + agas.CF4Flow() + agas.IbuFlow();
  fArFrac .Add(  agas.GetTime(), 100.0/losum* agas.ArFlow()  ,lastflag );
  fCF4Frac.Add(  agas.GetTime(), 100.0/losum* agas.CF4Flow() ,lastflag ); 
  fIbuFrac.Add(  agas.GetTime(), 100.0/losum* agas.IbuFlow() ,lastflag ); 
  fPatm   .Add(  agas.GetTime(), agas.Patm()*1000.0    ,lastflag );    
  fPD1    .Add(  agas.GetTime(), agas.PD1()     ,lastflag );     
  fPD2    .Add(  agas.GetTime(), agas.PD2()     ,lastflag );     
  fP      .Add(  agas.GetTime(), agas.P()       ,lastflag );     
  fPD4    .Add(  agas.GetTime(), agas.PD4()     ,lastflag );     
  	      
  fH2O    .Add(  agas.GetTime(), agas.H2O()     ,lastflag );     
  fO2     .Add(  agas.GetTime(), agas.O2()      ,lastflag );      
  fIbuMon .Add(  agas.GetTime(), agas.IbuMon()  ,lastflag );  
  fIbuAn  .Add(  agas.GetTime(), agas.IbuAn()   ,lastflag );   
  fCO2An  .Add(  agas.GetTime(), agas.CO2An()   ,lastflag );   
  fCF4An  .Add(  agas.GetTime(), agas.CF4An()   ,lastflag );   
  fTO2    .Add(  agas.GetTime(), agas.TO2()     ,lastflag );    
  fTTPC1  .Add(  agas.GetTime(), agas.TTPC1()   ,lastflag );    
  fTTPC2  .Add(  agas.GetTime(), agas.TTPC2()   ,lastflag );    
  fTTPC3  .Add(  agas.GetTime(), agas.TTPC3()   ,lastflag );    
  		 
  fDensity1  .Add(  agas.GetTime(), agas.P()/(agas.TTPC1()+273.15)   ,lastflag );    
  fDensity2  .Add(  agas.GetTime(), agas.P()/(agas.TTPC2()+273.15)   ,lastflag );    
  fDensity3  .Add(  agas.GetTime(), agas.P()/(agas.TTPC3()+273.15)   ,lastflag );    
  		 
  fD0TC1  .Add(  agas.GetTime(), agas.D0TC1() ,lastflag );    
  fD0TC2  .Add(  agas.GetTime(), agas.D0TC2() ,lastflag );    
  fD0TC3  .Add(  agas.GetTime(), agas.D0TC3() ,lastflag );    
  fD0TC4  .Add(  agas.GetTime(), agas.D0TC4() ,lastflag );    
  fD0TC5  .Add(  agas.GetTime(), agas.D0TC5() ,lastflag );    
  fP0TC1  .Add(  agas.GetTime(), agas.P0TC1() ,lastflag );    
  fP0TC2  .Add(  agas.GetTime(), agas.P0TC2() ,lastflag );    
  fP0TC3  .Add(  agas.GetTime(), agas.P0TC3() ,lastflag );    
  fP0TC4  .Add(  agas.GetTime(), agas.P0TC4() ,lastflag );    
  fP0TC5  .Add(  agas.GetTime(), agas.P0TC5() ,lastflag );    
  fP0TC6  .Add(  agas.GetTime(), agas.P0TC6() ,lastflag );    
  fDXTC1  .Add(  agas.GetTime(), agas.DXTC1() ,lastflag );    
  fDXTC2  .Add(  agas.GetTime(), agas.DXTC2() ,lastflag );    
  fDXTC3  .Add(  agas.GetTime(), agas.DXTC3() ,lastflag );    
  		 
  fVdIn    .Add( agas.GetTime(),  agas.VdIn()    ,lastflag );    
  fVdInUnc .Add( agas.GetTime(),  agas.VdInUnc() ,lastflag );    
  fVdOut   .Add( agas.GetTime(),  agas.VdOut()   ,lastflag );    
  fVdOutUnc.Add( agas.GetTime(),  agas.VdOutUnc(),lastflag );    
  
  fGainIn    .Add(  agas.GetTime(), agas.GainIn()    ,lastflag );    
  fGainInUnc .Add(  agas.GetTime(), agas.GainInUnc() ,lastflag );    
  fGainOut   .Add(  agas.GetTime(), agas.GainOut()   ,lastflag );    
  fGainOutUnc.Add(  agas.GetTime(), agas.GainOutUnc(),lastflag );    
  		  
  return;	  
}		  
		  
/// Print average and rms gas properties results to std::cout
void ITPCGasAvg::Print(){
  //(fTPAll[0])->PrintMap();

  std::cout<<"TPC Gas properties from "
	   << fTimeMin <<" to "
	   << fTimeMax << std::endl;
  fArFlow .Print("Ar  flow            (l/min):\t");  
  fCF4Flow.Print("CF4 flow            (l/min):\t");  
  fIbuFlow.Print("Ibu flow            (l/min):\t");

  fArFrac .Print("Ar  fraction            (%):\t");  
  fCF4Frac.Print("CF4 fraction            (%):\t");  
  fIbuFrac.Print("Ibu fraction            (%):\t");

  fPatm   .Print("Atmospheric Pressure (mbar):\t");
  fPD1    .Print("dP1 (gap-atmosphere)(mbarg):\t");
  fPD2    .Print("dP2 (inner-gap)     (mbarg):\t");
  fP      .Print("Inner Pressure       (mbar):\t");
  fPD4    .Print("dP4 (ibumon-atmos)  (mbarg):\t");

  fH2O   .Print("H2O                   (ppm):\t"  );
  fO2    .Print("O2                    (ppm):\t"  );
  fIbuMon.Print("Ibu Monitor             (%):\t"  );
  fIbuAn .Print("Ibu Analyzer (crap)     (%):\t"  );
  fCO2An .Print("CO2 Analyzer          (ppm):\t"  );
  fCF4An .Print("CF4 Analyzer            (%):\t"  );
  fTO2   .Print("O2 Analyzer Temp     (degC):\t"  );

  fTTPC1 .Print("TPC1 Avg Gas Man Temp(degC):\t"  );
  fTTPC2 .Print("TPC2 Avg Gas Man Temp(degC):\t"  );
  fTTPC3 .Print("TPC3 Avg Gas Man Temp(degC):\t"  );

  fTTPC1 .Print("TPC1 Gas Density   (mbar/K):\t"  );
  fTTPC2 .Print("TPC2 Gas Density   (mbar/K):\t"  );
  fTTPC3 .Print("TPC3 Gas Density   (mbar/K):\t"  );

  fD0TC1 .Print("Ibu flowmeter1 Temp    (degC):\t");  
  fD0TC2 .Print("Ibu flowmeter2 Temp    (degC):\t");  
  fD0TC3 .Print("Ibu purge supply Temp  (degC):\t");    
  fD0TC4 .Print("Ar supply line Temp    (degC):\t");    
  fD0TC5 .Print("Ibu norm supply Temp   (degC):\t");  
  									  
  fP0TC1 .Print("purifier 1 top Temp    (degC):\t");    
  fP0TC2 .Print("purifier 1 mid Temp    (degC):\t");    
  fP0TC3 .Print("purifier 1 bot Temp    (degC):\t");    
  fP0TC4 .Print("purifier 2 top Temp    (degC):\t");    
  fP0TC5 .Print("purifier 2 mid Temp    (degC):\t");    
  fP0TC6 .Print("purifier 2 bot Temp    (degC):\t");    
  									  
  fDXTC1 .Print("Ibu bottle 1 Temp      (degC):\t");    
  fDXTC2 .Print("Ibu bottle 2 Temp      (degC):\t");    
  fDXTC3 .Print("Ibu line Temp          (degC):\t");    

  fVdIn      .Print("Drift speed input  (cm/mus):\t");
  fVdInUnc   .Print("Drift speed in unc.(cm/mus):\t");
  fVdOut     .Print("Drift speed return (cm/mus):\t");
  fVdOutUnc  .Print("Drift speed ret unc(cm/mus):\t");

  fGainIn    .Print("Gain of input gas          :\t");
  fGainInUnc .Print("Unc. in Gain of input gas  :\t");
  fGainOut   .Print("Gain or return gas         :\t");
  fGainOutUnc.Print("Unc. in Gain of return gas :\t");

}

/// Plot gas properties
void ITPCGasAvg::Plot(){

  std::cout<<"TPC Gas properties from "
	   << fTimeMin <<" to "
	   << fTimeMax << std::endl;
  fArFlow .MakePlot("Arflow",false);  
  fCF4Flow.MakePlot("CF4flow",false);  
  fIbuFlow.MakePlot("Ibuflow",false);

  fArFrac .MakePlot("Arfraction",false);  
  fCF4Frac.MakePlot("CF4fraction",false);  
  fIbuFrac.MakePlot("Ibufraction",false);

  fPatm   .MakePlot("Patm",false);
  fPD1    .MakePlot("dP1",false);
  fPD2    .MakePlot("dP2",false);
  fP      .MakePlot("P",false);
  fPD4    .MakePlot("dP4",false);

  fH2O   .MakePlot("H2Oppm",false);
  fO2    .MakePlot("O2ppm",false);
  fIbuMon.MakePlot("IbuMon",false);

  fCO2An .MakePlot("CO2ppm",false);
  fCF4An .MakePlot("CF4",false);
  fTO2   .MakePlot("TO2degC",false);

  fTTPC1 .MakePlot("TTPC1degC",false);
  fTTPC2 .MakePlot("TTPC2degC",false);
  fTTPC3 .MakePlot("TTPC3degC",false);

  fDensity1 .MakePlot("TTPC1density",false);
  fDensity2 .MakePlot("TTPC2density",false);
  fDensity3 .MakePlot("TTPC3density",false);

  fD0TC1 .MakePlot("D0TC1",false);  
  fD0TC2 .MakePlot("D0TC2",false);  
  fD0TC3 .MakePlot("fD0TC3",false);    
  fD0TC4 .MakePlot("fD0TC4",false);    
  fD0TC5 .MakePlot("fD0TC5",false);  
  		    				  
  fP0TC1 .MakePlot("fP0TC1",false);    
  fP0TC2 .MakePlot("fP0TC2",false);    
  fP0TC3 .MakePlot("fP0TC3",false);    
  fP0TC4 .MakePlot("fP0TC4",false);    
  fP0TC5 .MakePlot("fP0TC5",false);    
  fP0TC6 .MakePlot("fP0TC6",false);    
  		    				  
  fDXTC1 .MakePlot("fDXTC1",false);    
  fDXTC2 .MakePlot("fDXTC2",false);    
  fDXTC3 .MakePlot("fDXTC3",false);    

  fVdIn      .MakePlot("VdIn",false);
  fVdInUnc   .MakePlot("fVdInUnc",false);
  fVdOut     .MakePlot("fVdOut",false);
  fVdOutUnc  .MakePlot("fVdOutUnc",false);
				     
  fGainIn    .MakePlot("fGainIn",false);
  fGainInUnc .MakePlot("fGainInUnc",false);
  fGainOut   .MakePlot("fGainOut",false);
  fGainOutUnc.MakePlot("fGainOutUnc",false);

  char aname[256];
  sprintf(aname,"ttpcgasavg_%ld_to_%ld.root",fTimeMin,fTimeMax);
  TFile fout(aname,"recreate");
  fArFlow .fGraph->Write();
  fCF4Flow.fGraph->Write();
  fIbuFlow.fGraph->Write();
  	     
  fArFrac .fGraph->Write();
  fCF4Frac.fGraph->Write();
  fIbuFrac.fGraph->Write();
  
  fPatm   .fGraph->Write();
  fPD1    .fGraph->Write();
  fPD2    .fGraph->Write();
  fP      .fGraph->Write();
  fPD4    .fGraph->Write();
  
  fH2O   .fGraph->Write();
  fO2    .fGraph->Write();
  fIbuMon.fGraph->Write();
  
  fCO2An .fGraph->Write();
  fCF4An .fGraph->Write();
  fTO2   .fGraph->Write();
  
  fTTPC1 .fGraph->Write();
  fTTPC2 .fGraph->Write();
  fTTPC3 .fGraph->Write();
  
  fDensity1 .fGraph->Write();
  fDensity2 .fGraph->Write();
  fDensity3 .fGraph->Write();
  
  fD0TC1 .fGraph->Write();
  fD0TC2 .fGraph->Write();
  fD0TC3 .fGraph->Write();
  fD0TC4 .fGraph->Write();
  fD0TC5 .fGraph->Write();
  
  fP0TC1 .fGraph->Write();
  fP0TC2 .fGraph->Write();
  fP0TC3 .fGraph->Write();
  fP0TC4 .fGraph->Write();
  fP0TC5 .fGraph->Write();
  fP0TC6 .fGraph->Write();
  
  fDXTC1 .fGraph->Write();
  fDXTC2 .fGraph->Write();
  fDXTC3 .fGraph->Write();
  
  fVdIn      .fGraph->Write();
  fVdInUnc   .fGraph->Write();
  fVdOut     .fGraph->Write();
  fVdOutUnc  .fGraph->Write();
  
  fGainIn    .fGraph->Write();
  fGainInUnc .fGraph->Write();
  fGainOut   .fGraph->Write();
  fGainOutUnc.fGraph->Write();
  fout.Close();
  return;
}