G4PreCompoundModel Class Reference

#include <G4PreCompoundModel.hh>

Inheritance diagram for G4PreCompoundModel:

Public Member Functions
	G4PreCompoundModel (G4ExcitationHandler *ptr=nullptr)
	~G4PreCompoundModel () override
G4HadFinalState *	ApplyYourself (const G4HadProjectile &thePrimary, G4Nucleus &theNucleus) override
G4ReactionProductVector *	DeExcite (G4Fragment &aFragment) override
void	BuildPhysicsTable (const G4ParticleDefinition &) override
void	InitialiseModel () override
void	ModelDescription (std::ostream &outFile) const override
void	DeExciteModelDescription (std::ostream &outFile) const override
	G4PreCompoundModel (const G4PreCompoundModel &)=delete
const G4PreCompoundModel &	operator= (const G4PreCompoundModel &right)=delete
G4bool	operator== (const G4PreCompoundModel &right) const =delete
G4bool	operator!= (const G4PreCompoundModel &right) const =delete
Public Member Functions inherited from G4VPreCompoundModel
	G4VPreCompoundModel (G4ExcitationHandler *ptr=nullptr, const G4String &modelName="PrecompoundModel")
virtual	~G4VPreCompoundModel ()
void	SetExcitationHandler (G4ExcitationHandler *ptr)
G4ExcitationHandler *	GetExcitationHandler () const
	G4VPreCompoundModel (const G4VPreCompoundModel &)=delete
const G4VPreCompoundModel &	operator= (const G4VPreCompoundModel &right)=delete
G4bool	operator== (const G4VPreCompoundModel &right) const =delete
G4bool	operator!= (const G4VPreCompoundModel &right) const =delete
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 &aTrack, G4Nucleus &targetNucleus)
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)
G4int	GetVerboseLevel () const
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
virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const
virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
	G4HadronicInteraction (const G4HadronicInteraction &right)=delete
const G4HadronicInteraction &	operator= (const G4HadronicInteraction &right)=delete
G4bool	operator== (const G4HadronicInteraction &right) const =delete
G4bool	operator!= (const G4HadronicInteraction &right) const =delete

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 64 of file G4PreCompoundModel.hh.

Constructor & Destructor Documentation

G4PreCompoundModel() [1/2]

G4PreCompoundModel::G4PreCompoundModel ( G4ExcitationHandler * ptr = nullptr )

explicit

Definition at line 74 of file G4PreCompoundModel.cc.

  : G4VPreCompoundModel(ptr,"PRECO")
{
  //G4cout << "### NEW PrecompoundModel " << this << G4endl;
  if(nullptr == ptr) { SetExcitationHandler(new G4ExcitationHandler()); }
 
  fNuclData = G4NuclearLevelData::GetInstance();
  proton = G4Proton::Proton();
  neutron = G4Neutron::Neutron();
  modelID = G4PhysicsModelCatalog::GetModelID("model_PRECO");
}

Referenced by G4PreCompoundModel(), operator!=(), operator=(), and operator==().

~G4PreCompoundModel()

G4PreCompoundModel::~G4PreCompoundModel ( )

override

Definition at line 88 of file G4PreCompoundModel.cc.

{
  delete theEmission;
  delete theTransition;
  delete GetExcitationHandler();
  delete theMultiFrag;
  theResult.Clear();
}

G4PreCompoundModel() [2/2]

G4PreCompoundModel::G4PreCompoundModel ( const G4PreCompoundModel & )

delete

Member Function Documentation

ApplyYourself()

G4HadFinalState * G4PreCompoundModel::ApplyYourself ( const G4HadProjectile & thePrimary, G4Nucleus & theNucleus )

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 149 of file G4PreCompoundModel.cc.

