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)
void	SetUserForwardSteppingAction (G4UserSteppingAction *anAction)
void	SetAdjointTrackingMode (G4bool aBool)
Public Member Functions inherited from G4UserSteppingAction
	G4UserSteppingAction ()
virtual	~G4UserSteppingAction ()
void	SetSteppingManagerPointer (G4SteppingManager *pValue)

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

(

)

Definition at line 44 of file G4AdjointSteppingAction.cc.

References G4AdjointCrossSurfChecker::GetInstance().

  : 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),
    theUserFwdSteppingAction(0)
{ 
  theG4AdjointCrossSurfChecker = G4AdjointCrossSurfChecker::GetInstance();
}

G4AdjointSteppingAction::~G4AdjointSteppingAction

(

)

Definition at line 54 of file G4AdjointSteppingAction.cc.

{;}

Member Function Documentation

G4bool G4AdjointSteppingAction::GetDidAdjParticleReachTheExtSource ( )

inline

Definition at line 79 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::GetDidAdjParticleReachTheExtSource(), and G4AdjointTrackingAction::PostUserTrackingAction().

{return did_adj_part_reach_ext_source;}

G4double G4AdjointSteppingAction::GetLastEkin ( )

inline

Definition at line 82 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointTrackingAction::PostUserTrackingAction(), and G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_ekin;}

G4ThreeVector G4AdjointSteppingAction::GetLastMomentum ( )

inline

Definition at line 80 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointTrackingAction::PostUserTrackingAction(), and G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_momentum;}

G4ParticleDefinition* G4AdjointSteppingAction::GetLastPartDef ( )

inline

Definition at line 85 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointTrackingAction::PostUserTrackingAction(), and G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_part_def;}

G4ThreeVector G4AdjointSteppingAction::GetLastPosition ( )

inline

Definition at line 81 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointTrackingAction::PostUserTrackingAction(), and G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_pos;}

G4double G4AdjointSteppingAction::GetLastWeight ( )

inline

Definition at line 83 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointTrackingAction::PostUserTrackingAction(), and G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_weight;}

void G4AdjointSteppingAction::SetAdjointTrackingMode ( G4bool  aBool )

inline

Definition at line 90 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointTrackingAction::PreUserTrackingAction().

{is_adjoint_tracking_mode =aBool;}

void G4AdjointSteppingAction::SetExtSourceEMax ( G4double  Emax )

inline

Definition at line 77 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::SetExtSourceEmax().

{ext_sourceEMax=Emax;} 

void G4AdjointSteppingAction::SetPrimWeight ( G4double  weight )

inline

Definition at line 84 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointTrackingAction::PreUserTrackingAction(), and G4AdjointSimManager::RegisterAdjointPrimaryWeight().

{prim_weight=weight;} 

void G4AdjointSteppingAction::SetStartEvent ( G4bool  aBool )

inline

Definition at line 78 of file G4AdjointSteppingAction.hh.

{start_event =aBool;}

void G4AdjointSteppingAction::SetUserAdjointSteppingAction ( G4UserSteppingAction *  anAction )

inline

Definition at line 86 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::SetAdjointSteppingAction().

                                                                              {
                                       theUserAdjointSteppingAction = anAction;}

void G4AdjointSteppingAction::SetUserForwardSteppingAction ( G4UserSteppingAction *  anAction )

inline

Definition at line 88 of file G4AdjointSteppingAction.hh.

                                                                              {
                                          theUserFwdSteppingAction = anAction;}

void G4AdjointSteppingAction::UserSteppingAction ( const G4Step *  aStep )

virtual

Reimplemented from G4UserSteppingAction.

Definition at line 59 of file G4AdjointSteppingAction.cc.

References G4AdjointCrossSurfChecker::CrossingOneOfTheRegisteredSurface(), fStopAndKill, fWorldBoundary, G4ParticleDefinition::GetBaryonNumber(), G4Track::GetDefinition(), G4Track::GetKineticEnergy(), G4Track::GetMomentum(), G4ParticleDefinition::GetParticleType(), G4Step::GetPostStepPoint(), G4StepPoint::GetStepStatus(), G4Step::GetTrack(), G4Track::GetWeight(), G4Track::SetTrackStatus(), and G4UserSteppingAction::UserSteppingAction().

{
  G4Track* aTrack =aStep->GetTrack();
  //forward tracking mode
   if(!is_adjoint_tracking_mode){
     //check if last adjoint did reach the external source
     if (!did_adj_part_reach_ext_source) {
        aTrack->SetTrackStatus(fStopAndKill);
        return;
     }
     if(theUserFwdSteppingAction)
       theUserFwdSteppingAction->UserSteppingAction(aStep);
     return;
   }


  //Apply first the user adjoint stepping action
  //---------------------------
  if (theUserAdjointSteppingAction) theUserAdjointSteppingAction->UserSteppingAction(aStep);


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