Added data members to the track segment data structure to make

[u/mrichter/AliRoot.git] / HLT / TPCLib / Ref / AliHLTTPCMemHandler.h

// @(#) $Id$

#ifndef ALIHLTTPC_MEMHANDLER_H
#define ALIHLTTPC_MEMHANDLER_H

#include "AliHLTTPCRootTypes.h"
#include "AliHLTTPCDigitData.h"

class AliHLTTPCSpacePointData;
class AliHLTTPCDigitRowData;
class AliHLTTPCTrackSegmentData;
class AliHLTTPCTrackArray;
class AliHLTTPCRandomPointData;

class AliHLTTPCMemHandler{
 private:
  
  Byte_t *fPt;//!
  UInt_t fSize;

  AliHLTTPCRandomDigitData **fDPt;//!
  AliHLTTPCRandomDigitData *fDigits;//!
  Bool_t IsRandom;
  Int_t fNRandom;
  Int_t fNGenerate;
  Int_t fNUsed;
  Int_t fNDigits;

  void Write(UInt_t *comp, UInt_t & index, UInt_t & subindex, UShort_t value);
  UShort_t Read(UInt_t *comp, UInt_t & index, UInt_t & subindex);
  UShort_t Test(UInt_t *comp, UInt_t index, UInt_t subindex); 
  
  void DigitizePoint(Int_t row,Int_t pad, Int_t time,Int_t charge);
  void QSort(AliHLTTPCRandomDigitData **a, Int_t first, Int_t last);
  Int_t ComparePoints(UInt_t row,UShort_t pad,UShort_t time);
  Int_t CompareDigits(AliHLTTPCRandomDigitData *a,AliHLTTPCRandomDigitData *b);
  void AddData(AliHLTTPCDigitData *data,UInt_t & ndata,
                      UInt_t row,UShort_t pad,UShort_t time,UShort_t charge);
  void AddRandom(AliHLTTPCDigitData *data,UInt_t & ndata);
  void MergeDataRandom(AliHLTTPCDigitData *data,UInt_t & ndata,
                      UInt_t row,UShort_t pad,UShort_t time,UShort_t charge);
  void AddDataRandom(AliHLTTPCDigitData *data,UInt_t & ndata,
                      UInt_t row,UShort_t pad,UShort_t time,UShort_t charge);

 protected:
  Int_t fRowMin;
  Int_t fRowMax;
  Int_t fSlice;
  Int_t fPatch;

  Int_t fEtaMinTimeBin[159]; //for ROI in eta only
  Int_t fEtaMaxTimeBin[159];
  
  FILE *fInBinary;//!
  FILE *fOutBinary;//!
  
 public:
  AliHLTTPCMemHandler();
  virtual ~AliHLTTPCMemHandler();
  
  void Reset(){CloseBinaryInput();CloseBinaryOutput();Free();}  
  void Init(Int_t s,Int_t p, Int_t *r=0);

  Bool_t SetBinaryInput(char *name);
  Bool_t SetBinaryInput(FILE *file);
  void CloseBinaryInput();
  
  Bool_t SetBinaryOutput(char *name);
  Bool_t SetBinaryOutput(FILE *file);
  void CloseBinaryOutput();

  //Random cluster
  void SetRandomCluster(Int_t maxnumber);
  void SetRandomSeed(UInt_t seed){srand(seed);}
  void SetRandomSeed();

  void ResetRandom(){fNDigits = 0; fNUsed = 0;}
  void Generate(Int_t row);
  void SetNGenerate(Int_t number){(number>fNRandom)?fNGenerate=fNRandom:fNGenerate = number;}

  void SetROI(Float_t *eta,Int_t *slice);
  void ResetROI();

  //Digit IO
  Bool_t Memory2Binary(UInt_t nrow,AliHLTTPCDigitRowData *data);
  Bool_t Binary2Memory(UInt_t & nrow,AliHLTTPCDigitRowData *data);

  Int_t Memory2CompMemory(UInt_t nrow,AliHLTTPCDigitRowData *data,UInt_t *comp);
  Int_t CompMemory2Memory(UInt_t nrow,AliHLTTPCDigitRowData *data,UInt_t *comp);
  Bool_t CompMemory2CompBinary(UInt_t nrow,UInt_t *comp, UInt_t size=0);
  Bool_t CompBinary2CompMemory(UInt_t & nrow,UInt_t *comp);

  virtual AliHLTTPCDigitRowData *CompBinary2Memory(UInt_t & nrow);
  virtual Bool_t Memory2CompBinary(UInt_t nrow,AliHLTTPCDigitRowData *data);
  
  UInt_t GetNRow(UInt_t *comp,UInt_t size);

  //Point IO
  Bool_t Memory2Binary(UInt_t npoint,AliHLTTPCSpacePointData *data);
  Bool_t Binary2Memory(UInt_t & npoint,AliHLTTPCSpacePointData *data);
  Bool_t Transform(UInt_t npoint,AliHLTTPCSpacePointData *data,Int_t slice);
  static void UpdateRowPointer(AliHLTTPCDigitRowData *&tempPt);
  
  //Track IO
  Bool_t Memory2Binary(UInt_t ntrack,AliHLTTPCTrackSegmentData *data);
  Bool_t Binary2Memory(UInt_t & ntrack,AliHLTTPCTrackSegmentData *data);
  Bool_t TrackArray2Binary(AliHLTTPCTrackArray *array);
  Bool_t Binary2TrackArray(AliHLTTPCTrackArray *array);
  Bool_t TrackArray2Memory(UInt_t & ntrack,AliHLTTPCTrackSegmentData *data,AliHLTTPCTrackArray *array);
  Bool_t Memory2TrackArray(UInt_t ntrack,AliHLTTPCTrackSegmentData *data,AliHLTTPCTrackArray *array);
  Bool_t Memory2TrackArray(UInt_t ntrack,AliHLTTPCTrackSegmentData *data,AliHLTTPCTrackArray *array,Int_t slice);
    
