#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)

Public Member Functions inherited from G4StopDeexcitationAlgorithm
	G4StopDeexcitationAlgorithm ()

virtual	~G4StopDeexcitationAlgorithm ()

virtual G4ReactionProductVector *	BreakUp (G4double A, G4double Z, G4double excitation, const G4ThreeVector &p)=0

Detailed Description

Definition at line 44 of file G4StopTheoDeexcitation.hh.

Constructor & Destructor Documentation

◆ G4StopTheoDeexcitation()

G4StopTheoDeexcitation::G4StopTheoDeexcitation ( )

Definition at line 53 of file G4StopTheoDeexcitation.cc.

                                               {
  G4HadronicDeprecate("G4StopTheoDeexcitation");
}

◆ ~G4StopTheoDeexcitation()

G4StopTheoDeexcitation::~G4StopTheoDeexcitation ( )

virtual

Definition at line 59 of file G4StopTheoDeexcitation.cc.

60{}

Member Function Documentation

◆ BreakUp()

G4ReactionProductVector * G4StopTheoDeexcitation::BreakUp	(	G4double	A,
		G4double	Z,
		G4double	excitation,
		const G4ThreeVector &	p
	)

virtual

Implements G4StopDeexcitationAlgorithm.

Definition at line 63 of file G4StopTheoDeexcitation.cc.

{
  G4ExcitationHandler theHandler;
 
  // MF and FB parameters modified by MGP to force evaporation 
  // Max A and Z values for use Fermi Breakup
  // theHandler.SetMaxAandZForFermiBreakUp(16, 10);
  //  theHandler.SetMaxAandZForFermiBreakUp(2, 1);
  // Min excitation energy (per nucleon) for use MultiFrag
 
  theHandler.SetMinEForMultiFrag(300*GeV);
 
  // Deexcite the nucleus 
 
  G4double atomicMass = G4NucleiProperties::GetNuclearMass(static_cast<G4int>(A),static_cast<G4int>(Z));
  G4double mass = atomicMass + excitation;
  G4double pMag = p.mag();
  G4LorentzVector initialMomentum(p.x(),p.y(),p.z(),std::sqrt(pMag*pMag + mass*mass));
  G4Fragment theExcitedNucleus(static_cast<G4int>(A),static_cast<G4int>(Z),initialMomentum);
 
  return theHandler.BreakItUp(theExcitedNucleus);
}