]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TOF/AliTOFDigitMap.cxx
bugfix: corrected calculation of slice and partition from track point Id
[u/mrichter/AliRoot.git] / TOF / AliTOFDigitMap.cxx
index 56c47786e9c11b9d9ec4104f132967878c42aef8..da9be4df8884db24481041af63fa9d88b06a4118 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.12  2007/02/20 15:57:00  decaro
+Raw data update: to read the TOF raw data defined in UNPACKED mode
+
 */
 
 ////////////////////////////////////////////////////////////////////////
@@ -46,48 +49,54 @@ $Log$
 ClassImp(AliTOFDigitMap)
 
 AliTOFDigitMap::AliTOFDigitMap():
-  fNSector(-1),
-  fNplate(-1),
-  fNstrip(-1),
-  fNpx(-1),
-  fNpz(-1),
+  fNSector(AliTOFGeometry::NSectors()),
+  fNplate(AliTOFGeometry::NPlates()),
+  fNstrip(AliTOFGeometry::NStripC()),
+  fNpx(AliTOFGeometry::NpadX()),
+  fNpz(AliTOFGeometry::NpadZ()),
   fMaxIndex(-1),
-  fDigitMap(0x0),
-  fTOFGeometry(new AliTOFGeometry())
+  fDigitMap(0x0)
 {
 //
 // Default ctor
 //
 
-  fNSector = AliTOFGeometry::NSectors();
-  fNplate = AliTOFGeometry::NPlates();
-  fNstrip = fTOFGeometry->NStripC();//fTOFGeometry->NMaxNstrip();
-  fNpx  = AliTOFGeometry::NpadX();
-  fNpz  = AliTOFGeometry::NpadZ();
   fMaxIndex=fNSector*fNplate*fNstrip*fNpx*fNpz;
-
   fDigitMap = new Int_t*[fMaxIndex];
+
   for (Int_t i=0; i<fMaxIndex; i++) fDigitMap[i] = new Int_t[kMaxDigitsPerPad];
   Clear();
 }
 
 ////////////////////////////////////////////////////////////////////////
-AliTOFDigitMap::AliTOFDigitMap(const AliTOFDigitMap & /*digitMap*/):
-  TObject(),
-  fNSector(-1),
-  fNplate(-1),
-  fNstrip(-1),
-  fNpx(-1),
-  fNpz(-1),
-  fMaxIndex(-1),
-  fDigitMap(0x0),
-  fTOFGeometry(new AliTOFGeometry())
+AliTOFDigitMap::AliTOFDigitMap(const AliTOFDigitMap & digitMap):
+  TObject(digitMap),
+  fNSector(digitMap.fNSector),
+  fNplate(digitMap.fNplate),
+  fNstrip(digitMap.fNstrip),
+  fNpx(digitMap.fNpx),
+  fNpz(digitMap.fNpz),
+  fMaxIndex(digitMap.fMaxIndex),
+  fDigitMap(0x0)
 {
 //
-// Dummy copy constructor
+// dummy copy constructor
 //
-  ;
+  
+  fMaxIndex=fNSector*fNplate*fNstrip*fNpx*fNpz;
+  fDigitMap = new Int_t*[fMaxIndex];
 
+  for (Int_t i=0; i<fMaxIndex; i++) fDigitMap[i] = new Int_t[kMaxDigitsPerPad];
+}
+
+////////////////////////////////////////////////////////////////////////
+AliTOFDigitMap &
+AliTOFDigitMap::operator=(const AliTOFDigitMap & /*digitMap*/)
+{
+//
+// dummy copy const
+//
+    return *this;
 }
 
  
@@ -99,9 +108,9 @@ AliTOFDigitMap::~AliTOFDigitMap()
 //
   if (fDigitMap) {
     for (Int_t i=0; i<fMaxIndex; i++)  delete[] fDigitMap[i];
+    delete [] fDigitMap;
   }
 
-  fTOFGeometry = 0;
 
 }
 
@@ -121,7 +130,7 @@ void AliTOFDigitMap::Clear(const Option_t*)
 }
 
 ////////////////////////////////////////////////////////////////////////
