G4AdjointSimManager Class Reference

#include <G4AdjointSimManager.hh>

Public Member Functions
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 ()
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	SetAdjointTrackingAction (G4UserTrackingAction *anAction)
void	SetAdjointRunAction (G4UserRunAction *anAction)
void	UseUserStackingActionInFwdTrackingPhase (G4bool aBool)

Static Public Member Functions
static G4AdjointSimManager *	GetInstance ()

Detailed Description

Definition at line 149 of file G4AdjointSimManager.hh.

Member Function Documentation

ConsiderParticleAsPrimary()

void G4AdjointSimManager::ConsiderParticleAsPrimary

(

const G4String &

particle_name

)

Definition at line 422 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->ConsiderParticleAsPrimary(particle_name);
}

Referenced by G4AdjointSimMessenger::SetNewValue().

DefineAdjointSourceOnTheExtSurfaceOfAVolume()

G4bool G4AdjointSimManager::DefineAdjointSourceOnTheExtSurfaceOfAVolume

(

const G4String &

volume_name

)

Definition at line 398 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;
}

Referenced by G4AdjointSimMessenger::SetNewValue().

DefineExtSourceOnTheExtSurfaceOfAVolume()

G4bool G4AdjointSimManager::DefineExtSourceOnTheExtSurfaceOfAVolume

(

const G4String &

volume_name

)

Definition at line 364 of file G4AdjointSimManager.cc.

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

Referenced by G4AdjointSimMessenger::SetNewValue().

DefineSphericalAdjointSource()

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

Definition at line 377 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;    
}

Referenced by G4AdjointSimMessenger::SetNewValue().

DefineSphericalAdjointSourceWithCentreAtTheCentreOfAVolume()

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

Definition at line 387 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;
}

Referenced by G4AdjointSimMessenger::SetNewValue().

DefineSphericalExtSource()

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

Definition at line 349 of file G4AdjointSimManager.cc.

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

Referenced by G4AdjointSimMessenger::SetNewValue().

DefineSphericalExtSourceWithCentreAtTheCentreOfAVolume()

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

Definition at line 356 of file G4AdjointSimManager.cc.

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

Referenced by G4AdjointSimMessenger::SetNewValue().

GetAdjointSimMode()

G4bool G4AdjointSimManager::GetAdjointSimMode ( )

inline

Definition at line 163 of file G4AdjointSimManager.hh.

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

GetAdjointSourceArea()

G4double G4AdjointSimManager::GetAdjointSourceArea ( )

inline

Definition at line 195 of file G4AdjointSimManager.hh.

{return area_of_the_adjoint_source;} 

Referenced by G4AdjointPrimaryGeneratorAction::GeneratePrimaries().

GetAdjointTrackingMode()

G4bool G4AdjointSimManager::GetAdjointTrackingMode ( )

inline

Definition at line 162 of file G4AdjointSimManager.hh.

{return  adjoint_tracking_mode;} //true if an adjoint track is being processed

GetCosthAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetCosthAtEndOfLastAdjointTrack ( )

inline

Definition at line 175 of file G4AdjointSimManager.hh.

{return last_cos_th;}

GetDidAdjParticleReachTheExtSource()

G4bool G4AdjointSimManager::GetDidAdjParticleReachTheExtSource

(

)

Definition at line 277 of file G4AdjointSimManager.cc.

{
  return theAdjointSteppingAction->GetDidAdjParticleReachTheExtSource();
}

Referenced by SetAdjointTrackingMode().

GetDirectionAtEndOfLastAdjointTrack()

G4ThreeVector G4AdjointSimManager::GetDirectionAtEndOfLastAdjointTrack ( )

inline

Definition at line 171 of file G4AdjointSimManager.hh.

{ return last_direction;}

GetEkinAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetEkinAtEndOfLastAdjointTrack ( )

inline

Definition at line 172 of file G4AdjointSimManager.hh.

{ return last_ekin;} 

GetEkinNucAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetEkinNucAtEndOfLastAdjointTrack ( )

inline

Definition at line 173 of file G4AdjointSimManager.hh.

{ return last_ekin_nuc;}

