G4NuclearRadii Class Reference

#include <G4NuclearRadii.hh>

Static Public Member Functions
static G4double	ExplicitRadius (G4int Z, G4int A)
static G4double	Radius (G4int Z, G4int A)
static G4double	RadiusRMS (G4int Z, G4int A)
static G4double	RadiusNNGG (G4int Z, G4int A)
static G4double	RadiusECS (G4int Z, G4int A)
static G4double	RadiusHNGG (G4int A)
static G4double	RadiusKNGG (G4int A)
static G4double	RadiusND (G4int A)
static G4double	RadiusCB (G4int Z, G4int A)
static G4double	ParticleRadius (const G4ParticleDefinition *)
static G4double	CoulombFactor (const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)
static G4double	CoulombFactor (G4int Z, G4int A, const G4ParticleDefinition *theParticle, G4double ekin)
static G4double	NeutronInelasticShape (G4int Z, G4double ekin)
static G4double	ProtonInelasticShape (G4int Z, G4double ekin)

Static Public Attributes
static G4Pow *	fG4pow = G4Pow::GetInstance()
static const G4double	r0 [93]

Detailed Description

Definition at line 43 of file G4NuclearRadii.hh.

Member Function Documentation

CoulombFactor() [1/2]

G4double G4NuclearRadii::CoulombFactor ( const G4ParticleDefinition * theParticle, const G4ParticleDefinition * nucleon, G4double ekin )

static

Definition at line 180 of file G4NuclearRadii.cc.

{
  G4double tR = 0.895*CLHEP::fermi;
  G4double pR = ParticleRadius(theParticle);
 
  G4double pZ = theParticle->GetPDGCharge()*fInvep;
  G4double tZ = nucleon->GetPDGCharge()*fInvep;
 
  G4double pM = theParticle->GetPDGMass(); 
  G4double tM = nucleon->GetPDGMass();
 
  G4double pElab = ekin + pM;
  G4double totTcm  = std::sqrt(pM*pM + tM*tM + 2.*pElab*tM) - pM -tM;
  
  G4double bC = fAlpha*pZ*tZ/(pR + tR);
  return (totTcm > bC) ? 1. - bC/totTcm : 0.0;
}

Referenced by G4HadronNucleonXsc::HadronNucleonXscNS(), and G4HadronNucleonXsc::KaonNucleonXscVG().

CoulombFactor() [2/2]

G4double G4NuclearRadii::CoulombFactor ( G4int Z, G4int A, const G4ParticleDefinition * theParticle, G4double ekin )

static

Definition at line 201 of file G4NuclearRadii.cc.

{
  G4double tR = RadiusCB(Z, A);
  G4double pR = ParticleRadius(theParticle);
 
  G4double pZ = theParticle->GetPDGCharge()*fInvep;
 
  G4double pM = theParticle->GetPDGMass(); 
  G4double tM = G4NucleiProperties::GetNuclearMass(A, Z);
 
  G4double pElab = ekin + pM;
  G4double totTcm  = std::sqrt(pM*pM + tM*tM + 2.*pElab*tM) - pM -tM;
  
  G4double bC = fAlpha*pZ*Z/(pR + tR);
  return (totTcm > bC) ? 1. - bC/totTcm : 0.0;
}

ExplicitRadius()

G4double G4NuclearRadii::ExplicitRadius ( G4int Z, G4int A )

static

Definition at line 51 of file G4NuclearRadii.cc.

{
  G4double R = 0.0;
  // Special rms radii for light nucleii
  if(Z <= 4) {
    if(A == 1)                { R = 0.895*CLHEP::fermi; }// p
    else if(A == 2)           { R = 2.13*CLHEP::fermi; }// d
    else if(Z == 1 && A == 3) { R = 1.80*CLHEP::fermi; }// t
    else if(Z == 2 && A == 3) { R = 1.96*CLHEP::fermi; }// He3
    else if(Z == 2 && A == 4) { R = 1.68*CLHEP::fermi; }// He4
    else if(Z == 3)           { R = 2.40*CLHEP::fermi; }// Li7
    else if(Z == 4)           { R = 2.51*CLHEP::fermi; }// Be9
  }
  return R;
}

Referenced by Radius(), RadiusCB(), RadiusNNGG(), and RadiusRMS().

NeutronInelasticShape()

