G4NeutronHPNames Class Reference

#include <G4NeutronHPNames.hh>

Public Member Functions
	G4NeutronHPNames ()
	G4NeutronHPNames (G4int maxOffSet)
	~G4NeutronHPNames ()
G4NeutronHPDataUsed	GetName (G4int A, G4int Z, G4String base, G4String rest, G4bool &active)
G4NeutronHPDataUsed	GetName (G4int A, G4int Z, G4int M, G4String base, G4String rest, G4bool &active)
G4String	GetName (G4int i)
void	SetMaxOffSet (G4int anOffset)
G4String	itoa (int current)

Data Fields
G4int	theMaxOffSet

Static Public Attributes
static const G4String	theString [100]

Detailed Description

Definition at line 38 of file G4NeutronHPNames.hh.

Constructor & Destructor Documentation

G4NeutronHPNames::G4NeutronHPNames ( )

inline

Definition at line 42 of file G4NeutronHPNames.hh.

References theMaxOffSet.

{theMaxOffSet = 5;}

G4NeutronHPNames::G4NeutronHPNames ( G4int  maxOffSet )

inline

Definition at line 43 of file G4NeutronHPNames.hh.

References theMaxOffSet.

{theMaxOffSet = maxOffSet;}

G4NeutronHPNames::~G4NeutronHPNames ( )

inline

Definition at line 44 of file G4NeutronHPNames.hh.

{}

Member Function Documentation

G4NeutronHPDataUsed G4NeutronHPNames::GetName ( G4int  A, G4int  Z, G4String  base, G4String  rest, G4bool &  active  )

inline

Definition at line 47 of file G4NeutronHPNames.hh.

References GetName().

Referenced by G4NeutronHPJENDLHEData::BuildPhysicsTable(), GetName(), G4NeutronHPIsoData::GetName(), G4NeutronHPFFFissionFS::Init(), G4NeutronHPFSFissionFS::Init(), G4NeutronHPFissionBaseFS::Init(), G4NeutronHPCaptureFS::Init(), G4NeutronHPElasticFS::Init(), G4NeutronHPIsoData::Init(), G4NeutronHPInelasticBaseFS::Init(), G4NeutronHPInelasticCompFS::Init(), G4FissionLibrary::Init(), and G4WendtFissionFragmentGenerator::InitializeANucleus().

{ G4int M = 0; return GetName( A, Z, M, base, rest, active); };

G4NeutronHPDataUsed G4NeutronHPNames::GetName	(	G4int	A,
		G4int	Z,
		G4int	M,
		G4String	base,
		G4String	rest,
		G4bool &	active
	)

Definition at line 66 of file G4NeutronHPNames.cc.

References G4cerr, G4cout, G4endl, G4NeutronHPDataUsed::GetA(), G4NeutronHPManager::GetDataStream(), G4NeutronHPManager::GetInstance(), G4NeutronHPDataUsed::GetName(), G4NeutronHPManager::GetVerboseLevel(), G4NeutronHPDataUsed::GetZ(), G4SandiaTable::GetZtoA(), G4NeutronHPDataUsed::IsThisNaturalAbundance(), itoa(), G4NeutronHPDataUsed::SetA(), G4NeutronHPDataUsed::SetM(), G4NeutronHPDataUsed::SetName(), G4NeutronHPDataUsed::SetNaturalAbundanceFlag(), G4NeutronHPDataUsed::SetZ(), theMaxOffSet, and theString.

