1 /**************************************************************************
2 * Copyright(c) 1998-2003, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
19 //This class contains all the necessary methods to create the Raw Data
20 //files (slides) for the ITS data challenges for:
26 //#include <Riostream.h>
27 #include <TClonesArray.h>
29 #include "AliITSdigit.h"
30 #include "AliITSDDLRawData.h"
31 #include "AliRawDataHeaderSim.h"
32 #include "AliITSRawStreamSPD.h"
33 #include "AliITSRawStreamSDD.h"
34 #include "AliITSDDLModuleMapSDD.h"
35 #include "AliITSRawStreamSSD.h"
36 #include "AliITSIntMap.h"
37 #include "AliBitPacking.h"
39 #include "AliFstream.h"
41 ClassImp(AliITSDDLRawData)
43 ////////////////////////////////////////////////////////////////////////////////////////
44 AliITSDDLRawData::AliITSDDLRawData():
52 ////////////////////////////////////////////////////////////////////////////////////////
54 AliITSDDLRawData::AliITSDDLRawData(const AliITSDDLRawData &source) :
56 fVerbose(source.fVerbose),
57 fIndex(source.fIndex),
58 fHalfStaveModule(source.fHalfStaveModule){
62 ////////////////////////////////////////////////////////////////////////////////////////
64 AliITSDDLRawData& AliITSDDLRawData::operator=(const AliITSDDLRawData &source){
66 this->fIndex=source.fIndex;
67 this->fHalfStaveModule=source.fHalfStaveModule;
68 this->fVerbose=source.fVerbose;
72 ////////////////////////////////////////////////////////////////////////////////////////
76 void AliITSDDLRawData::GetDigitsSSD(TClonesArray *ITSdigits,Int_t mod,Int_t modR,Int_t ddl,UInt_t *buf){
77 //This method packs the SSD digits in a proper 32 bits structure
78 // Revised by Enrico Fragiacomo
84 Int_t ndigits = ITSdigits->GetEntries();
89 ftxt.open("SSDdigits.txt",ios::app);
91 for (Int_t digit=0;digit<ndigits;digit++) {
92 digs = (AliITSdigit*)ITSdigits->UncheckedAt(digit);
93 iz=digs->GetCoord1(); // If iz==0, O side and if iz=1 N side
94 ix=digs->GetCoord2(); // Strip Number
95 is=digs->GetCompressedSignal(); // ADC Signal
96 // cout<<" Module:"<<mod-500<<" N/P side:"<<iz<<" Strip Number:"<<ix<<" Amplidute:"<<is-1<<endl;
97 if(is<0) is = 4096 + is;
99 ftxt<<"DDL:"<<ddl<<" Mod: "<<modR<<" N/P: "<<iz<<" Strip: "<<ix<<" Value: "<<is-1<<endl;
104 AliBitPacking::PackWord(word,baseWord,0,11);//ADC data
106 word = (iz==0) ? ix : 1535-ix ; // on N-side 1535-768 -> 0-767
107 AliBitPacking::PackWord(word,baseWord,12,22);//Strip Number
109 word = mod%12; // ADC-number (12 ADCs per AD module)
110 word += ( word<6 ) ? 0 : 2; // ADC range 0-5 and 8-13
111 AliBitPacking::PackWord(word,baseWord,24,27);//ADC Channel
113 word = mod/12+1; // AD-number (AD module index ranges 1-9)
114 AliBitPacking::PackWord(word,baseWord,28,31);//AD slot
116 buf[fIndex]=baseWord;
124 ////////////////////////////////////////////////////////////////////////////////////////
125 //Silicon Drift Detector
128 void AliITSDDLRawData::GetDigitsSDD(TClonesArray *ITSdigits,Int_t mod,Int_t modR,Int_t ddl,UInt_t *buf){
129 //This method packs the SDD digits in a proper 32 bits structure
135 Int_t ndigits = ITSdigits->GetEntries();
138 Int_t digarr[512][256];
139 for(Int_t i=0;i<512;i++){
140 for(Int_t j=0;j<256;j++){
144 //word to select the 12 carlos for the 12 modules
145 UInt_t carlosid=0x30000000+mod;
148 buf[fIndex]=carlosid;
154 Bool_t flag = kFALSE;
156 Int_t bitinfo1[4] = {3,8,3,7}; //vector with info on bit for timebin info
157 Int_t wordinfo1[4]= {0,0,0,0}; //vector with word info for timebin info
158 Int_t bitinfo2[2] = {3,18}; //vector with info on bit for EOR (end of row) info
159 Int_t wordinfo2[3]= {1,65593}; //vector with word info for anode info
161 /* for time bin info: word n bits meaning
162 0 3 next info is timebin
163 8 3 next word is 8 bit long
164 tb value 8 timebin value
165 n (2->7) 3 next info is n bit long
166 signal n signal value
168 for anode info: 1 3 next 18 bits are for EOR
169 increments the anode value
171 EOR 18 error codes + other info
176 ftxt.open("SDDdigits.txt",ios::app);
177 for (Int_t digit=0;digit<ndigits;digit++) {
178 digs = (AliITSdigit*)ITSdigits->UncheckedAt(digit);
179 iz=digs->GetCoord1(); // Anode
180 ix=digs->GetCoord2(); // Time
181 is=digs->GetCompressedSignal(); // ADC Signal
184 ftxt<<"DDL:"<<ddl<<" MID:"<<modR<<" An:"<<iz<<" T:"<<ix<<" A:"<<is<<endl;
185 if (is>255){Error("GetDigitsSDD", "bits words is needed)!!!");}
188 for(Int_t anode=0;anode<512;anode++){
191 AliBitPacking::PackWord(word2,baseWord,first,last);
204 for(Int_t tb=0;tb<256;tb++){
205 if(digarr[anode][tb]!=0){
208 AliBitPacking::PackWord(word2,baseWord,first,last);
214 //non lossy compression as it is done in Carlos
215 //(data are already 10to8bit compressed by AMBRA
217 /* if value < 8 value = value - (1 << 2) (word is 2 bit long)
218 if value < 16 value = value - (1 << 3) (word is 3 bit long)
219 if value < 32 value = value - (1 << 4) (word is 4 bit long)
220 if value < 64 value = value - (1 << 5) (word is 5 bit long)
221 if value <128 value = value - (1 << 6) (word is 6 bit long)
222 if value >=128value = value - (1 << 7) (word is 7 bit long)
225 if(digarr[anode][tb]<8){
228 wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]);
230 if(digarr[anode][tb]>=8 && digarr[anode][tb]<16){
233 wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]);
235 if(digarr[anode][tb]>=16 && digarr[anode][tb]<32){
238 wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]);
240 if(digarr[anode][tb]>=32 && digarr[anode][tb]<64){
243 wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]);
245 if(digarr[anode][tb]>=64 && digarr[anode][tb]<128){
248 wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]);
250 if(digarr[anode][tb]>=128){
253 wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]);
256 for(Int_t ie=0;ie<4;ie++){
260 AliBitPacking::PackWord(word2,baseWord,first,last);
265 last = first+bitinfo1[ie]-1;
266 if(first < 30 && last < 30){
267 AliBitPacking::PackWord(wordinfo1[ie],baseWord,first,last);
272 UInt_t w = AliBitPacking::UnpackWord(wordinfo1[ie],0,29-first);
273 AliBitPacking::PackWord(w,baseWord,first,29);
274 Int_t lb = 29-first+1;
275 diff = bitinfo1[ie]-lb;
276 word2 = AliBitPacking::UnpackWord(wordinfo1[ie],lb,lb+diff-1);
278 if(anode<256) word = 2;//channel 0 of carlos
279 else word = 3; //channel 1 of carlos
280 AliBitPacking::PackWord(word,baseWord,30,31);
282 buf[fIndex]=baseWord;
289 word2 = wordinfo1[ie];
292 if(anode<256) word = 2; //channel 0 of carlos
293 else word = 3; //channel 1 of carlos
294 AliBitPacking::PackWord(word,baseWord,30,31);
296 buf[fIndex]=baseWord;
309 for(Int_t i=0;i<2;i++){
312 AliBitPacking::PackWord(word2,baseWord,first,last);
321 if(first < 30 && last < 30){
322 AliBitPacking::PackWord(word,baseWord,first,last); //3 bit code =1 -> next 18 bits for EOR
328 UInt_t w = AliBitPacking::UnpackWord(word,0,29-first);
329 AliBitPacking::PackWord(w,baseWord,first,29);
330 Int_t lb = 29-first+1;
332 word2 = AliBitPacking::UnpackWord(word,lb,lb+diff-1);
334 if(anode<256) word = 2;
336 AliBitPacking::PackWord(word,baseWord,30,31);
338 buf[fIndex]=baseWord;
351 if(anode<256) word = 2;
353 AliBitPacking::PackWord(word,baseWord,30,31);
355 buf[fIndex]=baseWord;
374 ////////////////////////////////////////////////////////////////////////////////////////
378 void AliITSDDLRawData::GetDigitsSPD(TClonesArray *ITSdigits,Int_t mod,Int_t ddl, UInt_t *buf){
379 //This method packs the SPD digits in a proper 32 structure
380 //Since data is zero suppressed,the coordinates for the chip having zero digits
381 //doesn't get listed in the galice.root file. However the SPD format requires
382 //the empty chip to be written with chip header and chip trailer.
384 Int_t chipLow = AliITSRawStreamSPD::GetOnlineChipFromOffline(mod,0);
385 Int_t chipHigh = AliITSRawStreamSPD::GetOnlineChipFromOffline(mod,159);
386 if (chipLow>chipHigh) {chipLow -= 4; chipHigh += 4;}
387 UInt_t hs = AliITSRawStreamSPD::GetOnlineHSFromOffline(mod);
389 // create int map to later hold all digits sorted
390 AliITSIntMap* digMap = new AliITSIntMap();
393 Int_t chipHitCount=0; //Number of Hit in the current chip
394 Int_t previousChip=-1; //Previuos chip respect to the actual aone
395 Int_t ndigits = ITSdigits->GetEntries(); //number of digits in the current module
396 //cout<<" Number of digits in the current module:"<<ndigits<<" module:"<<mod<<endl;
403 ftxt.open("SPDdigits.txt",ios::app);
404 for (Int_t digit=0;digit<ndigits;digit++){
405 digs = (AliITSdigit*)ITSdigits->UncheckedAt(digit);
406 /*---------------------------------------------------------------------------
407 * Each module contains 5 read out chips of 256 rows and 32 columns.
408 * So, the cell number in Z direction varies from 0 to 159.
409 * ---------------------------------------------------------------------*/
410 Int_t iz=digs->GetCoord1(); // Cell number in Z direction
411 Int_t ix=digs->GetCoord2(); // Cell number in X direction
414 ftxt<<"DDL:"<<ddl<<" Mod:"<<mod<<" Row:"<<ix<<" Col:"<<iz<<endl;
415 UInt_t dummyDDL,dummyHS,chip,col,row;
416 AliITSRawStreamSPD::OfflineToOnline(mod,iz,ix,dummyDDL,dummyHS,chip,col,row);
418 // insert digit into map...
419 // (reverse order of cols and rows as in real raw data)
420 digMap->Insert(chip*256*32+(32-col)*256+(256-row),row);
424 UInt_t nrHits = digMap->GetNrEntries();
427 for (UInt_t nHit=0; nHit<nrHits; nHit++) {
428 Int_t key = digMap->GetKeyIndex(nHit);
430 Int_t col = 32 - (key%(256*32))/256;
431 Int_t row = digMap->GetValIndex(nHit);
433 if(previousChip==-1) { // first hit
434 //loop over chip without digits
435 //Even if there aren't digits for a given chip
436 //the chip header and the chip trailer are stored
437 for (Int_t i=chipLow; i<chip; i++) {
438 WriteChipHeader(i,hs,baseWord);
439 WriteChipTrailer(buf,0,baseWord);
441 WriteChipHeader(chip,hs,baseWord);
442 WriteHit(buf,row,col,baseWord);
447 if(previousChip!=(Int_t)chip) {
448 WriteChipTrailer(buf,chipHitCount,baseWord);
450 for(Int_t i=previousChip+1; i<chip; i++) {
451 WriteChipHeader(i,hs,baseWord);
452 WriteChipTrailer(buf,0,baseWord);
454 WriteChipHeader(chip,hs,baseWord);
458 WriteHit(buf,row,col,baseWord);
461 //Even if there aren't digits for a given chip
462 //the chip header and the chip trailer are stored
463 WriteChipTrailer(buf,chipHitCount,baseWord);
465 for(Int_t i=chip+1;i<=chipHigh;i++){
466 WriteChipHeader(i,hs,baseWord);
467 WriteChipTrailer(buf,0,baseWord);
471 //In this module there aren't digits but
472 //the chip header and chip trailer are stored anyway
473 for(Int_t i=chipLow; i<=chipHigh; i++){
474 WriteChipHeader(i,hs,baseWord);
475 WriteChipTrailer(buf,0,baseWord);
486 /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
488 Int_t AliITSDDLRawData::RawDataSPD(TBranch* branch){
489 //This method creates the Raw data files for SPD detectors
490 const Int_t kSize=21000; //256*32*5=40960 max number of digits per module
491 UInt_t buf[kSize]; //One buffer cell can contain 2 digits
494 TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
496 AliFstream* outfile; // logical name of the output file
497 AliRawDataHeaderSim header;
500 for(Int_t i=0;i<AliDAQ::NumberOfDdls("ITSSPD");i++){
501 strcpy(fileName,AliDAQ::DdlFileName("ITSSPD",i)); //The name of the output file.
502 outfile = new AliFstream(fileName);
503 //write Dummy DATA HEADER
504 UInt_t dataHeaderPosition=outfile->Tellp();
505 outfile->WriteBuffer((char*)(&header),sizeof(header));
506 //Loops over Modules of a particular DDL
507 for (Int_t mod=0; mod<AliITSRawStreamSPD::kModulesPerDDL; mod++){
508 Int_t moduleNumber = AliITSRawStreamSPD::GetModuleNumber(i, mod);
510 branch->GetEvent(moduleNumber);
511 //For each Module, buf contains the array of data words in Binary format
512 //fIndex gives the number of 32 bits words in the buffer for each module
513 GetDigitsSPD(digits,moduleNumber,i,buf);
514 outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t)));
515 for(Int_t i=0;i<(fIndex+1);i++){
521 //Write REAL DATA HEADER
522 UInt_t currentFilePosition=outfile->Tellp();
523 outfile->Seekp(dataHeaderPosition);
524 header.fSize=currentFilePosition-dataHeaderPosition;
525 outfile->WriteBuffer((char*)(&header),sizeof(header));
532 /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
534 Int_t AliITSDDLRawData::RawDataSSD(TBranch* branch){
536 //This method creates the Raw data files for SSD detectors
537 const Int_t kSize=1536;//768*2 Number of stripe * number of sides(N and P)
541 TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
543 AliFstream* outfile; // logical name of the output file
544 AliRawDataHeaderSim header;
547 for(Int_t i=0;i<AliDAQ::NumberOfDdls("ITSSSD");i++){
548 strcpy(fileName,AliDAQ::DdlFileName("ITSSSD",i)); //The name of the output file.
549 outfile = new AliFstream(fileName);
550 //write Dummy DATA HEADER
551 UInt_t dataHeaderPosition=outfile->Tellp();
552 outfile->WriteBuffer((char*)(&header),sizeof(header));
554 //Loops over Modules of a particular DDL
555 for (Int_t mod=0; mod<AliITSRawStreamSSD::kModulesPerDDL; mod++){
556 Int_t moduleNumber = AliITSRawStreamSSD::GetModuleNumber(i, mod);
557 if(moduleNumber!=-1){
559 branch->GetEvent(moduleNumber);
560 //For each Module, buf contains the array of data words in Binary format
561 //fIndex gives the number of 32 bits words in the buffer for each module
562 GetDigitsSSD(digits,mod,moduleNumber,i,buf);
563 outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t)));
568 //Write REAL DATA HEADER
569 UInt_t currentFilePosition=outfile->Tellp();
570 outfile->Seekp(dataHeaderPosition);
571 header.fSize=currentFilePosition-dataHeaderPosition;
572 header.SetAttribute(0); // valid data
573 outfile->WriteBuffer((char*)(&header),sizeof(header));
580 /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
582 Int_t AliITSDDLRawData::RawDataSDD(TBranch* branch, AliITSDDLModuleMapSDD* ddlsdd){
583 //This method creates the Raw data files for SDD detectors
584 const Int_t kSize=131072; //256*512
588 TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
590 AliFstream* outfile; // logical name of the output file
591 AliRawDataHeaderSim header;
592 UInt_t skippedword, carlosFooterWord,fifoFooterWord,jitterWord;
594 retcode = AliBitPacking::PackWord(0x3FFFFFFF,carlosFooterWord,0,31);
595 retcode = AliBitPacking::PackWord(0x3F1F1F1F,fifoFooterWord,0,31);
596 retcode = AliBitPacking::PackWord(0x7F00000E,jitterWord,0,31);
599 for(Int_t i=0;i<AliDAQ::NumberOfDdls("ITSSDD");i++){
600 strcpy(fileName,AliDAQ::DdlFileName("ITSSDD",i)); //The name of the output file.
601 outfile = new AliFstream(fileName);
602 //write Dummy DATA HEADER
603 UInt_t dataHeaderPosition=outfile->Tellp();
604 outfile->WriteBuffer((char*)(&header),sizeof(header));
607 //first 1 "dummy" word to be skipped
608 retcode = AliBitPacking::PackWord(0xFFFFFFFF,skippedword,0,31);
609 outfile->WriteBuffer((char*)(&skippedword),sizeof(skippedword));
611 //Loops over Modules of a particular DDL
612 for (Int_t mod=0; mod<AliITSRawStreamSDD::kModulesPerDDL; mod++){
613 Int_t moduleNumber = ddlsdd->GetModuleNumber(i, mod);
614 if(moduleNumber!=-1){
616 branch->GetEvent(moduleNumber);
618 //For each Module, buf contains the array of data words in Binary format
619 //fIndex gives the number of 32 bits words in the buffer for each module
620 // cout<<"MODULE NUMBER:"<<mapSDD[i][mod]<<endl;
621 GetDigitsSDD(digits,mod,moduleNumber,i,buf);
622 outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t)));
623 for(Int_t iw=0;iw<3;iw++) outfile->WriteBuffer((char*)(&carlosFooterWord),sizeof(carlosFooterWord));
627 // 12 words with FIFO footers (=4 FIFO x 3 3F1F1F1F words per DDL)
628 for(Int_t iw=0;iw<12;iw++) outfile->WriteBuffer((char*)(&fifoFooterWord),sizeof(fifoFooterWord));
630 outfile->WriteBuffer((char*)(&jitterWord),sizeof(jitterWord));
632 //Write REAL DATA HEADER
633 UInt_t currentFilePosition=outfile->Tellp();
634 outfile->Seekp(dataHeaderPosition);
635 header.fSize=currentFilePosition-dataHeaderPosition;
636 header.SetAttribute(0); // valid data
637 outfile->WriteBuffer((char*)(&header),sizeof(header));
644 /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
646 void AliITSDDLRawData::WriteChipHeader(Int_t ChipAddr,Int_t halfStave,UInt_t &BaseWord){
647 //This method writes a chip header
648 //cout<<"Chip: "<<ChipAddr<<" Half Stave module:"<<halfStave<<endl;
650 AliBitPacking::PackWord(ChipAddr,BaseWord,16,19);
651 // At the moment the event count is always 0 (bits 20-26)
652 AliBitPacking::PackWord(0,BaseWord,20,26);
653 AliBitPacking::PackWord(halfStave,BaseWord,27,29);
654 AliBitPacking::PackWord(0x1,BaseWord,30,31);
656 }//end WriteChipHeader
658 /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
660 void AliITSDDLRawData::WriteChipTrailer(UInt_t *buf,Int_t ChipHitCount,UInt_t &BaseWord){
661 //This method writes a chip trailer
663 if((ChipHitCount%2)!=0){
664 AliBitPacking::PackWord(0xC000,BaseWord,16,31);
666 AliBitPacking::PackWord(ChipHitCount,BaseWord,0,13);
667 AliBitPacking::PackWord(0x0,BaseWord,14,15);
669 buf[fIndex]=BaseWord;
672 }//end WriteChipTrailer
674 /////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
676 void AliITSDDLRawData::WriteHit(UInt_t *buf,Int_t RowAddr,Int_t HitAddr,UInt_t &BaseWord){
677 //This method writs an hit
679 AliBitPacking::PackWord(HitAddr,BaseWord,16,20);
680 AliBitPacking::PackWord(RowAddr,BaseWord,21,28);
681 AliBitPacking::PackWord(2,BaseWord,30,31);
684 AliBitPacking::PackWord(HitAddr,BaseWord,0,4);
685 AliBitPacking::PackWord(RowAddr,BaseWord,5,12);
686 AliBitPacking::PackWord(2,BaseWord,14,15);
688 buf[fIndex]=BaseWord;