G4AdjointSteppingAction Class Reference

#include <G4AdjointSteppingAction.hh>

Inheritance diagram for G4AdjointSteppingAction:

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

Private Attributes
G4double	ext_sourceEMax
G4AdjointCrossSurfChecker *	theG4AdjointCrossSurfChecker
G4bool	start_event
G4bool	did_adj_part_reach_ext_source
G4ThreeVector	last_momentum
G4ThreeVector	last_pos
G4double	last_ekin
G4double	last_weight
G4double	prim_weight
G4ParticleDefinition *	last_part_def
G4UserSteppingAction *	theUserAdjointSteppingAction
G4UserSteppingAction *	theUserFwdSteppingAction
G4bool	is_adjoint_tracking_mode
G4bool	is_adjoint_geantino_tracking_mode

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

Here is the call graph for this function:

~G4AdjointSteppingAction()

G4AdjointSteppingAction::~G4AdjointSteppingAction

(

)

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

Here is the caller graph for this function:

GetLastEkin()

G4double G4AdjointSteppingAction::GetLastEkin ( )

inline

Definition at line 82 of file G4AdjointSteppingAction.hh.

{return last_ekin;}

Here is the caller graph for this function:

GetLastMomentum()

G4ThreeVector G4AdjointSteppingAction::GetLastMomentum ( )

inline

Definition at line 80 of file G4AdjointSteppingAction.hh.

{return last_momentum;}

Here is the caller graph for this function:

GetLastPartDef()

G4ParticleDefinition* G4AdjointSteppingAction::GetLastPartDef ( )

inline

Definition at line 85 of file G4AdjointSteppingAction.hh.

{return last_part_def;}

Here is the caller graph for this function:

GetLastPosition()

G4ThreeVector G4AdjointSteppingAction::GetLastPosition ( )

inline

Definition at line 81 of file G4AdjointSteppingAction.hh.

{return last_pos;}

Here is the caller graph for this function:

GetLastWeight()

G4double G4AdjointSteppingAction::GetLastWeight ( )

inline

Definition at line 83 of file G4AdjointSteppingAction.hh.

{return last_weight;}

Here is the caller graph for this function:

SetAdjointGeantinoTrackingMode()

void G4AdjointSteppingAction::SetAdjointGeantinoTrackingMode ( G4bool  aBool )

inline

Definition at line 91 of file G4AdjointSteppingAction.hh.

{is_adjoint_geantino_tracking_mode =aBool;}

SetAdjointTrackingMode()

void G4AdjointSteppingAction::SetAdjointTrackingMode ( G4bool  aBool )

inline

Definition at line 90 of file G4AdjointSteppingAction.hh.

{is_adjoint_tracking_mode =aBool;}

Here is the caller graph for this function:

SetExtSourceEMax()

void G4AdjointSteppingAction::SetExtSourceEMax ( G4double  Emax )

inline

Definition at line 77 of file G4AdjointSteppingAction.hh.

{ext_sourceEMax=Emax;} 

Here is the caller graph for this function:

SetPrimWeight()

void G4AdjointSteppingAction::SetPrimWeight ( G4double  weight )

inline

Definition at line 84 of file G4AdjointSteppingAction.hh.

{prim_weight=weight;} 

Here is the caller graph for this function:

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

Here is the caller graph for this function:

SetUserForwardSteppingAction()

void G4AdjointSteppingAction::SetUserForwardSteppingAction ( G4UserSteppingAction *  anAction )

inline

Definition at line 88 of file G4AdjointSteppingAction.hh.

                                                                              {
                                          theUserFwdSteppingAction = anAction;}

Here is the caller graph for this function:

UserSteppingAction()

void G4AdjointSteppingAction::UserSteppingAction ( const G4Step *  aStep )

virtual

Reimplemented from G4UserSteppingAction.

Definition at line 59 of file G4AdjointSteppingAction.cc.

{
  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
  //G4cout<<aStep->GetPostStepPoint()->GetStepStatus()<<std::endl;
  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;

  }

}

Here is the call graph for this function:

Member Data Documentation

did_adj_part_reach_ext_source

G4bool G4AdjointSteppingAction::did_adj_part_reach_ext_source

private

Definition at line 98 of file G4AdjointSteppingAction.hh.

ext_sourceEMax

G4double G4AdjointSteppingAction::ext_sourceEMax

private

Definition at line 94 of file G4AdjointSteppingAction.hh.

is_adjoint_geantino_tracking_mode

G4bool G4AdjointSteppingAction::is_adjoint_geantino_tracking_mode

private

Definition at line 107 of file G4AdjointSteppingAction.hh.

is_adjoint_tracking_mode

G4bool G4AdjointSteppingAction::is_adjoint_tracking_mode

private

Definition at line 106 of file G4AdjointSteppingAction.hh.

last_ekin

G4double G4AdjointSteppingAction::last_ekin

private

Definition at line 100 of file G4AdjointSteppingAction.hh.

last_momentum

G4ThreeVector G4AdjointSteppingAction::last_momentum

private

Definition at line 99 of file G4AdjointSteppingAction.hh.

last_part_def

G4ParticleDefinition* G4AdjointSteppingAction::last_part_def

private

Definition at line 103 of file G4AdjointSteppingAction.hh.

last_pos

G4ThreeVector G4AdjointSteppingAction::last_pos

private

Definition at line 99 of file G4AdjointSteppingAction.hh.

last_weight

G4double G4AdjointSteppingAction::last_weight

private

Definition at line 101 of file G4AdjointSteppingAction.hh.

prim_weight

G4double G4AdjointSteppingAction::prim_weight

private

Definition at line 102 of file G4AdjointSteppingAction.hh.

start_event

G4bool G4AdjointSteppingAction::start_event

private

Definition at line 96 of file G4AdjointSteppingAction.hh.

theG4AdjointCrossSurfChecker

G4AdjointCrossSurfChecker* G4AdjointSteppingAction::theG4AdjointCrossSurfChecker

private

Definition at line 95 of file G4AdjointSteppingAction.hh.

theUserAdjointSteppingAction

G4UserSteppingAction* G4AdjointSteppingAction::theUserAdjointSteppingAction

private

Definition at line 104 of file G4AdjointSteppingAction.hh.

theUserFwdSteppingAction

G4UserSteppingAction* G4AdjointSteppingAction::theUserFwdSteppingAction

private

Definition at line 105 of file G4AdjointSteppingAction.hh.

---

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

• G4AdjointSteppingAction.hh
• G4AdjointSteppingAction.cc