X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=TPC%2FAliTPCDDLRawData.cxx;h=dcec79dfcc00d142ba41362912fce55b02b6bcc2;hb=108cb8b981fd4f1bf0b1ec2b66c46daacb58697f;hp=5c16b16dec78c32bd3b52c71fb7a72db9fa49ca6;hpb=30c1018ed2b0e5ba76bf85a7a2d6695a50bf0f4a;p=u%2Fmrichter%2FAliRoot.git

diff --git a/TPC/AliTPCDDLRawData.cxx b/TPC/AliTPCDDLRawData.cxx
index 5c16b16dec7..dcec79dfcc0 100644
--- a/TPC/AliTPCDDLRawData.cxx
+++ b/TPC/AliTPCDDLRawData.cxx
@@ -14,52 +14,56 @@
  **************************************************************************/
 /* $Id$ */
 
+
 //This class conteins all the methods to create raw data 
 //as par a given DDL.
 //It produces DDL with both compressed and uncompressed format.
 //For compression we use the optimized table wich needs 
 //to be provided.
 
-
-#include "TObjArray.h"
-#include "Riostream.h"
+#include <TObjArray.h>
+#include <TString.h>
+#include <TSystem.h>
+#include <Riostream.h>
 #include <stdio.h>
 #include <stdlib.h>
-#include "AliTPCCompression.h"
-#include "AliTPCBuffer160.h"
+#include "AliAltroBuffer.h"
+#include "AliTPCAltroMapping.h"
 #include "AliTPCDDLRawData.h"
+#include "AliDAQ.h"
 
 ClassImp(AliTPCDDLRawData)
 ////////////////////////////////////////////////////////////////////////////////////////
 