-Int_t AliTOFDigitMap::CheckedIndex(Int_t *vol) const
+Int_t AliTOFDigitMap::CheckedIndex(Int_t * const vol) const
 {
   //
   // Return checked index for vol
@@ -150,6 +159,13 @@ void AliTOFDigitMap::AddDigit(Int_t *vol, Int_t idigit)
   //
   // 0 means empty pad, we need to shift indeces by 1
 
+  if (fDigitMap[CheckedIndex(vol)][kMaxDigitsPerPad-1]!=0) {
+    AliDebug(1,Form("In the volume (Se%i, Pl%i, St%i, PadR%i, Pad%i) there is not more possibility to add other digits.", vol[0], vol[1], vol[2], vol[4], vol[3]));
+    AliDebug(1,Form("Then, the digit number %i will be not inserted in the digit map, i.e. it will be lost.", idigit));
+    AliDebug(1,Form("Please, check the possibility to increase the digit map size (succently set to %i)", kMaxDigitsPerPad));
+    return;
+  }
+
   for (Int_t slot=0; slot<kMaxDigitsPerPad; slot++) {
 
     if (fDigitMap[CheckedIndex(vol)][slot]==0) {
@@ -159,6 +175,7 @@ void AliTOFDigitMap::AddDigit(Int_t *vol, Int_t idigit)
     //else continue;
 
   }
+
 }
 
 ////////////////////////////////////////////////////////////////////////
@@ -192,9 +209,10 @@ Int_t AliTOFDigitMap::GetDigitIndex(Int_t *vol, Int_t label) const
     return -1;
   }
 
-  if (CheckedIndex(vol)==-1) return -1;
+  Int_t ci = CheckedIndex(vol);
+  if (ci==-1) return -1;
   
-  Int_t dummy = fDigitMap[CheckedIndex(vol)][label];
+  Int_t dummy = fDigitMap[ci][label];
   
   if (dummy>0) return dummy-1;
   else return -1;
@@ -218,37 +236,159 @@ FlagType AliTOFDigitMap::TestDigit(Int_t *vol) const
 }
 
 ////////////////////////////////////////////////////////////////////////
-AliTOFDigitMap & AliTOFDigitMap::operator = (const AliTOFDigitMap & /*rhs*/) 
+Int_t AliTOFDigitMap::GetFilledCellNumber() const
 {
-// Dummy assignment operator
-    return *this;
+  //
+  // Returns the number of filled cells of the TOF digit map
+  //
+
+  Int_t counter = 0;
+
+  for (Int_t index = 0; index < fMaxIndex; ++index)
+  {
+    for (Int_t label = 0; label < kMaxDigitsPerPad; ++label)
+    {
+      if (fDigitMap[index][label] > 0)
+      {
+       ++counter;
+       break;
+      }
+    }
+  }
+
+  return counter;
 }
+
 ////////////////////////////////////////////////////////////////////////
