GLO: Add mtc pt trending

Signed-off-by: Felix Schlepper <felix.schlepper@cern.ch>

Author: f3sch

Modules/GLO/glo-itstpc-mtch-qcmn-test.json

@@ -49,6 +49,7 @@
           "GID": "ITS-TPC,ITS",
           "verbose": "false",
           "doMTCRatios": "true",
+          "doMTCTrending": "true",
           "doK0QC": "true",
           "isSync": "true",
           "trackSourcesK0": "ITS,TPC,ITS-TPC,ITS-TPC-TOF,TPC-TOF,TPC-TRD,ITS-TPC-TRD,TPC-TRD-TOF,ITS-TPC-TOF,ITS-TPC-TRD-TOF"
@@ -131,16 +132,55 @@
             "title": "Fitted K0s mass",
             "varexp": "gloFitK0sMassSignal.mean:time",
             "option": "*L",
+            "graphYRange": "0.45:0.55",
             "graphAxisLabel": "K0s mass (GeV/c^{2}):time"
           },
           {
             "name": "glo_k0s_sigma_trending",
             "title": "Fitted K0s sigma",
             "varexp": "gloFitK0sMassSignal.sigma:time",
             "option": "*L",
+            "graphYRange": "0:0.01",
             "graphAxisLabel": "#sigma_{K0s mass (GeV/c^{2})}:time"
           }
         ]
+      },
+      "MTCTrend": {
+        "active": "true",
+        "className": "o2::quality_control::postprocessing::TrendingTask",
+        "moduleName": "QcGLO",
+        "detectorName": "GLO",
+        "resumeTrend": "false",
+        "initTrigger": [
+          "once"
+        ],
+        "updateTrigger": [
+          "newobject:qcdb:GLO/MO/ITSTPCMatchingTask/gloMTCTrendingObject"
+        ],
+        "stopTrigger": [
+          "usercontrol"
+        ],
+        "dataSources": [
+          {
+            "type": "repository",
+            "path": "GLO/MO/ITSTPCMatchingTask",
+            "names": [
+              "gloMTCTrendingObject"
+            ],
+            "reductorName": "o2::quality_control_modules::glo::MTCReductor",
+            "moduleName": "QcGLO"
+          }
+        ],
+        "plots": [
+          {
+            "name": "glo_mtc_pt_trending",
+            "title": "MTC #it{p}_{T} trending",
+            "varexp": "gloMTCTrendingObject.pt:time",
+            "option": "*L",
+            "graphYRange": "0:1.02",
+            "graphAxisLabel": "MTC:time"
+          }
+        ]
       }
     }
   }

Modules/GLO/include/GLO/ITSTPCMatchingTask.h

@@ -21,6 +21,7 @@
 #include "GLOQC/MatchITSTPCQC.h"
 #include "Common/TH1Ratio.h"
 #include "GLO/Helpers.h"
+#include "GLO/Reductors.h"
 
 #include <CommonConstants/PhysicsConstants.h>
 
@@ -56,6 +57,9 @@ class ITSTPCMatchingTask final : public TaskInterface
   std::unique_ptr<common::TH1FRatio> mEffPt;
   std::unique_ptr<common::TH1FRatio> mEffEta;
   std::unique_ptr<common::TH1FRatio> mEffPhi;
+  bool mDoMTCTrending{ false };
+  float mMTCTrendingPt{ 1.5 };
+  std::unique_ptr<TF1> mMTCTrendingObject;
 
   bool mDoK0s{ false };
   bool mPublishK0s3D{ false };

Modules/GLO/include/GLO/LinkDef.h

@@ -14,4 +14,5 @@
 #pragma link C++ class o2::quality_control_modules::glo::ITSTPCMatchingTask + ;
 #pragma link C++ class o2::quality_control_modules::glo::ITSTPCmatchingCheck + ;
 #pragma link C++ class o2::quality_control_modules::glo::K0sFitReductor + ;
