/********************************************************************************
* Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH *
* *
* This software is distributed under the terms of the *
* GNU Lesser General Public Licence version 3 (LGPL) version 3, *
* copied verbatim in the file "LICENSE" *
********************************************************************************/
Int_t sql_params_read()
{
// Create a Runtime Database singleton.
FairRuntimeDb* db = FairRuntimeDb::instance();
// Set the SQL IO as first input
FairParTSQLIo* inp = new FairParTSQLIo();
// Verbosity level
inp->SetVerbosity(1);
inp->open();
db->setFirstInput(inp);
// Create the container via the factory if not already created
FairDbTutPar* p1 = (FairDbTutPar*)(db->getContainer("TUTParDefault"));
FairDbTutPar* p2 = (FairDbTutPar*)(db->getContainer("TUTParAlternative"));
// Create a dummy runID using date in UTC from which
// corresponding parameters will be initialised
ValTimeStamp tStamp(2015,02,20,10,10,10);
UInt_t runId = tStamp.GetSec();
cout << "-I- looking for parameters at runID# " << runId << endl;
cout << "-I- corresponding time in runID (UTC) " << tStamp.Format("iso") << endl;
// Use the generated RunID to initialised the parameter
// using the SQL-based IO input
db->initContainers(runId);
// Get the container after initialisation
// from the RuntimeDB
FairDbTutPar* pp1 = (FairDbTutPar*)(db->getContainer("TUTParDefault"));
FairDbTutPar* pp2 = (FairDbTutPar*)(db->getContainer("TUTParAlternative"));
cout << endl;
cout << "-I- Reading Parameter data from SQL Database: \n" << endl;
cout << endl;
pp1->Print();
pp2->Print();
cout << endl;
if (db) delete db;
return 0;
}
Int_t calgeo_read()
{
// Create a Runtime Database singleton.
FairRuntimeDb* db = FairRuntimeDb::instance();
// Set the SQL IO as first input
FairParTSQLIo* inp = new FairParTSQLIo();
// Verbosity level
//inp->SetVerbosity(1);
inp->open();
db->setFirstInput(inp);
// Create in memory the relevant container
R3BCaloGeometryPar* par = (R3BCaloGeometryPar*)(db->getContainer("CaloGeometryPar"));
// Create a dummy runID using date in UTC from which
// corresponding parameters will be initialised
ValTimeStamp tStamp(2014,03,21,18,00,00);
UInt_t runId = tStamp.GetSec();
cout << "-I- looking for parameters at runID# " << runId << endl;
cout << "-I- corresponding time in runID (UTC) " << tStamp.Format("iso") << endl;
// Use the generated RunID to initialised the parameter
// using the SQL-based IO input
db->initContainers(runId);
// Get the container after initialisation
// from the RuntimeDB
R3BCaloGeometryPar* par = (R3BCaloGeometryPar*)(db->getContainer("CaloGeometryPar"));
cout << endl;
cout << "-I- Reading Parameter data from SQL Database: \n" << endl;
cout << endl;
par->Print();
cout << endl;
if (db) delete db;
return 0;
}
/********************************************************************************
* Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH *
* *
* This software is distributed under the terms of the *
* GNU Lesser General Public Licence version 3 (LGPL) version 3, *
* copied verbatim in the file "LICENSE" *
********************************************************************************/
Int_t sql_params_read_bin()
{
// ---- Load libraries -------------------------------------------------
gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
basiclibs();
gSystem->Load("libGenVector");
gSystem->Load("libGeoBase");
gSystem->Load("libFairDB");
gSystem->Load("libParBase");
gSystem->Load("libBase");
gSystem->Load("libMCStack");
gSystem->Load("libField");
gSystem->Load("libTutorial5");
// Create a Runtime Database singleton.
FairRuntimeDb* db = FairRuntimeDb::instance();
// Set the SQL IO as first input
FairParTSQLIo* inp = new FairParTSQLIo();
// Verbosity level
inp->SetVerbosity(1);
inp->open();
db->setFirstInput(inp);
// Create the container via the factory if not already created
FairDbTutParBin* p1 = (FairDbTutParBin*)(db->getContainer("TUTParBin"));
p1->Print();
// Create a dummy runID using date in UTC from which
// corresponding parameters will be initialised
ValTimeStamp tStamp(2015,02,20,10,07,48);
UInt_t runId = tStamp.GetSec();
cout << "-I- looking for parameters at runID# " << runId << endl;
cout << "-I- corresponding time in runID (UTC) " << tStamp.AsString("c") << endl;
// Use the generated RunID to initialised the parameter
// using the SQL-based IO input
db->initContainers(runId);
cout << endl;
cout << "-I- Initialisation from SQL Database:" << endl;
cout << endl;
// Get the container after initialisation
// from the RuntimeDB
FairDbTutParBin* pp1 = (FairDbTutParBin*)(db->getContainer("TUTParBin"));
pp1->Print();
cout << endl;
if (db) delete db;
return 0;
}
void run(TString runNumber)
{
TStopwatch timer;
timer.Start();
const Int_t nModules = 200;
TString dirIn1 = "/Volumes/Data/kresan/s438/data/";
TString dirIn2 = "/Volumes/Data/kresan/s438/tcal/";
TString dirOut = "/Volumes/Data/kresan/s438/tcal/";
TString inputFileName1 = dirIn1 + runNumber + "_tcal.root"; // name of input file
TString inputFileName2 = dirIn2 + runNumber + "_tcal_1.root"; // name of input file
TString parFileName = dirIn1 + "params_" + runNumber + "_tcal.root"; // name of parameter file
TString outputFileName = dirOut + runNumber + "_tcal_temp.root"; // name of output file
// Create analysis run -------------------------------------------------------
Int_t nFiles;
if (0 == runNumber.CompareTo("r122"))
{
nFiles = 22;
}
else if (0 == runNumber.CompareTo("r126"))
{
nFiles = 2;
}
// Create analysis run -------------------------------------------------------
FairRunAna* run = new FairRunAna();
run->SetInputFile(inputFileName1.Data());
run->AddFriend(inputFileName2.Data());
for (Int_t i = 1; i < nFiles; i++)
{
inputFileName1 = dirIn1 + runNumber + "_tcal_";
inputFileName1 += i;
inputFileName1 += ".root";
run->AddFile(inputFileName1.Data());
}
run->SetOutputFile(outputFileName.Data());
// ---------------------------------------------------------------------------
// ----- Runtime DataBase info -----------------------------------------------
FairRuntimeDb* rtdb = run->GetRuntimeDb();
FairParRootFileIo* parIo1 = new FairParRootFileIo();
parIo1->open(parFileName);
rtdb->setFirstInput(parIo1);
rtdb->setOutput(parIo1);
// Set the SQL IO as second input
FairParTSQLIo* inp = new FairParTSQLIo();
inp->SetVerbosity(1);
inp->open();
rtdb->setFirstInput(inp);
rtdb->saveOutput();
// ---------------------------------------------------------------------------
// Time calibration ----------------------------------------------------------
R3BLandTcal* landTcal = new R3BLandTcal("LandTcal", 1);
landTcal->SetNofModules(nModules, 40);
run->AddTask(landTcal);
R3BLosTcal* losTcal = new R3BLosTcal("LosTcal", 1);
losTcal->SetNofModules(20);
run->AddTask(losTcal);
// ---------------------------------------------------------------------------
// Initialize ----------------------------------------------------------------
run->Init();
FairLogger::GetLogger()->SetLogScreenLevel("INFO");
// ---------------------------------------------------------------------------
// Run -----------------------------------------------------------------------
run->Run(0, 100000);
// ---------------------------------------------------------------------------
timer.Stop();
Double_t rtime = timer.RealTime();
Double_t ctime = timer.CpuTime();
cout << endl << endl;
cout << "Macro finished succesfully." << endl;
cout << "Output file is " << outputFileName << endl;
cout << "Parameter file is " << parFileName << endl;
cout << "Real time " << rtime << " s, CPU time " << ctime << "s" << endl << endl;
}
