G4AdjointSimManager.cc

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
// G4AdjointCrossSurfChecker implementation
//
// --------------------------------------------------------------------
//   Class Name:   G4AdjointCrossSurfChecker
//   Author:       L. Desorgher, 2007-2009
//   Organisation: SpaceIT GmbH
//   Contract:     ESA contract 21435/08/NL/AT
//   Customer:     ESA/ESTEC
// --------------------------------------------------------------------
 
#include "G4AdjointSimManager.hh"
 
#include "G4AdjointCrossSurfChecker.hh"
#include "G4AdjointPrimaryGeneratorAction.hh"
#include "G4AdjointSimMessenger.hh"
#include "G4AdjointStackingAction.hh"
#include "G4AdjointSteppingAction.hh"
#include "G4AdjointTrackingAction.hh"
#include "G4ParticleTable.hh"
#include "G4PhysicsLogVector.hh"
#include "G4Run.hh"
#include "G4RunManager.hh"
#include "G4WorkerRunManager.hh"
#include "G4UserEventAction.hh"
#include "G4UserRunAction.hh"
#include "G4UserStackingAction.hh"
#include "G4UserSteppingAction.hh"
#include "G4UserTrackingAction.hh"
#include "G4VUserPrimaryGeneratorAction.hh"
 
// --------------------------------------------------------------------
//
G4ThreadLocal G4AdjointSimManager* G4AdjointSimManager::instance = nullptr;
 
// --------------------------------------------------------------------
//
G4AdjointSimManager::G4AdjointSimManager()
{
  theAdjointPrimaryGeneratorAction = new G4AdjointPrimaryGeneratorAction();
  theAdjointSteppingAction = new G4AdjointSteppingAction();
  theAdjointTrackingAction = new G4AdjointTrackingAction(theAdjointSteppingAction);
  theAdjointStackingAction = new G4AdjointStackingAction(theAdjointTrackingAction);
  theAdjointTrackingAction->SetListOfPrimaryFwdParticles(
  theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles());
 
  theMessenger = new G4AdjointSimMessenger(this);
}
 
// --------------------------------------------------------------------
//
G4AdjointSimManager::~G4AdjointSimManager()
{
  delete theAdjointRunAction;
  delete theAdjointPrimaryGeneratorAction;
  delete theAdjointSteppingAction;
  delete theAdjointEventAction;
  delete theAdjointTrackingAction;
  delete theAdjointStackingAction;
  delete theMessenger;
}
 
