Setting run number in CDB manager added. GRP now read from default cdb storage instea...
[u/mrichter/AliRoot.git] / FIT / AliFITRawReader.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id:  */
17
18 //____________________________________________________________________
19 //                                                                          
20 // FIT 
21 // Class for reading FIT RAW data in TOF data format
22 //
23 #include "AliFITRawReader.h"
24 #include "AliBitPacking.h"
25 #include "TBits.h"
26
27 #include <Riostream.h>
28 #include "TMath.h"
29 #include "TH1F.h"
30 #include "AliLog.h"
31
32 using std::cout;
33 using std::endl;
34 using std::ios_base;
35 ClassImp(AliFITRawReader)
36   
37   AliFITRawReader::AliFITRawReader (AliRawReader *rawReader)
38     :  TTask("FITRawReader","read FIT raw  data"),
39        fRawReader(rawReader),
40        fData(NULL),
41        fPosition(0),
42        fBunchID(0),
43        fPrintout(kFALSE)
44     
45 {
46   //
47   // create an object to read FIT raw digits
48   
49   fRawReader->Reset();
50   fRawReader->Select("FIT"); 
51   for ( Int_t k=0; k<500; k++)   fAllData[k] = -1;
52  
53 }
54 //_____________________________________________________________________________
55 AliFITRawReader::~AliFITRawReader ()
56 {
57   // 
58 }
59
60 //_____________________________________________________________________________
61
62 Bool_t  AliFITRawReader::Next()
63 {
64   // read the next raw digit
65   // returns kFALSE if there is no digit left
66 // allData array collect data from all channels in one :
67 // allData[0] - allData[79]   CFD channels side C
68 // allData[80] - allData[159] CFD channels side A
69 // allData[160] -   allData[239]  QT0 channels side C
70 // allData[240] -   allData[319]  QT0 channels side A
71 // allData[320] -   allData[399]  QT1 channels side C
72 // allData[400] -   allData[479]  QT1 channels side A
73
74   UInt_t word;
75   Int_t time=0,  itdc=0, ichannel=0, uu; 
76   Int_t numberOfWordsInTRM=0, iTRM=0;
77   Int_t tdcTime, koef;
78   Int_t trm_chain_header =  0x00000000;
79   Int_t  trm_chain_trailer[2] = { 0x10000000, 0x20000000};
80   UInt_t  filler =  0x70000000;
81   Bool_t correct=kTRUE;
82   Int_t header;
83   
84   for ( Int_t k=0; k<500; k++)   fAllData[k] = -1;
85   
86   
87   do {
88     if (!fRawReader->ReadNextData(fData)) return kFALSE;
89   } while (fRawReader->GetDataSize() == 0);
90   
91   fPosition = 0;
92   cout.setf( ios_base::hex, ios_base::basefield );
93   if(fPrintout)
94     cout<<" CDH :: BC ID "<< (fRawReader->GetBCID())<<
95       " Event size"<<fRawReader->GetDataSize()<<
96       " orbit ID "<<fRawReader->GetOrbitID()<< 
97       " event index "<<fRawReader->GetEventIndex()<<
98       " event type " <<fRawReader->GetType()<<endl;
99   //DRM header
100   for (Int_t i=0; i<6; i++) {
101     word = GetNextWord();
102     if(fPrintout && i==0) cout<<" DRM header:: event words "<<AliBitPacking::UnpackWord(word,4, 20)<<endl;;
103     if (fPrintout && i==4 ) cout<<" L0BC ID "<< AliBitPacking::UnpackWord(word,4, 15)<<endl;
104     header = AliBitPacking::UnpackWord(word,28,31);
105     if( header !=4 ) 
106       {
107         AliWarning(Form(" !!!! wrong  DRM header  %x!!!!", word));
108         fRawReader->AddFatalErrorLog(kWrongDRMHeader,Form("w=%x",word));
109         break;
110       }
111   }
112   for (Int_t ntrm=0; ntrm<4; ntrm++)
113     {
114       //TRMheader  
115       word = GetNextWord();
116       header = AliBitPacking::UnpackWord(word,28,31);
117       if ( header != 4 )
118         {
119           AliWarning(Form(" !!!! wrong TRM header  %x!!!!", word));
120           fRawReader->AddMajorErrorLog(kWrongTRMHeader,Form("w=%x",word));
121           break;
122         }
123       numberOfWordsInTRM=AliBitPacking::UnpackWord(word,4,16);
124       if(fPrintout) {
125         cout<<" TRM header :: event words "<<numberOfWordsInTRM;
126         cout<<" ACQ bits "<<AliBitPacking::UnpackWord(word,17,18);
127         cout<<" L bit "<<AliBitPacking::UnpackWord(word,19,19)<<endl;
128       }
129       iTRM=AliBitPacking::UnpackWord(word,0,3);
130       for( Int_t ichain=0; ichain<2; ichain++)
131         {
132           //chain header
133           word = GetNextWord();
134           uu = word & trm_chain_header;
135           if(uu != trm_chain_header) 
136             {
137               AliWarning(Form(" !!!! wrong CHAIN  0  header %x!!!!", word));
138               fRawReader->AddMajorErrorLog(kWrongChain0Header,Form("w=%x",word));
139               break;
140             }
141           fBunchID=AliBitPacking::UnpackWord(word,4,15);
142           if(fPrintout)
143             cout<<" chain "<< ichain<<" header:: BunchID  "<<fBunchID<<endl;
144           word = GetNextWord();
145           tdcTime =  AliBitPacking::UnpackWord(word,31,31);   
146           while(tdcTime==1)
147             {
148               correct = kTRUE;
149               itdc=AliBitPacking::UnpackWord(word,24,27);
150               ichannel=AliBitPacking::UnpackWord(word,21,23);
151               time=AliBitPacking::UnpackWord(word,0,20);
152               
153               koef = GetChannel(iTRM,itdc,ichain,ichannel);
154               if (koef != 0 && fPrintout) 
155                 cout<<"RawReader>> "<<"koef "<<koef<<" trm "<<iTRM<<" tdc "<<itdc<<" chain "<<ichain<<              " channel "<<ichannel<<" time "<<time<<endl;
156               if (koef ==-1 ){
157                 AliWarning(Form("Incorrect lookup table ! "));
158                 fRawReader->AddMajorErrorLog(kIncorrectLUT);
159                 correct=kFALSE;
160               }
161               if(correct)   fAllData[koef]=time; 
162               word = GetNextWord();
163                 tdcTime =  AliBitPacking::UnpackWord(word,31,31);   
164             }
165             
166           uu = word&trm_chain_trailer[ichain];
167           if(uu != trm_chain_trailer[ichain] )
168             {
169               AliWarning(Form(" !!!! wrong CHAIN %i trailer %x !!!!", ichain, word));
170               fRawReader->AddMajorErrorLog(kWrongChain0Trailer,Form("w=%x",word));
171               break;
172             }
173           if(fPrintout)
174             cout<<" trailer:: event counter "<< AliBitPacking::UnpackWord(word,16,27)<<endl;
175         }
176       
177         word = GetNextWord(); //TRM trailer
178         header = AliBitPacking::UnpackWord(word,28,31);
179         if( header !=5 )
180           {
181             AliWarning(Form(" !!!! wrong TRM GLOBAL trailer  %x!!!!", word));
182             fRawReader->AddMajorErrorLog(kWrongTRMTrailer,Form("w=%x",word));
183             break;
184           }
185         if(fPrintout)
186           cout<<"  TRM trailer :: event counter "<< AliBitPacking::UnpackWord(word,16,27)<<endl;
187     } //TRM loop
188   word = GetNextWord(); //
189   if (word == filler )  word = GetNextWord(); 
190   header = AliBitPacking::UnpackWord(word,28,31);
191   if( header !=5 )
192     {
193       AliWarning(Form(" !!!! wrong DRM GLOBAL trailer  %x!!!!", word));
194       // fRawReader->AddFatalErrorLog(kWrongDRMTrailer,Form("w=%x",word));
195     }
196   if(fPrintout)
197     cout<<" DRM trailer ::event counter "<< AliBitPacking::UnpackWord(word,4,15)<<endl;
198   cout.setf( ios_base::dec, ios_base::basefield );
199   
200   return kTRUE;
201 }
202 //_____________________________________________________________________________
203 Int_t AliFITRawReader::GetPosition()
204 {
205   // Sets the position in the
206   // input stream
207   if (((fRawReader->GetDataSize() * 8) % 32) != 0)
208     AliFatal(Form("Incorrect raw data size ! %d words are found !",fRawReader->GetDataSize()));
209   return (fRawReader->GetDataSize() * 8) / 32;
210 }
211 //_____________________________________________________________________________
212 UInt_t AliFITRawReader::GetNextWord()
213 {
214   // Read the next 32 bit word in backward direction
215   // The input stream access is given by fData and fPosition
216   
217   
218   //   fPosition--;
219   Int_t iBit = fPosition * 32;
220   Int_t iByte = iBit / 8;
221   
222   UInt_t word = 0;
223   word  = fData[iByte+3]<<24;
224   word |= fData[iByte+2]<<16;
225   word |= fData[iByte+1]<<8;
226   word |= fData[iByte];
227   fPosition++;
228   
229   return word;   
230 }
231
232 //_____________________________________________________________________________
233 UInt_t AliFITRawReader::GetChannel(Int_t iTRM, Int_t iTDC, Int_t iChain, Int_t ichannel)
234 {
235   return   iTRM*120 + iChain*60 + iTDC*4 +ichannel;
236 }