G4PhotonEvaporation Class Reference

#include <G4PhotonEvaporation.hh>

Inheritance diagram for G4PhotonEvaporation:

Public Member Functions
	G4PhotonEvaporation ()
virtual	~G4PhotonEvaporation ()
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakItUp (const G4Fragment &nucleus)
virtual G4FragmentVector *	BreakUp (const G4Fragment &nucleus)
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)
virtual void	SetEmissionStrategy (G4VEmissionProbability *probAlgorithm)
void	SetVerboseLevel (G4int verbose)
void	SetICM (G4bool)
void	RDMForced (G4bool)
void	SetMaxHalfLife (G4double)
void	SetTimeLimit (G4double value)
void	SetEOccupancy (G4ElectronOccupancy eOccupancy)
G4ElectronOccupancy	GetEOccupancy ()
G4int	GetVacantShellNumber ()
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="Anonymous")
virtual	~G4VEvaporationChannel ()
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)=0
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)=0
G4String	GetName () const
void	SetName (const G4String &aName)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Additional Inherited Members
Protected Attributes inherited from G4VEvaporationChannel
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 74 of file G4PhotonEvaporation.hh.

Constructor & Destructor Documentation

G4PhotonEvaporation()

G4PhotonEvaporation::G4PhotonEvaporation

(

)

Definition at line 72 of file G4PhotonEvaporation.cc.

  :_verbose(0),_myOwnProbAlgorithm (true),
   _eOccupancy(0), _vShellNumber(-1),_gammaE(0.)
{ 
  _probAlgorithm = new G4E1Probability;
 
  G4double timeLimit = DBL_MAX;
  char* env = getenv("G4AddTimeLimitToPhotonEvaporation"); 
  if(env) { timeLimit = 1.e-16*second; }
 
  G4DiscreteGammaDeexcitation* p = new G4DiscreteGammaDeexcitation();
  p->SetICM(false);
  p->SetTimeLimit(timeLimit);
 
  _discrDeexcitation = p;
  _contDeexcitation = new G4ContinuumGammaDeexcitation;
  _nucleus = 0;
}

~G4PhotonEvaporation()

G4PhotonEvaporation::~G4PhotonEvaporation ( )

virtual

Definition at line 91 of file G4PhotonEvaporation.cc.

{ 
  if(_myOwnProbAlgorithm) delete _probAlgorithm;
  delete _discrDeexcitation;
  delete _contDeexcitation;
}

Member Function Documentation

BreakItUp()

G4FragmentVector * G4PhotonEvaporation::BreakItUp ( const G4Fragment &  nucleus )

virtual

Definition at line 244 of file G4PhotonEvaporation.cc.

{
  // The same pointer of primary nucleus
  _nucleus = new G4Fragment(nucleus);
  _contDeexcitation->SetNucleus(_nucleus);
  _discrDeexcitation->SetNucleus(_nucleus);
 
  //G4cout << "G4PhotonEvaporation::BreakItUp:  " << nucleus << G4endl;
 
  // Do the whole gamma chain 
  G4FragmentVector* products = _contDeexcitation->DoChain();  
  if( !products ) { products = new G4FragmentVector; }
 
  // Products from continuum gamma transitions
  if (_verbose > 0) {
    G4cout << " = BreakItUp = " << products->size()
       << " gammas from ContinuumDeexcitation " << G4endl;
  }
 
  // Products from discrete gamma transitions
  G4FragmentVector* discrProducts = _discrDeexcitation->DoChain();
  if(discrProducts) {
    _eOccupancy = _discrDeexcitation->GetEO();
    _vShellNumber = _discrDeexcitation->GetVacantSN();
 
    // not sure if the following line is needed!
    _discrDeexcitation->SetVaccantSN(-1);
 
    if (_verbose > 0) {
      G4cout << " = BreakItUp = " << discrProducts->size() 
         << " gammas from DiscreteDeexcitation " << G4endl;
    }
    G4FragmentVector::iterator i;
    for (i = discrProducts->begin(); i != discrProducts->end(); ++i)
      {
    products->push_back(*i);
      }
    delete discrProducts;
  }
  // Add deexcited nucleus to products
  products->push_back(_nucleus);
 
  if (_verbose > 0) {
    G4cout << "*-*-* Photon evaporation: " << products->size() << G4endl;
  }
  return products;
}

Referenced by G4NeutronHPCaptureFS::ApplyYourself(), and G4NuclearDecayChannel::DecayIt().

BreakUp()

