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 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 169 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 190 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 43 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().

ParticleRadius()

G4double G4NuclearRadii::ParticleRadius ( const G4ParticleDefinition *  p )

static

Definition at line 158 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().

Radius()

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

static

Definition at line 59 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 148 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 CoulombFactor(), G4CoulombBarrier::G4CoulombBarrier(), G4FermiCoulombBarrier::G4FermiCoulombBarrier(), and G4CoulombBarrier::GetCoulombBarrier().

RadiusECS()

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

static

Definition at line 101 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 122 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 133 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 138 of file G4NuclearRadii.cc.

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

RadiusNNGG()

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

static

Definition at line 87 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 78 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 90 of file G4NuclearRadii.hh.

Referenced by Radius(), RadiusCB(), RadiusECS(), RadiusHNGG(), RadiusKNGG(), 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 91 of file G4NuclearRadii.hh.

Referenced by RadiusCB().

---

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

• G4NuclearRadii.hh
• G4NuclearRadii.cc