Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4CoulombBarrier Class Reference

#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 G4CoulombBarrieroperator= (const G4CoulombBarrier &right)=delete
 
- Public Member Functions inherited from G4VCoulombBarrier
 G4VCoulombBarrier (G4int anA, G4int aZ)
 
virtual ~G4VCoulombBarrier ()=default
 
void SetParameters (G4double rho, G4double r0)
 
 G4VCoulombBarrier (const G4VCoulombBarrier &right)=delete
 
const G4VCoulombBarrieroperator= (const G4VCoulombBarrier &right)=delete
 

Additional Inherited Members

- Protected Attributes inherited from G4VCoulombBarrier
G4Powg4calc
 
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.

40{
41 factor = CLHEP::elm_coupling*Z;
42 SetParameters(0.4*G4NuclearRadii::RadiusCB(Z, A), 1.5*CLHEP::fermi);
43}
const G4double A[17]
static G4double RadiusCB(G4int Z, G4int A)
G4VCoulombBarrier(G4int anA, G4int aZ)
void SetParameters(G4double rho, G4double r0)

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

55{
56 // Data comes from
57 // Dostrovsky, Fraenkel and Friedlander
58 // Physical Review, vol 116, num. 3 1959
59 //
60 // const G4int size = 5;
61 // const G4double Zlist[size] = {10.0, 20.0, 30.0, 50.0, 70.0};
62 // const G4double Kprot[size] = {0.42, 0.58, 0.68, 0.77, 0.80};
63 //
64 G4double res = 1.0;
65 if(theZ == 1) {
66 res = (aZ >= 70) ? 0.80 :
67 (((0.2357e-5*aZ) - 0.42679e-3)*aZ + 0.27035e-1)*aZ + 0.19025;
68 res += 0.06*(theA - 1);
69
70 } else if(theZ == 2 && theA <= 4) {
71 res = (aZ >= 70) ? 0.98 :
72 (((0.23684e-5*aZ) - 0.42143e-3)*aZ + 0.25222e-1)*aZ + 0.46699;
73 res += 0.12*(4 - theA);
74 }
75 return res;
76}
double G4double
Definition G4Types.hh:83

Referenced by G4GEMCoulombBarrier::GetCoulombBarrier().

◆ GetCoulombBarrier()

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

Implements G4VCoulombBarrier.

Reimplemented in G4GEMCoulombBarrier.

Definition at line 45 of file G4CoulombBarrier.cc.

47{
48 if(0 == theZ) { return 0.0; }
49 G4double cb = factor*ZRes/(G4NuclearRadii::RadiusCB(ZRes,ARes) + theRho);
50 if(U > 0.0) { cb /= (1.0 + std::sqrt( U/((2*ARes)*CLHEP::MeV) )); }
51 return cb;
52}

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

◆ operator=()

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

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