#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 Z)

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 450 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 481 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, filename);
    if(newman) { 
      res = true;
      G4cout << "G4NuclearLevelData::AddPrivateData for Z= " << Z
             << " A= " << A << " from <" << filename 
             << "> 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 578 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 G4Radioactivation::CalculateChainsFromParent(), G4GNASHTransitions::CalculateProbability(), LBE::ConstructGeneral(), G4CompetitiveFission::G4CompetitiveFission(), G4Evaporation::G4Evaporation(), G4EvaporationChannel::G4EvaporationChannel(), G4EvaporationLevelDensityParameter::G4EvaporationLevelDensityParameter(), G4FermiFragmentsPoolVI::G4FermiFragmentsPoolVI(), G4FissionBarrier::G4FissionBarrier(), G4FissionLevelDensityParameter::G4FissionLevelDensityParameter(), G4FissionLevelDensityParameterINCLXX::G4FissionLevelDensityParameterINCLXX(), G4FissionProbability::G4FissionProbability(), G4GEMChannel::G4GEMChannel(), G4GEMChannelVI::G4GEMChannelVI(), G4GEMProbability::G4GEMProbability(), G4GNASHTransitions::G4GNASHTransitions(), G4IonINCLXXPhysics::G4IonINCLXXPhysics(), G4IonPhysicsPHP::G4IonPhysicsPHP(), G4IonQMDPhysics::G4IonQMDPhysics(), G4LevelManager::G4LevelManager(), G4PhotonEvaporation::G4PhotonEvaporation(), G4PreCompoundEmission::G4PreCompoundEmission(), G4PreCompoundModel::G4PreCompoundModel(), G4PreCompoundTransitions::G4PreCompoundTransitions(), G4RadioactiveDecayPhysics::G4RadioactiveDecayPhysics(), G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp(), G4VEmissionProbability::G4VEmissionProbability(), G4VPreCompoundFragment::G4VPreCompoundFragment(), G4ExcitationHandler::Initialise(), G4NeutronRadCapture::InitialiseModel(), G4PhysListUtil::InitialiseParameters(), and G4RadioactiveDecayMessenger::SetNewValue().

◆ GetLevelDensity()

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

Definition at line 620 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(), G4GEMProbabilityVI::ComputeProbability(), G4EvaporationProbability::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 552 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 465 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(), and StreamLevels().

◆ GetLowEdgeLevelEnergy()

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

Definition at line 565 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 514 of file G4NuclearLevelData.cc.

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

◆ GetMaxLevelEnergy()

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

Definition at line 540 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 509 of file G4NuclearLevelData.cc.

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

◆ GetPairingCorrection() [1/2]

G4PairingCorrection * G4NuclearLevelData::GetPairingCorrection ( )

Definition at line 610 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 630 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 605 of file G4NuclearLevelData.cc.

{
  return fDeexPrecoParameters;
}

Referenced by LBE::ConstructGeneral(), G4FermiFragmentsPoolVI::G4FermiFragmentsPoolVI(), G4GNASHTransitions::G4GNASHTransitions(), G4IonINCLXXPhysics::G4IonINCLXXPhysics(), G4IonPhysicsPHP::G4IonPhysicsPHP(), G4IonQMDPhysics::G4IonQMDPhysics(), G4LevelReader::G4LevelReader(), G4PreCompoundEmission::G4PreCompoundEmission(), G4PreCompoundTransitions::G4PreCompoundTransitions(), G4RadioactiveDecayPhysics::G4RadioactiveDecayPhysics(), G4VPreCompoundFragment::G4VPreCompoundFragment(), G4ExcitationHandler::Initialise(), G4VEmissionProbability::Initialise(), G4PhotonEvaporation::Initialise(), G4Evaporation::InitialiseChannels(), G4NeutronRadCapture::InitialiseModel(), and G4PreCompoundModel::InitialiseModel().

◆ GetShellCorrection()

G4ShellCorrection * G4NuclearLevelData::GetShellCorrection ( )

Definition at line 615 of file G4NuclearLevelData.cc.

{
  return fShellCorrection;
}

Referenced by G4FissionBarrier::G4FissionBarrier().

◆ MaxLevelEnergy()

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

Definition at line 546 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 )