130 for(
size_t iReg = 0; iReg < numRegions; ++iReg) {
134 for(
size_t jMat = 0; jMat < numOfMat; ++jMat) {
182 if(0 > coupleIndex) {
return 0.0; }
201 if(0 > coupleIndex)
return 0.0;
205 if(tmax <= cutEnergy)
return 0.0;
226 if(0 > coupleIndex) {
return; }
234 if( maxEnergy < tmax) tmax = maxEnergy;
235 if( tmin >= tmax)
return;
240 G4double totalMomentum = sqrt(kineticEnergy*(totalEnergy + fMass));
250 if( deltaTkin <= 0. && fVerbose > 0)
252 G4cout<<
"G4PAIPhotModel::SampleSecondary e- deltaTkin = "<<deltaTkin<<
G4endl;
254 if( deltaTkin <= 0.)
return;
256 if( deltaTkin > tmax) deltaTkin = tmax;
268 kineticEnergy -= deltaTkin;
270 if( kineticEnergy <= 0. )
279 G4ThreeVector dir = totalMomentum*direction - deltaRay->GetMomentum();
280 direction = dir.unit();
283 vdp->push_back(deltaRay);
292 if( deltaTkin <= 0. )
294 G4cout<<
"G4PAIPhotonModel::SampleSecondary gamma deltaTkin = "<<deltaTkin<<
G4endl;
296 if( deltaTkin <= 0.)
return;
298 if( deltaTkin >= kineticEnergy )
300 deltaTkin = kineticEnergy;
304 G4double sintheta = sqrt((1.+costheta)*(1.-costheta));
308 G4double dirx = sintheta*cos(phi), diry = sintheta*sin(phi), dirz = costheta;
311 deltaDirection.rotateUz(direction);
313 if( kineticEnergy > 0.)
315 kineticEnergy -= deltaTkin;
330 vdp->push_back(photonRay);
344 if(0 > coupleIndex) {
return eloss; }
386 G4double etot = kineticEnergy + particleMass;
387 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*particleMass)/(etot*etot);
388 G4double siga = (1.0/beta2 - 0.5) * twopi_mc2_rcl2 * tmax * step
389 * electronDensity * q * q;
417 tmax = 2.0*electron_mass_c2*(gamma*gamma - 1.) /
418 (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