//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Waveform as recorded during the IROC ALICE Test
//      Scope: 
//
//
// Environment:
//      Software developed for the GEM-TPC detector
//
// Author List:
//      Felix Boehmer    TUM            (original author)
//
//
//-----------------------------------------------------------



#include "TpcWaveform.h"
#include <algorithm>
#include <map>
#include <cstdlib>

ClassImp(TpcWaveform)

TpcWaveform::TpcWaveform() 
{;}

TpcWaveform::TpcWaveform(long double timestamp, long int step, float tlevel, float tdelay,
			 float fstart, float fstop, float foffset, int nsamples)
  :
  fTimeStamp(timestamp),                
  fTimeStep(step),
  fTriggerLevel(tlevel),
  fTriggerDelay(tdelay),
  fFrameStart(fstart),
  fFrameStop(fstop),
  fFrameOffset(foffset),
  fNSamples(nsamples),
  fSamples(std::vector<float>(nsamples,0.)),
  fTimeStampStr("")
{;}

TpcWaveform::TpcWaveform(const char* timestamp, long int step, float tlevel, float tdelay,
			 float fstart, float fstop, float foffset, int nsamples)
  :
  fTimeStep(step),
  fTriggerLevel(tlevel),
  fTriggerDelay(tdelay),
  fFrameStart(fstart),
  fFrameStop(fstop),
  fFrameOffset(foffset),
  fNSamples(nsamples),
  fSamples(std::vector<float>(nsamples,0.))
{
  setTimeStamp(timestamp);
}


TpcWaveform::~TpcWaveform() 
{;}

void
TpcWaveform::setSamples(float* sArray) {
  for(unsigned int s=0; s<fNSamples; s++) 
    fSamples[s] = sArray[s];
}



float
TpcWaveform::getMaximum() const {
  std::vector<float> temp(fSamples);
  std::sort(temp.begin(),temp.end());
  return temp.back();    
}

float
TpcWaveform::getMinimum() const {
  std::vector<float> temp(fSamples);
  std::sort(temp.begin(),temp.end());
  return temp.front();    
}


int 
TpcWaveform::getMaxSample() const {
  std::map<float,unsigned int> temp;
  for(unsigned int s=0; s<fNSamples; s++) 
    temp[(fSamples[s])] = s;
  return temp.end()->second;
}

int 
TpcWaveform::getMinSample() const {
  std::map<float,unsigned int> temp;
  for(unsigned int s=0; s<fNSamples; s++) 
    temp[(fSamples[s])] = s;
  return temp.begin()->second;
}
  
void
TpcWaveform::setNSamples(int s) {
  fNSamples = s;
  fSamples.resize(s,0);
}

void
TpcWaveform::setTimeStamp(const char* ts) {
  fTimeStamp = strtold(ts,NULL);
  fTimeStampStr.assign(ts);
}