GetFwdParticleIndexAtEndOfLastAdjointTrack()

G4int G4AdjointSimManager::GetFwdParticleIndexAtEndOfLastAdjointTrack ( )

inline

Definition at line 178 of file G4AdjointSimManager.hh.

{return last_fwd_part_index;}

GetFwdParticleNameAtEndOfLastAdjointTrack()

const G4String & G4AdjointSimManager::GetFwdParticleNameAtEndOfLastAdjointTrack ( )

inline

Definition at line 176 of file G4AdjointSimManager.hh.

{return last_fwd_part_name;}

GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack()

G4int G4AdjointSimManager::GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack ( )

inline

Definition at line 177 of file G4AdjointSimManager.hh.

{return last_fwd_part_PDGEncoding;}

GetIDOfLastAdjParticleReachingExtSource()

G4int G4AdjointSimManager::GetIDOfLastAdjParticleReachingExtSource ( )

inline

Definition at line 169 of file G4AdjointSimManager.hh.

{return ID_of_last_particle_that_reach_the_ext_source;};                     

GetInstance()

G4AdjointSimManager * G4AdjointSimManager::GetInstance ( )

static

Definition at line 116 of file G4AdjointSimManager.cc.

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

Referenced by G4AdjointPrimaryGeneratorAction::GeneratePrimaries().

GetListOfPrimaryFwdParticles()

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

(

)

Definition at line 283 of file G4AdjointSimManager.cc.

{
  return theAdjointPrimaryGeneratorAction->GetListOfPrimaryFwdParticles();
}

GetNbEvtOfLastRun()

G4int G4AdjointSimManager::GetNbEvtOfLastRun ( )

inline

Definition at line 159 of file G4AdjointSimManager.hh.

{return nb_evt_of_last_run;}

GetNormalisationMode()

G4int G4AdjointSimManager::GetNormalisationMode ( )

inline

Definition at line 202 of file G4AdjointSimManager.hh.

{return normalisation_mode;};

GetNumberNucleonsInIon()

G4double G4AdjointSimManager::GetNumberNucleonsInIon ( )

inline

Definition at line 203 of file G4AdjointSimManager.hh.

{return nb_nuc;};

GetPositionAtEndOfLastAdjointTrack()

G4ThreeVector G4AdjointSimManager::GetPositionAtEndOfLastAdjointTrack ( )

inline

Definition at line 170 of file G4AdjointSimManager.hh.

{ return last_pos;}

GetPrimaryIonName()

const G4String & G4AdjointSimManager::GetPrimaryIonName

(

)

Definition at line 447 of file G4AdjointSimManager.cc.

{
  return theAdjointPrimaryGeneratorAction->GetPrimaryIonName();
}

GetWeightAtEndOfLastAdjointTrack()

G4double G4AdjointSimManager::GetWeightAtEndOfLastAdjointTrack ( )

inline

Definition at line 174 of file G4AdjointSimManager.hh.

{return last_weight;}

NeglectParticleAsPrimary()

void G4AdjointSimManager::NeglectParticleAsPrimary

(

const G4String &

particle_name

)

Definition at line 428 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->NeglectParticleAsPrimary(particle_name);
}

Referenced by G4AdjointSimMessenger::SetNewValue().

RegisterAdjointPrimaryWeight()

void G4AdjointSimManager::RegisterAdjointPrimaryWeight

(

G4double

aWeight

)

Definition at line 453 of file G4AdjointSimManager.cc.

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

Referenced by G4AdjointPrimaryGeneratorAction::GeneratePrimaries().

RegisterAtEndOfAdjointTrack()

void G4AdjointSimManager::RegisterAtEndOfAdjointTrack

(

)

Definition at line 289 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++;
}

Referenced by SetAdjointTrackingMode().

RunAdjointSimulation()

void G4AdjointSimManager::RunAdjointSimulation

(

G4int

nb_evt

)

Definition at line 123 of file G4AdjointSimManager.cc.

