defects fixed
authorbasanta <basanta@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 8 Oct 2010 16:23:17 +0000 (16:23 +0000)
committerbasanta <basanta@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 8 Oct 2010 16:23:17 +0000 (16:23 +0000)
PMD/AliPMDDDLRawData.cxx
PMD/AliPMDDDLRawData.h
PMD/AliPMDMappingData.cxx
PMD/AliPMDMappingData.h
PMD/CreateMappingCDB.C

index 6c0a1b9..53a6732 100644 (file)
 #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
@@ -45,6 +52,8 @@ AliPMDDDLRawData::AliPMDDDLRawData():
 //____________________________________________________________________________
 AliPMDDDLRawData::AliPMDDDLRawData(const AliPMDDDLRawData& ddlraw):
   TObject(ddlraw),
+  fDdlinfo(ddlraw.fDdlinfo),
+  fMapData(ddlraw.fMapData),
   fDigits(ddlraw.fDigits)
 {
   //Copy Constructor 
@@ -55,7 +64,9 @@ AliPMDDDLRawData & AliPMDDDLRawData::operator=(const AliPMDDDLRawData& ddlraw)
   //Assignment operator 
   if(this != &ddlraw)
     {
-      fDigits = ddlraw.fDigits;
+      fDdlinfo = ddlraw.fDdlinfo;
+      fMapData = ddlraw.fMapData;
+      fDigits  = ddlraw.fDigits;
     }
   return *this;
 }
@@ -87,7 +98,6 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD)
   AliDebug(1,Form("Number of modules inside treeD = %d",nmodules));
 
   const Int_t kDDL          = AliDAQ::NumberOfDdls("PMD");
-  Int_t modulePerDDL        = 0;
 
 
   AliRawDataHeaderSim header;
@@ -103,36 +113,22 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD)
 
   Char_t filename[80];
 
-  // open the ddl file info to know the modules per DDL
-  TString ddlinfofileName(gSystem->Getenv("ALICE_ROOT"));
-  ddlinfofileName += "/PMD/PMD_ddl_info.dat";
-
-  ifstream infileddl;
-  infileddl.open(ddlinfofileName.Data(), ios::in); // ascii file
-  if(!infileddl) AliError("Could not read the ddl info file");
-
-
-
+  Int_t modulePerDDL        = 0;
   Int_t mmodule = 0;
   Int_t ddlno;
-  Int_t modno;
-  Int_t *modulenoddl = 0x0;
+  Int_t modulenoddl[12] = {-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1};
 
   for(Int_t iddl = 0; iddl < kDDL; iddl++)
     {
-      if (infileddl.eof()) break;
-      infileddl >> ddlno >> modulePerDDL;
-
+      ddlno = iddl;
+      modulePerDDL = fDdlinfo->GetNoOfModulePerDdl(iddl);
       if (modulePerDDL == 0) continue;
-
-      modulenoddl = new Int_t [modulePerDDL];
-      for (Int_t im = 0; im < modulePerDDL; im++)
+      for (Int_t im = 0; im < 12; im++)
        {
-         infileddl >> modno;
-         modulenoddl[im] = modno;
+         modulenoddl[im] = fDdlinfo->GetModulesPerDdl(iddl,im);;
        }
 
-      strcpy(filename,AliDAQ::DdlFileName("PMD",iddl));
+      strncpy(filename,AliDAQ::DdlFileName("PMD",iddl),80);
       
       outfile = new AliFstream(filename);
       
@@ -149,16 +145,15 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD)
            }
        }
 
-      for(Int_t ium = 0; ium < modulePerDDL; ium++)
+      for(Int_t ium = 0; ium < 12; ium++)
        {
-         //if (iddl == 4 && ium == 6) mmodule = 42;
-
          // Extract energy deposition per cell and pack it
          // in a 32-bit word and returns all the total words
          // per one unit-module
          
          mmodule = modulenoddl[ium];
-         GetUMDigitsData(treeD, mmodule, iddl, modulePerDDL, contentsBus, busPatch);
+         if(mmodule == -1) continue;
+         GetUMDigitsData(treeD, mmodule, iddl, contentsBus, busPatch);
        }
 
       Int_t ij = 0;
