#include <G4GEMCoulombBarrierHE.hh>

Inheritance diagram for G4GEMCoulombBarrierHE:

Public Member Functions
	G4GEMCoulombBarrierHE (G4int anA, G4int aZ)

	~G4GEMCoulombBarrierHE ()

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

virtual G4double	BarrierPenetrationFactor (G4double) const

G4double	CalcCompoundRadius (G4int ARes) const

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

virtual	~G4VCoulombBarrier ()

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

G4int	GetA (void) const

G4int	GetZ (void) const

Detailed Description

Definition at line 39 of file G4GEMCoulombBarrierHE.hh.

Constructor & Destructor Documentation

◆ G4GEMCoulombBarrierHE()

G4GEMCoulombBarrierHE::G4GEMCoulombBarrierHE	(	G4int	anA,
		G4int	aZ
	)

Definition at line 37 of file G4GEMCoulombBarrierHE.cc.

                                                                :
  G4VCoulombBarrier(anA,aZ) 
{}

Referenced by G4GEMCoulombBarrierHE().

◆ ~G4GEMCoulombBarrierHE()

G4GEMCoulombBarrierHE::~G4GEMCoulombBarrierHE ( )

Definition at line 41 of file G4GEMCoulombBarrierHE.cc.

42{}

Member Function Documentation

◆ BarrierPenetrationFactor()

virtual G4double G4GEMCoulombBarrierHE::BarrierPenetrationFactor ( G4double ) const

inlinevirtual

Definition at line 59 of file G4GEMCoulombBarrierHE.hh.

60 {return 1.0;};

◆ CalcCompoundRadius()

G4double G4GEMCoulombBarrierHE::CalcCompoundRadius ( G4int ARes ) const

Definition at line 71 of file G4GEMCoulombBarrierHE.cc.

{
  G4Pow* g4pow = G4Pow::GetInstance();
  G4double AresOneThird = g4pow->Z13(ARes);
  G4double AejectOneThird = g4pow->Z13(GetA());
 
  G4double Result = 1.12*(AresOneThird + AejectOneThird) - 
    0.86*(AresOneThird+AejectOneThird)/(AresOneThird*AejectOneThird);
 
  return Result*fermi;
}

◆ GetCoulombBarrier()

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

virtual

Implements G4VCoulombBarrier.

Definition at line 44 of file G4GEMCoulombBarrierHE.cc.

{
  G4double Barrier = 0.0;
  if (ZRes > ARes || ARes < 1) {
    G4cout << "G4GEMCoulombBarrierHE::GetCoulombBarrier: "
       << "Wrong values for "
       << "residual nucleus A = " << ARes << " "
       << "and residual nucleus Z = " << ZRes << G4endl;
    throw G4HadronicException(__FILE__, __LINE__,"FATAL error");
  }
  if (GetZ() == 0) {
    Barrier = 0.0;   // If there is no charge there is neither barrier
  } else {
    G4double CompoundRadius = CalcCompoundRadius(ARes);
    Barrier = ( elm_coupling * GetZ() * ZRes)/(CompoundRadius+3.75*fermi);
    
    // Barrier penetration coeficient
    //    G4double K = BarrierPenetrationFactor(ZRes);
    //    Barrier *= K;
    
    Barrier /= (1.0 + std::sqrt(U/static_cast<G4double>(2*ARes)));
  }
  return Barrier;
}

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

Public Member Functions