G4FragmentVector * G4PhotonEvaporation::BreakUp ( const G4Fragment &  nucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 187 of file G4PhotonEvaporation.cc.

{
  //G4cout << "G4PhotonEvaporation::BreakUp" << G4endl;
  _nucleus = new G4Fragment(nucleus);
 
  _contDeexcitation->SetNucleus(_nucleus);
  _discrDeexcitation->SetNucleus(_nucleus);
  
  // Do one photon emission
 
  // Products from continuum gamma transitions
 
  G4FragmentVector* products = _contDeexcitation->DoTransition();  
  if( !products ) { products = new G4FragmentVector(); }
  else if(_verbose > 0) {
    G4cout << "G4PhotonEvaporation::BreakUp " << products->size() 
       << " gammas from ContinuesDeexcitation " << G4endl;
    G4cout << "   Residual: " << nucleus << G4endl;
  }
 
  if (0 == products->size())
    {
      // Products from discrete gamma transitions
      G4FragmentVector* discrProducts = _discrDeexcitation->DoTransition();
 
      if (discrProducts) {
    _eOccupancy = _discrDeexcitation->GetEO();
    _vShellNumber = _discrDeexcitation->GetVacantSN();
 
    // not sure if the following line is needed!
    _discrDeexcitation->SetVaccantSN(-1);
    //
    if (_verbose > 0) {
      G4cout << " = BreakUp = " << discrProducts->size() 
         << " gammas from DiscreteDeexcitation " 
         << G4endl;
      G4cout << "   Residual: " << nucleus << G4endl;
    }
    G4FragmentVector::iterator i;
    for (i = discrProducts->begin(); i != discrProducts->end(); ++i)
      {
        products->push_back(*i);
      }
    delete discrProducts;
      }
    }
  
  // Add deexcited nucleus to products
  products->push_back(_nucleus);
 
  if (_verbose > 0) {
    G4cout << "*-*-*-* Photon evaporation: " << products->size() << G4endl;
  }
 
  return products;
}

Referenced by G4NuclearDecayChannel::DecayIt().

BreakUpFragment()

G4FragmentVector * G4PhotonEvaporation::BreakUpFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 142 of file G4PhotonEvaporation.cc.

{
  //G4cout << "G4PhotonEvaporation::BreakUpFragment" << G4endl;
  // The same pointer of primary nucleus
  _nucleus = nucleus;
  _contDeexcitation->SetNucleus(_nucleus);
  _discrDeexcitation->SetNucleus(_nucleus);
 
  // Do the whole gamma chain 
  G4FragmentVector* products = _contDeexcitation->DoChain();  
  if( !products ) { products = new G4FragmentVector(); }
 
  if (_verbose > 0) {
    G4cout << "G4PhotonEvaporation::BreakUpFragment " << products->size() 
       << " gammas from ContinuumDeexcitation " << G4endl;
    G4cout << "   Residual: " << nucleus << G4endl;
  }
  // Products from discrete gamma transitions
  G4FragmentVector* discrProducts = _discrDeexcitation->DoChain();
  if(discrProducts) {
    _eOccupancy = _discrDeexcitation->GetEO();
    _vShellNumber = _discrDeexcitation->GetVacantSN();
 
    // not sure if the following line is needed!
    _discrDeexcitation->SetVaccantSN(-1);
 
    if (_verbose > 0) {
      G4cout << "G4PhotonEvaporation::BreakUpFragment " << discrProducts->size() 
         << " gammas from DiscreteDeexcitation " << G4endl;
      G4cout << "   Residual: " << nucleus << G4endl;
    }
    G4FragmentVector::iterator i;
    for (i = discrProducts->begin(); i != discrProducts->end(); ++i)
      {
    products->push_back(*i);
      }
    delete discrProducts;
  }
 
  if (_verbose > 0) {
    G4cout << "*-*-* Photon evaporation: " << products->size() << G4endl;
  }
  return products;
}

Referenced by G4NeutronRadCapture::ApplyYourself().

EmittedFragment()

G4Fragment * G4PhotonEvaporation::EmittedFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 98 of file G4PhotonEvaporation.cc.

{
  //G4cout << "G4PhotonEvaporation::EmittedFragment" << G4endl;
  _nucleus = nucleus;
  
  // Do one photon emission by the continues deexcitation  
  _contDeexcitation->SetNucleus(_nucleus);
  _contDeexcitation->Initialize();
 
  if(_contDeexcitation->CanDoTransition()) {  
    G4Fragment* gamma = _contDeexcitation->GenerateGamma();
    if(gamma) { 
      if (_verbose > 0) {
    G4cout << "G4PhotonEvaporation::EmittedFragment continium deex: "   
           << gamma << G4endl;
        G4cout << "   Residual: " << nucleus << G4endl;
      }
      return gamma; 
    }
  }
 
  // Do one photon emission by the discrete deexcitation 
  _discrDeexcitation->SetNucleus(_nucleus);
  _discrDeexcitation->Initialize();
 
  if(_discrDeexcitation->CanDoTransition()) {  
    G4Fragment* gamma = _discrDeexcitation->GenerateGamma();
    if(gamma) { 
      if (_verbose > 0) {
    G4cout << "G4PhotonEvaporation::EmittedFragment discrete deex: "   
           << gamma << G4endl;
        G4cout << "   Residual: " << nucleus << G4endl;
      }
      return gamma; 
    }
  }
 
  if (_verbose > 0) {
    G4cout << "G4PhotonEvaporation unable emit gamma: " 
       << nucleus << G4endl;
  }
  return 0;
}

GetEmissionProbability()

G4double G4PhotonEvaporation::GetEmissionProbability ( G4Fragment *  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 293 of file G4PhotonEvaporation.cc.

{
  _nucleus = theNucleus;
  G4double prob = 
    _probAlgorithm->EmissionProbability(*_nucleus,_nucleus->GetExcitationEnergy());
  return prob;
}

GetEOccupancy()

G4ElectronOccupancy G4PhotonEvaporation::GetEOccupancy ( )

inline

Definition at line 106 of file G4PhotonEvaporation.hh.

{ return _eOccupancy;} ;

GetVacantShellNumber()

G4int G4PhotonEvaporation::GetVacantShellNumber ( )

inline

Definition at line 108 of file G4PhotonEvaporation.hh.

{ return _vShellNumber;};

Referenced by G4NuclearDecayChannel::DecayIt().

RDMForced()

void G4PhotonEvaporation::RDMForced

(

G4bool

fromRDM

)

Definition at line 337 of file G4PhotonEvaporation.cc.

{
 (static_cast <G4DiscreteGammaDeexcitation*> (_discrDeexcitation))->SetRDM(fromRDM);
}

Referenced by G4NuclearDecayChannel::DecayIt().

SetEmissionStrategy()

void G4PhotonEvaporation::SetEmissionStrategy ( G4VEmissionProbability *  probAlgorithm )

virtual

Definition at line 302 of file G4PhotonEvaporation.cc.

{
 
  // CD - not sure about always wanting to delete this pointer....
 
  if(_myOwnProbAlgorithm) delete _probAlgorithm;
 
  _probAlgorithm = probAlgorithm;
 
  _myOwnProbAlgorithm = false;
}

SetEOccupancy()

void G4PhotonEvaporation::SetEOccupancy

(

G4ElectronOccupancy

eOccupancy

)

Definition at line 342 of file G4PhotonEvaporation.cc.

{
  _discrDeexcitation->SetEO(eo);
}

SetICM()

void G4PhotonEvaporation::SetICM

(

G4bool

ic

)

Definition at line 322 of file G4PhotonEvaporation.cc.

{
 (static_cast <G4DiscreteGammaDeexcitation*> (_discrDeexcitation))->SetICM(ic);
}

Referenced by G4NeutronHPCaptureFS::ApplyYourself(), and G4NuclearDecayChannel::DecayIt().

SetMaxHalfLife()

void G4PhotonEvaporation::SetMaxHalfLife

(

G4double

hl

)

Definition at line 327 of file G4PhotonEvaporation.cc.

{
 (static_cast <G4DiscreteGammaDeexcitation*> (_discrDeexcitation))->SetHL(hl);
}

Referenced by G4NuclearDecayChannel::DecayIt().

SetTimeLimit()

void G4PhotonEvaporation::SetTimeLimit

(

G4double

value

)

Definition at line 332 of file G4PhotonEvaporation.cc.

{
 _discrDeexcitation->SetTimeLimit(val);
}

SetVerboseLevel()

void G4PhotonEvaporation::SetVerboseLevel

(

G4int

verbose

)

Definition at line 315 of file G4PhotonEvaporation.cc.

{
  _verbose = verbose;
  _contDeexcitation->SetVerboseLevel(verbose);
  _discrDeexcitation->SetVerboseLevel(verbose);
}

Referenced by G4NuclearDecayChannel::DecayIt().

---

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

• G4PhotonEvaporation.hh
• G4PhotonEvaporation.cc