#include <G4NuclearLevelData.hh>

Public Member Functions
	~G4NuclearLevelData ()

const G4LevelManager *	GetLevelManager (G4int Z, G4int A)

G4bool	AddPrivateData (G4int Z, G4int A, const G4String &filename)

G4int	GetMinA (G4int Z) const

G4int	GetMaxA (G4int Z) const

G4double	GetMaxLevelEnergy (G4int Z, G4int A) const

G4float	MaxLevelEnergy (G4int Z, G4int A) const

G4double	GetLevelEnergy (G4int Z, G4int A, G4double energy)

G4double	GetLowEdgeLevelEnergy (G4int Z, G4int A, G4double energy)

G4double	FindLevel (G4int Z, G4int A, G4double resMass, G4double Mass, G4double partMass, G4double T)

G4DeexPrecoParameters *	GetParameters ()

G4PairingCorrection *	GetPairingCorrection ()

G4ShellCorrection *	GetShellCorrection ()

G4double	GetLevelDensity (G4int Z, G4int A, G4double U)

G4double	GetPairingCorrection (G4int Z, G4int A)

void	UploadNuclearLevelData (G4int Zlim)

void	StreamLevels (std::ostream &os, G4int Z, G4int A)

	G4NuclearLevelData (G4NuclearLevelData &)=delete

G4NuclearLevelData &	operator= (const G4NuclearLevelData &right)=delete

Static Public Member Functions
static G4NuclearLevelData *	GetInstance ()

Detailed Description

Definition at line 59 of file G4NuclearLevelData.hh.

Constructor & Destructor Documentation

◆ ~G4NuclearLevelData()

G4NuclearLevelData::~G4NuclearLevelData ( )

Definition at line 449 of file G4NuclearLevelData.cc.

{
  delete fLevelReader;
  delete fDeexPrecoParameters;
  delete fShellCorrection;
  delete fPairingCorrection;
  for(G4int Z=1; Z<ZMAX; ++Z) {
    size_t nn = (fLevelManagers[Z]).size();
    for(size_t j=0; j<nn; ++j) { 
      delete (fLevelManagers[Z])[j]; 
    }
  }
}

◆ G4NuclearLevelData()

G4NuclearLevelData::G4NuclearLevelData ( G4NuclearLevelData & )

delete

Member Function Documentation

◆ AddPrivateData()

G4bool G4NuclearLevelData::AddPrivateData	(	G4int	Z,
		G4int	A,
		const G4String &	filename )

Definition at line 480 of file G4NuclearLevelData.cc.

{
  G4bool res = false; 
  if(Z > 0 && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) { 
    G4AutoLock l(&nuclearLevelDataMutex);
    const G4LevelManager* newman = fLevelReader->MakeLevelManager(Z, A, fname);
    // if file is corrupted G4LevelReader should make the exception
    if(newman != nullptr) {
      res = true;
      if(0 < fDeexPrecoParameters->GetVerbose()) {
    G4cout << "G4NuclearLevelData::AddPrivateData for Z= " << Z
           << " A= " << A << " from <" << fname 
           << "> is done" << G4endl;
      }
      const G4int idx = A - AMIN[Z];
      delete (fLevelManagers[Z])[idx]; 
      (fLevelManagers[Z])[idx] = newman;
      (fLevelManagerFlags[Z])[idx] = true;
    }
    l.unlock();
  } else {
    G4ExceptionDescription ed;
    ed << "private nuclear level data for Z= " << Z << " A= " << A
       << " outside allowed limits ";
    G4Exception("G4NuclearLevelData::AddPrivateData","had0433",FatalException,
                ed,"Stop execution");
  }
  return res;
}

Referenced by G4RadioactiveDecayMessenger::SetNewValue().

◆ FindLevel()

G4double G4NuclearLevelData::FindLevel	(	G4int	Z,
		G4int	A,
		G4double	resMass,
		G4double	Mass,
		G4double	partMass,
		G4double	T )