G4double G4NuclearRadii::NeutronInelasticShape ( G4int Z, G4double ekin )

static

Definition at line 222 of file G4NuclearRadii.cc.

{
  G4double A = (Z < 100) ? aeff[Z] : aeff[100];
  G4double elog = G4Log(ekin/CLHEP::GeV)/llog10;
  G4double p3 = 0.6 + 13./A - 0.0005*A;
  G4double p4 = 7.2449 - 0.018242*A;
  G4double p5 = 1.36 + 1.8/A + 0.0005*A;
  G4double p6 = 1. + 200./A + 0.02*A;
  G4double p7 = 3.0 - (A - 70.)*(A - 200.)/11000.;
 
  G4double firstexp = G4Exp(-p4*(elog + p5));
  G4double secondexp = G4Exp(-p6*(elog + p7));
 
  return (1. + p3*firstexp/(1. + firstexp))/(1. + secondexp);
}

ParticleRadius()

G4double G4NuclearRadii::ParticleRadius ( const G4ParticleDefinition * p )

static

Definition at line 169 of file G4NuclearRadii.cc.

{
  G4double R = CLHEP::fermi;
  G4int pdg = std::abs(p->GetPDGEncoding());
  if(pdg == 2112 || pdg == 2212)   { R *= 0.895; }
  else if(pdg == 211)  { R *= 0.663; }
  else if(pdg == 321)  { R *= 0.340; }
  else { R *= 0.5; }
  return R;
}

Referenced by CoulombFactor(), and CoulombFactor().

ProtonInelasticShape()

G4double G4NuclearRadii::ProtonInelasticShape ( G4int Z, G4double ekin )

static

Definition at line 239 of file G4NuclearRadii.cc.

{
  G4double A = (Z < 100) ? aeff[Z] : aeff[100];
  G4double elog = G4Log(ekin/CLHEP::GeV)/llog10;
  G4double ff1 = 5.6  - 0.016*A; // slope of the drop at medium energies.
  G4double ff2 = 1.37 + 1.37/A;  // start of the slope.
  G4double ff3 = 0.8  + 18./A - 0.002*A;   // stephight
  G4double res = (1.0 + ff3*(1.0 - (1.0/(1+G4Exp(-ff1*(elog + ff2))))));
  ff1 = 8. - 8./A - 0.008*A; // slope of the rise
  ff2 = 2.34 - 5.4/A - 0.0028*A; // start of the rise
  res /= (1.0 + G4Exp(-ff1*(elog + ff2)));
  return res;
}

Radius()

G4double G4NuclearRadii::Radius ( G4int Z, G4int A )

static

Definition at line 67 of file G4NuclearRadii.cc.

{
  G4double R = ExplicitRadius(Z, A);
  if(0.0 == R) {
    if (A <= 50) {
      G4double y = 1.1;
      if( A <= 15)      { y = 1.26; }
      else if( A <= 20) { y = 1.19; }
      else if( A <= 30) { y = 1.12; }
      G4double x = fG4pow->Z13(A);
      R = y*(x - 1./x);
    } else {
      R = fG4pow->powZ(A, 0.27);
    }
    R *= CLHEP::fermi;
  }
  return R;
}

RadiusCB()

G4double G4NuclearRadii::RadiusCB ( G4int Z, G4int A )

static

Definition at line 159 of file G4NuclearRadii.cc.

{
  G4double R = ExplicitRadius(Z, A);
  if(0.0 == R) {
    G4int z = std::min(Z, 92);
    R = r0[z]*fG4pow->Z13(A)*CLHEP::fermi;
  }
  return R;
}

Referenced by G4FermiBreakUpUtil::CoulombBarrier(), CoulombFactor(), G4CoulombBarrier::G4CoulombBarrier(), and G4CoulombBarrier::GetCoulombBarrier().

RadiusECS()

G4double G4NuclearRadii::RadiusECS ( G4int Z, G4int A )

static

Definition at line 109 of file G4NuclearRadii.cc.

