plots to check dead dcols or regions at SOR and EOR

Author: iravasen

Modules/ITS/include/ITS/ITSClusterTask.h

@@ -74,18 +74,20 @@ class ITSClusterTask : public TaskInterface
 
   static constexpr int NLayer = 7;
   static constexpr int NLayerIB = 3;
+  static constexpr int NLayerOB = 4;
   static constexpr int NStavesIB = 12 + 16 + 20;
   static constexpr int NStavesOB = 24 + 30 + 42 + 48;
 
   std::vector<TObject*> mPublishedObjects;
-
   // Task
   TH1D* hTFCounter = nullptr;
 
   // Inner barrel
   TH1D* hClusterTopologySummaryIB[NLayer][48][9] = { { { nullptr } } };
   TH1D* hGroupedClusterSizeSummaryIB[NLayer][48][9] = { { { nullptr } } };
   TH1D* hClusterSizeSummaryIB[NLayer][48][9] = { { { nullptr } } };
+  TH2I* hDcolMapIBSor = nullptr;
+  TH2I* hDcolMapIBEor = nullptr;
 
   std::shared_ptr<TH2DRatio> hAverageClusterOccupancySummaryIB[NLayer];
   std::shared_ptr<TH2DRatio> hAverageClusterSizeSummaryIB[NLayer];
@@ -98,6 +100,9 @@ class ITSClusterTask : public TaskInterface
   std::shared_ptr<TH2DRatio> hAverageClusterOccupancySummaryOB[NLayer];
   std::shared_ptr<TH2DRatio> hAverageClusterSizeSummaryOB[NLayer];
 
+  TH2I* hDcolMapOBSor[NLayerOB] = { nullptr };
+  TH2I* hDcolMapOBEor[NLayerOB] = { nullptr };
+
   // Layer summary
   TH1L* hClusterSizeLayerSummary[NLayer] = { nullptr };
   TH1L* hClusterTopologyLayerSummary[NLayer] = { nullptr };
@@ -135,9 +140,13 @@ class ITSClusterTask : public TaskInterface
   int nRphiBinsOB = 1;
   int nZBinsOB = 1;
   static constexpr int NFlags = 4;
+  int nResetCycle = 3; // to reset plots every nResetCycle cycles
+  int nCycle = 0;      // cycle counter
+  int nCycleStop = 3;  // number of cycles after which to stop filling histos
 
   const int mOccUpdateFrequency = 100000;
-  const int mNLanes[4] = { 432, 864, 2520, 3816 }; // IB, ML, OL, TOTAL lane
+  const int mNLanes[4] = { 432, 864, 2520, 3816 };                      // IB, ML, OL, TOTAL lane
+  const int mChips[NLayer] = { 108, 144, 180, 2688, 3360, 8232, 9408 }; // # chips in L0, L2, L2, L3, L4, L5, L6
   int mDoPublish1DSummary = 0;
   int mNThreads = 1;
   int nBCbins = 103;

Modules/ITS/src/ITSClusterTask.cxx

