#include <G4LENDManager.hh>

Public Member Functions
G4GIDI_target *	GetLENDTarget (G4ParticleDefinition *, G4String, G4int iZ, G4int iA, G4int iM=0)

std::vector< G4String >	IsLENDTargetAvailable (G4ParticleDefinition *, G4int iZ, G4int iA, G4int iM=0)

G4int	GetNucleusEncoding (G4int iZ, G4int iA, G4int iM)

G4NistElementBuilder *	GetNistElementBuilder ()

G4int	GetVerboseLevel ()

G4bool	RequestChangeOfVerboseLevel (G4int)

G4double	GetExcitationEnergyOfExcitedIsomer (G4int, G4int, G4int)

Static Public Member Functions
static G4LENDManager *	GetInstance ()

Detailed Description

Definition at line 63 of file G4LENDManager.hh.

Member Function Documentation

◆ GetExcitationEnergyOfExcitedIsomer()

G4double G4LENDManager::GetExcitationEnergyOfExcitedIsomer	(	G4int	iZ,
		G4int	iA,
		G4int	iM )

Definition at line 407 of file G4LENDManager.cc.

{
   G4double EE = 0.0;
   G4int nucCode = GetNucleusEncoding( iZ , iA , iM );
   auto it = mExcitationEnergy.find( nucCode );
   if ( it != mExcitationEnergy.end() ) {
      EE = it->second;
   } else {
      if ( iM == 0 ) {
         G4cout << "G4LENDManager::GetExcitationEnergyOfExcitedIsomer is called for ground state (iM=0) nucleus" << G4endl;
      } else {
         G4cout << "Can not find excitation energy for Z = " << iZ << ", A = " << iA << ", M = " << iM << " and the energy set to 0." << G4endl;
      }
   }
   return EE;
}

Referenced by G4LENDElastic::ApplyYourself().

◆ GetInstance()

static G4LENDManager * G4LENDManager::GetInstance ( )

inlinestatic

Definition at line 77 of file G4LENDManager.hh.

      {
         if ( lend_manager == NULL) lend_manager = new G4LENDManager();
         return lend_manager;
      };

Referenced by G4LENDElastic::ApplyYourself(), G4LENDCrossSection::G4LENDCrossSection(), and G4LENDModel::G4LENDModel().

◆ GetLENDTarget()

G4GIDI_target * G4LENDManager::GetLENDTarget	(	G4ParticleDefinition *	proj,
		G4String	evaluation,
		G4int	iZ,
		G4int	iA,
		G4int	iM = 0 )

Definition at line 208 of file G4LENDManager.cc.

{
 
   G4GIDI_target* anLENDTarget = NULL;
 
   if ( iM > 9 ) { 
      throw G4HadronicException(__FILE__, __LINE__, "Requested isomer level of target is too high." ); 
   }
 
   G4int iTarg = GetNucleusEncoding( iZ , iA , iM );
 
   // Searching in current map
   for ( std::vector < lend_target >::iterator 
         it = v_lend_target.begin() ; it != v_lend_target.end() ; it++ )
   {
      if ( it->proj == proj && it->target_code == iTarg && it->evaluation == evaluation ) 
      {
         //find! 
         return it->target;
      }
   }
 
 
   if ( proj_lend_map.find ( proj ) == proj_lend_map.end() ) {
      G4cout << proj->GetParticleName() << " is not supported by this LEND library." << G4endl;
      return anLENDTarget; // return NULL 
   }
 
   G4GIDI* xlend = proj_lend_map.find ( proj ) -> second; 
 
   if ( xlend->isThisDataAvailable( evaluation, iZ, iA , iM ) )
   {
 
      if ( verboseLevel > 1 ) {
         G4cout << evaluation << " for " << ionTable->GetIonName( iZ , iA , 0 ) 
                 << " with Isomer level of " << iM  << " is exist in this LEND." << G4endl;
      }
 
      anLENDTarget = xlend->readTarget( evaluation , iZ , iA , iM );
 
      lend_target new_target; 
      new_target.lend = xlend; 
      new_target.target = anLENDTarget; 
      new_target.proj = proj;
      new_target.evaluation = evaluation;
      new_target.target_code = iTarg;
      
      v_lend_target.push_back( new_target );
 
//    found EXACT
      return anLENDTarget;
 
   }
   else 
   {
//    NO EXACT DATA (Evaluatino & Z,A,M)
                                                                        
      //Searching available evaluation and natural abundance data and give suggestions.  
      //
      if ( verboseLevel > 1 ) 
         G4cout << evaluation << " for " << ionTable->GetIonName( iZ , iA , 0 ) 
                << " with Isomer level of " << iM << " is not exist in this LEND." << G4endl;
 
      std::vector< std::string >* available =  xlend->getNamesOfAvailableLibraries( iZ, iA , iM );
      if ( available->size() > 0 ) {
//       EXACT Z,A,M but Evaluation is different 
         if ( verboseLevel > 1 ) 
         {
            G4cout << " However you can use following evaluation(s) for the target. " << G4endl;
 
            std::vector< std::string >::iterator its;
            for ( its = available->begin() ; its != available->end() ; its++ ) 
               G4cout << *its << G4endl;
 
            G4cout << G4endl;
         }
      } else if ( xlend->isThisDataAvailable( evaluation, iZ, 0 , iM ) ) {
//      
//    checking natural abundance data for Z
//
//       EXACT natural abundance data for the evaluation 
         if ( verboseLevel > 1 ) 
            G4cout << " However you can use natural abundance data for the target. " << G4endl;
      }
      else
      {
         std::vector< std::string >* available_nat =  xlend->getNamesOfAvailableLibraries( iZ, 0 , iM );
//
         if ( available_nat->size() > 0 ) {
//          EXACT natural abundance data for Z but differnet evaluation
            if ( verboseLevel > 1 ) {
               G4cout << " However you can use following evaluation(s) for natural abundace of the target. " << G4endl;
 
               std::vector< std::string >::iterator its;
               for ( its = available_nat->begin() ; its != available_nat->end() ; its++ ) 
                  G4cout << *its << G4endl;
               G4cout << G4endl;
            }
         }
         delete available_nat;
      }
      delete available;
//    return NULL if exact data is not available               
      return anLENDTarget; // return NULL   
   }
   
   return anLENDTarget; 
}

