// --------------------------------------------------------------------------
//
// Macro to create a geometry file for the TRD
// This file is needed to create the digi parameters as
// next step
//
// F. Uhlig 07.05.2010
//
// D. Emschermann - 26.05.2010 - fix TRD geo input file
//
// expecting the trd_v13q.root file in directory "geometry"
//
// --------------------------------------------------------------------------

void create_geometry_file_root_test(TString geoName = "trd_v15a_1e")
{
   TString inDir   = gSystem->Getenv("VMCWORKDIR");

   TString outDir  = "data/";
   TString outFile = outDir + "test.mc."+ geoName + ".root";
   TString parFile = outDir + "params.root";

   TString caveGeom   = "cave.geo";
   TString trdGeom    = inDir + "/geometry/trd/" + geoName + ".geo.root";
   TString trdGeomOut = outDir + "geofile_" + geoName + ".root";

   // Function needed for CTest runtime dependency
   TString depFile = Remove_CTest_Dependency_File(outDir, "create_geometry_file_root_test");

   FairRunSim* run = new FairRunSim();
   run->SetName("TGeant3");
   run->SetOutputFile(outFile);
   run->SetMaterials("media.geo");

   if (caveGeom != "") {
      FairModule* cave = new CbmCave("CAVE");
      cave->SetGeometryFileName(caveGeom);
      run->AddModule(cave);
   }
   if (trdGeom != "") {
      FairDetector* trd = new CbmTrd("TRD", kTRUE);
      trd->SetGeometryFileName(trdGeom);
      run->AddModule(trd);
   }

   FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
   run->SetGenerator(primGen);
   FairBoxGenerator* boxGen1 = new FairBoxGenerator(11, 20);
   boxGen1->SetPRange(2.,2.);
   boxGen1->SetPhiRange(0.,0.);
   boxGen1->SetThetaRange(0.,0.);
   boxGen1->SetCosTheta();
   boxGen1->Init();
   primGen->AddGenerator(boxGen1);

   // Constant Field
   CbmFieldConst* magField = new CbmFieldConst();
   magField->SetField(0, 14 ,0 ); // values are in kG
   magField->SetFieldRegion(-74, -39 ,-22 , 74, 39 , 160);
   run->SetField(magField);

   run->Init();
   run->Run(1);

   // Create separate file with the geomanager only
   run->CreateGeometryFile(trdGeomOut);

   // ------------------------------------------------------------------------

   cout << " Test passed" << endl;
   cout << " All ok " << endl;

  // Function needed for CTest runtime dependency
  Generate_CTest_Dependency_File(depFile);

}