// --------------------------------------------------------------------------
//
// Macro to create digi parameters for the TRD
// This macro needs a geometry file with the TRD geometry only.
//
// The way how the pad layout looks like has to be implemented in a task called
// in this macro. CbmTrdCreateDigiPar is the implementation used here.
//
// F.Uhlig 07.05.2010
//
// --------------------------------------------------------------------------

void create_digipar_root(TString geoName = "trd_v13q")
{
   TString inFile   = "data/test.mc." + geoName + ".root";
   TString geoFile  = "geofile_" + geoName + ".root";
   TString outFile  = "data/test.esd." + geoName + ".root";
   TString digiFile = geoName + ".digi.par.long";  // Digi Parameter Output File

   FairRunAna* run = new FairRunAna();
   run->SetInputFile(inFile);
   run->SetOutputFile(outFile);
   run->SetGeomFile(geoFile);

   FairRuntimeDb* rtdb = run->GetRuntimeDb();

   FairParAsciiFileIo* parIo2 = new FairParAsciiFileIo();
   parIo2->open(digiFile, "out");
   rtdb->setOutput(parIo2);

   CbmTrdCreateDigiPar* trdDigitizer = new CbmTrdCreateDigiPar();
   run->AddTask(trdDigitizer);

   rtdb->saveOutput();
   run->Init();
   rtdb->print();
   CbmTrdDigiPar* digiPar = (CbmTrdDigiPar*) rtdb->getContainer("CbmTrdDigiPar");
   digiPar->setChanged();
   digiPar->setInputVersion(run->GetRunId(), 1);
   rtdb->print();
   rtdb->saveOutput();
}