3 // Author: C. Loizides <loizides@ikf.uni-frankfurt.de>
4 //*-- Copyright © ALICE HLT Group
6 #include "AliL3StandardIncludes.h"
8 #include "AliL3RootTypes.h"
9 #include "AliL3Logging.h"
10 #include "AliL3Transform.h"
11 #include "AliL3MemHandler.h"
12 #include "AliL3DigitData.h"
14 #include "../RAW/AliRawReaderRoot.h"
15 #include "../RAW/AliRawReaderDate.h"
17 #include "AliL3DDLTPCRawStream.h"
18 #include "AliL3DDLRawReaderFile.h"
20 #include "AliL3DDLDataFileHandler.h"
26 /** \class AliL3DDLDataFileHandler
28 //_____________________________________________________________
29 // AliL3DDLDataFileHandler
31 // This class does converts from the DDL format of offline
32 // into the memory I/O handling of the HLT binary files.
34 // Examples: see ddl2binary in exa and the general
35 // AliL3MemHandler class description
40 ClassImp(AliL3DDLDataFileHandler)
42 AliL3DDLDataFileHandler::AliL3DDLDataFileHandler()
44 // default constructor
49 AliL3DDLDataFileHandler::~AliL3DDLDataFileHandler()
55 void AliL3DDLDataFileHandler::FreeAll()
57 // frees all heap memory
58 if(fReader) delete fReader;
59 if(fTPCStream) delete fTPCStream;
66 Bool_t AliL3DDLDataFileHandler::SetReaderInput(Char_t *name,Int_t event)
68 // sets the input of the reader
70 if(fReader) delete fReader;
71 if(fTPCStream) delete fTPCStream;
74 fReader=new AliRawReaderRoot(name,event);
77 fReader=new AliRawReaderDate((void *)name);
79 fTPCStream=new AliTPCRawStream(fReader);
84 Bool_t AliL3DDLDataFileHandler::SetReaderInput(Char_t *name, Bool_t add)
86 // sets the input of the reader
88 LOG(AliL3Log::kError,"AliL3DDLDataFileHandler::SetReaderInput","File Open")
89 <<"Reader ptr is already in use"<<ENDLOG;
93 fReader=new AliL3DDLRawReaderFile(name,add);
94 fTPCStream=new AliL3DDLTPCRawStream(fReader);
98 Bool_t AliL3DDLDataFileHandler::SetReaderInput(AliL3DDLRawReaderFile *rf)
100 // sets the input of the reader
102 LOG(AliL3Log::kError,"AliL3RawDataFileHandler::SetReaderInput","File Open")
103 <<"Reader ptr is already in use, delete it first"<<ENDLOG;
107 //Open the raw data file with given file.
109 fTPCStream=new AliL3DDLTPCRawStream(fReader);
115 void AliL3DDLDataFileHandler::CloseReaderInput()
117 // closes the input of the reader
119 LOG(AliL3Log::kWarning,"AliL3RawDataFileHandler::CloseReaderInput","File Close")
120 <<"Nothing to Close"<<ENDLOG;
131 Bool_t AliL3DDLDataFileHandler::IsDigit(Int_t /*i*/) const
139 AliL3DigitRowData * AliL3DDLDataFileHandler::DDLData2Memory(UInt_t &nrow,Int_t event)
141 // transfers the DDL data to the memory
143 if((fEvent>=0)&&(event!=fEvent)){
145 if(fReader) delete fReader;
146 if(fTPCStream) delete fTPCStream;
147 fReader=new AliRawReaderRoot(fFilename,event);
148 fTPCStream=new AliTPCRawStream(fReader);
151 AliL3DigitRowData *data = 0;
155 LOG(AliL3Log::kWarning,"AliL3DDLDataFileHandler::DDLData2Memory","File")
156 <<"No Input avalible: no object AliL3DDLRawReaderFile"<<ENDLOG;
160 Int_t nrows=fRowMax-fRowMin+1;
162 Int_t * ndigits = new Int_t[nrows];
163 UShort_t ***charges=new UShort_t**[nrows];
164 for(Int_t r=fRowMin;r<=fRowMax;r++){
165 Int_t lrow=r-fRowMin;
166 charges[lrow]=new UShort_t*[AliL3Transform::GetNPads(r)];
167 for(Int_t k=0;k<AliL3Transform::GetNPads(r);k++){
168 charges[lrow][k]=new UShort_t[AliL3Transform::GetNTimeBins()];
169 for(Int_t j=0;j<AliL3Transform::GetNTimeBins();j++) charges[lrow][k][j]=0;
173 Int_t ddlsToSearch=0;
174 Int_t ddls[9]={-1,-1,-1,-1,-1,-1,-1,-1,-1};
175 Int_t ddlid=-1,lddlid=-1;
176 for(Int_t r=fRowMin;r<=fRowMax;r++){
177 ndigits[r-fRowMin] = 0; //now digits on row
179 Int_t patch=AliL3Transform::GetPatch(r);
181 AliL3Transform::Slice2Sector(fSlice,r,sector,row);
183 if(sector<36) //taken from AliTPCBuffer160.cxx
184 ddlid=sector*2+patch;
186 ddlid=70+(sector-36)*4+patch;
188 if((lddlid!=ddlid-1)&&(r==30)){ //dont forget the split row on the last ddl
189 ddls[ddlsToSearch++]=ddlid-1;
193 if((lddlid==-1)||(ddlid!=lddlid)){
194 ddls[ddlsToSearch++]=ddlid;
197 if((r==90)||(r==139)){ //dont forget the split row on the next ddl
198 ddls[ddlsToSearch++]=ddlid+1;
203 // for(Int_t i=0;i<ddlsToSearch;i++) cout << ddls[i] <<endl;
205 if(ddls[0]>ddls[ddlsToSearch-1]) {
206 Int_t tempddl = ddls[0];
207 ddls[0] = ddls[ddlsToSearch-1];
208 ddls[ddlsToSearch-1] = tempddl;
212 fReader->Select(0,ddls[0],ddls[ddlsToSearch-1]);
215 fTPCStream->SetDDLID(ddls[i]); //ddl to read out
217 Int_t zerosup = AliL3Transform::GetZeroSup();
218 Int_t adcsat = AliL3Transform::GetADCSat();
222 while (fTPCStream->Next()){
224 if(fTPCStream->IsNewSector() || fTPCStream->IsNewRow()) {
225 Int_t sector=fTPCStream->GetSector();
226 Int_t row=fTPCStream->GetRow();
227 AliL3Transform::Sector2Slice(slice,srow,sector,row);
229 LOG(AliL3Log::kError,"AliL3DDLDataFileHandler::DDLDigits2Memory","Slice")
230 <<AliL3Log::kDec<<"Found slice "<<slice<<", expected "<<fSlice<<ENDLOG;
236 //test row criteria (patch boundaries)
237 if((srow<fRowMin)||(srow>fRowMax))continue;
239 Int_t pad=fTPCStream->GetPad();
240 if(fTPCStream->IsNewPad()) {
241 if((pad<0)||(pad>=AliL3Transform::GetNPads(srow))){
242 LOG(AliL3Log::kError,"AliL3DDLDataFileHandler::DDLDigits2Memory","Pad")
243 <<AliL3Log::kDec<<"Pad value out of bounds "<<pad<<" "
244 <<AliL3Transform::GetNPads(srow)<<ENDLOG;
249 Int_t time=fTPCStream->GetTime();
250 if((time<0)||(time>=AliL3Transform::GetNTimeBins())){
251 LOG(AliL3Log::kError,"AliL3DDLDataFileHandler::DDLDigits2Memory","Time")
252 <<AliL3Log::kDec<<"Time out of bounds "<<time<<" "
253 <<AliL3Transform::GetNTimeBins()<<ENDLOG;
258 UShort_t dig=fTPCStream->GetSignal();
259 if(dig <= zerosup) continue;
260 if(dig >= adcsat) dig = adcsat;
262 ndigits[lrow]++; //for this row only
263 ndigitcount++; //total number of digits to be published
265 charges[lrow][pad][time]=dig;
268 Int_t size = sizeof(AliL3DigitData)*ndigitcount
269 + nrows*sizeof(AliL3DigitRowData);
271 LOG(AliL3Log::kDebug,"AliL3DDLDataFileHandler::DDLDigits2Memory","Digits")
272 <<AliL3Log::kDec<<"Found "<<ndigitcount<<" Digits"<<ENDLOG;
274 data=(AliL3DigitRowData*) Allocate(size);
275 nrow = (UInt_t)nrows;
276 AliL3DigitRowData *tempPt = data;
278 for(Int_t r=fRowMin;r<=fRowMax;r++){
279 Int_t lrow=r-fRowMin;
281 tempPt->fNDigit = ndigits[lrow];
284 for(Int_t pad=0;pad<AliL3Transform::GetNPads(r);pad++){
285 for(Int_t time=0;time<AliL3Transform::GetNTimeBins();time++){
286 UShort_t dig=charges[lrow][pad][time];
289 if(localcount >= ndigits[lrow])
290 LOG(AliL3Log::kFatal,"AliL3DDLDataFileHandler::DDLDigits2Binary","Memory")
291 <<AliL3Log::kDec<<"Mismatch: localcount "<<localcount<<" ndigits "
292 <<ndigits[lrow]<<ENDLOG;
295 tempPt->fDigitData[localcount].fCharge=dig;
296 tempPt->fDigitData[localcount].fPad=pad;
297 tempPt->fDigitData[localcount].fTime=time;
299 tempPt->fDigitData[localcount].fTrackID[0] = 0;
300 tempPt->fDigitData[localcount].fTrackID[1] = 0;
301 tempPt->fDigitData[localcount].fTrackID[2] = 0;
307 if(localcount != ndigits[lrow])
308 LOG(AliL3Log::kFatal,"AliL3DDLDataFileHandler::DDLDigits2Binary","Memory")
309 <<AliL3Log::kDec<<"Mismatch: localcount "<<localcount<<" ndigits "
310 <<ndigits[lrow]<<ENDLOG;
313 Byte_t *tmp = (Byte_t*)tempPt;
314 Int_t size = sizeof(AliL3DigitRowData)
315 + ndigits[lrow]*sizeof(AliL3DigitData);
317 tempPt = (AliL3DigitRowData*)tmp;
320 //delete charge array
321 for(Int_t r=fRowMin;r<=fRowMax;r++){
322 Int_t lrow=r-fRowMin;
323 for(Int_t k=0;k<AliL3Transform::GetNPads(r);k++)
324 delete charges[lrow][k];
325 delete charges[lrow];
333 AliL3DigitRowData * AliL3DDLDataFileHandler::DDLData2Memory(UInt_t &nrow,Int_t event)
335 // transfers the DDL data to the memory
337 if((fEvent>=0)&&(event!=fEvent)){
339 if(fReader) delete fReader;
340 if(fTPCStream) delete fTPCStream;
341 fReader=new AliRawReaderRoot(fFilename,event);
342 fTPCStream=new AliTPCRawStream(fReader);
345 AliL3DigitRowData *data = 0;
349 LOG(AliL3Log::kWarning,"AliL3DDLDataFileHandler::DDLData2Memory","File")
350 <<"No Input avalible: no object AliL3DDLRawReaderFile"<<ENDLOG;
354 Int_t nrows=fRowMax-fRowMin+1;
356 Int_t ddlsToSearch=0;
357 Int_t ddls[9]={-1,-1,-1,-1,-1,-1,-1,-1,-1};
358 Int_t ddlid=-1,lddlid=-1;
359 for(Int_t r=fRowMin;r<=fRowMax;r++){
361 Int_t patch=AliL3Transform::GetPatch(r);
363 AliL3Transform::Slice2Sector(fSlice,r,sector,row);
365 if(sector<36) //taken from AliTPCBuffer160.cxx
366 ddlid=sector*2+patch;
368 ddlid=70+(sector-36)*4+patch;
370 if((lddlid==-1)||(ddlid!=lddlid)){
371 ddls[ddlsToSearch++]=ddlid;
376 // for(Int_t i=0;i<ddlsToSearch;i++) cout << ddls[i] <<endl;
380 fReader->Select(0,ddls[0],ddls[ddlsToSearch-1]);
383 fTPCStream->SetDDLID(ddls[i]); //ddl to read out
386 nrow = (UInt_t)nrows;
392 Bool_t AliL3DDLDataFileHandler::DDLData2CompBinary(Int_t event)
394 // transfers the DDL data to the memory and converts it
398 AliL3DigitRowData *digits=0;
399 digits = DDLData2Memory(ndigits,event);
400 out = Memory2CompBinary(ndigits,digits);