+#pragma link C++ class o2::quality_control_modules::glo::MTCReductor + ;
 #endif

Modules/GLO/include/GLO/Reductors.h

@@ -20,13 +20,25 @@ namespace o2::quality_control_modules::glo
 class K0sFitReductor final : public quality_control::postprocessing::ReductorTObject
 {
   void* getBranchAddress() final { return &mStats; };
-  const char* getBranchLeafList() final { return "mean/D:sigma"; };
-  void update(TObject* obj) override;
+  const char* getBranchLeafList() final { return "mean/F:sigma/F"; };
+  void update(TObject* obj) final;
 
  private:
   struct {
-    Double_t mean;
-    Double_t stddev;
+    Float_t mean;
+    Float_t sigma;
+  } mStats;
+};
+
+class MTCReductor final : public quality_control::postprocessing::ReductorTObject
+{
+  void* getBranchAddress() final { return &mStats; };
+  const char* getBranchLeafList() final { return "pt/F"; };
+  void update(TObject* obj) final;
+
+ private:
+  struct {
+    Float_t pt;
   } mStats;
 };
 

Modules/GLO/src/ITSTPCMatchingTask.cxx

@@ -55,7 +55,12 @@ void ITSTPCMatchingTask::initialize(o2::framework::InitContext& /*ctx*/)
   // Sync
   mIsSync = common::getFromConfig(mCustomParameters, "isSync", false);
   // MTC ratios
-  mDoMTCRatios = common::getFromConfig(mCustomParameters, "doMTCRatios", false);
+  if ((mDoMTCRatios = common::getFromConfig(mCustomParameters, "doMTCRatios", false))) {
+    if ((mDoMTCTrending = getFromConfig(mCustomParameters, "doMTCTrending", false))) {
+      mMTCTrendingPt = getFromConfig(mCustomParameters, "trendingPt", 1.5f);
+      mMTCTrendingObject.reset(new TF1("gloMTCTrendingObject", "pol0"));
+    }
+  }
   // K0s
   mMatchITSTPCQC.setDoK0QC((mDoK0s = getFromConfig(mCustomParameters, "doK0QC", false)));
   if (mIsSync && mDoK0s) {
@@ -131,6 +136,10 @@ void ITSTPCMatchingTask::endOfCycle()
       makeRatio(mEffPt, mMatchITSTPCQC.getFractionITSTPCmatch(gloqc::MatchITSTPCQC::ITS));
       getObjectsManager()->startPublishing(mEffPt.get(), PublicationPolicy::Once);
       getObjectsManager()->setDefaultDrawOptions(mEffPt->GetName(), "logx");
+      if (mDoMTCTrending) {
+        mMTCTrendingObject->SetParameter(0, mEffPt->GetBinContent(mEffPt->FindBin(mMTCTrendingPt)));
+        getObjectsManager()->startPublishing<true>(mMTCTrendingObject.get(), PublicationPolicy::Once);
+      }
 
       // Eta
       makeRatio(mEffEta, mMatchITSTPCQC.getFractionITSTPCmatchEta(gloqc::MatchITSTPCQC::ITS));

Modules/GLO/src/Reductors.cxx

@@ -25,7 +25,17 @@ void K0sFitReductor::update(TObject* obj)
   }
 
   mStats.mean = f->GetParameter(helpers::K0sFitter::Parameters::Mass);
-  mStats.stddev = f->GetParameter(helpers::K0sFitter::Parameters::Sigma);
+  mStats.sigma = f->GetParameter(helpers::K0sFitter::Parameters::Sigma);
+}
+
+void MTCReductor::update(TObject* obj)
+{
+  auto f = dynamic_cast<TF1*>(obj);
+  if (!f) {
+    return;
+  }
+
+  mStats.pt = f->GetParameter(0);
 }
 
 } // namespace o2::quality_control_modules::glo