gphysics.cc

Go to the documentation of this file.

// gphysics
#include "gphysics.h"
#include "gphysics_options.h"
#include "gphysicsConventions.h"
#include "gphysicsRuntime.h"
 
// gemc
#include "gutilities.h"
 
// geant4 version
#include "G4Version.hh"
 
// geant4
 
// From examples/extended/physicslists/extensibleFactory:
// The extensible factory also allows for the extension of lists by adding
// (using "+" as a separator) or replacing (using "_" as a separator)
// specific physics constructors.  These can be specified by
// pre-defining a short name (e.g. RADIO for G4RadioactiveDecayPhysics) or
// providing the full class name.
#include "G4PhysListFactoryAlt.hh"
 
// allow ourselves to give the user extra info about available physics constructors (ctors)
#include "G4PhysicsConstructorFactory.hh"
 
 
// Implementation note (non-Doxygen):
// The public API and semantics are documented in gphysics.h. This translation unit focuses on:
// - Reading options
// - Instantiating the requested reference physics list using the Geant4 extensible factory
// - Registering default add-on constructors
GPhysics::GPhysics(const std::shared_ptr<GOptions>& gopts) :
    GBase(gopts, GPHYSICS_LOGGER),
    physList(nullptr) {
    log->debug(CONSTRUCTOR, "GPhysics");
 
    bool        showPhys  = gopts->getSwitch("showPhysics");
    std::string gphysList = gopts->getRequiredScalarString("phys_list");
 
    if (showPhys) {
        // Diagnostic mode: print what Geant4 advertises and return without building a list.
        printAvailable();
        return;
    }
 
    // g4alt::G4PhysListFactoryAlt is the extensible factory
    // including the G4PhysListFactoryAlt.hh header and the line:
    //    using namespace g4alt;
    // would make this a drop-in replacement, but we'll list the explicit
    // namespace here just for clarity
    // The Geant4 factory syntax is whitespace-sensitive; normalize user input.
    std::string              g4physList = gutilities::removeAllSpacesFromString(gphysList);
 
    try {
        physList = new gphysics::RuntimePhysicsList(g4physList);
    } catch (const std::exception& error) {
        log->error(gphysics::ERR_PHYSLISTERROR, "%1", error.what());
    }
    if (const auto cut = gopts->getOptionalScalarDouble("production_cut")) {
        std::string error;
        if (!static_cast<gphysics::RuntimePhysicsList*>(physList)->setProductionCut(*cut, error)) {
            log->error(gphysics::ERR_PHYSLISTERROR, "%1", error);
        }
    }
 
    if (!physList) {
        log->error(gphysics::ERR_PHYSLISTERROR, "physics list <%1> could not be loaded.", gphysList);
    }
 
    log->info(2, "G4PhysListFactory: <%1> loaded.", g4physList);
}
 
 
// Implementation note (non-Doxygen):
// Prints Geant4 version metadata and queries Geant4 registries for the available lists/constructors.
void GPhysics::printAvailable() const {
    std::string g4ver = gutilities::replaceCharInStringWithChars(G4Version, "$", "");
 
    log->info(0, "Geant4 Version %1 %2", g4ver, G4Date);;
 
    g4alt::G4PhysListFactory factory;
    factory.PrintAvailablePhysLists();
 
    log->info(0, "Available Geant4 Physics Lists:");
 
    G4PhysicsConstructorRegistry* g4pctorFactory = G4PhysicsConstructorRegistry::Instance();
    g4pctorFactory->PrintAvailablePhysicsConstructors();
}