{ 
  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;
  } 
  
  //Replace the user defined actions by the adjoint actions
  //---------------------------------------------------------
  SetAdjointPrimaryRunAndStackingActions();
  SetRestOfAdjointActions();
  
  //Update the list of primaries
  //-----------------------------
  theAdjointPrimaryGeneratorAction->UpdateListOfPrimaryParticles();
 
  adjoint_sim_mode=true;
  
  ID_of_last_particle_that_reach_the_ext_source=0;
  
  //Make the run
  //------------
  
  nb_evt_of_last_run =nb_evt;
  G4RunManager::GetRunManager()->BeamOn(theAdjointPrimaryGeneratorAction->GetNbOfAdjointPrimaryTypes()*2*nb_evt);
 
  //Restore the user defined actions
  //--------------------------------
  ResetRestOfUserActions(); 
  ResetUserPrimaryRunAndStackingActions();
  adjoint_sim_mode=false; 
 
  /*
  //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();
  */  
}

Referenced by G4AdjointSimMessenger::SetNewValue().

SetAdjointEventAction()

void G4AdjointSimManager::SetAdjointEventAction

(

G4UserEventAction *

anAction

)

Definition at line 461 of file G4AdjointSimManager.cc.

{
  theAdjointEventAction = anAction;
}

SetAdjointRunAction()

void G4AdjointSimManager::SetAdjointRunAction

(

G4UserRunAction *

anAction

)

Definition at line 485 of file G4AdjointSimManager.cc.

{
  theAdjointRunAction=anAction;
} 

SetAdjointSourceEmax()

void G4AdjointSimManager::SetAdjointSourceEmax

(

G4double

Emax

)

Definition at line 416 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->SetEmax(Emax);
}

Referenced by G4AdjointSimMessenger::SetNewValue().

SetAdjointSourceEmin()

void G4AdjointSimManager::SetAdjointSourceEmin

(

G4double

Emin

)

Definition at line 410 of file G4AdjointSimManager.cc.

{
  theAdjointPrimaryGeneratorAction->SetEmin(Emin);
}

Referenced by G4AdjointSimMessenger::SetNewValue().

SetAdjointStackingAction()

void G4AdjointSimManager::SetAdjointStackingAction

(

G4UserStackingAction *

anAction

)

Definition at line 473 of file G4AdjointSimManager.cc.

{
  theAdjointStackingAction->SetUserAdjointStackingAction(anAction);
}

SetAdjointSteppingAction()

void G4AdjointSimManager::SetAdjointSteppingAction

(

G4UserSteppingAction *

anAction

)

Definition at line 467 of file G4AdjointSimManager.cc.

{
  theAdjointSteppingAction->SetUserAdjointSteppingAction(anAction);
}

SetAdjointTrackingAction()

void G4AdjointSimManager::SetAdjointTrackingAction

(

G4UserTrackingAction *

anAction

)

Definition at line 479 of file G4AdjointSimManager.cc.

{
  theAdjointTrackingAction=anAction;
}

SetAdjointTrackingMode()

void G4AdjointSimManager::SetAdjointTrackingMode

(

G4bool

aBool

)

Definition at line 254 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);
  } 
}

Referenced by G4AdjointPrimaryGeneratorAction::GeneratePrimaries().

SetExtSourceEmax()

void G4AdjointSimManager::SetExtSourceEmax

(

G4double

Emax

)

Definition at line 371 of file G4AdjointSimManager.cc.

{
  theAdjointSteppingAction->SetExtSourceEMax(Emax);
}

Referenced by G4AdjointSimMessenger::SetNewValue().

SetNormalisationMode()

void G4AdjointSimManager::SetNormalisationMode ( G4int  n )

inline

Definition at line 201 of file G4AdjointSimManager.hh.

{normalisation_mode=n;};

SetPrimaryIon()

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

Definition at line 441 of file G4AdjointSimManager.cc.

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

UseUserStackingActionInFwdTrackingPhase()

void G4AdjointSimManager::UseUserStackingActionInFwdTrackingPhase ( G4bool  aBool )

inline

Definition at line 215 of file G4AdjointSimManager.hh.

{use_user_StackingAction=aBool;}

---

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

• G4AdjointSimManager.hh
• G4AdjointSimManager.cc