{
  G4double R=0.;
  const G4double c[3]={0.77329745, 1.38206072, 30.28295235};
  const G4double c1=c[0];
  const G4double c2=c[1];
  const G4double c3=c[2];
 
  // Special rms radii for light nuclei
  if (A <= 30) {
    G4double vn = 0.5*A + fG4pow->powN(0.028*A,2) - fG4pow->powN(0.011*A,3);
    G4double dev = vn - (A-Z);
    R = c1*fG4pow->Z13(A) + c2/fG4pow->Z13(A) + c3*dev*dev/(A*A);
  } else if (A<=50){
    G4double y = 1.1; 
    G4double x = fG4pow->Z13(A);
    R = y*(x - 1./x);
  }
  return R*CLHEP::fermi;
}

RadiusHNGG()

G4double G4NuclearRadii::RadiusHNGG ( G4int A )

static

Definition at line 130 of file G4NuclearRadii.cc.

{
  G4double R = CLHEP::fermi;
  if(A > 20) {
    R *= 1.08*fG4pow->Z13(A)*(0.8 + 0.2*G4Exp(-(G4double)(A - 20)/20.));
  } else {
    R *= 1.08*fG4pow->Z13(A)*(1.0 + 0.1*G4Exp(-(G4double)(A - 20)/20.));
  }
  return R;
}

Referenced by G4ComponentGGHadronNucleusXsc::ComputeCrossSections().

RadiusKNGG()

G4double G4NuclearRadii::RadiusKNGG ( G4int A )

static

Definition at line 141 of file G4NuclearRadii.cc.

{
  return 1.3*CLHEP::fermi*fG4pow->Z13(A);
}

Referenced by G4ComponentGGHadronNucleusXsc::ComputeCrossSections().

RadiusND()

G4double G4NuclearRadii::RadiusND ( G4int A )

static

Definition at line 146 of file G4NuclearRadii.cc.

{
  G4double R = CLHEP::fermi;
  if (1 == A) {
    R *= 0.895;
  } else {
    R *= fG4pow->Z13(A);
    if (A <= 3.) { R *= 0.8; }
    else { R *= 1.7; }
  }
  return R;
}

RadiusNNGG()

G4double G4NuclearRadii::RadiusNNGG ( G4int Z, G4int A )

static

Definition at line 95 of file G4NuclearRadii.cc.

{
  G4double R = ExplicitRadius(Z, A);
  if(0.0 == R) {
    if(A > 20) {
      R = 1.08*fG4pow->Z13(A)*(0.85 + 0.15*G4Exp(-(G4double)(A - 21)/40.));
    } else {
      R = 1.08*fG4pow->Z13(A)*(1.0 + 0.3*G4Exp(-(G4double)(A - 21)/10.));
    }
    R *= CLHEP::fermi;
  }
  return R;
}

RadiusRMS()

G4double G4NuclearRadii::RadiusRMS ( G4int Z, G4int A )

static

Definition at line 86 of file G4NuclearRadii.cc.

{
  G4double R = ExplicitRadius(Z, A);
  if(0.0 == R) {
    R = 1.24*fG4pow->powZ(A, 0.28)*CLHEP::fermi;
  }
  return R;
}

Member Data Documentation

fG4pow

G4Pow * G4NuclearRadii::fG4pow = G4Pow::GetInstance()

static

Definition at line 94 of file G4NuclearRadii.hh.

Referenced by Radius(), RadiusCB(), RadiusECS(), RadiusHNGG(), RadiusKNGG(), RadiusND(), RadiusNNGG(), and RadiusRMS().

r0

const G4double G4NuclearRadii::r0

static

Initial value:

= {
 1.2,
 1.3, 1.3, 1.3, 1.3,1.17,1.54,1.65,1.71, 1.7,1.75, 
 1.7,1.57,1.53, 1.4, 1.3,1.30,1.44, 1.4, 1.4, 1.4, 
 1.4, 1.4,1.46, 1.4, 1.4,1.46,1.55, 1.5,1.38,1.48, 
 1.4, 1.4, 1.4,1.46, 1.4, 1.4, 1.4, 1.4, 1.4,1.45, 
 1.4, 1.4, 1.4, 1.4, 1.4, 1.4,1.45,1.48, 1.4,1.52, 
1.46, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.5, 
 1.4, 1.4, 1.4, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.4, 
 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3,1.33,1.43, 
 1.3,1.32,1.34, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 
 1.3, 1.3}

Definition at line 253 of file G4NuclearRadii.hh.

Referenced by RadiusCB().

---

The documentation for this class was generated from the following files:

• G4NuclearRadii.hh
• G4NuclearRadii.cc