bug fixed
[u/mrichter/AliRoot.git] / PMD / AliPMDDDLRawData.cxx
index b752b39..faca049 100644 (file)
 #include <TTree.h>
 #include <TBranch.h>
 #include <TMath.h>
+#include <TString.h>
+#include <TSystem.h>
 
+#include "AliLog.h"
+#include "AliRawDataHeaderSim.h"
+#include "AliBitPacking.h"
 #include "AliPMDdigit.h"
+#include "AliPMDBlockHeader.h"
+#include "AliPMDDspHeader.h"
+#include "AliPMDPatchBusHeader.h"
+#include "AliPMDRawStream.h"
+#include "AliPMDddlinfoData.h"
+#include "AliPMDMappingData.h"
 #include "AliPMDDDLRawData.h"
-
+#include "AliDAQ.h"
+#include "AliFstream.h"
+#include "AliCDBManager.h"
+#include "AliCDBStorage.h"
+#include "AliCDBEntry.h"
 
 ClassImp(AliPMDDDLRawData)
 
 AliPMDDDLRawData::AliPMDDDLRawData():
+  fDdlinfo(GetDdlinfoData()),
+  fMapData(GetMappingData()),
   fDigits(new TClonesArray("AliPMDdigit", 1000))
 {
   // Default Constructor
   //
+
+}
+//____________________________________________________________________________
+AliPMDDDLRawData::AliPMDDDLRawData(const AliPMDDDLRawData& ddlraw):
+  TObject(ddlraw),
+  fDdlinfo(ddlraw.fDdlinfo),
+  fMapData(ddlraw.fMapData),
+  fDigits(ddlraw.fDigits)
+{
+  //Copy Constructor 
+}
+//____________________________________________________________________________
+AliPMDDDLRawData & AliPMDDDLRawData::operator=(const AliPMDDDLRawData& ddlraw)
+{
+  //Assignment operator 
+  if(this != &ddlraw)
+    {
+      fDdlinfo = ddlraw.fDdlinfo;
+      fMapData = ddlraw.fMapData;
+      fDigits  = ddlraw.fDigits;
+    }
+  return *this;
 }
 //____________________________________________________________________________
 
@@ -41,437 +80,778 @@ AliPMDDDLRawData::~AliPMDDDLRawData()
 }
 
 //____________________________________________________________________________
