Geant4  10.02.p03
G4CoulombBarrier Class Reference

#include <G4CoulombBarrier.hh>

Inheritance diagram for G4CoulombBarrier:
Collaboration diagram for G4CoulombBarrier:

Public Member Functions

 G4CoulombBarrier ()
 
 G4CoulombBarrier (G4int anA, G4int aZ)
 
virtual ~G4CoulombBarrier ()
 
G4double GetCoulombBarrier (G4int ARes, G4int ZRes, G4double U) 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

 G4CoulombBarrier (const G4CoulombBarrier &right)
 
const G4CoulombBarrieroperator= (const G4CoulombBarrier &right)
 
G4bool operator== (const G4CoulombBarrier &right) const
 
G4bool operator!= (const G4CoulombBarrier &right) const
 
virtual G4double BarrierPenetrationFactor (G4double) const
 
G4double CalcCompoundRadius (const G4double ZRes) const
 

Detailed Description

Definition at line 43 of file G4CoulombBarrier.hh.

Constructor & Destructor Documentation

◆ G4CoulombBarrier() [1/3]

G4CoulombBarrier::G4CoulombBarrier ( )

Definition at line 42 of file G4CoulombBarrier.cc.

42  : G4VCoulombBarrier(1,0)
43 {}

◆ G4CoulombBarrier() [2/3]

G4CoulombBarrier::G4CoulombBarrier ( G4int  anA,
G4int  aZ 
)

Definition at line 45 of file G4CoulombBarrier.cc.

46  : G4VCoulombBarrier(anA,aZ)
47 {}

◆ ~G4CoulombBarrier()

G4CoulombBarrier::~G4CoulombBarrier ( )
virtual

Definition at line 49 of file G4CoulombBarrier.cc.

50 {}

◆ G4CoulombBarrier() [3/3]

G4CoulombBarrier::G4CoulombBarrier ( const G4CoulombBarrier right)
private

Member Function Documentation

◆ BarrierPenetrationFactor()

G4double G4CoulombBarrier::BarrierPenetrationFactor ( G4double  ) const
privatevirtual

Reimplemented in G4ProtonCoulombBarrier, G4He3CoulombBarrier, G4AlphaCoulombBarrier, G4DeuteronCoulombBarrier, and G4TritonCoulombBarrier.

Definition at line 52 of file G4CoulombBarrier.cc.

53 {
54  return 1.0;
55 }
Here is the caller graph for this function:

◆ CalcCompoundRadius()

G4double G4CoulombBarrier::CalcCompoundRadius ( const G4double  ZRes) const
inlineprivate

Definition at line 63 of file G4CoulombBarrier.hh.

64  {
65  return 2.173*CLHEP::fermi*(1.0+0.006103*static_cast<G4double>(GetZ())*ZRes)/
66  (1.0+0.009443*static_cast<G4double>(GetZ())*ZRes);
67  }
static const double fermi
Definition: SystemOfUnits.h:82
G4int GetZ(void) const
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:

◆ GetCoulombBarrier()

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

New coulomb Barrier according to original Dostrovski's paper

Implements G4VCoulombBarrier.

Definition at line 57 of file G4CoulombBarrier.cc.

60 {
61  G4double Barrier = 0.0;
62  if (GetA() == 1 && GetZ() == 0) {
63  Barrier = 0.0; // Neutron Coulomb Barrier is 0
64  } else {
65 
66  // JMQ: old coulomb barrier commented since it does not agree with
67  // Dostrovski's prescription
68  // and too low barriers are obtained (for protons at least)
69 
71  G4double rho=1.2*fermi;
72  if(GetA()==1 && GetZ()==1){ rho=0.0;}
73 
74  const G4double RN=1.5*fermi;
75  // VI cleanup
76  Barrier=elm_coupling*(GetZ()*ZRes)/(RN * G4Pow::GetInstance()->Z13(ARes) + rho);
77 
78  // Barrier penetration coeficient
80 
81  Barrier *= K;
82 
83  // JMQ : the following statement has unknown origin and dimensionally
84  // is meaningless( energy divided by mass number in argument of sqrt function).
85  // Energy dependence of Coulomb barrier penetrability should be included in
86  // proper way (if needed..)
87  // Barrier /= (1.0 + std::sqrt(U/(2.0*static_cast<G4double>(ARes))));
88  //
89  }
90  return Barrier;
91 }
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
int elm_coupling
Definition: hepunit.py:286
G4int GetA(void) const
G4double Z13(G4int Z) const
Definition: G4Pow.hh:127
G4int GetZ(void) const
double G4double
Definition: G4Types.hh:76
static const double fermi
Definition: G4SIunits.hh:102
virtual G4double BarrierPenetrationFactor(G4double) const
Here is the call graph for this function:

◆ operator!=()

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

◆ operator=()

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

◆ operator==()

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

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