  //Memory Allocation
  UInt_t GetAllocatedSize(){return fSize;}  
  UInt_t GetFileSize();
  UInt_t GetMemorySize(UInt_t nrow,UInt_t *comp);
  UInt_t GetCompMemorySize(UInt_t nrow,AliHLTTPCDigitRowData *data);
  UInt_t GetRandomSize();

  Byte_t *Allocate(UInt_t size);
  Byte_t *Allocate();  // allocate size of Binary Input File
  Byte_t *Allocate(AliHLTTPCTrackArray *array);
  Byte_t *GetDataPointer(UInt_t &size) {size = fSize; return fPt;}
  FILE *GetFilePointer() {return fInBinary;}
  void   Free();
  
  //Getters:
  Int_t GetRowMin(){return fRowMin;}
  Int_t GetRowMax(){return fRowMax;}
  Int_t GetSlice(){return fSlice;}
  Int_t GetPatch(){return fPatch;}
  
  //virtual functions:
  virtual void FreeDigitsTree() {return;}
  virtual Bool_t SetAliInput(char *name){return 0;}
  virtual void CloseAliInput(){return;} 
  virtual Bool_t IsDigit(Int_t i=0){return 0;}
  virtual Bool_t SetMCOutput(char *name){return 0;}
  virtual Bool_t SetMCOutput(FILE *file){return 0;}
  virtual void CloseMCOutput(){return;}
  virtual Bool_t AliDigits2Binary(Int_t event=0,Bool_t altro=kFALSE){return 0;}
  virtual Bool_t AliDigits2CompBinary(Int_t event=0,Bool_t altro=kFALSE){return 0;}  
  virtual AliHLTTPCDigitRowData *AliDigits2Memory(UInt_t & nrow,Int_t event=0){return 0;}
  virtual AliHLTTPCDigitRowData *AliAltroDigits2Memory(UInt_t & nrow,Int_t event=0,Bool_t eventmerge=kFALSE){return 0;}
  virtual void AliDigits2RootFile(AliHLTTPCDigitRowData *rowPt,Char_t *new_digitsfile){return;}
  virtual Bool_t AliPoints2Binary(Int_t eventn=0){return 0;}
  virtual AliHLTTPCSpacePointData *AliPoints2Memory(UInt_t & npoint,Int_t eventn=0){return 0;}

  //AliHLTTPCRawDataFileHandler
  virtual Bool_t SetRawInput(Char_t *name){return 0;}
  virtual Bool_t SetRawInput(std::ifstream *file){return 0;}
  virtual void CloseRawInput(){} 
  virtual Int_t ReadRawInput(){return 0;}
  virtual Short_t** GetRawData(Int_t &channels, Int_t & timebins){return 0;}

  virtual Bool_t SetRawOutput(Char_t *name){return 0;}
  virtual Bool_t SetRawOutput(std::ofstream *file){return 0;}
  virtual void CloseRawOutput(){} 
  virtual Bool_t SaveRawOutput(){return 0;}

  virtual Bool_t SetMappingFile(Char_t *name){return 0;}
  virtual Bool_t SetMappingFile(FILE *file){return 0;}
  virtual void CloseMappingFile(){} 
  virtual Int_t ReadMappingFile(){return 0;}
  
  virtual Bool_t SetRawPedestalsInput(Char_t *name){return 0;}
  virtual Bool_t SetRawPedestalsInput(std::ifstream *file){return 0;}
  virtual void CloseRawPedestalsInput(){} 
  virtual Int_t ReadRawPedestalsInput(){return 0;}

  virtual AliHLTTPCDigitRowData* RawData2Memory(UInt_t &nrow,Int_t event=-1){return 0;}
  virtual Bool_t RawData2CompMemory(Int_t event=-1){return 0;}

  //AliHLTTPCDDLDataFileHandler
#ifdef use_newio
  virtual Bool_t SetReaderInput(Char_t *name,Int_t event=0){return 0;}
#else
  virtual Bool_t SetReaderInput(Char_t *name,Bool_t add=kTRUE){return 0;}
#endif
  virtual void CloseReaderInput(){};

  virtual AliHLTTPCDigitRowData* DDLData2Memory(UInt_t &nrow,Int_t event=-1){return 0;}
  virtual Bool_t DDLData2CompBinary(Int_t event=-1){return 0;}

  ClassDef(AliHLTTPCMemHandler,1) // Memory handler class
};

inline Int_t  AliHLTTPCMemHandler::ComparePoints(UInt_t row,UShort_t pad,UShort_t time){
  if(fNUsed>=fNDigits) return -2;

  if(pad==fDPt[fNUsed]->fPad&&time==fDPt[fNUsed]->fTime) return 0;

  if(pad<fDPt[fNUsed]->fPad) return -1;
  if(pad==fDPt[fNUsed]->fPad&&time<fDPt[fNUsed]->fTime)  return -1;

  return 1;
}

inline Int_t AliHLTTPCMemHandler::CompareDigits(AliHLTTPCRandomDigitData *a,AliHLTTPCRandomDigitData *b){
  if(a->fPad==b->fPad && a->fTime == b->fTime) return 0;

  if(a->fPad<b->fPad) return -1;
  if(a->fPad==b->fPad && a->fTime<b->fTime) return -1;
  
  return 1;
}

#endif