//SSD
#include <stdlib.h>
+//#include <Riostream.h>
#include <TClonesArray.h>
#include <TTree.h>
#include "AliITSdigit.h"
#include "AliITSDDLRawData.h"
-#include "AliRawDataHeader.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"
ClassImp(AliITSDDLRawData)
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;
}
for (Int_t digit=0;digit<ndigits;digit++) {
digs = (AliITSdigit*)ITSdigits->UncheckedAt(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:"<<mod-500<<" N/P side:"<<iz<<" Strip Number:"<<ix<<" Amplidute:"<<is-1<<endl;
+ if(is<0) is = 4096 + is;
if (fVerbose==2)
ftxt<<"DDL:"<<ddl<<" Mod: "<<modR<<" N/P: "<<iz<<" Strip: "<<ix<<" Value: "<<is-1<<endl;
+
baseWord=0;
- word=is-1;
- AliBitPacking::PackWord(word,baseWord,0,9);//ADC data
- word=ix;
- AliBitPacking::PackWord(word,baseWord,10,19);//Strip Number
- word=iz;
- AliBitPacking::PackWord(word,baseWord,20,20);//ADC Channel ID (N or P side)
- word=mod;
- AliBitPacking::PackWord(word,baseWord,21,31);//ADC module ID
+
+ word=is;
+ AliBitPacking::PackWord(word,baseWord,0,11);//ADC data
+
+ word = (iz==0) ? ix : 1535-ix ; // on N-side 1535-768 -> 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
}
}
//word to select the 12 carlos for the 12 modules
- UInt_t carlosid=0;
- if(mod==0) carlosid=805306368;
- if(mod==1) carlosid=805306369;
- if(mod==2) carlosid=805306370;
- if(mod==3) carlosid=805306371;
- if(mod==4) carlosid=805306372;
- if(mod==5) carlosid=805306373;
- if(mod==6) carlosid=805306374;
- if(mod==7) carlosid=805306375;
- if(mod==8) carlosid=805306376;
- if(mod==9) carlosid=805306377;
- if(mod==10) carlosid=805306378;
- if(mod==11) carlosid=805306379;
+ UInt_t carlosid=0x30000000+mod;
fIndex++;
buf[fIndex]=carlosid;
ix=digs->GetCoord2(); // Time
is=digs->GetCompressedSignal(); // ADC Signal
digarr[iz][ix]=is;
- if (fVerbose==2)
+ if (fVerbose==2)
ftxt<<"DDL:"<<ddl<<" MID:"<<modR<<" An:"<<iz<<" T:"<<ix<<" A:"<<is<<endl;
if (is>255){Error("GetDigitsSDD", "bits words is needed)!!!");}
}
for(Int_t anode=0;anode<512;anode++){
- if(flag){
+ 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++){
//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.
- //The index of the half stave is calculated as (mod/2).
- 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 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:"<<ndigits<<" module:"<<mod<<endl;
+
AliITSdigit *digs;
- fHalfStaveModule++; //It's a private variable used to distinguish between the firs
- //and the second module of an Half Stave Module
ofstream ftxt;
if(ndigits){
//loop over digits
digs = (AliITSdigit*)ITSdigits->UncheckedAt(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;
+ 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:"<<ddl<<" Mod:"<<mod<<" Row:"<<ix<<" Col:"<<iz<<endl;
- hitRow=iz-chipNo*32;
- if(fHalfStaveModule){
- chipNo+=5;
- fHalfStaveModule=-1;
- }//end if
- if(previousChip==-1){
+ UInt_t dummyDDL,dummyHS,chip,col,row;
+ AliITSRawStreamSPD::OfflineToOnline(mod,iz,ix,dummyDDL,dummyHS,chip,col,row);
+
+ // insert digit into map...
+ // (reverse order of cols and rows as in real raw data)
+ digMap->Insert(chip*256*32+(31-col)*256+(255-row),row);
+ }
+ }
+
+ UInt_t nrHits = digMap->GetNrEntries();
+ if (nrHits>0) {
+ Int_t chip = 0;
+ for (UInt_t nHit=0; nHit<nrHits; nHit++) {
+ Int_t key = digMap->GetKeyIndex(nHit);
+ chip = key/(256*32);
+ Int_t col = 31 - (key%(256*32))/256;
+ Int_t row = digMap->GetValIndex(nHit);
+
+ if(previousChip==-1) { // first hit
//loop over chip without digits
//Even if there aren't digits for a given chip
//the chip header and the chip trailer are stored
- for(Int_t i=0;i<(iz/32);i++){
- if(chipNo>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);
+ for (Int_t i=chipLow; i<chip; i++) {
+ WriteChipHeader(i,hs,baseWord);
+ WriteChipTrailer(buf,0,baseWord);
+ }
+ WriteChipHeader(chip,hs,baseWord);
+ WriteHit(buf,row,col,baseWord);
chipHitCount++;
- WriteHit(buf,ix,hitRow,baseWord);
- previousChip=chipNo;
+ previousChip=chip;
}//end if
else{
- if(previousChip!=chipNo){
+ if(previousChip!=(Int_t)chip) {
WriteChipTrailer(buf,chipHitCount,baseWord);
chipHitCount=0;
- for(Int_t i=previousChip+1;i<chipNo;i++){
- WriteChipHeader(i,(mod/2),baseWord);
+ for(Int_t i=previousChip+1; i<chip; i++) {
+ WriteChipHeader(i,hs,baseWord);
WriteChipTrailer(buf,0,baseWord);
- chipHitCount=0;
}//end for
- WriteChipHeader(chipNo,(mod/2),baseWord);
- previousChip=chipNo;
+ WriteChipHeader(chip,hs,baseWord);
+ previousChip=chip;
}//end if
chipHitCount++;
- WriteHit(buf,ix,hitRow,baseWord);
+ WriteHit(buf,row,col,baseWord);
}//end else
}//end for
//Even if there aren't digits for a given chip
//the chip header and the chip trailer are stored
- Int_t end=4;
- if(chipNo>4)end+=5;
WriteChipTrailer(buf,chipHitCount,baseWord);
chipHitCount=0;
- for(Int_t i=chipNo+1;i<=end;i++){
- WriteChipHeader(i,(mod/2),baseWord);
+ for(Int_t i=chip+1;i<=chipHigh;i++){
+ WriteChipHeader(i,hs,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;
+ //the chip header and chip trailer are stored anyway
+ for(Int_t i=chipLow; i<=chipHigh; i++){
+ WriteChipHeader(i,hs,baseWord);
+ WriteChipTrailer(buf,0,baseWord);
}//end for
}//end else
+
+ delete digMap;
+
if(fVerbose==2)
ftxt.close();
return;
TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
char fileName[15];
- ofstream outfile; // logical name of the output file
- AliRawDataHeader header;
+ AliFstream* outfile; // logical name of the output file
+ AliRawDataHeaderSim header;
//loop over DDLs
for(Int_t i=0;i<AliDAQ::NumberOfDdls("ITSSPD");i++){
strcpy(fileName,AliDAQ::DdlFileName("ITSSPD",i)); //The name of the output file.
-#ifndef __DECCXX
- outfile.open(fileName,ios::binary);
-#else
- outfile.open(fileName);
-#endif
+ outfile = new AliFstream(fileName);
//write Dummy DATA HEADER
- UInt_t dataHeaderPosition=outfile.tellp();
- outfile.write((char*)(&header),sizeof(header));
+ UInt_t dataHeaderPosition=outfile->Tellp();
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
//Loops over Modules of a particular DDL
for (Int_t mod=0; mod<AliITSRawStreamSPD::kModulesPerDDL; mod++){
Int_t moduleNumber = AliITSRawStreamSPD::GetModuleNumber(i, mod);
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(digits,mod,i,buf);
- outfile.write((char *)buf,((fIndex+1)*sizeof(UInt_t)));
+ GetDigitsSPD(digits,moduleNumber,i,buf);
+ outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t)));
for(Int_t i=0;i<(fIndex+1);i++){
buf[i]=0;
}//end for
}//end for
//Write REAL DATA HEADER
- UInt_t currentFilePosition=outfile.tellp();
- outfile.seekp(dataHeaderPosition);
+ UInt_t currentFilePosition=outfile->Tellp();
+ outfile->Seekp(dataHeaderPosition);
header.fSize=currentFilePosition-dataHeaderPosition;
- outfile.write((char*)(&header),sizeof(header));
- outfile.close();
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
+ delete outfile;
}//end for
return 0;
Int_t AliITSDDLRawData::RawDataSSD(TBranch* branch){
- //This method creates the Raw data files for SSD detectors
+ //This method creates the Raw data files for SSD detectors
const Int_t kSize=1536;//768*2 Number of stripe * number of sides(N and P)
UInt_t buf[kSize];
fIndex=-1;
TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
char fileName[15];
- ofstream outfile; // logical name of the output file
- AliRawDataHeader header;
+ AliFstream* outfile; // logical name of the output file
+ AliRawDataHeaderSim header;
//loop over DDLs
for(Int_t i=0;i<AliDAQ::NumberOfDdls("ITSSSD");i++){
strcpy(fileName,AliDAQ::DdlFileName("ITSSSD",i)); //The name of the output file.
-#ifndef __DECCXX
- outfile.open(fileName,ios::binary);
-#else
- outfile.open(fileName);
-#endif
+ outfile = new AliFstream(fileName);
//write Dummy DATA HEADER
- UInt_t dataHeaderPosition=outfile.tellp();
- outfile.write((char*)(&header),sizeof(header));
+ UInt_t dataHeaderPosition=outfile->Tellp();
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
//Loops over Modules of a particular DDL
for (Int_t mod=0; mod<AliITSRawStreamSSD::kModulesPerDDL; mod++){
//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(digits,mod,moduleNumber,i,buf);
- outfile.write((char *)buf,((fIndex+1)*sizeof(UInt_t)));
+ outfile->WriteBuffer((char *)buf,((fIndex+1)*sizeof(UInt_t)));
fIndex=-1;
}//end if
}//end for
//Write REAL DATA HEADER
- UInt_t currentFilePosition=outfile.tellp();
- outfile.seekp(dataHeaderPosition);
+ UInt_t currentFilePosition=outfile->Tellp();
+ outfile->Seekp(dataHeaderPosition);
header.fSize=currentFilePosition-dataHeaderPosition;
header.SetAttribute(0); // valid data
- outfile.write((char*)(&header),sizeof(header));
- outfile.close();
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
+ delete outfile;
}//end for
return 0;
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
-Int_t AliITSDDLRawData::RawDataSDD(TBranch* branch){
+Int_t AliITSDDLRawData::RawDataSDD(TBranch* branch, AliITSDDLModuleMapSDD* ddlsdd){
//This method creates the Raw data files for SDD detectors
const Int_t kSize=131072; //256*512
UInt_t buf[kSize];
TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
char fileName[15];
- ofstream outfile; // logical name of the output file
- AliRawDataHeader header;
- UInt_t skippedword = AliBitPacking::PackWord(2,skippedword,0,31);
+ AliFstream* outfile; // logical name of the output file
+ AliRawDataHeaderSim header;
+ UInt_t skippedword, carlosFooterWord,fifoFooterWord,jitterWord;
+ Bool_t retcode;
+ retcode = AliBitPacking::PackWord(0x3FFFFFFF,carlosFooterWord,0,31);
+ retcode = AliBitPacking::PackWord(0x3F1F1F1F,fifoFooterWord,0,31);
+ retcode = AliBitPacking::PackWord(0x7F00000E,jitterWord,0,31);
//loop over DDLs
for(Int_t i=0;i<AliDAQ::NumberOfDdls("ITSSDD");i++){
strcpy(fileName,AliDAQ::DdlFileName("ITSSDD",i)); //The name of the output file.
-#ifndef __DECCXX
- outfile.open(fileName,ios::binary);
-#else
- outfile.open(fileName);
-#endif
+ outfile = new AliFstream(fileName);
//write Dummy DATA HEADER
- UInt_t dataHeaderPosition=outfile.tellp();
- outfile.write((char*)(&header),sizeof(header));
+ UInt_t dataHeaderPosition=outfile->Tellp();
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
+
+
+ //first 1 "dummy" word to be skipped
+ retcode = AliBitPacking::PackWord(0xFFFFFFFF,skippedword,0,31);
+ outfile->WriteBuffer((char*)(&skippedword),sizeof(skippedword));
- //first 9 "dummy" words to be skipped
- for(Int_t iw=0;iw<9;iw++){
- outfile.write((char*)&skippedword,sizeof(skippedword));
- }
-
//Loops over Modules of a particular DDL
for (Int_t mod=0; mod<AliITSRawStreamSDD::kModulesPerDDL; mod++){
- Int_t moduleNumber = AliITSRawStreamSDD::GetModuleNumber(i, mod);
+ Int_t moduleNumber = ddlsdd->GetModuleNumber(i, mod);
if(moduleNumber!=-1){
digits->Clear();
branch->GetEvent(moduleNumber);
//fIndex gives the number of 32 bits words in the buffer for each module
// cout<<"MODULE NUMBER:"<<mapSDD[i][mod]<<endl;
GetDigitsSDD(digits,mod,moduleNumber,i,buf);
- outfile.write((char *)buf,((fIndex+1)*sizeof(UInt_t)));
+ 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
+ // 12 words with FIFO footers (=4 FIFO x 3 3F1F1F1F words per DDL)
+ 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);
+ UInt_t currentFilePosition=outfile->Tellp();
+ outfile->Seekp(dataHeaderPosition);
header.fSize=currentFilePosition-dataHeaderPosition;
header.SetAttribute(0); // valid data
- outfile.write((char*)(&header),sizeof(header));
-
- outfile.close();
+ outfile->WriteBuffer((char*)(&header),sizeof(header));
+ delete outfile;
}//end for
return 0;