#include <Doxymodules_field.h>

Inheritance diagram for F04SteppingAction:

Public Member Functions
	F04SteppingAction ()=default

	~F04SteppingAction () override=default

void	UserSteppingAction (const G4Step *) override

Private Attributes
G4LogicalVolume *	fTargetVolume = nullptr

G4LogicalVolume *	fTestPlaneVolume = nullptr

Detailed Description

Definition at line 119 of file Doxymodules_field.h.

Constructor & Destructor Documentation

◆ F04SteppingAction()

F04SteppingAction::F04SteppingAction ( )

default

◆ ~F04SteppingAction()

F04SteppingAction::~F04SteppingAction ( )

overridedefault

Member Function Documentation

◆ UserSteppingAction()

void F04SteppingAction::UserSteppingAction ( const G4Step * theStep )

override

Definition at line 43 of file F04SteppingAction.cc.

{
  G4Track* theTrack = theStep->GetTrack();
 
  // Get pointers to test volumes (only once)
  if ( ! fTargetVolume ) {
    fTargetVolume
      = G4LogicalVolumeStore::GetInstance()->GetVolume("Target");
    fTestPlaneVolume
     = G4LogicalVolumeStore::GetInstance()->GetVolume("TestPlane");
  }
 
  if (theTrack->GetParentID()==0) {
    //This is a primary track
    G4LogicalVolume* theVolume
      = theStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
    if (theVolume != fTargetVolume) {
       theTrack->SetTrackStatus(fStopAndKill);
       return;
    }
  }
 
//  G4ParticleDefinition* particleType = theTrack->GetDefinition();
 
  // check if it is entering the test volume
  G4StepPoint* thePrePoint = theStep->GetPreStepPoint();
  G4VPhysicalVolume* thePrePV = thePrePoint->GetPhysicalVolume();
  G4LogicalVolume* thePreLV = thePrePV->GetLogicalVolume();
 
  G4LogicalVolume* thePostLV = nullptr;
  G4StepPoint* thePostPoint = theStep->GetPostStepPoint();
 
  if (thePostPoint) {
     G4VPhysicalVolume* thePostPV = thePostPoint->GetPhysicalVolume();
     if (thePostPV) thePostLV = thePostPV->GetLogicalVolume();
  }
 
  if (thePostLV == fTestPlaneVolume &&
      thePreLV  != fTestPlaneVolume) {
 
//     G4double x = theTrack->GetPosition().x();
//     G4double y = theTrack->GetPosition().y();
 
     G4ThreeVector theMomentumDirection = theTrack->
                          GetDynamicParticle()->GetMomentumDirection();
//     G4double theTotalMomentum = theTrack->GetDynamicParticle()->
//                                                    GetTotalMomentum();
 
     // then kill the track
     theTrack->SetTrackStatus(fStopAndKill);
     return;
 
  }
 
//  G4double z = theTrack->GetPosition().z();
 
  auto  trackInformation =
                      (F04UserTrackInformation*)theTrack->GetUserInformation();
 
  if (trackInformation->GetTrackStatusFlag() != reverse) {
     if (thePreLV  != fTargetVolume ) {
        if ( theTrack->      GetMomentumDirection().z()>0.0 &&
             theTrack->GetVertexMomentumDirection().z()<0.0     )
        {
             trackInformation->SetTrackStatusFlag(reverse);
        }
     }
  }
 
  // check if it is alive
  if (theTrack->GetTrackStatus() == fAlive) { return; }
 
  if (thePostPoint->GetProcessDefinedStep() != nullptr) {
     if (thePostPoint->GetProcessDefinedStep()->GetProcessName() != "Decay")
                                                                       return;
  }
 
}

Member Data Documentation

◆ fTargetVolume

G4LogicalVolume* F04SteppingAction::fTargetVolume = nullptr

private

Definition at line 51 of file F04SteppingAction.hh.

◆ fTestPlaneVolume

G4LogicalVolume* F04SteppingAction::fTestPlaneVolume = nullptr