@@ -366,13 +361,11 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD)
       delete outfile;
     } // DDL Loop over
 
-  delete [] modulenoddl;
-  infileddl.close();
 }
 //____________________________________________________________________________
+
 void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule,
-                                      Int_t ddlno, Int_t modulePerDDL,
-                                      Int_t *contentsBus,
+                                      Int_t ddlno,  Int_t *contentsBus,
                                       UInt_t busPatch[][1536])
 {
   // Retrieves digits data UnitModule by UnitModule
@@ -385,9 +378,9 @@ void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule,
   Int_t  det = 0, smn = 0, irow = 0, icol = 0;
   Int_t  parity = 0;
   
-  Int_t totPatchBus = 0, bPatchBus = 0, ePatchBus = 0;
-  Int_t ibus = 0, totmcm = 0, rows = 0, cols = 0, rowe = 0, cole = 0;
-  Int_t moduleno = 0;
+  //  Int_t totPatchBus = 0, bPatchBus = 0, ePatchBus = 0;
+  //  Int_t ibus = 0, totmcm = 0, rows = 0, cols = 0, rowe = 0, cole = 0;
+  //  Int_t moduleno = 0;
   Int_t busno = 0;
   Int_t patchBusNo[kMaxBus], mcmperBus[kMaxBus];
   Int_t startRowBus[kMaxBus], startColBus[kMaxBus];
@@ -405,104 +398,34 @@ void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule,
       endColBus[i]   = -1;
     }
 
+  // Fetch the DDL mapping info from the mapping database
 
-  TString fileName(gSystem->Getenv("ALICE_ROOT"));
-
-  if(ddlno == 0)
-    {
-      fileName += "/PMD/PMD_Mapping_ddl0.dat";
-    }
-  else if(ddlno == 1)
-    {
-      fileName += "/PMD/PMD_Mapping_ddl1.dat";
-    }
-  else if(ddlno == 2)
-    {
-      fileName += "/PMD/PMD_Mapping_ddl2.dat";
-    }
-  else if(ddlno == 3)
-    {
-      fileName += "/PMD/PMD_Mapping_ddl3.dat";
-    }
-  else if(ddlno == 4)
-    {
-      fileName += "/PMD/PMD_Mapping_ddl4.dat";
-    }
-  else if(ddlno == 5)
-    {
-      fileName += "/PMD/PMD_Mapping_ddl5.dat";
-    }
-
-  ifstream infile;
-  infile.open(fileName.Data(), ios::in); // ascii file
-  if(!infile)
-    AliError(Form("Could not read the mapping file for DDL No = %d",ddlno));
-
-  for (Int_t im = 0; im < modulePerDDL; im++)
-    {
-      infile >> moduleno;
-      infile >> totPatchBus >> bPatchBus >> ePatchBus;
-      
-      if (totPatchBus == 0) continue;    // BKN
+  DdlMapping(ddlno, imodule, beginPatchBus, endPatchBus,
+            patchBusNo, mcmperBus, startRowBus, endRowBus,
+            startColBus, endColBus);
 
-      if (moduleno == imodule)
-       {
-         beginPatchBus = bPatchBus;
-         endPatchBus   = ePatchBus;
-       }
-      
-      for(Int_t i=0; i<totPatchBus; i++)
-       {
-         infile >> ibus >> totmcm >> rows >> rowe >> cols >> cole;
-
-         if (moduleno == imodule)
-           {
-             patchBusNo[ibus]   = ibus;
-             mcmperBus[ibus]    = totmcm;
-             startRowBus[ibus]  = rows;
-             startColBus[ibus]  = cols;
-             endRowBus[ibus]    = rowe;
-             endColBus[ibus]    = cole;
-           }
-       }
-
-    }
-
-  infile.close();
-
-  // Read if some chains are off
-  TString rchainName(gSystem->Getenv("ALICE_ROOT"));
-  rchainName += "/PMD/PMD_removed_chains.dat";
-
-  ifstream rchainfile;
-  rchainfile.open(rchainName.Data(), ios::in); // ascii file
-  if(!rchainfile)AliError("Could not read the removed cahins file");
+  // 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];
 
