Add cluster error and shape parameterization class + routine to fit parameters

[u/mrichter/AliRoot.git] / TPC / AliTPCDDLRawData.C

#if !defined(__CINT__) 

#include "AliTPCDDLRawData.h"
#include "AliTPCCompression.h"
#endif


void AliTPCDDLRawData(Int_t eventNumber=0, Int_t LDCsNumber=12){
  AliTPCDDLRawData *util=new AliTPCDDLRawData();
  AliTPCCompression *u=new AliTPCCompression();
  TStopwatch timer;
  static const Int_t NumTable=5;

  util->SetVerbose(1);
  u->SetVerbose(1);  
  
  //The Altro File "AltroFormatDDL.dat" is built from "AliTPCDDL.dat"
  //util->RawDataAltro();
  
  /*
  //The file "AltroFormatDDL.dat" is converted in a txt file "AltroFormatDDL.txt"
  //that is used for debugging
  u->ReadAltroFormat("AltroFormatDDL.txt","AltroFormatDDL.dat");
  */

  /*
  //TABLES CREATION 
  //Tables are created and stored in as sequence of binary files
  u->CreateTables("AltroFormatDDL.dat",NumTable);
  */


  while (eventNumber<=0){
    cout<<"Insert the event number:";
    cin>>eventNumber;
    cout<<endl;
  }

  //SLICE CREATION
  //Slices are built here
  timer.Start();
  util->RawData(LDCsNumber,eventNumber);
  timer.Stop();
  timer.Print();

  /*
  //SLICE CHECKING
  //An Altro File is created from the slides
  cout<<"slice control"<<endl;
  util->RawDataAltroDecode(LDCsNumber,eventNumber,0);
  ///The Altro file AltroDDLRecomposed.dat is converted in a txt file AltroDDLRecomposed.txt
  //This file must be equal to the ones created above.
  u->ReadAltroFormat("AltroDDLRecomposed.txt","AltroDDLRecomposed.dat");
  */
  

  //SLICE COMPRESSION
  cout<<"Slice Compression"<<endl;
  //Slices are compressed here using the tables created above or an optimized set of tables 
  //(Tables file for Huffman coding are required)
  timer.Start();
  util->RawDataCompDecompress(LDCsNumber,eventNumber,0);
  timer.Stop();
  timer.Print();

  /*  
  //SLICE DECOMPRESSION
  timer.Start();
  util->RawDataCompDecompress(LDCsNumber,eventNumber,1);
  timer.Stop();
  timer.Print();
  */
  
  /*
  //SLICE DECOMPRESSED CHECKING  
  //A new Altro file is created from the decompressed slides
  util->RawDataAltroDecode(LDCsNumber,eventNumber,1);
  //Convertion of the Altro file AltroDDLRecomposedDec.dat in a txt file AltroDDLRecomposedDec.txt
  //Useful for debugging
  u->ReadAltroFormat("AltroDDLRecomposedDec.txt","AltroDDLRecomposedDec.dat");
  */
  delete util;
  delete u;
  return;
}