27 #include "CLHEP/Units/PhysicalConstants.h"
90 G4cout <<
" G4LEPTSIonisationModel::SampleSecondaries( SetProposedKineticEnergy " << P1KinEn <<
" " << P0KinEn <<
" - " << Energylost <<
G4endl;
110 P3Dir /= P3Dir.mag();
113 fvect->push_back(e3);
119 G4double cp0 = std::sqrt(P0KinEn*(P0KinEn + 2.*CLHEP::electron_mass_c2));
120 G4double cp1 = std::sqrt(P1KinEn*(P1KinEn + 2.*CLHEP::electron_mass_c2));
123 P2Dir.rotateUz(P0Dir);
125 fvect->push_back(e2);
G4LEPTSIonisationModel(const G4String &modelName="G4LEPTSIonisationModel")
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
std::map< const G4Material *, G4double > theIonisPot
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4double GetMeanFreePath(const G4Material *mate, const G4ParticleDefinition *aParticle, G4double kineticEnergy)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double tmax=DBL_MAX)
G4GLOB_DLL std::ostream G4cout
G4double SampleEnergyLoss(const G4Material *aMaterial, G4double eMin, G4double eMax)
const G4ThreeVector & GetMomentumDirection() const
void BuildPhysicsTable(const G4ParticleDefinition &aParticleType)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
~G4LEPTSIonisationModel()
static G4Electron * Electron()
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4ThreeVector SampleNewDirection(const G4Material *aMaterial, G4ThreeVector Dir, G4double e, G4double el)
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4double theLowestEnergyLimit
void SetDeexcitationFlag(G4bool val)
G4ParticleChangeForGamma * fParticleChangeForGamma
G4ParticleChangeForGamma * GetParticleChangeForGamma()
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
const G4Material * GetMaterial() const
std::map< const G4Material *, G4double > theIonisPotInt