F04SteppingAction Class Reference

#include <F04SteppingAction.hh>

Inheritance diagram for F04SteppingAction:

Public Member Functions
	F04SteppingAction ()
virtual	~F04SteppingAction ()
virtual void	UserSteppingAction (const G4Step *)
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 43 of file F04SteppingAction.hh.

Constructor & Destructor Documentation

F04SteppingAction::F04SteppingAction

(

)

Definition at line 45 of file F04SteppingAction.cc.

 : G4UserSteppingAction(),
   fSteppingMessenger(0),
   fTargetVolume(0),
   fTestPlaneVolume(0)
{
  fSteppingMessenger = new F04SteppingActionMessenger();
}

F04SteppingAction::~F04SteppingAction ( )

virtual

Definition at line 56 of file F04SteppingAction.cc.

{
  delete fSteppingMessenger ;
}

Member Function Documentation

void F04SteppingAction::UserSteppingAction ( const G4Step *  theStep )

virtual

Reimplemented from G4UserSteppingAction.

Definition at line 63 of file F04SteppingAction.cc.

References fAlive, fStopAndKill, G4LogicalVolumeStore::GetInstance(), G4VPhysicalVolume::GetLogicalVolume(), G4Track::GetMomentumDirection(), G4Track::GetParentID(), G4StepPoint::GetPhysicalVolume(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4StepPoint::GetProcessDefinedStep(), G4VProcess::GetProcessName(), G4Step::GetTrack(), G4Track::GetTrackStatus(), F04UserTrackInformation::GetTrackStatusFlag(), G4Track::GetUserInformation(), G4Track::GetVertexMomentumDirection(), G4LogicalVolumeStore::GetVolume(), reverse, G4Track::SetTrackStatus(), F04UserTrackInformation::SetTrackStatusFlag(), and CLHEP::Hep3Vector::z().

{
  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 = 0;
  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();

  F04UserTrackInformation* 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() != 0) {
     if (thePostPoint->GetProcessDefinedStep()->GetProcessName() != "Decay")
                                                                       return;
  }

}

---

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

• F04SteppingAction.hh
• F04SteppingAction.cc