readtofSimResParams.C

// macro to read tofSimResTables from database
//
// Frank Geurts Nov. 24, 2020

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

void readtofSimResParams(string ZReadTime = "2029-12-31 23:59:59")
{
  const int mNTray = 120;
  const int mNCellsPerTray = 192;
  const int mNVPD = 19;
  //double params[mNTray][mNCellsPerTray];

  //-- 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());

// --- BTOF
  StDbTable* tofSimResParams = configNode->addDbTable("tofSimResParams");
  dbManager->fetchDbTable(tofSimResParams);

  cout<<"--- BTOF table ---"<< endl;
  cout<<"version : " << tofSimResParams->getVersion()<<endl;
  cout<<"begin data/time :" << tofSimResParams->getBeginDateTime()<<endl;
  cout<<"end data/time   :" << tofSimResParams->getEndDateTime()<<endl;

  tofSimResParams_st* tofTable = static_cast<tofSimResParams_st*>(tofSimResParams->GetTable());
  if(!tofTable) {
    cout << " No BTOF table? arrghhhhhh! " << endl;
    return;
  }

// --- VPD
  StDbTable* vpdSimParams = configNode->addDbTable("vpdSimParams");
  dbManager->fetchDbTable(vpdSimParams);

  cout<<"--- VPD table ---"<< endl;
  cout<<"version : " << vpdSimParams->getVersion()<<endl;
  cout<<"begin data/time :" << vpdSimParams->getBeginDateTime()<<endl;
  cout<<"end data/time   :" << vpdSimParams->getEndDateTime()<<endl;

  vpdSimParams_st* vpdTable = static_cast<vpdSimParams_st*>(vpdSimParams->GetTable());
  if(!vpdTable) {
    cout << " No VPD table? arrghhhhhh! " << endl;
    return;
  }


  //--- store BTOF data in text file
  cout <<"--- Read from BTOF table ---"<<endl;
  Int_t nRows = tofSimResParams->GetNRows();
  cout << " NRows = " << nRows << endl;

  ofstream outData;
  outData.open("tofSimResParams.dat");

  for ( int row = 0; row < nRows; row++){
    //cout << "row " << setw(6) << row << ": ";
    for ( int i = 0; i < mNTray; i++ ){ //  nTrays
      //cout << ".";
             for ( int j = 0; j < mNCellsPerTray; j++ ){
         size_t index = i * mNTray + j;
                //params[i][j] =  tofTable[0].resolution[index];
                 outData  << setw(6) << i << setw(6) << j << setw(6) << tofTable[0].resolution[index]  << endl;
       }
     }
     //cout << endl;
   }
   outData.close();


// --- Store VPD data in text file
  cout <<"--- Read from VPD table ---"<<endl;
  nRows = vpdSimParams->GetNRows();
  cout << " NRows = " << nRows << endl;

  outData.open("vpdSimParams.dat");

  for ( int row = 0; row < nRows; row++){
    //cout << "row " << setw(6) << row << ": ";
    for ( int i = 0; i < 2*mNVPD; i++ ){ //  nTubes
      //params.tubeId = table->tubeID[i];
      //params.singleTubeRes = table->tubeRes[i];
      //params.tubeStatusFlag = table->tubeStatusFlag[i];
      //params.tubeTriggerFlag = table->tubeTriggerFlag[i];
      //mSimParams[table->tubeID[i]] = params;
      outData << setw(6) << i << setw(6) << vpdTable->tubeRes[i] << endl;

     }
     //cout << endl;
  }
    outData.close();

}