#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

virtual G4double	GetCoulombBarrier (G4int ARes, G4int ZRes, G4double U=0.0) const =0

virtual G4double	BarrierPenetrationFactor (G4int aZ) const

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
G4Pow *	g4calc

G4int	theA

G4int	theZ

G4double	theR0

G4double	theRho = 0.0

G4double	factor = 0.0

Detailed Description

Definition at line 38 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);
}

◆ ~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.

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)*CLHEP::MeV) )); }
  return cb;
}

Referenced by G4EvaporationChannel::ComputeProbability(), and G4EvaporationChannel::GetEmissionProbability().

◆ operator=()

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