G4SingleScatter Class Reference

#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 ( )

inline

Definition at line 197 of file G4ScreenedNuclearRecoil.hh.

{ }

virtual G4SingleScatter::~G4SingleScatter ( )

inlinevirtual

Definition at line 200 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 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:

• source/geant4.10.03.p02/examples/extended/electromagnetic/TestEm7/include/G4ScreenedNuclearRecoil.hh
• source/geant4.10.03.p02/examples/extended/electromagnetic/TestEm7/src/G4ScreenedNuclearRecoil.cc