Added new library libFMDutil. This library contains utility classes that
[u/mrichter/AliRoot.git] / FMD / scripts / ShowRaw.C
1 //
2 // $Id$
3 //
4 // Script to read a raw data file, and dump it to std::cout 
5 //
6 #include <iomanip>
7
8 void 
9 ShowRaw(Int_t det=2,  bool verbose=false, Int_t event=0) 
10 {
11   TString file(Form("raw%d/FMD_%d.ddl", event, AliFMD::kBaseDDL + det - 1));
12
13   std::cout << "Reading raw data file " << file << std::endl;
14   
15   TH1* h = new TH1F("rawData", "Raw Data", 128, 0, 1024);
16   
17   
18   // This method creates a text file containing the same information
19   // stored in an Altro file. The information in the text file is
20   // organized pad by pad and and for each pad it consists in a
21   // sequence of bunches (Bunch length +2, Time bin of the last
22   // amplitude sample in the bunch, amplitude values) It is used
23   // mainly //for debugging
24
25   AliAltroBuffer buff(file.Data(),0);
26   Int_t numWords,padNum,rowNum,secNum=0;
27   Int_t value = 0;
28   Int_t zero  = 0;
29   // if (!buff.ReadDataHeader()) {
30   // std::cout<< file << " isn't a valid data file!" << std::endl;
31   // }
32   
33   while(buff.ReadTrailerBackward(numWords,padNum,rowNum,secNum)){
34     if (verbose) 
35       std::cout << "Ring: " << (secNum == 0 ? 'I' : 'O') 
36                 << " Sector: " << std::setw(2) << rowNum 
37                 << " Strip:  " << std::setw(3) << padNum 
38                 << " Words:  " << std::setw(4) << numWords << std::endl;
39     if (numWords == 0) zero++;
40     if (numWords % 4){
41       if (verbose) 
42         std::cout << "Skipping trailer of " 
43                   << (4 - numWords % 4) << " words" << std::endl;
44       for(Int_t j = 0; j < (4 - numWords % 4); j++)
45         value=buff.GetNextBackWord(); 
46     }//end if
47     for(Int_t i = 0; i <numWords; i++) {
48       value=buff.GetNextBackWord();
49       if (verbose) {
50         std::cout << std::setw(5) <<  value << std::flush;
51         if (i % 16 == 15) std::cout << std::endl;
52       }
53       h->Fill(value);
54     }//end for
55     if (verbose)
56       std::cout << std::endl;
57     if (zero > 1) {
58       std::cout << "Error: Read zero channels - should not happen" 
59                 << std::endl;
60       break;
61     }
62   }//end while
63
64   TCanvas* c = new TCanvas("raw", "Raw Data");
65   c->SetFillColor(0);
66   c->SetLogy();
67   c->SetBorderMode(0);
68   h->Draw();
69   return;
70 }
71 //____________________________________________________________________
72 //
73 // EOF
74 //