Definition at line 579 of file G4NuclearLevelData.cc.

{
  G4double T = ekin;
  G4double E2 = (Mass - partMass)*(Mass - partMass);
  G4double Eex = std::sqrt(E2 - 2.0*Mass*T) - resMass;
  if(Eex <= GetMaxLevelEnergy(Z, A)) {
    if(Eex <= 0.0) { 
      Eex = 0.0;
      T = (E2 - resMass*resMass)*0.5/Mass; 
    } else {
      const G4LevelManager* man = GetLevelManager(Z, A);
      if(man != nullptr) {
        size_t idx = man->NearestLevelIndex(Eex, man->NumberOfTransitions());
        for(;;) {
          Eex = (G4double)man->LevelEnergy(idx);
          T = (E2 - (resMass+Eex)*(resMass+Eex))*0.5/Mass;
          if(T >= 0.0 || 0 == idx) { break; }
          --idx; 
        }
      }
    }
    T = std::max(T, 0.0);
  }
  return T;
}

◆ GetInstance()

G4NuclearLevelData * G4NuclearLevelData::GetInstance ( )

static

Definition at line 423 of file G4NuclearLevelData.cc.

{
  if (nullptr == theInstance)  { 
    G4AutoLock l(&nuclearLevelDataMutex);
    if (nullptr == theInstance)  { 
      static G4NuclearLevelData theData;
      theInstance = &theData; 
    }
    l.unlock();
  }
  return theInstance;
}   

Referenced by G4NeutronRadCaptureHP::BuildPhysicsTable(), G4Radioactivation::CalculateChainsFromParent(), G4GNASHTransitions::CalculateProbability(), LBE::ConstructGeneral(), G4CompetitiveFission::G4CompetitiveFission(), G4Evaporation::G4Evaporation(), G4EvaporationChannel::G4EvaporationChannel(), G4EvaporationLevelDensityParameter::G4EvaporationLevelDensityParameter(), G4FissionBarrier::G4FissionBarrier(), G4FissionLevelDensityParameter::G4FissionLevelDensityParameter(), G4FissionLevelDensityParameterINCLXX::G4FissionLevelDensityParameterINCLXX(), G4FissionProbability::G4FissionProbability(), G4GEMChannel::G4GEMChannel(), G4GEMChannelVI::G4GEMChannelVI(), G4GEMProbability::G4GEMProbability(), G4GNASHTransitions::G4GNASHTransitions(), G4IonINCLXXPhysics::G4IonINCLXXPhysics(), G4IonPhysicsPHP::G4IonPhysicsPHP(), G4IonQMDPhysics::G4IonQMDPhysics(), G4LevelManager::G4LevelManager(), G4LightIonQMDPhysics::G4LightIonQMDPhysics(), G4PhotonEvaporation::G4PhotonEvaporation(), G4PreCompoundEmission::G4PreCompoundEmission(), G4PreCompoundModel::G4PreCompoundModel(), G4PreCompoundTransitions::G4PreCompoundTransitions(), G4RadioactiveDecayPhysics::G4RadioactiveDecayPhysics(), G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp(), G4VEmissionProbability::G4VEmissionProbability(), G4VPreCompoundFragment::G4VPreCompoundFragment(), G4ExcitationHandler::Initialise(), G4FermiBreakUpVI::Initialise(), G4FermiFragmentsPoolVI::Initialise(), G4NeutronFissionVI::InitialiseModel(), G4NeutronRadCapture::InitialiseModel(), G4PhysListUtil::InitialiseParameters(), G4RadioactiveDecayMessenger::SetNewValue(), and G4RadioactiveDecay::StreamInfo().

◆ GetLevelDensity()

G4double G4NuclearLevelData::GetLevelDensity	(	G4int	Z,
		G4int	A,
		G4double	U )

Definition at line 621 of file G4NuclearLevelData.cc.

