#include <G4StopTheoDeexcitation.hh>
Inheritance diagram for G4StopTheoDeexcitation:
Public Member Functions | |
G4StopTheoDeexcitation () | |
virtual | ~G4StopTheoDeexcitation () |
virtual G4ReactionProductVector * | BreakUp (G4double A, G4double Z, G4double excitation, const G4ThreeVector &p) |
Definition at line 44 of file G4StopTheoDeexcitation.hh.
G4StopTheoDeexcitation::G4StopTheoDeexcitation | ( | ) |
Definition at line 53 of file G4StopTheoDeexcitation.cc.
References G4HadronicDeprecate.
00053 { 00054 G4HadronicDeprecate("G4StopTheoDeexcitation"); 00055 }
G4StopTheoDeexcitation::~G4StopTheoDeexcitation | ( | ) | [virtual] |
G4ReactionProductVector * G4StopTheoDeexcitation::BreakUp | ( | G4double | A, | |
G4double | Z, | |||
G4double | excitation, | |||
const G4ThreeVector & | p | |||
) | [virtual] |
Implements G4StopDeexcitationAlgorithm.
Definition at line 63 of file G4StopTheoDeexcitation.cc.
References G4ExcitationHandler::BreakItUp(), G4NucleiProperties::GetNuclearMass(), and G4ExcitationHandler::SetMinEForMultiFrag().
00065 { 00066 G4ExcitationHandler theHandler; 00067 00068 // MF and FB parameters modified by MGP to force evaporation 00069 // Max A and Z values for use Fermi Breakup 00070 // theHandler.SetMaxAandZForFermiBreakUp(16, 10); 00071 // theHandler.SetMaxAandZForFermiBreakUp(2, 1); 00072 // Min excitation energy (per nucleon) for use MultiFrag 00073 00074 theHandler.SetMinEForMultiFrag(300*GeV); 00075 00076 // Deexcite the nucleus 00077 00078 G4double atomicMass = G4NucleiProperties::GetNuclearMass(static_cast<G4int>(A),static_cast<G4int>(Z)); 00079 G4double mass = atomicMass + excitation; 00080 G4double pMag = p.mag(); 00081 G4LorentzVector initialMomentum(p.x(),p.y(),p.z(),std::sqrt(pMag*pMag + mass*mass)); 00082 G4Fragment theExcitedNucleus(static_cast<G4int>(A),static_cast<G4int>(Z),initialMomentum); 00083 00084 return theHandler.BreakItUp(theExcitedNucleus); 00085 }