-  Int_t rows1 = 0, rowe1 = 0, cols1 = 0, cole1 = 0;
-  Int_t rows2 = 0, rowe2 = 0, cols2 = 0, cole2 = 0;
-
-  for (Int_t im = 0; im < 48; im++)
+  for (Int_t idet = 0; idet < 2; idet++)
     {
-      rchainfile >> det >> smn >> rows1 >> rowe1 >> cols1 >> cole1
-                >> rows2 >> rowe2 >> cols2 >> cole2;
-      
-      srowoff1[det][smn] = rows1;
-      erowoff1[det][smn] = rowe1;
-      scoloff1[det][smn] = cols1;
-      ecoloff1[det][smn] = cole1;
-      srowoff2[det][smn] = rows2;
-      erowoff2[det][smn] = rowe2;
-      scoloff2[det][smn] = cols2;
-      ecoloff2[det][smn] = cole2;
+      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);
+       }
     }
 
-  rchainfile.close();
-
   treeD->GetEntry(imodule); 
   Int_t nentries = fDigits->GetLast();
   Int_t totword = nentries+1;
@@ -881,5 +804,60 @@ Int_t AliPMDDDLRawData::ComputeParity(UInt_t baseword)
   Int_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++)
+    {
+      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);
+    }
 
+}
+//____________________________________________________________________________
+
+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;
+}
+//____________________________________________________________________________
+AliPMDMappingData* AliPMDDDLRawData::GetMappingData() const
+{
+  // 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;
+}
 //____________________________________________________________________________
index 76f6b47..313d9cd 100644 (file)
 
 class TClonesArray;
 class TTree;
-
+class AliPMDddlinfoData;
+class AliPMDMappingData;
 class AliPMDdigit;
+class AliCDBManager;
+class AliCDBStorage;
+class AliCDBEntry;
 
 class AliPMDDDLRawData:public TObject
 {
@@ -29,8 +33,7 @@ class AliPMDDDLRawData:public TObject
 
   void WritePMDRawData(TTree *treeD);
   void GetUMDigitsData(TTree *treeD, Int_t imodule, Int_t ddlno,
-                      Int_t modulePerDDL, Int_t *contentsBus,
-                      UInt_t busPatch[][1536]);
+                      Int_t *contentsBus, UInt_t busPatch[][1536]);
   void TransformS2H(Int_t smn, Int_t &irow, Int_t &icol);
   void GetMCMCh(Int_t imodule, Int_t row, Int_t col,
                Int_t beginPatchBus, Int_t endPatchBus,
@@ -38,14 +41,27 @@ class AliPMDDDLRawData:public TObject
                Int_t *startRowBus, Int_t *startColBus,
                Int_t *endRowBus, Int_t *endColBus,
                Int_t & busno, UInt_t &mcmno, UInt_t &chno);
+  void 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[]);
+
+
+  AliPMDddlinfoData  *GetDdlinfoData() const;
+  AliPMDMappingData  *GetMappingData() const;
+
 
  protected:
 
+  AliPMDddlinfoData  *fDdlinfo;    //! ddl info data
+  AliPMDMappingData  *fMapData;    //! Mapping data
+
   Int_t ComputeParity(UInt_t baseword);
 
   TClonesArray *fDigits;    //! List of digits
 
-  ClassDef(AliPMDDDLRawData,9)    // To make RAW Data
+  ClassDef(AliPMDDDLRawData,10)    // To make RAW Data
 };
 #endif
 
