Geant4  10.02.p03
G4GEMCoulombBarrierHE Class Reference

#include <G4GEMCoulombBarrierHE.hh>

Inheritance diagram for G4GEMCoulombBarrierHE:
Collaboration 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
 

Private Member Functions

 G4GEMCoulombBarrierHE ()
 
 G4GEMCoulombBarrierHE (const G4GEMCoulombBarrierHE &right)
 
const G4GEMCoulombBarrierHEoperator= (const G4GEMCoulombBarrierHE &right)
 
G4bool operator== (const G4GEMCoulombBarrierHE &right) const
 
G4bool operator!= (const G4GEMCoulombBarrierHE &right) const
 

Detailed Description

Definition at line 39 of file G4GEMCoulombBarrierHE.hh.

Constructor & Destructor Documentation

◆ G4GEMCoulombBarrierHE() [1/3]

G4GEMCoulombBarrierHE::G4GEMCoulombBarrierHE ( G4int  anA,
G4int  aZ 
)

Definition at line 37 of file G4GEMCoulombBarrierHE.cc.

37  :
38  G4VCoulombBarrier(anA,aZ)
39 {}

◆ ~G4GEMCoulombBarrierHE()

G4GEMCoulombBarrierHE::~G4GEMCoulombBarrierHE ( )

Definition at line 41 of file G4GEMCoulombBarrierHE.cc.

42 {}

◆ G4GEMCoulombBarrierHE() [2/3]

G4GEMCoulombBarrierHE::G4GEMCoulombBarrierHE ( )
private

◆ G4GEMCoulombBarrierHE() [3/3]

G4GEMCoulombBarrierHE::G4GEMCoulombBarrierHE ( const G4GEMCoulombBarrierHE right)
private

Member Function Documentation

◆ BarrierPenetrationFactor()

virtual G4double G4GEMCoulombBarrierHE::BarrierPenetrationFactor ( G4double  ) const
inlinevirtual

Definition at line 59 of file G4GEMCoulombBarrierHE.hh.

60  {return 1.0;};
Here is the call graph for this function:

◆ CalcCompoundRadius()

G4double G4GEMCoulombBarrierHE::CalcCompoundRadius ( G4int  ARes) const

Definition at line 71 of file G4GEMCoulombBarrierHE.cc.

72 {
73  G4Pow* g4pow = G4Pow::GetInstance();
74  G4double AresOneThird = g4pow->Z13(ARes);
75  G4double AejectOneThird = g4pow->Z13(GetA());
76 
77  G4double Result = 1.12*(AresOneThird + AejectOneThird) -
78  0.86*(AresOneThird+AejectOneThird)/(AresOneThird*AejectOneThird);
79 
80  return Result*fermi;
81 }
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
Definition: G4Pow.hh:56
G4int GetA(void) const
G4double Z13(G4int Z) const
Definition: G4Pow.hh:127
double G4double
Definition: G4Types.hh:76
static const double fermi
Definition: G4SIunits.hh:102
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetCoulombBarrier()

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

Implements G4VCoulombBarrier.

Definition at line 44 of file G4GEMCoulombBarrierHE.cc.

46 {
47  G4double Barrier = 0.0;
48  if (ZRes > ARes || ARes < 1) {
49  G4cout << "G4GEMCoulombBarrierHE::GetCoulombBarrier: "
50  << "Wrong values for "
51  << "residual nucleus A = " << ARes << " "
52  << "and residual nucleus Z = " << ZRes << G4endl;
53  throw G4HadronicException(__FILE__, __LINE__,"FATAL error");
54  }
55  if (GetZ() == 0) {
56  Barrier = 0.0; // If there is no charge there is neither barrier
57  } else {
58  G4double CompoundRadius = CalcCompoundRadius(ARes);
59  Barrier = ( elm_coupling * GetZ() * ZRes)/(CompoundRadius+3.75*fermi);
60 
61  // Barrier penetration coeficient
62  // G4double K = BarrierPenetrationFactor(ZRes);
63  // Barrier *= K;
64 
65  Barrier /= (1.0 + std::sqrt(U/static_cast<G4double>(2*ARes)));
66  }
67  return Barrier;
68 }
int elm_coupling
Definition: hepunit.py:286
G4GLOB_DLL std::ostream G4cout
G4double CalcCompoundRadius(G4int ARes) const
G4int GetZ(void) const
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76
static const double fermi
Definition: G4SIunits.hh:102
Here is the call graph for this function:

◆ operator!=()

G4bool G4GEMCoulombBarrierHE::operator!= ( const G4GEMCoulombBarrierHE right) const
private

◆ operator=()

const G4GEMCoulombBarrierHE& G4GEMCoulombBarrierHE::operator= ( const G4GEMCoulombBarrierHE right)
private

◆ operator==()

G4bool G4GEMCoulombBarrierHE::operator== ( const G4GEMCoulombBarrierHE right) const
private

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