TPC splines: add possibility to merge specific sector

Author: cbmsw

Detectors/TPC/calibration/include/TPCCalibration/TPCFastSpaceChargeCorrectionHelper.h

@@ -119,11 +119,11 @@ class TPCFastSpaceChargeCorrectionHelper
   /// \param mainCorrection main correction
   /// \param scale scaling factor for the main correction
   /// \param additionalCorrections vector of pairs of additional corrections and their scaling factors
-  /// \param prn printout flag
+  /// \param sector when >=0 perform the merge only for the sector
   /// \return main correction merged with additional corrections
   void mergeCorrections(
     o2::gpu::TPCFastSpaceChargeCorrection& mainCorrection, float scale,
-    const std::vector<std::pair<const o2::gpu::TPCFastSpaceChargeCorrection*, float>>& additionalCorrections, bool prn);
+    const std::vector<std::pair<const o2::gpu::TPCFastSpaceChargeCorrection*, float>>& additionalCorrections, int sector = -1);
 
   /// how far the voxel mean is allowed to be outside of the voxel (1.1 means 10%)
   void setVoxelMeanValidityRange(double range)

Detectors/TPC/calibration/src/TPCFastSpaceChargeCorrectionHelper.cxx

@@ -1032,7 +1032,7 @@ void TPCFastSpaceChargeCorrectionHelper::initInverse(std::vector<o2::gpu::TPCFas
 
 void TPCFastSpaceChargeCorrectionHelper::mergeCorrections(
   o2::gpu::TPCFastSpaceChargeCorrection& mainCorrection, float mainScale,
-  const std::vector<std::pair<const o2::gpu::TPCFastSpaceChargeCorrection*, float>>& additionalCorrections, bool /*prn*/)
+  const std::vector<std::pair<const o2::gpu::TPCFastSpaceChargeCorrection*, float>>& additionalCorrections, int SelectedSector)
 {
   /// merge several corrections
 
@@ -1041,7 +1041,15 @@ void TPCFastSpaceChargeCorrectionHelper::mergeCorrections(
 
   const auto& geo = mainCorrection.getGeometry();
 
-  for (int sector = 0; sector < geo.getNumberOfSectors(); sector++) {
+  int firstSector = 0;
+  int lastSector = geo.getNumberOfSectors() - 1;
+  if (SelectedSector >= 0) {
+    assert(SelectedSector < geo.getNumberOfSectors());
+    firstSector = SelectedSector;
+    lastSector = SelectedSector;
+  }
+
+  for (int sector = firstSector; sector <= lastSector; sector++) {
 
     auto myThread = [&](int iThread) {
       for (int row = iThread; row < geo.getNumberOfRows(); row += mNthreads) {