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 ()

---

Detailed Description

Definition at line 74 of file G4PhotonEvaporation.hh.

---

Constructor & Destructor Documentation

G4PhotonEvaporation::G4PhotonEvaporation

(

)

Definition at line 72 of file G4PhotonEvaporation.cc.

References DBL_MAX, G4DiscreteGammaDeexcitation::SetICM(), and G4VGammaDeexcitation::SetTimeLimit().

  :_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

(

)

[virtual]

Definition at line 91 of file G4PhotonEvaporation.cc.

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

---

Member Function Documentation

G4FragmentVector * G4PhotonEvaporation::BreakItUp

(

const G4Fragment &

nucleus

)

[virtual]

Definition at line 244 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::DoChain(), G4cout, G4endl, G4VGammaDeexcitation::GetEO(), G4VGammaDeexcitation::GetVacantSN(), G4VGammaDeexcitation::SetNucleus(), and G4VGammaDeexcitation::SetVaccantSN().

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

{
  // 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;
}

G4FragmentVector * G4PhotonEvaporation::BreakUp

(

const G4Fragment &

nucleus

)

[virtual]

Implements G4VEvaporationChannel.

Definition at line 187 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::DoTransition(), G4cout, G4endl, G4VGammaDeexcitation::GetEO(), G4VGammaDeexcitation::GetVacantSN(), G4VGammaDeexcitation::SetNucleus(), and G4VGammaDeexcitation::SetVaccantSN().

Referenced by G4NuclearDecayChannel::DecayIt().

{
  //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;
}

G4FragmentVector * G4PhotonEvaporation::BreakUpFragment

(

G4Fragment *

theNucleus

)

[virtual]

Reimplemented from G4VEvaporationChannel.

Definition at line 142 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::DoChain(), G4cout, G4endl, G4VGammaDeexcitation::GetEO(), G4VGammaDeexcitation::GetVacantSN(), G4VGammaDeexcitation::SetNucleus(), and G4VGammaDeexcitation::SetVaccantSN().

Referenced by G4NeutronRadCapture::ApplyYourself().

{
  //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;
}

G4Fragment * G4PhotonEvaporation::EmittedFragment

(

G4Fragment *

theNucleus

)

[virtual]

Reimplemented from G4VEvaporationChannel.

Definition at line 98 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::CanDoTransition(), G4cout, G4endl, G4VGammaDeexcitation::GenerateGamma(), G4VGammaDeexcitation::Initialize(), and G4VGammaDeexcitation::SetNucleus().

{
  //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;
}

G4double G4PhotonEvaporation::GetEmissionProbability

(

G4Fragment *

theNucleus

)

[virtual]

Implements G4VEvaporationChannel.

Definition at line 293 of file G4PhotonEvaporation.cc.

References G4VEmissionProbability::EmissionProbability(), and G4Fragment::GetExcitationEnergy().

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

G4ElectronOccupancy G4PhotonEvaporation::GetEOccupancy

(

)

[inline]

Definition at line 106 of file G4PhotonEvaporation.hh.

{ return _eOccupancy;} ;

G4int G4PhotonEvaporation::GetVacantShellNumber

(

)

[inline]

Definition at line 108 of file G4PhotonEvaporation.hh.

Referenced by G4NuclearDecayChannel::DecayIt().

{ return _vShellNumber;};

void G4PhotonEvaporation::RDMForced

(

G4bool

)

Definition at line 337 of file G4PhotonEvaporation.cc.

Referenced by G4NuclearDecayChannel::DecayIt().

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

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;
}

void G4PhotonEvaporation::SetEOccupancy

(

G4ElectronOccupancy

eOccupancy

)

Definition at line 342 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::SetEO().

{
  _discrDeexcitation->SetEO(eo);
}

void G4PhotonEvaporation::SetICM

(

G4bool

)

Definition at line 322 of file G4PhotonEvaporation.cc.

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

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

void G4PhotonEvaporation::SetMaxHalfLife

(

G4double

)

Definition at line 327 of file G4PhotonEvaporation.cc.

Referenced by G4NuclearDecayChannel::DecayIt().

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

void G4PhotonEvaporation::SetTimeLimit

(

G4double

value

)

Definition at line 332 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::SetTimeLimit().

{
 _discrDeexcitation->SetTimeLimit(val);
}

void G4PhotonEvaporation::SetVerboseLevel

(

G4int

verbose

)

Definition at line 315 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::SetVerboseLevel().

Referenced by G4NuclearDecayChannel::DecayIt().

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

---

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

• G4PhotonEvaporation.hh
• G4PhotonEvaporation.cc

---