102 mass = proton_mass_c2;
146 if(
IsMaster() && mass < GeV && part->GetParticleName() !=
"GenericIon") {
178 for (
G4int j=0; j<nelm; ++j) {
180 if(iz < Z) { Z =
iz; }
204 if(kinEnergy <= 0.0) {
return cross; }
216 if(costmin > costmax) {
235 std::vector<G4DynamicParticle*>* fvect,
267 if(costmin > costmax) {
270 G4double ratio = ecross/(cross + ecross);
288 newDirection.rotateUz(direction);
296 /(targetMass + (
mass + kinEnergy)*(1.0 - cost));
299 if(trec > kinEnergy) { trec = kinEnergy; }
313 newDirection*sqrt(finalT*(2*
mass + finalT))).unit();
315 fvect->push_back(newdp);
327 if(edep < 0.0) { edep = 0.0; }
const G4ParticleDefinition * particle
static G4double GetNuclearMass(const G4double A, const G4double Z)
void SetTargetMass(G4double value)
std::vector< G4Element * > G4ElementVector
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
const std::vector< G4double > * pCuts
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
const G4MaterialCutsCouple * currentCouple
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
const G4Material * currentMaterial
G4ParticleDefinition * GetDefinition() const
G4double GetMomentumSquare() const
const G4ElementVector * GetElementVector() const
static G4NistManager * Instance()
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4double SetupTarget(G4int Z, G4double cut=DBL_MAX)
G4int currentMaterialIndex
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4IonTable * GetIonTable() const
const G4MaterialCutsCouple * CurrentCouple() const
double A(double temperature)
void ProposeNonIonizingEnergyDeposit(G4double anEnergyPart)
const G4ThreeVector & GetMomentumDirection() const
G4double SetupKinematic(G4double kinEnergy, const G4Material *mat)
static G4Proton * Proton()
void DefineMaterial(const G4MaterialCutsCouple *)
G4int SelectIsotopeNumber(const G4Element *)
G4eCoulombScatteringModel(G4bool combined=true)
std::vector< G4EmElementSelector * > * GetElementSelectors()
virtual G4double MinPrimaryEnergy(const G4Material *, const G4ParticleDefinition *, G4double)
G4double ComputeElectronCrossSection(G4double CosThetaMin, G4double CosThetaMax)
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
static G4ParticleTable * GetParticleTable()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double ComputeNuclearCrossSection(G4double CosThetaMin, G4double CosThetaMax)
G4double PolarAngleLimit() const
G4double GetAtomicMassAmu(const G4String &symb) const
const G4ParticleDefinition * theProton
G4NistManager * fNistManager
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
void SetProposedKineticEnergy(G4double proposedKinEnergy)
size_t GetNumberOfElements() const
G4WentzelOKandVIxSection * wokvi
void SetupParticle(const G4ParticleDefinition *)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax)
virtual ~G4eCoulombScatteringModel()
G4ParticleChangeForGamma * fParticleChange
const G4Element * SelectRandomAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
void Initialise(const G4ParticleDefinition *, G4double CosThetaLim)
G4ParticleChangeForGamma * GetParticleChangeForGamma()
G4ThreeVector & SampleSingleScattering(G4double CosThetaMin, G4double CosThetaMax, G4double elecRatio=0.0)