X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSDDLRawData.cxx;h=1d3eac594f0de0e090241fb65a0f735c331954c3;hb=6e05de7534f432aa964cafe641d20ea9797d0e28;hp=0afd33d4a227e07f20acb77aa30db6c5822cb980;hpb=7f1a504ba1f7576118d5a67e5a8eece12b64f2cb;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSDDLRawData.cxx b/ITS/AliITSDDLRawData.cxx index 0afd33d4a22..1d3eac594f0 100644 --- a/ITS/AliITSDDLRawData.cxx +++ b/ITS/AliITSDDLRawData.cxx @@ -23,45 +23,57 @@ //SSD #include -#include +//#include #include #include -#include -#include "AliITS.h" -#include "AliITSgeom.h" -#include "AliITSdigitSPD.h" -#include "AliITSdigitSDD.h" -#include "AliITSdigitSSD.h" +#include "AliITSdigit.h" #include "AliITSDDLRawData.h" - +#include "AliRawDataHeaderSim.h" +#include "AliITSRawStreamSPD.h" +#include "AliITSRawStreamSDD.h" +#include "AliITSDDLModuleMapSDD.h" +#include "AliITSRawStreamSSD.h" +#include "AliITSIntMap.h" +#include "AliBitPacking.h" +#include "AliDAQ.h" +#include "AliFstream.h" +#include "AliITSFOSignalsSPD.h" + +using std::ofstream; +using std::ios; +using std::endl; ClassImp(AliITSDDLRawData) //////////////////////////////////////////////////////////////////////////////////////// -AliITSDDLRawData::AliITSDDLRawData(){ +AliITSDDLRawData::AliITSDDLRawData(): +fVerbose(0), +fIndex(-1), +fHalfStaveModule(-1), +fSDDRawFormat(7){ //Default constructor - fIndex=-1; - fHalfStaveModule=-1; - fVerbose=0; + } //////////////////////////////////////////////////////////////////////////////////////// AliITSDDLRawData::AliITSDDLRawData(const AliITSDDLRawData &source) : - TObject(source){ + TObject(source), +fVerbose(source.fVerbose), +fIndex(source.fIndex), +fHalfStaveModule(source.fHalfStaveModule), +fSDDRawFormat(source.fSDDRawFormat){ //Copy Constructor - this->fIndex=source.fIndex; - this->fHalfStaveModule=source.fHalfStaveModule; - this->fVerbose=source.fVerbose; - return; } //////////////////////////////////////////////////////////////////////////////////////// AliITSDDLRawData& AliITSDDLRawData::operator=(const AliITSDDLRawData &source){ //Assigment operator - this->fIndex=source.fIndex; - this->fHalfStaveModule=source.fHalfStaveModule; - this->fVerbose=source.fVerbose; + if(this==&source) return *this; + fIndex=source.fIndex; + fHalfStaveModule=source.fHalfStaveModule; + fVerbose=source.fVerbose; + fSDDRawFormat=source.fSDDRawFormat; return *this; } @@ -71,6 +83,7 @@ AliITSDDLRawData& AliITSDDLRawData::operator=(const AliITSDDLRawData &source){ void AliITSDDLRawData::GetDigitsSSD(TClonesArray *ITSdigits,Int_t mod,Int_t modR,Int_t ddl,UInt_t *buf){ //This method packs the SSD digits in a proper 32 bits structure + // Revised by Enrico Fragiacomo Int_t ix; Int_t iz; Int_t is; @@ -85,21 +98,29 @@ void AliITSDDLRawData::GetDigitsSSD(TClonesArray *ITSdigits,Int_t mod,Int_t modR } for (Int_t digit=0;digitUncheckedAt(digit); - iz=digs->GetCoord1(); // If iz==0, N side and if iz=1 P side - ix=digs->GetCoord2(); // Strip Numbar + iz=digs->GetCoord1(); // If iz==0, O side and if iz=1 N side + ix=digs->GetCoord2(); // Strip Number is=digs->GetCompressedSignal(); // ADC Signal // cout<<" Module:"<4095) is = 4095; if (fVerbose==2) ftxt<<"DDL:"< 0-767 + AliBitPacking::PackWord(word,baseWord,12,22);//Strip Number + + word = mod%12; // ADC-number (12 ADCs per AD module) + word += ( word<6 ) ? 0 : 2; // ADC range 0-5 and 8-13 + AliBitPacking::PackWord(word,baseWord,24,27);//ADC Channel + + word = mod/12+1; // AD-number (AD module index ranges 1-9) + AliBitPacking::PackWord(word,baseWord,28,31);//AD slot fIndex++; buf[fIndex]=baseWord; }//end for @@ -113,18 +134,97 @@ void AliITSDDLRawData::GetDigitsSSD(TClonesArray *ITSdigits,Int_t mod,Int_t modR //Silicon Drift Detector // +void AliITSDDLRawData::GetDigitsSDDCompressed(TClonesArray *ITSdigits, Int_t mod, UInt_t *buf){ +//This method packs the SDD digits in the compressed format (32 bit per digit) +// see AliITSRawStreamSDDCompressed for details on the dta format + + UInt_t dataWord=0; + Int_t ndigits = ITSdigits->GetEntries(); + AliITSdigit *digs; + if(ndigits){ + for (Int_t digit=0;digitUncheckedAt(digit); + Int_t iz=digs->GetCoord1(); // Anode + Int_t ix=digs->GetCoord2(); // Time + Int_t is=digs->GetCompressedSignal(); // ADC Signal - 8 bit + dataWord=mod<<27; + Int_t sid=0; + if(iz>=256){ + sid=1; + iz-=256; + } + dataWord+=sid<<26; + dataWord+=iz<<18; + dataWord+=ix<<10; + UInt_t adcEncoded=0; + Int_t shift=0; + if(is < 8) shift=2; + else if(is<16) shift=3; + else if(is<32) shift=4; + else if(is<64) shift=5; + else if(is<128) shift=6; + else shift=7; + adcEncoded=shift+((is-(1<GetEntries(); AliITSdigit *digs; ofstream ftxt; + Int_t digarr[512][256]; + for(Int_t i=0;i<512;i++){ + for(Int_t j=0;j<256;j++){ + digarr[i][j]=0; + } + } + //word to select the 12 carlos for the 12 modules + UInt_t carlosid=0x30000000+mod; + + fIndex++; + buf[fIndex]=carlosid; + Int_t first=0; + Int_t last=0; + Int_t diff=0; + Int_t nbit=0; + UInt_t word2=0; + Bool_t flag = kFALSE; + baseWord=0; + Int_t bitinfo1[4] = {3,8,3,7}; //vector with info on bit for timebin info + Int_t wordinfo1[4]= {0,0,0,0}; //vector with word info for timebin info + Int_t bitinfo2[2] = {3,18}; //vector with info on bit for EOR (end of row) info + Int_t wordinfo2[3]= {1,65593}; //vector with word info for anode info + + /* for time bin info: word n bits meaning + 0 3 next info is timebin + 8 3 next word is 8 bit long + tb value 8 timebin value + n (2->7) 3 next info is n bit long + signal n signal value + + for anode info: 1 3 next 18 bits are for EOR + increments the anode value + + EOR 18 error codes + other info + */ + if(ndigits){ - //cout<<"Mudule "<GetCoord1(); // Anode ix=digs->GetCoord2(); // Time is=digs->GetCompressedSignal(); // ADC Signal + digarr[iz][ix]=is; if (fVerbose==2) ftxt<<"DDL:"<255){cout<<"WARNING (!) bits words is needed)!!!\n";} - baseWord=0; - /* - //10 bits words for amplitude value - word=is; - PackWord(baseWord,word,0,9);//ADC data - word=ix; - PackWord(baseWord,word,10,17);//Time bucket - word=iz; - PackWord(baseWord,word,18,26);//Anode Number - word=mod; - PackWord(baseWord,word,27,31);//Module number - */ + if (is>255){Error("GetDigitsSDD", "bits words is needed)!!!");} + } - //8bits words for amplitude value - word=is; - PackWord(baseWord,word,0,7);//ADC data - word=ix; - PackWord(baseWord,word,8,15);//Time bucket - word=iz; - PackWord(baseWord,word,16,24);//Anode Number - word=mod; - PackWord(baseWord,word,25,31);//Module number - - fIndex++; - buf[fIndex]=baseWord; - }//end for - }//end if + for(Int_t anode=0;anode<512;anode++){ + if(flag){ + last = first+diff-1; + AliBitPacking::PackWord(word2,baseWord,first,last); + flag = kFALSE; + first = last+1; + diff=0; + } + if(anode == 256){ + last = 0; + first = 0; + flag = kFALSE; + diff = 0; + word2=0; + } + + for(Int_t tb=0;tb<256;tb++){ + if(digarr[anode][tb]!=0){ + if(flag){ + last = first+diff-1; + AliBitPacking::PackWord(word2,baseWord,first,last); + flag = kFALSE; + first = last+1; + diff=0; + } + wordinfo1[1] = tb; + //non lossy compression as it is done in Carlos + //(data are already 10to8bit compressed by AMBRA + + /* if value < 8 value = value - (1 << 2) (word is 2 bit long) + if value < 16 value = value - (1 << 3) (word is 3 bit long) + if value < 32 value = value - (1 << 4) (word is 4 bit long) + if value < 64 value = value - (1 << 5) (word is 5 bit long) + if value <128 value = value - (1 << 6) (word is 6 bit long) + if value >=128value = value - (1 << 7) (word is 7 bit long) + + */ + //if(digarr[anode][tb]<4) continue; // not write <4 cnts above tL + if(digarr[anode][tb]<8){ + bitinfo1[3] = 2; + wordinfo1[2] = 2; + wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]); + } + if(digarr[anode][tb]>=8 && digarr[anode][tb]<16){ + bitinfo1[3] = 3; + wordinfo1[2] = 3; + wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]); + } + if(digarr[anode][tb]>=16 && digarr[anode][tb]<32){ + bitinfo1[3] = 4; + wordinfo1[2] = 4; + wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]); + } + if(digarr[anode][tb]>=32 && digarr[anode][tb]<64){ + bitinfo1[3] = 5; + wordinfo1[2] = 5; + wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]); + } + if(digarr[anode][tb]>=64 && digarr[anode][tb]<128){ + bitinfo1[3] = 6; + wordinfo1[2] = 6; + wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]); + } + if(digarr[anode][tb]>=128){ + bitinfo1[3] = 7; + wordinfo1[2] = 7; + wordinfo1[3] = digarr[anode][tb]-(1 << bitinfo1[3]); + } + + for(Int_t ie=0;ie<4;ie++){ + + if(flag){ + last = first+diff-1; + AliBitPacking::PackWord(word2,baseWord,first,last); + flag = kFALSE; + first = last+1; + diff=0; + } + last = first+bitinfo1[ie]-1; + if(first < 30 && last < 30){ + AliBitPacking::PackWord(wordinfo1[ie],baseWord,first,last); + first = last+1; + } + else{ + if(first<=29){ + UInt_t w = AliBitPacking::UnpackWord(wordinfo1[ie],0,29-first); + AliBitPacking::PackWord(w,baseWord,first,29); + Int_t lb = 29-first+1; + diff = bitinfo1[ie]-lb; + word2 = AliBitPacking::UnpackWord(wordinfo1[ie],lb,lb+diff-1); + flag = kTRUE; + if(anode<256) word = 2;//channel 0 of carlos + else word = 3; //channel 1 of carlos + AliBitPacking::PackWord(word,baseWord,30,31); + fIndex++; + buf[fIndex]=baseWord; + first=0; + last = 0; + baseWord=0; + word = 0; + } + else{ + word2 = wordinfo1[ie]; + diff = bitinfo1[ie]; + flag = kTRUE; + if(anode<256) word = 2; //channel 0 of carlos + else word = 3; //channel 1 of carlos + AliBitPacking::PackWord(word,baseWord,30,31); + fIndex++; + buf[fIndex]=baseWord; + first=0; + last=0; + baseWord=0; + word = 0; + } + } + } + + }//END IF + + }//end loop on tb + + for(Int_t i=0;i<2;i++){ + if(flag){ + last = first+diff-1; + AliBitPacking::PackWord(word2,baseWord,first,last); + flag = kFALSE; + first = last+1; + diff=0; + } + + word = wordinfo2[i]; + nbit = bitinfo2[i]; + last = first+nbit-1; + if(first < 30 && last < 30){ + AliBitPacking::PackWord(word,baseWord,first,last); //3 bit code =1 -> next 18 bits for EOR + first = last+1; + } + + else{ + if(first<=29){ + UInt_t w = AliBitPacking::UnpackWord(word,0,29-first); + AliBitPacking::PackWord(w,baseWord,first,29); + Int_t lb = 29-first+1; + diff = nbit-lb; + word2 = AliBitPacking::UnpackWord(word,lb,lb+diff-1); + flag = kTRUE; + if(anode<256) word = 2; + else word = 3; + AliBitPacking::PackWord(word,baseWord,30,31); + fIndex++; + buf[fIndex]=baseWord; + first=0; + last = 0; + baseWord=0; + if(anode==255){ + flag=kFALSE; + word2=0; + } + } + else{ + word2 = word; + diff = nbit; + flag = kTRUE; + if(anode<256) word = 2; + else word = 3; + AliBitPacking::PackWord(word,baseWord,30,31); + fIndex++; + buf[fIndex]=baseWord; + first=0; + last=0; + baseWord=0; + if(anode==255){ + flag=kFALSE; + word2=0; + } + } + } + } + } //end for + + } if(fVerbose==2) ftxt.close(); return; @@ -172,673 +429,318 @@ void AliITSDDLRawData::GetDigitsSDD(TClonesArray *ITSdigits,Int_t mod,Int_t modR //PIXEL // -void AliITSDDLRawData::GetDigitsSPD(TClonesArray *ITSdigits,Int_t mod,Int_t ddl, UInt_t *buf){ +void AliITSDDLRawData::GetDigitsSPD(TClonesArray *ITSdigits,Int_t mod,Int_t ddl, UInt_t *buf, AliITSFOSignalsSPD* foSignals){ //This method packs the SPD digits in a proper 32 structure //Since data is zero suppressed,the coordinates for the chip having zero digits //doesn't get listed in the galice.root file. However the SPD format requires //the empty chip to be written with chip header and chip trailer. - Int_t ix; - Int_t iz; - Int_t chipNo=0; + + Int_t chipLow = AliITSRawStreamSPD::GetOnlineChipFromOffline(mod,0); + Int_t chipHigh = AliITSRawStreamSPD::GetOnlineChipFromOffline(mod,159); + + if (chipLow>chipHigh) {chipLow -= 4; chipHigh += 4;} + UInt_t eq = AliITSRawStreamSPD::GetOnlineEqIdFromOffline(mod); + UInt_t hs = AliITSRawStreamSPD::GetOnlineHSFromOffline(mod); + + // create int map to later hold all digits sorted + AliITSIntMap* digMap = new AliITSIntMap(); + UInt_t baseWord=0; - UInt_t hitRow=0; - Int_t chipHitCount=0; //Number of Hit in the current chip - Int_t previousChip=-1; //Previuos chip respect to the actual aone + Int_t ndigits = ITSdigits->GetEntries(); //number of digits in the current module //cout<<" Number of digits in the current module:"<UncheckedAt(digit); + if (fVerbose==2) ftxt.open("SPDdigits.txt",ios::app); + for (Int_t digit=0; digitUncheckedAt(digit); /*--------------------------------------------------------------------------- * Each module contains 5 read out chips of 256 rows and 32 columns. - * So, the cell number in Z direction varies from 0 to 159. Therefore, - * to get the chip address (0 to 4), we need to divide column number by 32. + * So, the cell number in Z direction varies from 0 to 159. * ---------------------------------------------------------------------*/ - iz=digs->GetCoord1(); // Cell number in Z direction - ix=digs->GetCoord2(); // Cell number in X direction - chipNo=iz/32; - if(fVerbose==2) - ftxt<<"DDL:"<4) - WriteChipHeader(i+5,(mod/2),baseWord); - else - WriteChipHeader(i,(mod/2),baseWord); - WriteChipTrailer(buf,chipHitCount,baseWord); - chipHitCount=0; - }//end for - WriteChipHeader(chipNo,(mod/2),baseWord); - chipHitCount++; - WriteHit(buf,ix,hitRow,baseWord); - previousChip=chipNo; - }//end if - else{ - if(previousChip!=chipNo){ - WriteChipTrailer(buf,chipHitCount,baseWord); - chipHitCount=0; - for(Int_t i=previousChip+1;i4)end+=5; - WriteChipTrailer(buf,chipHitCount,baseWord); - chipHitCount=0; - for(Int_t i=chipNo+1;i<=end;i++){ - WriteChipHeader(i,(mod/2),baseWord); - WriteChipTrailer(buf,0,baseWord); - chipHitCount=0; - }//end for - }//end if - else{ - //In this module there aren't digits but - //the chip header and chip trailer are store anyway - if(fHalfStaveModule){ - chipNo=5; - fHalfStaveModule=-1; - }//end if - for(Int_t i=0;i<5;i++){ - WriteChipHeader(chipNo+i,(mod/2),baseWord); - WriteChipTrailer(buf,chipHitCount,baseWord); - chipHitCount=0; - }//end for - }//end else - if(fVerbose==2) - ftxt.close(); - return; -}//end GetDigitsSPD + Int_t iz=digs->GetCoord1(); // Cell number in Z direction + Int_t ix=digs->GetCoord2(); // Cell number in X direction + + if(fVerbose==2) ftxt<<"DDL:"< sum){ - cout<<"WARNING::Word to be filled is not within desired length"<Insert(chip*256*32+(31-col)*256+(255-row),row); + } } - offSet=sum; - offSet<<=StartBit; - BaseWord=BaseWord|offSet; - //The Word to be filled is shifted to the position StartBit - //and the remaining Left and Right bits are filled with 1 - sum=(UInt_t)TMath::Power(2,StartBit)-1; - dummyWord=0xFFFFFFFF<GetNrEntries(); + for (UInt_t nHit=0; nHit>=StartBit; + Int_t key = digMap->GetKeyIndex(nHit); + chip = key/(256*32); + Int_t col = 31 - (key%(256*32))/256; + Int_t row = digMap->GetValIndex(nHit); + + // add trailer for previous chip (if there was one...) + if (chip>previousChip && previousChip>chipLow-1) { + WriteChipTrailer(buf, chipHitCount, foSignals->GetSignal(eq,hs,previousChip), baseWord); + } + + // add headers/trailers for chips without hits (if any) + for (Int_t ch=previousChip+1; chGetSignal(eq,hs,ch), baseWord); + } + + // if new chip, add header + if (chip>previousChip) { + WriteChipHeader(chip, hs, baseWord); + chipHitCount = 0; + previousChip = chip; + } + + chipHitCount++; + + // add pixel hit + WriteHit(buf,row,col,baseWord); + + } + + // add trailer for last chip (if there was one...) + if (chip>chipLow-1) { + WriteChipTrailer(buf, chipHitCount, foSignals->GetSignal(eq,hs,chip), baseWord); + } + + // add REMAINING headers/trailers for chips without hits (if any) + for (Int_t ch=chip+1; ch<=chipHigh; ch++) { + WriteChipHeader(ch, hs, baseWord); + WriteChipTrailer(buf, 0, foSignals->GetSignal(eq,hs,ch), baseWord); + } + // _______________________________________________________________________ + + + delete digMap; + + if(fVerbose==2) + ftxt.close(); return; -} +}//end GetDigitsSPD ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////// -Int_t AliITSDDLRawData::RawDataSPD(AliITS *ITS,TTree *TD ,Int_t LDCsNumber,Int_t eventNumber){ +Int_t AliITSDDLRawData::RawDataSPD(TBranch* branch, AliITSFOSignalsSPD* foSignals){ //This method creates the Raw data files for SPD detectors - //Silicon Pixel Detector - const Int_t kDDLsNumber=20; // Number of DDL in SPD - const Int_t kModulesPerDDL=12; // Number of modules in each DDL - //Row ==> DDLs - //Column ==> Modules - Int_t mapSPD[kDDLsNumber][kModulesPerDDL]={{ 0, 1, 4, 5, 80, 81, 84, 85, 88, 89, 92, 93}, - { 2, 3, 6, 7, 82, 83, 86, 87, 90, 91, 94, 95}, - { 8, 9,12,13, 96, 97,100,101,104,105,108,109}, - {10,11,14,15, 98, 99,102,103,106,107,110,111}, - {16,17,20,21,112,113,116,117,120,121,124,125}, - {18,19,22,23,114,115,118,119,122,123,126,127}, - {24,25,28,29,128,129,132,133,136,137,140,141}, - {26,27,30,31,130,131,134,135,138,139,142,143}, - {32,33,36,37,144,145,148,149,152,153,156,157}, - {34,35,38,39,146,147,150,151,154,155,158,159}, - {40,41,44,45,160,161,164,165,168,169,172,173}, - {42,43,46,47,162,163,166,167,170,171,174,175}, - {48,49,52,53,176,177,180,181,184,185,188,189}, - {50,51,54,55,178,179,182,183,186,187,190,191}, - {56,57,60,61,192,193,196,197,200,201,204,205}, - {58,59,62,63,194,195,198,199,202,203,206,207}, - {64,65,68,69,208,209,212,213,216,217,220,221}, - {66,67,70,71,210,211,214,215,218,219,222,223}, - {72,73,76,77,224,225,228,229,232,233,236,237}, - {74,75,78,79,226,227,230,231,234,235,238,239}}; - Int_t ddlsPerFile=kDDLsNumber/LDCsNumber; - if(kDDLsNumber%LDCsNumber)ddlsPerFile++; - if(fVerbose) - cout<<"Number of DDL per File: "<DigitsAddress(subd); - - Int_t nbytes = 0; - char fileName[15]; - - UInt_t miniHeaderPosition=0; //variable used to store the position of the Mini Header inside a file - ofstream outfile; // logical name of the output file - Int_t flag=0; // 0==> Uncompressed data 1==>Compressed Data - Int_t detector=1; // 1==>ITS (Pixel) 2==>ITS (Drift) 3==>ITS (Strip) 0==>TPC ...... - UInt_t size=0; // size of the data block that follows the mini header - Int_t magicWord=0x123456; // Magic word used to distinguish between data and garbage - sprintf(fileName,"Ev%dSPDslice%d",eventNumber,sliceNumber); //The name of the output file. There are as many slides as the number of LDC -#ifndef __DECCXX - outfile.open(fileName,ios::binary); -#else - outfile.open(fileName); -#endif - UInt_t miniHeader[3]; - Int_t miniHeaderSize=sizeof(UInt_t)*3; - Int_t version=1; //version of the mini header + + TClonesArray*& digits = * (TClonesArray**) branch->GetAddress(); + TString fileName; + AliFstream* outfile; // logical name of the output file + AliRawDataHeaderSim header; + //loop over DDLs - for(Int_t i=0;iTellp(); + outfile->WriteBuffer((char*)(&header),sizeof(header)); //Loops over Modules of a particular DDL - for (Int_t mod=0; modResetDigits(); - nbytes += TD->GetEvent(mapSPD[i][mod]); + for (Int_t mod=0; modClear(); + branch->GetEvent(moduleNumber); //For each Module, buf contains the array of data words in Binary format //fIndex gives the number of 32 bits words in the buffer for each module - GetDigitsSPD(digitsInITS,mapSPD[i][mod],i,buf); - outfile.write((char *)buf,((fIndex+1)*sizeof(UInt_t))); + GetDigitsSPD(digits, moduleNumber, ddl, buf, foSignals); + outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t))); for(Int_t i=0;i<(fIndex+1);i++){ buf[i]=0; }//end for fIndex=-1; }//end for - //Write REAL MINI HEADER - UInt_t currentFilePosition=outfile.tellp(); - outfile.seekp(miniHeaderPosition); - size=currentFilePosition-miniHeaderPosition-miniHeaderSize; - - outfile.write((char*)(&size),sizeof(UInt_t)); - outfile.seekp(currentFilePosition); - if(countDDL==ddlsPerFile){ - outfile.close(); - sliceNumber++; - sprintf(fileName,"Ev%dSPDslice%d",eventNumber,sliceNumber); - if(i!=(kDDLsNumber-1)) -#ifndef __DECCXX - outfile.open(fileName,ios::binary); -#else - outfile.open(fileName); -#endif - countDDL=0; - }//end if + //Write REAL DATA HEADER + UInt_t currentFilePosition=outfile->Tellp(); + outfile->Seekp(dataHeaderPosition); + header.fSize=currentFilePosition-dataHeaderPosition; + outfile->WriteBuffer((char*)(&header),sizeof(header)); + delete outfile; }//end for - outfile.close(); + return 0; } ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////// -Int_t AliITSDDLRawData::RawDataSSD(AliITS *ITS,TTree *TD ,Int_t LDCsNumber,Int_t eventNumber){ - //This method creates the Raw data files for SSD detectors - //Strip detector - const Int_t kDDLsNumber=16; // Number of DDL in SSD - const Int_t kModulesPerDDL=109; // Number of modules in each DDL - //DDL from 32 to 47 (16 DDL) - //Row ==> DDLs - //Column ==> Modules - Int_t mapSSD[kDDLsNumber][kModulesPerDDL]={ - //104 - //DDL[32][]= - { 500, 501, 502, 503, 504, 505, 506, 507, 508, 509, 510, - 522, 523, 524, 525, 526, 527, 528, 529, 530, 531, 532, - 1204,1205,1206,1207,1208,1209,1210,1211,1212,1213,1214, - 1226,1227,1228,1229,1230,1231,1232,1233,1234,1235,1236, - 1248,1249,1250,1251,1252,1253,1254,1255,1256,1257,1258,1259, - 2098,2099,2100,2101,2102,2103,2104,2105,2106,2107,2108,2109, - 2123,2124,2125,2126,2127,2128,2129,2130,2131,2132,2133,2134, - 2148,2149,2150,2151,2152,2153,2154,2155,2156,2157,2158,2159, - 2173,2174,2175,2176,2177,2178,2179,2180,2181,2182,2183,2184,-1,-1,-1,-1,-1}, - //93 - //DDL[33][]= - { 544, 545, 546, 547, 548, 549, 550, 551, 552, 553, 554, - 566, 567, 568, 569, 570, 571, 572, 573, 574, 575, 576, - 588, 589, 590, 591, 592, 593, 594, 595, 596, 597, 598, - 610, 611, 612, 613, 614, 615, 616, 617, 618, 619, 620, - 1273,1274,1275,1276,1277,1278,1279,1280,1281,1282,1283,1284, - 1298,1299,1300,1301,1302,1303,1304,1305,1306,1307,1308,1309, - 1323,1324,1325,1326,1327,1328,1329,1330,1331,1332,1333,1334, - 1348,1349,1350,1351,1352,1353,1354,1355,1356,1357,1358,1359, - 1373,1374,1375,1376,1377,1378,1379,1380,1381,1382,1383,1384,-1,-1,-1,-1,-1}, - //103 - //DDL[34][]= - { 632, 633, 634, 635, 636, 637, 638, 639, 640, 641, 642, - 654, 655, 656, 657, 658, 659, 660, 661, 662, 663, 664, - 676, 677, 678, 679, 680, 681, 682, 683, 684, 685, 686, - 698, 699, 700, 701, 702, 703, 704, 705, 706, 707, 708, - 720, 721, 722, 723, 724, 725, 726, 727, 728, 729, 730, - 1398,1399,1400,1401,1402,1403,1404,1405,1406,1407,1408,1409, - 1423,1424,1425,1426,1427,1428,1429,1430,1431,1432,1433,1434, - 1448,1449,1450,1451,1452,1453,1454,1455,1456,1457,1458,1459, - 1473,1474,1475,1476,1477,1478,1479,1480,1481,1482,1483,1484,-1,-1,-1,-1,-1,-1}, - //104 - //DDL[35][]= - { 742, 743, 744, 745, 746, 747, 748, 749, 750, 751, 752, - 764, 765, 766, 767, 768, 769, 770, 771, 772, 773, 774, - 786, 787, 788, 789, 790, 791, 792, 793, 794, 795, 796, - 808, 809, 810, 811, 812, 813, 814, 815, 816, 817, 818, - 1498,1499,1500,1501,1502,1503,1504,1505,1506,1507,1508,1509, - 1523,1524,1525,1526,1527,1528,1529,1530,1531,1532,1533,1534, - 1548,1549,1550,1551,1552,1553,1554,1555,1556,1557,1558,1559, - 1573,1574,1575,1576,1577,1578,1579,1580,1581,1582,1583,1584, - 1598,1599,1600,1601,1602,1603,1604,1605,1606,1607,1608,1609,-1,-1,-1,-1,-1}, - //104 - //DDL[36[]= - { 830, 831, 832, 833, 834, 835, 836, 837, 838, 839, 840, - 852, 853, 854, 855, 856, 857, 858, 859, 860, 861, 862, - 874, 875, 876, 877, 878, 879, 880, 881, 882, 883, 884, - 896, 897, 898, 899, 900, 901, 902, 903, 904, 905, 906, - 1623,1624,1625,1626,1627,1628,1629,1630,1631,1632,1633,1634, - 1648,1649,1650,1651,1652,1653,1654,1655,1656,1657,1658,1659, - 1673,1674,1675,1676,1677,1678,1679,1680,1681,1682,1683,1684, - 1698,1699,1700,1701,1702,1703,1704,1705,1706,1707,1708,1709, - 1723,1724,1725,1726,1727,1728,1729,1730,1731,1732,1733,1734,-1,-1,-1,-1,-1}, - //104 - //DDL[37][]= - { 918, 919, 920, 921, 922, 923, 924, 925, 926, 927, 928, - 940, 941, 942, 943, 944, 945, 946, 947, 948, 949, 950, - 962, 963, 964, 965, 966, 967, 968, 969, 970, 971, 972, - 984, 985, 986, 987, 988, 989, 990, 991, 992, 993, 994, - 1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759, - 1773,1774,1775,1776,1777,1778,1779,1780,1781,1782,1783,1784, - 1798,1799,1800,1801,1802,1803,1804,1805,1806,1807,1808,1809, - 1823,1824,1825,1826,1827,1828,1829,1830,1831,1832,1833,1834, - 1848,1849,1850,1851,1852,1853,1854,1855,1856,1857,1858,1859,-1,-1,-1,-1,-1}, - //103 - //DDL[38][]= - {1006,1007,1008,1009,1010,1011,1012,1013,1014,1015,1016, - 1028,1029,1030,1031,1032,1033,1034,1035,1036,1037,1038, - 1050,1051,1052,1053,1054,1055,1056,1057,1058,1059,1060, - 1072,1073,1074,1075,1076,1077,1078,1079,1080,1081,1082, - 1094,1095,1096,1097,1098,1099,1100,1101,1102,1103,1104, - 1873,1874,1875,1876,1877,1878,1879,1880,1881,1882,1883,1884, - 1898,1899,1900,1901,1902,1903,1904,1905,1906,1907,1908,1909, - 1923,1924,1925,1926,1927,1928,1929,1930,1931,1932,1933,1934, - 1948,1949,1950,1951,1952,1953,1954,1955,1956,1957,1958,1959,-1,-1,-1,-1,-1,-1}, - //104 - //DDL[39][]= - {1116,1117,1118,1119,1120,1121,1122,1123,1124,1125,1126, - 1138,1139,1140,1141,1142,1143,1144,1145,1146,1147,1148, - 1160,1161,1162,1163,1164,1165,1166,1167,1168,1169,1170, - 1182,1183,1184,1185,1186,1187,1188,1189,1190,1191,1192, - 1973,1974,1975,1976,1977,1978,1979,1980,1981,1982,1983,1984, - 1998,1999,2000,2001,2002,2003,2004,2005,2006,2007,2008,2009, - 2023,2024,2025,2026,2027,2028,2029,2030,2031,2032,2033,2034, - 2048,2049,2050,2051,2052,2053,2054,2055,2056,2057,2058,2059, - 2073,2074,2075,2076,2077,2078,2079,2080,2081,2082,2083,2084,-1,-1,-1,-1,-1}, - //109 - //DDL[40][]= - { 511, 512, 513, 514, 515, 516, 517, 518, 519, 520, 521, - 533, 534, 535, 536, 537, 538, 539, 540, 541, 542, 543, - 1215,1216,1217,1218,1219,1220,1221,1222,1223,1224,1225, - 1237,1238,1239,1240,1241,1242,1243,1244,1245,1246,1247, - 1260,1261,1262,1263,1264,1265,1266,1267,1268,1269,1270,1271,1272, - 2110,2111,2112,2113,2114,2115,2116,2117,2118,2119,2120,2121,2122, - 2135,2136,2137,2138,2139,2140,2141,2142,2143,2144,2145,2146,2147, - 2160,2161,2162,2163,2164,2165,2166,2167,2168,2169,2170,2171,2172, - 2185,2186,2187,2188,2189,2190,2191,2192,2193,2194,2195,2196,2197}, - //109 - //DDL[41][]= - { 555, 556, 557, 558, 559, 560, 561, 562, 563, 564, 565, - 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 587, - 599, 600, 601, 602, 603, 604, 605, 606, 607, 608, 609, - 621, 622, 623, 624, 625, 626, 627, 628, 629, 630, 631, - 1285,1286,1287,1288,1289,1290,1291,1292,1293,1294,1295,1296,1297, - 1310,1311,1312,1313,1314,1315,1316,1317,1318,1319,1320,1321,1322, - 1335,1336,1337,1338,1339,1340,1341,1342,1343,1344,1345,1346,1347, - 1360,1361,1362,1363,1364,1365,1366,1367,1368,1369,1370,1371,1372, - 1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397}, - //107 - //DDL[42][]= - { 643, 644, 645, 646, 647, 648, 649, 650, 651, 652, 653, - 665, 666, 667, 668, 669, 670, 671, 672, 673, 674, 675, - 687, 688, 689, 690, 691, 692, 693, 694, 695, 696, 697, - 709, 710, 711, 712, 713, 714, 715, 716, 717, 718, 719, - 731, 732, 733, 734, 735, 736, 737, 738, 739, 740, 741, - 1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422, - 1435,1436,1437,1438,1439,1440,1441,1442,1443,1444,1445,1446,1447, - 1460,1461,1462,1463,1464,1465,1466,1467,1468,1469,1470,1471,1472, - 1485,1486,1487,1488,1489,1490,1491,1492,1493,1494,1495,1496,1497,-1,-1}, - //109 - //DDL[43][]= - { 753, 754, 755, 756, 757, 758, 759, 760, 761, 762, 763, - 775, 776, 777, 778, 779, 780, 781, 782, 783, 784, 785, - 797, 798, 799, 800, 801, 802, 803, 804, 805, 806, 807, - 819, 820, 821, 822, 823, 824, 825, 826, 827, 828, 829, - 1510,1511,1512,1513,1514,1515,1516,1517,1518,1519,1520,1521,1522, - 1535,1536,1537,1538,1539,1540,1541,1542,1543,1544,1545,1546,1547, - 1560,1561,1562,1563,1564,1565,1566,1567,1568,1569,1570,1571,1572, - 1585,1586,1587,1588,1589,1590,1591,1592,1593,1594,1595,1596,1597, - 1610,1611,1612,1613,1614,1615,1616,1617,1618,1619,1620,1621,1622}, - //109 - //DDL[44][]= - { 841, 842, 843, 844, 845, 846, 847, 848, 849, 850, 851, - 863, 864, 865, 866, 867, 868, 869, 870, 871, 872, 873, - 885, 886, 887, 888, 889, 890, 891, 892, 893, 894, 895, - 907, 908, 909, 910, 911, 912, 913, 914, 915, 916, 917, - 1635,1636,1637,1638,1639,1640,1641,1642,1643,1644,1645,1646,1647, - 1660,1661,1662,1663,1664,1665,1666,1667,1668,1669,1670,1671,1672, - 1685,1686,1687,1688,1689,1690,1691,1692,1693,1694,1695,1696,1697, - 1710,1711,1712,1713,1714,1715,1716,1717,1718,1719,1720,1721,1722, - 1735,1736,1737,1738,1739,1740,1741,1742,1743,1744,1745,1746,1747}, - //109 - //DDL[45][]= - {929, 930, 931, 932, 933, 934, 935, 936, 937, 938, 939, - 951, 952, 953, 954, 955, 956, 957, 958, 959, 960, 961, - 973, 974, 975, 976, 977, 978, 979, 980, 981, 982, 983, - 995, 996, 997, 998, 999,1000,1001,1002,1003,1004,1005, - 1760,1761,1762,1763,1764,1765,1766,1767,1768,1769,1770,1771,1772, - 1785,1786,1787,1788,1789,1790,1791,1792,1793,1794,1795,1796,1797, - 1810,1811,1812,1813,1814,1815,1816,1817,1818,1819,1820,1821,1822, - 1835,1836,1837,1838,1839,1840,1841,1842,1843,1844,1845,1846,1847, - 1860,1861,1862,1863,1864,1865,1866,1867,1868,1869,1870,1871,1872}, - //109 - //DDL[46][]= - {1017,1018,1019,1020,1021,1022,1023,1024,1025,1026,1027, - 1039,1040,1041,1042,1043,1044,1045,1046,1047,1048,1049, - 1061,1062,1063,1064,1065,1066,1067,1068,1069,1070,1071, - 1083,1084,1085,1086,1087,1088,1089,1090,1091,1092,1093, - 1105,1106,1107,1108,1109,1110,1111,1112,1113,1114,1115, - 1885,1886,1887,1888,1889,1890,1891,1892,1893,1894,1895,1896,1897, - 1910,1911,1912,1913,1914,1915,1916,1917,1918,1919,1920,1921,1922, - 1935,1936,1937,1938,1939,1940,1941,1942,1943,1944,1945,1946,1947, - 1960,1961,1962,1963,1964,1965,1966,1967,1968,1969,1970,1971,1972}, - //109 - //DDL[47][]= - {1127,1128,1129,1130,1131,1132,1133,1134,1135,1136,1137, - 1149,1150,1151,1152,1153,1154,1155,1156,1157,1158,1159, - 1171,1172,1173,1174,1175,1176,1177,1178,1179,1180,1181, - 1193,1194,1195,1196,1197,1198,1199,1200,1201,1202,1203, - 1985,1986,1987,1988,1989,1990,1991,1992,1993,1994,1995,1996,1997, - 2010,2011,2012,2013,2014,2015,2016,2017,2018,2019,2020,2021,2022, - 2035,2036,2037,2038,2039,2040,2041,2042,2043,2044,2045,2046,2047, - 2060,2061,2062,2063,2064,2065,2066,2067,2068,2069,2070,2071,2072, - 2085,2086,2087,2088,2089,2090,2091,2092,2093,2094,2095,2096,2097} }; - - - Int_t ddlsPerFile=kDDLsNumber/LDCsNumber; - if(kDDLsNumber%LDCsNumber)ddlsPerFile++; - if(fVerbose) - cout<<"Number of DDL per File: "<DigitsAddress(subd); - char fileName[15]; - - UInt_t miniHeaderPosition=0; //variable used to store the position of the Mini Header inside the file - ofstream outfile; // logical name of the output file - Int_t flag=0; // 0==> Uncompressed data 1==>Compressed Data - Int_t detector=3; // 1==>ITS (Pixel) 2==>ITS (Drift) 3==>ITS (Strip) 0==>TPC ...... - UInt_t size=0; // size of the data block that follows the mini header - Int_t magicWord=0x123456; // Magic word used to distinguish between data and garbage - sprintf(fileName,"Ev%dSSDslice%d",eventNumber,sliceNumber); //The name of the output file. There are as many slides as the number of LDC -#ifndef __DECCXX - outfile.open(fileName,ios::binary); -#else - outfile.open(fileName); -#endif - UInt_t miniHeader[3]; - Int_t miniHeaderSize=sizeof(UInt_t)*3; - Int_t version=1; //version of the mini header - //loop over DDLs - - for(Int_t i=0;iGetAddress(); + TString fileName; + AliFstream* outfile; // logical name of the output file + AliRawDataHeaderSim header; + + //loop over DDLs + for(Int_t i=0;iTellp(); + outfile->WriteBuffer((char*)(&header),sizeof(header)); //Loops over Modules of a particular DDL - for (Int_t mod=0; modResetDigits(); - nbytes += TD->GetEvent(mapSSD[i][mod]); + for (Int_t mod=0; modClear(); + branch->GetEvent(moduleNumber); //For each Module, buf contains the array of data words in Binary format //fIndex gives the number of 32 bits words in the buffer for each module - GetDigitsSSD(digitsInITS,mod,mapSSD[i][mod],i,buf); - outfile.write((char *)buf,((fIndex+1)*sizeof(UInt_t))); - for(Int_t i=0;i<(fIndex+1);i++){ - buf[i]=0; - }//end for + GetDigitsSSD(digits,mod,moduleNumber,i,buf); + outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t))); fIndex=-1; }//end if }//end for - //Write REAL MINI HEADER - UInt_t currentFilePosition=outfile.tellp(); - outfile.seekp(miniHeaderPosition); - size=currentFilePosition-miniHeaderPosition-miniHeaderSize; - - outfile.write((char*)(&size),sizeof(UInt_t)); - outfile.seekp(currentFilePosition); - if(countDDL==ddlsPerFile){ - outfile.close(); - sliceNumber++; - sprintf(fileName,"Ev%dSSDslice%d",eventNumber,sliceNumber); - if(i!=(kDDLsNumber-1)) -#ifndef __DECCXX - outfile.open(fileName,ios::binary); -#else - outfile.open(fileName); -#endif - countDDL=0; - }//end if + + //Write REAL DATA HEADER + UInt_t currentFilePosition=outfile->Tellp(); + outfile->Seekp(dataHeaderPosition); + header.fSize=currentFilePosition-dataHeaderPosition; + header.SetAttribute(0); // valid data + outfile->WriteBuffer((char*)(&header),sizeof(header)); + delete outfile; }//end for - outfile.close(); + return 0; } ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////// -Int_t AliITSDDLRawData::RawDataSDD(AliITS *ITS,TTree *TD ,Int_t LDCsNumber,Int_t eventNumber){ +Int_t AliITSDDLRawData::RawDataSDD(TBranch* branch, const AliITSDDLModuleMapSDD* ddlsdd){ //This method creates the Raw data files for SDD detectors - //Silicon Drift detector - const Int_t kDDLsNumber=12; // Number of DDL in SPD - const Int_t kModulesPerDDL=22; // Number of modules in each DDL - //Row ==> DDLs - //Column ==> Modules - Int_t mapSDD[kDDLsNumber][kModulesPerDDL]= {{240,241,242,246,247,248,252,253,254,258,259,260,264,265,266,270,271,272,276,277,278,-1}, - {243,244,245,249,250,251,255,256,257,261,262,263,267,268,269,273,274,275,279,280,281,-1}, - {282,283,284,288,289,290,294,295,296,300,301,302,306,307,308,312,313,314,318,319,320,-1}, - {285,286,287,291,292,293,297,298,299,303,304,305,309,310,311,315,316,317,321,322,323,-1}, - {324,325,326,327,332,333,334,335,340,341,342,343,348,349,350,351,356,357,358,359,364,365}, - {328,329,330,331,336,337,338,339,344,345,346,347,352,353,354,355,360,361,362,363,368,369}, - {366,367,372,373,374,375,380,381,382,383,388,389,390,391,396,397,398,399,404,405,406,407}, - {370,371,376,377,378,379,384,385,386,387,392,393,394,395,400,401,402,403,408,409,410,411}, - {412,413,414,415,420,421,422,423,428,429,430,431,436,437,438,439,444,445,446,447,452,453}, - {416,417,418,419,424,425,426,427,432,433,434,435,440,441,442,443,448,449,450,451,456,457}, - {454,455,460,461,462,463,468,469,470,471,476,477,478,479,484,485,486,487,492,493,494,495}, - {458,459,464,465,466,467,472,473,474,475,480,481,482,483,488,489,490,491,496,497,498,499}}; - - Int_t ddlsPerFile=kDDLsNumber/LDCsNumber; - if(kDDLsNumber%LDCsNumber)ddlsPerFile++; - if(fVerbose) - cout<<"Number of DDL per File: "<DigitsAddress(subd); - char fileName[15]; - - UInt_t miniHeaderPosition=0; // variable used to store the position of the Mini Header inside a file - ofstream outfile; // logical name of the output file - Int_t flag=0; // 0==> Uncompressed data 1==>Compressed Data - Int_t detector=2; // 1==>ITS (Pixel) 2==>ITS (Drift) 3==>ITS (Strip) 0==>TPC ...... - UInt_t size=0; // size of the data block that follows the mini header - Int_t magicWord=0x123456; // Magic word used to distinguish between data and garbage - sprintf(fileName,"Ev%dSDDslice%d",eventNumber,sliceNumber); //The name of the output file. There are as many slides as the number of LDC -#ifndef __DECCXX - outfile.open(fileName,ios::binary); -#else - outfile.open(fileName); -#endif - UInt_t miniHeader[3]; - Int_t miniHeaderSize=sizeof(UInt_t)*3; - Int_t version=1; //version of the mini header - //loop over DDLs - for(Int_t i=0;iGetAddress(); + TString fileName; + AliFstream* outfile; // logical name of the output file + AliRawDataHeaderSim header; + + if(fSDDRawFormat!=0){ + for(Int_t ibit=0; ibit<8; ibit++) header.SetAttribute(ibit); + }else{ + for(Int_t ibit=0; ibit<5; ibit++) header.SetAttribute(ibit); + for(Int_t ibit=5; ibit<8; ibit++) header.ResetAttribute(ibit); + } + UInt_t skippedword=0; + UInt_t carlosFooterWord=0; + UInt_t fifoFooterWord=0; + UInt_t jitterWord=0; + Bool_t retcode; + retcode = AliBitPacking::PackWord(0x3FFFFFFF,carlosFooterWord,0,31); + retcode = AliBitPacking::PackWord(0x3F1F1F1F,fifoFooterWord,0,31); + if(fSDDRawFormat!=0) retcode = AliBitPacking::PackWord(0x7F000000,jitterWord,0,31); + else retcode = AliBitPacking::PackWord(0x80000000,jitterWord,0,31); + + //loop over DDLs + for(Int_t i=0;iTellp(); + outfile->WriteBuffer((char*)(&header),sizeof(header)); + + + //first 1 "dummy" word to be skipped + if(fSDDRawFormat!=0){ + retcode = AliBitPacking::PackWord(0xFFFFFFFF,skippedword,0,31); + outfile->WriteBuffer((char*)(&skippedword),sizeof(skippedword)); + } + //Loops over Modules of a particular DDL - for (Int_t mod=0; modResetDigits(); - nbytes += TD->GetEvent(mapSDD[i][mod]); + for (Int_t mod=0; modGetModuleNumber(i, mod); + if(moduleNumber!=-1){ + digits->Clear(); + branch->GetEvent(moduleNumber); + //For each Module, buf contains the array of data words in Binary format //fIndex gives the number of 32 bits words in the buffer for each module // cout<<"MODULE NUMBER:"<WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t))); + }else{ + GetDigitsSDD(digits,mod,moduleNumber,i,buf); + outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t))); + for(Int_t iw=0;iw<3;iw++) outfile->WriteBuffer((char*)(&carlosFooterWord),sizeof(carlosFooterWord)); + } fIndex=-1; }//end if }//end for - - //Write REAL MINI HEADER - UInt_t currentFilePosition=outfile.tellp(); - outfile.seekp(miniHeaderPosition); - size=currentFilePosition-miniHeaderPosition-miniHeaderSize; - - outfile.write((char*)(&size),sizeof(UInt_t)); - outfile.seekp(currentFilePosition); - if(countDDL==ddlsPerFile){ - outfile.close(); - sliceNumber++; - sprintf(fileName,"Ev%dSDDslice%d",eventNumber,sliceNumber); - if(i!=(kDDLsNumber-1)) -#ifndef __DECCXX - outfile.open(fileName,ios::binary); -#else - outfile.open(fileName); -#endif - countDDL=0; - }//end if + // 12 words with FIFO footers (=4 FIFO x 3 3F1F1F1F words per DDL) + if(fSDDRawFormat!=0){ + for(Int_t iw=0;iw<12;iw++) outfile->WriteBuffer((char*)(&fifoFooterWord),sizeof(fifoFooterWord)); + } + outfile->WriteBuffer((char*)(&jitterWord),sizeof(jitterWord)); + //Write REAL DATA HEADER + UInt_t currentFilePosition=outfile->Tellp(); + outfile->Seekp(dataHeaderPosition); + header.fSize=currentFilePosition-dataHeaderPosition; + header.SetAttribute(0); // valid data + outfile->WriteBuffer((char*)(&header),sizeof(header)); + delete outfile; }//end for - outfile.close(); + return 0; } ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////// -void AliITSDDLRawData::WriteChipHeader(Int_t ChipAddr,Int_t EventCnt,UInt_t &BaseWord){ +void AliITSDDLRawData::WriteChipHeader(Int_t ChipAddr,Int_t halfStave,UInt_t &BaseWord){ //This method writes a chip header - //cout<<"Chip: "<