// --------------------------------------------------------------------
//
G4AdjointSimManager* G4AdjointSimManager::GetInstance()
{
  if (instance == nullptr) instance = new G4AdjointSimManager;
  return instance;
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::RunAdjointSimulation(G4int nb_evt)
{
 
  if (welcome_message) {
    G4cout << "****************************************************************" << std::endl;
    G4cout << "*** Geant4 Reverse/Adjoint Monte Carlo mode                  ***" << std::endl;
    G4cout << "*** Author:       L.Desorgher                                ***" << std::endl;
    G4cout << "*** Company:      SpaceIT GmbH, Bern, Switzerland            ***" << std::endl;
    G4cout << "*** Sponsored by: ESA/ESTEC contract contract 21435/08/NL/AT ***" << std::endl;
    G4cout << "****************************************************************" << std::endl;
    welcome_message = false;
  }
 
  // Switch to adjoint simulation mode
  //---------------------------------------------------------
  SwitchToAdjointSimulationMode();
 
  // Make the run
  //------------
  nb_evt_of_last_run = nb_evt;
  if (G4Threading::G4GetThreadId() < 0){
    G4RunManager::GetRunManager()->BeamOn(
    G4int(nb_evt * theAdjointPrimaryGeneratorAction->GetNbOfAdjointPrimaryTypes()));
  }
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetRestOfAdjointActions()
{ if (G4Threading::G4GetThreadId() == -1) return;
  G4RunManager* theRunManager = G4RunManager::GetRunManager();
 
  if (!user_action_already_defined) DefineUserActions();
 
  // Replace the user action by the adjoint actions
  //-------------------------------------------------
 
  theRunManager->G4RunManager::SetUserAction(theAdjointEventAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointSteppingAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointTrackingAction);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SwitchToAdjointSimulationMode()
{
  // Replace the user defined actions by the adjoint actions
  //---------------------------------------------------------
  SetAdjointActions();
 
  // Update the list of primaries
  //-----------------------------
  if (G4Threading::G4GetThreadId() != -1)  theAdjointPrimaryGeneratorAction->UpdateListOfPrimaryParticles();
  adjoint_sim_mode = true;
  ID_of_last_particle_that_reach_the_ext_source = 0;
 
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::BackToFwdSimulationMode()
{
  // Restore the user defined actions
  //--------------------------------
  ResetUserActions();
  adjoint_sim_mode = false;
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointActions()
{
  auto theRunManager = G4RunManager::GetRunManager();
 
 
  if (!user_action_already_defined) DefineUserActions();
 
 
  // Replace the user action by the adjoint actions
  //-------------------------------------------------
  theRunManager->G4RunManager::SetUserAction(this);
  if (G4RunManager::GetRunManager()->GetRunManagerType() == G4RunManager::sequentialRM ||
          G4RunManager::GetRunManager()->GetRunManagerType() == G4RunManager::workerRM) {
  theRunManager->G4RunManager::SetUserAction(theAdjointPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointStackingAction);
  if (use_user_StackingAction)
    theAdjointStackingAction->SetUserFwdStackingAction(fUserStackingAction);
  else
    theAdjointStackingAction->SetUserFwdStackingAction(nullptr);
  theRunManager->G4RunManager::SetUserAction(theAdjointEventAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointSteppingAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointTrackingAction);
  if (use_user_TrackingAction)
    theAdjointTrackingAction->SetUserForwardTrackingAction(fUserTrackingAction);
  else
    theAdjointTrackingAction->SetUserForwardTrackingAction(nullptr);
}
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointPrimaryRunAndStackingActions()
{
  G4RunManager* theRunManager = G4RunManager::GetRunManager();
 
  if (!user_action_already_defined) DefineUserActions();
 
  // Replace the user action by the adjoint actions
  //-------------------------------------------------
 
  theRunManager->G4RunManager::SetUserAction(theAdjointRunAction);
  if (G4Threading::G4GetThreadId() != -1) {
  theRunManager->G4RunManager::SetUserAction(theAdjointPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointStackingAction);
  if (use_user_StackingAction)
    theAdjointStackingAction->SetUserFwdStackingAction(fUserStackingAction);
  else
    theAdjointStackingAction->SetUserFwdStackingAction(nullptr);
}
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::ResetUserActions()
{
  G4RunManager* theRunManager = G4RunManager::GetRunManager();
 
  // Restore the user defined actions
  //-------------------------------
  theRunManager->G4RunManager::SetUserAction(fUserRunAction);
  if (G4Threading::G4GetThreadId() != -1) {
  theRunManager->G4RunManager::SetUserAction(fUserEventAction);
  theRunManager->G4RunManager::SetUserAction(fUserSteppingAction);
  theRunManager->G4RunManager::SetUserAction(fUserTrackingAction);
  theRunManager->G4RunManager::SetUserAction(fUserPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(fUserStackingAction);
  }
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::ResetRestOfUserActions()
{ 
  if (G4Threading::G4GetThreadId() == -1) return;
  G4RunManager* theRunManager = G4RunManager::GetRunManager();
 
  // Restore the user defined actions
  //-------------------------------
 
  theRunManager->G4RunManager::SetUserAction(fUserEventAction);
  theRunManager->G4RunManager::SetUserAction(fUserSteppingAction);
  theRunManager->G4RunManager::SetUserAction(fUserTrackingAction);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::ResetUserPrimaryRunAndStackingActions()
{
  G4RunManager* theRunManager = G4RunManager::GetRunManager();
  // Restore the user defined actions
  //-------------------------------
  theRunManager->G4RunManager::SetUserAction(fUserRunAction);
  if (G4Threading::G4GetThreadId() != -1) {
  theRunManager->G4RunManager::SetUserAction(fUserPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(fUserStackingAction);
  }
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::DefineUserActions()
{
  G4RunManager* theRunManager = G4RunManager::GetRunManager();
  fUserRunAction =
    const_cast<G4UserRunAction*>(theRunManager->GetUserRunAction());
 if (G4Threading::G4GetThreadId() != -1) {
  fUserTrackingAction = const_cast<G4UserTrackingAction*>(theRunManager->GetUserTrackingAction());
  fUserEventAction = const_cast<G4UserEventAction*>(theRunManager->GetUserEventAction());
  fUserSteppingAction = const_cast<G4UserSteppingAction*>(theRunManager->GetUserSteppingAction());
  theAdjointSteppingAction->SetUserForwardSteppingAction(fUserSteppingAction);
  fUserPrimaryGeneratorAction =
    const_cast<G4VUserPrimaryGeneratorAction*>(theRunManager->GetUserPrimaryGeneratorAction());
  fUserStackingAction = const_cast<G4UserStackingAction*>(theRunManager->GetUserStackingAction());
  user_action_already_defined = true;
 }
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::GetAdjointTrackingMode()
{
  return theAdjointTrackingAction->GetIsAdjointTrackingMode();
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointTrackingMode(G4bool aBool)  // could be removed
{
  adjoint_tracking_mode = aBool;
 
  if (adjoint_tracking_mode) {
    SetRestOfAdjointActions();
    theAdjointStackingAction->SetAdjointMode(true);
    theAdjointStackingAction->SetKillTracks(false);
  }
  else {
    ResetRestOfUserActions();
    theAdjointStackingAction->SetAdjointMode(false);
    if (GetDidAdjParticleReachTheExtSource()) {
      theAdjointStackingAction->SetKillTracks(false);
      RegisterAtEndOfAdjointTrack();
    }
    else
      theAdjointStackingAction->SetKillTracks(true);
  }
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::GetDidAdjParticleReachTheExtSource()
{
  return (GetNbOfAdointTracksReachingTheExternalSurface() > 0);
}
 
// --------------------------------------------------------------------
//
std::vector<G4ParticleDefinition*>* G4AdjointSimManager::GetListOfPrimaryFwdParticles()
{
  return theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles();
}
 
// --------------------------------------------------------------------
//
std::size_t G4AdjointSimManager::GetNbOfPrimaryFwdParticles()
{
  return theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles()->size();
}
 
// --------------------------------------------------------------------
//
G4ThreeVector G4AdjointSimManager::GetPositionAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetPositionAtEndOfLastAdjointTrack(i);
}
 
// --------------------------------------------------------------------
//
G4ThreeVector G4AdjointSimManager::GetDirectionAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetDirectionAtEndOfLastAdjointTrack(i);
}
 
// --------------------------------------------------------------------
//
G4double G4AdjointSimManager::GetEkinAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetEkinAtEndOfLastAdjointTrack(i);
}
 
// --------------------------------------------------------------------
//
G4double G4AdjointSimManager::GetEkinNucAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetEkinNucAtEndOfLastAdjointTrack(i);
}
 
// --------------------------------------------------------------------
//
G4double G4AdjointSimManager::GetWeightAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetWeightAtEndOfLastAdjointTrack(i);
}
 
// --------------------------------------------------------------------
//
G4double G4AdjointSimManager::GetCosthAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetCosthAtEndOfLastAdjointTrack(i);
}
 
// --------------------------------------------------------------------
//
const G4String& G4AdjointSimManager::GetFwdParticleNameAtEndOfLastAdjointTrack()
{
  return theAdjointTrackingAction->GetFwdParticleNameAtEndOfLastAdjointTrack();
}
 
// --------------------------------------------------------------------
//
G4int G4AdjointSimManager::GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack(i);
}
 
// --------------------------------------------------------------------
//
G4int G4AdjointSimManager::GetFwdParticleIndexAtEndOfLastAdjointTrack(std::size_t i)
{
  return theAdjointTrackingAction->GetLastFwdParticleIndex(i);
}
 
// --------------------------------------------------------------------
//
std::size_t G4AdjointSimManager::GetNbOfAdointTracksReachingTheExternalSurface()
{
  return theAdjointTrackingAction->GetNbOfAdointTracksReachingTheExternalSurface();
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::ClearEndOfAdjointTrackInfoVectors()
{
  theAdjointTrackingAction->ClearEndOfAdjointTrackInfoVectors();
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::RegisterAtEndOfAdjointTrack()
{
  last_pos = theAdjointSteppingAction->GetLastPosition();
  last_direction = theAdjointSteppingAction->GetLastMomentum();
  last_direction /= last_direction.mag();
  last_cos_th = last_direction.z();
  G4ParticleDefinition* aPartDef = theAdjointSteppingAction->GetLastPartDef();
 
  last_fwd_part_name = aPartDef->GetParticleName();
 
  last_fwd_part_name.erase(0, 4);
 
  last_fwd_part_PDGEncoding =
    G4ParticleTable::GetParticleTable()->FindParticle(last_fwd_part_name)->GetPDGEncoding();
 
  std::vector<G4ParticleDefinition*>* aList =
    theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles();
  last_fwd_part_index = -1;
  G4int i = 0;
  while (i < (G4int)aList->size() && last_fwd_part_index < 0) {
    if ((*aList)[i]->GetParticleName() == last_fwd_part_name) last_fwd_part_index = i;
    ++i;
  }
 
  last_ekin = theAdjointSteppingAction->GetLastEkin();
  last_ekin_nuc = last_ekin;
  if (aPartDef->GetParticleType() == "adjoint_nucleus") {
    nb_nuc = static_cast<G4double>(aPartDef->GetBaryonNumber());
    last_ekin_nuc /= nb_nuc;
  }
 
  last_weight = theAdjointSteppingAction->GetLastWeight();
 
  last_pos_vec.push_back(last_pos);
  last_direction_vec.push_back(last_direction);
  last_ekin_vec.push_back(last_ekin);
  last_ekin_nuc_vec.push_back(last_ekin_nuc);
  last_cos_th_vec.push_back(last_cos_th);
  last_weight_vec.push_back(last_weight);
  last_fwd_part_PDGEncoding_vec.push_back(last_fwd_part_PDGEncoding);
  last_fwd_part_index_vec.push_back(last_fwd_part_index);
  ID_of_last_particle_that_reach_the_ext_source++;
  ID_of_last_particle_that_reach_the_ext_source_vec.push_back(
    ID_of_last_particle_that_reach_the_ext_source);
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::DefineSphericalExtSource(G4double radius, G4ThreeVector pos)
{
  G4double area;
  return G4AdjointCrossSurfChecker::GetInstance()->AddaSphericalSurface("ExternalSource", radius,
                                                                        pos, area);
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::DefineSphericalExtSourceWithCentreAtTheCentreOfAVolume(
  G4double radius, const G4String& volume_name)
{
  G4double area;
  G4ThreeVector center;
  return G4AdjointCrossSurfChecker::GetInstance()
    ->AddaSphericalSurfaceWithCenterAtTheCenterOfAVolume("ExternalSource", radius, volume_name,
                                                         center, area);
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::DefineExtSourceOnTheExtSurfaceOfAVolume(const G4String& volume_name)
{
  G4double area;
  return G4AdjointCrossSurfChecker::GetInstance()->AddanExtSurfaceOfAvolume("ExternalSource",
                                                                            volume_name, area);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetExtSourceEmax(G4double Emax)
{
  theAdjointSteppingAction->SetExtSourceEMax(Emax);
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::DefineSphericalAdjointSource(G4double radius, G4ThreeVector pos)
{
  G4double area;
  G4bool aBool = G4AdjointCrossSurfChecker::GetInstance()->AddaSphericalSurface("AdjointSource",
                                                                                radius, pos, area);
  theAdjointPrimaryGeneratorAction->SetSphericalAdjointPrimarySource(radius, pos);
  area_of_the_adjoint_source = area;
  return aBool;
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::DefineSphericalAdjointSourceWithCentreAtTheCentreOfAVolume(
  G4double radius, const G4String& volume_name)
{
  G4double area;
  G4ThreeVector center;
  G4bool aBool =
    G4AdjointCrossSurfChecker::GetInstance()->AddaSphericalSurfaceWithCenterAtTheCenterOfAVolume(
      "AdjointSource", radius, volume_name, center, area);
  theAdjointPrimaryGeneratorAction->SetSphericalAdjointPrimarySource(radius, center);
  area_of_the_adjoint_source = area;
  return aBool;
}
 
// --------------------------------------------------------------------
//
G4bool G4AdjointSimManager::DefineAdjointSourceOnTheExtSurfaceOfAVolume(const G4String& volume_name)
{
  G4double area;
  G4bool aBool = G4AdjointCrossSurfChecker::GetInstance()->AddanExtSurfaceOfAvolume(
    "AdjointSource", volume_name, area);
  area_of_the_adjoint_source = area;
  if (aBool) {
    theAdjointPrimaryGeneratorAction->SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume(volume_name);
  }
  return aBool;
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointSourceEmin(G4double Emin)
{
  theAdjointPrimaryGeneratorAction->SetEmin(Emin);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointSourceEmax(G4double Emax)
{
  theAdjointPrimaryGeneratorAction->SetEmax(Emax);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::ConsiderParticleAsPrimary(const G4String& particle_name)
{
  theAdjointPrimaryGeneratorAction->ConsiderParticleAsPrimary(particle_name);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::NeglectParticleAsPrimary(const G4String& particle_name)
{
  theAdjointPrimaryGeneratorAction->NeglectParticleAsPrimary(particle_name);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetPrimaryIon(G4ParticleDefinition* adjointIon,
                                        G4ParticleDefinition* fwdIon)
{
  theAdjointPrimaryGeneratorAction->SetPrimaryIon(adjointIon, fwdIon);
}
 
// --------------------------------------------------------------------
//
const G4String& G4AdjointSimManager::GetPrimaryIonName()
{
  return theAdjointPrimaryGeneratorAction->GetPrimaryIonName();
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::RegisterAdjointPrimaryWeight(G4double aWeight)
{
  theAdjointPrimaryWeight = aWeight;
  theAdjointSteppingAction->SetPrimWeight(aWeight);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointEventAction(G4UserEventAction* anAction)
{
  theAdjointEventAction = anAction;
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointSteppingAction(G4UserSteppingAction* anAction)
{
  theAdjointSteppingAction->SetUserAdjointSteppingAction(anAction);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointStackingAction(G4UserStackingAction* anAction)
{
  theAdjointStackingAction->SetUserAdjointStackingAction(anAction);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetAdjointRunAction(G4UserRunAction* anAction)
{
  theAdjointRunAction = anAction;
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetNbOfPrimaryFwdGammasPerEvent(G4int nb)
{
  theAdjointPrimaryGeneratorAction->SetNbPrimaryFwdGammasPerEvent(nb);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetNbAdjointPrimaryGammasPerEvent(G4int nb)
{
  theAdjointPrimaryGeneratorAction->SetNbAdjointPrimaryGammasPerEvent(nb);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::SetNbAdjointPrimaryElectronsPerEvent(G4int nb)
{
  theAdjointPrimaryGeneratorAction->SetNbAdjointPrimaryElectronsPerEvent(nb);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::BeginOfRunAction(const G4Run* aRun)
{  if (!adjoint_sim_mode){
    if(fUserRunAction  != nullptr) fUserRunAction->BeginOfRunAction(aRun);
   }
   else {
    if (theAdjointRunAction  != nullptr) theAdjointRunAction->BeginOfRunAction(aRun);
   }
  
  //fUserRunAction->BeginOfRunAction(aRun);
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::EndOfRunAction(const G4Run* aRun)
{
  if (!adjoint_sim_mode) {
    if (fUserRunAction != nullptr) fUserRunAction->EndOfRunAction(aRun);
  }
  else if (theAdjointRunAction != nullptr)
    theAdjointRunAction->EndOfRunAction(aRun);
 
  BackToFwdSimulationMode();
 
}
 
// --------------------------------------------------------------------
//
G4ParticleDefinition* G4AdjointSimManager::GetLastGeneratedFwdPrimaryParticle()
{
  return theAdjointPrimaryGeneratorAction->GetLastGeneratedFwdPrimaryParticle();
}
 
// --------------------------------------------------------------------
//
void G4AdjointSimManager::ResetDidOneAdjPartReachExtSourceDuringEvent()
{
  theAdjointSteppingAction->ResetDidOneAdjPartReachExtSourceDuringEvent();
}