masked halfchambers for the HC map (#1820)

* QCtasks multiple fixes

* QCtasks multiple fixes_1

* maskedchambers

* removedthez-lab

Author: maroozm

Modules/TRD/include/TRD/TrackletsTask.h

@@ -71,7 +71,7 @@ class TrackletsTask final : public TaskInterface
   TH1F* mTrackletPosition = nullptr;
   TH1F* mTrackletPositionRaw = nullptr;
   TH1F* mTrackletsPerEvent = nullptr;
-  TH2F* mTrackletsPerEvent2D = nullptr;
+  TH2F* mTrackletsPerHC2D = nullptr;
   std::array<TH2F*, 18> moHCMCMn;
   TH1F* mTrackletSlopen = nullptr;
   TH1F* mTrackletSlopeRawn = nullptr;

Modules/TRD/src/TrackletsTask.cxx

@@ -148,29 +148,57 @@ void TrackletsTask::buildHistograms()
   getObjectsManager()->startPublishing(mTrackletPositionRaw);
   mTrackletsPerEvent = new TH1F("trackletsperevent", "Number of Tracklets per event;Tracklets in Event;Counts", 25000, 0, 25000);
   getObjectsManager()->startPublishing(mTrackletsPerEvent);
-  mTrackletsPerEvent2D = new TH2F("trackletsperevent2d", "Tracklets distribution in half-chambers;Sector_Side;Stack_Side;Number of tracklets per half-chamber", 36, 0, 36, 30, 0, 30);
-  getObjectsManager()->startPublishing(mTrackletsPerEvent2D);
-  getObjectsManager()->setDefaultDrawOptions("trackletsperevent2d", "COLZ");
-  getObjectsManager()->setDisplayHint(mTrackletsPerEvent2D->GetName(), "logz");
-  mTrackletsPerEvent2D->GetXaxis()->SetTitle("Sector_Side");
-  mTrackletsPerEvent2D->GetXaxis()->CenterTitle(kTRUE);
-  mTrackletsPerEvent2D->GetYaxis()->SetTitle("Stack_Layer");
-  mTrackletsPerEvent2D->GetYaxis()->CenterTitle(kTRUE);
+  mTrackletsPerHC2D = new TH2F("trackletsperHC2D", "Tracklets distribution in half-chambers;Sector_Side;Stack_Side", 36, 0, 36, 30, 0, 30);
+  mTrackletsPerHC2D->SetStats(0);
+  mTrackletsPerHC2D->GetXaxis()->SetTitle("Sector_Side");
+  mTrackletsPerHC2D->GetXaxis()->CenterTitle(kTRUE);
+  mTrackletsPerHC2D->GetYaxis()->SetTitle("Stack_Layer");
+  mTrackletsPerHC2D->GetYaxis()->CenterTitle(kTRUE);
   for (int s = 0; s < o2::trd::constants::NSTACK; ++s) {
     for (int l = 0; l < o2::trd::constants::NLAYER; ++l) {
       std::string label = fmt::format("{0}_{1}", s, l);
       int pos = s * o2::trd::constants::NLAYER + l + 1;
-      mTrackletsPerEvent2D->GetYaxis()->SetBinLabel(pos, label.c_str());
+      mTrackletsPerHC2D->GetYaxis()->SetBinLabel(pos, label.c_str());
     }
   }
   for (int sm = 0; sm < o2::trd::constants::NSECTOR; ++sm) {
     for (int side = 0; side < 2; ++side) {
       std::string label = fmt::format("{0}_{1}", sm, side == 0 ? "A" : "B");
       int pos = sm * 2 + side + 1;
-      mTrackletsPerEvent2D->GetXaxis()->SetBinLabel(pos, label.c_str());
+      mTrackletsPerHC2D->GetXaxis()->SetBinLabel(pos, label.c_str());
     }
   }
 
+  // Mask known inactive halfchambers in the active chamber map
+  TLine* line[6];
+  std::pair<int, int> x, y;
+  for (int iHC = 0; iHC < o2::trd::constants::NCHAMBER * 2; ++iHC) {
+    if (mChamberStatus != nullptr) {
+      if (mChamberStatus->isMasked(iHC)) {
+        int stackLayer = o2::trd::HelperMethods::getStack(iHC / 2) * o2::trd::constants::NLAYER + o2::trd::HelperMethods::getLayer(iHC / 2);
+        int sectorSide = (iHC / o2::trd::constants::NHCPERSEC) * 2 + (iHC % 2);
+        x.first = sectorSide;
+        x.second = sectorSide + 1;
+        y.first = stackLayer;
+        y.second = stackLayer + 1;
+
+        line[0] = new TLine(x.first, y.first, x.second, y.second);
+        line[1] = new TLine(x.second, y.first, x.first, y.second);
+        line[2] = new TLine(x.first, y.first, x.second, y.first);
+        line[3] = new TLine(x.first, y.second, x.second, y.second);
+        line[4] = new TLine(x.first, y.first, x.first, y.second);
+        line[5] = new TLine(x.second, y.first, x.second, y.second);
+        for (int i = 0; i < 6; ++i) {
+          line[i]->SetLineColor(kBlack);
+          mTrackletsPerHC2D->GetListOfFunctions()->Add(line[i]);
+        }
+      }
+    }
+  }
+  getObjectsManager()->startPublishing(mTrackletsPerHC2D);
+  getObjectsManager()->setDefaultDrawOptions("trackletsperHC2D", "COLZ");
+  getObjectsManager()->setDisplayHint(mTrackletsPerHC2D->GetName(), "logz");
+
   for (Int_t sm = 0; sm < o2::trd::constants::NSECTOR; ++sm) {
     std::string label = fmt::format("TrackletHCMCMnoise_{0}", sm);
     std::string title = fmt::format("MCM in Tracklets data stream for sector {0} noise in;mcm in rob in layer;ROB in stack", sm);
@@ -413,7 +441,7 @@ void TrackletsTask::monitorData(o2::framework::ProcessingContext& ctx)
           int hcid = tracklets[currenttracklet].getHCID();
           int stackLayer = o2::trd::HelperMethods::getStack(hcid / 2) * o2::trd::constants::NLAYER + o2::trd::HelperMethods::getLayer(hcid / 2);
           int sectorSide = (hcid / o2::trd::constants::NHCPERSEC) * 2 + (hcid % 2);
-          mTrackletsPerEvent2D->Fill(sectorSide, stackLayer);
+          mTrackletsPerHC2D->Fill(sectorSide, stackLayer);
           // 8 rob x 16 mcm each per chamber
           //  5 stack(y), 6 layers(x)
           //  y=stack_rob, x=layer_mcm
@@ -488,7 +516,7 @@ void TrackletsTask::reset()
   mTrackletPosition->Reset();
   mTrackletPositionRaw->Reset();
   mTrackletsPerEvent->Reset();
-  mTrackletsPerEvent2D->Reset();
+  mTrackletsPerHC2D->Reset();
   for (auto h : moHCMCMn) {
     h->Reset();
   }