#include "AliSTARTRawReader.h"
#include "AliSTARTRawData.h"
#include "AliSTARTdigit.h"
#include "AliBitPacking.h"

#include <Riostream.h>
#include "TMath.h"
#include "TH1F.h"
#include "TArrayI.h"
#include "AliLog.h"
 
ClassImp(AliSTARTRawReader)
  
  AliSTARTRawReader::AliSTARTRawReader (AliRawReader *rawReader, TTree* tree)
    :  TTask("STARTRawReader","read raw START data"),
       fTree(tree),
       fRawReader(rawReader)
{
  //
// create an object to read STARTraw digits
  AliDebug(1,"Start ");

 
}
 AliSTARTRawReader::~AliSTARTRawReader ()
{
  // 
}

Bool_t  AliSTARTRawReader::Next()
{
// read the next raw digit
// returns kFALSE if there is no digit left
  AliBitPacking *pack ;

  UInt_t word, unpackword; 
  Int_t time, adc, pmt;
  TArrayI *timeTDC1 = new TArrayI(24);
  TArrayI * chargeTDC1 = new TArrayI(24);
  TArrayI *timeTDC2 = new TArrayI(24);
  TArrayI *chargeTDC2 = new TArrayI(24);

  fRawReader->Reset();
  fRawReader->Select(13,0,1);
 
  
  if (!fRawReader->ReadHeader()){
    Error("ReadSTARTRaw","Couldn't read header");
    return kFALSE;
  }
  AliSTARTdigit *fDigits = new AliSTARTdigit(); 
  fTree->Branch("START","AliSTARTdigit",&fDigits,400,1);
  if (!fRawReader->ReadNextInt(word)) return kFALSE;
   
  for (Int_t i=0; i<24; i++) //time LED
    {
      word=0;
      unpackword=0;
      fRawReader->ReadNextInt(word);
       unpackword=pack->UnpackWord(word,0,8);
      pmt=unpackword; 
     
      word=0;
      unpackword=0;
      fRawReader->ReadNextInt(word);
      unpackword=pack->UnpackWord(word,8,31);
      time=unpackword;
      timeTDC1->AddAt(time,i);
      word=0;
      unpackword=0;
    }
  cout<<endl;
  for (Int_t i=0; i<24; i++) //time CFD
    {
      word=0;
      unpackword=0;
      fRawReader->ReadNextInt(word);
      unpackword=pack->UnpackWord(word,0,8);
      pmt=unpackword;
      word=0;
      unpackword=0;
      fRawReader->ReadNextInt(word);
      unpackword=pack->UnpackWord(word,8,31);
      time=unpackword;
      timeTDC2->AddAt(time,pmt-24);
      word=0;
      unpackword=0;
    } 
  cout<<endl;

  for (Int_t i=0; i<24; i++)
    {
      //  QTC
      word=0;
      unpackword=0;
      fRawReader->ReadNextInt(word);
      unpackword=pack->UnpackWord(word,0,8);
      pmt=unpackword;
      word=0;
      unpackword=0;
      fRawReader->ReadNextInt(word);
      unpackword= pack->UnpackWord(word,8,31);
      adc=unpackword; //
      chargeTDC1->AddAt(adc,pmt-48);

      word=0;
      unpackword=0;
    }
  cout<<endl;
  
   for (Int_t i=0; i<24; i++)   //QTC amplified
     {
      word=0;
      unpackword=0;
    
      fRawReader->ReadNextInt(word);
      unpackword=pack->UnpackWord(word,0,8);
      pmt=unpackword;
      word=0;
      unpackword=0;
      fRawReader->ReadNextInt(word);
      unpackword= pack->UnpackWord(word,8,31);
      adc=unpackword;  
      chargeTDC2->AddAt(adc,pmt-72);
     }
   cout<<endl;
   fDigits->SetTime(*timeTDC2);
   fDigits->SetADC(*chargeTDC1);
 
   fDigits->SetTimeAmp(*timeTDC1);
   fDigits->SetADCAmp(*chargeTDC2);

   word=0;
   unpackword=0;
    
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,0,8);
   pmt=unpackword;

   word=0;
   unpackword=0;
   
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,8,31);
   time=unpackword;
   fDigits->SetMeanTime(time);   
    
   // Best time right &left  
   word=0;
   unpackword=0;
   
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,0,8);
   pmt=unpackword;

   word=0;
   unpackword=0;

   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,8,31);
   time=unpackword;
   fDigits->SetTimeBestRight(time);   
 

   // best time left 
   word=0;
   unpackword=0;
     
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,0,8);
   pmt=unpackword;

   word=0;
   unpackword=0;
   
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,8,31);
   time=unpackword;
   fDigits->SetTimeBestLeft(time);   
   

   // best time differece  
    word=0;
   unpackword=0;
   
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,0,8);
   pmt=unpackword;

   word=0;
   unpackword=0;
   
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,8,31);
   time=unpackword;
   fDigits->SetDiffTime(time);   
 
 //  multiplicity 
   word=0;
   unpackword=0;
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,0,8);
   pmt=unpackword;
   word=0;
   unpackword=0;
   fRawReader->ReadNextInt(word);
   unpackword=pack->UnpackWord(word,8,31);
   time=unpackword;
   fDigits->SetSumMult(time);   
      
    fTree->Fill();

    return kTRUE;
}