G4AdjointSimManager Class Reference

#include <G4AdjointSimManager.hh>

Inheritance diagram for G4AdjointSimManager:

Collaboration diagram for G4AdjointSimManager:

Public Member Functions
virtual void	BeginOfRunAction (const G4Run *aRun)
virtual void	EndOfRunAction (const G4Run *aRun)
void	RunAdjointSimulation (G4int nb_evt)
G4int	GetNbEvtOfLastRun ()
void	SetAdjointTrackingMode (G4bool aBool)
G4bool	GetAdjointTrackingMode ()
G4bool	GetAdjointSimMode ()
G4bool	GetDidAdjParticleReachTheExtSource ()
void	RegisterAtEndOfAdjointTrack ()
void	RegisterAdjointPrimaryWeight (G4double aWeight)
G4int	GetIDOfLastAdjParticleReachingExtSource ()
G4ThreeVector	GetPositionAtEndOfLastAdjointTrack ()
G4ThreeVector	GetDirectionAtEndOfLastAdjointTrack ()
G4double	GetEkinAtEndOfLastAdjointTrack ()
G4double	GetEkinNucAtEndOfLastAdjointTrack ()
G4double	GetWeightAtEndOfLastAdjointTrack ()
G4double	GetCosthAtEndOfLastAdjointTrack ()
const G4String &	GetFwdParticleNameAtEndOfLastAdjointTrack ()
G4int	GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack ()
G4int	GetFwdParticleIndexAtEndOfLastAdjointTrack ()
std::vector< G4ParticleDefinition * > *	GetListOfPrimaryFwdParticles ()
size_t	GetNbOfPrimaryFwdParticles ()
G4bool	DefineSphericalExtSource (G4double radius, G4ThreeVector pos)
G4bool	DefineSphericalExtSourceWithCentreAtTheCentreOfAVolume (G4double radius, const G4String &volume_name)
G4bool	DefineExtSourceOnTheExtSurfaceOfAVolume (const G4String &volume_name)
void	SetExtSourceEmax (G4double Emax)
G4bool	DefineSphericalAdjointSource (G4double radius, G4ThreeVector pos)
G4bool	DefineSphericalAdjointSourceWithCentreAtTheCentreOfAVolume (G4double radius, const G4String &volume_name)
G4bool	DefineAdjointSourceOnTheExtSurfaceOfAVolume (const G4String &volume_name)
void	SetAdjointSourceEmin (G4double Emin)
void	SetAdjointSourceEmax (G4double Emax)
G4double	GetAdjointSourceArea ()
void	ConsiderParticleAsPrimary (const G4String &particle_name)
void	NeglectParticleAsPrimary (const G4String &particle_name)
void	SetPrimaryIon (G4ParticleDefinition *adjointIon, G4ParticleDefinition *fwdIon)
const G4String &	GetPrimaryIonName ()
void	SetNormalisationMode (G4int n)
G4int	GetNormalisationMode ()
G4double	GetNumberNucleonsInIon ()
void	SetAdjointEventAction (G4UserEventAction *anAction)
void	SetAdjointSteppingAction (G4UserSteppingAction *anAction)
void	SetAdjointStackingAction (G4UserStackingAction *anAction)
void	SetAdjointRunAction (G4UserRunAction *anAction)
void	UseUserStackingActionInFwdTrackingPhase (G4bool aBool)
void	UseUserTrackingActionInFwdTrackingPhase (G4bool aBool)
void	SetNbOfPrimaryFwdGammasPerEvent (G4int)
void	SwitchToAdjointSimulationMode ()
void	BackToFwdSimulationMode ()
Public Member Functions inherited from G4UserRunAction
	G4UserRunAction ()
virtual	~G4UserRunAction ()
virtual G4Run *	GenerateRun ()
void	SetMaster (G4bool val=true)
G4bool	IsMaster () const

Static Public Member Functions
static G4AdjointSimManager *	GetInstance ()

Private Member Functions
void	SetRestOfAdjointActions ()
void	SetAdjointPrimaryRunAndStackingActions ()
void	SetAdjointActions ()
void	ResetRestOfUserActions ()
void	ResetUserPrimaryRunAndStackingActions ()
void	ResetUserActions ()
void	DefineUserActions ()
	G4AdjointSimManager ()
	~G4AdjointSimManager ()