-AliTPCDDLRawData::AliTPCDDLRawData(const AliTPCDDLRawData &source){
+AliTPCDDLRawData::AliTPCDDLRawData(const AliTPCDDLRawData &source):
+  TObject(source),
+  fVerbose(0)
+{
   // Copy Constructor
+  fVerbose=source.fVerbose;
   return;
 }
 
 AliTPCDDLRawData& AliTPCDDLRawData::operator=(const AliTPCDDLRawData &source){
   //Assigment operator
+  fVerbose=source.fVerbose;
   return *this;
 }
 
 
 ////////////////////////////////////////////////////////////////////////////
-void AliTPCDDLRawData::RawData(Int_t LDCsNumber){
-  //Raw data slides generation
+void AliTPCDDLRawData::RawData(const char* inputFileName){
+  //Raw data generation
   //Number of DDL=2*36+4*36=216
   //2 DDL for each inner sector
   //4 DDL for each outer sector
-  Int_t ddlPerFile=216/LDCsNumber;
-  Int_t offset=1;
-  if (216%LDCsNumber) ddlPerFile++;
-  cout<<"Number of DDL per slide: "<<ddlPerFile<<endl;
   ifstream f;
 #ifndef __DECCXX
-  f.open("AliTPCDDL.dat",ios::binary);
+  f.open(inputFileName,ios::binary);
 #else
-  f.open("AliTPCDDL.dat");
+  f.open(inputFileName);
 #endif
-  if(!f){cout<<"File doesn't exist !!"<<endl;return;}
+  if(!f){Error("RawData", "File doesn't exist !!");return;}
   struct DataPad{
     Int_t Sec;
     Int_t SubSec;
@@ -70,35 +74,51 @@ void AliTPCDDLRawData::RawData(Int_t LDCsNumber){
   };
   DataPad data;
 
-  //AliTPCBuffer160 is used in write mode to generate AltroFormat.dat file
-  Int_t sliceNumber=1;
+  //AliAltroBuffer is used in write mode to generate raw data file
   char  filename[15];
-  sprintf(filename,"TPCslice%d",sliceNumber); 
-  cout<<"   Creating "<<filename<<endl;
-  AliTPCBuffer160 *buffer=new AliTPCBuffer160(filename,1);
-
-  ULong_t count=0;
+  Int_t ddlNumber=0;
+  AliAltroBuffer *buffer=NULL;
   Int_t pSecNumber=-1;  //Previous Sector number
   Int_t pRowNumber=-1;  //Previous Row number  
   Int_t pPadNumber=-1;  //Previous Pad number
   Int_t pTimeBin=-1;    //Previous Time-Bin
   Int_t pSubSector=-1;  //Previous Sub Sector
   Int_t bunchLength=0;
-  Int_t countDDL=0;
   Int_t nwords=0;
-  ULong_t numPackets=0;
+  UInt_t numPackets=0;
+
+  TString path = gSystem->Getenv("ALICE_ROOT");
+  path += "/TPC/mapping/Patch";
+  TString path2;
+  AliTPCAltroMapping *mapping[6];
+  for(Int_t i = 0; i < 6; i++) {
+    path2 = path;
+    path2 += i;
+    path2 += ".data";
+    mapping[i] = new AliTPCAltroMapping(path2.Data());
+  }
+
+
   while (f.read((char*)(&data),sizeof(data))){
-    count++;
     if (pPadNumber==-1){
       pSecNumber=data.Sec;
       pRowNumber=data.Row;
       pPadNumber=data.Pad;
       pTimeBin=data.Time;
       pSubSector=data.SubSec;
+
+      if(data.Sec<36)
+	ddlNumber=data.Sec*2+data.SubSec;
+      else
+	ddlNumber=72+(data.Sec-36)*4+data.SubSec;
+      strcpy(filename,AliDAQ::DdlFileName("TPC",ddlNumber));
+      Int_t patchIndex = data.SubSec;
+      if(data.Sec>=36) patchIndex += 2;
+      buffer=new AliAltroBuffer(filename,mapping[patchIndex]);
       //size magic word sector number sub-sector number 0 for TPC 0 for uncompressed
-      buffer->WriteMiniHeader(0,pSecNumber,pSubSector,0,0);//Dummy;
+      buffer->WriteDataHeader(kTRUE,kFALSE);//Dummy;
       bunchLength=1;
-      buffer->FillBuffer(data.Dig-offset);
+      buffer->FillBuffer(data.Dig);
       nwords++;
     }//end if
     else{
@@ -119,25 +139,21 @@ void AliTPCDDLRawData::RawData(Int_t LDCsNumber){
 	  nwords=0;
 
 	  if(pSubSector!=data.SubSec){
-	    countDDL++;
-	    if(countDDL==ddlPerFile){
-	      //size magic word sector number sub-sector number 0 for TPC 0 for uncompressed
-	      buffer->Flush();
-	      buffer->WriteMiniHeader(1,pSecNumber,pSubSector,0,0);
-	      //cout<<"Mini header for DDL:"<<PSecNumber<<" Sub-sec:"<<PSubSector<<endl;
-	      delete buffer;
-	      sliceNumber++;
-	      sprintf(filename,"TPCslice%d",sliceNumber);
-	      cout<<"   Creating "<<filename<<endl;
-	      buffer=new AliTPCBuffer160(filename,1);
-	      buffer->WriteMiniHeader(0,data.Sec,data.SubSec,0,0);//Dummy;
-	      countDDL=0;
-	    }//end if
-	    else{
-	      buffer->Flush();
-	      buffer->WriteMiniHeader(1,pSecNumber,pSubSector,0,0);
-	      buffer->WriteMiniHeader(0,data.Sec,data.SubSec,0,0);//Dummy;
-	    }
+	    //size magic word sector number sub-sector number 0 for TPC 0 for uncompressed
+	    buffer->Flush();
+	    buffer->WriteDataHeader(kFALSE,kFALSE);
+	    //cout<<"Data header for DDL:"<<PSecNumber<<" Sub-sec:"<<PSubSector<<endl;
+	    delete buffer;
+
+	    if(data.Sec<36)
+	      ddlNumber=data.Sec*2+data.SubSec;
+	    else
+	      ddlNumber=72+(data.Sec-36)*4+data.SubSec;
+	    strcpy(filename,AliDAQ::DdlFileName("TPC",ddlNumber));
+	    Int_t patchIndex = data.SubSec;
+	    if(data.Sec>=36) patchIndex += 2;
+	    buffer=new AliAltroBuffer(filename,mapping[patchIndex]);
+	    buffer->WriteDataHeader(kTRUE,kFALSE);//Dummy;
 	    pSubSector=data.SubSec;
 	  }//end if
 	}//end if
@@ -148,262 +164,24 @@ void AliTPCDDLRawData::RawData(Int_t LDCsNumber){
 	pSecNumber=data.Sec;
       }//end else
       pTimeBin=data.Time;
-      buffer->FillBuffer(data.Dig-offset);
-      nwords++;
-    }//end else
-  }//end while
-  buffer->FillBuffer(pTimeBin);
-  buffer->FillBuffer(bunchLength+2);
-  nwords+=2;
-  buffer->WriteTrailer(nwords,pPadNumber,pRowNumber,pSecNumber);
-  //write the  M.H.
-  buffer->Flush();
-  buffer->WriteMiniHeader(1,pSecNumber,pSubSector,0,0);
-  //cout<<"Mini header for D D L:"<<pSecNumber<<" Sub-sec:"<<pSubSector<<endl;
-  delete buffer;
-  cout<<"Number of digits: "<<count<<endl;
-  f.close();
-  return;
-}
-////////////////////////////////////////////////////////////////////////////
-////////////////////////////////////////////////////////////////////////////
-
-
-Int_t AliTPCDDLRawData::RawDataCompDecompress(Int_t LDCsNumber,Int_t Comp){
-  //This method is used to compress and decompress the slides
-  static const Int_t kNumTables=5;
-  char filename[20];
-  char dest[20];
-  fstream f;
-  ULong_t size=0;
-  //Int_t MagicWord,DDLNumber,SecNumber,SubSector,Detector;
-  Int_t flag=0;
-  for(Int_t i=1;i<=LDCsNumber;i++){
-    if(!Comp){
-      sprintf(filename,"TPCslice%d",i);
-      sprintf(dest,"TPCslice%d.comp",i);
-    }
-    else{
-      sprintf(filename,"TPCslice%d.comp",i);
-      sprintf(dest,"TPCslice%d.decomp",i);
-    }
-#ifndef __DECCXX
-    f.open(filename,ios::binary|ios::in);
-#else
-    f.open(filename,ios::in);
-#endif
-    if(!f){cout<<"File doesn't exist \n";exit(1);}
-    cout<<filename<<"  "<<dest<<endl;
-    ofstream fdest;
-#ifndef __DECCXX
-    fdest.open(dest,ios::binary);
-#else
-    fdest.open(dest);
-#endif
-    //loop over the DDL block 
-    //Each block contains a Mini Header followed by raw data (ALTRO FORMAT)
-    //The number of block is ceil(216/LDCsNumber)
-    ULong_t miniHeader[3];
-    //here the Mini Header is read
-    while( (f.read((char*)(miniHeader),sizeof(ULong_t)*3)) ){
-      size=miniHeader[0];
-      //Int_t dim=sizeof(ULong_t)+sizeof(Int_t)*5;
-      //cout<<" Sec "<<SecNumber<<" SubSector "<<SubSector<<" size "<<size<<endl;
-      //open the temporay File
-      ofstream fo;
-      char temp[15]="TempFile";
-#ifndef __DECCXX
-      fo.open(temp,ios::binary);
-#else
-      fo.open(temp);
-#endif
-      Int_t car=0;
-      for(ULong_t j=0;j<size;j++){
-	f.read((char*)(&car),1);
-	fo.write((char*)(&car),1);
-      }//end for
-      fo.close();
-      //The temp file is compressed or decompressed
-      AliTPCCompression *util = new AliTPCCompression();
-      if(!Comp)
-	util->CompressDataOptTables(kNumTables,temp,"TempCompDecomp");
-      else
-	util->DecompressDataOptTables(kNumTables,temp,"TempCompDecomp");
-      delete util;
-      //the temp compressed file is open and copied to the final file fdest
-      ifstream fi;
-#ifndef __DECCXX
-      fi.open("TempCompDecomp",ios::binary);
-#else
-      fi.open("TempCompDecomp");
-#endif
-      fi.seekg(0,ios::end);
-      size=fi.tellg();
-      fi.seekg(0);
-      //The Mini Header is updated (size and Compressed flag) 
-      //and written into the output file
-      miniHeader[0]=size;
-      if(!Comp)
-	flag=1;
-      else
-	flag=0;
-      ULong_t aux=0xFFFF;
-      aux<<=16;
-      aux|=flag;
-      aux|=0xFF;
-      miniHeader[2]=miniHeader[2]&aux;
-      fdest.write((char*)(miniHeader),sizeof(ULong_t)*3);
-      //The compressem temp file is copied into the output file fdest
-      for(ULong_t j=0;j<size;j++){
-	fi.read((char*)(&car),1);
-	fdest.write((char*)(&car),1);
-      }//end for
-      fi.close();
-    }//end while
-    f.clear();
-    f.close();
-    fdest.close();
-    remove("TempFile");
-    remove("TempCompDecomp");
-  }//end for
-  return 0;
-}
-
-/////////////////////////////////////////////////////////////////////////////////
-void AliTPCDDLRawData::RawDataAltro()const{
-  //This method is used to build the Altro format from AliTPCDDL.dat
-  //It is used to debug the code and create the tables used in the compresseion phase
-  Int_t offset=1;
-  ifstream f;
-#ifndef __DECCXX
-  f.open("AliTPCDDL.dat",ios::binary);
-#else
-  f.open("AliTPCDDL.dat");
-#endif
-  if(!f){cout<<"File doesn't exist !!"<<endl;return;}
-  struct DataPad{
-    Int_t Sec;
-    Int_t SubSec;
-    Int_t Row;
-    Int_t Pad;
-    Int_t Dig;
-    Int_t Time;
-  };
-  DataPad data;
-
-  //AliTPCBuffer160 is used in write mode to generate AltroFormat.dat file
-  char  filename[30]="AltroFormatDDL.dat";
-  cout<<"   Creating "<<filename<<endl;
-  AliTPCBuffer160 *buffer=new AliTPCBuffer160(filename,1);
-
-  ULong_t count=0;
-  Int_t pSecNumber=-1;  //Previous Sector number
-  Int_t pRowNumber=-1;  //Previous Row number  
-  Int_t pPadNumber=-1;  //Previous Pad number
-  Int_t pTimeBin=-1;    //Previous Time-Bin
-  Int_t bunchLength=0;
-  Int_t nwords=0;
-  ULong_t numPackets=0;
-  while (f.read((char*)(&data),sizeof(data))){
-    count++;
-    if (pPadNumber==-1){
-      pSecNumber=data.Sec;
-      pRowNumber=data.Row;
-      pPadNumber=data.Pad;
-      pTimeBin=data.Time;
-      bunchLength=1;
-      buffer->FillBuffer(data.Dig-offset);
-      nwords++;
-    }//end if
-    else{
-      if ( (data.Time==(pTimeBin+1)) &&
-	   (pPadNumber==data.Pad) &&
-	   (pRowNumber==data.Row) &&
-	   (pSecNumber==data.Sec)){
-	bunchLength++;
-      }//end if
-      else{
-	buffer->FillBuffer(pTimeBin);
-	buffer->FillBuffer(bunchLength+2);
-	nwords+=2;
-	if ((pPadNumber!=data.Pad)||(pRowNumber!=data.Row)||(pSecNumber!=data.Sec)){
-	  //Trailer is formatted and inserted!!
-	  buffer->WriteTrailer(nwords,pPadNumber,pRowNumber,pSecNumber);
-	  numPackets++;
-	  nwords=0;
-	}//end if
-	
-	bunchLength=1;
-	pPadNumber=data.Pad;
-	pRowNumber=data.Row;
-	pSecNumber=data.Sec;
-      }//end else
-      pTimeBin=data.Time;
-      buffer->FillBuffer(data.Dig-offset);
+      buffer->FillBuffer(data.Dig);
       nwords++;
     }//end else
   }//end while
-  buffer->FillBuffer(pTimeBin);
-  buffer->FillBuffer(bunchLength+2);
-  nwords+=2;
-  buffer->WriteTrailer(nwords,pPadNumber,pRowNumber,pSecNumber);
-  delete buffer;
-  cout<<"Number of digits: "<<count<<endl;
-  f.close(); 
-  return;
-}
+  if (buffer) {
+    buffer->FillBuffer(pTimeBin);
+    buffer->FillBuffer(bunchLength+2);
+    nwords+=2;
+    buffer->WriteTrailer(nwords,pPadNumber,pRowNumber,pSecNumber);
+    //write the  D.H.
+    buffer->Flush();
+    buffer->WriteDataHeader(kFALSE,kFALSE);
+    //cout<<"Data header for D D L:"<<pSecNumber<<" Sub-sec:"<<pSubSector<<endl;
+    delete buffer;
+  }
+
+  for(Int_t i = 0; i < 6; i++) delete mapping[i];
 
-/////////////////////////////////////////////////////////////////////////
-void AliTPCDDLRawData::RawDataAltroDecode(Int_t LDCsNumber,Int_t Comp){
-  //This method merges the slides in only one file removing at the same 
-  //time all the mini headers. The file so obtained must be Altro format
-  //complaiant.
-  //It is used mainly in the debugging phase 
-  char filename[15];
-  char dest[30];
-  fstream f;
-  if(!Comp)
-    sprintf(dest,"AltroDDLRecomposed.dat");
-  else
-    sprintf(dest,"AltroDDLRecomposedDec.dat");
-  ofstream fdest;
-
-#ifndef __DECCXX
-  fdest.open(dest,ios::binary);
-#else
-  fdest.open(dest);
-#endif
-  ULong_t size=0;
-  //Int_t MagicWord,DDLNumber,SecNumber,SubSector,Detector,flag=0;
-  for(Int_t i=1;i<=LDCsNumber;i++){
-    if(!Comp)
-      sprintf(filename,"TPCslice%d",i);  
-    else
-      sprintf(filename,"TPCslice%d.decomp",i);  
-#ifndef __DECCXX
-    f.open(filename,ios::binary|ios::in);
-#else
-    f.open(filename,ios::in);
-#endif
-    if(!f){cout<<"The file doesn't exist"<<endl;exit(1);}
-    //loop over the DDL block 
-    //Each block contains a Mini Header followed by raw data (ALTRO FORMAT)
-    //The number of block is ceil(216/LDCsNumber)
-    ULong_t miniHeader[3];
-    //here the Mini Header is read
-    while( (f.read((char*)(miniHeader),sizeof(ULong_t)*3)) ){
-      //cout<<"Mini header dimension "<<miniHeader[0]<<endl;
-      Int_t car=0;
-      size=miniHeader[0];
-      for(ULong_t j=0;j<size;j++){
-	f.read((char*)(&car),1);
-	fdest.write((char*)(&car),1);
-      }//end for
-    }//end while
-    f.clear();
-    f.close();
-  }//end for
-  fdest.close();
+  f.close();
   return;
 }
-