-void AliPMDDDLRawData::WritePMDRawData(TTree *treeD, Int_t evtno)
+void AliPMDDDLRawData::WritePMDRawData(TTree *treeD)
 {
   // write digits into raw data format
 
-  ofstream outfile;
+  AliFstream *outfile;
 
   TBranch *branch = treeD->GetBranch("PMDDigit");
-  if (!branch) return;
+  if (!branch)
+    {
+      AliError("PMD Digit branch not found");
+      return;
+    }
   branch->SetAddress(&fDigits);  
   
   Int_t   nmodules = (Int_t) treeD->GetEntries();
-  //  cout << " nmodules = " << nmodules << endl;
+  AliDebug(1,Form("Number of modules inside treeD = %d",nmodules));
 
-  const Int_t kSize         = 4608;
-  const Int_t kDDL          = 6;
-  Int_t modulePerDDL        = 0;
+  const Int_t kDDL          = AliDAQ::NumberOfDdls("PMD");
 
 
+  AliRawDataHeaderSim header;
   UInt_t sizeRawData = 0;
-  UInt_t magicWord   = 0x123456;
-  UInt_t detectorID  = 0;
-  UInt_t ddlID       = 0;
-  Int_t  flag        = 0;
-  Int_t  version     = 1;
-  UInt_t mHSize      = 3*sizeof(UInt_t);
-  Int_t  totword     = 0;
-
-  UInt_t miniHeader[3];
+  
+  const Int_t kbusSize = 51;
+  const Int_t kSize = 1536;
   UInt_t buffer[kSize];
 
+  UInt_t busPatch[kbusSize][1536];
+
+  Int_t contentsBus[kbusSize];
+
   Char_t filename[80];
 
+  Int_t modulePerDDL        = 0;
   Int_t mmodule = 0;
+  Int_t ddlno;
+  Int_t modulenoddl[12] = {-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1};
+
   for(Int_t iddl = 0; iddl < kDDL; iddl++)
     {
-      sprintf(filename,"Ev%dPMDddl%d.dat",evtno,iddl);
-#ifndef __DECCXX
-      outfile.open(filename,ios::binary);
-#else
-      outfile.open(filename);
-#endif
-      
-      if (iddl < 4)
-       {
-         modulePerDDL = 6;
-         mmodule = 6*iddl;
-         detectorID = 0;
-       }
-      else if (iddl == 4)
+      ddlno = iddl;
+      modulePerDDL = fDdlinfo->GetNoOfModulePerDdl(iddl);
+      if (modulePerDDL == 0) continue;
+      for (Int_t im = 0; im < 12; im++)
        {
-         modulePerDDL = 12;
-         mmodule = 24;
-         detectorID = 1;
+         modulenoddl[im] = fDdlinfo->GetModulesPerDdl(iddl,im);;
        }
-      else if (iddl == 5)
-       {
-         modulePerDDL = 12;
-         mmodule = 36;
-         detectorID = 1;
-       }
-
-      miniHeader[0] = sizeRawData;
-      PackWord(0,23,magicWord,miniHeader[1]);
-      PackWord(24,31,detectorID,miniHeader[1]);
-      ddlID = iddl;
-      PackWord(0,15,ddlID,miniHeader[2]);
-      PackWord(16,23,flag,miniHeader[2]);
-      PackWord(24,31,version,miniHeader[2]);
-
-
-      // Write the Dummy Mini Header into the file
-      Int_t bHPosition = outfile.tellp();
-      outfile.write((char*)(miniHeader),mHSize);
 
+      strncpy(filename,AliDAQ::DdlFileName("PMD",iddl),79);
+      
+      outfile = new AliFstream(filename);
+      
+      // Write the Dummy Data Header into the file
+      Int_t bHPosition = outfile->Tellp();
+      outfile->WriteBuffer((char*)(&header),sizeof(header));
 
-      for(Int_t ium = 0; ium < modulePerDDL; ium++)
+      for (Int_t ibus = 0; ibus < kbusSize; ibus++)
        {
-         
-         for (Int_t i = 0; i < kSize; i++)
+         contentsBus[ibus] = 0;
+         for (Int_t ich = 0; ich < kSize; ich++)
            {
-             buffer[i] = 0;
+             busPatch[ibus][ich] = 0;
            }
+       }
+
+      for(Int_t ium = 0; ium < 12; ium++)
+       {
          // Extract energy deposition per cell and pack it
          // in a 32-bit word and returns all the total words
          // per one unit-module
          
-         GetUMDigitsData(treeD, mmodule, ium, iddl, totword, buffer);
-         
-         outfile.write((char*)buffer,totword*sizeof(UInt_t));
-
-         mmodule++;
+         mmodule = modulenoddl[ium];
+         if(mmodule == -1) continue;
+         GetUMDigitsData(treeD, mmodule, iddl, contentsBus, busPatch);
+       }
 
+      Int_t ij = 0;
+      Int_t dsp[10];
+      Int_t dspBus[10];
+      for (Int_t i = 0; i < 10; i++)
+       {
+         dsp[i] = 0;
+         dspBus[i] = 0;
+         for (Int_t ibus=0; ibus < 5; ibus++)
+           {
+             ij++;
+             if (contentsBus[ij] > 0)
+               {
+                 dsp[i] += contentsBus[ij];
+               }
+             dspBus[i]++;
+           }
+         // Add the patch Bus header to the DSP contents
+         dsp[i] += 4*dspBus[i];
        }
 
-      // Write real mini header
-      // take the pointer to the beginning of the mini header
-      // write the total number of words per ddl and bring the
-      // pointer to the current file position and close it
-      UInt_t cFPosition = outfile.tellp();
-      sizeRawData = cFPosition - bHPosition - mHSize;
-      outfile.seekp(bHPosition);
-      outfile.write((char*)(&sizeRawData),sizeof(UInt_t));
-      outfile.seekp(cFPosition);
+      Int_t dspBlockARDL    = 0;
+      Int_t dspBlockBRDL    = 0;
+      Int_t remainder       = 0;
 
-      outfile.close();
-    } // DDL Loop over
 
+      for (Int_t i = 0; i < 5; i++)
+       {
+         Int_t ieven = 2*i;
+         Int_t iodd  = 2*i + 1;
+         if (dsp[ieven] > 0)
+           {
+             dspBlockARDL += dsp[ieven];
+             remainder = dsp[ieven]%2;
+             if (remainder == 1)
+               {
+                 dspBlockARDL++;
+               }
+           }
+         if (dsp[iodd] > 0)
+           {
+             dspBlockBRDL += dsp[iodd];
+             remainder = dsp[iodd]%2;
+             if (remainder == 1)
+               {
+                 dspBlockBRDL++;
+               }
+           }
+       }
+
+      dspBlockARDL += 50;
+      dspBlockBRDL += 50;
 
-}
-//____________________________________________________________________________
+      // Start writing the DDL file
 
