#include <XPhononStackingAction.hh>

Inheritance diagram for XPhononStackingAction:

Collaboration diagram for XPhononStackingAction:

Public Member Functions
	XPhononStackingAction ()

virtual	~XPhononStackingAction ()

virtual G4ClassificationOfNewTrack	ClassifyNewTrack (const G4Track *aTrack)

Public Member Functions inherited from G4UserStackingAction
	G4UserStackingAction ()

virtual	~G4UserStackingAction ()

void	SetStackManager (G4StackManager *value)

virtual void	NewStage ()

virtual void	PrepareNewEvent ()

Additional Inherited Members
Protected Attributes inherited from G4UserStackingAction
G4StackManager *	stackManager

Detailed Description

Definition at line 39 of file XPhononStackingAction.hh.

Constructor & Destructor Documentation

◆ XPhononStackingAction()

XPhononStackingAction::XPhononStackingAction ( )

Definition at line 50 of file XPhononStackingAction.cc.

50 {;}

◆ ~XPhononStackingAction()

XPhononStackingAction::~XPhononStackingAction ( )

virtual

Definition at line 54 of file XPhononStackingAction.cc.

54 {;}

Member Function Documentation

◆ ClassifyNewTrack()

G4ClassificationOfNewTrack XPhononStackingAction::ClassifyNewTrack ( const G4Track * aTrack )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 59 of file XPhononStackingAction.cc.

                                                              {
   G4ClassificationOfNewTrack classification = fUrgent;
 
   if (aTrack->GetParentID() == 0) {
     //Obtain LatticeManager for phonon dynamics
     G4LatticeManager* LM = G4LatticeManager::GetLatticeManager();
     
     G4int pol = G4PhononPolarization::Get(aTrack->GetDefinition());
     
     //Compute random wave-vector (override whatever ParticleGun did)
     G4ThreeVector Ran = G4RandomDirection();
     
     //Store wave-vector as track information
     G4PhononTrackMap* theKmap = G4PhononTrackMap::GetPhononTrackMap();
     theKmap->SetK(aTrack, Ran);
     
     //Compute direction of propagation from wave vector
     G4ThreeVector momentumDir = LM->MapKtoVDir(aTrack->GetVolume(), pol, Ran);
     
     //Compute true velocity of propagation
     G4double velocity = LM->MapKtoV(aTrack->GetVolume(), pol, Ran);
     
     //cast to non-const pointer so we can set the velocity
     G4Track* theTrack = const_cast<G4Track*>(aTrack);
 
     theTrack->SetMomentumDirection(momentumDir);
     theTrack->SetVelocity(velocity);
     theTrack->UseGivenVelocity(true);
   }
   
   return classification; 
 }

Here is the call graph for this function:

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