Put 0.1 as minimal value of the gas gain (Raphaelle)
[u/mrichter/AliRoot.git] / T0 / AliT0RawReader.cxx
index f146111..3aefb2d 100644 (file)
 //                                                                          
 // T0 
 // Class for reading T0 RAW data in TOF data format
-// Alla.Maevskaya@cern.ch
 //
-//____________________________________________________________________
-
 #include "AliT0RawReader.h"
 #include "AliT0Parameters.h"
 #include "AliBitPacking.h"
 #include "TBits.h"
 
-//#include <Riostream.h>
+#include <Riostream.h>
+#include "TMath.h"
+#include "TH1F.h"
+#include "TArrayI.h"
 #include "AliLog.h"
  
 ClassImp(AliT0RawReader)
@@ -39,7 +39,10 @@ ClassImp(AliT0RawReader)
        fData(NULL),
        fPosition(0),
        fParam(NULL),
-       fIsOnline(isOnline)
+       fIsOnline(isOnline),
+       fBunchID(0),
+       fPrintout(kFALSE)
+
 {
   //
 // create an object to read T0raw digits
@@ -58,6 +61,15 @@ ClassImp(AliT0RawReader)
 {
   // 
 }
+/*
+AliT0RawReader::AliT0RawReader(const AliT0RawReader& o): TTask(o),
+     fRawReader(rawReader),
+       fData(NULL),
+       fPosition(0)
+{
+  //
+}
+*/
 
 
 Bool_t  AliT0RawReader::Next()
@@ -88,54 +100,57 @@ Bool_t  AliT0RawReader::Next()
 //  allData[48]  mean (T0) signal  
 // allData[49]   time difference (vertex)
 
   UInt_t word;
   Int_t time=0,  itdc=0, ichannel=0, uu; 
   Int_t numberOfWordsInTRM=0, iTRM=0;
   Int_t tdcTime, koef,hit=0;
   Int_t koefhits[110];
-  Int_t trmChainHeader =  0x00000000;
-  Int_t  trmChainTrailer =  0x10000000;
+  Int_t trm_chain_header =  0x00000000;
+  Int_t  trm_chain_trailer =  0x10000000;
   
-  Int_t  filler =  0x70000000;
+  UInt_t  filler =  0x70000000;
   Bool_t correct=kTRUE;
   Int_t header;
 
-   Int_t fNTRM = fParam->GetNumberOfTRMs();
-   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 = 0;
-    //    cout.setf( ios_base::hex, ios_base::basefield );
-    
-    //DRM header
+    }
+  }
+  do {
+    if (!fRawReader->ReadNextData(fData)) return kFALSE;
+  } while (fRawReader->GetDataSize() == 0);
+  
+  fPosition = 0;
+  //  cout.setf( ios_base::hex, ios_base::basefield );
+  if(fPrintout)
+    cout<<" CDH :: BC ID "<< (fRawReader->GetBCID())<<
+      " Event size"<<fRawReader->GetDataSize()<<
+      " orbit ID "<<fRawReader->GetOrbitID()<< 
+      " event index "<<fRawReader->GetEventIndex()<<
+      " event type " <<fRawReader->GetType()<<endl;
+  //DRM header
     for (Int_t i=0; i<6; i++) {
       word = GetNextWord();
-      if(i==0 &&  AliBitPacking::UnpackWord(word,4,20) ==7) fNTRM=0;
-      //      cout<<" DRM header "<<word<<endl;
+      if(fPrintout && i==0) cout<<" DRM header:: event words "<<AliBitPacking::UnpackWord(word,4, 20);
+      //  cout<<i<<" DRM header "<<word<<endl;
+      if (fPrintout && i==4 ) cout<<" L0BC ID "<< AliBitPacking::UnpackWord(word,4, 15)<<endl;
       header = AliBitPacking::UnpackWord(word,28,31);
-      if( header !=4 )
+      if( header !=4 ) 
        {
          AliWarning(Form(" !!!! wrong  DRM header  %x!!!!", word));
          fRawReader->AddFatalErrorLog(kWrongDRMHeader,Form("w=%x",word));
          break;
        }
     }
