HLT/exa/read.C

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