#include <G4CoulombBarrier.hh>

Inheritance diagram for G4CoulombBarrier:

Public Member Functions
	G4CoulombBarrier (G4int anA, G4int aZ)

	~G4CoulombBarrier () override=default

G4double	GetCoulombBarrier (G4int ARes, G4int ZRes, G4double U) const override

G4double	BarrierPenetrationFactor (G4int aZ) const override

	G4CoulombBarrier (const G4CoulombBarrier &right)=delete

const G4CoulombBarrier &	operator= (const G4CoulombBarrier &right)=delete

Public Member Functions inherited from G4VCoulombBarrier
	G4VCoulombBarrier (G4int anA, G4int aZ)

virtual	~G4VCoulombBarrier ()=default

void	SetParameters (G4double rho, G4double r0)

	G4VCoulombBarrier (const G4VCoulombBarrier &right)=delete

const G4VCoulombBarrier &	operator= (const G4VCoulombBarrier &right)=delete

Additional Inherited Members
Protected Attributes inherited from G4VCoulombBarrier
G4int	theA

G4int	theZ

G4double	theRho {0.0}

G4double	factor {0.0}

Detailed Description

Definition at line 37 of file G4CoulombBarrier.hh.

Constructor & Destructor Documentation

◆ G4CoulombBarrier() [1/2]

G4CoulombBarrier::G4CoulombBarrier	(	G4int	anA,
		G4int	aZ )

explicit

Definition at line 38 of file G4CoulombBarrier.cc.

  : G4VCoulombBarrier(A, Z)
{
  factor = CLHEP::elm_coupling*Z;
  SetParameters(0.4*G4NuclearRadii::RadiusCB(Z, A), 1.5*CLHEP::fermi);
}

Referenced by G4CoulombBarrier(), G4GEMCoulombBarrier::G4GEMCoulombBarrier(), and operator=().

◆ ~G4CoulombBarrier()

G4CoulombBarrier::~G4CoulombBarrier ( )

overridedefault

◆ G4CoulombBarrier() [2/2]

G4CoulombBarrier::G4CoulombBarrier ( const G4CoulombBarrier & right )

delete

Member Function Documentation

◆ BarrierPenetrationFactor()

G4double G4CoulombBarrier::BarrierPenetrationFactor ( G4int aZ ) const

overridevirtual

Reimplemented from G4VCoulombBarrier.

Definition at line 54 of file G4CoulombBarrier.cc.

{
  // Data comes from 
  // Dostrovsky, Fraenkel and Friedlander
  // Physical Review, vol 116, num. 3 1959
  // 
  // const G4int size = 5;
  // const G4double Zlist[size] = {10.0, 20.0, 30.0, 50.0, 70.0};
  // const G4double Kprot[size] = {0.42, 0.58, 0.68, 0.77, 0.80};
  // 
  G4double res = 1.0;
  if(theZ == 1) {
    res = (aZ >= 70) ? 0.80 :
    (((0.2357e-5*aZ) - 0.42679e-3)*aZ + 0.27035e-1)*aZ + 0.19025;
    res += 0.06*(theA - 1);
 
  } else if(theZ == 2 && theA <= 4) {
    res = (aZ >= 70) ? 0.98 :
    (((0.23684e-5*aZ) - 0.42143e-3)*aZ + 0.25222e-1)*aZ + 0.46699;
    res += 0.12*(4 - theA);
  }
  return res;
}

Referenced by G4GEMCoulombBarrier::GetCoulombBarrier().

◆ GetCoulombBarrier()

G4double G4CoulombBarrier::GetCoulombBarrier	(	G4int	ARes,
		G4int	ZRes,
		G4double	U ) const

overridevirtual

Implements G4VCoulombBarrier.

Reimplemented in G4GEMCoulombBarrier.

Definition at line 45 of file G4CoulombBarrier.cc.

{
  if (0 == theZ) { return 0.0; }
  G4double cb = factor*ZRes/(G4NuclearRadii::RadiusCB(ZRes,ARes) + theRho);
  if (U > 0.0) { cb /= (1.0 + std::sqrt( U/((2*(ARes + theA))*CLHEP::MeV) )); }
  return cb;
}

◆ operator=()

const G4CoulombBarrier & G4CoulombBarrier::operator= ( const G4CoulombBarrier & right )