-void AliPMDDDLRawData::ReadPMDRawData(Int_t evtno)
-{
-  // reads the raw data
-  ifstream infile;
+      AliPMDBlockHeader blHeader;
+      AliPMDDspHeader   dspHeader;
+      AliPMDPatchBusHeader pbusHeader;
 
-  Char_t filename[80];
+      const Int_t kblHLen   = blHeader.GetHeaderLength();
+      const Int_t kdspHLen  = dspHeader.GetHeaderLength();
+      const Int_t kpbusHLen = pbusHeader.GetHeaderLength();
 
-  const Int_t kDDL = 6;
-  Int_t imodule = 0;
+      UInt_t dspRDL = 0;
+      UInt_t dspBlockHeaderWord[8];
+      UInt_t dspHeaderWord[10];
+      UInt_t patchBusHeaderWord[4];
+      Int_t  iskip[5];
+      UInt_t ddlEndWord[2] = {0xDEADFACE, 0xDEADFACE};
 
-  for(Int_t iddl = 0; iddl < kDDL; iddl++)
-    {
-      sprintf(filename,"Ev%dPMDddl%d.dat",evtno,iddl);
-#ifndef __DECCXX
-      infile.open(filename,ios::binary);
-#else
-      infile.open(filename);
-#endif
-      
-      Int_t  ium;
-      Int_t  irow;
-      Int_t  icol;
-      UInt_t baseword;
-      UInt_t miniHeader[3];
-      UInt_t sizeRawData;
-      UInt_t mHSize = 3*sizeof(UInt_t);
-
-      infile.read((char*)(miniHeader), mHSize);
-      
-      sizeRawData = miniHeader[0];
-      Int_t totword = sizeRawData/4;
-      
-      UInt_t adc;
-      UInt_t chno;
-      UInt_t mcmno;
-      
-      for (Int_t iword = 0; iword < totword; iword++)
+      for (Int_t iblock = 0; iblock < 2; iblock++)
        {
-         infile.read((char*)(&baseword),sizeof(UInt_t));
-         
-         UnpackWord(0,11,adc,baseword);
-         UnpackWord(12,17,chno,baseword);
-         UnpackWord(18,28,mcmno,baseword);
+         // DSP Block Header
          
-         GetRowCol(iddl, mcmno, chno, ium, irow, icol);
-         if (iddl < 4)
+         for (Int_t i=0; i<kblHLen; i++)
            {
-             imodule = iddl*6 + ium;
+             dspBlockHeaderWord[i] = 0;
            }
-         else if (iddl == 4)
+         if (iblock == 0)
            {
-             imodule = 24 + ium;
+             dspBlockHeaderWord[1] = (UInt_t) (dspBlockARDL + kblHLen);
+             dspBlockHeaderWord[2] = (UInt_t) dspBlockARDL;
            }
-         else if (iddl == 5)
+         else if (iblock == 1)
            {
-             imodule = 36 + ium;
+             dspBlockHeaderWord[1] = (UInt_t) (dspBlockBRDL + kblHLen);
+             dspBlockHeaderWord[2] = (UInt_t) dspBlockBRDL;
            }
 
-         fprintf(fpw3,"%d %d %d %d %d %d\n",imodule,irow,icol,mcmno,chno,adc);
-         
+         outfile->WriteBuffer((char*)dspBlockHeaderWord,kblHLen*sizeof(UInt_t));
+
+         if (iblock == 0)
+           {
+             iskip[0] = 0;
+             iskip[1] = 10;
+             iskip[2] = 20;
+             iskip[3] = 30;
+             iskip[4] = 40;
+           }
+         else if (iblock == 1)
+           {
+             iskip[0] = 5;
+             iskip[1] = 15;
+             iskip[2] = 25;
+             iskip[3] = 35;
+             iskip[4] = 45;
+           }
+
+         for (Int_t idsp = 0; idsp < 5; idsp++)
+           {
+             // DSP Header
+             Int_t dspno = 0;
+             if (iblock == 0)
+               {
+                 dspno = 2*idsp;
+                 dspRDL = (UInt_t) dsp[dspno];
+               }
+             else if (iblock == 1)
+               {
+                 dspno = 2*idsp + 1;
+                 dspRDL = (UInt_t) dsp[dspno];
+               }
+
+             for (Int_t i=0; i<kdspHLen; i++)
+               {
+                 dspHeaderWord[i] = 0;
+               }
+             remainder = dspRDL%2;
+             if (remainder == 1) dspRDL++;
+
+             dspHeaderWord[1] = dspRDL + kdspHLen;
+             dspHeaderWord[2] = dspRDL;
+             dspHeaderWord[3] = dspno;
+             if (remainder == 1) dspHeaderWord[8] = 1; // setting the padding word
+
+
+             outfile->WriteBuffer((char*)dspHeaderWord,kdspHLen*sizeof(UInt_t));
+
+             for (Int_t ibus = 0; ibus < 5; ibus++)
+               {
+                 // Patch Bus Header
+
+                 Int_t busno = iskip[idsp] + ibus + 1;
+                 Int_t patchbusRDL = contentsBus[busno];
+
+                 if (patchbusRDL > 0)
+                   {
+                     patchBusHeaderWord[0] = 0;
+                     patchBusHeaderWord[1] = (UInt_t) (patchbusRDL + kpbusHLen);
+                     patchBusHeaderWord[2] = (UInt_t) patchbusRDL;
+                     patchBusHeaderWord[3] = (UInt_t) busno;
+                   }
+                 else if (patchbusRDL == 0)
+                   {
+                     patchBusHeaderWord[0] = 0;
+                     patchBusHeaderWord[1] = (UInt_t) kpbusHLen;
+                     patchBusHeaderWord[2] = (UInt_t) 0;
+                     patchBusHeaderWord[3] = (UInt_t) busno;
+                   }
+
+
+                 outfile->WriteBuffer((char*)patchBusHeaderWord,4*sizeof(UInt_t));
+
+                 for (Int_t iword = 0; iword < patchbusRDL; iword++)
+                   {
+                     buffer[iword] = busPatch[busno][iword];
+                   }
+                 
+                 outfile->WriteBuffer((char*)buffer,patchbusRDL*sizeof(UInt_t));
+
+               } // End of patch bus loop
+
+
+             // Adding a padding word if the total words odd
+             if (remainder == 1)
+               {
+                 UInt_t paddingWord = dspHeader.GetDefaultPaddingWord();
+                 outfile->WriteBuffer((char*)(&paddingWord),sizeof(UInt_t));
+               }
+           }
        }
 
-    } // ddl loop
+      // Write two extra word at the end of each DDL file
+      outfile->WriteBuffer((char*)ddlEndWord,2*sizeof(UInt_t));
+
+      // Write real data header
+      // take the pointer to the beginning of the data header
+      // write the total number of words per ddl and bring the
+      // pointer to the current file position and close it
+      UInt_t cFPosition = outfile->Tellp();
+      sizeRawData = cFPosition - bHPosition - sizeof(header);
+
+      header.fSize = cFPosition - bHPosition;
+      header.SetAttribute(0);  // valid data
+      outfile->Seekp(bHPosition);
+      outfile->WriteBuffer((char*)(&header),sizeof(header));
+      outfile->Seekp(cFPosition);
+
+      delete outfile;
+    } // DDL Loop over
+
 }
 //____________________________________________________________________________
