#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 43 of file G4CoulombBarrier.cc.

  : G4VCoulombBarrier(A, Z)
{
  factor = CLHEP::elm_coupling*Z;
  theRho = 0.5*G4NuclearRadii::RadiusCB(Z, A);
}

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 58 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 50 of file G4CoulombBarrier.cc.

{
  G4double cb = factor*ZRes/(G4NuclearRadii::RadiusCB(ZRes,ARes) + theRho);
  cb /= (1.0 + std::sqrt( U / ((ARes + theA)*e0)));
  return cb;
}

◆ operator=()

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