◆ GetNistElementBuilder()

G4NistElementBuilder * G4LENDManager::GetNistElementBuilder ( )

inline

Definition at line 87 of file G4LENDManager.hh.

87{ return nistElementBuilder; };

Referenced by G4LENDCrossSection::create_used_target_map(), and G4LENDModel::create_used_target_map().

◆ GetNucleusEncoding()

G4int G4LENDManager::GetNucleusEncoding	(	G4int	iZ,
		G4int	iA,
		G4int	iM )

Definition at line 343 of file G4LENDManager.cc.

{   
   G4int value = ionTable->GetNucleusEncoding( iZ , iA ); // Ground State
                                                     // G4double E=0.0, G4int J=0);
   value += iM;
   return value;
}

Referenced by G4LENDCapture::ApplyYourself(), G4LENDElastic::ApplyYourself(), G4LENDFission::ApplyYourself(), G4LENDInelastic::ApplyYourself(), G4LENDModel::ApplyYourself(), G4LENDCrossSection::create_used_target_map(), G4LENDModel::create_used_target_map(), GetExcitationEnergyOfExcitedIsomer(), G4LENDCrossSection::GetIsoCrossSection(), GetLENDTarget(), G4LENDCombinedModel::HasData(), G4LENDGammaModel::HasData(), and G4LENDCrossSection::IsIsoApplicable().

◆ GetVerboseLevel()

G4int G4LENDManager::GetVerboseLevel ( )

inline

Definition at line 89 of file G4LENDManager.hh.

89{ return verboseLevel; };

Referenced by G4LENDCrossSection::DumpLENDTargetInfo(), G4LENDModel::DumpLENDTargetInfo(), and G4LENDModel::returnUnchanged().

◆ IsLENDTargetAvailable()

std::vector< G4String > G4LENDManager::IsLENDTargetAvailable	(	G4ParticleDefinition *	proj,
		G4int	iZ,
		G4int	iA,
		G4int	iM = 0 )

Definition at line 318 of file G4LENDManager.cc.

{
 
   std::vector< G4String > vEvaluation; 
   if ( proj_lend_map.find ( proj ) == proj_lend_map.end() ) 
   {
      G4cout << proj->GetParticleName() << " is not supported by this LEND." << G4endl;
      return vEvaluation; // return empty 
   }
 
   G4GIDI* xlend = proj_lend_map.find ( proj ) -> second; 
   std::vector< std::string >* available =  xlend->getNamesOfAvailableLibraries( iZ, iA , iM );
 
   if ( available->size() > 0 ) {
      std::vector< std::string >::iterator its;
      for ( its = available->begin() ; its != available->end() ; its++ ) 
         vEvaluation.push_back ( *its );
   }
   delete available;
 
   return vEvaluation;
}

◆ RequestChangeOfVerboseLevel()

G4bool G4LENDManager::RequestChangeOfVerboseLevel ( G4int newValue )

Definition at line 391 of file G4LENDManager.cc.

{
   G4bool result=false;
   if ( newValue >= verboseLevel) 
   {
      verboseLevel = newValue;
      result=true;
   }
   else
   {
      G4cout << "Since other LEND model or cross section have set the higher verbose level (" << verboseLevel << ") in LENDManager, you cannot change the value now." << G4endl; 
   }
 
   return result;
}