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)

---

Detailed Description

Definition at line 69 of file G4AdjointSteppingAction.hh.

---

Constructor & Destructor Documentation

G4AdjointSteppingAction::G4AdjointSteppingAction

(

)

Definition at line 43 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)
{ 
  theG4AdjointCrossSurfChecker = G4AdjointCrossSurfChecker::GetInstance();
}

G4AdjointSteppingAction::~G4AdjointSteppingAction

(

)

Definition at line 52 of file G4AdjointSteppingAction.cc.

{;}

---

Member Function Documentation

G4bool G4AdjointSteppingAction::GetDidAdjParticleReachTheExtSource

(

)

[inline]

Definition at line 79 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::GetDidAdjParticleReachTheExtSource().

{return did_adj_part_reach_ext_source;}

G4double G4AdjointSteppingAction::GetLastEkin

(

)

[inline]

Definition at line 82 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_ekin;}

G4ThreeVector G4AdjointSteppingAction::GetLastMomentum

(

)

[inline]

Definition at line 80 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_momentum;}

G4ParticleDefinition* G4AdjointSteppingAction::GetLastPartDef

(

)

[inline]

Definition at line 85 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_part_def;}

G4ThreeVector G4AdjointSteppingAction::GetLastPosition

(

)

[inline]

Definition at line 81 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_pos;}

G4double G4AdjointSteppingAction::GetLastWeight

(

)

[inline]

Definition at line 83 of file G4AdjointSteppingAction.hh.

Referenced by G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

{return last_weight;}

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 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::UserSteppingAction

(

const G4Step *

)

[virtual]

Reimplemented from G4UserSteppingAction.

Definition at line 57 of file G4AdjointSteppingAction.cc.

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

{
  //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

---