#ifndef _LEVEL1DATA_H_
#define _LEVEL1DATA_H_

#include <string>
#include <vector>


#include "DatabaseHandles.h"

typedef std::vector<std::vector<double> > Vector2D;

class CMarkup;

class Level1Data
{
public:
	Level1Data() {}
	virtual ~Level1Data() {}
	
	
	void Init(CMarkup* xml, const std::string& section="level1");

        bool AverageList(const std::vector<int>& events, char* l1_ver) ;

    bool Read(const int EventNo, char* l1_ver );


	bool ReadSummary(const int  EventNo, char* l1_ver); 
			
	void ReadScans(const int  EventNo);

//	void ReadSignals(const int  EventNo);
	void ReadGridded(const int  EventNo);

	void ReadZPT(const int  EventNo);

	void ReadEphem(const int EventNo);
		
//	void ReadLockDown(const int  EventNo);
	
	void ReadDriftData(const int EventNo);
	void ReadMSIS(const int EventNo) ;

//	void ReadSolarExtent(const int EventNo);
		
	int EventNo() const { return _EventNo; }
	int OrbitNo() const { return _OrbitNo; }
	std::string Mode() const { return _Mode; }
	std::string StartTime() const { return _StartTime; }
	std::string EndTime() const { return _EndTime; }
	std::vector<double> ScanElev() const { return _ElevAng; }
	std::vector<double> ScanData(const int band) const {return _ScanData.at(band-1); }
    std::vector<double> GetTransmittance(const int band) const { return _DetectorData.at(band-1) ; }
    //std::vector<double> GetSignal(const int band) const { return _DetData.at(band-1) ; }
    // RGS Changere,  with signal averages,  we want to use Transmissions rather than signal for computed transmissions (see poor man FOV for instance)
    //  force the 'GetSignal' function to return Transmission * vexo
    std::vector<double> GetSignal(const int band) const { 
       std::vector<double>  CurrentSignal = _DetectorData.at(band-1);
       double  Currentexo    = _exoV.at(band-1);
       for(int i = 0; i<CurrentSignal.size();i++){CurrentSignal[i] = CurrentSignal[i]*Currentexo;}
       return  CurrentSignal ; 
       }

    void SetSignal(const int band, const std::vector<double>& v ) { _DetData[band-1] = v; } 

	std::vector<double> GetImpactAlts() { return _TPAlt; } 
	std::vector<double> GetViewAngle() { return _ViewAngle; }
	std::vector<double> LockDownElev() { return _LockDownElev; }
	std::vector<double> SolarExtent() { return _SolarExtent; }
	std::vector<double> SunSensorRefraction() { return _SSRefraction; }
std::vector<double> SunSensorTransmission(const int iregion) { return _SSTransmission.at(iregion) ; }
	std::vector<double> GetTimes() { return _Times; }
	std::vector<float> Alt() { return _Altitude; }
	std::vector<float> Pres() {return _Pressure; }
	std::vector<float> Temp() { return _Temperature; }
std::vector<float> msisAlt() { return _msisAlt; }
std::vector<float> msisO() { return _msisO; }
std::vector<float> msisO2() { return _msisO2; }
std::vector<float> msisN2() { return _msisN2; }

	double DriftVariance( const int band) const { return _noise.at(band-1) ;} 
	double GetExoCounts( const int band) const { return _exoV.at(band-1); }
 	double DetectorTemp(const int band) const {return _FilterTemps.at(band-1); }
	double Attenuator(const int band) const { return _Attenuators.at(band-1); }
    
    double ThermalChiSq(const int band) const { return _ThermalChiSq.at(band-1); }

	float RadiusCurvature() const { return _CurvatureRadius; }  // need to get from database
	float EarthSunDistance() const { return _EarthSunDistance; }
	float LOS_Bearing() const { return _LOS_Bearing; }
	float Lat83() { return _Lat; }
	float Lon83() { return _Lon; }
	float RegOffset() { return _RegOffset; }
//
	void SetLockDownElev(const std::vector<double>& v) { _LockDownElev = v; }
    void Set_TP_Alts(const std::vector<double>& v) { _TPAlt = v; }

protected:
// from summary info
	int _EventNo, _OrbitNo;
        
//	std::string _topIndex, _lastIndex;
	std::string _StartTime, _EndTime;   // mysql format  yyyy-mm-dd hh:mm:ss
	std::string _Mode;   //r or s
// from read scan data
	std::vector<double> _ElevAng;   // Solar Scan Elevation Data
	Vector2D _ScanData;			  // Solar to Limb function for bands 1-16
// All the signals on an  gridded altitude
	std::vector<double> _Times; // seconds since _StartTime
	Vector2D _DetectorData;	  // bands   1-16  Transmission and difference signals (17-24)
    Vector2D _DetData;
    Vector2D _SSTransmission;
	std::vector<double> _exoV; // bands 1-16 exo values  (counts)
	std::vector<double> _noise; // noise values for signals std dev
// The Z, P, T Profile
	std::vector<float> _Altitude, _Temperature, _Pressure,
	_msisAlt,
	_msisO,
	_msisO2,
	_msisN2;
// Ephemeris Data  Tangent Point lats and longs Doppler shifts etc.
	std::vector<double> 
		_ViewAngle,
		_TPAlt,
		_TPLat,
		_TPLon,
		_SSRefraction;
      // _SSTransmission;
;
//  radius of curvature
	float
		_CurvatureRadius,
		_EarthSunDistance,
		_LOS_Bearing; 
// location interpolated to 83km   Needs implemented
	float _Lat, _Lon,
		_RegOffset;
// Filter Temperatures
	std::vector<double> _FilterTemps;  // (K) bands 1-16 filter temperatures
// attenuator settings
	std::vector<double> _Attenuators;  //band 1-16
    std::vector<double> _ThermalChiSq;    // NO Thermal Correction
	
// Lockdown of Band 3
	std::vector<double> _LockDownElev;   // radians
	std::vector<double> _LockDownAz;	//Radians
	std::vector<double> _SolarExtent;  //radians
// database info
	gatsDBpp::GATS_DBptr _dbconn;	
	bool _blob;
// this is string data
	std::string _summaryTable, _mergedTPtable,
		_transmissionTable,_signalTable, _statsTable , _solarTable, _solarAngleTable,
		_ephemTable, _timeTable, _lockdownTable, _altsTable,
		_driftTable,_solarDatTable,_miscTable,_msisTable;

};

extern Level1Data gLevel1Data;

#endif //_LEVEL1DATA_H_