-void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule, Int_t ium,
-                                      Int_t ddlno, Int_t & totword,
-                                      UInt_t *buffer)
+
+void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule,
+                                      Int_t ddlno,  Int_t *contentsBus,
+                                      UInt_t busPatch[][1536])
 {
+  // Retrieves digits data UnitModule by UnitModule
 
-  UInt_t dataword, baseword;
+  const Int_t kMaxBus = 51;
 
-  UInt_t mcmno, chno;
-  UInt_t adc;
-  Int_t  irownew = 0;
-  Int_t  icolnew = 0;
-  Int_t  det, smn, irow, icol;
+  UInt_t baseword = 0;
+  UInt_t mcmno = 0, chno = 0;
+  UInt_t adc = 0;
+  Int_t  det = 0, smn = 0, irow = 0, icol = 0;
+  UInt_t  parity = 0;
+  
+  Int_t busno = 0;
+  Int_t patchBusNo[kMaxBus], mcmperBus[kMaxBus];
+  Int_t startRowBus[kMaxBus], startColBus[kMaxBus];
+  Int_t endRowBus[kMaxBus], endColBus[kMaxBus];
+
+  Int_t beginPatchBus = -1;
+  Int_t endPatchBus   = -1;
+  for(Int_t i = 0; i < kMaxBus; i++)
+    {
+      patchBusNo[i]  = -1;
+      mcmperBus[i]   = -1;
+      startRowBus[i] = -1;
+      startColBus[i] = -1;
+      endRowBus[i]   = -1;
+      endColBus[i]   = -1;
+    }
+
+  // Fetch the DDL mapping info from the mapping database
 
+  DdlMapping(ddlno, imodule, beginPatchBus, endPatchBus,
+            patchBusNo, mcmperBus, startRowBus, endRowBus,
+            startColBus, endColBus);
+
+  // Read if some chains are off from the ddlinfo database
+
+  Int_t srowoff1[2][24], erowoff1[2][24];
+  Int_t scoloff1[2][24], ecoloff1[2][24];
+  Int_t srowoff2[2][24], erowoff2[2][24];
+  Int_t scoloff2[2][24], ecoloff2[2][24];
+
+  for (Int_t idet = 0; idet < 2; idet++)
+    {
+      for (Int_t im = 0; im < 24; im++)
+       {
+         srowoff1[idet][im] = fDdlinfo->GetStartRowA(idet,im);
+         erowoff1[idet][im] = fDdlinfo->GetEndRowA(idet,im);
+         scoloff1[idet][im] = fDdlinfo->GetStartColA(idet,im);
+         ecoloff1[idet][im] = fDdlinfo->GetEndColA(idet,im);
+         srowoff2[idet][im] = fDdlinfo->GetStartRowB(idet,im);
+         erowoff2[idet][im] = fDdlinfo->GetEndRowB(idet,im);
+         scoloff2[idet][im] = fDdlinfo->GetStartColB(idet,im);
+         ecoloff2[idet][im] = fDdlinfo->GetEndColB(idet,im);
+       }
+    }
 
   treeD->GetEntry(imodule); 
   Int_t nentries = fDigits->GetLast();
-  totword = nentries+1;
+  Int_t totword = nentries+1;
+
+  AliPMDdigit *pmddigit = 0x0;
 
   for (Int_t ient = 0; ient < totword; ient++)
     {
-      fPMDdigit = (AliPMDdigit*)fDigits->UncheckedAt(ient);
+      pmddigit = (AliPMDdigit*)fDigits->UncheckedAt(ient);
       
-      det    = fPMDdigit->GetDetector();
-      smn    = fPMDdigit->GetSMNumber();
-      irow   = fPMDdigit->GetRow();
-      icol   = fPMDdigit->GetColumn();
-      adc    = (UInt_t) fPMDdigit->GetADC();
+      det    = pmddigit->GetDetector();
+      smn    = pmddigit->GetSMNumber();
+      irow   = pmddigit->GetRow();
+      icol   = pmddigit->GetColumn();
+      Float_t aadc = pmddigit->GetADC();
+      if (aadc < 0.) aadc = 0.;
+      adc    = (UInt_t) aadc;
+
+      TransformS2H(smn,irow,icol);
+      
+      // remove the non-existence channels
+
+      //printf("%d %d %d %d\n",det,smn,irow,icol);
+      //printf("--- %d   %d   %d   %d\n",srowoff[det][smn],erowoff[det][smn],
+      //     scoloff[det][smn],ecoloff[det][smn]);
 
-      if(smn < 12)
+      if (irow >= srowoff1[det][smn] && irow <= erowoff1[det][smn])
        {
-         irownew = icol;
-         icolnew = irow;
+         if (icol >= scoloff1[det][smn] && icol <= ecoloff1[det][smn])
+           {
+             continue;
+           }
        }
-      else if( smn >= 12 && smn < 24)
+      if (irow >= srowoff2[det][smn] && irow <= erowoff2[det][smn])
        {
-         irownew = irow;
-         icolnew = icol;
+         if (icol >= scoloff2[det][smn] && icol <= ecoloff2[det][smn])
+           {
+             continue;
+           }
        }
 
 
-
-      GetMCMCh(ddlno, ium, irownew, icolnew, mcmno, chno);
-
-      fprintf(fpw2,"%d %d %d %d %d %d\n",imodule,irownew,icolnew,mcmno,chno,adc);
+      GetMCMCh(imodule, irow, icol, beginPatchBus, endPatchBus,
+              mcmperBus, startRowBus, startColBus,
+              endRowBus, endColBus, busno, mcmno, chno);
 
       baseword = 0;
-      dataword = adc;
-      PackWord(0, 11,dataword,baseword);
-      dataword = chno;
-      PackWord(12,17,dataword,baseword);
-      dataword = mcmno;
-      PackWord(18,28,dataword,baseword);
-      dataword = 0;
-      PackWord(29,29,dataword,baseword);
-      dataword = 0;
-      PackWord(30,30,dataword,baseword);
-      dataword = 1;
-      PackWord(31,31,dataword,baseword);
-      
-      buffer[ient] = baseword;
-      
+      AliBitPacking::PackWord(adc,baseword,0,11);
+      AliBitPacking::PackWord(chno,baseword,12,17);
+      AliBitPacking::PackWord(mcmno,baseword,18,28);
+      AliBitPacking::PackWord(0,baseword,29,30);
+      parity = ComputeParity(baseword);      // generate the parity bit
+      AliBitPacking::PackWord(parity,baseword,31,31);
 
+      Int_t jj = contentsBus[busno];
+      busPatch[busno][jj] = baseword;
 
-      //      fprintf(fpw1,"%d %d %d %d %d\n",det,smn,irow,icol,adc);
-      
+      contentsBus[busno]++;
     }
 
-
 }