@@ -40,6 +40,8 @@ ITSClusterTask::ITSClusterTask() : TaskInterface() {}
 
 ITSClusterTask::~ITSClusterTask()
 {
+  delete hDcolMapIBSor;
+  delete hDcolMapIBEor;
   delete hTFCounter;
   delete hEmptyLaneFractionGlobal;
   delete hClusterVsBunchCrossing;
@@ -53,6 +55,8 @@ ITSClusterTask::~ITSClusterTask()
       delete hClusterCenterMap[iLayer];
       delete hLongClustersPerChip[iLayer];
       delete hMultPerChipWhenLongClusters[iLayer];
+      delete hDcolMapOBSor[iLayer - NLayerIB];
+      delete hDcolMapOBEor[iLayer - NLayerIB];
     }
 
     else {
@@ -123,6 +127,15 @@ void ITSClusterTask::startOfActivity(const Activity& /*activity*/)
 void ITSClusterTask::startOfCycle()
 {
   ILOG(Debug, Devel) << "startOfCycle" << ENDM;
+
+  if (nCycle % nResetCycle == 0) {
+    hDcolMapIBEor->Reset();
+    for (int i = 0; i < NLayerOB; i++) {
+      hDcolMapOBEor[i]->Reset();
+    }
+  }
+
+  nCycle++;
 }
 
 void ITSClusterTask::monitorData(o2::framework::ProcessingContext& ctx)
@@ -248,6 +261,12 @@ void ITSClusterTask::monitorData(o2::framework::ProcessingContext& ctx)
         hAverageClusterSizeSummaryIB[lay]->getNum()->Fill(chip, sta, (double)npix);
         hAverageClusterSizeSummaryIB[lay]->getDen()->Fill(chip, sta, 1.);
         hClusterCenterMap[lay]->Fill(cluster.getCol(), cluster.getRow());
+
+        if (nCycle < nCycleStop) {
+          hDcolMapIBSor->Fill((int)(cluster.getCol() / 2), ChipID);
+        }
+        hDcolMapIBEor->Fill((int)(cluster.getCol() / 2), ChipID);
+
         if (mDoPublish1DSummary == 1) {
           hClusterTopologySummaryIB[lay][sta][chip]->Fill(ClusterID);
         }
@@ -266,6 +285,12 @@ void ITSClusterTask::monitorData(o2::framework::ProcessingContext& ctx)
           hGroupedClusterSizeLayerSummary[lay]->Fill(npix);
         }
       } else {
+
+        if (nCycle < nCycleStop) {
+          hDcolMapOBSor[lay - NLayerIB]->Fill((int)(cluster.getCol() / 32), ChipID);
+        }
+        hDcolMapOBEor[lay - NLayerIB]->Fill((int)(cluster.getCol() / 32), ChipID);
+
         hAverageClusterOccupancySummaryOB[lay]->getNum()->Fill(lane, sta, 1. / (mNChipsPerHic[lay] / mNLanePerHic[lay])); // 14 To have occupation per chip -> 7 because we're considering lanes
         hAverageClusterSizeSummaryOB[lay]->getNum()->Fill(lane, sta, (double)npix);
         hAverageClusterSizeSummaryOB[lay]->getDen()->Fill(lane, sta, 1);
@@ -435,6 +460,7 @@ void ITSClusterTask::reset()
 {
   ILOG(Debug, Devel) << "Resetting the histograms" << ENDM;
   hTFCounter->Reset();
+  hDcolMapIBSor->Reset();
   hClusterVsBunchCrossing->Reset();
   hEmptyLaneFractionGlobal->Reset("ICES");
   mGeneralOccupancy->Reset();
@@ -463,6 +489,7 @@ void ITSClusterTask::reset()
         }
       }
     } else {
+      hDcolMapOBSor[iLayer - NLayerIB]->Reset();
       hLongClustersPerStave[iLayer - NLayerIB]->Reset();
       hAverageClusterOccupancySummaryOB[iLayer]->Reset();
       hAverageClusterSizeSummaryOB[iLayer]->Reset();
@@ -485,6 +512,12 @@ void ITSClusterTask::reset()
 
 void ITSClusterTask::createAllHistos()
 {
+  hDcolMapIBSor = new TH2I("hDcolMapIBSor", "Double column hitmap IB chips - SOR;Dcol;ChipID", 512, -0.5, 511.5, mNLanes[0], -0.5, mNLanes[0] - 0.5);
+  addObject(hDcolMapIBSor);
+
+  hDcolMapIBEor = new TH2I("hDcolMapIBEor", "Double column hitmap IB chips - EOR;Dcol;ChipID", 512, -0.5, 511.5, mNLanes[0], -0.5, mNLanes[0] - 0.5);
+  addObject(hDcolMapIBEor);
+
   hTFCounter = new TH1D("TFcounter", "TFcounter", 1, 0, 1);
   hTFCounter->SetTitle("TF counter");
   addObject(hTFCounter);
@@ -535,6 +568,19 @@ void ITSClusterTask::createAllHistos()
       addObject(hLongClustersPerStave[iLayer - NLayerIB]);
       formatAxes(hLongClustersPerStave[iLayer - NLayerIB], "Stave", "number of long clusters", 1, 1.10);
       hLongClustersPerStave[iLayer - NLayerIB]->SetStats(0);
+
+      // dcols maps
+      int chipStart = 0, chipStop = 0;
+      for (int j = 0; j < iLayer; j++) {
+        chipStart += mChips[j];
+      }
+      chipStop = chipStart + mChips[iLayer];
+
+      hDcolMapOBSor[iLayer - NLayerIB] = new TH2I(Form("hDcolMapL%dSor", iLayer), Form("Double column hitmap L%d chips - SOR;RegionID;ChipID", iLayer), 32, -0.5, 31.5, mChips[iLayer], chipStart - 0.5, chipStop - 0.5);
+      addObject(hDcolMapOBSor[iLayer - NLayerIB]);
+
+      hDcolMapOBEor[iLayer - NLayerIB] = new TH2I(Form("hDcolMapL%dEor", iLayer), Form("Double column hitmap L%d chips - EOR;RegionID;ChipID", iLayer), 32, -0.5, 31.5, mChips[iLayer], chipStart - 0.5, chipStop - 0.5);
+      addObject(hDcolMapOBEor[iLayer - NLayerIB]);
     }
 
     hClusterSizeLayerSummary[iLayer] = new TH1L(Form("Layer%d/AverageClusterSizeSummary", iLayer), Form("Layer%dAverageClusterSizeSummary", iLayer), 128 * 128, 0, 128 * 128);
@@ -712,6 +758,9 @@ void ITSClusterTask::getJsonParameters()
   std::string LayerConfig = o2::quality_control_modules::common::getFromConfig<std::string>(mCustomParameters, "layer", "0000000");
   mDoPublishDetailedSummary = o2::quality_control_modules::common::getFromConfig<int>(mCustomParameters, "publishDetailedSummary", mDoPublishDetailedSummary);
 
+  nCycleStop = o2::quality_control_modules::common::getFromConfig<int>(mCustomParameters, "nCycleStop", nCycleStop);
+  nResetCycle = o2::quality_control_modules::common::getFromConfig<int>(mCustomParameters, "nResetCycle", nResetCycle);
+
   for (int ilayer = 0; ilayer < NLayer; ilayer++) {
     if (LayerConfig[ilayer] != '0') {
       mEnableLayers[ilayer] = 1;