G4SingleScatter Class Reference

#include <G4ScreenedNuclearRecoil.hh>

Inheritance 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 178 of file G4ScreenedNuclearRecoil.hh.

Constructor & Destructor Documentation

G4SingleScatter::G4SingleScatter ( )

inline

Definition at line 181 of file G4ScreenedNuclearRecoil.hh.

{ }

virtual G4SingleScatter::~G4SingleScatter ( )

inlinevirtual

Definition at line 184 of file G4ScreenedNuclearRecoil.hh.

{}

Member Function Documentation

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

virtual

Implements G4ScreenedCollisionStage.

Definition at line 661 of file G4ScreenedNuclearRecoil.cc.

References G4CoulombKinematicsInfo::a2, G4ParticleChange::AddSecondary(), python.hepunit::amu_c2, G4CoulombKinematicsInfo::cosZeta, G4CoulombKinematicsInfo::eRecoil, G4UniformRand, G4Track::GetDynamicParticle(), G4ParticleChange::GetEnergy(), G4ScreenedNuclearRecoil::GetKinematics(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentum(), G4DynamicParticle::GetMomentumDirection(), G4ScreenedNuclearRecoil::GetParticleChange(), G4ScreenedNuclearRecoil::GetValidCollision(), G4ParticleChange::ProposeEnergy(), G4ParticleChange::ProposeMomentumDirection(), G4CoulombKinematicsInfo::recoilIon, CLHEP::Hep3Vector::rotateUz(), G4VParticleChange::SetNumberOfSecondaries(), G4CoulombKinematicsInfo::sinZeta, and CLHEP::Hep3Vector::unit().

                                                                                                                      {
        
        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*CLHEP::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);   
        }
}

---

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

• G4ScreenedNuclearRecoil.hh
• G4ScreenedNuclearRecoil.cc