{
  const G4ParticleDefinition* primary = thePrimary.GetDefinition();
  if(primary != neutron && primary != proton) {
    G4ExceptionDescription ed;
    ed << "G4PreCompoundModel is used for ";
    if(primary) { ed << primary->GetParticleName(); }
    G4Exception("G4PreCompoundModel::ApplyYourself()","had0033",FatalException,
        ed,"");
    return nullptr;
  }
 
  G4int Zp = 0;
  G4int Ap = 1;
  if(primary == proton) { Zp = 1; }
 
  G4double timePrimary=thePrimary.GetGlobalTime();
  G4int A = theNucleus.GetA_asInt();
  G4int Z = theNucleus.GetZ_asInt();
 
  //G4cout << "### G4PreCompoundModel::ApplyYourself: A= " << A << " Z= " << Z
  //     << " Ap= " << Ap << " Zp= " << Zp << G4endl;
 
  // 4-Momentum
  G4LorentzVector p = thePrimary.Get4Momentum();
  G4double mass = G4NucleiProperties::GetNuclearMass(A, Z);
  p += G4LorentzVector(0.0,0.0,0.0,mass);
  //G4cout << "Primary 4-mom " << p << "  mass= " << mass << G4endl;
 
  // prepare fragment
  G4Fragment anInitialState(A + Ap, Z + Zp, p);
  anInitialState.SetNumberOfExcitedParticle(2, 1);
  anInitialState.SetNumberOfHoles(1,0);
  anInitialState.SetCreationTime(thePrimary.GetGlobalTime());
  anInitialState.SetCreatorModelID(modelID);
 
  // call excitation handler
  G4ReactionProductVector* result = DeExcite(anInitialState);
 
  // fill particle change
  theResult.Clear();
  theResult.SetStatusChange(stopAndKill);
  for(auto const & prod : *result) {
    G4DynamicParticle * aNewDP = new G4DynamicParticle(prod->GetDefinition(),
                                                       prod->GetTotalEnergy(),
                                                       prod->GetMomentum());
    G4HadSecondary aNew = G4HadSecondary(aNewDP);
    G4double time = std::max(prod->GetFormationTime(), 0.0);
    aNew.SetTime(timePrimary + time);
    aNew.SetCreatorModelID(prod->GetCreatorModelID());
    delete prod;
    theResult.AddSecondary(aNew);
  }
  delete result;
  return &theResult;
}

BuildPhysicsTable()

void G4PreCompoundModel::BuildPhysicsTable ( const G4ParticleDefinition & )

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 99 of file G4PreCompoundModel.cc.

{
  InitialiseModel();
}

DeExcite()

G4ReactionProductVector * G4PreCompoundModel::DeExcite ( G4Fragment & aFragment )

overridevirtual

Implements G4VPreCompoundModel.

Definition at line 209 of file G4PreCompoundModel.cc.

{
  if (!isInitialised) { InitialiseModel(); }
  if (usePrecoInterface) { return fInterface->DeExcite(aFragment); }
 
  G4ReactionProductVector* result = new G4ReactionProductVector();
 
  G4int Z = aFragment.GetZ_asInt();
  G4int A = aFragment.GetA_asInt();
  G4double U = aFragment.GetExcitationEnergy();
 
  if (1 < fVerbose) {
    G4cout << "### G4PreCompoundModel::DeExcite Z=" << Z << " A=" << A
       << " U(MeV)=" << U << " Umin=" << fLowLimitExc*A << G4endl;
    G4cout << aFragment << G4endl;
  }  
  
  // Conditions to skip pre-compound and perform equilibrium emission 
  if (!isActive || (Z < minZ && A < minA) || U < fLowLimitExc*A ||
      U > fHighLimitExc*A || 0 < aFragment.GetNumberOfLambdas() ||
      0 == Z || Z >= A) {
    PerformEquilibriumEmission(aFragment, result);
    return result;
  }
  
  // apply multi-fragmentation above threshold
  G4FragmentVector* theTempResult;
  if (U > A*fMinEForMultiFrag) {
    theTempResult = theMultiFrag->BreakItUp(aFragment);
  }
  else {
    theTempResult = new G4FragmentVector();
    theTempResult->push_back(new G4Fragment(aFragment));
  }
 
  // apply pre-compound on the list of fragments 
  for (auto & ptr : *theTempResult) {
    DoIt(result, *ptr);
    delete ptr;
  }
  delete theTempResult;
  return result;
}

