]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TOF/AliTOFDDLRawData.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Fix for coverity and compilation warnings
[u/mrichter/AliRoot.git] / TOF / AliTOFDDLRawData.cxx
CommitLineData
7e6dce66 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
16/*
d0eb8f39 17$Log$
96f01799 18Revision 1.19 2007/06/22 11:37:47 cvetan
19Fixes in order to write correct raw-data on big-endian platforms (Marco)
20
08f92f14 21Revision 1.18 2007/05/21 13:26:19 decaro
22Correction on matching_window control and bug fixing (R.Preghenella)
23
018d19cf 24Revision 1.17 2007/05/10 09:29:34 hristov
25Last moment fixes and changes from v4-05-Release (Silvia)
26
d321691a 27Revision 1.16 2007/05/03 09:07:22 decaro
28Double digit in the same TDC channel. Wrong sequence during the raw data writing in unpacked mode: solved
29
a403781d 30Revision 1.15 2007/04/23 16:51:39 decaro
31Digits-to-raw_data conversion: correction for a more real description (A.De Caro, R.Preghenella)
32
515faf5d 33Revision 1.14 2007/03/28 10:50:33 decaro
34Rounding off problem in rawData coding/decoding: solved
35
d11fc181 36Revision 1.13 2007/02/20 15:57:00 decaro
37Raw data update: to read the TOF raw data defined in UNPACKED mode
38
15ec34b9 39Revision 1.12 2006/08/22 13:29:42 arcelli
40removal of effective c++ warnings (C.Zampolli)
41
655e379f 42Revision 1.11 2006/08/10 14:46:54 decaro
43TOF raw data format: updated version
44
d0eb8f39 45Revision 1.10.1 2006/06/28 A.De Caro
46 Update TOF raw data format
47 according to the final version
48 (see ALICE internal note in preparation
49 'ALICE TOF raw data format')
50
571dda3d 51Revision 0.02 2005/7/25 A.De Caro
52 Update number of bits allocated for time-of-flight
53 and 'charge' measurements
54
7e6dce66 55Revision 0.01 2004/6/11 A.De Caro, S.B.Sellitto, R.Silvestri
56 First implementation: global methods RawDataTOF
57 GetDigits
7e6dce66 58*/
571dda3d 59
d0eb8f39 60////////////////////////////////////////////////////////////////////
61// //
62// This class contains the methods to create the Raw Data files //
63// for the TOF detector starting from the Digits. //
64// In this implementation, we defined the structure //
65// of the ALICE-TOF raw data (according to the //
66// ALICE technical note, in preparation) //
67// starting from the TOF digit format. //
68// //
69////////////////////////////////////////////////////////////////////
7e6dce66 70
0e46b9ae 71#include "Riostream.h"
72
73#include "TBranch.h"
74#include "TClonesArray.h"
75#include "TMath.h"
d0eb8f39 76#include "TRandom.h"
43f77f2d 77
571dda3d 78#include "AliBitPacking.h"
d0eb8f39 79#include "AliDAQ.h"
0e46b9ae 80#include "AliLog.h"
d0eb8f39 81//#include "AliRawDataHeader.h"
82#include "AliRawDataHeaderSim.h"
32ead898 83#include "AliFstream.h"
571dda3d 84
0e46b9ae 85#include "AliTOFDDLRawData.h"
d0eb8f39 86#include "AliTOFDigitMap.h"
7e6dce66 87#include "AliTOFdigit.h"
0e46b9ae 88#include "AliTOFGeometry.h"
571dda3d 89#include "AliTOFRawStream.h"
32ead898 90//#include "AliTOFCableLengthMap.h"
d0eb8f39 91
92extern TRandom *gRandom;
7e6dce66 93
94ClassImp(AliTOFDDLRawData)
95
d0eb8f39 96//---------------------------------------------------------------------------
655e379f 97AliTOFDDLRawData::AliTOFDDLRawData():
6e4ef37b 98 TObject(),
655e379f 99 fVerbose(0),
100 fIndex(-1),
32ead898 101 fPackedAcquisition(kFALSE),
515faf5d 102 fFakeOrphaneProduction(kFALSE),
103 fMatchingWindow(8192),
655e379f 104 fTOFdigitMap(new AliTOFDigitMap()),
177b002e 105 fTOFdigitArray(0x0)
106 //fTOFrawStream(AliTOFRawStream()),
107 //fTOFCableLengthMap(new AliTOFCableLengthMap())
7e6dce66 108{
109 //Default constructor
d3c7bfac 110}
d0eb8f39 111//----------------------------------------------------------------------------
112AliTOFDDLRawData::AliTOFDDLRawData(const AliTOFDDLRawData &source) :
655e379f 113 TObject(source),
8a190ba2 114 fVerbose(source.fVerbose),
115 fIndex(source.fIndex),
116 fPackedAcquisition(source.fPackedAcquisition),
117 fFakeOrphaneProduction(source.fFakeOrphaneProduction),
118 fMatchingWindow(source.fMatchingWindow),
119 fTOFdigitMap(source.fTOFdigitMap),
177b002e 120 fTOFdigitArray(source.fTOFdigitArray)
121 //fTOFrawStream(source.fTOFrawStream),
122 //fTOFCableLengthMap(source.fTOFCableLengthMap)
655e379f 123 {
7e6dce66 124 //Copy Constructor
7e6dce66 125 return;
126}
127
d0eb8f39 128//----------------------------------------------------------------------------
129AliTOFDDLRawData& AliTOFDDLRawData::operator=(const AliTOFDDLRawData &source) {
7e6dce66 130 //Assigment operator
8a190ba2 131
132 if (this == &source)
133 return *this;
134
135 TObject::operator=(source);
136 fIndex=source.fIndex;
137 fVerbose=source.fVerbose;
138 fPackedAcquisition=source.fPackedAcquisition;
139 fFakeOrphaneProduction=source.fFakeOrphaneProduction;
140 fMatchingWindow=source.fMatchingWindow;
141 fTOFdigitMap=source.fTOFdigitMap;
142 fTOFdigitArray=source.fTOFdigitArray;
177b002e 143 //fTOFrawStream=source.fTOFrawStream;
144 //fTOFCableLengthMap=source.fTOFCableLengthMap;
7e6dce66 145 return *this;
146}
147
d321691a 148//----------------------------------------------------------------------------
149AliTOFDDLRawData::~AliTOFDDLRawData()
150{
151 delete fTOFdigitMap;
177b002e 152 //delete fTOFrawStream;
153 //delete fTOFCableLengthMap;
d321691a 154}
d0eb8f39 155//----------------------------------------------------------------------------
156Int_t AliTOFDDLRawData::RawDataTOF(TBranch* branch)
7e6dce66 157{
571dda3d 158 //
159 // This method creates the Raw data files for TOF detector
160 //
161
6e4ef37b 162 const Int_t kSize = 5000; // max number of digits per DDL file times 2
571dda3d 163
7e6dce66 164 UInt_t buf[kSize];
7e6dce66 165
6e4ef37b 166 // To clear the digit indices map for each event
167 fTOFdigitMap->Clear();
168
177b002e 169 //fTOFrawStream->Clear();
170
d0eb8f39 171 fIndex = -1;
172
173 fTOFdigitArray = * (TClonesArray**) branch->GetAddress();
7e6dce66 174
abfcded3 175 const Int_t kDimension=16;
176 char fileName[kDimension-1];
6e4ef37b 177 AliFstream* outfile; // logical name of the output file
d0eb8f39 178
179 //AliRawDataHeader header;
180 AliRawDataHeaderSim header;
181
571dda3d 182 UInt_t sizeRawData = 0;
7e6dce66 183
d0eb8f39 184 branch->GetEvent();
185
186 GetDigits();
187
3e2f0097 188 //if (!fPackedAcquisition) fMatchingWindow = 2097152;//AdC
189
d0eb8f39 190 Int_t jj = -1;
191 Int_t iDDL = -1;
192 Int_t nDDL = -1;
193 Int_t nTRM = 0;
194 Int_t iChain = -1;
195
196 UInt_t nWordsPerTRM = 0;
197
7e6dce66 198 //loop over TOF DDL files
d0eb8f39 199 for (nDDL=0; nDDL<AliDAQ::NumberOfDdls("TOF"); nDDL++) {
571dda3d 200
de40f9c1 201 strncpy(fileName,AliDAQ::DdlFileName("TOF",nDDL),kDimension); //The name of the output file
571dda3d 202
08f92f14 203 outfile = new AliFstream(fileName);
177b002e 204 //iDDL = fTOFrawStream->GetDDLnumberPerSector(nDDL);
205 iDDL = AliTOFRawStream::GetDDLnumberPerSector(nDDL);
d0eb8f39 206
207 // write Dummy DATA HEADER
08f92f14 208 UInt_t dataHeaderPosition = outfile->Tellp();
209 outfile->WriteBuffer((char*)(&header),sizeof(header));
7e6dce66 210
d0eb8f39 211 // DRM section: trailer
212 MakeDRMtrailer(buf);
7e6dce66 213
d0eb8f39 214 // LTM section
7e6dce66 215 fIndex++;
d0eb8f39 216 buf[fIndex] = MakeFiller();
217 MakeLTMtrailer(buf);
218 MakeLTMdata(buf);
219 MakeLTMheader(buf);
220
221 // loop on TRM number per DRM
222 for (nTRM=AliTOFGeometry::NTRM(); nTRM>=3; nTRM--) {
223
224 nWordsPerTRM = 0;
225
226 // the slot number 3 of the even DRM (i.e. right) doesn't contain TDC digit data
227 if (iDDL%2==0 && nTRM==3) continue;
228
229 // TRM global trailer
230 MakeTRMtrailer(buf);
231
232 // loop on TRM chain number per TRM
233 for (iChain=AliTOFGeometry::NChain()-1; iChain>=0; iChain--) {
234
235 // TRM chain trailer
236 MakeTRMchainTrailer(iChain, buf);
237 nWordsPerTRM++;
238
239 // TRM TDC digits
240 MakeTDCdigits(nDDL, nTRM, iChain, buf, nWordsPerTRM);
241
242 // TRM chain header
243 MakeTRMchainHeader(nTRM, iChain, buf);
244 nWordsPerTRM++;
245
246 } // end loop on iChain
7e6dce66 247
d0eb8f39 248 // TRM global header
249 MakeTRMheader(nTRM, buf);
250
251 // TRM filler in case where TRM data number is odd
252 if (nWordsPerTRM%2!=0) MakeTRMfiller(buf, nWordsPerTRM);
253
254 } // end loop on nTRM
255
256 // DRM section: in
257 MakeDRMheader(nDDL, buf);
258
259 ReverseArray(buf, fIndex+1);
7e6dce66 260
08f92f14 261 outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t)));
7e6dce66 262
d0eb8f39 263 for (jj=0; jj<(fIndex+1); jj++) buf[jj]=0;
264 fIndex = -1;
7e6dce66 265
266 //Write REAL DATA HEADER
08f92f14 267 UInt_t currentFilePosition = outfile->Tellp();
571dda3d 268 sizeRawData = currentFilePosition - dataHeaderPosition - sizeof(header);
d0eb8f39 269 header.fSize = currentFilePosition - dataHeaderPosition;
7e6dce66 270 header.SetAttribute(0); // valid data
08f92f14 271 outfile->Seekp(dataHeaderPosition);
272 outfile->WriteBuffer((char*)(&header),sizeof(header));
273 outfile->Seekp(currentFilePosition);
571dda3d 274
08f92f14 275 delete outfile;
7e6dce66 276
d0eb8f39 277 } //end loop on DDL file number
278
279 return 0;
280
281}
282
283//----------------------------------------------------------------------------
284void AliTOFDDLRawData::GetDigits()
285{
286 //
287 // Fill the TOF volumes' map with the TOF digit indices
288 //
289
290 Int_t vol[5] = {-1,-1,-1,-1,-1};
291
292 Int_t digit = -1;
293 Int_t ndigits = fTOFdigitArray->GetEntries();
96b13296 294 AliDebug(2, Form(" Number of read digits = %d",ndigits));
d0eb8f39 295 AliTOFdigit *digs;
296
297 // loop on TOF digits
298 for (digit=0; digit<ndigits; digit++) {
299 digs = (AliTOFdigit*)fTOFdigitArray->UncheckedAt(digit);
300
301 vol[0] = digs->GetSector(); // Sector Number (0-17)
302 vol[1] = digs->GetPlate(); // Plate Number (0-4)
303 vol[2] = digs->GetStrip(); // Strip Number (0-14/18)
304 vol[3] = digs->GetPadx(); // Pad Number in x direction (0-47)
305 vol[4] = digs->GetPadz(); // Pad Number in z direction (0-1)
306
307 fTOFdigitMap->AddDigit(vol, digit);
308
309 } // close loop on digit del TOF
96b13296 310 AliDebug(2,Form(" Number of mapped digits = %d",fTOFdigitMap->GetFilledCellNumber()));
d0eb8f39 311}
312
313//----------------------------------------------------------------------------
314void AliTOFDDLRawData::MakeDRMheader(Int_t nDDL, UInt_t *buf)
315{
316 //
317 // DRM global header
318 //
319
177b002e 320 //Int_t iDDL = fTOFrawStream->GetDDLnumberPerSector(nDDL);
321 Int_t iDDL = AliTOFRawStream::GetDDLnumberPerSector(nDDL);
d0eb8f39 322
177b002e 323 //Int_t iSector = fTOFrawStream->GetSectorNumber(nDDL);
324 Int_t iSector = AliTOFRawStream::GetSectorNumber(nDDL);
d0eb8f39 325
326 UInt_t baseWord=0;
327 UInt_t word;
328
329 // DRM event CRC
330 baseWord=0;
331 word = 1; // 0001 -> DRM data are coming from the VME slot number 1
332 AliBitPacking::PackWord(word,baseWord, 0, 3);
333 word = 0; // event CRC
334 AliBitPacking::PackWord(word,baseWord, 4,19);
335 word = 0; // reserved for future use
336 AliBitPacking::PackWord(word,baseWord,20,27);
337 word = 4; // 0100 -> DRM header ID
338 AliBitPacking::PackWord(word,baseWord,28,31);
339 fIndex++;
340 buf[fIndex]=baseWord;
341
73620dbf 342 // DRM status header 4
343 baseWord=0;
344 word = 1; // 0001 -> DRM data are coming from the VME slot number 1
345 AliBitPacking::PackWord(word,baseWord, 0, 3);
346 word = 0; // temperature
347 AliBitPacking::PackWord(word,baseWord, 4,13);
348 word = 0; // zero
349 AliBitPacking::PackWord(word,baseWord, 14,14);
350 word = 0; // ACK
351 AliBitPacking::PackWord(word,baseWord, 15,15);
352 word = 0; // Sens AD
353 AliBitPacking::PackWord(word,baseWord, 16,18);
354 word = 0; // zero
355 AliBitPacking::PackWord(word,baseWord, 19,19);
356 word = 0; // reserved for future use
357 AliBitPacking::PackWord(word,baseWord, 20,27);
358 word = 4; // 0100 -> DRM header ID
359 AliBitPacking::PackWord(word,baseWord,28,31);
360 fIndex++;
361 buf[fIndex]=baseWord;
362
d0eb8f39 363 // DRM status header 3
364 baseWord=0;
365 word = 1; // 0001 -> DRM data are coming from the VME slot number 1
366 AliBitPacking::PackWord(word,baseWord, 0, 3);
73620dbf 367 word = 0; // L0 BCID
368 AliBitPacking::PackWord(word,baseWord, 4,15);
369 word = 0; // Run Time info
370 AliBitPacking::PackWord(word,baseWord, 16,27);
d0eb8f39 371 word = 4; // 0100 -> DRM header ID
372 AliBitPacking::PackWord(word,baseWord,28,31);
373 fIndex++;
374 buf[fIndex]=baseWord;
375
376 // DRM status header 2
377 baseWord=0;
378 word = 1; // 0001 -> DRM data are coming from the VME slot number 1
379 AliBitPacking::PackWord(word,baseWord, 0, 3);
380
381 if (iDDL%2==1) {
382 word = 2047; // enable ID: [00000000000;11111111111] for odd
383 // (i.e. right) crates
384 AliBitPacking::PackWord(word,baseWord, 4,14);
385 } else {
386 word = 2045; // enable ID: [00000000000;11111111101] for even
387 // (i.e. left) crates
388 AliBitPacking::PackWord(word,baseWord, 4,14);
389 }
390
391 word = 0; //
392 AliBitPacking::PackWord(word,baseWord,15,15);
393 word = 0; // fault ID
73620dbf 394 AliBitPacking::PackWord(word,baseWord,16,26);
395 word = 0; // RTO
396 AliBitPacking::PackWord(word,baseWord,27,27);
d0eb8f39 397 word = 4; // 0100 -> DRM header ID
398 AliBitPacking::PackWord(word,baseWord,28,31);
399 fIndex++;
400 buf[fIndex]=baseWord;
7e6dce66 401
d0eb8f39 402 // DRM status header 1
403 baseWord=0;
404 word = 1; // 0001 -> DRM data are coming from the VME slot number 1
405 AliBitPacking::PackWord(word,baseWord, 0, 3);
406
407 if (iDDL%2==1) {
408 word = 2047; // slot ID: [00000000000;11111111111] for odd
409 // (i.e. right) crates
410 AliBitPacking::PackWord(word,baseWord, 4,14);
411 } else {
412 word = 2045; // slot ID: [00000000000;11111111101] for even
413 // (i.e. left) crates
414 AliBitPacking::PackWord(word,baseWord, 4,14);
415 }
416
417 word = 1; // LHC clock status: 1/0
418 AliBitPacking::PackWord(word,baseWord,15,15);
73620dbf 419 word = 0; // Vers ID
420 AliBitPacking::PackWord(word,baseWord,16,20);
421 word = 0; // DRMH size
422 AliBitPacking::PackWord(word,baseWord,21,24);
d0eb8f39 423 word = 0; // reserved for future use
73620dbf 424 AliBitPacking::PackWord(word,baseWord,25,27);
d0eb8f39 425 word = 4; // 0100 -> DRM header ID
426 AliBitPacking::PackWord(word,baseWord,28,31);
427 fIndex++;
428 buf[fIndex]=baseWord;
429
430 // DRM global header
431 baseWord=0;
432 word = 1; // 0001 -> DRM data are coming from the VME slot number 1
433 AliBitPacking::PackWord(word,baseWord, 0, 3);
434 word = fIndex+1 + 1; // event words
435 AliBitPacking::PackWord(word,baseWord, 4,20);
436 word = iDDL; // crate ID [0;3]
437 AliBitPacking::PackWord(word,baseWord,21,22);
438 word = iSector; // sector ID [0;17]
439 AliBitPacking::PackWord(word,baseWord,23,27);
440 word = 4; // 0100 -> DRM header ID
441 AliBitPacking::PackWord(word,baseWord,28,31);
442 fIndex++;
443 buf[fIndex]=baseWord;
444
7e6dce66 445}
446
d0eb8f39 447//----------------------------------------------------------------------------
448void AliTOFDDLRawData::MakeDRMtrailer(UInt_t *buf)
449{
450 //
451 // DRM global trailer
452 //
453
454 UInt_t baseWord;
455 UInt_t word;
456
457 baseWord=0;
458 word = 1; // 0001 -> DRM data are coming from the VME slot number 1
459 AliBitPacking::PackWord(word,baseWord, 0, 3);
460 word = 0; // local event counter
461 AliBitPacking::PackWord(word,baseWord, 4,15);
462 word = 0; // reserved for future use
463 AliBitPacking::PackWord(word,baseWord,16,27);
464 word = 5; // 0101 -> DRM trailer ID
465 AliBitPacking::PackWord(word,baseWord,28,31);
466 fIndex++;
467 buf[fIndex]=baseWord;
468
469}
470
471//----------------------------------------------------------------------------
472void AliTOFDDLRawData::MakeLTMheader(UInt_t *buf)
473{
474 //
475 // LTM header
476 //
477
478 UInt_t baseWord;
479 UInt_t word;
480
481 baseWord=0;
482 word = 2; // 0010 -> LTM data are coming from the VME slot number 2
483 AliBitPacking::PackWord(word,baseWord, 0, 3);
484 word = 35; // event words
485 AliBitPacking::PackWord(word,baseWord, 4,16);
486 word = 0; // crc error
487 AliBitPacking::PackWord(word,baseWord,17,17);
488 word = 0; // fault
489 AliBitPacking::PackWord(word,baseWord,18,23);
490 word = 0;
491 AliBitPacking::PackWord(word,baseWord,24,27);
492 word = 4; // 0100 -> LTM header ID
493 AliBitPacking::PackWord(word,baseWord,28,31);
494 fIndex++;
495 buf[fIndex]=baseWord;
496
497}
498
499//----------------------------------------------------------------------------
500void AliTOFDDLRawData::MakeLTMdata(UInt_t *buf)
501{
502 //
503 // LTM data
504 //
505
506 UInt_t baseWord;
507 UInt_t word;
96b13296 508
d0eb8f39 509 baseWord=0;
96b13296 510 word = 0;
511 AliBitPacking::PackWord(word,baseWord,0,9);
512 word = 0;
d0eb8f39 513 AliBitPacking::PackWord(word,baseWord,10,19);
96b13296 514 word = 0;
d0eb8f39 515 AliBitPacking::PackWord(word,baseWord,20,29);
516 word = 0;
96b13296 517 AliBitPacking::PackWord(word,baseWord,30,30);
518 word = 0;
519 AliBitPacking::PackWord(word,baseWord,31,31);
520
521
522 // OR45, OR46, OR47
523 fIndex++;
524 buf[fIndex]=baseWord;
525
526 // OR42, OR43, OR44
527 fIndex++;
528 buf[fIndex]=baseWord;
529
530 // OR39, OR40, OR41
531 fIndex++;
532 buf[fIndex]=baseWord;
533
534 // OR36, OR37, OR38
535 fIndex++;
536 buf[fIndex]=baseWord;
537
538 // OR33, OR34, OR35
539 fIndex++;
540 buf[fIndex]=baseWord;
541
542 // OR30, OR31, OR32
543 fIndex++;
544 buf[fIndex]=baseWord;
545
546 // OR27, OR28, OR29
547 fIndex++;
548 buf[fIndex]=baseWord;
549
550 // OR24, OR25, OR26
551 fIndex++;
552 buf[fIndex]=baseWord;
553
554 // OR21, OR22, OR23
555 fIndex++;
556 buf[fIndex]=baseWord;
557
558 // OR18, OR19, OR20
559 fIndex++;
560 buf[fIndex]=baseWord;
d0eb8f39 561
96b13296 562 // OR15, OR16, OR17
d0eb8f39 563 fIndex++;
564 buf[fIndex]=baseWord;
565
96b13296 566 // OR12, OR12, OR24
d0eb8f39 567 fIndex++;
568 buf[fIndex]=baseWord;
569
96b13296 570 // OR9, OR10, OR11
d0eb8f39 571 fIndex++;
572 buf[fIndex]=baseWord;
573
96b13296 574 // OR6, OR7, OR8
d0eb8f39 575 fIndex++;
576 buf[fIndex]=baseWord;
577
96b13296 578 // OR3, OR4, OR5
d0eb8f39 579 fIndex++;
580 buf[fIndex]=baseWord;
581
96b13296 582 // OR0, OR1, OR2
583 fIndex++;
584 buf[fIndex]=baseWord;
585
586
587
d0eb8f39 588 baseWord=0;
96b13296 589 word = 100; // Local temperature in LTM11 -> 4 X 25 degree (environment temperature)
d0eb8f39 590 AliBitPacking::PackWord(word,baseWord, 0, 9);
96b13296 591 word = 100; // Local temperature in LTM10 -> 4 X 25 degree (environment temperature)
d0eb8f39 592 AliBitPacking::PackWord(word,baseWord,10,19);
96b13296 593 word = 100; // Local temperature in LTM9 -> 4 X 25 degree (environment temperature)
d0eb8f39 594 AliBitPacking::PackWord(word,baseWord,20,29);
595 word = 0;
96b13296 596 AliBitPacking::PackWord(word,baseWord,30,30);
597 word = 0;
598 AliBitPacking::PackWord(word,baseWord,31,31);
599
600 fIndex++;
601 buf[fIndex]=baseWord;
602
603 // Local temperature in LTM8, LMT7, LTM6 -> 4 X 25 degree (environment temperature)
604 fIndex++;
605 buf[fIndex]=baseWord;
606
607 // Local temperature in LTM5, LMT4, LTM3 -> 4 X 25 degree (environment temperature)
608 fIndex++;
609 buf[fIndex]=baseWord;
610
611 // Local temperature in LTM2, LMT1, LTM0 -> 4 X 25 degree (environment temperature)
612 fIndex++;
613 buf[fIndex]=baseWord;
614
615
d0eb8f39 616
96b13296 617 // Local temperature in T7, T6, T5 -> 4 X 25 degree (environment temperature)
d0eb8f39 618 fIndex++;
619 buf[fIndex]=baseWord;
620
96b13296 621 // Local temperature in T4, T3, T2 -> 4 X 25 degree (environment temperature)
d0eb8f39 622 fIndex++;
623 buf[fIndex]=baseWord;
624
96b13296 625 // Local temperature in T1, T0 -> 4 X 25 degree (environment temperature)
626 // Local temperature in VTH15 -> Thereshould voltage for FEAC15
627 fIndex++;
628 buf[fIndex]=baseWord;
629
630 baseWord=0;
631 word = 0; // VTH13 -> Thereshould voltage for FEAC16
632 AliBitPacking::PackWord(word,baseWord, 0, 9);
633 word = 0; // VTH14 -> Thereshould voltage for FEAC14
634 AliBitPacking::PackWord(word,baseWord,10,19);
635 word = 0; // GND-FEAC7 -> Voltage drop between GND and FEAC7
636 AliBitPacking::PackWord(word,baseWord,20,29);
637 word = 0;
638 AliBitPacking::PackWord(word,baseWord,30,30);
639 word = 0;
640 AliBitPacking::PackWord(word,baseWord,31,31);
641
d0eb8f39 642 fIndex++;
643 buf[fIndex]=baseWord;
644
645 // VTH11 -> Thereshould voltage for FEAC11
d0eb8f39 646 // VTH12 -> Thereshould voltage for FEAC12
96b13296 647 // GND-FEAC6 -> Voltage drop between GND and FEAC6
d0eb8f39 648 fIndex++;
649 buf[fIndex]=baseWord;
650
96b13296 651 // VTH9 -> Thereshould voltage for FEAC11
652 // VTH10 -> Thereshould voltage for FEAC12
653 // GND-FEAC5 -> Voltage drop between GND and FEAC6
d0eb8f39 654 fIndex++;
655 buf[fIndex]=baseWord;
656
96b13296 657 // VTH7 -> Thereshould voltage for FEAC11
658 // VTH8 -> Thereshould voltage for FEAC12
659 // GND-FEAC4 -> Voltage drop between GND and FEAC6
d0eb8f39 660 fIndex++;
661 buf[fIndex]=baseWord;
662
96b13296 663 // VTH5 -> Thereshould voltage for FEAC11
664 // VTH6 -> Thereshould voltage for FEAC12
665 // GND-FEAC3 -> Voltage drop between GND and FEAC6
d0eb8f39 666 fIndex++;
667 buf[fIndex]=baseWord;
668
96b13296 669 // VTH3 -> Thereshould voltage for FEAC11
670 // VTH4 -> Thereshould voltage for FEAC12
671 // GND-FEAC2 -> Voltage drop between GND and FEAC6
d0eb8f39 672 fIndex++;
673 buf[fIndex]=baseWord;
674
96b13296 675 // VTH1 -> Thereshould voltage for FEAC11
676 // VTH2 -> Thereshould voltage for FEAC12
677 // GND-FEAC1 -> Voltage drop between GND and FEAC6
d0eb8f39 678 fIndex++;
679 buf[fIndex]=baseWord;
680
96b13296 681 // LV15
682 // VTH0 -> Thereshould voltage for FEAC12
683 // GND-FEAC0 -> Voltage drop between GND and FEAC6
d0eb8f39 684 fIndex++;
685 buf[fIndex]=baseWord;
686
96b13296 687 // LV15
688 // VTH0 -> Thereshould voltage for FEAC12
689 // GND-FEAC0 -> Voltage drop between GND and FEAC6
d0eb8f39 690 fIndex++;
691 buf[fIndex]=baseWord;
692
96b13296 693 // LV12
694 // LV13
695 // LV14
d0eb8f39 696 fIndex++;
697 buf[fIndex]=baseWord;
698
96b13296 699 // LV9
700 // LV10
701 // LV11
d0eb8f39 702 fIndex++;
703 buf[fIndex]=baseWord;
704
96b13296 705 // LV6
706 // LV7
707 // LV8
d0eb8f39 708 fIndex++;
709 buf[fIndex]=baseWord;
710
96b13296 711 // LV3
712 // LV4
713 // LV5
d0eb8f39 714 fIndex++;
715 buf[fIndex]=baseWord;
716
96b13296 717 // LV0
718 // LV1
719 // LV2
d0eb8f39 720 fIndex++;
721 buf[fIndex]=baseWord;
722
723
96b13296 724
d0eb8f39 725 baseWord=0;
96b13296 726 word = 0; // PDL45 -> Delay Line setting for PDL41
d0eb8f39 727 AliBitPacking::PackWord(word,baseWord, 0, 7);
96b13296 728 word = 0; // PDL46 -> Delay Line setting for PDL42
d0eb8f39 729 AliBitPacking::PackWord(word,baseWord, 8,15);
96b13296 730 word = 0; // PDL47 -> Delay Line setting for PDL43
d0eb8f39 731 AliBitPacking::PackWord(word,baseWord,16,23);
96b13296 732 word = 0; // PDL48 -> Delay Line setting for PDL44
d0eb8f39 733 AliBitPacking::PackWord(word,baseWord,24,31);
734 fIndex++;
735 buf[fIndex]=baseWord;
736
d0eb8f39 737 // Delay Line setting for PDL37, PDL38, PDL39, PDL40
738 fIndex++;
739 buf[fIndex]=baseWord;
740
741 // Delay Line setting for PDL33, PDL34, PDL35, PDL36
742 fIndex++;
743 buf[fIndex]=baseWord;
744
745 // Delay Line setting for PDL29, PDL30, PDL31, PDL32
746 fIndex++;
747 buf[fIndex]=baseWord;
748
749 // Delay Line setting for PDL25, PDL26, PDL27, PDL28
750 fIndex++;
751 buf[fIndex]=baseWord;
752
753 // Delay Line setting for PDL21, PDL22, PDL23, PDL24
754 fIndex++;
755 buf[fIndex]=baseWord;
756
757 // Delay Line setting for PDL17, PDL18, PDL19, PDL20
758 fIndex++;
759 buf[fIndex]=baseWord;
760
761 // Delay Line setting for PDL13, PDL14, PDL15, PDL16
762 fIndex++;
763 buf[fIndex]=baseWord;
764
765 // Delay Line setting for PDL9, PDL10, PDL11, PDL12
766 fIndex++;
767 buf[fIndex]=baseWord;
768
769 // Delay Line setting for PDL5, PDL6, PDL7, PDL8
770 fIndex++;
771 buf[fIndex]=baseWord;
772
773 // Delay Line setting for PDL1, PDL2, PDL3, PDL4
774 fIndex++;
775 buf[fIndex]=baseWord;
776
777}
778
779//----------------------------------------------------------------------------
780void AliTOFDDLRawData::MakeLTMtrailer(UInt_t *buf)
781{
782 //
783 // LTM trailer
784 //
785
786 UInt_t baseWord;
787 UInt_t word;
788
789 baseWord=0;
790 word = 2; // 0010 -> LTM data are coming from the VME slot number 2
791 AliBitPacking::PackWord(word,baseWord, 0, 3);
792 word = 0; // event crc
793 AliBitPacking::PackWord(word,baseWord, 4,15);
794 word = 0; // event number
795 AliBitPacking::PackWord(word,baseWord,16,27);
796 word = 5; // 0101 -> LTM trailer ID
797 AliBitPacking::PackWord(word,baseWord,28,31);
798 fIndex++;
799 buf[fIndex]=baseWord;
800
801}
802
803//----------------------------------------------------------------------------
804void AliTOFDDLRawData::MakeTRMheader(Int_t nTRM, UInt_t *buf)
805{
806 //
807 // TRM header for the TRM number nTRM [ 3;12]
808 //
809
810 if (nTRM<3 || nTRM>12) {
811 AliWarning(Form(" TRM number is out of the right range [3;12] (nTRM = %3i",nTRM));
812 return;
813 }
814
815 UInt_t baseWord;
816 UInt_t word;
817
818 baseWord = 0;
819 word = nTRM; // TRM data coming from the VME slot number nTRM
820 AliBitPacking::PackWord(word,baseWord, 0, 3);
821 word = 0; // event words
822 AliBitPacking::PackWord(word,baseWord, 4,16);
15ec34b9 823
824 if (fPackedAcquisition)
825 word = 0; // ACQuisition mode: [0;3] see document
826 else
827 word = 3; // ACQuisition mode: [0;3] see document
d0eb8f39 828 AliBitPacking::PackWord(word,baseWord,17,18);
829 word = 0; // description of a SEU inside LUT tables for INL compensation;
830 // the data are unaffected
831 AliBitPacking::PackWord(word,baseWord,19,19);
832 word = 0; // Must Be Zero (MBZ)
833 AliBitPacking::PackWord(word,baseWord,20,27);
834 word = 4; // 0100 -> TRM header ID
835 AliBitPacking::PackWord(word,baseWord,28,31);
836 fIndex++;
837 buf[fIndex]=baseWord;
838
839}
840
841//----------------------------------------------------------------------------
842void AliTOFDDLRawData::MakeTRMtrailer(UInt_t *buf)
843{
844 //
845 // TRM trailer
846 //
847
848 UInt_t baseWord;
849 UInt_t word;
850
851 baseWord=0;
852 word = 15; // 1111 -> TRM trailer ID 1
853 AliBitPacking::PackWord(word,baseWord, 0, 3);
854 word = 0; // event CRC
855 AliBitPacking::PackWord(word,baseWord, 4,15);
856 word = 0; // local event counter == DRM local event counter
857 AliBitPacking::PackWord(word,baseWord,16,27);
858 word = 5; // 0101 -> TRM trailer ID 2
859 AliBitPacking::PackWord(word,baseWord,28,31);
860 fIndex++;
861 buf[fIndex]=baseWord;
862
863}
864
865//----------------------------------------------------------------------------
866void AliTOFDDLRawData::MakeTRMchainHeader(Int_t nTRM, Int_t iChain,
867 UInt_t *buf)
868{
869 //
870 // TRM chain header
871 //
872
873 UInt_t baseWord;
874 UInt_t word;
875
876 if (nTRM<3 || nTRM>12) {
877 AliWarning(Form(" TRM number is out of the right range [3;12] (nTRM = %3i", nTRM));
878 return;
879 }
880
881 if (iChain<0 || iChain>1) {
882 AliWarning(Form(" Chain number is out of the right range [0;1] (iChain = %3i", iChain));
883 return;
884 }
885
886 baseWord=0;
887 word = nTRM; // TRM data coming from the VME slot ID nTRM
888 AliBitPacking::PackWord(word,baseWord, 0, 3);
889 word = 0; // bunch ID
890 AliBitPacking::PackWord(word,baseWord, 4,15);
891 word = 100; // PB24 temperature -> 4 X 25 degree (environment temperature)
892 AliBitPacking::PackWord(word,baseWord,16,23);
893 word = (Int_t)(5 * gRandom->Rndm()); // PB24 ID [0;4]
894 AliBitPacking::PackWord(word,baseWord,24,26);
895 word = 0; // TS
896 AliBitPacking::PackWord(word,baseWord,27,27);
897 switch (iChain) {
898 case 0:
899 word = 0; // 0000 -> TRM chain 0 ID
900 break;
901 case 1:
902 word = 2; // 0010 -> TRM chain 1 ID
903 break;
904 }
905 AliBitPacking::PackWord(word,baseWord,28,31);
906 fIndex++;
907 buf[fIndex]=baseWord;
908
909}
910
911//----------------------------------------------------------------------------
912void AliTOFDDLRawData::MakeTRMfiller(UInt_t *buf, UInt_t nWordsPerTRM)
913{
914 //
915 // TRM filler
916 //
917
918 Int_t jj = -1;
919
920 fIndex++;
921 for (jj=fIndex; jj>fIndex-(Int_t)nWordsPerTRM; jj--) {
922 buf[jj] = buf[jj-1];
923 }
924
925 buf[fIndex-nWordsPerTRM] = MakeFiller();
926
927}
928
929//----------------------------------------------------------------------------
3407e19b 930UInt_t AliTOFDDLRawData::MakeFiller() const
d0eb8f39 931{
932 //
933 // Filler word definition: to make even the number of words per TRM/LTM
934 //
935
936 UInt_t baseWord;
937 UInt_t word;
938
939 baseWord=0;
940 word = 0; // 0000 -> filler ID 1
941 AliBitPacking::PackWord(word,baseWord, 0, 3);
942 word = 0; // MBZ
943 AliBitPacking::PackWord(word,baseWord, 4,27);
944 word = 7; // 0111 -> filler ID 2
945 AliBitPacking::PackWord(word,baseWord, 28,31);
946
947 return baseWord;
948
949}
950
951//----------------------------------------------------------------------------
952void AliTOFDDLRawData::MakeTRMchainTrailer(Int_t iChain, UInt_t *buf)
953{
954 //
955 // TRM chain trailer
956 //
957
958 if (iChain<0 || iChain>1) {
959 AliWarning(Form(" Chain number is out of the right range [0;1] (iChain = %3i", iChain));
960 return;
961 }
962
963 UInt_t baseWord;
964 UInt_t word;
965
966 baseWord=0;
967 word = 0; // status
968 AliBitPacking::PackWord(word,baseWord, 0, 3);
969 word = 0; // MBZ
970 AliBitPacking::PackWord(word,baseWord, 4,15);
971 word = 0; // event counter
972 AliBitPacking::PackWord(word,baseWord,16,27);
973 switch (iChain) {
974 case 0:
975 word = 1; // 0001 -> TRM chain 0 trailer ID
976 break;
977 case 1:
978 word = 3; // 0011 -> TRM chain 1 trailer ID
979 break;
980 }
981 AliBitPacking::PackWord(word,baseWord,28,31);
982 fIndex++;
983 buf[fIndex]=baseWord;
984
985}
986
987//----------------------------------------------------------------------------
988void AliTOFDDLRawData::MakeTDCdigits(Int_t nDDL, Int_t nTRM, Int_t iChain,
989 UInt_t *buf, UInt_t &nWordsPerTRM)
990{
991 //
992 // TRM TDC digit
993 //
994
96f01799 995 const Double_t kOneMoreFilledCell = 1./(AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors());
996 Double_t percentFilledCells = Double_t(fTOFdigitMap->GetFilledCellNumber())/(AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors());
15ec34b9 997
d0eb8f39 998 if (nDDL<0 || nDDL>71) {
999 AliWarning(Form(" DDL number is out of the right range [0;71] (nDDL = %3i", nDDL));
1000 return;
1001 }
1002
1003 if (nTRM<3 || nTRM>12) {
1004 AliWarning(Form(" TRM number is out of the right range [3;12] (nTRM = %3i", nTRM));
1005 return;
1006 }
1007
1008 if (iChain<0 || iChain>1) {
1009 AliWarning(Form(" Chain number is out of the right range [0;1] (iChain = %3i", iChain));
1010 return;
1011 }
1012
15ec34b9 1013 Int_t psArray[1000];
1014 UInt_t localBuffer[1000];
1015 Int_t localIndex = -1;
1016
d0eb8f39 1017 Int_t iDDL = nDDL%AliTOFGeometry::NDDL();
1018
1019 Int_t volume[5] = {-1, -1, -1, -1, -1};
1020 Int_t indexDigit[3] = {-1, -1, -1};
1021
1022 Int_t totCharge = -1;
1023 Int_t timeOfFlight = -1;
1024
15ec34b9 1025 Int_t trailingSpurious = -1;
1026 Int_t leadingSpurious = -1;
1027
d0eb8f39 1028 AliTOFdigit *digs;
1029
1030 UInt_t baseWord=0;
1031 UInt_t word=0;
1032
1033 Int_t jj = -1;
1034 Int_t nTDC = -1;
1035 Int_t iCH = -1;
1036
a403781d 1037 //Int_t numberOfMeasuresPerChannel = 0;
1038 //Int_t maxMeasuresPerChannelInTDC = 0;
515faf5d 1039
1040 Bool_t outOut = HeadOrTail();
1041
d0eb8f39 1042 ofstream ftxt;
1043
1044 if (fVerbose==2) ftxt.open("TOFdigits.txt",ios::app);
1045
1046 for (jj=0; jj<5; jj++) volume[jj] = -1;
1047
1048 // loop on TDC number
1049 for (nTDC=AliTOFGeometry::NTdc()-1; nTDC>=0; nTDC--) {
1050
1051 // the DRM odd (i.e. left) slot number 3 doesn't contain TDC digit data
1052 // for TDC numbers 3-14
96b13296 1053 if (iDDL%2==1 && nTRM==3 && (Int_t)(nTDC/3)!=0) continue;
d0eb8f39 1054
1055 // loop on TDC channel number
1056 for (iCH=AliTOFGeometry::NCh()-1; iCH>=0; iCH--) {
1057
a403781d 1058 //numberOfMeasuresPerChannel = 0;
515faf5d 1059
96b13296 1060 for (Int_t aa=0; aa<5; aa++) volume[aa]=-1;
177b002e 1061 AliTOFRawStream::EquipmentId2VolumeId(nDDL, nTRM, iChain, nTDC, iCH, volume);
d0eb8f39 1062
96b13296 1063 AliDebug(3,Form(" volume -> %2d %1d %2d %2d %1d",volume[0],volume[1],volume[2],volume[3],volume[4]));
1064
d0eb8f39 1065 if (volume[0]==-1 || volume[1]==-1 || volume[2]==-1 ||
1066 volume[3]==-1 || volume[4]==-1) continue;
d0eb8f39 1067
96b13296 1068 AliDebug(3,Form(" ====== %2d %1d %2d %2d %1d",volume[0],volume[1],volume[2],volume[3],volume[4]));
1069
d0eb8f39 1070 for (jj=0; jj<3; jj++) indexDigit[jj] = -1;
1071
1072 fTOFdigitMap->GetDigitIndex(volume, indexDigit);
1073
15ec34b9 1074 if (indexDigit[0]<0) {
1075
515faf5d 1076 trailingSpurious = Int_t(2097152*gRandom->Rndm());
1077 leadingSpurious = Int_t(2097152*gRandom->Rndm());
15ec34b9 1078
515faf5d 1079 if ( fFakeOrphaneProduction &&
1080 ( ( fPackedAcquisition && percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) ) ||
1081 (!fPackedAcquisition && percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) ) ) ) {
15ec34b9 1082
1083 percentFilledCells+=kOneMoreFilledCell;
1084
1085 Int_t dummyPS = 0;
1086
515faf5d 1087 if (outOut) {
15ec34b9 1088 word = trailingSpurious; // trailing edge measurement
1089 dummyPS = 2;
1090 }
1091 else {
1092 word = leadingSpurious; // leading edge measurement
1093 dummyPS = 1;
1094 }
1095
1096 if (fVerbose==2) {
1097 if (nDDL<10) ftxt << " " << nDDL;
1098 else ftxt << " " << nDDL;
1099 if (nTRM<10) ftxt << " " << nTRM;
1100 else ftxt << " " << nTRM;
1101 ftxt << " " << iChain;
1102 if (nTDC<10) ftxt << " " << nTDC;
1103 else ftxt << " " << nTDC;
1104 ftxt << " " << iCH;
1105 if (volume[0]<10) ftxt << " -> " << volume[0];
1106 else ftxt << " -> " << volume[0];
1107 ftxt << " " << volume[1];
1108 if (volume[2]<10) ftxt << " " << volume[2];
1109 else ftxt << " " << volume[2];
1110 ftxt << " " << volume[4];
1111 if (volume[3]<10) ftxt << " " << volume[3];
1112 else ftxt << " " << volume[3];
515faf5d 1113 ftxt << " " << -1;
1114 if (word<10) ftxt << " " << word;
1115 else if (word>=10 && word<100) ftxt << " " << word;
1116 else if (word>=100 && word<1000) ftxt << " " << word;
1117 else if (word>=1000 && word<10000) ftxt << " " << word;
1118 else if (word>=10000 && word<100000) ftxt << " " << word;
1119 else if (word>=100000 && word<1000000) ftxt << " " << word;
1120 else ftxt << " " << word;
1121 ftxt << " " << dummyPS << endl;
15ec34b9 1122 }
1123
1124 AliBitPacking::PackWord(word,baseWord, 0,20);
1125 word = iCH; // TDC channel ID [0;7]
1126 AliBitPacking::PackWord(word,baseWord,21,23);
1127 word = nTDC; // TDC ID [0;14]
1128 AliBitPacking::PackWord(word,baseWord,24,27);
1129 word = 0; // error flag
1130 AliBitPacking::PackWord(word,baseWord,28,28);
1131 word = dummyPS; // Packing Status [0;3]
1132 AliBitPacking::PackWord(word,baseWord,29,30);
1133 word = 1; // TRM TDC digit ID
1134 AliBitPacking::PackWord(word,baseWord,31,31);
1135
1136 localIndex++;
1137 localBuffer[localIndex]=baseWord;
1138 psArray[localIndex]=dummyPS;
1139
1140 nWordsPerTRM++;
1141 baseWord=0;
1142
515faf5d 1143 } // if ( fFakeOrphaneProduction && ( ( fPackedAcquisition && percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) ) or ... ) )
15ec34b9 1144 } // if (indexDigit[0]<0)
1145
d0eb8f39 1146 for (jj=0; jj<3;jj++) {
1147
1148 if (indexDigit[jj]<0) continue;
a403781d 1149
96b13296 1150 AliDebug(3,Form(" ====== %2d %1d %2d %2d %1d -> %1d %d",volume[0],volume[1],volume[2],volume[3],volume[4],jj,indexDigit[jj]));
1151
d0eb8f39 1152 digs = (AliTOFdigit*)fTOFdigitArray->UncheckedAt(indexDigit[jj]);
1153
1154 if (digs->GetSector()!=volume[0] ||
1155 digs->GetPlate() !=volume[1] ||
1156 digs->GetStrip() !=volume[2] ||
1157 digs->GetPadx() !=volume[3] ||
6e4ef37b 1158 digs->GetPadz() !=volume[4]) AliWarning(Form(" --- ERROR --- %2i (%2i) %1i (%1i) %2i (%2i) %2i (%2i) %1i (%1i)",
1159 digs->GetSector(), volume[0],
1160 digs->GetPlate(), volume[1],
1161 digs->GetStrip(), volume[2],
1162 digs->GetPadx(), volume[3],
1163 digs->GetPadz(), volume[4])
1164 );
d0eb8f39 1165
96b13296 1166 timeOfFlight = (Int_t)(digs->GetTdc());
1167 /*+
cfceaf15 1168 fTOFCableLengthMap->GetCableTimeShiftBin(nDDL, nTRM, iChain, nTDC)*/;
d11fc181 1169
cfceaf15 1170 //if (timeOfFlight>=fMatchingWindow+1000) continue; // adc
96b13296 1171 //if (timeOfFlight>=fMatchingWindow) continue; // adc
1172 if (digs->GetTdcND()>=fMatchingWindow) {
1173 AliDebug(2,"Out of matching window.");
1174 continue; // adc
1175 }
1176 else AliDebug(2,"In matching window");
515faf5d 1177
a403781d 1178 //numberOfMeasuresPerChannel++;
1179
515faf5d 1180 // totCharge = (Int_t)digs->GetAdc(); //Use realistic ToT, for Standard production with no miscalibration/Slewing it == fAdC in digit (see AliTOFDigitizer)
1181 totCharge = (Int_t)(digs->GetToT());
d0eb8f39 1182 // temporary control
d11fc181 1183 if (totCharge<0) totCharge = 0;//TMath::Abs(totCharge);
d0eb8f39 1184
15ec34b9 1185 if (fPackedAcquisition) {
1186
d0eb8f39 1187 if (fVerbose==2) {
1188 if (nDDL<10) ftxt << " " << nDDL;
1189 else ftxt << " " << nDDL;
1190 if (nTRM<10) ftxt << " " << nTRM;
1191 else ftxt << " " << nTRM;
1192 ftxt << " " << iChain;
1193 if (nTDC<10) ftxt << " " << nTDC;
1194 else ftxt << " " << nTDC;
1195 ftxt << " " << iCH;
1196 if (volume[0]<10) ftxt << " -> " << volume[0];
1197 else ftxt << " -> " << volume[0];
1198 ftxt << " " << volume[1];
1199 if (volume[2]<10) ftxt << " " << volume[2];
1200 else ftxt << " " << volume[2];
1201 ftxt << " " << volume[4];
1202 if (volume[3]<10) ftxt << " " << volume[3];
1203 else ftxt << " " << volume[3];
515faf5d 1204 if (totCharge<10) ftxt << " " << totCharge;
1205 else if (totCharge>=10 && totCharge<100) ftxt << " " << totCharge;
1206 else ftxt << " " << totCharge;
1207 if (timeOfFlight<10) ftxt << " " << timeOfFlight << endl;
1208 else if (timeOfFlight>=10 && timeOfFlight<100) ftxt << " " << timeOfFlight << endl;
1209 else if (timeOfFlight>=100 && timeOfFlight<1000) ftxt << " " << timeOfFlight << endl;
1210 else ftxt << " " << timeOfFlight << endl;
d0eb8f39 1211 }
1212
515faf5d 1213 word = timeOfFlight%8192; // time-of-fligth measurement
d0eb8f39 1214 AliBitPacking::PackWord(word,baseWord, 0,12);
1215
515faf5d 1216 if (totCharge>=256) totCharge = 255;
d0eb8f39 1217 word = totCharge; // time-over-threshould measurement
1218 AliBitPacking::PackWord(word,baseWord,13,20);
1219
1220 word = iCH; // TDC channel ID [0;7]
1221 AliBitPacking::PackWord(word,baseWord,21,23);
1222 word = nTDC; // TDC ID [0;14]
1223 AliBitPacking::PackWord(word,baseWord,24,27);
1224 word = 0; // error flag
1225 AliBitPacking::PackWord(word,baseWord,28,28);
15ec34b9 1226 word = 0; // Packing Status [0;3]
1227 AliBitPacking::PackWord(word,baseWord,29,30);
1228 word = 1; // TRM TDC digit ID
1229 AliBitPacking::PackWord(word,baseWord,31,31);
1230
1231 localIndex++;
1232 localBuffer[localIndex]=baseWord;
1233
1234 nWordsPerTRM++;
1235 baseWord=0;
1236
515faf5d 1237 if ( fFakeOrphaneProduction &&
1238 percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) ) {
15ec34b9 1239
1240 percentFilledCells+=kOneMoreFilledCell;
1241
515faf5d 1242 trailingSpurious = Int_t(2097152*gRandom->Rndm());
1243 leadingSpurious = Int_t(2097152*gRandom->Rndm());
15ec34b9 1244
1245 Int_t dummyPS = 0;
1246
515faf5d 1247 if (outOut) {
15ec34b9 1248 word = trailingSpurious; // trailing edge measurement
1249 dummyPS = 2;
1250 }
1251 else {
1252 word = leadingSpurious; // leading edge measurement
1253 dummyPS = 1;
1254 }
1255
1256 if (fVerbose==2) {
1257 if (nDDL<10) ftxt << " " << nDDL;
1258 else ftxt << " " << nDDL;
1259 if (nTRM<10) ftxt << " " << nTRM;
1260 else ftxt << " " << nTRM;
1261 ftxt << " " << iChain;
1262 if (nTDC<10) ftxt << " " << nTDC;
1263 else ftxt << " " << nTDC;
1264 ftxt << " " << iCH;
1265 if (volume[0]<10) ftxt << " -> " << volume[0];
1266 else ftxt << " -> " << volume[0];
1267 ftxt << " " << volume[1];
1268 if (volume[2]<10) ftxt << " " << volume[2];
1269 else ftxt << " " << volume[2];
1270 ftxt << " " << volume[4];
1271 if (volume[3]<10) ftxt << " " << volume[3];
1272 else ftxt << " " << volume[3];
515faf5d 1273 ftxt << " " << -1;
1274 if (word<10) ftxt << " " << word;
1275 else if (word>=10 && word<100) ftxt << " " << word;
1276 else if (word>=100 && word<1000) ftxt << " " << word;
1277 else if (word>=1000 && word<10000) ftxt << " " << word;
1278 else if (word>=10000 && word<100000) ftxt << " " << word;
1279 else if (word>=100000 && word<1000000) ftxt << " " << word;
1280 else ftxt << " " << word;
1281 ftxt << " " << dummyPS << endl;
15ec34b9 1282 }
1283
1284 AliBitPacking::PackWord(word,baseWord, 0,20);
1285 word = iCH; // TDC channel ID [0;7]
1286 AliBitPacking::PackWord(word,baseWord,21,23);
1287 word = nTDC; // TDC ID [0;14]
1288 AliBitPacking::PackWord(word,baseWord,24,27);
1289 word = 0; // error flag
1290 AliBitPacking::PackWord(word,baseWord,28,28);
1291 word = dummyPS; // Packing Status [0;3]
1292 AliBitPacking::PackWord(word,baseWord,29,30);
1293 word = 1; // TRM TDC digit ID
1294 AliBitPacking::PackWord(word,baseWord,31,31);
1295
1296 localIndex++;
1297 localBuffer[localIndex]=baseWord;
1298 psArray[localIndex]=dummyPS;
1299
1300 nWordsPerTRM++;
1301 baseWord=0;
1302
515faf5d 1303 } // if ( fFakeOrphaneProduction && percentFilledCells<0.12 && gRandom->Rndm()<(0.12-percentFilledCells) )
15ec34b9 1304
1305
1306 } // if (fPackedAcquisition)
1307 else { // if (!fPackedAcquisition)
1308
515faf5d 1309 if ( fFakeOrphaneProduction &&
1310 percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && outOut ) {
15ec34b9 1311
1312 percentFilledCells+=kOneMoreFilledCell;
1313
515faf5d 1314 trailingSpurious = Int_t(2097152*gRandom->Rndm());
15ec34b9 1315 word = trailingSpurious;
1316 Int_t dummyPS = 2;
1317
1318 if (fVerbose==2) {
1319 if (nDDL<10) ftxt << " " << nDDL;
1320 else ftxt << " " << nDDL;
1321 if (nTRM<10) ftxt << " " << nTRM;
1322 else ftxt << " " << nTRM;
1323 ftxt << " " << iChain;
1324 if (nTDC<10) ftxt << " " << nTDC;
1325 else ftxt << " " << nTDC;
1326 ftxt << " " << iCH;
1327 if (volume[0]<10) ftxt << " -> " << volume[0];
1328 else ftxt << " -> " << volume[0];
1329 ftxt << " " << volume[1];
1330 if (volume[2]<10) ftxt << " " << volume[2];
1331 else ftxt << " " << volume[2];
1332 ftxt << " " << volume[4];
1333 if (volume[3]<10) ftxt << " " << volume[3];
1334 else ftxt << " " << volume[3];
515faf5d 1335 ftxt << " " << -1;
1336 if (word<10) ftxt << " " << word;
1337 else if (word>=10 && word<100) ftxt << " " << word;
1338 else if (word>=100 && word<1000) ftxt << " " << word;
1339 else if (word>=1000 && word<10000) ftxt << " " << word;
1340 else if (word>=10000 && word<100000) ftxt << " " << word;
1341 else if (word>=100000 && word<1000000) ftxt << " " << word;
1342 else ftxt << " " << word;
1343 ftxt << " " << dummyPS << endl;
15ec34b9 1344 }
1345
1346 AliBitPacking::PackWord(word,baseWord, 0,20);
1347 word = iCH; // TDC channel ID [0;7]
1348 AliBitPacking::PackWord(word,baseWord,21,23);
1349 word = nTDC; // TDC ID [0;14]
1350 AliBitPacking::PackWord(word,baseWord,24,27);
1351 word = 0; // error flag
1352 AliBitPacking::PackWord(word,baseWord,28,28);
1353 word = dummyPS; // Packing Status [0;3]
1354 AliBitPacking::PackWord(word,baseWord,29,30);
1355 word = 1; // TRM TDC digit ID
1356 AliBitPacking::PackWord(word,baseWord,31,31);
1357
1358 localIndex++;
1359 localBuffer[localIndex]=baseWord;
1360 psArray[localIndex]=dummyPS;
1361
1362 nWordsPerTRM++;
1363 baseWord=0;
1364
515faf5d 1365 } // if ( fFakeOrphaneProduction && percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && outOut )
15ec34b9 1366
1367
515faf5d 1368 word = (timeOfFlight + Int_t(totCharge*AliTOFGeometry::ToTBinWidth()/AliTOFGeometry::TdcBinWidth()))%2097152; // trailing edge measurement
15ec34b9 1369
1370 if (fVerbose==2) {
1371 if (nDDL<10) ftxt << " " << nDDL;
1372 else ftxt << " " << nDDL;
1373 if (nTRM<10) ftxt << " " << nTRM;
1374 else ftxt << " " << nTRM;
1375 ftxt << " " << iChain;
1376 if (nTDC<10) ftxt << " " << nTDC;
1377 else ftxt << " " << nTDC;
1378 ftxt << " " << iCH;
1379 if (volume[0]<10) ftxt << " -> " << volume[0];
1380 else ftxt << " -> " << volume[0];
1381 ftxt << " " << volume[1];
1382 if (volume[2]<10) ftxt << " " << volume[2];
1383 else ftxt << " " << volume[2];
1384 ftxt << " " << volume[4];
1385 if (volume[3]<10) ftxt << " " << volume[3];
1386 else ftxt << " " << volume[3];
515faf5d 1387 ftxt << " " << -1;
1388 if (word<10) ftxt << " " << word;
1389 else if (word>=10 && word<100) ftxt << " " << word;
1390 else if (word>=100 && word<1000) ftxt << " " << word;
1391 else if (word>=1000 && word<10000) ftxt << " " << word;
1392 else if (word>=10000 && word<100000) ftxt << " " << word;
1393 else if (word>=100000 && word<1000000) ftxt << " " << word;
1394 else ftxt << " " << word;
1395 ftxt << " " << 2 << endl;
15ec34b9 1396 }
1397
1398 AliBitPacking::PackWord(word,baseWord, 0,20);
1399
1400 word = iCH; // TDC channel ID [0;7]
1401 AliBitPacking::PackWord(word,baseWord,21,23);
1402 word = nTDC; // TDC ID [0;14]
1403 AliBitPacking::PackWord(word,baseWord,24,27);
1404 word = 0; // error flag
1405 AliBitPacking::PackWord(word,baseWord,28,28);
1406 word = 2; // Packing Status [0;3]
d0eb8f39 1407 AliBitPacking::PackWord(word,baseWord,29,30);
1408 word = 1; // TRM TDC digit ID
1409 AliBitPacking::PackWord(word,baseWord,31,31);
15ec34b9 1410
1411 localIndex++;
1412 localBuffer[localIndex]=baseWord;
1413 psArray[localIndex]=2;
d0eb8f39 1414
1415 nWordsPerTRM++;
1416 baseWord=0;
1417
515faf5d 1418 word = timeOfFlight%2097152; // leading edge measurement
15ec34b9 1419
1420 if (fVerbose==2) {
1421 if (nDDL<10) ftxt << " " << nDDL;
1422 else ftxt << " " << nDDL;
1423 if (nTRM<10) ftxt << " " << nTRM;
1424 else ftxt << " " << nTRM;
1425 ftxt << " " << iChain;
1426 if (nTDC<10) ftxt << " " << nTDC;
1427 else ftxt << " " << nTDC;
1428 ftxt << " " << iCH;
1429 if (volume[0]<10) ftxt << " -> " << volume[0];
1430 else ftxt << " -> " << volume[0];
1431 ftxt << " " << volume[1];
1432 if (volume[2]<10) ftxt << " " << volume[2];
1433 else ftxt << " " << volume[2];
1434 ftxt << " " << volume[4];
1435 if (volume[3]<10) ftxt << " " << volume[3];
1436 else ftxt << " " << volume[3];
515faf5d 1437 ftxt << " " << -1;
1438 if (word<10) ftxt << " " << word;
1439 else if (word>=10 && word<100) ftxt << " " << word;
1440 else if (word>=100 && word<1000) ftxt << " " << word;
1441 else if (word>=1000 && word<10000) ftxt << " " << word;
1442 else if (word>=10000 && word<100000) ftxt << " " << word;
1443 else if (word>=100000 && word<1000000) ftxt << " " << word;
1444 else ftxt << " " << word;
1445 ftxt << " " << 1 << endl;
15ec34b9 1446 }
1447
1448 AliBitPacking::PackWord(word,baseWord, 0,20);
1449
1450 word = iCH; // TDC channel ID [0;7]
1451 AliBitPacking::PackWord(word,baseWord,21,23);
1452 word = nTDC; // TDC ID [0;14]
1453 AliBitPacking::PackWord(word,baseWord,24,27);
1454 word = 0; // error flag
1455 AliBitPacking::PackWord(word,baseWord,28,28);
1456 word = 1; // Packing Status [0;3]
1457 AliBitPacking::PackWord(word,baseWord,29,30);
1458 word = 1; // TRM TDC digit ID
1459 AliBitPacking::PackWord(word,baseWord,31,31);
1460
1461 localIndex++;
1462 localBuffer[localIndex]=baseWord;
1463 psArray[localIndex]=1;
1464
1465 nWordsPerTRM++;
1466 baseWord=0;
1467
1468
515faf5d 1469 if ( fFakeOrphaneProduction &&
1470 percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && !outOut ) {
15ec34b9 1471
1472 percentFilledCells+=kOneMoreFilledCell;
1473
515faf5d 1474 leadingSpurious = Int_t(2097152*gRandom->Rndm());
15ec34b9 1475 word = leadingSpurious;
1476 Int_t dummyPS = 1;
1477
1478 if (fVerbose==2) {
1479 if (nDDL<10) ftxt << " " << nDDL;
1480 else ftxt << " " << nDDL;
1481 if (nTRM<10) ftxt << " " << nTRM;
1482 else ftxt << " " << nTRM;
1483 ftxt << " " << iChain;
1484 if (nTDC<10) ftxt << " " << nTDC;
1485 else ftxt << " " << nTDC;
1486 ftxt << " " << iCH;
1487 if (volume[0]<10) ftxt << " -> " << volume[0];
1488 else ftxt << " -> " << volume[0];
1489 ftxt << " " << volume[1];
1490 if (volume[2]<10) ftxt << " " << volume[2];
1491 else ftxt << " " << volume[2];
1492 ftxt << " " << volume[4];
1493 if (volume[3]<10) ftxt << " " << volume[3];
1494 else ftxt << " " << volume[3];
515faf5d 1495 ftxt << " " << -1;
1496 if (word<10) ftxt << " " << word;
1497 else if (word>=10 && word<100) ftxt << " " << word;
1498 else if (word>=100 && word<1000) ftxt << " " << word;
1499 else if (word>=1000 && word<10000) ftxt << " " << word;
1500 else if (word>=10000 && word<100000) ftxt << " " << word;
1501 else if (word>=100000 && word<1000000) ftxt << " " << word;
1502 else ftxt << " " << word;
1503 ftxt << " " << dummyPS << endl;
15ec34b9 1504 }
1505
1506 AliBitPacking::PackWord(word,baseWord, 0,20);
1507 word = iCH; // TDC channel ID [0;7]
1508 AliBitPacking::PackWord(word,baseWord,21,23);
1509 word = nTDC; // TDC ID [0;14]
1510 AliBitPacking::PackWord(word,baseWord,24,27);
1511 word = 0; // error flag
1512 AliBitPacking::PackWord(word,baseWord,28,28);
1513 word = dummyPS; // Packing Status [0;3]
1514 AliBitPacking::PackWord(word,baseWord,29,30);
1515 word = 1; // TRM TDC digit ID
1516 AliBitPacking::PackWord(word,baseWord,31,31);
1517
1518 localIndex++;
1519 localBuffer[localIndex]=baseWord;
1520 psArray[localIndex]=dummyPS;
1521
1522 nWordsPerTRM++;
1523 baseWord=0;
1524
515faf5d 1525 } // if ( fFakeOrphaneProduction && percentFilledCells<0.24 && gRandom->Rndm()<(0.24-percentFilledCells) && !outOut )
15ec34b9 1526
1527
1528 } // if (!fPackedAcquisition)
1529
d0eb8f39 1530 } //end loop on digits in the same volume
1531
a403781d 1532 //if (numberOfMeasuresPerChannel>maxMeasuresPerChannelInTDC)
1533 //maxMeasuresPerChannelInTDC = numberOfMeasuresPerChannel;
1534
d0eb8f39 1535 } // end loop on TDC channel number
1536
e8285af5 1537 //AliDebug(3,Form(" TDC number %2i: numberOfMeasuresPerChannel = %2i --- maxMeasuresPerChannelInTDC = %2i ", nTDC, numberOfMeasuresPerChannel, maxMeasuresPerChannelInTDC));
a403781d 1538
515faf5d 1539 if (localIndex==-1) continue;
15ec34b9 1540
515faf5d 1541 if (fPackedAcquisition) {
15ec34b9 1542
3c609b5c 1543 for (jj=0; jj<=localIndex; jj++) {
15ec34b9 1544 fIndex++;
1545 buf[fIndex] = localBuffer[jj];
515faf5d 1546 localBuffer[jj] = 0;
1547 psArray[jj] = -1;
15ec34b9 1548 }
515faf5d 1549
15ec34b9 1550 }
515faf5d 1551 else {
a403781d 1552 /*
515faf5d 1553 if (maxMeasuresPerChannelInTDC = 1) {
1554
1555 for (Int_t jj=0; jj<=localIndex; jj++) {
1556 if (psArray[jj]==2) {
1557 fIndex++;
1558 buf[fIndex] = localBuffer[jj];
1559 localBuffer[jj] = 0;
1560 psArray[jj] = -1;
1561 }
1562 }
1563 for (Int_t jj=0; jj<=localIndex; jj++) {
1564 if (psArray[jj]==1) {
1565 fIndex++;
1566 buf[fIndex] = localBuffer[jj];
1567 localBuffer[jj] = 0;
1568 psArray[jj] = -1;
1569 }
1570 }
1571
1572 } // if (maxMeasuresPerChannelInTDC = 1)
1573 else if (maxMeasuresPerChannelInTDC>1) {
1574
96b13296 1575 AliDebug(3,Form(" In the TOF DDL %2i, TRM %2i, TDC %2i, chain %1i, the maximum number of t.o.f. good measurements per channel is %2i",
515faf5d 1576 nDDL, nTRM, iChain, nTDC, iCH, maxMeasuresPerChannelInTDC));
a403781d 1577 */
3c609b5c 1578 for (jj=0; jj<=localIndex; jj++) {
515faf5d 1579 fIndex++;
1580 buf[fIndex] = localBuffer[jj];
1581 localBuffer[jj] = 0;
1582 psArray[jj] = -1;
1583 }
1584
a403781d 1585 //} // else if (maxMeasuresPerChannelInTDC>1)
515faf5d 1586
1587 } // else (!fPackedAcquisition)
1588
1589 localIndex = -1;
1590
a403781d 1591 //maxMeasuresPerChannelInTDC = 0;
1592
515faf5d 1593 } // end loop on TDC number
15ec34b9 1594
d0eb8f39 1595
1596 if (fVerbose==2) ftxt.close();
1597
1598}
1599
1600//----------------------------------------------------------------------------
1601void AliTOFDDLRawData::ReverseArray(UInt_t a[], Int_t n) const
1602{
1603 //
1604 // Reverses the n elements of array a
1605 //
1606
1607 Int_t ii, temp;
1608
1609 for (ii=0; ii<n/2; ii++) {
1610 temp = a[ii];
1611 a[ii] = a[n-ii-1];
1612 a[n-ii-1] = temp;
1613 }
1614
1615 return;
1616
1617}
15ec34b9 1618
1619//----------------------------------------------------------------------------
1620Bool_t AliTOFDDLRawData::HeadOrTail() const
1621{
1622 //
1623 // Returns the result of a 'pitch and toss'
1624 //
1625
1626 Double_t dummy = gRandom->Rndm();
1627
1628 if (dummy<0.5) return kFALSE;
1629 else return kTRUE;
1630
1631}