G4StatMFFragment Class Reference

#include <G4StatMFFragment.hh>

Public Member Functions
	G4StatMFFragment (G4int anA, G4int aZ)
virtual	~G4StatMFFragment ()
G4bool	operator== (const G4StatMFFragment &right) const
G4bool	operator!= (const G4StatMFFragment &right) const
G4double	GetCoulombEnergy (void) const
G4double	GetEnergy (const G4double T) const
G4double	GetInvLevelDensity (void) const
G4int	GetA (void) const
G4int	GetZ (void) const
void	SetPosition (const G4ThreeVector aPosition)
G4ThreeVector	GetPosition (void)
void	SetMomentum (const G4ThreeVector aMomentum)
G4ThreeVector	GetMomentum (void)
G4Fragment *	GetFragment (const G4double T)
G4double	GetNuclearMass (void)

Detailed Description

Definition at line 41 of file G4StatMFFragment.hh.

Constructor & Destructor Documentation

G4StatMFFragment::G4StatMFFragment ( G4int  anA, G4int  aZ  )

inline

Definition at line 45 of file G4StatMFFragment.hh.

                                          :
    theA(anA),theZ(aZ),
    _position(0.0,0.0,0.0),
    _momentum(0.0,0.0,0.0)
    {}

virtual G4StatMFFragment::~G4StatMFFragment ( )

inlinevirtual

Definition at line 53 of file G4StatMFFragment.hh.

{};

Member Function Documentation

G4int G4StatMFFragment::GetA ( void  ) const

inline

Definition at line 77 of file G4StatMFFragment.hh.

{return theA;}

G4double G4StatMFFragment::GetCoulombEnergy

(

void

)

const

Definition at line 67 of file G4StatMFFragment.cc.

References python.hepunit::elm_coupling, G4StatMFParameters::GetKappaCoulomb(), and G4StatMFParameters::Getr0().

Referenced by GetEnergy(), and SumCoulombEnergy::operator()().

{
    if (theZ <= 0.1) return 0.0;
    G4double Coulomb = (3./5.)*(elm_coupling*theZ*theZ)*
    std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.)/
    (G4StatMFParameters::Getr0()*std::pow(theA,1./3.));
                        
    return Coulomb;
}

G4double G4StatMFFragment::GetEnergy

(

const G4double

T

)

const

Definition at line 78 of file G4StatMFFragment.cc.

References G4StatMFParameters::DBetaDT(), G4cerr, G4endl, G4StatMFParameters::GetBeta0(), GetCoulombEnergy(), G4StatMFParameters::GetCriticalTemp(), GetInvLevelDensity(), and G4NucleiProperties::GetMassExcess().

{
    if (theA < 1 || theZ < 0 || theZ > theA) {
    G4cerr << "G4StatMFFragment::GetEnergy: A = " << theA 
           << ", Z = " << theZ << G4endl;
    throw G4HadronicException(__FILE__, __LINE__, 
        "G4StatMFFragment::GetEnergy: Wrong values for A and Z!");
    }
    G4double BulkEnergy = G4NucleiProperties::GetMassExcess(static_cast<G4int>(theA),
                                static_cast<G4int>(theZ));
    
    if (theA < 4) return BulkEnergy - GetCoulombEnergy();
    
    G4double SurfaceEnergy;
    if (G4StatMFParameters::DBetaDT(T) == 0.0) SurfaceEnergy = 0.0;
    else SurfaceEnergy = (5./2.)*std::pow(theA,2.0/3.0)*T*T*
         G4StatMFParameters::GetBeta0()/
         (G4StatMFParameters::GetCriticalTemp()*
          G4StatMFParameters::GetCriticalTemp());
                     
                     
    G4double ExchangeEnergy = theA*T*T/GetInvLevelDensity();
    if (theA != 4) ExchangeEnergy += SurfaceEnergy;       
    
    return  BulkEnergy + ExchangeEnergy - GetCoulombEnergy();       
    
}

G4Fragment * G4StatMFFragment::GetFragment

(

const G4double

T

)

Definition at line 118 of file G4StatMFFragment.cc.

References GetNuclearMass(), and CLHEP::Hep3Vector::mag2().

{
    G4double U = CalcExcitationEnergy(T);
    
    G4double M = GetNuclearMass();

    G4LorentzVector FourMomentum(_momentum,std::sqrt(_momentum.mag2()+(M+U)*(M+U)));

    G4Fragment * theFragment = new G4Fragment(static_cast<G4int>(theA),static_cast<G4int>(theZ),FourMomentum);

    return theFragment;
}

G4double G4StatMFFragment::GetInvLevelDensity

(

void

)

const

Definition at line 107 of file G4StatMFFragment.cc.

References G4StatMFParameters::GetEpsilon0().

Referenced by GetEnergy().

{
    // Calculate Inverse Density Level
    // Epsilon0*(1 + 3 /(Af - 1))
    if (theA == 1) return 0.0;
    else return
       G4StatMFParameters::GetEpsilon0()*(1.0+3.0/(theA - 1.0));
}

G4ThreeVector G4StatMFFragment::GetMomentum ( void  )

inline

Definition at line 87 of file G4StatMFFragment.hh.

{return _momentum;}

G4double G4StatMFFragment::GetNuclearMass ( void  )

inline

Definition at line 91 of file G4StatMFFragment.hh.

References G4IonTable::GetIonMass(), G4ParticleTable::GetIonTable(), and G4ParticleTable::GetParticleTable().

Referenced by GetFragment().

    {return G4ParticleTable::GetParticleTable()->GetIonTable()
                           ->GetIonMass(theZ, theA);}

G4ThreeVector G4StatMFFragment::GetPosition ( void  )

inline

Definition at line 83 of file G4StatMFFragment.hh.

{return _position;}

G4int G4StatMFFragment::GetZ ( void  ) const

inline

Definition at line 79 of file G4StatMFFragment.hh.

{return theZ;}

G4bool G4StatMFFragment::operator!=

(

const G4StatMFFragment &

right

)

const

Definition at line 59 of file G4StatMFFragment.cc.

{
//  throw G4HadronicException(__FILE__, __LINE__, "G4StatMFFragment::operator!= meant to not be accessable");
    return true;
}

G4bool G4StatMFFragment::operator==

(

const G4StatMFFragment &

right

)

const

Definition at line 52 of file G4StatMFFragment.cc.

{
//  throw G4HadronicException(__FILE__, __LINE__, "G4StatMFFragment::operator== meant to not be accessable");
    return false;
}

void G4StatMFFragment::SetMomentum ( const G4ThreeVector  aMomentum )

inline

Definition at line 85 of file G4StatMFFragment.hh.

{_momentum = aMomentum;}

void G4StatMFFragment::SetPosition ( const G4ThreeVector  aPosition )

inline

Definition at line 81 of file G4StatMFFragment.hh.

{_position = aPosition;}

---

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

• G4StatMFFragment.hh
• G4StatMFFragment.cc