G4GEMCoulombBarrier Class Reference

#include <G4GEMCoulombBarrier.hh>

Inheritance diagram for G4GEMCoulombBarrier:

Public Member Functions
	G4GEMCoulombBarrier (G4int anA, G4int aZ)
virtual	~G4GEMCoulombBarrier ()
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 37 of file G4GEMCoulombBarrier.hh.

Constructor & Destructor Documentation

G4GEMCoulombBarrier::G4GEMCoulombBarrier	(	G4int	anA,
		G4int	aZ
	)

Definition at line 36 of file G4GEMCoulombBarrier.cc.

References G4GEMCoulombBarrier().

Referenced by G4GEMCoulombBarrier().

                                                            :
  G4VCoulombBarrier(anA,aZ) 
{}

G4GEMCoulombBarrier::~G4GEMCoulombBarrier ( )

virtual

Definition at line 40 of file G4GEMCoulombBarrier.cc.

{}

Member Function Documentation

virtual G4double G4GEMCoulombBarrier::BarrierPenetrationFactor ( G4double  ) const

inlinevirtual

Reimplemented in G4AlphaGEMCoulombBarrier, and G4DeuteronGEMCoulombBarrier.

Definition at line 56 of file G4GEMCoulombBarrier.hh.

  {return 1.0;};

G4double G4GEMCoulombBarrier::CalcCompoundRadius

(

G4int

ARes

)

const

Definition at line 71 of file G4GEMCoulombBarrier.cc.

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

{      
  G4Pow* g4pow = G4Pow::GetInstance();
  G4double AresOneThird = g4pow->Z13(ARes);
  G4int A = GetA();
  G4double AejectOneThird = g4pow->Z13(A);

  G4double Result = 0.0;
  if(A == 1){
    Result = 1.7* AresOneThird;

  } else if (A <= 4){
    Result = 1.7* AresOneThird + 1.2;

  } else {
    Result = 1.12*(AresOneThird + AejectOneThird) - 
      0.86*(AresOneThird+AejectOneThird)/(AresOneThird*AejectOneThird)+3.75;
  }
  return Result*fermi;
}

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

virtual

Implements G4VCoulombBarrier.

Definition at line 43 of file G4GEMCoulombBarrier.cc.

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

{
  G4double Barrier = 0.0;
  if (ZRes > ARes || ARes < 1) {
    G4cout << "G4GEMCoulombBarrier::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;
      
    // Barrier penetration coeficient
    if(GetA() <= 4) { Barrier *= BarrierPenetrationFactor(G4double(ZRes)); }
  
    //JMQ 200709 effective decrease  of barrier with E* (Barashenkov)
    // (not inclued in original Furihata's formulation)
    Barrier /= (1.0 + std::sqrt(U/(static_cast<G4double>(2*ARes))));
  }
  return Barrier;
}

---

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

• G4GEMCoulombBarrier.hh
• G4GEMCoulombBarrier.cc