131 for(
size_t iReg = 0; iReg < numRegions; ++iReg) {
135 for(
size_t jMat = 0; jMat < numOfMat; ++jMat) {
183 if(0 > coupleIndex) {
return 0.0; }
202 if(0 > coupleIndex)
return 0.0;
206 if(tmax <= cutEnergy)
return 0.0;
227 if(0 > coupleIndex) {
return; }
235 if( maxEnergy < tmax) tmax = maxEnergy;
236 if( tmin >= tmax)
return;
241 G4double totalMomentum = sqrt(kineticEnergy*(totalEnergy + fMass));
251 if( deltaTkin <= 0. && fVerbose > 0)
253 G4cout<<
"G4PAIPhotModel::SampleSecondary e- deltaTkin = "<<deltaTkin<<
G4endl;
255 if( deltaTkin <= 0.)
return;
257 if( deltaTkin > tmax) deltaTkin = tmax;
269 kineticEnergy -= deltaTkin;
271 if( kineticEnergy <= 0. )
280 G4ThreeVector dir = totalMomentum*direction - deltaRay->GetMomentum();
281 direction = dir.unit();
284 vdp->push_back(deltaRay);
293 if( deltaTkin <= 0. )
295 G4cout<<
"G4PAIPhotonModel::SampleSecondary gamma deltaTkin = "<<deltaTkin<<
G4endl;
297 if( deltaTkin <= 0.)
return;
299 if( deltaTkin >= kineticEnergy )
301 deltaTkin = kineticEnergy;
305 G4double sintheta = sqrt((1.+costheta)*(1.-costheta));
309 G4double dirx = sintheta*cos(phi), diry = sintheta*sin(phi), dirz = costheta;
312 deltaDirection.rotateUz(direction);
314 if( kineticEnergy > 0.)
316 kineticEnergy -= deltaTkin;
331 vdp->push_back(photonRay);
345 if(0 > coupleIndex) {
return eloss; }
387 G4double etot = kineticEnergy + particleMass;
388 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*particleMass)/(etot*etot);
389 G4double siga = (1.0/beta2 - 0.5) * twopi_mc2_rcl2 * tmax * step
390 * electronDensity * q * q;
418 tmax = 2.0*electron_mass_c2*(gamma*gamma - 1.) /
419 (1. + 2.0*gamma*ratio + ratio*ratio);
const G4ParticleDefinition * fParticle
G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kinEnergy)
G4double SamplePostStepPhotonTransfer(G4int coupleIndex, G4double scaledTkin) const
virtual ~G4PAIPhotModel()
G4double LowEnergyLimit() const
G4ParticleChangeForLoss * GetParticleChangeForLoss()
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
G4double HighEnergyLimit() const
G4PAIPhotData * fModelData
G4PAIPhotData * GetPAIPhotData()
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
G4double GetPlasmonRatio(G4int coupleIndex, G4double scaledTkin) const
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
static G4MaterialTable * GetMaterialTable()
G4double SampleAlongStepPhotonTransfer(G4int coupleIndex, G4double kinEnergy, G4double scaledTkin, G4double stepFactor) const
std::vector< G4Material * > G4MaterialTable
G4VEmAngularDistribution * GetAngularDistribution()
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy)
std::vector< const G4MaterialCutsCouple * > fMaterialCutsCoupleVector
void SetMomentumDirection(const G4ThreeVector &aDirection)
G4ParticleDefinition * GetDefinition() const
const G4String & GetParticleName() const
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
void Initialise(const G4MaterialCutsCouple *, G4double cut, G4PAIPhotModel *)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
static const G4double reg
std::vector< const G4Region * > fPAIRegionVector
const G4MaterialCutsCouple * CurrentCouple() const
G4GLOB_DLL std::ostream G4cout
G4int FindCoupleIndex(const G4MaterialCutsCouple *)
G4double GetElectronDensity() const
const G4ParticleDefinition * fPositron
const G4ThreeVector & GetMomentumDirection() const
G4double DEDXPerVolume(G4int coupleIndex, G4double scaledTkin, G4double cut) const
G4double GetCharge() const
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4MaterialCutsCouple * FindCouple(G4Material *mat)
std::vector< G4EmElementSelector * > * GetElementSelectors()
static size_t GetNumberOfMaterials()
G4ParticleChangeForLoss * fParticleChange
void SetProposedMomentumDirection(const G4ThreeVector &dir)
void SetKineticEnergy(G4double aEnergy)
G4PAIPhotModel(const G4ParticleDefinition *p=0, const G4String &nam="PAI")
static G4Positron * Positron()
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
virtual G4double SampleFluctuations(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double, G4double, G4double)
G4double SamplePostStepPlasmonTransfer(G4int coupleIndex, G4double scaledTkin) const
virtual G4double Dispersion(const G4Material *, const G4DynamicParticle *, G4double, G4double)
void SetAngularDistribution(G4VEmAngularDistribution *)
G4double CrossSectionPerVolume(G4int coupleIndex, G4double scaledTkin, G4double tcut, G4double tmax) const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double SampleAlongStepPlasmonTransfer(G4int coupleIndex, G4double kinEnergy, G4double scaledTkin, G4double stepFactor) const
const G4ParticleDefinition * fElectron
static G4Electron * Electron()
void DefineForRegion(const G4Region *r)
void SetParticle(const G4ParticleDefinition *p)
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
static const double eplus
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
const G4Element * SelectRandomAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
const G4Material * GetMaterial() const