-    // cout<<" nTRMs "<<fNTRM<<endl;
-      for (Int_t ntrm=0; ntrm< fNTRM; ntrm++)
+    for (Int_t ntrm=0; ntrm< fNTRM; ntrm++)
       {
        //TRMheader  
        word = GetNextWord();
+       //      cout<<" TRM "<<word<<endl;
        header = AliBitPacking::UnpackWord(word,28,31);
-       Int_t event = AliBitPacking::UnpackWord(word,4,15);
-       //      cout<<" TRM header "<<word<<" event "<<event<<endl;
        if ( header != 4 )
          {
            AliWarning(Form(" !!!! wrong TRM header  %x!!!!", word));
@@ -143,35 +158,43 @@ Bool_t  AliT0RawReader::Next()
            break;
          }
        numberOfWordsInTRM=AliBitPacking::UnpackWord(word,4,16);
-       //cout<<" numberOfWordsInTRM "<<numberOfWordsInTRM<<endl;
+       if(fPrintout) {
+         cout<<" TRM header :: event words "<<numberOfWordsInTRM;
+         cout<<" ACQ bits "<<AliBitPacking::UnpackWord(word,17,18);
+         cout<<" L bit "<<AliBitPacking::UnpackWord(word,19,19)<<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 & trmChainHeader;
-           if(uu != trmChainHeader) 
+           // cout<<" chain header "<<word<<endl;
+           uu = word & trm_chain_header;
+           if(uu != trm_chain_header) 
              {
                AliWarning(Form(" !!!! wrong CHAIN  0  header %x!!!!", word));
                fRawReader->AddMajorErrorLog(kWrongChain0Header,Form("w=%x",word));
                break;
              }
+           fBunchID=AliBitPacking::UnpackWord(word,4,15);
+           if(fPrintout)
+             cout<<" chain "<< ichain<<" header:: BunchID  "<<fBunchID;
            word = GetNextWord();
-           // cout<<" next "<<word<<endl;
+           //   cout<<" next word "<<word<<endl;
            tdcTime =  AliBitPacking::UnpackWord(word,31,31);   
+           //      for (; tdcTime==1; tdcTime) 
            while(tdcTime==1)
              {
-               //                      cout<<" packed "<<word<<endl;
+               correct = kTRUE;
                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;
+                 //                              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);
@@ -180,26 +203,28 @@ Bool_t  AliT0RawReader::Next()
                if(correct){
                  hit=koefhits[koef];
                  if(hit>5) {  
-                   AliWarning(Form("Too many hits for %i channel  ! ",koef)); 
-                   break; 
+                   AliDebug(10,Form("Too many hits for %i channel - %i ! ",koef,koefhits[koef] )); 
+                   hit=5;
                  }
                  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&trmChainTrailer;
-           if(uu != trmChainTrailer )
+           //   cout<<" trailer :: "<<word<<endl;
+           uu = word&trm_chain_trailer;
+           if(uu != trm_chain_trailer )
              {
                AliWarning(Form(" !!!! wrong CHAIN 0 trailer %x !!!!", word));
                fRawReader->AddMajorErrorLog(kWrongChain0Trailer,Form("w=%x",word));
                break;
              }
+           if(fPrintout)
+             cout<<" trailer:: event counter "<< AliBitPacking::UnpackWord(word,16,27)<<endl;
          }
             
        word = GetNextWord(); //TRM trailer
@@ -211,9 +236,11 @@ Bool_t  AliT0RawReader::Next()
            fRawReader->AddMajorErrorLog(kWrongTRMTrailer,Form("w=%x",word));
            break;
          }
+       if(fPrintout)
+         cout<<"  TRM trailer :: event counter "<< AliBitPacking::UnpackWord(word,16,27)<<endl;
       } //TRM loop
     word = GetNextWord(); //
-    // cout<<" after TRM trailer "<<word<<endl;
+    //  cout<<" after TRM trailer "<<word<<endl;
     if (word == filler )  word = GetNextWord(); 
      header = AliBitPacking::UnpackWord(word,28,31);
      if( header !=5 )
@@ -221,8 +248,10 @@ Bool_t  AliT0RawReader::Next()
         AliWarning(Form(" !!!! wrong DRM GLOBAL trailer  %x!!!!", word));
         fRawReader->AddFatalErrorLog(kWrongDRMTrailer,Form("w=%x",word));
       }
-     // cout.setf( ios_base::dec, ios_base::basefield );
-         
+     if(fPrintout)
+       cout<<" DRM trailer ::event counter "<< AliBitPacking::UnpackWord(word,4,15)<<endl;
+     cout.setf( ios_base::dec, ios_base::basefield );
+    
      return kTRUE;
 }
 //_____________________________________________________________________________