{
  if(fDeexPrecoParameters->GetLevelDensityFlag()) {
    return A*fDeexPrecoParameters->GetLevelDensity();
  }
  const G4LevelManager* man = GetLevelManager(Z, A);
  return (man) ? man->LevelDensity(U) 
    : 0.058025*A*(1.0 + 5.9059/fG4calc->Z13(A));
}

Referenced by G4GNASHTransitions::CalculateProbability(), G4PreCompoundTransitions::CalculateProbability(), G4EvaporationProbability::ComputeProbability(), G4GEMProbabilityVI::ComputeProbability(), G4GEMProbabilityVI::ComputeTotalProbability(), G4PreCompoundModel::DeExcite(), G4PhotonEvaporation::GetEmissionProbability(), G4EvaporationLevelDensityParameter::LevelDensityParameter(), G4FissionLevelDensityParameter::LevelDensityParameter(), G4FissionLevelDensityParameterINCLXX::LevelDensityParameter(), G4PreCompoundIon::ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), G4HETCNeutron::SampleKineticEnergy(), and G4EvaporationProbability::TotalProbability().

◆ GetLevelEnergy()

G4double G4NuclearLevelData::GetLevelEnergy	(	G4int	Z,
		G4int	A,
		G4double	energy )

Definition at line 553 of file G4NuclearLevelData.cc.

{
  G4double e = energy;
  if(e <= GetMaxLevelEnergy(Z, A)) {
    const G4LevelManager* man = GetLevelManager(Z, A);
    if(man != nullptr) { 
      e = (G4double)man->NearestLevelEnergy(e, man->NumberOfTransitions()); 
    }
  }
  return e;
}

◆ GetLevelManager()

const G4LevelManager * G4NuclearLevelData::GetLevelManager	(	G4int	Z,
		G4int	A )

Definition at line 464 of file G4NuclearLevelData.cc.

{
  if(Z < 1 || Z >= ZMAX || A < AMIN[Z] || A > AMAX[Z]) { return nullptr; } 
  const G4int idx = A - AMIN[Z];
  if( !(fLevelManagerFlags[Z])[idx] ) {
    G4AutoLock l(&nuclearLevelDataMutex);
    if( !(fLevelManagerFlags[Z])[idx] ) {
      (fLevelManagers[Z])[idx] = fLevelReader->CreateLevelManager(Z, A);
      (fLevelManagerFlags[Z])[idx] = true;
    }
    l.unlock();
  }
  return (fLevelManagers[Z])[idx];
}

Referenced by G4Radioactivation::CalculateChainsFromParent(), FindLevel(), G4GEMChannelVI::G4GEMChannelVI(), GetLevelDensity(), GetLevelEnergy(), GetLowEdgeLevelEnergy(), G4FermiFragmentsPoolVI::Initialise(), and StreamLevels().

◆ GetLowEdgeLevelEnergy()

G4double G4NuclearLevelData::GetLowEdgeLevelEnergy	(	G4int	Z,
		G4int	A,
		G4double	energy )

Definition at line 566 of file G4NuclearLevelData.cc.

{
  G4double e = GetMaxLevelEnergy(Z, A);
  if(energy < e) { 
    const G4LevelManager* man = GetLevelManager(Z, A);
    if(man != nullptr) { 
      e = (G4double)man->NearestLowEdgeLevelEnergy(energy); 
    }
  }
  return e;
}

◆ GetMaxA()

G4int G4NuclearLevelData::GetMaxA ( G4int Z ) const

Definition at line 515 of file G4NuclearLevelData.cc.

{
  return (Z >= 0 && Z < ZMAX) ? AMAX[Z] : 0; 
}

Referenced by G4FermiFragmentsPoolVI::Initialise().

◆ GetMaxLevelEnergy()

G4double G4NuclearLevelData::GetMaxLevelEnergy	(	G4int	Z,
		G4int	A ) const

Definition at line 541 of file G4NuclearLevelData.cc.

{
  return (0 < Z && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) ?
    (G4double)(LEVELMAX[LEVELIDX[Z] + A - AMIN[Z]]) : 0.0;
}

