#include <G4CascadeDeexciteBase.hh>

Inheritance diagram for G4CascadeDeexciteBase:

Public Member Functions
	G4CascadeDeexciteBase (const char *name)

virtual	~G4CascadeDeexciteBase ()

virtual void	setVerboseLevel (G4int verbose=0)

Public Member Functions inherited from G4VCascadeDeexcitation
	G4VCascadeDeexcitation (const char *name)

virtual	~G4VCascadeDeexcitation ()

virtual void	collide (G4InuclParticle bullet, G4InuclParticle target, G4CollisionOutput &globalOutput)

virtual void	deExcite (const G4Fragment &fragment, G4CollisionOutput &output)=0

Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const char *name, G4int verbose=0)

virtual	~G4VCascadeCollider ()

Protected Member Functions
virtual G4bool	explosion (const G4Fragment &target) const

virtual G4bool	explosion (G4int A, G4int Z, G4double excitation) const

virtual G4bool	validateOutput (const G4Fragment &target, G4CollisionOutput &output)

virtual G4bool	validateOutput (const G4Fragment &target, const std::vector< G4InuclElementaryParticle > &particles)

virtual G4bool	validateOutput (const G4Fragment &target, const std::vector< G4InuclNuclei > &fragments)

void	getTargetData (const G4Fragment &target)

const G4Fragment &	makeFragment (G4LorentzVector mom, G4int A, G4int Z, G4double EX=0.)

const G4Fragment &	makeFragment (G4int A, G4int Z, G4double EX=0.)

Protected Member Functions inherited from G4VCascadeCollider
virtual void	setName (const char *name)

Protected Attributes
G4CascadeCheckBalance *	balance

G4int	A

G4int	Z

G4LorentzVector	PEX

G4double	EEXS

G4Fragment	aFragment

Protected Attributes inherited from G4VCascadeCollider
const char *	theName

G4int	verboseLevel

Detailed Description

Definition at line 45 of file G4CascadeDeexciteBase.hh.

Constructor & Destructor Documentation

G4CascadeDeexciteBase::G4CascadeDeexciteBase ( const char * name )

Definition at line 49 of file G4CascadeDeexciteBase.cc.

References balance, and G4CascadeParameters::checkConservation().

   : G4VCascadeDeexcitation(name), balance(0), A(0), Z(0), EEXS(0) {
   if (G4CascadeParameters::checkConservation())
     balance = new G4CascadeCheckBalance(name);
 }

G4CascadeDeexciteBase::~G4CascadeDeexciteBase ( )

virtual

Definition at line 55 of file G4CascadeDeexciteBase.cc.

References balance.

                                               {
   delete balance;
 }

Member Function Documentation

G4bool G4CascadeDeexciteBase::explosion ( const G4Fragment & target ) const

protectedvirtual

Definition at line 110 of file G4CascadeDeexciteBase.cc.

References G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), and G4Fragment::GetZ_asInt().

Referenced by G4CascadeDeexcitation::deExcite(), and G4PreCompoundDeexcitation::deExcite().

                                                                         {
   return explosion(fragment.GetA_asInt(), fragment.GetZ_asInt(),
            fragment.GetExcitationEnergy());     // in MeV
 }

G4bool G4CascadeDeexciteBase::explosion	(	G4int	A,
		G4int	Z,
		G4double	excitation
	)		const

protectedvirtual

Definition at line 115 of file G4CascadeDeexciteBase.cc.

References G4InuclSpecialFunctions::bindingEnergy(), G4cout, G4endl, G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                 {
   if (verboseLevel) G4cout << " >>> " << theName << "::explosion ?" << G4endl;
 
   const G4int a_cut = 20;
   const G4double be_cut = 3.0;
 
   // Neutron balls, or small fragments with high excitations can explode
   return ((fragA <= a_cut || fragZ==0) && 
       (excitation >= be_cut * bindingEnergy(fragA,fragZ))
       );
 }

void G4CascadeDeexciteBase::getTargetData ( const G4Fragment & target )

protected

Definition at line 67 of file G4CascadeDeexciteBase.cc.

References A, EEXS, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetMomentum(), G4Fragment::GetZ_asInt(), python.hepunit::GeV, PEX, and Z.

Referenced by G4NonEquilibriumEvaporator::deExcite(), G4BigBanger::deExcite(), and G4EquilibriumEvaporator::deExcite().

                                                                   {
   A = target.GetA_asInt();
   Z = target.GetZ_asInt();
   PEX = target.GetMomentum()/GeV;   // Convert from G4 to Bertini units
   EEXS = target.GetExcitationEnergy();
 }

const G4Fragment & G4CascadeDeexciteBase::makeFragment	(	G4LorentzVector	mom,
		G4int	A,
		G4int	Z,
		G4double	EX = `0.`
	)

protected

Definition at line 87 of file G4CascadeDeexciteBase.cc.

