diff --git a/macro/mcbm/mcbm_hadron_analysis.C b/macro/mcbm/mcbm_hadron_analysis.C index 3dfdce70122e82fa6e30eb7541222f7a3385975e..f2a76053eb99ca79e876abee305bea7cdb5c2ef3 100644 --- a/macro/mcbm/mcbm_hadron_analysis.C +++ b/macro/mcbm/mcbm_hadron_analysis.C @@ -37,26 +37,32 @@ void mcbm_hadron_analysis(Int_t nEvents = 10, TString setupName = "mcbm_beam_202 parInput1->open(ParFile.Data()); rtdb->setFirstInput(parInput1); - // ------ Geometry setup ------------------------------------------------ - std::cout << '\n'; - TString setupFile = srcDir + "/geometry/setup/setup_" + setupName + ".C"; - TString setupFunct = "setup_"; - setupFunct = setupFunct + setupName + "()"; - std::cout << "-I- Loading macro " << setupFile << '\n'; - gROOT->LoadMacro(setupFile); - gROOT->ProcessLine(setupFunct); + // ----- Load the geometry setup ------------------------------------- + std::cout << std::endl; + // std::cout << "-I- " << myName << ": Loading setup " << setup << std::endl; CbmSetup* setup = CbmSetup::Instance(); - // setup->RemoveModule(ECbmModuleId::kTrd); + setup->LoadSetup(setupName); + TString geoTag; // ------------------------------------------------------------------------ + // ----- MuCh digi parameters initialization -------------------------------------- if (setup->IsActive(ECbmModuleId::kMuch)) { // Parameter file name - TString geoTag; setup->GetGeoTag(ECbmModuleId::kMuch, geoTag); Int_t muchFlag = (geoTag.Contains("mcbm") ? 1 : 0); TString parFile = gSystem->Getenv("VMCWORKDIR"); - parFile += "/parameters/much/much_" + geoTag(0, 4) + "_digi_sector.root"; + + if (muchFlag) { + std::cout << geoTag << std::endl; + parFile = parFile + "/parameters/much/much_" + geoTag + "_digi_sector.root"; + std::cout << "Using parameter file " << parFile << std::endl; + } + else { + std::cout << geoTag(0, 4) << std::endl; + parFile = parFile + "/parameters/much/much_" + geoTag(0, 4) + "_digi_sector.root"; + std::cout << "Using parameter file " << parFile << std::endl; + } // Initialization of the geometry scheme auto muchGeoScheme = CbmMuchGeoScheme::Instance();