G4FermiDecayProbability Class Reference

#include <G4FermiDecayProbability.hh>

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

Detailed Description

Definition at line 39 of file G4FermiDecayProbability.hh.

Constructor & Destructor Documentation

G4FermiDecayProbability::G4FermiDecayProbability ( )

explicit

Definition at line 36 of file G4FermiDecayProbability.cc.

{}

G4FermiDecayProbability::~G4FermiDecayProbability

(

)

Definition at line 39 of file G4FermiDecayProbability.cc.

{}

Member Function Documentation

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

Definition at line 43 of file G4FermiDecayProbability.cc.

{
  G4double prob = 0.0;
  G4double mass1 = f1->GetTotalEnergy();
  G4double mass2 = f2->GetTotalEnergy();
  //G4cout << "ComputeProbability M1= " << mass1 << " M2= " << mass2 << G4endl;
  G4double ekin = etot - mass1 - mass2
    - f1->GetCoulombBarrier(f2->GetA(), f2->GetZ(), 0.0); 
  //G4cout << " Ekin= " << ekin << G4endl;
  if(ekin <= 0.0) { return prob; }

  // 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; 
}

Here is the call graph for this function:

---

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

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