G4INCL::HFB Namespace Reference

Functions
G4double	getRadiusParameterHFB (const ParticleType t, const G4int A, const G4int Z)
	Get the radius and diffuseness parameters from HFB calculations. More...
G4double	getSurfaceDiffusenessHFB (const ParticleType t, const G4int A, const G4int Z)
void	initialize ()

Function Documentation

getRadiusParameterHFB()

G4double G4INCL::HFB::getRadiusParameterHFB	(	const ParticleType	t,
		const G4int	A,
		const G4int	Z
	)

Get the radius and diffuseness parameters from HFB calculations.

Definition at line 130 of file G4INCLHFB.cc.

                                                                                      {
        // HFB calculations
      G4double r0=0.;
        if(t==Neutron)
         if(radiusN[Z][A]>0.)r0=radiusN[Z][A];
        if(t==Proton)
         if(radiusP[Z][A]>0.)r0=radiusP[Z][A];
      return r0;
    }

References A, G4INCL::Neutron, G4INCL::Proton, G4INCL::anonymous_namespace{G4INCLHFB.cc}::radiusN, G4INCL::anonymous_namespace{G4INCLHFB.cc}::radiusP, and Z.

Referenced by G4INCL::ParticleTable::getRadiusParameter(), and G4INCL::ParticleTable::getSurfaceDiffuseness().

getSurfaceDiffusenessHFB()

G4double G4INCL::HFB::getSurfaceDiffusenessHFB	(	const ParticleType	t,
		const G4int	A,
		const G4int	Z
	)

Definition at line 140 of file G4INCLHFB.cc.

                                                                                         {
        // HFB calculations
      G4double a=0.;
        if(t==Neutron)
          if(diffusenessN[Z][A]>0.)a=diffusenessN[Z][A];
        if(t==Proton)
          if(diffusenessP[Z][A]>0.)a=diffusenessP[Z][A];
      return a;
    }

References A, G4INCL::anonymous_namespace{G4INCLHFB.cc}::diffusenessN, G4INCL::anonymous_namespace{G4INCLHFB.cc}::diffusenessP, G4INCL::Neutron, G4INCL::Proton, and Z.

Referenced by G4INCL::ParticleTable::getRadiusParameter(), and G4INCL::ParticleTable::getSurfaceDiffuseness().

initialize()

void G4INCL::HFB::initialize

(

)

Definition at line 81 of file G4INCLHFB.cc.

                        {
#else
      void initialize(const std::string &path) {
#endif
 
      // Clear the existing tables, if any
      cleanTable();
 
#ifdef INCLXX_IN_GEANT4_MODE
       if(!std::getenv("G4INCLDATA")) {
        G4ExceptionDescription ed;
        ed << " Data missing: set environment variable G4INCLDATA\n"
           << " to point to the directory containing data files needed\n"
           << " by the INCL++ model" << G4endl;
           G4Exception("G4INCLDataFile::readData()","table_radius_hfb.dat",
                FatalException, ed);
      }
      G4String dataPath0(std::getenv("G4INCLDATA"));
      G4String dataPath(dataPath0 + "/table_radius_hfb.dat");
#else
      // File name
      std::string dataPath(path + "/table_radius_hfb.dat");
      INCL_DEBUG("Reading radius and diffuseness parameters from file " << dataPath << '\n');
#endif
 
      // Open the file stream
      std::ifstream hfbTableIn(dataPath.c_str());
      if(!hfbTableIn.good()) {
        std::cerr << "Cannot open " << dataPath << " data file." << '\n';
        std::abort();
        return;
      }
 
      // read the file
      G4int z, a, nbnuclei=0;
      G4double rp, rn, dp, dn;
      while(hfbTableIn.good()) { /* Loop checking, 22.01.2018, J.L. Rodriguez */
        hfbTableIn >> z >> a >> rp >> rn >> dp >> dn;
        radiusP[z][a] = rp;
        radiusN[z][a] = rn;
        diffusenessP[z][a] = dp;
        diffusenessN[z][a] = dn;
        nbnuclei++;
      }
      hfbTableIn.close();
      INCL_DEBUG("Read " << nbnuclei << " nuclei" << '\n');
 
    }

References G4INCL::anonymous_namespace{G4INCLHFB.cc}::cleanTable(), G4INCL::anonymous_namespace{G4INCLHFB.cc}::diffusenessN, G4INCL::anonymous_namespace{G4INCLHFB.cc}::diffusenessP, FatalException, G4endl, G4Exception(), INCL_DEBUG, initialize(), G4INCL::anonymous_namespace{G4INCLHFB.cc}::radiusN, and G4INCL::anonymous_namespace{G4INCLHFB.cc}::radiusP.

Referenced by initialize(), and G4INCL::ParticleTable::initialize().