storetofTDIGOnTray.C

// macro to upload tofr5 INL tables to database
//
// based on
//  http://www.star.bnl.gov/STAR/comp/db/StoreDbTable.cc.html
//
// Jing Liu, 02/18/2005 
//

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

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


void storetofTDIGOnTray()
//int main(int argc, char *argv[])

{

  const Int_t NBOARD = 8;
// year8
//  const Int_t NTRAY = 5;
// year9
  const Int_t NTRAY = 94;
  const Int_t NVPDTRAY = 2;
  const Int_t NMAX = 120;

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

  //----------------------------------------
//  TString ZStoreTime = "2007-11-01 00:00:03";
  TString ZStoreTime = "2008-11-01 00:00:00";

  //-- add table to the container with descriptor given by Database
  StDbTable* tofTDIGOnTray = configNode->addDbTable("tofTDIGOnTray");

  //-- fill structures & store times
  tofTDIGOnTray_st *tdig = new tofTDIGOnTray_st[NMAX+NVPDTRAY];

// year8
/*
  Short_t trayId[NTRAY+NVPDTRAY] = {76, 77, 78, 79, 80, 121, 122};
  Short_t tdigId[NTRAY+NVPDTRAY][NBOARD] = {
    36, 26, 27, 29, 25, 23, 22, 21,
    30, 31, 33, 35, 28, 32, 34, 37,
    74, 63, 49, 72, 64, 53, 59, 52,
    70, 51, 75, 65, 66, 50, 71, 68,
    46, 62, 72, 69, 61, 58, 42, 43,
    89, 91,  0,  0, 84, 82,  0,  0,
    85, 79,  0,  0, 80, 83,  0,  0
  };
*/
// year9
  Short_t validtrayId[NTRAY+NVPDTRAY] =
                          { 1,  2,  3,  4,  5,  6,  7,  8,  9, 10,
                           11, 12,         15, 16, 17, 18, 19, 20,
                           21, 22, 23, 24, 25, 26, 27, 28, 29, 30,
                           31, 32, 33, 34, 35, 36, 37, 38, 39, 40,
                           41,         44, 45, 46, 47, 48, 49, 50,
                           51, 52, 53, 54, 55, 56, 57, 58, 59, 60,
                           61, 62, 63, 64, 65, 66, 67, 68, 69, 70,
                           71, 72,         75, 76, 77, 78,

                                                   97, 98, 99, 100,
                          101,           104, 105, 106, 107, 108, 109, 110,
                          111, 112, 113, 114, 115, 116, 117, 118, 119, 120,
                          121, 122};
  Short_t trayId[NMAX] = { 68, 58, 63,  6, 71, 67, 61, 50, 57, 60,
                           47, 46,  0,  0, 64, 72, 62, 42, 14, 11,
                           28,  7, 25, 33, 19, 17, 37, 21, 30, 35,
                            8, 20, 31, 27, 12,  9, 10, 29, 34, 13,
                           16,  0,  0, 36, 39, 23, 32, 15, 24, 18,
                           59, 66, 70, 38, 41, 45, 52, 54, 53, 55,
                           84, 75, 94, 91, 95, 48, 65, 22, 40, 49,
                            5,  2,  0,  0,  1,  4,  3, 51,  0,  0,
                            0,  0,  0,  0,  0,  0,  0,  0,  0,  0,
                            0,  0,  0,  0,  0,  0, 77, 89, 87, 90,
                           93,  0,  0, 73, 76, 83, 69, 80, 82, 74,
                          123, 79, 43, 86, 85,122, 81, 92, 96, 78};
                              // B-122;  C-123  renamed

  Short_t tdigId[NMAX][NBOARD];
  Short_t tdigId_Vpd[NVPDTRAY][NBOARD] = { 89, 91, 1102, 0, 885, 82, 0, 0,
                                           85, 79, 1101, 0,  80, 83, 0, 0};

  ifstream infile("TDIGMap");
  Short_t tId[NTRAY], boardId[NTRAY][NBOARD];
  for(int i=0;i<NTRAY;i++) {
    for(int j=0;j<NBOARD;j++) {
      int trayNum, iboard, boardNum;
      infile >> trayNum >> iboard >> boardNum;
      cout << trayNum << " " << boardNum << endl;
      tId[i] = (Short_t)trayNum;
      boardId[i][iboard-1] = (Short_t)boardNum;
    }
  }
  infile.close();

  for(int i=0;i<NMAX;i++) {
    for(int j=0;j<NBOARD;j++) {
      tdigId[i][j] = 0;
    }

    Short_t thisTray = trayId[i];
    int index = -1;
    for(int j=0;j<NTRAY;j++) {
      if(thisTray==tId[j]) {
        index = j;
        break;
      }
    }

    cout << "tray = " << thisTray << " index = " << index << endl;
    if(index<0||index>=NTRAY) continue;
    for(int j=0;j<NBOARD;j++) {
      tdigId[i][j] = boardId[index][j];
    }
  }

/*
  {
    89, 91,  0,  0, 84, 82,  0,  0,
    85, 79,  0,  0, 80, 83,  0,  0
  };
*/

  for(int i=0;i<NMAX+NVPDTRAY;i++) {
    tdig[i].trayId = i+1;
    for(int j=0;j<NBOARD;j++) {
      if(i<NMAX) {
        tdig[i].tdigId[j] = tdigId[i][j];
      } else {
        tdig[i].tdigId[j] = tdigId_Vpd[i-NMAX][j];
      }
    }
  }

  ofstream outData("testTDIGMap.dat");
  for(int i=0;i<NMAX+NVPDTRAY;i++) {
    for(int j=0;j<NBOARD;j++) {
      outData << setw(5) << tdig[i].tdigId[j];
    }
    outData << endl;
  }
  outData.close();

  //- store data in table
  tofTDIGOnTray->SetTable((char*)tdig, NMAX+NVPDTRAY);
  //- set store time
  dbManager->setStoreTime(ZStoreTime.Data());
  //- store table in dBase
  cout<<" here "<<endl;
  dbManager->storeDbTable(tofTDIGOnTray);
  cout<<"uploaded"<<endl;
      //    return 0;
}