CaToolsDebugger.cxx

Go to the documentation of this file.

/* Copyright (C) 2022 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergey Gorbunov, Sergei Zharko [committer] */

 
/*
root -l CAdebug.root
TNtuple *n = (TNtuple*) gDirectory->FindObjectAny("triplets")
n->Draw("chi2")
*/
 
#include "CaToolsDebugger.h"
 
#include "AlgoFairloggerCompat.h"
#include "TDirectory.h"
#include "TFile.h"
#include "TNtuple.h"
 
#include <cstdarg>
#include <iostream>
#include <memory>
#include <string>
 
using cbm::ca::tools::Debugger;
 
class DebuggerInitialiser {
 public:
  DebuggerInitialiser()
  {
    std::shared_ptr<Debugger> instance(new Debugger("CaDebug.root"));
    Debugger::SetInstance(instance);
  }
};
 
static DebuggerInitialiser debuggerInitialiser;
 
 
void Debugger::Write()
{
  fIsWritten          = true;
  TDirectory* currDir = gDirectory;
  if (fFile) {
    fFile->cd();
    for (unsigned int i = 0; i < fNtuples.size(); i++) {
      fNtuples[i]->Write();
    }
  }
  gDirectory = currDir;
}
 
Debugger::~Debugger()
{
  if (!fIsWritten && fFile) {
    LOG(error) << "CaDebugger: you forgot to Write()!!";
  }
}
 
void Debugger::AddNtuple(const char* name, const char* varlist)
{
 
  if (!fFile) {
    fFile = new TFile(fFileName.c_str(), "RECREATE");
  }
  assert(fFile);
  if (GetNtupleIndex(name) >= 0) return;
  TDirectory* currDir = gDirectory;
  fFile->cd();
  fNtuples.push_back(new TNtuple(name, name, varlist));
  gDirectory = currDir;
}
 
int Debugger::GetNtupleIndex(const char* name)
{
  int ind = -1;
  for (unsigned int i = 0; i < fNtuples.size(); i++) {
    if (strcmp(name, fNtuples[i]->GetName()) == 0) {
      ind = i;
      break;
    }
  }
  return ind;
}
 
void Debugger::FillNtuple(const char* name, float v[])
{
  int iNtuple = GetNtupleIndex(name);
  if (iNtuple < 0) {
    std::cerr << "CaDebugger: Ntuple (" << name << ") doesn't exist" << std::endl;
  }
  fNtuples[iNtuple]->Fill(v);
}