mcr/calib/tofino/PrintTofCalib.C

/*
======================================================================
  PrintTofCalib(const TString &fname, const TString &output)
   Dump the calibration results from a root file produced by FillTofCalib.
======================================================================
*/
#include "TFndHdt.h"

using namespace FROOT;

Int_t start_run, stop_run;
Int_t start_time, stop_time;
Float_t ped_e[K_N_TOFINO_SLABS];
Float_t ped_p[K_N_TOFINO_SLABS];
Float_t bhapeak_e[K_N_TOFINO_SLABS];
Float_t bhapeak_p[K_N_TOFINO_SLABS];
Float_t eq_e[K_N_TOFINO_SLABS];
Float_t eq_p[K_N_TOFINO_SLABS];

void PrintTofCalib(const TString &fname, const TString &output){
  if (fname==output){
    cout << ".root file and output file must be different!" << endl;
    return;
  }
  TFile *file = new TFile(fname.Data());
  if (!file) return;

  std::fstream cFstrm(output.Data(), std::ios::out);
  if (cFstrm.fail()) return;

  TTree *tofcalib;
  tofcalib=(TTree*) file->Get("CalibTree");

  tofcalib->SetBranchAddress("StartRun", &start_run);
  tofcalib->SetBranchAddress("StartTime", &start_time);
  tofcalib->SetBranchAddress("StopRun", &stop_run);
  tofcalib->SetBranchAddress("StopTime", &stop_time);
  tofcalib->SetBranchAddress("PedE", ped_e);
  tofcalib->SetBranchAddress("PedP", ped_p);
  tofcalib->SetBranchAddress("EqE", eq_e);
  tofcalib->SetBranchAddress("EqP", eq_p);

  tofcalib->GetEntry(0);


//    for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
//      TThread::Printf("%f ", ped_e[sl]);
//    }
//    for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
//      TThread::Printf("%f ", ped_p[sl]);
//    }
//    for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
//      TThread::Printf("%f ", eq_e[sl]);
//    }
//    for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
//      TThread::Printf("%f ", eq_p[sl]);
//    }

  TDatime start_time_dt, stop_time_dt;
  start_time_dt.Set(start_time-60);
  stop_time_dt.Set(stop_time);

  TString start_time_st=start_time_dt.AsString();
  TString stop_time_st=stop_time_dt.AsString();
 
  cFstrm << start_run << " " << start_time-60 << " \"" << start_time_st << "\" "
         << stop_run << " " << stop_time << " \"" << stop_time_st << "\"";

  cout << start_run << " " << start_time-60 << " \"" << start_time_st << "\" "
       << stop_run << " " << stop_time << " \"" << stop_time_st << "\"";

 for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
    cFstrm << " " << ped_e[sl];
    cout << " " << ped_e[sl];
  }
  for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
    cFstrm << " " << ped_p[sl];
    cout << " " << ped_p[sl];
  }
  for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
    cFstrm << " " << eq_e[sl];
    cout << " " << eq_e[sl];
  }
  for (Int_t sl=0; sl<K_N_TOFINO_SLABS; sl++){
    cFstrm << " " << eq_p[sl];
    cout << " " << eq_p[sl];
  }
  cFstrm << endl;
  cout << endl;
}

void PrintTofCalib(const TString &fname_part){
  TString fname, output;
  fname=fname_part+".root";
  output=fname_part+".txt";
  PrintTofCalib(fname, output);
}

---