#include <G4AblaInterface.hh>

Inheritance diagram for G4AblaInterface:

Collaboration diagram for G4AblaInterface:

Public Member Functions
	G4AblaInterface ()

virtual	~G4AblaInterface ()

virtual G4ReactionProductVector *	DeExcite (G4Fragment &aFragment)

virtual G4HadFinalState *	ApplyYourself (G4HadProjectile const &, G4Nucleus &)

virtual void	ModelDescription (std::ostream &outFile) const

virtual void	DeExciteModelDescription (std::ostream &outFile) const

Public Member Functions inherited from G4VPreCompoundModel
	G4VPreCompoundModel (G4ExcitationHandler *ptr=0, const G4String &modelName="PrecompoundModel")

virtual	~G4VPreCompoundModel ()

void	SetExcitationHandler (G4ExcitationHandler *ptr)

G4ExcitationHandler *	GetExcitationHandler () const

Public Member Functions inherited from G4HadronicInteraction
	G4HadronicInteraction (const G4String &modelName="HadronicModel")

virtual	~G4HadronicInteraction ()

virtual G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)

virtual G4bool	IsApplicable (const G4HadProjectile &, G4Nucleus &)

G4double	GetMinEnergy () const

G4double	GetMinEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMinEnergy (G4double anEnergy)

void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)

G4double	GetMaxEnergy () const

G4double	GetMaxEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMaxEnergy (const G4double anEnergy)

void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)

const G4HadronicInteraction *	GetMyPointer () const

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

const G4String &	GetModelName () const

void	DeActivateFor (const G4Material *aMaterial)

void	ActivateFor (const G4Material *aMaterial)

void	DeActivateFor (const G4Element *anElement)

void	ActivateFor (const G4Element *anElement)

G4bool	IsBlocked (const G4Material *aMaterial) const

G4bool	IsBlocked (const G4Element *anElement) const

void	SetRecoilEnergyThreshold (G4double val)

G4double	GetRecoilEnergyThreshold () const

G4bool	operator== (const G4HadronicInteraction &right) const

G4bool	operator!= (const G4HadronicInteraction &right) const

virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const

virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const

void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

Private Member Functions
G4ReactionProduct *	toG4Particle (G4int A, G4int Z, G4double kinE, G4double px, G4double py, G4double pz) const
	Convert an Abla particle to a G4ReactionProduct. More...

G4ParticleDefinition *	toG4ParticleDefinition (G4int A, G4int Z) const
	Convert A and Z to a G4ParticleDefinition. More...

Private Attributes
G4VarNtp *	ablaResult

G4Volant *	volant

G4Abla *	theABLAModel

G4long	eventNumber

Additional Inherited Members
Protected Member Functions inherited from G4HadronicInteraction
void	SetModelName (const G4String &nam)

G4bool	IsBlocked () const

void	Block ()

Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState	theParticleChange

G4int	verboseLevel

G4double	theMinEnergy

G4double	theMaxEnergy

G4bool	isBlocked

Detailed Description

Definition at line 49 of file G4AblaInterface.hh.

Constructor & Destructor Documentation

◆ G4AblaInterface()

G4AblaInterface::G4AblaInterface ( )

Definition at line 49 of file G4AblaInterface.cc.

                                  :
   G4VPreCompoundModel(NULL, "ABLA"),
   ablaResult(new G4VarNtp),
   volant(new G4Volant),
   theABLAModel(new G4Abla(volant, ablaResult)),
   eventNumber(0)
 {
   theABLAModel->initEvapora();
 }

Here is the call graph for this function:

◆ ~G4AblaInterface()

G4AblaInterface::~G4AblaInterface ( )

virtual

Definition at line 59 of file G4AblaInterface.cc.

                                   {
   delete volant;
   delete ablaResult;
   delete theABLAModel;
 }

Member Function Documentation

◆ ApplyYourself()

virtual G4HadFinalState* G4AblaInterface::ApplyYourself	(	G4HadProjectile const &	,
		G4Nucleus &
	)

inlinevirtual

Implements G4VPreCompoundModel.

Definition at line 56 of file G4AblaInterface.hh.

                                                                                {
     return NULL;
   }

Here is the call graph for this function:

◆ DeExcite()

G4ReactionProductVector * G4AblaInterface::DeExcite ( G4Fragment & aFragment )

virtual

Implements G4VPreCompoundModel.