index 85dcdc5..216a067 100644 (file)
@@ -48,10 +48,19 @@ AliPMDMappingData::AliPMDMappingData(const AliPMDMappingData& mapda) :
   Reset();
   for(Int_t iddl = 0; iddl < kDdl; iddl++)
     {
+      for(Int_t imod = 0; imod < 48; imod++)
+       {
+         fBeginPatchBus[iddl][imod] = mapda.GetBeginPatchBus(iddl,imod);
+         fEndPatchBus[iddl][imod]   = mapda.GetEndPatchBus(iddl,imod);
+       }
       for(Int_t ibus = 0; ibus < kBus; ibus++)
        {
-         fModuleNo[iddl][ibus]=
-           mapda.GetModuleNo(iddl,ibus);
+         fModuleNo[iddl][ibus]    = mapda.GetModuleNo(iddl,ibus);
+         fMcmperBus[iddl][ibus]   = mapda.GetMcmperBus(iddl,ibus);
+         fStartRowBus[iddl][ibus] = mapda.GetStartRowBus(iddl,ibus);
+         fEndRowBus[iddl][ibus]   = mapda.GetEndRowBus(iddl,ibus);
+         fStartColBus[iddl][ibus] = mapda.GetStartColBus(iddl,ibus);
+         fEndColBus[iddl][ibus]   = mapda.GetEndColBus(iddl,ibus);
        }
     }
 
@@ -66,10 +75,19 @@ AliPMDMappingData &AliPMDMappingData::operator =(const AliPMDMappingData& mapda)
 
   for(Int_t iddl = 0; iddl < kDdl; iddl++)
     {
+      for(Int_t imod = 0; imod < 48; imod++)
+       {
+         fBeginPatchBus[iddl][imod] = mapda.GetBeginPatchBus(iddl,imod);
+         fEndPatchBus[iddl][imod]   = mapda.GetEndPatchBus(iddl,imod);
+       }
       for(Int_t ibus = 0; ibus < kBus; ibus++)
        {
-         fModuleNo[iddl][ibus]=
-           mapda.GetModuleNo(iddl,ibus);
+         fModuleNo[iddl][ibus]    = mapda.GetModuleNo(iddl,ibus);
+         fMcmperBus[iddl][ibus]   = mapda.GetMcmperBus(iddl,ibus);
+         fStartRowBus[iddl][ibus] = mapda.GetStartRowBus(iddl,ibus);
+         fEndRowBus[iddl][ibus]   = mapda.GetEndRowBus(iddl,ibus);
+         fStartColBus[iddl][ibus] = mapda.GetStartColBus(iddl,ibus);
+         fEndColBus[iddl][ibus]   = mapda.GetEndColBus(iddl,ibus);
        }
     }
 
