Add functions for specific pad selections

* edge pads
* stack boundary rows
* cross region

Author: wiechula

Detectors/TPC/base/include/TPCBase/Mapper.h

@@ -396,6 +396,14 @@ class Mapper
 
   bool isOutOfSector(GlobalPosition3D posEle, const Sector& sector, const float margin = 0.f) const;
 
+  static bool isEdgePad(int rowInSector, int padInRow);
+  static bool isFirstOrLastRowInStack(int rowInSector);
+  static bool isBelowSpacerCross(int rowInSector, int padInRow);
+  static bool isHighCouplingPad(int rowInSector, int padInRow)
+  {
+    return isEdgePad(rowInSector, padInRow) || isFirstOrLastRowInStack(rowInSector) || isBelowSpacerCross(rowInSector, padInRow);
+  }
+
   static constexpr unsigned short getNumberOfIROCs() { return 36; }
   static constexpr unsigned short getNumberOfOROCs() { return 36; }
   static constexpr unsigned short getPadsInIROC() { return mPadsInIROC; }
@@ -523,6 +531,7 @@ class Mapper
   static constexpr unsigned int GLOBALPADOFFSET[NREGIONS]{0, 1200, 2400, 3840, 5280, 6720, 8160, 9760, 11360, 12960};                                                 ///< offset of number of pads for region
   static constexpr unsigned int ROWSPERREGION[NREGIONS]{17, 15, 16, 15, 18, 16, 16, 14, 13, 12};                                                                      ///< number of pad rows for region
   static constexpr unsigned int ROWOFFSET[NREGIONS]{0, 17, 32, 48, 63, 81, 97, 113, 127, 140};                                                                        ///< offset to calculate local row from global row
+  static constexpr unsigned int ROWOFFSETSTACK[4]{0, 63, 97, 127};                                                                                                    ///< offset to calculate local row from global row
   static constexpr float REGIONAREA[NREGIONS]{374.4f, 378.f, 453.6f, 470.88f, 864.f, 864.f, 1167.36f, 1128.96f, 1449.6f, 1456.8f};                                    ///< volume of each region in cm^2
   static constexpr float INVPADAREA[NREGIONS]{1 / 0.312f, 1 / 0.315f, 1 / 0.315f, 1 / 0.327f, 1 / 0.6f, 1 / 0.6f, 1 / 0.7296f, 1 / 0.7056f, 1 / 0.906f, 1 / 0.9105f}; ///< inverse size of the pad area padwidth*padLength
   static constexpr unsigned REGION[PADROWS] = {
@@ -542,7 +551,7 @@ class Mapper
     {0, 0, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4},             // region 7
     {0, 0, 1, 1, 2, 2, 3, 3, 3, 4, 4, 5, 5},                // region 8
     {0, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5, 5}                    // region 9
-  };                                                        ///< additional pads per row compared to first row
+  }; ///< additional pads per row compared to first row
   const inline static std::vector<unsigned int> OFFSETCRULOCAL[NREGIONS]{
     {0, 66, 132, 198, 266, 334, 402, 472, 542, 612, 684, 756, 828, 902, 976, 1050, 1124},         // region 0
     {0, 76, 152, 228, 306, 384, 462, 542, 622, 702, 784, 866, 948, 1032, 1116},                   // region 1
@@ -554,7 +563,7 @@ class Mapper
     {0, 110, 220, 332, 444, 556, 670, 784, 898, 1014, 1130, 1246, 1364, 1482},                    // region 7
     {0, 118, 236, 356, 476, 598, 720, 844, 968, 1092, 1218, 1344, 1472},                          // region 8
     {0, 128, 258, 388, 520, 652, 784, 918, 1052, 1188, 1324, 1462}                                // region 9
-  };                                                                                              ///< row offset in cru for given local pad row
+  }; ///< row offset in cru for given local pad row
   const inline static std::vector<unsigned int> PADSPERROW[NREGIONS]{
     {66, 66, 66, 68, 68, 68, 70, 70, 70, 72, 72, 72, 74, 74, 74, 74, 76},      // region 0
     {76, 76, 76, 78, 78, 78, 80, 80, 80, 82, 82, 82, 84, 84, 84},              // region 1
@@ -566,7 +575,7 @@ class Mapper
     {110, 110, 112, 112, 112, 114, 114, 114, 116, 116, 116, 118, 118, 118},    // region 7
     {118, 118, 120, 120, 122, 122, 124, 124, 124, 126, 126, 128, 128},         // region 8
     {128, 130, 130, 132, 132, 132, 134, 134, 136, 136, 138, 138}               // region 9
-  };                                                                           ///< number of pads per row in region
+  }; ///< number of pads per row in region
   static constexpr unsigned int OFFSETCRUGLOBAL[PADROWS]{
     0, 66, 132, 198, 266, 334, 402, 472, 542, 612, 684, 756, 828, 902, 976, 1050, 1124,         // region 0
     0, 76, 152, 228, 306, 384, 462, 542, 622, 702, 784, 866, 948, 1032, 1116,                   // region 1
@@ -578,7 +587,7 @@ class Mapper
     0, 110, 220, 332, 444, 556, 670, 784, 898, 1014, 1130, 1246, 1364, 1482,                    // region 7
     0, 118, 236, 356, 476, 598, 720, 844, 968, 1092, 1218, 1344, 1472,                          // region 8
     0, 128, 258, 388, 520, 652, 784, 918, 1052, 1188, 1324, 1462                                // region 9
-  };                                                                                            ///< row offset in cru for given global pad row
+  }; ///< row offset in cru for given global pad row
 
   static constexpr unsigned int LinksPerRegionPerEndpoint[NREGIONS][NENDPOINTS]{
     {8, 7},   // region 0
@@ -591,7 +600,7 @@ class Mapper
     {10, 10}, // region 7
     {10, 10}, // region 8
     {10, 10}, // region 9
-  };          ///< number of links per region per end point
+  }; ///< number of links per region per end point
 
  private:
   Mapper(const std::string& mappingDir);

Detectors/TPC/base/src/Mapper.cxx

@@ -298,5 +298,41 @@ void Mapper::setTraceLengths(std::string_view inputFile, std::vector<float>& len
   }
 }
 
+bool Mapper::isEdgePad(int rowInSector, int padInRow)
+{
+  const auto& mapper = instance();
+  return (padInRow == 0) || (padInRow == mapper.getNumberOfPadsInRowSector(rowInSector) - 1);
+}
+
+bool Mapper::isFirstOrLastRowInStack(int rowInSector)
+{
+  if (rowInSector == 0 || rowInSector == PADROWS - 1) {
+    return true;
+  }
+
+  const auto& mapper = instance();
+  for (int i = 1; i < 4; ++i) {
+    if (rowInSector == ROWOFFSETSTACK[i] || rowInSector == ROWOFFSETSTACK[i] - 1) {
+      return true;
+    }
+  }
+  return false;
+}
+
+bool Mapper::isBelowSpacerCross(int rowInSector, int padInRow)
+{
+  static std::vector<bool> ROWSBELOWCROSS{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+  if (ROWSBELOWCROSS[rowInSector]) {
+    return true;
+  }
+
+  const auto& mapper = instance();
+  const auto padCenter = mapper.getNumberOfPadsInRowSector(rowInSector) / 2;
+  if (padInRow == padCenter || padInRow == padCenter - 1) {
+    return true;
+  }
+  return false;
+}
+
 } // namespace tpc
 } // namespace o2