readtofTotbCorr.C

// macro to read tofr5 MAP table from database
//
// based on
//  http://www.star.bnl.gov/STAR/comp/db/StoreDbTable.cc.html
//
// Xin Dong, 04/14/2010
//


// #include "StDbLib/StDbManager.hh"
// #include "StDbLib/StDbConfigNode.hh"
// #include "StDbLib/StDbTable.h"
// #include "StDbLib/StDbDefs.hh"

#include <iostream>
#include <fstream>
#include <string>
#include "iomanip.h"
using namespace std;

void readtofTotbCorr(string ZReadTime = "2029-12-31 23:59:59")
{
  const int mNTray = 120;
  const int mNTDIG = 8;
  const int mNVPD = 19;

  //-- load dBase and Table definition libraries
  gSystem->Load("St_base");
  gSystem->Load("StChain");
  gSystem->Load("StUtilities");
  gSystem->Load("St_Tables.so");

  gSystem->Load("StDbLib.so");
  gSystem->Load("libStDb_Tables.so");

  //-- get the singleton manager
  StDbManager* dbManager = StDbManager::Instance();

  //-- connect to the db & get an empty container
  StDbConfigNode* configNode = dbManager->initConfig("Calibrations_tof");

  dbManager->setRequestTime(ZReadTime.c_str());

  StDbTable* tofTotbCorr = configNode->addDbTable("tofTotbCorr");

  dbManager->fetchDbTable(tofTotbCorr);

  cout<<tofTotbCorr->getVersion()<<endl;
  //cout<<tofTotbCorr->getTableName()<<endl;
  cout<<tofTotbCorr->getBeginDateTime()<<endl;
  cout<<tofTotbCorr->getEndDateTime()<<endl;


  tofTotbCorr_st* totcorr = static_cast<tofTotbCorr_st*>(tofTotbCorr->GetTable());

  if(!totcorr) {
    cout << " ahhhhhh! " << endl;
    return;
  }

  cout<<"Read out from DataBase-------------->"<<endl;

  Int_t nRows = tofTotbCorr->GetNRows();
  cout << " NRows = " << nRows << endl;
  //if(nRows!=mNTray*mNTDIG) {
  //  cout << " NRows doesn't match !!!" << endl;
  //}

  ofstream outData;
  outData.open("totCorr_read.dat");
  //for(int i=0;i<mNTray*mNTDIG;i++) {
  for(int i=0;i<nRows;i++) {
    outData << setw(6) << totcorr[i].trayId << setw(6) << totcorr[i].moduleId << setw(6) << totcorr[i].cellId << setw(6) << totcorr[i].tdcId << endl;
    for(int j=0;j<60;j++) {
      if(totcorr[i].tot[j]>1.e-4)
       outData << setw(15) << totcorr[i].tot[j];
    }
    outData << endl;
    for(int j=0;j<60;j++) {
      if(totcorr[i].tot[j]>1.e-4)
        outData << setw(15) << totcorr[i].corr[j];
    }
    outData << endl;
  }
  outData.close();

}