Referenced by FindLevel(), GetLevelEnergy(), and GetLowEdgeLevelEnergy().

◆ GetMinA()

G4int G4NuclearLevelData::GetMinA ( G4int Z ) const

Definition at line 510 of file G4NuclearLevelData.cc.

{
  return (Z >= 0 && Z < ZMAX) ? AMIN[Z] : 0;
}

Referenced by G4FermiFragmentsPoolVI::Initialise().

◆ GetPairingCorrection() [1/2]

G4PairingCorrection * G4NuclearLevelData::GetPairingCorrection ( )

Definition at line 611 of file G4NuclearLevelData.cc.

{
  return fPairingCorrection;
}

Referenced by G4GEMProbabilityVI::ComputeTotalProbability(), G4CompetitiveFission::G4CompetitiveFission(), G4FissionProbability::G4FissionProbability(), G4GEMChannelVI::G4GEMChannelVI(), G4GEMProbability::GetCoulombBarrier(), G4EvaporationChannel::GetEmissionProbability(), G4GEMChannel::GetEmissionProbability(), and G4EvaporationProbability::TotalProbability().

◆ GetPairingCorrection() [2/2]

G4double G4NuclearLevelData::GetPairingCorrection	(	G4int	Z,
		G4int	A )

Definition at line 631 of file G4NuclearLevelData.cc.

{
  if(fDeexPrecoParameters->GetLevelDensityFlag()) {
    return fPairingCorrection->GetPairingCorrection(A, Z);
  }
  G4int N = A - Z;
  const G4double par = 12.*CLHEP::MeV;
  G4double x = (A <= 36) ? 6.0 : std::sqrt((G4double)A);
  return (2 - Z + (Z/2)*2 - N + (N/2)*2)*par/x;
}

◆ GetParameters()

G4DeexPrecoParameters * G4NuclearLevelData::GetParameters ( )

Definition at line 606 of file G4NuclearLevelData.cc.

{
  return fDeexPrecoParameters;
}

Referenced by G4NeutronRadCaptureHP::BuildPhysicsTable(), LBE::ConstructGeneral(), G4GNASHTransitions::G4GNASHTransitions(), G4IonINCLXXPhysics::G4IonINCLXXPhysics(), G4IonPhysicsPHP::G4IonPhysicsPHP(), G4IonQMDPhysics::G4IonQMDPhysics(), G4LightIonQMDPhysics::G4LightIonQMDPhysics(), G4PreCompoundEmission::G4PreCompoundEmission(), G4PreCompoundTransitions::G4PreCompoundTransitions(), G4RadioactiveDecayPhysics::G4RadioactiveDecayPhysics(), G4VPreCompoundFragment::G4VPreCompoundFragment(), G4ExcitationHandler::Initialise(), G4FermiBreakUpVI::Initialise(), G4FermiFragmentsPoolVI::Initialise(), G4PhotonEvaporation::Initialise(), G4VEmissionProbability::Initialise(), G4Evaporation::InitialiseChannels(), G4NeutronFissionVI::InitialiseModel(), G4NeutronRadCapture::InitialiseModel(), G4PreCompoundModel::InitialiseModel(), and G4RadioactiveDecay::StreamInfo().

◆ GetShellCorrection()

G4ShellCorrection * G4NuclearLevelData::GetShellCorrection ( )

Definition at line 616 of file G4NuclearLevelData.cc.

{
  return fShellCorrection;
}

Referenced by G4FissionBarrier::G4FissionBarrier().

◆ MaxLevelEnergy()

G4float G4NuclearLevelData::MaxLevelEnergy	(	G4int	Z,
		G4int	A ) const

Definition at line 547 of file G4NuclearLevelData.cc.

{
  return (0 < Z && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) ?
    (LEVELMAX[LEVELIDX[Z] + A - AMIN[Z]]) : 0.0f;
}

◆ operator=()

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