+
 //____________________________________________________________________________
-void AliPMDDDLRawData::GetMCMCh(Int_t ddlno, Int_t um,
-                               Int_t row, Int_t col,
-                               UInt_t &mcmno, UInt_t &chno)
+void AliPMDDDLRawData::TransformS2H(Int_t smn, Int_t &irow, Int_t &icol)
 {
-  // This part will be modified once the final track layout on the PCB is
-  // designed. This will only change the coordinate of the individual cell
-  UInt_t ch[16][4] = { {3, 2, 29, 28},
-                      {5, 1, 31, 27},
-                      {6, 0, 30, 26},
-                      {7, 4, 25, 24},
-                      {8, 9, 20, 23},
-                      {10, 14, 16, 22},
-                      {11, 15, 17, 21},
-                      {12, 13, 18, 19},
-                      {35, 34, 61, 60},
-                      {37, 33, 63, 59},
-                      {38, 32, 62, 58},
-                      {39, 36, 57, 56},
-                      {40, 41, 52, 55},
-                      {42, 46, 48, 54},
-                      {43, 47, 49, 53},
-                      {44, 45, 50, 51} };
-  
-  if (ddlno == 0 || ddlno == 1)
-    {
-      // PRE plane, SU Mod = 1, 2
-      Int_t irownew = row%16;
-      Int_t icolnew = col%4;
-      Int_t irowdiv = row/16;
-      Int_t icoldiv = col/4;
-
-      mcmno = 72*um + 6*icoldiv + irowdiv;
-      chno  = ch[irownew][icolnew];
-    }
-  else if (ddlno == 2 || ddlno == 3)
-    {
-      // PRE plane,  SU Mod = 3, 4
-      Int_t irownew = row%16;
-      Int_t icolnew = col%4;
-      Int_t irowdiv = row/16;
-      Int_t icoldiv = col/4;
-      
-      mcmno = 72*um + 3*icoldiv + irowdiv;
-      chno  = ch[irownew][icolnew];
-    }
-  else if (ddlno == 4)
+  // Does the Software to Hardware coordinate transformation
+  //
+
+  Int_t  irownew = 0;
+  Int_t  icolnew = 0;
+
+  // First in digits we have all dimension 48x96
+  // Transform into the realistic one, i.e, For SM 0&1 96(row)x48(col)
+  // and for SM 2&3 48(row)x96(col)
+  // 
+  if(smn < 12)
     {
-      // CPV plane,  SU Mod = 1, 2
-      Int_t irownew = row%16;
-      Int_t icolnew = col%4;
-      Int_t irowdiv = row/16;
-      Int_t icoldiv = col/4;
-
-      mcmno = 72*um + 6*icoldiv + irowdiv;
-      chno  = ch[irownew][icolnew];
+      irownew = icol;
+      icolnew = irow;
     }
-  else if (ddlno == 5)
+  else if( smn >= 12 && smn < 24)
     {
-      // CPV plane,  SU Mod = 3, 4
-      Int_t irownew = row%16;
-      Int_t icolnew = col%4;
-      Int_t irowdiv = row/16;
-      Int_t icoldiv = col/4;
-      
-      mcmno = 72*um + 3*icoldiv + irowdiv;
-      chno  = ch[irownew][icolnew];
+      irownew = irow;
+      icolnew = icol;
     }
 
+  irow = irownew;
+  icol = icolnew;
+
 }
+
+
 //____________________________________________________________________________
 
