]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - T0/AliT0RawReader.cxx
fix bug in OnlineLookup
[u/mrichter/AliRoot.git] / T0 / AliT0RawReader.cxx
index dd1e72a959307c1c8c1a24c94a7b1e90dab72106..dbdff015dda6df3d63faaeac8cc6db577633aafd 100644 (file)
  
 ClassImp(AliT0RawReader)
   
-  AliT0RawReader::AliT0RawReader (AliRawReader *rawReader)
+  AliT0RawReader::AliT0RawReader (AliRawReader *rawReader, Bool_t isOnline)
     :  TTask("T0RawReader","read raw T0 data"),
        fRawReader(rawReader),
        fData(NULL),
-       fPosition(0)
+       fPosition(0),
+       fParam(NULL),
+       fIsOnline(isOnline)
 {
   //
 // create an object to read T0raw digits
@@ -23,13 +25,27 @@ ClassImp(AliT0RawReader)
  
   fRawReader->Reset();
   fRawReader->Select("T0");
+  fParam = AliT0Parameters::Instance();  
+  if (fIsOnline)
+    fParam->InitIfOnline();
+  else 
+    fParam->Init();
  
 }
  AliT0RawReader::~AliT0RawReader ()
 {
   // 
 }
+/*
+AliT0RawReader::AliT0RawReader(const AliT0RawReader& o): TTask(o),
+     fRawReader(rawReader),
+       fData(NULL),
+       fPosition(0)
+{
+  //
+}
+*/
+
 
 Bool_t  AliT0RawReader::Next()
 {
@@ -59,94 +75,139 @@ Bool_t  AliT0RawReader::Next()
 //  allData[48]  mean (T0) signal  
 // allData[49]   time difference (vertex)
 
-
- //  if (fDigits == 0x0) fDigits = new AliT0digit(); 
-  // fTree->Branch("T0","AliT0digit",&fDigits,405,1);
  
   UInt_t word;
-  Int_t time=0,  itdc=0, ichannel=0; 
+  Int_t time=0,  itdc=0, ichannel=0, uu
   Int_t numberOfWordsInTRM=0, iTRM=0;
-  Int_t tdcTime, koef,hit, meanTime, timeDiff ;
-
-
+  Int_t tdcTime, koef,hit=0;
+  Int_t koefhits[110];
+    Int_t  TRM_CHAIN_0_HEADER =  0x00000000;
+  Int_t  TRM_CHAIN_0_TRAILER =  0x10000000;
+  
+  Int_t  FILLER =  0x70000000;
+  Bool_t correct=kTRUE;
+  Int_t header;
 
-  AliT0Parameters* param = AliT0Parameters::Instance();   //-->Zhenya
-  param->Init();
- for ( Int_t k=0; k<110; k++) {
+   Int_t fNTRM = fParam->GetNumberOfTRMs();
+   for ( Int_t k=0; k<110; k++) {
+    koefhits[k]=0;
     for ( Int_t jj=0; jj<5; jj++) {
       fAllData[k][jj]=0;
-    }
-  }
-    do {
-    if (!fRawReader->ReadNextData(fData)) return kFALSE;
-  } while (fRawReader->GetDataSize() == 0);
-  
-  //  fPosition = GetPosition();
-  fPosition = 0;
-
-  //DRM header
-  for (Int_t i=0; i<4; i++) {
-    word = GetNextWord();
-  }
-  //TRMheader  
-   word = GetNextWord();
-   numberOfWordsInTRM=AliBitPacking::UnpackWord(word,4,16);
-   iTRM=AliBitPacking::UnpackWord(word,0,3);
-
-   //chain header
-   Int_t ichain=0;
-   word = GetNextWord();
-  
-   for (Int_t i=0; i<numberOfWordsInTRM; i++) {
-     word = GetNextWord();
-     tdcTime =  AliBitPacking::UnpackWord(word,31,31);   
-
-     if ( tdcTime == 1)
-       {
-        itdc=AliBitPacking::UnpackWord(word,24,27);
-        ichannel=AliBitPacking::UnpackWord(word,21,23);
-        time=AliBitPacking::UnpackWord(word,0,20);
-        //  koef = itdc*4 + ichannel/2;
-        koef = param->GetChannel(iTRM,itdc,ichain,ichannel);
-        //      cout<<" RawReader::Next ::"<<iTRM<<"  "<<itdc<<" "<<ichain<<" "<<ichannel<<" "<<  koef<<" "<<time<<endl;
-        if(fAllData[koef][0] == 0)  fAllData[koef][0]=time;  // yield only 1st particle
-         
-       }
+     }
    }
-   word = GetNextWord(); //chain trailer
-   word = GetNextWord(); //TRM trailer
-     
-  //TRMheader  
-   word = GetNextWord();
-   numberOfWordsInTRM=AliBitPacking::UnpackWord(word,4,16);
-   iTRM=AliBitPacking::UnpackWord(word,0,3);
-   
-   //chain header
-   word = GetNextWord();
-   
-   for (Int_t iword=0; iword<numberOfWordsInTRM; iword++) {
-     word = GetNextWord();
-     tdcTime =  AliBitPacking::UnpackWord(word,31,31);   
-
-     if ( tdcTime == 1)
+    do {
+      if (!fRawReader->ReadNextData(fData)) return kFALSE;
+    } while (fRawReader->GetDataSize() == 0);
+    
+    fPosition = 0;
+     cout.setf( ios_base::hex, ios_base::basefield );
+    
+    //DRM header
+    for (Int_t i=0; i<6; i++) {
+      word = GetNextWord();
+      header = AliBitPacking::UnpackWord(word,28,31);
+      if( header !=4 )
+       {
+         AliWarning(Form(" !!!! wrong  DRM header  %x!!!!", word));
+         fRawReader->AddFatalErrorLog(kWrongDRMHeader,Form("w=%x",word));
+         break;
+       }
+    }
+    
+    for (Int_t ntrm=0; ntrm< fNTRM; ntrm++)
+      {
+       //TRMheader  
+       word = GetNextWord();
+       header = AliBitPacking::UnpackWord(word,28,31);
+//     cout<<" TRM header "<<word<<endl;
+       if ( header != 4 )
+         {
+           AliWarning(Form(" !!!! wrong TRM header  %x!!!!", word));
+           fRawReader->AddMajorErrorLog(kWrongTRMHeader,Form("w=%x",word));
+           break;
+         }
+       numberOfWordsInTRM=AliBitPacking::UnpackWord(word,4,16);
+//     cout<<" numberOfWordsInTRM "<<numberOfWordsInTRM<<endl;
+       iTRM=AliBitPacking::UnpackWord(word,0,3);
+       for( Int_t ichain=0; ichain<2; ichain++)
+         {
+           //chain header
+           word = GetNextWord();
+//         cout<<" chain header "<<word<<endl;
+           uu = word & TRM_CHAIN_0_HEADER;
+           if(uu != TRM_CHAIN_0_HEADER) 
+             {
+               AliWarning(Form(" !!!! wrong CHAIN  0  header %x!!!!", word));
+               fRawReader->AddMajorErrorLog(kWrongChain0Header,Form("w=%x",word));
+               break;
+             }
+           word = GetNextWord();
+//         cout<<" next "<<word<<endl;
+           tdcTime =  AliBitPacking::UnpackWord(word,31,31);   
+//         for (; tdcTime==1; tdcTime) 
+       while(tdcTime==1)
+             {
+//                     cout<<" packed "<<word<<endl;
+               itdc=AliBitPacking::UnpackWord(word,24,27);
+               ichannel=AliBitPacking::UnpackWord(word,21,23);
+               time=AliBitPacking::UnpackWord(word,0,20);
+               
+               koef = fParam->GetChannel(iTRM,itdc,ichain,ichannel);
+                       if (koef != 0 ) 
+//               cout<<"RawReader>> "<<"koef "<<koef<<" trm "<<iTRM<<
+//                 " tdc "<<itdc<<" chain "<<ichain<<
+//                 " channel "<<ichannel<<" time "<<time<<endl;
+               if (koef ==-1 ){
+                 AliWarning(Form("Incorrect lookup table ! "));
+                 fRawReader->AddMajorErrorLog(kIncorrectLUT);
+                 correct=kFALSE;
+               }
+               if(correct){
+                 hit=koefhits[koef];
+                 if(hit>5) {  
+                   AliWarning(Form("Too many hits for %i channel  ! ",koef)); 
+                   break; 
+                 }
+                 fAllData[koef][hit]=time; 
+                 koefhits[koef]++;
+               }
+                 word = GetNextWord();
+//               cout<<" next word in cycle "<<word<<endl;
+                 tdcTime =  AliBitPacking::UnpackWord(word,31,31);   
+               
+           }
+           
+//         cout<<" chain trailer "<<word<<endl;
+           uu = word&TRM_CHAIN_0_TRAILER;
+           if(uu != TRM_CHAIN_0_TRAILER )
+             {
+               AliWarning(Form(" !!!! wrong CHAIN 0 trailer %x !!!!", word));
+                 fRawReader->AddMajorErrorLog(kWrongChain0Trailer,Form("w=%x",word));
+               break;
+             }
+         }
+            
+       word = GetNextWord(); //TRM trailer
+//     cout<<" TRM trailer "<<word<<endl;
+       header = AliBitPacking::UnpackWord(word,28,31);
+       if( header !=5 )
+         {
+           AliWarning(Form(" !!!! wrong TRM GLOBAL trailer  %x!!!!", word));
+           fRawReader->AddMajorErrorLog(kWrongTRMTrailer,Form("w=%x",word));
+           break;
+         }
+      } //TRM loop
+    word = GetNextWord(); //
+    if (word == FILLER )  word = GetNextWord(); 
+     header = AliBitPacking::UnpackWord(word,28,31);
+     if( header !=5 )
        {
-        itdc=AliBitPacking::UnpackWord(word,24,27);
-        ichannel=AliBitPacking::UnpackWord(word,21,23);
-        time=AliBitPacking::UnpackWord(word,0,20);
-        //      koef = itdc*4 + ichannel/2;
-        koef = param->GetChannel(iTRM,itdc,ichain,ichannel);
-          if(fAllData[koef][0] == 0)    fAllData[koef][0]=time;
-                //      if(allData[koef+55] == 0) allData[koef+55]=time; // yield only 1st particle
-       }
-   }
-   meanTime = fAllData[49][0];  // T0 !!!!!!
-   timeDiff = fAllData[50][0];
-
-   word = GetNextWord();
-   word = GetNextWord();
-   return kTRUE;
+        AliWarning(Form(" !!!! wrong DRM GLOBAL trailer  %x!!!!", word));
+        fRawReader->AddFatalErrorLog(kWrongDRMTrailer,Form("w=%x",word));
+      }
+     cout.setf( ios_base::dec, ios_base::basefield );
+    
+     return kTRUE;
 }
 //_____________________________________________________________________________
 Int_t AliT0RawReader::GetPosition()