References aFragment, G4cout, G4endl, G4InuclNuclei::getNucleiMass(), python.hepunit::GeV, G4Fragment::SetMomentum(), G4Fragment::SetNumberOfExcitedParticle(), G4Fragment::SetNumberOfHoles(), CLHEP::HepLorentzVector::setVectM(), G4Fragment::SetZandA_asInt(), G4VCascadeCollider::theName, CLHEP::HepLorentzVector::vect(), and G4VCascadeCollider::verboseLevel.

Referenced by G4NonEquilibriumEvaporator::deExcite(), G4EquilibriumEvaporator::deExcite(), and makeFragment().

                                               {
   if (verboseLevel>2) {
     G4cout << " >>> " << theName << "::makeFragment " << mom << " " << fragA
        << " " << fragZ << " " << EX << G4endl;
   }
 
   // Adjust four-momentum so that mass is nucleus + excitation
   G4double mass =
     G4InuclNuclei::getNucleiMass(fragA,fragZ) + EX/GeV;
   mom.setVectM(mom.vect(), mass);
 
   // Overwrite previous fragment contents, zeroing out excitons
   aFragment.SetZandA_asInt(fragZ, fragA);
   aFragment.SetMomentum(mom*GeV);           // Bertini uses GeV!
   aFragment.SetNumberOfHoles(0,0);
   aFragment.SetNumberOfExcitedParticle(0,0);
 
   return aFragment;
 }

const G4Fragment & G4CascadeDeexciteBase::makeFragment	(	G4int	A,
		G4int	Z,
		G4double	EX = `0.`
	)

protected

Definition at line 82 of file G4CascadeDeexciteBase.cc.

References makeFragment().

                                                                          {
   return makeFragment(zero, fragA, fragZ, EX);
 }

void G4CascadeDeexciteBase::setVerboseLevel ( G4int verbose = 0 )

virtual

Reimplemented from G4VCascadeCollider.

Reimplemented in G4EquilibriumEvaporator, G4PreCompoundDeexcitation, and G4CascadeDeexcitation.

Definition at line 59 of file G4CascadeDeexciteBase.cc.

References balance, and G4VCascadeCollider::setVerboseLevel().

Referenced by G4CascadeDeexcitation::setVerboseLevel(), G4PreCompoundDeexcitation::setVerboseLevel(), and G4EquilibriumEvaporator::setVerboseLevel().

                                                          {
   G4VCascadeDeexcitation::setVerboseLevel(verbose);
   if (balance) balance->setVerboseLevel(verbose);
 }

G4bool G4CascadeDeexciteBase::validateOutput	(	const G4Fragment &	target,
		G4CollisionOutput &	output
	)

protectedvirtual

Definition at line 131 of file G4CascadeDeexciteBase.cc.

References balance, G4CascadeCheckBalance::collide(), G4cout, G4endl, G4CascadeCheckBalance::okay(), G4VCascadeCollider::setVerboseLevel(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by G4NonEquilibriumEvaporator::deExcite(), G4BigBanger::deExcite(), and G4EquilibriumEvaporator::deExcite().

                                                     {
   if (!balance) return true;        // Skip checks unless requested
 
   if (verboseLevel > 1)
     G4cout << " >>> " << theName << "::validateOutput" << G4endl;
 
   balance->setVerboseLevel(verboseLevel);
   balance->collide(target, output);
   return balance->okay();           // Returns false if violations
 }

G4bool G4CascadeDeexciteBase::validateOutput	(	const G4Fragment &	target,
		const std::vector< G4InuclElementaryParticle > &	particles
	)

protectedvirtual

Definition at line 143 of file G4CascadeDeexciteBase.cc.

References balance, G4CascadeCheckBalance::collide(), G4cout, G4endl, G4CascadeCheckBalance::okay(), G4VCascadeCollider::setVerboseLevel(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                                     {
   if (!balance) return true;        // Skip checks unless requested
 
   if (verboseLevel > 1)
     G4cout << " >>> " << theName << "::validateOutput" << G4endl;
 
   balance->setVerboseLevel(verboseLevel);
   balance->collide(target, particles);
   return balance->okay();           // Returns false if violations
 }

G4bool G4CascadeDeexciteBase::validateOutput	(	const G4Fragment &	target,
		const std::vector< G4InuclNuclei > &	fragments
	)

protectedvirtual

Definition at line 155 of file G4CascadeDeexciteBase.cc.

References balance, G4CascadeCheckBalance::collide(), G4cout, G4endl, G4CascadeCheckBalance::okay(), G4VCascadeCollider::setVerboseLevel(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                         {
   if (!balance) return true;        // Skip checks unless requested
 
   if (verboseLevel > 1)
     G4cout << " >>> " << theName << "::validateOutput" << G4endl;
 
   balance->setVerboseLevel(verboseLevel);
   balance->collide(target, fragments);
   return balance->okay();           // Returns false if violations
 }