G4AdjointSteppingAction Class Reference

#include <G4AdjointSteppingAction.hh>

Inheritance diagram for G4AdjointSteppingAction:

Public Member Functions
	G4AdjointSteppingAction ()
	~G4AdjointSteppingAction ()
void	UserSteppingAction (const G4Step *)
void	SetExtSourceEMax (G4double Emax)
void	SetStartEvent (G4bool aBool)
G4bool	GetDidAdjParticleReachTheExtSource ()
G4ThreeVector	GetLastMomentum ()
G4ThreeVector	GetLastPosition ()
G4double	GetLastEkin ()
G4double	GetLastWeight ()
void	SetPrimWeight (G4double weight)
G4ParticleDefinition *	GetLastPartDef ()
void	SetUserAdjointSteppingAction (G4UserSteppingAction *anAction)
Public Member Functions inherited from G4UserSteppingAction
	G4UserSteppingAction ()
virtual	~G4UserSteppingAction ()
void	SetSteppingManagerPointer (G4SteppingManager *pValue)
virtual void	UserSteppingAction (const G4Step *)

Additional Inherited Members
Protected Attributes inherited from G4UserSteppingAction
G4SteppingManager *	fpSteppingManager

Detailed Description

Definition at line 69 of file G4AdjointSteppingAction.hh.

Constructor & Destructor Documentation

G4AdjointSteppingAction()

G4AdjointSteppingAction::G4AdjointSteppingAction

(

)

Definition at line 43 of file G4AdjointSteppingAction.cc.

  : ext_sourceEMax(0.), start_event(false),
    did_adj_part_reach_ext_source(false), last_ekin(0.), last_weight(0.),
    prim_weight(0.), last_part_def(0), theUserAdjointSteppingAction(0)
{ 
  theG4AdjointCrossSurfChecker = G4AdjointCrossSurfChecker::GetInstance();
}

~G4AdjointSteppingAction()

G4AdjointSteppingAction::~G4AdjointSteppingAction

(

)

Definition at line 52 of file G4AdjointSteppingAction.cc.

{;}

Member Function Documentation

GetDidAdjParticleReachTheExtSource()

G4bool G4AdjointSteppingAction::GetDidAdjParticleReachTheExtSource ( )

inline

Definition at line 79 of file G4AdjointSteppingAction.hh.

{return did_adj_part_reach_ext_source;}

Referenced by G4AdjointSimManager::GetDidAdjParticleReachTheExtSource().

GetLastEkin()

G4double G4AdjointSteppingAction::GetLastEkin ( )

inline

Definition at line 82 of file G4AdjointSteppingAction.hh.

{return last_ekin;}

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

GetLastMomentum()

G4ThreeVector G4AdjointSteppingAction::GetLastMomentum ( )

inline

Definition at line 80 of file G4AdjointSteppingAction.hh.

{return last_momentum;}

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

GetLastPartDef()

G4ParticleDefinition * G4AdjointSteppingAction::GetLastPartDef ( )

inline

Definition at line 85 of file G4AdjointSteppingAction.hh.

{return last_part_def;}

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

GetLastPosition()

G4ThreeVector G4AdjointSteppingAction::GetLastPosition ( )

inline

Definition at line 81 of file G4AdjointSteppingAction.hh.

{return last_pos;}

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

GetLastWeight()

G4double G4AdjointSteppingAction::GetLastWeight ( )

inline

Definition at line 83 of file G4AdjointSteppingAction.hh.

{return last_weight;}

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

SetExtSourceEMax()

void G4AdjointSteppingAction::SetExtSourceEMax ( G4double  Emax )

inline

Definition at line 77 of file G4AdjointSteppingAction.hh.

{ext_sourceEMax=Emax;} 

Referenced by G4AdjointSimManager::SetExtSourceEmax().

SetPrimWeight()

void G4AdjointSteppingAction::SetPrimWeight ( G4double  weight )

inline

Definition at line 84 of file G4AdjointSteppingAction.hh.

