Some of the coding violations corrected
[u/mrichter/AliRoot.git] / TPC / AliTPCDDLRawData.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-2003, 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 /* $Id$ */
16
17
18 //This class conteins all the methods to create raw data 
19 //as par a given DDL.
20 //It produces DDL with both compressed and uncompressed format.
21 //For compression we use the optimized table wich needs 
22 //to be provided.
23
24 #include <TObjArray.h>
25 #include <TString.h>
26 #include <TSystem.h>
27 #include <Riostream.h>
28 #include <stdio.h>
29 #include <stdlib.h>
30 //#include "AliAltroBuffer.h"
31 #include "AliAltroBufferV3.h"
32 #include "AliTPCAltroMapping.h"
33 #include "AliTPCDDLRawData.h"
34 #include "AliDAQ.h"
35
36 ClassImp(AliTPCDDLRawData)
37 ////////////////////////////////////////////////////////////////////////////////////////
38
39 AliTPCDDLRawData::AliTPCDDLRawData(const AliTPCDDLRawData &source):
40   TObject(source),
41   fVerbose(0)
42 {
43   // Copy Constructor
44   fVerbose=source.fVerbose;
45   return;
46 }
47
48 AliTPCDDLRawData& AliTPCDDLRawData::operator=(const AliTPCDDLRawData &source){
49   //Assigment operator
50   fVerbose=source.fVerbose;
51   return *this;
52 }
53
54
55 ////////////////////////////////////////////////////////////////////////////
56 void AliTPCDDLRawData::RawData(const char* inputFileName){
57   //Raw data generation
58   //Number of DDL=2*36+4*36=216
59   //2 DDL for each inner sector
60   //4 DDL for each outer sector
61   ifstream f;
62 #ifndef __DECCXX
63   f.open(inputFileName,ios::binary);
64 #else
65   f.open(inputFileName);
66 #endif
67   if(!f){Error("RawData", "File doesn't exist !!");return;}
68   struct DataPad{
69     Int_t Sec;
70     Int_t SubSec;
71     Int_t Row;
72     Int_t Pad;
73     Int_t Dig;
74     Int_t Time;
75   };
76   DataPad data;
77
78   //AliAltroBuffer is used in write mode to generate raw data file
79   char  filename[101];
80   Int_t ddlNumber=0;
81   AliAltroBuffer *buffer=NULL;
82   Int_t pSecNumber=-1;  //Previous Sector number
83   Int_t pRowNumber=-1;  //Previous Row number  
84   Int_t pPadNumber=-1;  //Previous Pad number
85   Int_t pTimeBin=-1;    //Previous Time-Bin
86   Int_t pSubSector=-1;  //Previous Sub Sector
87   Int_t bunchLength=0;
88   Int_t nwords=0;
89   UInt_t numPackets=0;
90
91   TString path = gSystem->Getenv("ALICE_ROOT");
92   path += "/TPC/mapping/Patch";
93   TString path2;
94   AliTPCAltroMapping *mapping[6];
95   for(Int_t i = 0; i < 6; i++) {
96     path2 = path;
97     path2 += i;
98     path2 += ".data";
99     mapping[i] = new AliTPCAltroMapping(path2.Data());
100   }
101
102
103   while (f.read((char*)(&data),sizeof(data))){
104     if (pPadNumber==-1){
105       pSecNumber=data.Sec;
106       pRowNumber=data.Row;
107       pPadNumber=data.Pad;
108       pTimeBin=data.Time;
109       pSubSector=data.SubSec;
110
111       if(data.Sec<36)
112         ddlNumber=data.Sec*2+data.SubSec;
113       else
114         ddlNumber=72+(data.Sec-36)*4+data.SubSec;
115       strncpy(filename,AliDAQ::DdlFileName("TPC",ddlNumber),100);
116       Int_t patchIndex = data.SubSec;
117       if(data.Sec>=36) patchIndex += 2;
118       //buffer=new AliAltroBuffer(filename,mapping[patchIndex]);
119       buffer=new AliAltroBufferV3(filename,mapping[patchIndex]);
120       //size magic word sector number sub-sector number 0 for TPC 0 for uncompressed
121       buffer->WriteDataHeader(kTRUE,kFALSE);//Dummy;
122       bunchLength=1;
123       buffer->FillBuffer(data.Dig);
124       nwords++;
125     }//end if
126     else{
127       if ( (data.Time==(pTimeBin+1)) &&
128            (pPadNumber==data.Pad) &&
129            (pRowNumber==data.Row) &&
130            (pSecNumber==data.Sec) &&
131            (pSubSector==data.SubSec)){
132         bunchLength++;
133       }//end if
134       else{
135         buffer->FillBuffer(pTimeBin);
136         buffer->FillBuffer(bunchLength+2);
137         nwords+=2;
138         if ((pPadNumber!=data.Pad)||(pRowNumber!=data.Row)||(pSecNumber!=data.Sec)||(pSubSector!=data.SubSec)){
139           //Trailer is formatted and inserted!!
140           buffer->WriteTrailer(nwords,pPadNumber,pRowNumber,pSecNumber);
141           numPackets++;
142           nwords=0;
143
144           if(pSecNumber!=data.Sec || pSubSector!=data.SubSec){
145             //size magic word sector number sub-sector number 0 for TPC 0 for uncompressed
146             buffer->Flush();
147             buffer->WriteDataHeader(kFALSE,kFALSE);
148             //cout<<"Data header for DDL:"<<PSecNumber<<" Sub-sec:"<<PSubSector<<endl;
149             delete buffer;
150
151             if(data.Sec<36)
152               ddlNumber=data.Sec*2+data.SubSec;
153             else
154               ddlNumber=72+(data.Sec-36)*4+data.SubSec;
155             strncpy(filename,AliDAQ::DdlFileName("TPC",ddlNumber),100);
156             Int_t patchIndex = data.SubSec;
157             if(data.Sec>=36) patchIndex += 2;
158             // buffer=new AliAltroBuffer(filename,mapping[patchIndex]);
159             buffer=new AliAltroBufferV3(filename,mapping[patchIndex]);
160             buffer->WriteDataHeader(kTRUE,kFALSE);//Dummy;
161             pSubSector=data.SubSec;
162           }//end if
163         }//end if
164         
165         bunchLength=1;
166         pPadNumber=data.Pad;
167         pRowNumber=data.Row;
168         pSecNumber=data.Sec;
169       }//end else
170       pTimeBin=data.Time;
171       buffer->FillBuffer(data.Dig);
172       nwords++;
173     }//end else
174   }//end while
175   if (buffer) {
176     buffer->FillBuffer(pTimeBin);
177     buffer->FillBuffer(bunchLength+2);
178     nwords+=2;
179     buffer->WriteTrailer(nwords,pPadNumber,pRowNumber,pSecNumber);
180     //write the  D.H.
181     buffer->Flush();
182     buffer->WriteDataHeader(kFALSE,kFALSE);
183     //cout<<"Data header for D D L:"<<pSecNumber<<" Sub-sec:"<<pSubSector<<endl;
184     delete buffer;
185   }
186
187   for(Int_t i = 0; i < 6; i++) delete mapping[i];
188
189   f.close();
190   return;
191 }