-void AliPMDDDLRawData::GetRowCol(Int_t ddlno, UInt_t mcmno, UInt_t chno,
-                                Int_t &um, Int_t &row, Int_t &col)
+void AliPMDDDLRawData::GetMCMCh(Int_t imodule, Int_t row, Int_t col,
+                               Int_t beginPatchBus, Int_t endPatchBus,
+                               Int_t *mcmperBus,
+                               Int_t *startRowBus, Int_t *startColBus,
+                               Int_t *endRowBus, Int_t *endColBus,
+                               Int_t & busno, UInt_t &mcmno, UInt_t &chno)
 {
-  UInt_t ch[64] = { 9, 5, 1, 0, 13, 4, 8, 12,
-                   16, 17, 20, 24, 28, 29, 21, 25,
-                   22, 26, 30, 31, 18, 27, 23, 19,
-                   15, 14, 11, 7, 3, 2, 10, 6,
-                   41, 37, 33, 32, 45, 36, 40, 44,
-                   48, 49, 52, 56, 60, 61, 53, 57,
-                   54, 58, 62, 63, 50, 59, 55, 51,
-                   47, 46, 43, 39, 35, 34, 42, 38 };
-
-  if (ddlno == 1 || ddlno == 2)
-    {
-      um  = mcmno/72;
-      Int_t mcmnonew = mcmno - 72*um;
-      Int_t rowcol  = ch[chno];
-      Int_t irownew = rowcol/4;
-      Int_t icolnew = rowcol%4;
-      
-      Int_t remmcm  = mcmnonew%6;
-      Int_t divmcm  = mcmnonew/6;
-      
-      row = 16*remmcm + irownew;
-      col =  4*divmcm + icolnew;
-    }
-  else   if (ddlno == 3 || ddlno == 4)
-    {
-      um  = mcmno/72;
-      Int_t mcmnonew = mcmno - 72*um;
-      Int_t rowcol  = ch[chno];
-      Int_t irownew = rowcol/4;
-      Int_t icolnew = rowcol%4;
-      
-      Int_t remmcm  = mcmnonew%3;
-      Int_t divmcm  = mcmnonew/3;
-      
-      row = 16*remmcm + irownew;
-      col =  4*divmcm + icolnew;
-    }
-  else if (ddlno == 4)
-    {
-      um  = mcmno/144;
-      Int_t mcmnonew = mcmno - 72*um;
-      Int_t rowcol  = ch[chno];
-      Int_t irownew = rowcol/4;
-      Int_t icolnew = rowcol%4;
-      
-      Int_t remmcm  = mcmnonew%6;
-      Int_t divmcm  = mcmnonew/6;
-      
-      row = 16*remmcm + irownew;
-      col =  4*divmcm + icolnew;
-    }
-  else if (ddlno == 5)
+  // This converts row col to hardware channel number
+  // This is the final version of mapping supplied by Mriganka
+
+    UInt_t iCh[16][4];
+
+    static const UInt_t kChDdl01[16][4] = { {6, 4, 5, 7},
+                                          {10, 2, 1, 9},
+                                          {12, 0, 3, 11},
+                                          {14, 8, 13, 15},
+                                          {16, 18, 23, 17},
+                                          {20, 28, 31, 19},
+                                          {22, 30, 29, 21},
+                                          {24, 26, 27, 25},
+                                          {38, 36, 37, 39},
+                                          {42, 34, 33, 41},
+                                          {44, 32, 35, 43},
+                                          {46, 40, 45, 47},
+                                          {48, 50, 55, 49},
+                                          {52, 60, 63, 51},
+                                          {54, 62, 61, 53},
+                                          {56, 58, 59, 57} };
+
+
+    static const UInt_t kChDdl23[16][4] = { {57, 59, 58, 56},
+                                           {53, 61, 62, 54},
+                                           {51, 63, 60, 52},
+                                           {49, 55, 50, 48},
+                                           {47, 45, 40, 46},
+                                           {43, 35, 32, 44},
+                                           {41, 33, 34, 42},
+                                           {39, 37, 36, 38},
+                                           {25, 27, 26, 24},
+                                           {21, 29, 30, 22},
+                                           {19, 31, 28, 20},
+                                           {17, 23, 18, 16},
+                                           {15, 13, 8, 14},
+                                           {11, 3, 0, 12},
+                                           {9, 1, 2, 10},
+                                           {7, 5, 4, 6} };
+    
+    
+    static const UInt_t kChDdl41[16][4] = { {56, 58, 59, 57},
+                                          {54, 62, 61, 53},
+                                          {52, 60, 63, 51},
+                                          {48, 50, 55, 49},
+                                          {46, 40, 45, 47},
+                                          {44, 32, 35, 43},
+                                          {42, 34, 33, 41},
+                                          {38, 36, 37, 39},
+                                          {24, 26, 27, 25},
+                                          {22, 30, 29, 21},
+                                          {20, 28, 31, 19},
+                                          {16, 18, 23, 17},
+                                          {14, 8, 13, 15},
+                                          {12, 0, 3, 11},
+                                          {10, 2, 1, 9},
+                                          {6, 4, 5, 7} };
+
+
+    static const UInt_t kChDdl42[16][4] = { {7, 5, 4, 6},
+                                           {9, 1, 2, 10},
+                                           {11, 3, 0, 12},
+                                           {15, 13, 8, 14},
+                                           {17, 23, 18, 16},
+                                           {19, 31, 28, 20},
+                                           {21, 29, 30, 22},
+                                           {25, 27, 26, 24},
+                                           {39, 37, 36, 38},
+                                           {41, 33, 34, 42},
+                                           {43, 35, 32, 44},
+                                           {47, 45, 40, 46},
+                                           {49, 55, 50, 48},
+                                           {51, 63, 60, 52},
+                                           {53, 61, 62, 54},
+                                           {57, 59, 58, 56} };
+
+
+    static const UInt_t kChDdl51[16][4] = { {7, 5, 4, 6},
+                                           {9, 1, 2, 10},
+                                           {11, 3, 0, 12},
+                                           {15, 13, 8, 14},
+                                           {17, 23, 18, 16},
+                                           {19, 31, 28, 20},
+                                           {21, 29, 30, 22},
+                                           {25, 27, 26, 24},
+                                           {39, 37, 36, 38},
+                                           {41, 33, 34, 42},
+                                           {43, 35, 32, 44},
+                                           {47, 45, 40, 46},
+                                           {49, 55, 50, 48},
+                                           {51, 63, 60, 52},
+                                           {53, 61, 62, 54},
+                                           {57, 59, 58, 56} };
+    
+
+
+    static const UInt_t kChDdl52[16][4] = { {56, 58, 59, 57},
+                                           {54, 62, 61, 53},
+                                           {52, 60, 63, 51},
+                                           {48, 50, 55, 49},
+                                           {46, 40, 45, 47},
+                                           {44, 32, 35, 43},
+                                           {42, 34, 33, 41},
+                                           {38, 36, 37, 39},
+                                           {24, 26, 27, 25},
+                                           {22, 30, 29, 21},
+                                           {20, 28, 31, 19},
+                                           {16, 18, 23, 17},
+                                           {14, 8, 13, 15},
+                                           {12, 0, 3, 11},
+                                           {10, 2, 1, 9},
+                                           {6, 4, 5, 7} };
+    
+    
+    for (Int_t i = 0; i < 16; i++)
+      {
+       for (Int_t j = 0; j < 4; j++)
+         {
+           
+           if(imodule < 6)                    iCh[i][j] = kChDdl01[i][j];
+           if(imodule >= 6 && imodule <= 11)  iCh[i][j] = kChDdl01[i][j];
+           if(imodule >= 12 && imodule <= 17) iCh[i][j] = kChDdl23[i][j];
+           if(imodule >= 18 && imodule <= 23) iCh[i][j] = kChDdl23[i][j];
+
+           if(imodule >= 24 && imodule <= 29) iCh[i][j] = kChDdl41[i][j];
+           if(imodule >= 42 && imodule <= 47) iCh[i][j] = kChDdl42[i][j];
+           if(imodule >= 36 && imodule <= 41) iCh[i][j] = kChDdl51[i][j];
+           if(imodule >= 30 && imodule <= 35) iCh[i][j] = kChDdl52[i][j];
+           
+         }
+      }
+
+
+  Int_t irownew = row%16;
+  Int_t icolnew = col%4;
+  
+  chno  = iCh[irownew][icolnew];
+  
+  
+  for (Int_t ibus = beginPatchBus; ibus <= endPatchBus; ibus++)
     {
-      um  = mcmno/144;
-      Int_t mcmnonew = mcmno - 72*um;
-      Int_t rowcol  = ch[chno];
-      Int_t irownew = rowcol/4;
-      Int_t icolnew = rowcol%4;
-      
-      Int_t remmcm  = mcmnonew%3;
-      Int_t divmcm  = mcmnonew/3;
-      
-      row = 16*remmcm + irownew;
-      col =  4*divmcm + icolnew;
-    }
+      Int_t srow = startRowBus[ibus];
+      Int_t erow = endRowBus[ibus];
+      Int_t scol = startColBus[ibus];
+      Int_t ecol = endColBus[ibus];
+      Int_t tmcm = mcmperBus[ibus];
 
+      if ((row >= srow && row <= erow) && (col >= scol && col <= ecol))
+       {
+         busno = ibus;
+         
+         // Find out the MCM Number
+         //
+
+         if (imodule < 6)                  mcmno = (col-scol)/4 + 1;
+         if (imodule >= 6 && imodule < 12) mcmno = (col-scol)/4 + 1;
+
+         if (imodule >= 12 && imodule < 18)
+           {
+             icolnew = (col - scol)/4;
+             mcmno = tmcm - icolnew;
+           }         
+         if (imodule >= 18 && imodule < 24)
+           {
+             icolnew = (col - scol)/4;
+             mcmno = tmcm - icolnew;
+           }         
+
+         // DDL = 4
+         if (imodule >= 24 && imodule < 30)
+           {
+
+             //if (tmcm == 24)
+             Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus];
+             if(rowdiff > 16)
+               {
+                 Int_t midrow = srow + 16;
+                 if(row >= srow && row < midrow)
+                   {
+                     mcmno = 12 + (col-scol)/4 + 1;
+                   }
+                 else if(row >= midrow && row <= erow)
+                   
+                   {
+                     mcmno = (col-scol)/4 + 1;
+                   }
+               }
+             else if (rowdiff < 16)
+               {
+                 mcmno = (col-scol)/4 + 1;
+               }
+           
+           }         
+         if (imodule >= 42 && imodule < 48)
+           {
+             Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus];
+             if(rowdiff > 16)
+               {
+                 Int_t midrow = srow + 16;
+                 if (row >= midrow && row <= erow)
+                   {
+                     mcmno = 12 + (ecol -col)/4 + 1;
+                   }
+                 else if (row >= srow && row < midrow)
+                   {
+                     mcmno = (ecol - col)/4 + 1;
+                   }
+               }
+             else if (rowdiff < 16)
+               {
+                 mcmno = (ecol - col)/4 + 1;
+               }
+           }         
+
+         // DDL = 5
+         if (imodule >= 30 && imodule < 36)
+           {
+             // CPV plane,  SU Mod = 1, 2 : ddl = 5
+             
+             //if(tmcm == 24)
+             Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus];
+             if(rowdiff > 16)
+               {
+                 Int_t midrow = srow + 16;
+                 if(row >= srow && row < midrow)
+                   {
+                     mcmno = 12 + (col-scol)/4 + 1;
+                   }
+                 else if(row >= midrow && row <= erow)
+                   {
+                     mcmno = (col-scol)/4 + 1;
+                   }
+               }
+             else if(rowdiff < 16)
+               {
+                 mcmno = (col-scol)/4 + 1;
+               }
+             
+           }
+         if (imodule >= 36 && imodule < 42)
+           {
+             Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus];
+             if(rowdiff > 16)
+               {
+                 Int_t midrow = srow + 16;
+                 if (row >= midrow && row <= erow)
+                   {
+                     mcmno = 12 + (ecol - col)/4 + 1;
+                   }
+                 else if (row >= srow && row < midrow)
+                   {
+                     mcmno = (ecol - col)/4 + 1;
+                   }
+               }
+             else if (rowdiff < 16)
+               {
+                 mcmno = (ecol - col)/4 + 1;
+               }
+           }
+
+       }
+    }
+} 
 
