G4PreCompoundModel Class Reference

#include <G4PreCompoundModel.hh>

Inheritance diagram for G4PreCompoundModel:

Collaboration diagram for G4PreCompoundModel:

Public Member Functions
	G4PreCompoundModel (G4ExcitationHandler *ptr=0)
virtual	~G4PreCompoundModel ()
virtual G4HadFinalState *	ApplyYourself (const G4HadProjectile &thePrimary, G4Nucleus &theNucleus)
virtual G4ReactionProductVector *	DeExcite (G4Fragment &aFragment)
virtual void	ModelDescription (std::ostream &outFile) const
virtual void	DeExciteModelDescription (std::ostream &outFile) const
void	UseHETCEmission ()
void	UseDefaultEmission ()
void	UseGNASHTransition ()
void	UseDefaultTransition ()
void	SetOPTxs (G4int opt)
void	UseSICB ()
void	UseNGB ()
void	UseSCO ()
void	UseCEMtr ()
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
void	PerformEquilibriumEmission (const G4Fragment &aFragment, G4ReactionProductVector *theResult) const
	G4PreCompoundModel (const G4PreCompoundModel &)
const G4PreCompoundModel &	operator= (const G4PreCompoundModel &right)
G4bool	operator== (const G4PreCompoundModel &right) const
G4bool	operator!= (const G4PreCompoundModel &right) const

Private Attributes
G4PreCompoundEmission *	theEmission
G4VPreCompoundTransitions *	theTransition
const G4ParticleDefinition *	proton
const G4ParticleDefinition *	neutron
G4double	fLevelDensity
G4bool	useHETCEmission
G4bool	useGNASHTransition
G4int	OPTxs
G4bool	useSICB
G4bool	useNGB
G4bool	useSCO
G4bool	useCEMtr
G4int	maxZ
G4int	maxA
G4HadFinalState	theResult

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 = 0

)

Definition at line 70 of file G4PreCompoundModel.cc.

  : G4VPreCompoundModel(ptr,"PRECO"), useHETCEmission(false), 
    useGNASHTransition(false), OPTxs(3), useSICB(false), 
    useNGB(false), useSCO(false), useCEMtr(true), maxZ(3), maxA(5) 
{
  if(!ptr) { SetExcitationHandler(new G4ExcitationHandler()); }
  G4PreCompoundParameters param;

  // 12/pi2 factor is used in real computation
  fLevelDensity = param.GetLevelDensity()*12.0/CLHEP::pi2;

  theEmission = new G4PreCompoundEmission();
  if(useHETCEmission) { theEmission->SetHETCModel(); }
  else { theEmission->SetDefaultModel(); }
  theEmission->SetOPTxs(OPTxs);
  theEmission->UseSICB(useSICB);

  if(useGNASHTransition) { theTransition = new G4GNASHTransitions; }
  else { theTransition = new G4PreCompoundTransitions(); }
  theTransition->UseNGB(useNGB);
  theTransition->UseCEMtr(useCEMtr);

  proton = G4Proton::Proton();
  neutron = G4Neutron::Neutron();
}

Here is the call graph for this function:

~G4PreCompoundModel()

G4PreCompoundModel::~G4PreCompoundModel ( )

virtual

Definition at line 98 of file G4PreCompoundModel.cc.

{
  delete theEmission;
  delete theTransition;
  delete GetExcitationHandler();
}

Here is the call graph for this function:

G4PreCompoundModel() [2/2]

G4PreCompoundModel::G4PreCompoundModel ( const G4PreCompoundModel &  )

private

Member Function Documentation

ApplyYourself()

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

virtual

Implements G4VPreCompoundModel.

