Writing TOF calibration to CCDB

Author: chiarazampolli

Detectors/GlobalTracking/include/GlobalTracking/CalibTOF.h

@@ -60,7 +60,7 @@ class CalibTOF
 
   ///< calibrate using the provided input
   void run(int flag, int sector = -1);
-  void fillOutput();
+  void fillOutput(int flag);
 
   ///< perform all initializations
   void init();
@@ -102,6 +102,12 @@ class CalibTOF
   float getNsigmaFractionProblematicCut() const { return mNsigmaFractionProblematicCut; }
   float getNsigmaSigmaProblematicCut() const { return mNsigmaSigmaProblematicCut; }
 
+  void setFillCCDB(bool flag) { mFillCCDB = flag; }
+  bool getFillCCDB() const { return mFillCCDB; }
+
+  void setCCDBpath(std::string path) { mCCDBpath = path; }
+  std::string getCCDBpath() const { return mCCDBpath; }
+  
  private:
   Int_t mDebugMode = 0; ///< >0= time slewing extra plot, >1= problematic fits stored
 
@@ -129,7 +135,9 @@ class CalibTOF
   // Data members
 
   bool mInitDone = false; ///< flag init already done
-
+  bool mFillCCDB = false; ///< flag init already doneto decide whether to fill or not the CCDB
+  std::string mCCDBpath = "http://ccdb-test.cern.ch:8080"; ///< path to CCDB 
+  
   ///========== Parameters to be set externally, e.g. from CCDB ====================
 
   // to be done later

Detectors/GlobalTracking/src/CalibTOF.cxx

@@ -16,6 +16,7 @@
 #include <TFile.h>
 #include "DataFormatsParameters/GRPObject.h"
 #include "ReconstructionDataFormats/PID.h"
+#include "CCDB/CcdbApi.h"
 
 #include "GlobalTracking/CalibTOF.h"
 
@@ -25,6 +26,7 @@
 #include "TRandom.h"
 
 using namespace o2::globaltracking;
+using CcdbApi = o2::ccdb::CcdbApi;
 
 ClassImp(CalibTOF);
 
@@ -259,9 +261,23 @@ void CalibTOF::run(int flag, int sector)
 }
 
 //______________________________________________
-void CalibTOF::fillOutput()
+void CalibTOF::fillOutput(int flag)
 {
   mOutputTree->Fill();
+
+  if (mFillCCDB) {
+    CcdbApi api;
+    api.init(mCCDBpath.c_str()); // or http://localhost:8080 for a local installation
+    
+    if (flag & kLHCphase) {
+      std::map<std::string, std::string> metadataLHCphase; // can be empty
+      api.storeAsTFileAny(mLHCphaseObj, "TOF/LHCphase", metadataLHCphase, (ulong)mMinTimestamp*1000, (ulong)mMaxTimestamp*1000); // we use as validity the timestamps that we got from the input for the calibration; but we need to convert to ms for the CCDB (at least for now that we use an integer for the timestamp)
+    }
+    if (flag & kChannelOffset || flag & kChannelTimeSlewing) {
+      std::map<std::string, std::string> metadataChannelCalib; // can be empty
+      api.storeAsTFileAny(mTimeSlewingObj, "TOF/ChannelCalib", metadataChannelCalib, (ulong)mMinTimestamp*1000); // contains both offset and time slewing; we use as validity the START ONLY timestamp that we got from the input for the calibration; but we need to convert to ms for the CCDB (at least for now that we use an integer for the timestamp), END is default
+    }
+  }
 }
 
 //______________________________________________

macro/run_calib_tof.C

@@ -34,14 +34,15 @@
 #include "err.h"
 
 void run_calib_tof(std::string path = "./", std::string outputfile = "o2calparams_tof.root",
-                   std::string inputfileCalib = "o2calibration_tof.root")
+                   std::string inputfileCalib = "o2calibration_tof.root", int calibType = o2::globaltracking::CalibTOF::kLHCphase)
 {
   bool onlymerge = false; // set to true if you have already the outputs from forked processes and you want only merge
 
   TString namefile(outputfile);
   namefile.ReplaceAll(".root", "");
 
-  const int ninstance = 4;
+  int ninstance = 4;
+  if (calibType == o2::globaltracking::CalibTOF::kLHCphase) ninstance = 1;
   o2::globaltracking::CalibTOF calib;
   calib.setDebugMode(1);
 
@@ -71,7 +72,7 @@ void run_calib_tof(std::string path = "./", std::string outputfile = "o2calparam
   // to be generalized for more than 2 forks
   int counter = 0;
 
-  pid_t pids[ninstance];
+  pid_t* pids = new pid_t[ninstance];
   int n = ninstance;
   /* Start children. */
   TStopwatch timerTot;
@@ -91,14 +92,15 @@ void run_calib_tof(std::string path = "./", std::string outputfile = "o2calparam
         cout << "Child process: Value returned by fork() = " << pids[i] << endl;
 
         // only for the first child
-        if (i == 0)
+        if (i == 0 && calibType & o2::globaltracking::CalibTOF::kLHCphase) 
           calib.run(o2::globaltracking::CalibTOF::kLHCphase);
-        for (int sect = i; sect < 18; sect += ninstance)
-          //calib.run(o2::globaltracking::CalibTOF::kChannelTimeSlewing, sect);
-          calib.run(o2::globaltracking::CalibTOF::kChannelOffset, sect);
-        calib.fillOutput();
-        outFile.cd();
-        outTree.Write();
+	if (calibType > o2::globaltracking::CalibTOF::kLHCphase) { // we do more than only the LHCphase
+	  for (int sect = i; sect < 18; sect += ninstance)
+	    calib.run((~o2::globaltracking::CalibTOF::kLHCphase) & calibType, sect); // we switch off the LHCphase, if it was enabled in the calibType
+	}
+	calib.fillOutput();
+	outFile.cd();
+	outTree.Write();
         if (i == 0)
           calib.getLHCphaseHisto()->Write();
         calib.getChTimeSlewingHistoAll()->Write();
@@ -117,6 +119,8 @@ void run_calib_tof(std::string path = "./", std::string outputfile = "o2calparam
     }
   }
 
+  delete [] pids;
+  
   timerTot.Stop();
   Printf("Time to run the calibration was:");
   timerTot.Print();
@@ -142,7 +146,7 @@ void run_calib_tof(std::string path = "./", std::string outputfile = "o2calparam
   timerTot.Print();
 
   outFile.cd();
-  calib.fillOutput();
+  calib.fillOutput(calibType);
   outTree.Write();
   outFile.Close();
 }