{prim_weight=weight;} 

Referenced by G4AdjointSimManager::RegisterAdjointPrimaryWeight().

SetStartEvent()

void G4AdjointSteppingAction::SetStartEvent ( G4bool  aBool )

inline

Definition at line 78 of file G4AdjointSteppingAction.hh.

{start_event =aBool;}

SetUserAdjointSteppingAction()

void G4AdjointSteppingAction::SetUserAdjointSteppingAction ( G4UserSteppingAction *  anAction )

inline

Definition at line 86 of file G4AdjointSteppingAction.hh.

{theUserAdjointSteppingAction = anAction;}

Referenced by G4AdjointSimManager::SetAdjointSteppingAction().

UserSteppingAction()

void G4AdjointSteppingAction::UserSteppingAction ( const G4Step *  aStep )

virtual

Reimplemented from G4UserSteppingAction.

Definition at line 57 of file G4AdjointSteppingAction.cc.

{
  //Apply first the user adjoint stepping action
  //---------------------------
  if (theUserAdjointSteppingAction) theUserAdjointSteppingAction->UserSteppingAction(aStep);
 
  G4Track* aTrack =aStep->GetTrack();
  G4double nb_nuc=1.;
  G4ParticleDefinition* thePartDef = aTrack->GetDefinition();
 
  if (thePartDef->GetParticleType() == "adjoint_nucleus"){
    nb_nuc=double(thePartDef->GetBaryonNumber());
  }
  //Kill conditions for adjoint particles reaching the maximum energy
  //-----------------------------------------------------------------
  if(aTrack->GetKineticEnergy() >= ext_sourceEMax*nb_nuc){
    aTrack->SetTrackStatus(fStopAndKill);
    did_adj_part_reach_ext_source=false;
    return;
  }
 
  G4double weight_factor = aTrack->GetWeight()/prim_weight;
 
  if ( (weight_factor>0 && weight_factor<=0) || weight_factor<= 1e-290 || weight_factor>1.e200)
  {
    //std::cout<<"Weight_factor problem! Value = "<<weight_factor<<std::endl;
    aTrack->SetTrackStatus(fStopAndKill);
    did_adj_part_reach_ext_source=false;
    return; 
  }
  
  
  //Kill conditions for surface crossing
  //--------------------------------------
  
  G4String surface_name;
  G4double cos_to_surface;
  G4bool GoingIn;
  G4ThreeVector crossing_pos;
  if (theG4AdjointCrossSurfChecker->CrossingOneOfTheRegisteredSurface(aStep, surface_name, crossing_pos, cos_to_surface, GoingIn) ){
    
    //G4cout<<"Test_step11"<<std::endl;
    //G4cout<<surface_name<<std::endl;
    if (surface_name == "ExternalSource") {
        //Registering still needed
        did_adj_part_reach_ext_source=true;
        aTrack->SetTrackStatus(fStopAndKill);
        //now register the adjoint particles reaching the external surface
        last_momentum =aTrack->GetMomentum();
        last_ekin=aTrack->GetKineticEnergy();
        last_weight = aTrack->GetWeight();
        last_part_def = aTrack->GetDefinition();
        last_pos = crossing_pos;
            return;
    }   
    else if (surface_name == "AdjointSource" && GoingIn) {
        did_adj_part_reach_ext_source=false;
        aTrack->SetTrackStatus(fStopAndKill);
        return;
    }  
  }
  //Check for reaching out of world
  if (aStep->GetPostStepPoint()->GetStepStatus() == fWorldBoundary) {
      did_adj_part_reach_ext_source=true;
      last_momentum =aTrack->GetMomentum();
      last_ekin=aTrack->GetKineticEnergy();
      last_weight = aTrack->GetWeight();
      last_part_def = aTrack->GetDefinition();
      last_pos = crossing_pos;
        return;
 
  }
  
}

---

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

• G4AdjointSteppingAction.hh
• G4AdjointSteppingAction.cc