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