RunManager.cc

Go to the documentation of this file.

// $Header: /cvs/lcd/slic/src/RunManager.cc,v 1.17 2010/01/08 23:18:56 jeremy Exp $
#include "RunManager.hh"

// lcdd
#include "LCDDParser.hh"
#include "LCDDProcessor.hh"

// slic
#include "EventAction.hh"
#include "EventSourceManager.hh"
#include "ModuleRegistry.hh"
#include "PhysicsListManager.hh"
#include "PrimaryGeneratorAction.hh"
#include "RunAction.hh"
#include "SteppingAction.hh"
#include "TrackingAction.hh"

// geant4
#include "G4StateManager.hh"
#include "G4UImanager.hh"

namespace slic
{

    RunManager::RunManager() :
        Module("RunManager", false),
        m_userActionsInitialized(false)
    {}

    RunManager::~RunManager()
    {}

    void RunManager::initializeUserActions()
    {
        SetUserAction(new PrimaryGeneratorAction);
        SetUserAction(new RunAction);
        SetUserAction(new EventAction);
        SetUserAction(new TrackingAction);
        SetUserAction(new SteppingAction);
        m_userActionsInitialized = true;
    }

    void RunManager::Initialize()
    {
        // This makes sure that physics initialization occurs before other user actions.
        G4RunManager::Initialize();

        // Initialize user actions here to avoid ordering problems.
        if (!m_userActionsInitialized)
            initializeUserActions();

        // Initialize the event generation manager.
        EventSourceManager::instance();

        // Setup the default event source.
        EventSourceManager::instance()->setupEventSource();

        // Print list of registered modules.
#ifdef SLIC_LOG
        ModuleRegistry::instance()->print();
#endif
    }

    void RunManager::InitializePhysics()
    {
        // Initialize the physics list.
        PhysicsListManager::instance()->initializePhysicsList();

        // Call the G4RunManager's intitialization method.
        G4RunManager::InitializePhysics();

        // Check if the LCDD subsystem got some limits.
        LCDDProcessor* lcdd = LCDDProcessor::instance();
        PhysicsListManager* pmgr = PhysicsListManager::instance();
        if ( lcdd->getLimitSetsBegin() != lcdd->getLimitSetsEnd() ) 
        {
            pmgr->enableLimits(true);            
        }

        // Enable physics limits, if necessary.
        if ( pmgr->enableLimits() )  
        {
            pmgr->setupUserLimitsProcesses();
        }
    }

    void RunManager::InitializeGeometry()
    {
        if ( !LCDDParser::instance()->isValidSetup() ) 
        {
            G4Exception("FATAL ERROR: Current LCDD geometry setup is not valid.");
        }

        G4RunManager::InitializeGeometry();
    }

    void RunManager::BeamOn(G4int n_event, const char* macroFile, G4int n_select)
    {
        // Attempt to initialize Geant4 if not in idle state.
        if ( G4StateManager::GetStateManager()->GetCurrentState() != G4State_Idle ) 
        {
            this->Initialize();
        }

        m_numberOfEventsToRun = n_event;

        G4RunManager::BeamOn(n_event, macroFile, n_select);
    }

    int RunManager::getNumberOfEventsToRun()
    {
        return m_numberOfEventsToRun;
    }
}

---