131 if(0 == numRegions) {
133 "no G4Regions are registered for the PAI model - World is used");
135 ->GetRegion(
"DefaultRegionForTheWorld",
false));
140 G4cout <<
"G4PAIModel is defined for " << numRegions <<
" regions "
142 G4cout <<
" total number of materials " << numOfMat <<
G4endl;
144 for(
size_t iReg = 0; iReg<numRegions; ++iReg) {
148 for(
size_t jMat = 0; jMat<numOfMat; ++jMat) {
162 << mat->
GetName() <<
"> CoupleIndex= "
165 <<
" cutsize= " << cuts.size() <<
G4endl;
170 for(
size_t i=0; i<
n; ++i) {
223 if(0 > coupleIndex) {
return 0.0; }
247 if(0 > coupleIndex) {
return 0.0; }
250 if(tmax <= cutEnergy) {
return 0.0; }
274 if(0 > coupleIndex) {
return; }
280 if(maxEnergy < tmax) { tmax = maxEnergy; }
281 if(tmin >= tmax) {
return; }
286 G4double totalMomentum = sqrt(kineticEnergy*(totalEnergy+fMass));
294 if( !(deltaTkin <= 0.) && !(deltaTkin > 0)) {
295 G4cout<<
"G4PAIModel::SampleSecondaries; deltaKIn = "<<deltaTkin/
keV
296 <<
" keV "<<
" Escaled(MeV)= " << scaledTkin <<
G4endl;
299 if( deltaTkin <= 0.) {
return; }
301 if( deltaTkin > tmax) { deltaTkin = tmax; }
314 kineticEnergy -= deltaTkin;
316 direction = dir.
unit();
320 vdp->push_back(deltaRay);
331 if(0 > coupleIndex) {
return eloss; }
369 G4double etot = kineticEnergy + particleMass;
370 G4double beta2 = kineticEnergy*(kineticEnergy + 2.0*particleMass)/(etot*etot);
372 * electronDensity * q * q;
389 (1. + 2.0*gamma*ratio + ratio*
ratio);