Referenced by ApplyYourself().

DeExciteModelDescription()

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

overridevirtual

Implements G4VPreCompoundModel.

Definition at line 384 of file G4PreCompoundModel.cc.

{
  outFile << "description of precompound model as used with DeExcite()" << "\n";
}

InitialiseModel()

void G4PreCompoundModel::InitialiseModel ( )

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 106 of file G4PreCompoundModel.cc.

{
  if (isInitialised) { return; }
  isInitialised = true;
 
  G4DeexPrecoParameters* param = fNuclData->GetParameters();
  fVerbose = param->GetVerbose();
 
  if (param->PrecoDummy() || eDeexcitation == param->GetPreCompoundType()) {
    isActive = false;
  }
  else if (ePrecoInterface == param->GetPreCompoundType()) {
    usePrecoInterface = true;
    fInterface = new G4PreCompoundInterface();
  }
  else {
    theEmission = new G4PreCompoundEmission();
    if (param->UseHETC()) { theEmission->SetHETCModel(); }
    theEmission->SetOPTxs(param->GetPrecoModelType());
 
    if (param->UseGNASH()) { theTransition = new G4GNASHTransitions; }
    else { theTransition = new G4PreCompoundTransitions(); }
    theTransition->UseNGB(param->NeverGoBack());
    theTransition->UseCEMtr(param->UseCEM());
  }
  if (isActive) {
    theMultiFrag = new G4StatMF(); 
    fLowLimitExc = param->GetPrecoLowEnergy();
    fHighLimitExc = param->GetPrecoHighEnergy();
    fMinEForMultiFrag = param->GetMinExPerNucleounForMF();
 
    useSCO = param->UseSoftCutoff();
 
    minZ = param->GetMinZForPreco();
    minA = param->GetMinAForPreco();
  }
  
  GetExcitationHandler()->Initialise();
}

Referenced by BuildPhysicsTable(), and DeExcite().

ModelDescription()

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

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 358 of file G4PreCompoundModel.cc.

{
  outFile 
    << "The GEANT4 precompound model is considered as an extension of the\n"
    <<  "hadron kinetic model. It gives a possibility to extend the low energy range\n"
    <<  "of the hadron kinetic model for nucleon-nucleus inelastic collision and it \n"
    <<  "provides a ”smooth” transition from kinetic stage of reaction described by the\n"
    <<  "hadron kinetic model to the equilibrium stage of reaction described by the\n"
    <<  "equilibrium deexcitation models.\n"
    <<  "The initial information for calculation of pre-compound nuclear stage\n"
    <<  "consists of the atomic mass number A, charge Z of residual nucleus, its\n"
    <<  "four momentum P0, excitation energy U and number of excitons n, which equals\n"
    <<  "the sum of the number of particles p (from them p_Z are charged) and the number of\n"
    <<  "holes h.\n"
    <<  "At the preequilibrium stage of reaction, we follow the exciton model approach in ref. [1],\n"
    <<  "taking into account the competition among all possible nuclear transitions\n"
    <<  "with ∆n = +2, −2, 0 (which are defined by their associated transition probabilities) and\n"
    <<  "the emission of neutrons, protons, deuterons, thritium and helium nuclei (also defined by\n"
    <<  "their associated emission  probabilities according to exciton model)\n"
    <<  "\n"
    <<  "[1] K.K. Gudima, S.G. Mashnik, V.D. Toneev, Nucl. Phys. A401 329 (1983)\n"
    <<  "\n";
}

operator!=()

G4bool G4PreCompoundModel::operator!= ( const G4PreCompoundModel & right ) const

delete

operator=()

const G4PreCompoundModel & G4PreCompoundModel::operator= ( const G4PreCompoundModel & right )

delete

operator==()

G4bool G4PreCompoundModel::operator== ( const G4PreCompoundModel & right ) const

delete

---

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

• G4PreCompoundModel.hh
• G4PreCompoundModel.cc