{

   G4int verboseLevel = G4NeutronHPManager::GetInstance()->GetVerboseLevel();

   //G4cout << Z << " " << A << " " << M << " " << base << " " << rest << G4endl;

   //Excited isomer indicator
   std::stringstream ss;
   G4String sM;
   if ( M > 0 ) 
   {
      ss << "m";
      ss << M;
      ss >> sM;
      ss.clear();
   }

   G4NeutronHPDataUsed result;
   aFlag = true;
if(getenv("NeutronHPNames")) G4cout << "Names::GetName entered for Z = " << Z << ", A = " << A <<G4endl;

    G4int myA = A;
    G4int myZ = Z;

    if(Z>92.5&&!getenv("AllowForHeavyElements") ) 
    {
      //G4cerr << "Please contact Hans-Peter.Wellisch@cern.ch"<<G4endl;
      G4cerr << "Please contact Geant4 Hadron Group Coordinator"<<G4endl;
      throw G4HadronicException(__FILE__, __LINE__, "G4NeutronHPNames::GetName - data with Z>92 are not provided");
    }

    G4String * theName = 0;
    G4String theFileName("");

//    G4int inc = 1;

    G4int flip_Z = 1;
    G4int delta_Z = 0;

    G4int flip_A = 1;
    G4int delta_A = 0;
    
    //std::ifstream * check = new std::ifstream(".dummy");
    std::istringstream* check = NULL;
    G4bool first = true;
if(getenv("NeutronHPNames"))  G4cout << "entered GetName!!!"<<G4endl;
    do   
    {
       aFlag = true;
       G4String * biff = new G4String(); // delete here as theName
       *biff = base+"/CrossSection/"+itoa(myZ)+"_"+itoa(myA)+sM+"_"+theString[myZ-1];
      
       if(theName!=0) delete theName;
       theName = biff;
       result.SetName(*theName);
       result.SetA(myA);
       result.SetZ(myZ);
       result.SetM(M);
if(getenv("NeutronHPNames")) G4cout <<"HPWD 1 "<<*theName<<G4endl;

     // T.K. debug for memory leak
     if ( check != NULL ) {
        //check->close();
        delete check;
     } 

       //check = new std::ifstream(*theName);
       check = new std::istringstream(std::ios::in);
       G4NeutronHPManager::GetInstance()->GetDataStream(*theName,*check);
       if ( !(*check) ) 
       {
      //check->close();
      delete check;
          check = 0;
          aFlag = false;
          if ( first )
          {
             aFlag = true;
             first = false;
             biff = new G4String(); // delete here as theName
             *biff = base+"/CrossSection/"+itoa(myZ)+"_"+"nat"+"_"+theString[myZ-1];
             if(theName!=0) delete theName;
             theName = biff;
if(getenv("NeutronHPNames"))    G4cout <<"HPWD 2 "<<*theName<<G4endl;
             result.SetName(*theName);
             G4double natA = myZ/G4SandiaTable::GetZtoA(myZ);
             result.SetA(natA);
             result.SetZ(myZ);
             result.SetM(M);
             //check = new std::ifstream(*theName);
             check = new std::istringstream(std::ios::in);
             G4NeutronHPManager::GetInstance()->GetDataStream(*theName,*check);
             if ( !(*check) ) 
             {
                //check->close();
            delete check;
                check = 0;
                aFlag = false;
             }
             else
             {
                biff = new G4String(); // delete here as theName
                if(theName!=0) delete theName;
                *biff = base+"/"+rest+"/"+itoa(myZ)+"_"+"nat"+"_"+theString[myZ-1];  
                if ( rest=="/CrossSection" ) *biff = base+rest+"/"+itoa(myZ)+"_"+"nat"+"_"+theString[myZ-1];  
                theName = biff;
if(getenv("NeutronHPNames"))    G4cout <<"HPWD 3 "<<*theName<<G4endl;
                result.SetName(*theName);
                natA = myZ/G4SandiaTable::GetZtoA(myZ);
                result.SetA(natA);
                result.SetZ(myZ);
                result.SetM(M);
                result.SetNaturalAbundanceFlag();
             }
          }
       }
       else
       {
// 070706 T. Koi Modified 
/*
          biff = new G4String(); // delete here as theName
          *biff = base+"/"+rest+itoa(myZ)+"_"+itoa(myA)+"_"+theString[myZ-1];  
          if(theName!=0) delete theName;
          theName = biff;
if(getenv("NeutronHPNames"))    G4cout <<"HPWD 4 "<<*theName<<G4endl;
          result.SetName(*theName);
          result.SetA(myA);
          result.SetZ(myZ);
*/

          G4double tmpA = myA;
          //std::ifstream* file = NULL;
          std::istringstream* file = NULL;
          G4String fileName;

          if ( rest == "/CrossSection" )
          {

             //fileName = base+"/"+rest+"/"+itoa(myZ)+"_"+itoa(myA)+sM+"_"+theString[myZ-1];
             fileName = base+rest+"/"+itoa(myZ)+"_"+itoa(myA)+sM+"_"+theString[myZ-1];
if(getenv("NeutronHPNames"))    G4cout <<"HPWD 4a "<<*theName<<G4endl;

          }
          else
          {

// For FS
             fileName = base+"/"+rest+"/"+itoa(myZ)+"_"+itoa(myA)+sM+"_"+theString[myZ-1];
             //file = new std::ifstream(fileName);
             file = new std::istringstream(std::ios::in);
             G4NeutronHPManager::GetInstance()->GetDataStream(fileName,*file);

             if ( *file )
             {

// isotope FS
if(getenv("NeutronHPNames"))    G4cout <<"HPWD 4b1 "<<*theName<<G4endl;
             }
             else
             {

// _nat_ FS
                fileName  = base+"/"+rest+"/"+itoa(myZ)+"_"+"nat"+"_"+theString[myZ-1];

                delete file;
                //file = new std::ifstream(fileName);
                file = new std::istringstream(std::ios::in);
                G4NeutronHPManager::GetInstance()->GetDataStream(fileName,*file);
                if ( *file )
                {

// FS neither isotope nor _nat_
if(getenv("NeutronHPNames"))    G4cout <<"HPWD 4b2a "<<*theName<<G4endl;
                   G4double natA = myZ/G4SandiaTable::GetZtoA(myZ);
                   tmpA = natA;
                }
                else
                {
if(getenv("NeutronHPNames"))    G4cout <<"HPWD 4b2c "<<*theName<<G4endl;
                   fileName="INVALID";
                }
             }

             delete file;

          }

          result.SetName(fileName);
          result.SetA(tmpA);
          result.SetZ(myZ);
          result.SetM(M);

       }

       do 
       {
//        if (std::abs(myZ-Z)>theMaxOffSet||myZ==0||myA==0)
          if ( delta_Z > theMaxOffSet )
          { 
             //if ( inc > 0 )
             //{
             //   inc*= -1;
             //   myZ = Z;
             //   myA = A;
             //}
             //else
             //{
                G4cout <<"G4NeutronHPNames: Sorry, this material does not come near to any data."<<G4endl;
                G4cout <<"G4NeutronHPNames: Please make sure G4NEUTRONHPDATA points to the" << G4endl;
                G4cout <<"                  directory, the neutron scattering data are located in." << G4endl;
                G4cout << "G4NeutronHPNames: The material was: A="<<A<<", Z="<<Z<<G4endl;
                //throw G4HadronicException(__FILE__, __LINE__, "In case the data sets are at present not available in the neutron data library, please contact Hans-Peter.Wellisch@cern.ch");
                throw G4HadronicException(__FILE__, __LINE__, "In case the data sets are at present not available in the neutron data library, please contact Hadron Group Coordinator");
                delete theName;
                theFileName = "";
                return result;
             //}
          }

          //if ( std::abs( myA - A ) > theMaxOffSet )
          if ( delta_A > 2*theMaxOffSet )
          {
             delta_A = 0;
             flip_A = 1;

             first = true;

             if ( flip_Z > 0 ) 
             {
                delta_Z +=1; 
             }
             myZ = Z + flip_Z * delta_Z;
             flip_Z *= -1;
             
             myA = A;
             if ( myZ > 100 ) 
             {
                myZ = 100;
             }
             if ( myZ < 1 ) 
             {
                myZ = 1;
             }
              
//             myZ += inc;
          }
          else
          {
             if ( flip_A > 0 )
             {
                delta_A += 1;
             }
             myA = A + flip_A * delta_A; 
             flip_A *= -1;

             if ( myA < 1 ) 
             {
                myA = 1;
             }
              
//             myA += inc;
          }

       }
       while( myZ == 0 || myA == 0 );  // No meaning 

    }
    while((!check) || (!(*check)));

    if(getenv("NeutronHPNamesLogging") || getenv("NeutronHPNames")) 
    {
      G4cout << "Names::GetName: last theName proposal = "<< G4endl;
      G4cout << *theName <<" "<<A<<" "<<Z<<" "<<result.GetName()<<G4endl;
    }

// administration and anouncement for lacking of exact data in NDL 
    if ( Z != result.GetZ() || A != result.GetA() )
    {
       if ( rest == "/CrossSection" )
       {
          G4String reac = base;
          G4String dir = getenv("G4NEUTRONHPDATA"); 
          reac.erase ( 0 , dir.length() );
          if ( getenv ( "G4NEUTRONHP_SKIP_MISSING_ISOTOPES" ) && !( Z == result.GetZ() && result.IsThisNaturalAbundance() ) )
          {
             if ( verboseLevel > 0 ) {
                G4cout << "NeutronHP: " << reac << " file for Z = " << Z << ", A = " << A << " is not found and CrossSection set to 0." << G4endl;
             }
             G4String new_name = base+"/"+rest+"/"+"0_0_Zero";  
             result.SetName( new_name );
          }
          else
          { 
             //080901 Add protection that deuteron data do not selected for hydrogen and so on by T. Koi
             if ( ( reac.find("Inelastic") != reac.size() && ( (Z == 1 && A == 1) || (Z == 2 && A == 4) ) ) 
               || ( reac.find("Capture")   != reac.size() && ( (Z == 1 && A == 3) || (Z == 2 && A == 4) ) )  
               || ( reac.find("Fission")   != reac.size() && ( (Z == 88 && A == 224) || (Z == 88 && A == 225) || (Z == 89 && A == 225) || (Z == 88 && A == 226) ) ) ) 
                   
             {
                G4String new_name = base+"/"+rest+"/"+"0_0_Zero";
                result.SetName( new_name );
             }
             else
             {
                if ( verboseLevel > 0 ) {
                   G4cout << "NeutronHP: " << reac << " file for Z = " << Z << ", A = " << A << " is not found and NeutronHP will use " << result.GetName() << G4endl;
                }
             }
          }
       }
    }

    delete theName;
    if(aFlag)
    {
      //check->close();
      delete check;
      check = NULL;
    }
    return result;
  }

