exb-trend

Author: maroozm

Diff for: Modules/TRD/calibParamTRD.json

@@ -6,13 +6,13 @@
         "host": "ccdb-test.cern.ch:8080"
       },
       "Activity": {
-        "number": ""
+        "number": "526486"
       },
       "conditionDB": {
         "url": "alice-ccdb.cern.ch"
       },
       "postprocessing": {
-      "matchAnyRunNumber": "true"
+        "matchAnyRunNumber": "true"
       }
     },
     "postprocessing": {
@@ -26,34 +26,57 @@
           {
             "type": "condition",
             "path": "TRD/Calib",
-            "names": [ "CalVdriftExB" ],
+            "names": [
+              "CalVdriftExB"
+            ],
             "reductorName": "o2::quality_control_modules::trd::CalibReductorTRD",
             "moduleName": "QcTRD"
           }
         ],
         "plots": [
           {
             "name": "vdrift_trend",
-            "title": "vDrift trend",
+            "title": "vDrift mean trend",
             "varexp": "CalVdriftExB.vdriftmean:time",
             "selection": "",
             "option": "*L",
             "graphYRange": "0:3.5",
             "graphXRange": "",
             "graphAxisLabel": "vdrift:time",
             "graphErrors": "0:CalVdriftExB.vdrifterr"
+          },
+          {
+            "name": "vdrift_trend2",
+            "title": "vDrift trend",
+            "varexp": "CalVdriftExB.vdrift[]:runNumber",
+            "selection": "",
+            "option": "*L",
+            "graphYRange": "-2:3.5",
+            "graphXRange": "",
+            "graphAxisLabel": "vdrift:time"
+          },
+          {
+            "name": "exb_trend",
+            "title": "ExB mean trend",
+            "varexp": "CalVdriftExB.exbmean:time",
+            "selection": "",
+            "option": "*L",
+            "graphYRange": "-2.0:2.0",
+            "graphXRange": "",
+            "graphAxisLabel": "ExB:time",
+            "graphErrors": "0:CalVdriftExB.exberr"
           }
         ],
         "initTrigger": [
           "userorcontrol"
         ],
         "updateTrigger": [
-          "foreachobject:ccdb:TRD/Calib/CalVdriftExB"
+          "newobject:ccdb:TRD/Calib/CalVdriftExB"
         ],
         "stopTrigger": [
           "userorcontrol"
         ]
       }
     }
   }
-}
+}

Diff for: Modules/TRD/include/TRD/CalibReductorTRD.h

@@ -11,7 +11,7 @@
 
 ///
 /// \file   CalibReductorTRD.h
-/// \author
+/// \author Salman Malik
 ///
 
 #ifndef QUALITYCONTROL_CALIBREDUCTORTRD_H
@@ -38,8 +38,11 @@ class CalibReductorTRD : public quality_control::postprocessing::ReductorConditi
 
  private:
   struct {
+    Float_t vdrift[o2::trd::constants::MAXCHAMBER];
     Float_t vdriftmean;
     Float_t vdrifterr;
+    Float_t exbmean;
+    Float_t exberr;
   } mStats;
 };
 

Diff for: Modules/TRD/src/CalibReductorTRD.cxx

@@ -11,7 +11,7 @@
 
 ///
 /// \file   CalibReductorTRD.cxx
-/// \author
+/// \author Salman Malik
 ///
 
 #include "TRD/CalibReductorTRD.h"
@@ -28,21 +28,29 @@ void* CalibReductorTRD::getBranchAddress()
 
 const char* CalibReductorTRD::getBranchLeafList()
 {
-  return "vdriftmean/F:vdrifterr/F";
+  return "vdrift[540]/F:vdriftmean/F:vdrifterr/F:exbmean/F:exberr/F";
 }
 
 bool CalibReductorTRD::update(ConditionRetriever& retriever)
 {
   if (auto retvdrift = retriever.retrieve<o2::trd::CalVdriftExB>()) {
-    double sum = 0;
-    double sumSq = 0;
+    double sumVdrift = 0, sumSqVdrift = 0;
+    double sumExb = 0, sumSqExb = 0;
+
     for (int i = 0; i < o2::trd::constants::MAXCHAMBER; i++) {
-      double vdrift = retvdrift->getVdrift(i);
-      sum += vdrift;
-      sumSq += vdrift * vdrift;
+      mStats.vdrift[i] = retvdrift->getVdrift(i);
+      sumVdrift += retvdrift->getVdrift(i);
+      sumSqVdrift += retvdrift->getVdrift(i) * retvdrift->getVdrift(i);
+      sumExb += retvdrift->getExB(i);
+      sumSqExb += retvdrift->getExB(i) * retvdrift->getExB(i);
     }
-    mStats.vdriftmean = sum / o2::trd::constants::MAXCHAMBER;
-    mStats.vdrifterr = std::sqrt(sumSq / o2::trd::constants::MAXCHAMBER - mStats.vdriftmean * mStats.vdriftmean);
+
+    mStats.vdriftmean = sumVdrift / o2::trd::constants::MAXCHAMBER;
+    mStats.vdrifterr = std::sqrt(sumSqVdrift / o2::trd::constants::MAXCHAMBER - mStats.vdriftmean * mStats.vdriftmean);
+
+    mStats.exbmean = sumExb / o2::trd::constants::MAXCHAMBER;
+    mStats.exberr = std::sqrt(sumSqExb / o2::trd::constants::MAXCHAMBER - mStats.exbmean * mStats.exbmean);
+
     return true;
   }
   return false;