-}
 //____________________________________________________________________________
 
-void AliPMDDDLRawData::PackWord(UInt_t startbit, UInt_t stopbit, 
-                               UInt_t dataword, UInt_t &packedword)
+UInt_t AliPMDDDLRawData::ComputeParity(UInt_t baseword)
 {
-  UInt_t bitLength  = stopbit - startbit + 1;
-  UInt_t bitContent = (UInt_t) (TMath::Power(2,bitLength) - 1);
-  if(bitContent < dataword)
-    {
-      cout << " *** ERROR *** bitContent is less than the dataword" << endl;
-      return;
-    }
-  UInt_t packedBits = 0;
-  if (packedword != 0)
-    packedBits = (UInt_t) (TMath::Log(packedword)/TMath::Log(2));
+  // Generate the parity bit
 
-  UInt_t counter;
-  if (packedBits <= stopbit)
+  UInt_t count = 0;
+  for(Int_t j=0; j<29; j++)
     {
-      counter   = 31 - stopbit;
+      if (baseword & 0x01 ) count++;
+      baseword >>= 1;
     }
-  else
+  UInt_t parity = count%2;
+  return parity;
+}
+//____________________________________________________________________________
+void AliPMDDDLRawData::DdlMapping(Int_t iddl, Int_t imodule,
+                                 Int_t &beginPatchBus, Int_t &endPatchBus,
+                                 Int_t patchBusNo[], Int_t mcmperBus[],
+                                 Int_t startRowBus[], Int_t endRowBus[],
+                                 Int_t startColBus[], Int_t endColBus[])
+{
+  // DDL Mapping fetching from mapping database
+
+  beginPatchBus = fMapData->GetBeginPatchBus(iddl,imodule);
+  endPatchBus   = fMapData->GetEndPatchBus(iddl,imodule);
+
+  for(Int_t ibus = beginPatchBus; ibus < endPatchBus+1; ibus++)
     {
-      counter   = 31 - packedBits;
+      patchBusNo[ibus]   = ibus;
+      mcmperBus[ibus]    = fMapData->GetMcmperBus(iddl,ibus);
+      startRowBus[ibus]  = fMapData->GetStartRowBus(iddl,ibus);
+      startColBus[ibus]  = fMapData->GetStartColBus(iddl,ibus);
+      endRowBus[ibus]    = fMapData->GetEndRowBus(iddl,ibus);
+      endColBus[ibus]    = fMapData->GetEndColBus(iddl,ibus);
     }
-  UInt_t dummyword = 0xFFFFFFFF;
-  dummyword >>= counter;
-  UInt_t lword = dataword << startbit;
-  UInt_t nword = lword | packedword;
-  packedword = dummyword & nword;
 
+}
+//____________________________________________________________________________
 