Definition at line 65 of file G4AblaInterface.cc.

                                                                         {
   volant->clear();
   ablaResult->clear();
 
   const G4int ARem = aFragment.GetA_asInt();
   const G4int ZRem = aFragment.GetZ_asInt();
   const G4double nuclearMass = aFragment.GetGroundStateMass() / MeV;
   const G4double eStarRem = aFragment.GetExcitationEnergy() / MeV;
   const G4double jRem = aFragment.GetAngularMomentum().mag() / hbar_Planck;
   const G4LorentzVector &pRem = aFragment.GetMomentum();
   const G4double eTotRem = pRem.e();
   const G4double eKinRem = (eTotRem - pRem.invariantMass()) / MeV;
   const G4double pxRem = pRem.x() / MeV;
   const G4double pyRem = pRem.y() / MeV;
   const G4double pzRem = pRem.z() / MeV;
 
   eventNumber++;
 
   theABLAModel->breakItUp(ARem, ZRem,
                           nuclearMass,
                           eStarRem,
                           jRem,
                           eKinRem,
                           pxRem,
                           pyRem,
                           pzRem,
                           eventNumber);
 
   G4ReactionProductVector *result = new G4ReactionProductVector;
 
   for(int j = 0; j < ablaResult->ntrack; ++j) { // Copy ABLA result to the EventInfo
     G4ReactionProduct *product = toG4Particle(ablaResult->avv[j],
                                               ablaResult->zvv[j],
                                               ablaResult->enerj[j],
                                               ablaResult->plab[j]*std::sin(ablaResult->tetlab[j]*pi/180.0)*std::cos(ablaResult->philab[j]*pi/180.0),
                                               ablaResult->plab[j]*std::sin(ablaResult->tetlab[j]*pi/180.0)*std::sin(ablaResult->philab[j]*pi/180.0),
                                               ablaResult->plab[j]*std::cos(ablaResult->tetlab[j]*pi/180.0));
     if(product)
       result->push_back(product);
   }
   return result;
 }

Here is the call graph for this function:

◆ DeExciteModelDescription()

void G4AblaInterface::DeExciteModelDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VPreCompoundModel.

Definition at line 147 of file G4AblaInterface.cc.

                                                                         {
    outFile << "ABLA V3 is a statistical model for nuclear de-excitation. It simulates\n"
      << "evaporation of neutrons, protons and alpha particles, as well as fission\n"
      << "where applicable. The code included in Geant4 is a C++ translation of the\n"
      << "original Fortran code. More details about the physics are available in the\n"
      << "the Geant4 Physics Reference Manual and in the reference articles.\n\n"
      << "References: A.R. Junghans et al., Nucl. Phys. A629 (1998) 635;\n"
      << "            J. Benlliure et al., Nucl. Phys. A628 (1998) 458.\n\n"; }

Here is the caller graph for this function:

◆ ModelDescription()

void G4AblaInterface::ModelDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4HadronicInteraction.

Definition at line 143 of file G4AblaInterface.cc.

                                                                 {
    outFile << "ABLA V3 does not provide an implementation of the ApplyYourself method!\n\n";
 }

Here is the caller graph for this function:

◆ toG4Particle()

G4ReactionProduct * G4AblaInterface::toG4Particle	(	G4int	A,
		G4int	Z,
		G4double	kinE,
		G4double	px,
		G4double	py,
		G4double	pz
	)		const

private

Convert an Abla particle to a G4ReactionProduct.

Definition at line 126 of file G4AblaInterface.cc.

                                                                                  {
   G4ParticleDefinition *def = toG4ParticleDefinition(A, Z);
   if(def == 0) { // Check if we have a valid particle definition
     return 0;
   }
   const double energy = kinE * MeV;
   const G4ThreeVector momentum(px, py, pz);
   const G4ThreeVector momentumDirection = momentum.unit();
   G4DynamicParticle p(def, momentumDirection, energy);
   G4ReactionProduct *r = new G4ReactionProduct(def);
   (*r) = p;
   return r;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ toG4ParticleDefinition()

G4ParticleDefinition * G4AblaInterface::toG4ParticleDefinition	(	G4int	A,
		G4int	Z
	)		const

private

Convert A and Z to a G4ParticleDefinition.

Definition at line 108 of file G4AblaInterface.cc.

                                                                                     {
   if     (A == 1 && Z == 1)  return G4Proton::Proton();
   else if(A == 1 && Z == 0)  return G4Neutron::Neutron();
   else if(A == 0 && Z == 1)  return G4PionPlus::PionPlus();
   else if(A == 0 && Z == -1) return G4PionMinus::PionMinus();
   else if(A == 0 && Z == 0)  return G4PionZero::PionZero();
   else if(A == 2 && Z == 1)  return G4Deuteron::Deuteron();
   else if(A == 3 && Z == 1)  return G4Triton::Triton();
   else if(A == 3 && Z == 2)  return G4He3::He3();
   else if(A == 4 && Z == 2)  return G4Alpha::Alpha();
   else if(A > 0 && Z > 0 && A >= Z) { // Returns ground state ion definition
     return G4IonTable::GetIonTable()->GetIon(Z, A, 0);
   } else { // Error, unrecognized particle
     G4cout << "Can't convert particle with A=" << A << ", Z=" << Z << " to G4ParticleDefinition, trouble ahead" << G4endl;
     return 0;
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

◆ ablaResult

G4VarNtp* G4AblaInterface::ablaResult

private

Definition at line 64 of file G4AblaInterface.hh.

◆ eventNumber

G4long G4AblaInterface::eventNumber

private

Definition at line 67 of file G4AblaInterface.hh.

◆ theABLAModel

G4Abla* G4AblaInterface::theABLAModel

private

Definition at line 66 of file G4AblaInterface.hh.

◆ volant

G4Volant* G4AblaInterface::volant

private

Definition at line 65 of file G4AblaInterface.hh.

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4AblaInterface()

◆ ~G4AblaInterface()

Member Function Documentation

◆ ApplyYourself()

◆ DeExcite()

◆ DeExciteModelDescription()

◆ ModelDescription()

◆ toG4Particle()

◆ toG4ParticleDefinition()

Member Data Documentation

◆ ablaResult

◆ eventNumber

◆ theABLAModel

◆ volant