HLT/exa/read.C

   1 // $Id$
   2
   3 void read(int min=0,int max=35)
   4 {
   5
   6   for(int slice=0; slice<35; slice++)
   7     {
   8       char fname[256];
   9       //sprintf(fname,"/prog/alice/data/Rawdata/PileUp/digits_%d_0.raw",slice);
  10       sprintf(fname,"/prog/alice/data/Rawdata/1_patch/pp/test_pileup/digits_%d_0.raw",slice);
  11       //sprintf(fname,"digits_%d_0.raw",slice);
  12       file = new AliL3FileHandler();
  13       if(!file->SetBinaryInput(fname))
  14         {
  15           cerr<<"Error opening file "<<fname<<endl;
  16           return;
  17         }
  18       int row[2]={0,175};
  19       file->Init(slice,0,row);
  20       UInt_t size;
  21       char name[256];
  22       AliL3DigitRowData *data = file->CompBinary2Memory(size);
  23
  24       for(Int_t r=0; r<175; r++)
  25         {
  26           UInt_t padrow=data->fRow;
  27           AliL3DigitData *dPt = (AliL3DigitData*)data->fDigitData;
  28           cout<<"padrow "<<padrow<<" ndigits "<<data->fNDigit<<endl;
  29
  30           for(Int_t d=0; d<data->fNDigit; d++)
  31             {
  32               if(d>0 && dPt[d]->fPad == dPt[d-1]->fPad && dPt[d]->fTime == dPt[d-1]->fTime)
  33                 cout<<"Slice "<<slice<<" padrow "<<padrow<<" pad "<<(int)dPt[d]->fPad<<" time "
  34                     <<(int)dPt[d]->fTime<<" charge "<<(int)dPt[d]->fCharge<<endl;
  35             }
  36
  37           file->UpdateRowPointer(data);
  38         }
  39
  40       file->CloseBinaryInput();
  41       delete file;
  42     }
  43 }