Private Attributes
G4AdjointSimMessenger *	theMessenger
bool	user_action_already_defined
G4UserRunAction *	fUserRunAction
G4UserEventAction *	fUserEventAction
G4VUserPrimaryGeneratorAction *	fUserPrimaryGeneratorAction
G4UserTrackingAction *	fUserTrackingAction
G4UserSteppingAction *	fUserSteppingAction
G4UserStackingAction *	fUserStackingAction
bool	use_user_StackingAction
bool	use_user_TrackingAction
G4UserRunAction *	theAdjointRunAction
G4UserEventAction *	theAdjointEventAction
G4AdjointPrimaryGeneratorAction *	theAdjointPrimaryGeneratorAction
G4AdjointTrackingAction *	theAdjointTrackingAction
G4AdjointSteppingAction *	theAdjointSteppingAction
G4AdjointStackingAction *	theAdjointStackingAction
G4bool	adjoint_tracking_mode
G4bool	adjoint_sim_mode
G4ThreeVector	last_pos
G4ThreeVector	last_direction
G4double	last_ekin
G4double	last_ekin_nuc
G4double	last_cos_th
G4String	last_fwd_part_name
G4int	last_fwd_part_PDGEncoding
G4int	last_fwd_part_index
G4double	last_weight
G4int	ID_of_last_particle_that_reach_the_ext_source
G4int	nb_evt_of_last_run
G4int	normalisation_mode
G4double	area_of_the_adjoint_source
G4double	nb_nuc
G4double	theAdjointPrimaryWeight
G4bool	welcome_message

Static Private Attributes
static G4ThreadLocal G4AdjointSimManager *	instance = 0

Additional Inherited Members
Protected Attributes inherited from G4UserRunAction
G4bool	isMaster

Detailed Description

Definition at line 150 of file G4AdjointSimManager.hh.

Constructor & Destructor Documentation

G4AdjointSimManager()

G4AdjointSimManager::G4AdjointSimManager ( )

private

Definition at line 70 of file G4AdjointSimManager.cc.

                                        :
   fUserRunAction(0), fUserEventAction(0),fUserPrimaryGeneratorAction(0),
   fUserTrackingAction(0), fUserSteppingAction(0), fUserStackingAction(0),
   theAdjointRunAction(0), theAdjointEventAction(0),
   adjoint_tracking_mode(false),last_ekin(0),last_ekin_nuc(0),
   last_cos_th(0),last_fwd_part_PDGEncoding(0),last_fwd_part_index(0),
   last_weight(0), ID_of_last_particle_that_reach_the_ext_source(0),
   nb_evt_of_last_run(0),area_of_the_adjoint_source(0),theAdjointPrimaryWeight(0)
{ 
 //Create adjoint actions;
 //----------------------
 theAdjointPrimaryGeneratorAction = new G4AdjointPrimaryGeneratorAction(); 
 theAdjointSteppingAction = new G4AdjointSteppingAction();
 theAdjointStackingAction = new G4AdjointStackingAction();
 theAdjointTrackingAction = new G4AdjointTrackingAction(theAdjointSteppingAction);
 theAdjointTrackingAction->SetListOfPrimaryFwdParticles(
 theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles());
 //Create messenger
 //----------------
  theMessenger = new G4AdjointSimMessenger(this);
  
  user_action_already_defined=false;
  use_user_StackingAction = false; 
  use_user_TrackingAction =true;

  adjoint_sim_mode = false;
  
  normalisation_mode=3;
  
  nb_nuc=1.;
  
  welcome_message =true;
  
  //Define user action and set this class instance  as RunAction
  //----------------
  DefineUserActions();
  G4RunManager* theRunManager = G4RunManager::GetRunManager();

  theRunManager->G4RunManager::SetUserAction(this);
/*
#ifdef G4MULTITHREADED

 if (theRunManager->GetRunManagerType() == G4RunManager::workerRM){
     G4cout<<"Here"<<std::endl;
     //G4MTAdjointSimManager::GetInstance()->RegisterLocalManager(this);
     G4cout<<"Here1"<<std::endl;
 }
#endif
*/
}

Here is the call graph for this function:

Here is the caller graph for this function:

~G4AdjointSimManager()

G4AdjointSimManager::~G4AdjointSimManager ( )

private

Definition at line 122 of file G4AdjointSimManager.cc.

{ 
  if (theAdjointRunAction) delete theAdjointRunAction;
  if (theAdjointPrimaryGeneratorAction) delete theAdjointPrimaryGeneratorAction;
  if (theAdjointSteppingAction) delete theAdjointSteppingAction;
  if (theAdjointEventAction) delete theAdjointEventAction;
  if (theAdjointTrackingAction) delete theAdjointTrackingAction;
  if (theAdjointStackingAction) delete theAdjointStackingAction;
  if (theMessenger) delete theMessenger;
}

Member Function Documentation

BackToFwdSimulationMode()

void G4AdjointSimManager::BackToFwdSimulationMode

(

)

Definition at line 221 of file G4AdjointSimManager.cc.

{  //Restore the user defined actions
   //--------------------------------
   ResetUserActions();
   adjoint_sim_mode=false;
}