Definition at line 108 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 0;
  }

  G4int Zp = 0;
  G4int Ap = 1;
  if(primary == proton) { Zp = 1; }

  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());
  
  // call excitation handler
  G4ReactionProductVector * result = DeExcite(anInitialState);

  // fill particle change
  theResult.Clear();
  theResult.SetStatusChange(stopAndKill);
  for(G4ReactionProductVector::iterator i= result->begin(); 
      i != result->end(); ++i)
    {
      G4DynamicParticle * aNew = 
    new G4DynamicParticle((*i)->GetDefinition(),
                  (*i)->GetTotalEnergy(),
                  (*i)->GetMomentum());
      delete (*i);
      theResult.AddSecondary(aNew);
    }
  delete result;
  
  //return the filled particle change
  return &theResult;
}

Here is the call graph for this function:

DeExcite()

G4ReactionProductVector * G4PreCompoundModel::DeExcite ( G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundModel.

Definition at line 166 of file G4PreCompoundModel.cc.

{
  G4ReactionProductVector * Result = new G4ReactionProductVector;
  G4double Eex = aFragment.GetExcitationEnergy();
  G4int Z = aFragment.GetZ_asInt(); 
  G4int A = aFragment.GetA_asInt(); 

  //G4cout << "### G4PreCompoundModel::DeExcite" << G4endl;
  //G4cout << aFragment << G4endl;
 
  // Perform Equilibrium Emission 
  if ((Z < maxZ && A < maxA) || Eex < MeV /*|| Eex > 3.*MeV*A*/) {
    PerformEquilibriumEmission(aFragment, Result);
    return Result;
  }
  
  // main loop  
  G4int count = 0;
  const G4int countmax = 1000;
  for (;;) {
    //G4cout << "### PreCompound loop over fragment" << G4endl;
    //G4cout << aFragment << G4endl;
    G4int EquilibriumExcitonNumber = 
      G4lrint(std::sqrt(aFragment.GetExcitationEnergy()
            *aFragment.GetA_asInt()*fLevelDensity));
    //   
    //    G4cout<<"Neq="<<EquilibriumExcitonNumber<<G4endl;
    //
    // J. M. Quesada (Jan. 08)  equilibrium hole number could be used as preeq.
    // evap. delimiter (IAEA report)
    
    // Loop for transitions, it is performed while there are 
    // preequilibrium transitions.
    G4bool ThereIsTransition = false;
    
    //        G4cout<<"----------------------------------------"<<G4endl;
    //        G4double NP=aFragment.GetNumberOfParticles();
    //        G4double NH=aFragment.GetNumberOfHoles();
    //        G4double NE=aFragment.GetNumberOfExcitons();
    //        G4cout<<" Ex. Energy="<<aFragment.GetExcitationEnergy()<<G4endl;
    //   G4cout<<"N. excitons="<<NE<<"  N. Part="<<NP<<"N. Holes ="<<NH<<G4endl;
    do {
      ++count;
      //G4cout<<"transition number .."<<count
      //      <<" n ="<<aFragment.GetNumberOfExcitons()<<G4endl;
      G4bool go_ahead = false;
      // soft cutoff criterium as an "ad-hoc" solution to force 
      // increase in  evaporation  
      G4int test = aFragment.GetNumberOfExcitons();
      if (test <= EquilibriumExcitonNumber) { go_ahead=true; }

      //J. M. Quesada (Apr. 08): soft-cutoff switched off by default
      if (useSCO && go_ahead)
    {
      G4double x = G4double(test)/G4double(EquilibriumExcitonNumber) - 1;
      if( G4UniformRand() < 1.0 -  G4Exp(-x*x/0.32) ) { go_ahead = false; }
    } 
        
      // JMQ: WARNING:  CalculateProbability MUST be called prior to Get!! 
      // (O values would be returned otherwise)
      G4double TotalTransitionProbability = 
    theTransition->CalculateProbability(aFragment);
      G4double P1 = theTransition->GetTransitionProb1();
      G4double P2 = theTransition->GetTransitionProb2();
      G4double P3 = theTransition->GetTransitionProb3();
      //G4cout<<"#0 P1="<<P1<<" P2="<<P2<<"  P3="<<P3<<G4endl;
      
      //J.M. Quesada (May 2008) Physical criterium (lamdas)  PREVAILS over 
      //                        approximation (critical exciton number)
      //V.Ivanchenko (May 2011) added check on number of nucleons
      //                        to send a fragment to FermiBreakUp
      if(!go_ahead || P1 <= P2+P3 || 
     (aFragment.GetZ_asInt() < maxZ && aFragment.GetA_asInt() < maxA) )
    {
      //G4cout<<"#4 EquilibriumEmission"<<G4endl; 
      PerformEquilibriumEmission(aFragment,Result);
      return Result;
    }
      else 
    {
      //
      // Check if number of excitons is greater than 0
      // else perform equilibrium emission
      if (aFragment.GetNumberOfExcitons() <= 0) 
        {
          PerformEquilibriumEmission(aFragment,Result);
          return Result;
        }
        
      G4double TotalEmissionProbability = 
        theEmission->GetTotalProbability(aFragment);
      //
      //G4cout<<"#1 TotalEmissionProbability="<<TotalEmissionProbability
      // <<" Nex= " <<aFragment.GetNumberOfExcitons()<<G4endl;
      //J.M.Quesada (May 08) this has already been done in order to decide  
      //                     what to do (preeq-eq) 
      // Sum of all probabilities
      G4double TotalProbability = TotalEmissionProbability 
        + TotalTransitionProbability;
            
      // Select subprocess
      if (TotalProbability*G4UniformRand() > TotalEmissionProbability) 
        {
          //G4cout<<"#2 Transition"<<G4endl; 
          // It will be transition to state with a new number of excitons
          ThereIsTransition = true;     
          // Perform the transition
          theTransition->PerformTransition(aFragment);
        } 
      else 
        {
          //G4cout<<"#3 Emission"<<G4endl; 
          // It will be fragment emission
          ThereIsTransition = false;
          Result->push_back(theEmission->PerformEmission(aFragment));
        }
    }
      // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
    } while (ThereIsTransition);   // end of do loop

    // stop if too many iterations
    if(count >= countmax) {
      G4ExceptionDescription ed;
      ed << "G4PreCompoundModel loop over " << countmax << " iterations; "
     << "current G4Fragment: \n" << aFragment;
      G4Exception("G4PreCompoundModel::DeExcite()","had0034",JustWarning,
          ed,"");
      PerformEquilibriumEmission(aFragment, Result);
      return Result;
    }
  } // end of for (;;) loop
  return Result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DeExciteModelDescription()

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

virtual

Reimplemented from G4VPreCompoundModel.

Definition at line 386 of file G4PreCompoundModel.cc.

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

ModelDescription()

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

virtual

Reimplemented from G4HadronicInteraction.

Definition at line 363 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, deutrons, 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

private

operator=()

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

private

operator==()

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

private

PerformEquilibriumEmission()

void G4PreCompoundModel::PerformEquilibriumEmission ( const G4Fragment &  aFragment, G4ReactionProductVector *  theResult  ) const

inlineprivate

Definition at line 137 of file G4PreCompoundModel.hh.

{
  G4ReactionProductVector* theEquilibriumResult = 
    GetExcitationHandler()->BreakItUp(aFragment);
  Result->insert(Result->end(),theEquilibriumResult->begin(), theEquilibriumResult->end());
  delete theEquilibriumResult;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetOPTxs()

void G4PreCompoundModel::SetOPTxs

(

G4int

opt

)

Definition at line 332 of file G4PreCompoundModel.cc.

{ 
  OPTxs = opt; 
  theEmission->SetOPTxs(OPTxs);
}

Here is the call graph for this function:

UseCEMtr()

void G4PreCompoundModel::UseCEMtr

(

)

Definition at line 354 of file G4PreCompoundModel.cc.

{ 
  useCEMtr = true; 
}

UseDefaultEmission()

void G4PreCompoundModel::UseDefaultEmission

(

)

Definition at line 310 of file G4PreCompoundModel.cc.

{ 
  useHETCEmission = false; 
  theEmission->SetDefaultModel();
}

Here is the call graph for this function:

UseDefaultTransition()

void G4PreCompoundModel::UseDefaultTransition

(

)

Definition at line 324 of file G4PreCompoundModel.cc.

                                              { 
  useGNASHTransition = false; 
  delete theTransition;
  theTransition = new G4PreCompoundTransitions();
  theTransition->UseNGB(useNGB);
  theTransition->UseCEMtr(useCEMtr);
}

Here is the call graph for this function:

UseGNASHTransition()

void G4PreCompoundModel::UseGNASHTransition

(

)

Definition at line 316 of file G4PreCompoundModel.cc.

                                            { 
  useGNASHTransition = true; 
  delete theTransition;
  theTransition = new G4GNASHTransitions;
  theTransition->UseNGB(useNGB);
  theTransition->UseCEMtr(useCEMtr);
}

Here is the call graph for this function:

UseHETCEmission()

void G4PreCompoundModel::UseHETCEmission

(

)

Definition at line 304 of file G4PreCompoundModel.cc.

{ 
  useHETCEmission = true; 
  theEmission->SetHETCModel();
}

Here is the call graph for this function:

UseNGB()

void G4PreCompoundModel::UseNGB

(

)

Definition at line 344 of file G4PreCompoundModel.cc.

{ 
  useNGB = true; 
}

UseSCO()

void G4PreCompoundModel::UseSCO

(

)

Definition at line 349 of file G4PreCompoundModel.cc.

{ 
  useSCO = true; 
}

UseSICB()

void G4PreCompoundModel::UseSICB

(

)

Definition at line 338 of file G4PreCompoundModel.cc.

{ 
  useSICB = true; 
  theEmission->UseSICB(useSICB);
}

Here is the call graph for this function:

Member Data Documentation

fLevelDensity

G4double G4PreCompoundModel::fLevelDensity

private

Definition at line 115 of file G4PreCompoundModel.hh.

maxA

G4int G4PreCompoundModel::maxA

private

Definition at line 130 of file G4PreCompoundModel.hh.

maxZ

G4int G4PreCompoundModel::maxZ

private

Definition at line 129 of file G4PreCompoundModel.hh.

neutron

const G4ParticleDefinition* G4PreCompoundModel::neutron

private

Definition at line 113 of file G4PreCompoundModel.hh.

OPTxs

G4int G4PreCompoundModel::OPTxs

private

Definition at line 121 of file G4PreCompoundModel.hh.

proton

const G4ParticleDefinition* G4PreCompoundModel::proton

private

Definition at line 112 of file G4PreCompoundModel.hh.

theEmission

G4PreCompoundEmission* G4PreCompoundModel::theEmission

private

Definition at line 109 of file G4PreCompoundModel.hh.

theResult

G4HadFinalState G4PreCompoundModel::theResult

private

Definition at line 132 of file G4PreCompoundModel.hh.

theTransition

G4VPreCompoundTransitions* G4PreCompoundModel::theTransition

private

Definition at line 110 of file G4PreCompoundModel.hh.

useCEMtr

G4bool G4PreCompoundModel::useCEMtr

private

Definition at line 127 of file G4PreCompoundModel.hh.

useGNASHTransition

G4bool G4PreCompoundModel::useGNASHTransition

private

Definition at line 118 of file G4PreCompoundModel.hh.

useHETCEmission

G4bool G4PreCompoundModel::useHETCEmission

private

Definition at line 117 of file G4PreCompoundModel.hh.

useNGB

G4bool G4PreCompoundModel::useNGB

private

Definition at line 125 of file G4PreCompoundModel.hh.

useSCO

G4bool G4PreCompoundModel::useSCO

private

Definition at line 126 of file G4PreCompoundModel.hh.

useSICB

G4bool G4PreCompoundModel::useSICB

private

Definition at line 124 of file G4PreCompoundModel.hh.

---

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

• G4PreCompoundModel.hh
• G4PreCompoundModel.cc