#include <G4ScreenedNuclearRecoil.hh>

Inheritance diagram for G4SingleScatter:

Collaboration diagram for G4SingleScatter:

Public Member Functions
	G4SingleScatter ()

virtual void	DoCollisionStep (class G4ScreenedNuclearRecoil *master, const class G4Track &aTrack, const class G4Step &aStep)

virtual	~G4SingleScatter ()

Public Member Functions inherited from G4ScreenedCoulombCrossSectionInfo
	G4ScreenedCoulombCrossSectionInfo ()

	~G4ScreenedCoulombCrossSectionInfo ()

Public Member Functions inherited from G4ScreenedCollisionStage
virtual	~G4ScreenedCollisionStage ()

Additional Inherited Members
Static Public Member Functions inherited from G4ScreenedCoulombCrossSectionInfo
static const char *	CVSHeaderVers ()

static const char *	CVSFileVers ()

Detailed Description

Definition at line 193 of file G4ScreenedNuclearRecoil.hh.

Constructor & Destructor Documentation

◆ G4SingleScatter()

G4SingleScatter::G4SingleScatter ( )

inline

Definition at line 197 of file G4ScreenedNuclearRecoil.hh.

197 { }

◆ ~G4SingleScatter()

virtual G4SingleScatter::~G4SingleScatter ( )

inlinevirtual

Definition at line 200 of file G4ScreenedNuclearRecoil.hh.

200 {}

Member Function Documentation

◆ DoCollisionStep()

void G4SingleScatter::DoCollisionStep	(	class G4ScreenedNuclearRecoil *	master,
		const class G4Track &	aTrack,
		const class G4Step &	aStep
	)

virtual

Implements G4ScreenedCollisionStage.

Definition at line 752 of file G4ScreenedNuclearRecoil.cc.

                                                                             {
         
         if(!master->GetValidCollision()) return;
 
         G4CoulombKinematicsInfo &kin=master->GetKinematics();
         G4ParticleChange &aParticleChange=master->GetParticleChange();
 
         const G4DynamicParticle* incidentParticle = aTrack.GetDynamicParticle();
         G4double incidentEnergy = incidentParticle->GetKineticEnergy();
         G4double eRecoil=kin.eRecoil;
         
         G4double azimuth=G4UniformRand()*(2.0*pi);
         G4double sa=std::sin(azimuth);
         G4double ca=std::cos(azimuth);
 
         G4ThreeVector recoilMomentumDirection(kin.sinZeta*ca, 
                                               kin.sinZeta*sa, kin.cosZeta);
         G4ParticleMomentum incidentDirection = 
           incidentParticle->GetMomentumDirection();
         recoilMomentumDirection=
           recoilMomentumDirection.rotateUz(incidentDirection);
         G4ThreeVector recoilMomentum=
           recoilMomentumDirection*std::sqrt(2.0*eRecoil*kin.a2*amu_c2);
         
         if(aParticleChange.GetEnergy() != 0.0) { 
           // DoKinematics hasn't stopped it!
           G4ThreeVector beamMomentum=
             incidentParticle->GetMomentum()-recoilMomentum;
           aParticleChange.ProposeMomentumDirection(beamMomentum.unit()) ;
           aParticleChange.ProposeEnergy(incidentEnergy-eRecoil);
         }
         
         if(kin.recoilIon) {
                 G4DynamicParticle* recoil = 
                   new G4DynamicParticle (kin.recoilIon,
                                 recoilMomentumDirection,eRecoil) ;
                 
                 aParticleChange.SetNumberOfSecondaries(1);
                 aParticleChange.AddSecondary(recoil);   
         }
 }

Here is the call graph for this function:

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