+AliPMDddlinfoData* AliPMDDDLRawData::GetDdlinfoData() const
+{
+  // The run number will be centralized in AliCDBManager,
+  // you don't need to set it here!
+  AliCDBEntry  *entry = AliCDBManager::Instance()->Get("PMD/Calib/Ddlinfo");
+  
+  if(!entry) AliFatal("ddlinfo object retrieval failed!");
+  
+  AliPMDddlinfoData *ddlinfo = 0;
+  if (entry) ddlinfo = (AliPMDddlinfoData*) entry->GetObject();
+  
+  if (!ddlinfo)  AliFatal("No ddl info data from  database !");
+  
+  return ddlinfo;
 }
 //____________________________________________________________________________
-void AliPMDDDLRawData::UnpackWord(UInt_t startbit, UInt_t stopbit, 
-                               UInt_t &dataword, UInt_t packedword)
+AliPMDMappingData* AliPMDDDLRawData::GetMappingData() const
 {
-  UInt_t bitLength  = stopbit - startbit + 1;
-  UInt_t bitContent = (UInt_t) (TMath::Power(2,bitLength) - 1);
-  bitContent <<= startbit;
-  dataword = packedword & bitContent;
-  dataword >>= startbit;
-
+  // The run number will be centralized in AliCDBManager,
+  // you don't need to set it here!
+  AliCDBEntry  *entry = AliCDBManager::Instance()->Get("PMD/Calib/Mapping");
+  
+  if(!entry) AliFatal("Mapping object retrieval failed!");
+  
+  AliPMDMappingData *mapda = 0;
+  if (entry) mapda = (AliPMDMappingData*) entry->GetObject();
+  
+  if (!mapda)  AliFatal("No mapping data from  database !");
+  
+  return mapda;
 }
 //____________________________________________________________________________
-