G4DiscreteGammaDeexcitation Class Reference

#include <G4DiscreteGammaDeexcitation.hh>

Inheritance diagram for G4DiscreteGammaDeexcitation:

Public Member Functions
	G4DiscreteGammaDeexcitation ()
virtual	~G4DiscreteGammaDeexcitation ()
virtual G4VGammaTransition *	CreateTransition ()
virtual G4bool	CanDoTransition ()
void	SetICM (G4bool hl)
void	SetRDM (G4bool hl)
void	SetHL (G4double hl)
Public Member Functions inherited from G4VGammaDeexcitation
	G4VGammaDeexcitation ()
virtual	~G4VGammaDeexcitation ()
G4FragmentVector *	DoTransition ()
G4FragmentVector *	DoChain ()
G4Fragment *	GenerateGamma ()
G4Fragment *	GetNucleus ()
void	SetNucleus (G4Fragment *nucleus)
void	SetVerboseLevel (G4int verbose)
void	Initialize ()
void	SetEO (G4ElectronOccupancy eo)
void	SetVaccantSN (G4int val)
G4ElectronOccupancy	GetEO ()
G4int	GetVacantSN ()
void	SetTimeLimit (G4double value)

Additional Inherited Members
Protected Member Functions inherited from G4VGammaDeexcitation
void	Update ()
Protected Attributes inherited from G4VGammaDeexcitation
G4VGammaTransition *	_transition
G4int	_verbose

Detailed Description

Definition at line 68 of file G4DiscreteGammaDeexcitation.hh.

Constructor & Destructor Documentation

G4DiscreteGammaDeexcitation::G4DiscreteGammaDeexcitation

(

)

Definition at line 67 of file G4DiscreteGammaDeexcitation.cc.

References G4NuclearLevelStore::GetInstance().

                                                        : 
  _nucleusZ(0), _nucleusA(0), _max_hl(1e-6*second), _icm(false),
  _rdm(false), _levelManager(0)
{
  _tolerance = CLHEP::keV;
  G4NuclearLevelStore::GetInstance();
}

G4DiscreteGammaDeexcitation::~G4DiscreteGammaDeexcitation ( )

virtual

Definition at line 75 of file G4DiscreteGammaDeexcitation.cc.

{}

Member Function Documentation

G4bool G4DiscreteGammaDeexcitation::CanDoTransition ( )

virtual

Implements G4VGammaDeexcitation.

Definition at line 130 of file G4DiscreteGammaDeexcitation.cc.

References G4VGammaDeexcitation::_transition, G4VGammaDeexcitation::_verbose, G4cout, G4endl, G4Fragment::GetExcitationEnergy(), G4VGammaDeexcitation::GetNucleus(), G4NuclearLevel::HalfLife(), G4NuclearLevelManager::MaxLevelEnergy(), G4NuclearLevelManager::MinLevelEnergy(), and G4NuclearLevelManager::NearestLevel().

{

  G4bool canDo = true;
    
  if (_transition == 0) {
    canDo = false;
    
    if (_verbose > 0)
      G4cout 
    << "G4DiscreteGammaDeexcitation::CanDoTransition - Null transition " 
    << G4endl;
  } 
  if (canDo)  {
    //remove the   _nucleusZ>98 limit L. Desorgher Sept. 2011
    //if (_nucleusZ<2 || _nucleusA<3 || _nucleusZ>98)
     if (_nucleusZ<2 || _nucleusA<3 )
      {
    canDo = false;
    if (_verbose > 0) 
      G4cout 
        << "G4DiscreteGammaDeexcitation::CanDoTransition - n/p/H/>Cf"
        << G4endl;
      }
  }

  G4Fragment* nucleus = GetNucleus();
  G4double excitation = nucleus->GetExcitationEnergy();
  //G4cout << "G4DiscreteGammaDeexcitation::CanDoTransition: " << nucleus << G4endl;

  if (canDo) {
    if (excitation <= _tolerance) {
      canDo = false;
      if (_verbose > 0) {
    G4cout 
      << "G4DiscreteGammaDeexcitation::CanDoTransition -  Excitation <= 0" 
      << excitation << "  " << excitation - _tolerance
      << G4endl;
      }
    } else { 
      if (excitation > _levelManager->MaxLevelEnergy() + _tolerance) { canDo = false; }
      //if (excitation < _levelManager->MinLevelEnergy() - _tolerance) canDo = false;  
      // The following is a protection to avoid looping in case of elements with very low
      // ensdf levels
      //if (excitation < _levelManager->MinLevelEnergy() * 0.9) canDo = false;  
  
      if (_verbose > 0) {
    G4cout << "G4DiscreteGammaDeexcitation::CanDoTransition -  Excitation " 
           << excitation << ", Min-Max are " 
           << _levelManager->MinLevelEnergy() << " "
           << _levelManager->MaxLevelEnergy() << G4endl;
      }
    }
  }
  
  if (canDo) {
    const G4NuclearLevel* level = _levelManager->NearestLevel(excitation);  
    if (!level) { 
      canDo = false; 
 
    } else {
      if (level->HalfLife() > _max_hl && !_rdm ) { canDo = false; }
      
      if (_verbose > 0) {
    G4cout << "G4DiscreteGammaDeexcitation::CanDoTransition -  Halflife " 
           << level->HalfLife() << ", Calling from RDM " 
           << (_rdm ? " True " : " False ")  << ", Max-HL = " <<  _max_hl 
           << G4endl;
      }
    }
  }
  if (_verbose > 0) {
    G4cout <<"G4DiscreteGammaDeexcitation::CanDoTransition - CanDo: " 
       <<  (canDo ? " True " : " False ")  << G4endl; 
  }
  
  return canDo;
      
}