@@ -83,6 +101,14 @@ AliPMDMappingData::~AliPMDMappingData()
 // ----------------------------------------------------------------- //
 void AliPMDMappingData::Reset()
 {
+  for(Int_t i = 0; i < 6; i++)
+    {
+      for(Int_t j = 0; j < 48; j++)
+       {
+         fBeginPatchBus[i][j] = -1;
+         fEndPatchBus[i][j] = -1;
+       }
+    }
 
   for(Int_t iddl = 0; iddl < kDdl; iddl++)
     {
@@ -99,6 +125,16 @@ void AliPMDMappingData::Reset()
 
 }
 // ----------------------------------------------------------------- //
+Int_t AliPMDMappingData:: GetBeginPatchBus(Int_t iddl, Int_t imod) const
+{
+  return fBeginPatchBus[iddl][imod];
+}
+// ----------------------------------------------------------------- //
+Int_t AliPMDMappingData:: GetEndPatchBus(Int_t iddl, Int_t imod) const
+{
+  return fEndPatchBus[iddl][imod];
+}
+// ----------------------------------------------------------------- //
 Int_t AliPMDMappingData:: GetModuleNo(Int_t iddl, Int_t ibus) const
 {
   return fModuleNo[iddl][ibus];
@@ -129,6 +165,13 @@ Int_t AliPMDMappingData:: GetEndColBus(Int_t iddl, Int_t ibus) const
   return fEndColBus[iddl][ibus];
 }
 // ----------------------------------------------------------------- //
+void AliPMDMappingData::SetPatchBus(Int_t iddl, Int_t imod,
+                                   Int_t bpatchbus, Int_t epatchbus)
+{
+  fBeginPatchBus[iddl][imod] = bpatchbus;
+  fEndPatchBus[iddl][imod]   = epatchbus;
+}
+// ----------------------------------------------------------------- //
 void AliPMDMappingData::SetModuleNo(Int_t iddl, Int_t ibus, Int_t modno)
 {
   fModuleNo[iddl][ibus] = modno;
@@ -161,6 +204,11 @@ void AliPMDMappingData::Print(Option_t *) const
 
   for(Int_t iddl = 0; iddl < kDdl; iddl++)
     {
+      for(Int_t imod = 0; imod < 48; imod++)
+       {
+         printf("%d %d %d %d \n",iddl, imod, fBeginPatchBus[iddl][imod],
+                fEndPatchBus[iddl][imod]);
+       }
       for(Int_t ibus = 0; ibus < kBus; ibus++)
        {
          printf("%d %d %d %d %d %d %d %d\n",iddl, ibus,
index 6f4771f..ac12763 100644 (file)
@@ -18,6 +18,8 @@ class AliPMDMappingData: public TNamed
   virtual ~AliPMDMappingData();
   void  Reset();
 
+  Int_t GetBeginPatchBus(Int_t iddl, Int_t imod) const;
+  Int_t GetEndPatchBus(Int_t iddl, Int_t imod) const;
   Int_t GetModuleNo(Int_t iddl, Int_t ibus) const;
   Int_t GetMcmperBus(Int_t iddl, Int_t ibus) const;
   Int_t GetStartRowBus(Int_t iddl, Int_t ibus) const;
@@ -25,6 +27,8 @@ class AliPMDMappingData: public TNamed
   Int_t GetStartColBus(Int_t iddl, Int_t ibus) const;
   Int_t GetEndColBus(Int_t iddl, Int_t ibus) const;
 
+  
+  void  SetPatchBus(Int_t iddl, Int_t imod, Int_t bpatchbus, Int_t epatchbus);
   void  SetModuleNo(Int_t iddl, Int_t ibus, Int_t modno);
   void  SetMcmperBus(Int_t iddl, Int_t ibus, Int_t totmcm);
   void  SetRowBus(Int_t iddl, Int_t ibus, Int_t rows, Int_t rowe);
@@ -39,6 +43,8 @@ class AliPMDMappingData: public TNamed
        kBus = 51    // Modules of patch bus
       };
 
+  Int_t fBeginPatchBus[6][48];
+  Int_t fEndPatchBus[6][48];
   Int_t fModuleNo[kDdl][kBus];
   Int_t fMcmperBus[kDdl][kBus];
   Int_t fStartRowBus[kDdl][kBus];
@@ -47,6 +53,6 @@ class AliPMDMappingData: public TNamed
   Int_t fEndColBus[kDdl][kBus];
 
 
-  ClassDef(AliPMDMappingData,1) // calibration class for gainfactors
+  ClassDef(AliPMDMappingData,2) // calibration class for gainfactors
 };
 #endif
index 6a4c4f7..f276e44 100644 (file)
@@ -32,7 +32,7 @@ void CreateMappingCDB()
        {
          infile >> moduleno;
          infile >> totPatchBus >> bPatchBus >> ePatchBus;
-         
+         mapda->SetPatchBus(iddl,moduleno,bPatchBus,ePatchBus);
          if (totPatchBus == 0) continue;
          
          for(Int_t i=0; i<totPatchBus; i++)