G4GEMCoulombBarrierHE Class Reference

#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 ()
G4int	GetA (void) const
G4int	GetZ (void) const

Detailed Description

Definition at line 39 of file G4GEMCoulombBarrierHE.hh.

Constructor & Destructor Documentation

G4GEMCoulombBarrierHE::G4GEMCoulombBarrierHE	(	G4int	anA,
		G4int	aZ
	)

Definition at line 37 of file G4GEMCoulombBarrierHE.cc.

References G4GEMCoulombBarrierHE().

Referenced by G4GEMCoulombBarrierHE().

                                                                :
  G4VCoulombBarrier(anA,aZ) 
{}

G4GEMCoulombBarrierHE::~G4GEMCoulombBarrierHE

(

)

Definition at line 41 of file G4GEMCoulombBarrierHE.cc.

{}

Member Function Documentation

virtual G4double G4GEMCoulombBarrierHE::BarrierPenetrationFactor ( G4double  ) const

inlinevirtual

Definition at line 59 of file G4GEMCoulombBarrierHE.hh.

  {return 1.0;};

G4double G4GEMCoulombBarrierHE::CalcCompoundRadius

(

G4int

ARes

)

const

Definition at line 71 of file G4GEMCoulombBarrierHE.cc.

References python.hepunit::fermi, G4VCoulombBarrier::GetA(), G4Pow::GetInstance(), and G4Pow::Z13().

{
  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;
}

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

virtual

Implements G4VCoulombBarrier.

Definition at line 44 of file G4GEMCoulombBarrierHE.cc.

References python.hepunit::elm_coupling, python.hepunit::fermi, G4cout, and G4endl.

{
  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:

• G4GEMCoulombBarrierHE.hh
• G4GEMCoulombBarrierHE.cc