G4String G4NeutronHPNames::GetName

(

G4int

i

)

Definition at line 63 of file G4NeutronHPNames.cc.

References theString.

{ return theString[i]; }

G4String G4NeutronHPNames::itoa ( int  current )

inline

Definition at line 56 of file G4NeutronHPNames.hh.

Referenced by GetName().

  {
    const char theDigits[11] = "0123456789";
    G4String result;
    int digit;
    do
    {
      digit = current-10*(current/10);
      result=theDigits[digit]+result;
      current/=10;
    }
    while(current!=0);
    return result;
  }

void G4NeutronHPNames::SetMaxOffSet ( G4int  anOffset )

inline

Definition at line 50 of file G4NeutronHPNames.hh.

References theMaxOffSet.

{ theMaxOffSet = anOffset; }

Field Documentation

G4int G4NeutronHPNames::theMaxOffSet

Definition at line 55 of file G4NeutronHPNames.hh.

Referenced by G4NeutronHPNames(), GetName(), and SetMaxOffSet().

const G4String G4NeutronHPNames::theString

static

Initial value:

= {"Hydrogen", "Helium",
 "Lithium", "Berylium", "Boron", "Carbon", "Nitrogen", "Oxygen", "Fluorine",
 "Neon", "Sodium", "Magnesium", "Aluminum", "Silicon", "Phosphorous", 
 "Sulfur", "Chlorine", "Argon", "Potassium", "Calcium", "Scandium",
 "Titanium", "Vanadium", "Chromium", "Manganese", "Iron", "Cobalt", "Nickel",
 "Copper", "Zinc", "Gallium", "Germanium", "Arsenic", "Selenium", "Bromine",
 "Krypton", "Rubidium", "Strontium", "Yttrium", "Zirconium", "Niobium",
 "Molybdenum", "Technetium", "Ruthenium", "Rhodium", "Palladium", "Silver",
 "Cadmium", "Indium", "Tin", "Antimony", "Tellurium", "Iodine", "Xenon",
 "Cesium", "Barium", "Lanthanum", "Cerium", "Praseodymium", "Neodymium",
 "Promethium", "Samarium", "Europium", "Gadolinium", "Terbium", "Dysprosium",
 "Holmium", "Erbium", "Thulium", "Ytterbium", "Lutetium", "Hafnium",
 "Tantalum", "Tungsten", "Rhenium", "Osmium", "Iridium", "Platinium", "Gold",
 "Mercury", "Thallium", "Lead", "Bismuth", "Polonium", "Astatine", "Radon", 
 "Francium", "Radium", "Actinium", "Thorium", "Protactinium", "Uranium", 
 "Neptunium", "Plutonium", "Americium", "Curium", "Berkelium", "Californium",
 "Einsteinium","Fermium"}

Definition at line 54 of file G4NeutronHPNames.hh.

Referenced by GetName(), and G4FissionProductYieldDist::MakeIsotopeName().

---

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

• G4NeutronHPNames.hh
• G4NeutronHPNames.cc