#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 60 of file G4NuclearLevelData.hh.

Constructor & Destructor Documentation

◆ ~G4NuclearLevelData()

G4NuclearLevelData::~G4NuclearLevelData ( )

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

{
  G4bool res = false; 
  if(Z > 0 && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) { 
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&nuclearLevelDataMutex);
#endif
    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;
    }
#ifdef G4MULTITHREADED
    G4MUTEXUNLOCK(&nuclearLevelDataMutex);
#endif
  } 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 G4RadioactiveDecayBaseMessenger::SetNewValue(), and G4RadioactiveDecaymessenger::SetNewValue().

◆ FindLevel()

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

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

{
  if (!theInstance)  { 
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&nuclearLevelDataMutex);
    if (!theInstance)  { 
#endif
      static G4NuclearLevelData theData;
      theInstance = &theData; 
#ifdef G4MULTITHREADED
    }
    G4MUTEXUNLOCK(&nuclearLevelDataMutex);
#endif
  }
  return theInstance;
}   

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4RadioactiveDecay::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(), G4RadioactiveDecayBaseMessenger::SetNewValue(), and G4RadioactiveDecaymessenger::SetNewValue().

◆ GetLevelDensity()

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

Definition at line 634 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 566 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 467 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] ) {
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&nuclearLevelDataMutex);
    if( !(fLevelManagerFlags[Z])[idx] ) {
#endif
      (fLevelManagers[Z])[idx] = fLevelReader->CreateLevelManager(Z, A);
      (fLevelManagerFlags[Z])[idx] = true;
#ifdef G4MULTITHREADED
    }
    G4MUTEXUNLOCK(&nuclearLevelDataMutex);
#endif
  }
  return (fLevelManagers[Z])[idx];
}

Referenced by G4Evaporation::BreakFragment(), G4UnstableFragmentBreakUp::BreakUpChain(), G4Radioactivation::CalculateChainsFromParent(), G4RadioactiveDecay::CalculateChainsFromParent(), FindLevel(), G4GEMChannelVI::G4GEMChannelVI(), GetLevelDensity(), GetLevelEnergy(), GetLowEdgeLevelEnergy(), and StreamLevels().

◆ GetLowEdgeLevelEnergy()

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

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

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

◆ GetMaxLevelEnergy()

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

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

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

◆ GetPairingCorrection() [1/2]

G4PairingCorrection * G4NuclearLevelData::GetPairingCorrection ( )

Definition at line 624 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 644 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 619 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 629 of file G4NuclearLevelData.cc.

{
  return fShellCorrection;
}

Referenced by G4FissionBarrier::G4FissionBarrier(), and G4LevelManager::G4LevelManager().

◆ MaxLevelEnergy()

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

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