G4FermiDecayProbability Class Reference

#include <G4FermiDecayProbability.hh>

Public Member Functions
G4double	ComputeProbability (G4int Z, G4int A, G4int spin, G4double TotalE, const G4FermiFragment *f1, const G4FermiFragment *f2) const
	G4FermiDecayProbability ()
	~G4FermiDecayProbability ()

Private Member Functions
	G4FermiDecayProbability (const G4FermiDecayProbability &right)=delete
G4bool	operator!= (const G4FermiDecayProbability &right) const =delete
const G4FermiDecayProbability &	operator= (const G4FermiDecayProbability &right)=delete
G4bool	operator== (const G4FermiDecayProbability &right) const =delete

Detailed Description

Definition at line 38 of file G4FermiDecayProbability.hh.

Constructor & Destructor Documentation

G4FermiDecayProbability() [1/2]

G4FermiDecayProbability::G4FermiDecayProbability ( )

explicit

Definition at line 35 of file G4FermiDecayProbability.cc.

{}

~G4FermiDecayProbability()

G4FermiDecayProbability::~G4FermiDecayProbability

(

)

Definition at line 38 of file G4FermiDecayProbability.cc.

{}

G4FermiDecayProbability() [2/2]

G4FermiDecayProbability::G4FermiDecayProbability ( const G4FermiDecayProbability &  right )

privatedelete

Member Function Documentation

ComputeProbability()

G4double G4FermiDecayProbability::ComputeProbability	(	G4int	Z,
		G4int	A,
		G4int	spin,
		G4double	TotalE,
		const G4FermiFragment *	f1,
		const G4FermiFragment *	f2
	)		const

Definition at line 42 of file G4FermiDecayProbability.cc.

{
  G4double prob = 0.0;
  G4double mass1 = f1->GetTotalEnergy();
  G4double mass2 = f2->GetTotalEnergy();
  G4double bCouloumb = f1->GetCoulombBarrier(f2->GetA(), f2->GetZ(), 0.0);
  if(etot <= mass1 + mass2 + bCouloumb) { return prob; }
  
  //G4cout << "ComputeProbability M1= " << mass1 << " M2= " << mass2 << G4endl;
  G4double ekin = etot - mass1 - mass2;
 
  // mass factors
  G4double massFactor = mass1*mass2/(mass1 + mass2);
  massFactor *= std::sqrt(massFactor);
 
  // Spin factor S_n
  G4double S_n = 1.0;
  if(spin >= 0) {
    G4int spin1 = f1->GetSpin();
    G4int spin2 = f2->GetSpin();
    if(spin1 >= 0 && spin2 >= 0) {
      S_n = (spin1+1)*(spin2+1); 
    }
  }
    
  // Permutation Factor G_n
  // search for identical fragments
  G4double G_n = (f1 == f2) ? 0.5 : 1.0;
 
  prob = A*massFactor*S_n*G_n*std::sqrt(ekin);
 
  //G4cout << "prob= " << prob << " Coeff= " << Coeff << G4endl;
  return prob; 
}

References A, G4FermiFragment::GetA(), G4FermiFragment::GetCoulombBarrier(), G4FermiFragment::GetSpin(), G4FermiFragment::GetTotalEnergy(), and G4FermiFragment::GetZ().

Referenced by G4FermiFragmentsPoolVI::Initialise(), and G4FermiBreakUpVI::SampleDecay().

operator!=()

G4bool G4FermiDecayProbability::operator!= ( const G4FermiDecayProbability &  right ) const

privatedelete

operator=()

const G4FermiDecayProbability & G4FermiDecayProbability::operator= ( const G4FermiDecayProbability &  right )

privatedelete

operator==()

G4bool G4FermiDecayProbability::operator== ( const G4FermiDecayProbability &  right ) const

privatedelete

---

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

• source/processes/hadronic/models/de_excitation/fermi_breakup/include/G4FermiDecayProbability.hh
• source/processes/hadronic/models/de_excitation/fermi_breakup/src/G4FermiDecayProbability.cc