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