G4VGammaTransition * G4DiscreteGammaDeexcitation::CreateTransition ( )

virtual

Implements G4VGammaDeexcitation.

Definition at line 78 of file G4DiscreteGammaDeexcitation.cc.

References G4VGammaDeexcitation::_verbose, G4NuclearLevel::Energy(), G4cout, G4endl, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4NuclearLevelStore::GetInstance(), G4NuclearLevelStore::GetManager(), G4VGammaDeexcitation::GetNucleus(), G4Fragment::GetZ_asInt(), G4NuclearLevelManager::IsValid(), G4NuclearLevelManager::NearestLevel(), and G4DiscreteGammaTransition::SetICM().

{
  G4Fragment* nucleus = GetNucleus();
  G4int A = nucleus->GetA_asInt();
  G4int Z = nucleus->GetZ_asInt();
  //  _verbose =2;
  //  G4cout << "G4DiscreteGammaDeexcitation::CreateTransition: " << nucleus << G4endl;
  if (_nucleusA != A || _nucleusZ != Z) 
    {
      _nucleusA = A;
      _nucleusZ = Z;
      _levelManager = G4NuclearLevelStore::GetInstance()->GetManager(Z,A);
    }

  if (_levelManager->IsValid()) 
    { 
      if (_verbose > 1)
    {
      G4cout 
        << "G4DiscreteGammaDeexcitation::CreateTransition - (A,Z) is valid " 
        << G4endl;
    }
    
      G4double excitation = nucleus->GetExcitationEnergy();
      const G4NuclearLevel* level =_levelManager->NearestLevel(excitation);
    
      if (level != 0)  
    { 
      if (_verbose > 0) {
        G4cout 
          << "G4DiscreteGammaDeexcitation::CreateTransition - Created from level energy " 
          << level->Energy() << ", excitation is " 
          << excitation << G4endl;
      }
      G4DiscreteGammaTransition* dtransit = new G4DiscreteGammaTransition(*level,Z,A);
      dtransit->SetICM(_icm);  
      return dtransit;
    }
      else 
    { 
      if (_verbose > 0) {
        G4cout 
          << "G4DiscreteGammaDeexcitation::CreateTransition - No transition created from "
          << excitation << " within tolerance " << _tolerance << G4endl;
      }
      return 0; 
    }
    }
  return 0;
}

void G4DiscreteGammaDeexcitation::SetHL ( G4double  hl )

inline

Definition at line 89 of file G4DiscreteGammaDeexcitation.hh.

{ _max_hl = hl; };

void G4DiscreteGammaDeexcitation::SetICM ( G4bool  hl )

inline

Definition at line 85 of file G4DiscreteGammaDeexcitation.hh.

Referenced by G4PhotonEvaporation::G4PhotonEvaporation().

{ _icm = hl; };

void G4DiscreteGammaDeexcitation::SetRDM ( G4bool  hl )

inline

Definition at line 87 of file G4DiscreteGammaDeexcitation.hh.

{ _rdm = hl; };

---

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

• G4DiscreteGammaDeexcitation.hh
• G4DiscreteGammaDeexcitation.cc