Here is the call graph for this function:

Here is the caller graph for this function:

BeginOfRunAction()

void G4AdjointSimManager::BeginOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

Definition at line 620 of file G4AdjointSimManager.cc.

{

 if (!adjoint_sim_mode){
  if(fUserRunAction) fUserRunAction->BeginOfRunAction(aRun);
 }
 else {
  if (theAdjointRunAction) theAdjointRunAction->BeginOfRunAction(aRun);
 }
}

Here is the call graph for this function:

ConsiderParticleAsPrimary()

void G4AdjointSimManager::ConsiderParticleAsPrimary

(

const G4String &

particle_name

)

Definition at line 550 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->ConsiderParticleAsPrimary(particle_name);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefineAdjointSourceOnTheExtSurfaceOfAVolume()

G4bool G4AdjointSimManager::DefineAdjointSourceOnTheExtSurfaceOfAVolume

(

const G4String &

volume_name

)

Definition at line 526 of file G4AdjointSimManager.cc.

{
   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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefineExtSourceOnTheExtSurfaceOfAVolume()

G4bool G4AdjointSimManager::DefineExtSourceOnTheExtSurfaceOfAVolume

(

const G4String &

volume_name

)

Definition at line 492 of file G4AdjointSimManager.cc.

{
   G4double area;
   return G4AdjointCrossSurfChecker::GetInstance()->AddanExtSurfaceOfAvolume( "ExternalSource", volume_name,area);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefineSphericalAdjointSource()

G4bool G4AdjointSimManager::DefineSphericalAdjointSource	(	G4double	radius,
		G4ThreeVector	pos
	)

Definition at line 505 of file G4AdjointSimManager.cc.

{   
   G4double area;
   G4bool aBool = G4AdjointCrossSurfChecker::GetInstance()->AddaSphericalSurface("AdjointSource", radius, pos, area); 
   theAdjointPrimaryGeneratorAction->SetSphericalAdjointPrimarySource(radius, pos); 
   area_of_the_adjoint_source=area;
   return aBool;    
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefineSphericalAdjointSourceWithCentreAtTheCentreOfAVolume()

G4bool G4AdjointSimManager::DefineSphericalAdjointSourceWithCentreAtTheCentreOfAVolume	(	G4double	radius,
		const G4String &	volume_name
	)

Definition at line 515 of file G4AdjointSimManager.cc.

{
   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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefineSphericalExtSource()

G4bool G4AdjointSimManager::DefineSphericalExtSource	(	G4double	radius,
		G4ThreeVector	pos
	)

Definition at line 477 of file G4AdjointSimManager.cc.

{   
   G4double area;
   return G4AdjointCrossSurfChecker::GetInstance()->AddaSphericalSurface("ExternalSource", radius, pos, area);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefineSphericalExtSourceWithCentreAtTheCentreOfAVolume()

G4bool G4AdjointSimManager::DefineSphericalExtSourceWithCentreAtTheCentreOfAVolume	(	G4double	radius,
		const G4String &	volume_name
	)

Definition at line 484 of file G4AdjointSimManager.cc.

{
   G4double area;
   G4ThreeVector center;
   return G4AdjointCrossSurfChecker::GetInstance()->AddaSphericalSurfaceWithCenterAtTheCenterOfAVolume( "ExternalSource", radius, volume_name,center, area);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefineUserActions()

void G4AdjointSimManager::DefineUserActions ( )

private

Definition at line 308 of file G4AdjointSimManager.cc.

{ 
   G4RunManager* theRunManager =  G4RunManager::GetRunManager();
   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() );
   fUserRunAction= const_cast<G4UserRunAction*>( theRunManager->GetUserRunAction() );
   fUserStackingAction= const_cast<G4UserStackingAction* >( theRunManager->GetUserStackingAction() );
   user_action_already_defined=true;    
}

Here is the call graph for this function:

Here is the caller graph for this function:

EndOfRunAction()

void G4AdjointSimManager::EndOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

Definition at line 632 of file G4AdjointSimManager.cc.

{if (!adjoint_sim_mode){
  if(fUserRunAction) fUserRunAction->EndOfRunAction(aRun);
 }
 else if (theAdjointRunAction) theAdjointRunAction->EndOfRunAction(aRun);
/*
#ifdef G4MULTITHREADED
 if (G4RunManager::GetRunManager()->GetRunManagerType() == G4RunManager::workerRM){
  if (adjoint_sim_mode) BackToFwdSimulationMode();
 }
#endif
*/

}

Here is the call graph for this function:

GetAdjointSimMode()

G4bool G4AdjointSimManager::GetAdjointSimMode ( )

inline

Definition at line 166 of file G4AdjointSimManager.hh.

{return  adjoint_sim_mode;} //true if an adjoint simulation is running

Here is the call graph for this function:

Here is the caller graph for this function:

GetAdjointSourceArea()

G4double G4AdjointSimManager::GetAdjointSourceArea ( )

inline

Definition at line 199 of file G4AdjointSimManager.hh.

{return area_of_the_adjoint_source;} 

Here is the call graph for this function:

Here is the caller graph for this function:

GetAdjointTrackingMode()

G4bool G4AdjointSimManager::GetAdjointTrackingMode

(

)

Definition at line 322 of file G4AdjointSimManager.cc.

                                                  {
  return  theAdjointTrackingAction->GetIsAdjointTrackingMode();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetCosthAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetCosthAtEndOfLastAdjointTrack

(

)

Definition at line 395 of file G4AdjointSimManager.cc.

                                                             {
 return theAdjointTrackingAction->GetCosthAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetDidAdjParticleReachTheExtSource()

G4bool G4AdjointSimManager::GetDidAdjParticleReachTheExtSource

(

)

Definition at line 350 of file G4AdjointSimManager.cc.

{
  return theAdjointSteppingAction->GetDidAdjParticleReachTheExtSource();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetDirectionAtEndOfLastAdjointTrack()

G4ThreeVector G4AdjointSimManager::GetDirectionAtEndOfLastAdjointTrack

(

)

Definition at line 375 of file G4AdjointSimManager.cc.

                                                                      {
   return theAdjointTrackingAction->GetDirectionAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetEkinAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetEkinAtEndOfLastAdjointTrack

(

)

Definition at line 380 of file G4AdjointSimManager.cc.

                                                            {
 return theAdjointTrackingAction->GetEkinAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetEkinNucAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetEkinNucAtEndOfLastAdjointTrack

(

)

Definition at line 385 of file G4AdjointSimManager.cc.

                                                               {
 return theAdjointTrackingAction->GetEkinNucAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetFwdParticleIndexAtEndOfLastAdjointTrack()

G4int G4AdjointSimManager::GetFwdParticleIndexAtEndOfLastAdjointTrack

(

)

Definition at line 411 of file G4AdjointSimManager.cc.

{return theAdjointTrackingAction->GetLastFwdParticleIndex();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetFwdParticleNameAtEndOfLastAdjointTrack()

const G4String & G4AdjointSimManager::GetFwdParticleNameAtEndOfLastAdjointTrack

(

)

Definition at line 400 of file G4AdjointSimManager.cc.

{return theAdjointTrackingAction->GetFwdParticleNameAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack()

G4int G4AdjointSimManager::GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack

(

)

Definition at line 405 of file G4AdjointSimManager.cc.

                                                                           {
 return theAdjointTrackingAction->GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetIDOfLastAdjParticleReachingExtSource()

G4int G4AdjointSimManager::GetIDOfLastAdjParticleReachingExtSource ( )

inline

Definition at line 172 of file G4AdjointSimManager.hh.

{return ID_of_last_particle_that_reach_the_ext_source;};

Here is the call graph for this function:

GetInstance()

G4AdjointSimManager * G4AdjointSimManager::GetInstance ( void  )

static

Definition at line 134 of file G4AdjointSimManager.cc.

{
  if (instance == 0) instance = new G4AdjointSimManager;
  return instance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetListOfPrimaryFwdParticles()

std::vector< G4ParticleDefinition * > * G4AdjointSimManager::GetListOfPrimaryFwdParticles

(

)

Definition at line 356 of file G4AdjointSimManager.cc.

{
  return theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetNbEvtOfLastRun()

G4int G4AdjointSimManager::GetNbEvtOfLastRun ( )

inline

Definition at line 162 of file G4AdjointSimManager.hh.

{return nb_evt_of_last_run;}

Here is the call graph for this function:

Here is the caller graph for this function:

GetNbOfPrimaryFwdParticles()

size_t G4AdjointSimManager::GetNbOfPrimaryFwdParticles

(

)

Definition at line 362 of file G4AdjointSimManager.cc.

{
  return theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles()->size();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetNormalisationMode()

G4int G4AdjointSimManager::GetNormalisationMode ( )

inline

Definition at line 206 of file G4AdjointSimManager.hh.

{return normalisation_mode;};

GetNumberNucleonsInIon()

G4double G4AdjointSimManager::GetNumberNucleonsInIon ( )

inline

Definition at line 207 of file G4AdjointSimManager.hh.

{return nb_nuc;};

Here is the call graph for this function:

GetPositionAtEndOfLastAdjointTrack()

G4ThreeVector G4AdjointSimManager::GetPositionAtEndOfLastAdjointTrack

(

)

Definition at line 369 of file G4AdjointSimManager.cc.

                                                                      {
   return theAdjointTrackingAction->GetPositionAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetPrimaryIonName()

const G4String & G4AdjointSimManager::GetPrimaryIonName

(

)

Definition at line 575 of file G4AdjointSimManager.cc.

{
  return theAdjointPrimaryGeneratorAction->GetPrimaryIonName();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetWeightAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetWeightAtEndOfLastAdjointTrack

(

)

Definition at line 390 of file G4AdjointSimManager.cc.

                                                              {
 return theAdjointTrackingAction->GetWeightAtEndOfLastAdjointTrack();
}

Here is the call graph for this function:

Here is the caller graph for this function:

NeglectParticleAsPrimary()

void G4AdjointSimManager::NeglectParticleAsPrimary

(

const G4String &

particle_name

)

Definition at line 556 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->NeglectParticleAsPrimary(particle_name);
}

Here is the call graph for this function:

Here is the caller graph for this function:

RegisterAdjointPrimaryWeight()

void G4AdjointSimManager::RegisterAdjointPrimaryWeight

(

G4double

aWeight

)

Definition at line 581 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryWeight = aWeight;
  theAdjointSteppingAction->SetPrimWeight(aWeight);
} 

Here is the call graph for this function:

Here is the caller graph for this function:

RegisterAtEndOfAdjointTrack()

void G4AdjointSimManager::RegisterAtEndOfAdjointTrack

(

)

Definition at line 417 of file G4AdjointSimManager.cc.

{
  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.remove(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;
  size_t i=0;
  while(i<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=double(aPartDef->GetBaryonNumber());
    last_ekin_nuc /=nb_nuc;
  }

  last_weight = theAdjointSteppingAction->GetLastWeight(); 
  
  /* G4PhysicsLogVector* theWeightVector=0;
  if (last_fwd_part_name =="e-")  theWeightVector=electron_last_weight_vector;
  else if (last_fwd_part_name =="gamma") theWeightVector=gamma_last_weight_vector;
  else if (last_fwd_part_name =="proton") theWeightVector=proton_last_weight_vector;
  
  if (theWeightVector){

    size_t ind =  size_t(std::log10(last_weight/theAdjointPrimaryWeight)*10. + 200);
    G4double low_val =theWeightVector->GetLowEdgeEnergy(ind);
    G4bool aBool = true;
    G4double bin_weight = theWeightVector->GetValue(low_val, aBool)+1.;
    theWeightVector->PutValue(ind, bin_weight);
  }
  */
  /*if ((last_weight/theAdjointPrimaryWeight)>1.) last_weight*=1000. ;
  else if ( (last_weight/theAdjointPrimaryWeight)>0.1) last_weight*=100. ;
  else if ( (last_weight/theAdjointPrimaryWeight)>0.01) last_weight*=10. ;*/
  
  
  //G4cout <<"Last Weight "<<last_weight<<'\t'<<theAdjointPrimaryWeight<<'\t'<<last_weight/theAdjointPrimaryWeight<<std::endl;
  /*if (last_weight/theAdjointPrimaryWeight >10.) {
    G4cout<<"Warning a weight increase by a factor : "<<last_weight/theAdjointPrimaryWeight<<std::endl;
  }
  */
 
   ID_of_last_particle_that_reach_the_ext_source++;
}

Here is the call graph for this function:

Here is the caller graph for this function:

ResetRestOfUserActions()

void G4AdjointSimManager::ResetRestOfUserActions ( )

private

Definition at line 283 of file G4AdjointSimManager.cc.

{
  G4RunManager* theRunManager =  G4RunManager::GetRunManager();
  
  //Restore the user defined actions
  //-------------------------------
 
  theRunManager->G4RunManager::SetUserAction(fUserEventAction);
  theRunManager->G4RunManager::SetUserAction(fUserSteppingAction);
  theRunManager->G4RunManager::SetUserAction(fUserTrackingAction);
}

Here is the call graph for this function:

Here is the caller graph for this function:

ResetUserActions()

void G4AdjointSimManager::ResetUserActions ( )

private

Definition at line 268 of file G4AdjointSimManager.cc.

{
  G4RunManager* theRunManager =  G4RunManager::GetRunManager();

  //Restore the user defined actions
  //-------------------------------

  theRunManager->G4RunManager::SetUserAction(fUserEventAction);
  theRunManager->G4RunManager::SetUserAction(fUserSteppingAction);
  theRunManager->G4RunManager::SetUserAction(fUserTrackingAction);
  theRunManager->G4RunManager::SetUserAction(fUserPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(fUserStackingAction);
}

Here is the call graph for this function:

Here is the caller graph for this function:

ResetUserPrimaryRunAndStackingActions()

void G4AdjointSimManager::ResetUserPrimaryRunAndStackingActions ( )

private

Definition at line 297 of file G4AdjointSimManager.cc.

{ 
  G4RunManager* theRunManager =  G4RunManager::GetRunManager();
  //Restore the user defined actions
  //-------------------------------
  theRunManager->G4RunManager::SetUserAction(fUserRunAction);
  theRunManager->G4RunManager::SetUserAction(fUserPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(fUserStackingAction);
}

Here is the call graph for this function:

RunAdjointSimulation()

void G4AdjointSimManager::RunAdjointSimulation

(

G4int

nb_evt

)

Definition at line 141 of file G4AdjointSimManager.cc.

{ if (G4RunManager::GetRunManager()->GetRunManagerType() != G4RunManager::sequentialRM) return; //only for sequential mode
  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;
  G4RunManager::GetRunManager()->BeamOn(nb_evt*theAdjointPrimaryGeneratorAction->GetNbOfAdjointPrimaryTypes());
  //G4RunManager::GetRunManager()->BeamOn(theAdjointPrimaryGeneratorAction->GetNbOfAdjointPrimaryTypes()*2*nb_evt);

  //Back to Fwd Simulation Mode
  //--------------------------------
  BackToFwdSimulationMode();

  /*
  //Register the weight vector
  //--------------------------
  std::ofstream FileOutputElectronWeight("ElectronWeight.txt", std::ios::out);
  FileOutputElectronWeight<<std::setiosflags(std::ios::scientific);
  FileOutputElectronWeight<<std::setprecision(6);
  G4bool aBool = electron_last_weight_vector->Store(FileOutputElectronWeight, true);
  FileOutputElectronWeight.close();
  
  std::ofstream FileOutputProtonWeight("ProtonWeight.txt", std::ios::out);
  FileOutputProtonWeight<<std::setiosflags(std::ios::scientific);
  FileOutputProtonWeight<<std::setprecision(6);
  aBool = proton_last_weight_vector->Store(FileOutputProtonWeight, true);
  FileOutputProtonWeight.close();
  
  std::ofstream FileOutputGammaWeight("GammaWeight.txt", std::ios::out);
  FileOutputGammaWeight<<std::setiosflags(std::ios::scientific);
  FileOutputGammaWeight<<std::setprecision(6);
  aBool = gamma_last_weight_vector->Store(FileOutputGammaWeight, true);
  FileOutputGammaWeight.close();
  */  
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAdjointActions()

void G4AdjointSimManager::SetAdjointActions ( )

private

Definition at line 230 of file G4AdjointSimManager.cc.

{  
  G4RunManager* theRunManager =  G4RunManager::GetRunManager();
  
  if (!user_action_already_defined) DefineUserActions();

 //Replace the user action by the adjoint actions
 //------------------------------------------------- 
  
  theRunManager->G4RunManager::SetUserAction(theAdjointPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointStackingAction);
  if (use_user_StackingAction)  theAdjointStackingAction->SetUserFwdStackingAction(fUserStackingAction);
  else theAdjointStackingAction->SetUserFwdStackingAction(0);
  theRunManager->G4RunManager::SetUserAction(theAdjointEventAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointSteppingAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointTrackingAction);
  if (use_user_TrackingAction)  theAdjointTrackingAction->SetUserForwardTrackingAction(fUserTrackingAction);
  else theAdjointTrackingAction->SetUserForwardTrackingAction(0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAdjointEventAction()

void G4AdjointSimManager::SetAdjointEventAction

(

G4UserEventAction *

anAction

)

Definition at line 589 of file G4AdjointSimManager.cc.

{
  theAdjointEventAction = anAction;
}

Here is the caller graph for this function:

SetAdjointPrimaryRunAndStackingActions()

void G4AdjointSimManager::SetAdjointPrimaryRunAndStackingActions ( )

private

Definition at line 251 of file G4AdjointSimManager.cc.

{  
  G4RunManager* theRunManager =  G4RunManager::GetRunManager();
  
  if (!user_action_already_defined) DefineUserActions();
  
 //Replace the user action by the adjoint actions
 //------------------------------------------------- 
  
  theRunManager->G4RunManager::SetUserAction(theAdjointRunAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointPrimaryGeneratorAction);
  theRunManager->G4RunManager::SetUserAction(theAdjointStackingAction);
  if (use_user_StackingAction)  theAdjointStackingAction->SetUserFwdStackingAction(fUserStackingAction);
  else theAdjointStackingAction->SetUserFwdStackingAction(0);
}

Here is the call graph for this function:

SetAdjointRunAction()

void G4AdjointSimManager::SetAdjointRunAction

(

G4UserRunAction *

anAction

)

Definition at line 608 of file G4AdjointSimManager.cc.

{
  theAdjointRunAction=anAction;
} 

Here is the caller graph for this function:

SetAdjointSourceEmax()

void G4AdjointSimManager::SetAdjointSourceEmax

(

G4double

Emax

)

Definition at line 544 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->SetEmax(Emax);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAdjointSourceEmin()

void G4AdjointSimManager::SetAdjointSourceEmin

(

G4double

Emin

)

Definition at line 538 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->SetEmin(Emin);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAdjointStackingAction()

void G4AdjointSimManager::SetAdjointStackingAction

(

G4UserStackingAction *

anAction

)

Definition at line 601 of file G4AdjointSimManager.cc.

{
  theAdjointStackingAction->SetUserAdjointStackingAction(anAction);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAdjointSteppingAction()

void G4AdjointSimManager::SetAdjointSteppingAction

(

G4UserSteppingAction *

anAction

)

Definition at line 595 of file G4AdjointSimManager.cc.

{
  theAdjointSteppingAction->SetUserAdjointSteppingAction(anAction);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAdjointTrackingMode()

void G4AdjointSimManager::SetAdjointTrackingMode

(

G4bool

aBool

)

Definition at line 327 of file G4AdjointSimManager.cc.

{
  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);
  } 
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetExtSourceEmax()

void G4AdjointSimManager::SetExtSourceEmax

(

G4double

Emax

)

Definition at line 499 of file G4AdjointSimManager.cc.

{
  theAdjointSteppingAction->SetExtSourceEMax(Emax);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetNbOfPrimaryFwdGammasPerEvent()

void G4AdjointSimManager::SetNbOfPrimaryFwdGammasPerEvent

(

G4int

nb

)

Definition at line 614 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->SetNbPrimaryFwdGammasPerEvent(nb);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetNormalisationMode()

void G4AdjointSimManager::SetNormalisationMode ( G4int  n )

inline

Definition at line 205 of file G4AdjointSimManager.hh.

{normalisation_mode=n;};

SetPrimaryIon()

void G4AdjointSimManager::SetPrimaryIon	(	G4ParticleDefinition *	adjointIon,
		G4ParticleDefinition *	fwdIon
	)

Definition at line 569 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->SetPrimaryIon(adjointIon, fwdIon);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetRestOfAdjointActions()

void G4AdjointSimManager::SetRestOfAdjointActions ( )

private

Definition at line 192 of file G4AdjointSimManager.cc.

{
  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);

}

Here is the call graph for this function:

Here is the caller graph for this function:

SwitchToAdjointSimulationMode()

void G4AdjointSimManager::SwitchToAdjointSimulationMode

(

)

Definition at line 208 of file G4AdjointSimManager.cc.

{ //Replace the user defined actions by the adjoint actions
  //---------------------------------------------------------
  SetAdjointActions();

  //Update the list of primaries
  //-----------------------------
  theAdjointPrimaryGeneratorAction->UpdateListOfPrimaryParticles();
  adjoint_sim_mode=true;
  ID_of_last_particle_that_reach_the_ext_source=0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

UseUserStackingActionInFwdTrackingPhase()

void G4AdjointSimManager::UseUserStackingActionInFwdTrackingPhase ( G4bool  aBool )

inline

Definition at line 218 of file G4AdjointSimManager.hh.

{use_user_StackingAction=aBool;}

UseUserTrackingActionInFwdTrackingPhase()

void G4AdjointSimManager::UseUserTrackingActionInFwdTrackingPhase ( G4bool  aBool )

inline

Definition at line 219 of file G4AdjointSimManager.hh.

{use_user_TrackingAction=aBool;}

Here is the call graph for this function:

Member Data Documentation

adjoint_sim_mode

G4bool G4AdjointSimManager::adjoint_sim_mode

private

Definition at line 288 of file G4AdjointSimManager.hh.

adjoint_tracking_mode

G4bool G4AdjointSimManager::adjoint_tracking_mode

private

Definition at line 287 of file G4AdjointSimManager.hh.

area_of_the_adjoint_source

G4double G4AdjointSimManager::area_of_the_adjoint_source

private

Definition at line 307 of file G4AdjointSimManager.hh.

fUserEventAction

G4UserEventAction* G4AdjointSimManager::fUserEventAction

private

Definition at line 268 of file G4AdjointSimManager.hh.

fUserPrimaryGeneratorAction

G4VUserPrimaryGeneratorAction* G4AdjointSimManager::fUserPrimaryGeneratorAction

private

Definition at line 269 of file G4AdjointSimManager.hh.

fUserRunAction

G4UserRunAction* G4AdjointSimManager::fUserRunAction

private

Definition at line 267 of file G4AdjointSimManager.hh.

fUserStackingAction

G4UserStackingAction* G4AdjointSimManager::fUserStackingAction

private

Definition at line 272 of file G4AdjointSimManager.hh.

fUserSteppingAction

G4UserSteppingAction* G4AdjointSimManager::fUserSteppingAction

private

Definition at line 271 of file G4AdjointSimManager.hh.

fUserTrackingAction

G4UserTrackingAction* G4AdjointSimManager::fUserTrackingAction

private

Definition at line 270 of file G4AdjointSimManager.hh.

ID_of_last_particle_that_reach_the_ext_source

G4int G4AdjointSimManager::ID_of_last_particle_that_reach_the_ext_source

private

Definition at line 300 of file G4AdjointSimManager.hh.

instance

G4ThreadLocal G4AdjointSimManager * G4AdjointSimManager::instance = 0

staticprivate

Definition at line 237 of file G4AdjointSimManager.hh.

last_cos_th

G4double G4AdjointSimManager::last_cos_th

private

Definition at line 295 of file G4AdjointSimManager.hh.

last_direction

G4ThreeVector G4AdjointSimManager::last_direction

private

Definition at line 293 of file G4AdjointSimManager.hh.

last_ekin

G4double G4AdjointSimManager::last_ekin

private

Definition at line 294 of file G4AdjointSimManager.hh.

last_ekin_nuc

G4double G4AdjointSimManager::last_ekin_nuc

private

Definition at line 294 of file G4AdjointSimManager.hh.

last_fwd_part_index

G4int G4AdjointSimManager::last_fwd_part_index

private

Definition at line 298 of file G4AdjointSimManager.hh.

last_fwd_part_name

G4String G4AdjointSimManager::last_fwd_part_name

private

Definition at line 296 of file G4AdjointSimManager.hh.

last_fwd_part_PDGEncoding

G4int G4AdjointSimManager::last_fwd_part_PDGEncoding

private

Definition at line 297 of file G4AdjointSimManager.hh.

last_pos

G4ThreeVector G4AdjointSimManager::last_pos

private

Definition at line 292 of file G4AdjointSimManager.hh.

last_weight

G4double G4AdjointSimManager::last_weight

private

Definition at line 299 of file G4AdjointSimManager.hh.

nb_evt_of_last_run

G4int G4AdjointSimManager::nb_evt_of_last_run

private

Definition at line 302 of file G4AdjointSimManager.hh.

nb_nuc

G4double G4AdjointSimManager::nb_nuc

private

Definition at line 308 of file G4AdjointSimManager.hh.

normalisation_mode

G4int G4AdjointSimManager::normalisation_mode

private

Definition at line 303 of file G4AdjointSimManager.hh.

theAdjointEventAction

G4UserEventAction* G4AdjointSimManager::theAdjointEventAction

private

Definition at line 279 of file G4AdjointSimManager.hh.

theAdjointPrimaryGeneratorAction

G4AdjointPrimaryGeneratorAction* G4AdjointSimManager::theAdjointPrimaryGeneratorAction

private

Definition at line 280 of file G4AdjointSimManager.hh.

theAdjointPrimaryWeight

G4double G4AdjointSimManager::theAdjointPrimaryWeight

private

Definition at line 309 of file G4AdjointSimManager.hh.

theAdjointRunAction

G4UserRunAction* G4AdjointSimManager::theAdjointRunAction

private

Definition at line 278 of file G4AdjointSimManager.hh.

theAdjointStackingAction

G4AdjointStackingAction* G4AdjointSimManager::theAdjointStackingAction

private

Definition at line 283 of file G4AdjointSimManager.hh.

theAdjointSteppingAction

G4AdjointSteppingAction* G4AdjointSimManager::theAdjointSteppingAction

private

Definition at line 282 of file G4AdjointSimManager.hh.

theAdjointTrackingAction

G4AdjointTrackingAction* G4AdjointSimManager::theAdjointTrackingAction

private

Definition at line 281 of file G4AdjointSimManager.hh.

theMessenger

G4AdjointSimMessenger* G4AdjointSimManager::theMessenger

private

Definition at line 262 of file G4AdjointSimManager.hh.

use_user_StackingAction

bool G4AdjointSimManager::use_user_StackingAction

private

Definition at line 273 of file G4AdjointSimManager.hh.

use_user_TrackingAction

bool G4AdjointSimManager::use_user_TrackingAction

private

Definition at line 274 of file G4AdjointSimManager.hh.

user_action_already_defined

bool G4AdjointSimManager::user_action_already_defined

private

Definition at line 266 of file G4AdjointSimManager.hh.

welcome_message

G4bool G4AdjointSimManager::welcome_message

private

Definition at line 317 of file G4AdjointSimManager.hh.

---

The documentation for this class was generated from the following files:

• G4AdjointSimManager.hh
• G4AdjointSimManager.cc