/////////////////////////////////////////////////
//    TSTELEvent IMPLEMENTATION 
//    Apr 2011
//    ver. 1.0 event unpacker
//    Taken from N. Kurs 05/06/2008 
/////////////////////////////////////////////////
#include "TSTELEvent.h"
#include "Riostream.h"
#include "TROOT.h"
#include "TH1.h"
#include "TH2.h"
#include "TList.h"
#include <stdio.h>

TSTELEvent *gstelevent = 0;

//Constructor
TSTELEvent::TSTELEvent()
{ 
  ifstream cfile;
  if (gstelevent==0) gstelevent = this;
/*03062011
  fX1raw  = MakeTH1('I',"X1raw","X1raw",2048,0,2048);
  fX2raw  = MakeTH1('I',"X2raw","X2raw",2048,0,2048);
  fY1raw  = MakeTH1('I',"Y1raw","Y1raw",2048,0,2048);
  fY2raw  = MakeTH1('I',"Y2raw","Y2raw",2048,0,2048);
  fTS1raw = MakeTH1('I',"TS1raw","TS1raw",2048,0,2048);
  fTS2raw = MakeTH1('I',"TS2raw","TS2raw",2048,0,2048);
  fQS1raw = MakeTH1('I',"QS1raw","QS1raw",2048,0,2048);
  fQS2raw = MakeTH1('I',"QS2raw","QS2raw",2048,0,2048);
  fXY     = MakeTH2('I',"XY", "XY",200,-20,20,200,-20,20);
  fTS1    = MakeTH1('I',"TS1","TS1",40,-1,1);
  fTS2    = MakeTH1('I',"TS2","TS2",40,-1,1);
03062011*/
//  cout << "open stelz file" << endl;
  cfile.open("SteCalib.dat");
  cfile >> CalFacX;
  cfile >> CalOffX;
  cfile >> CalFacY;
  cfile >> CalOffY;
  cfile >> CalFacTS1;
  cfile >> CalOffTS1;
  cfile >> CalFacTS2;
  cfile >> CalOffTS2;
  cfile.close();
  //  cout << "close stelz file" << endl;
}
//Destructor 
TSTELEvent::~TSTELEvent()
{
// cout<<"TSTELEvent destructor...done"<<endl;
 gstelevent = 0;
}


int TSTELEvent::Calibrate()  
{
  X = CalFacX*(X1raw-X2raw) + CalOffX;
  Y = CalFacY*(Y1raw-Y2raw) + CalOffY;
  TS1 = CalFacTS1*TS1raw + CalOffTS1;
  TS2 = CalFacTS2*TS2raw + CalOffTS2;

  return 0;
}
/*03062011
int TSTELEvent::FillHistos()  
{

 fX1raw->Fill(X1raw);
 fX2raw->Fill(X2raw);
 fY1raw->Fill(Y1raw);
 fY2raw->Fill(Y2raw);
 fTS1raw->Fill(TS1raw);
 fTS2raw->Fill(TS2raw);
 fQS1raw->Fill(QS1raw);
 fQS2raw->Fill(QS2raw);
 fXY->Fill(X,Y);
 fTS1->Fill(TS1);
 fTS2->Fill(TS2);
  
 return 0;
}03062011*/