-Int_t AliTOFDigitMap::GetFilledCellNumber() const
+Bool_t AliTOFDigitMap::StripDigitCheck(Int_t iSector, Int_t iPlate, Int_t iStrip) const
 {
   //
-  // Returns the number of filled cells of the TOF digit map
+  // Returns:
+  //           kFALSE if the strip doesn't contain digits
+  //           kTRUE  if the strip contains at least one digit
+  //
+
+  Int_t volume[5] = {iSector, iPlate, iStrip, -1, -1};
+  Bool_t counter = kFALSE;
+
+  for (Int_t iPadX=0; iPadX<fNpx; iPadX++)
+    for (Int_t iPadZ=0; iPadZ<fNpz; iPadZ++)
+      {
+       volume[3] = iPadX;
+       volume[4] = iPadZ;
+       for (Int_t label=0; label<kMaxDigitsPerPad; label++) {
+         if (GetDigitIndex(volume, label)>=0) {
+           counter = kTRUE;
+           break;
+         }
+       }
+      }
+
+  return counter;
+
+}
+
+////////////////////////////////////////////////////////////////////////
+Int_t AliTOFDigitMap::DigitInStrip(Int_t iSector, Int_t iPlate, Int_t iStrip) const
+{
+  //
+  // Returns number of digits in the strip iStrip,
+  //         in the plate iPlate of the sector iSector
+  //
+
+  Int_t volume[5] = {iSector, iPlate, iStrip, -1, -1};
+  Int_t counter = 0;
+
+  for (Int_t iPadX=0; iPadX<fNpx; iPadX++)
+    for (Int_t iPadZ=0; iPadZ<fNpz; iPadZ++)
+      for (Int_t label=0; label<kMaxDigitsPerPad; label++) {
+       volume[3] = iPadX;
+       volume[4] = iPadZ;
+       if (GetDigitIndex(volume, label)>=0)
+         counter++;
+      }
+
+  return counter;
+
+}
+
+////////////////////////////////////////////////////////////////////////
+Int_t AliTOFDigitMap::FilledCellsInStrip(Int_t iSector, Int_t iPlate, Int_t iStrip) const
+{
+  //
+  // Returns number of digits in the strip iStrip,
+  //         in the plate iPlate of the sector iSector
+  //
+
+  Int_t volume[5] = {iSector, iPlate, iStrip, -1, -1};
+  Int_t counter = 0;
+
+  for (Int_t iPadX=0; iPadX<fNpx; iPadX++)
+    for (Int_t iPadZ=0; iPadZ<fNpz; iPadZ++) {
+      volume[3] = iPadX;
+      volume[4] = iPadZ;
+      if (GetDigitIndex(volume, 0)>=0)
+       counter++;
+    }
+
+  return counter;
+
+}
+
+////////////////////////////////////////////////////////////////////////
+void AliTOFDigitMap::ResetDigitNumber(Int_t *vol, Int_t dig)
+{
+  //
+  // Reset digit into pad vol
+  //
+
+  for (Int_t slot=0; slot<kMaxDigitsPerPad; slot++) {
+    if (fDigitMap[CheckedIndex(vol)][slot]-1==dig) {
+      fDigitMap[CheckedIndex(vol)][slot] = 0;
+    }
+  }
+
+}
+
+////////////////////////////////////////////////////////////////////////
+void AliTOFDigitMap::ResetDigit(Int_t *vol, Int_t dig)
+{
+  //
+  // Reset digit into pad vol
+  //
+  // 0 means empty pad, we need to shift indeces by 1
+
+  fDigitMap[CheckedIndex(vol)][dig] = 0;
+
+}
+
+////////////////////////////////////////////////////////////////////////
+void AliTOFDigitMap::ResetDigit(Int_t *vol)
+{
+  //
+  // Reset digit into pad vol
+  //
+  // 0 means empty pad, we need to shift indices by 1
+
+  for (Int_t slot=0; slot<kMaxDigitsPerPad; slot++)
+    fDigitMap[CheckedIndex(vol)][slot] = 0;
+
+}
+
+////////////////////////////////////////////////////////////////////////
+Int_t AliTOFDigitMap::GetNumberOfDigits(Int_t *vol)
+{
+  //
+  // Returns the number of digit
+  //   into pad volume vol
+  //
+  // 0 means empty pad
   //
 
-  Int_t volume[5] = {-1, -1, -1, -1, -1};
   Int_t counter = 0;
 
-  for (Int_t iSector=0; iSector<fNSector; iSector++)
-    for (Int_t iPlate=0; iPlate<fNplate; iPlate++)
-      for (Int_t iStrip=0; iStrip<fNstrip; iStrip++)
-       for (Int_t iPadX=0; iPadX<fNpx; iPadX++)
-         for (Int_t iPadZ=0; iPadZ<fNpz; iPadZ++)
-           {
-
-             volume[0] = iSector;
-             volume[1] = iPlate;
-             volume[2] = iStrip;
-             volume[3] = iPadX;
-             volume[4] = iPadZ;
-
-             //if (CheckedIndex(volume)!=-1) counter++;
-             if (GetDigitIndex(volume, 0)>0) counter++;
-           }
+  for (Int_t slot=0; slot<kMaxDigitsPerPad; slot++)
+    if (GetDigitIndex(vol, slot)>=0) counter++;
 
   return counter;