G4INCL::NuclearPotential Namespace Reference

Data Structures
class	INuclearPotential
class	NuclearPotentialConstant
class	NuclearPotentialEnergyIsospin
class	NuclearPotentialEnergyIsospinSmooth
class	NuclearPotentialIsospin

Functions
INuclearPotential const *	createPotential (const PotentialType type, const G4int theA, const G4int theZ, const G4bool pionPotential)
	Create an INuclearPotential object. More...
void	clearCache ()
	Clear the INuclearPotential cache. More...

Function Documentation

void G4INCL::NuclearPotential::clearCache

(

)

Clear the INuclearPotential cache.

Definition at line 98 of file G4INCLINuclearPotential.cc.

Referenced by G4INCL::INCL::~INCL().

                      {
      if(nuclearPotentialCache) {
        for(std::map<long,INuclearPotential const *>::const_iterator i = nuclearPotentialCache->begin(), e=nuclearPotentialCache->end(); i!=e; ++i)
          delete i->second;
        nuclearPotentialCache->clear();
        delete nuclearPotentialCache;
        nuclearPotentialCache = NULL;
      }
    }

INuclearPotential const * G4INCL::NuclearPotential::createPotential	(	const PotentialType	type,
		const G4int	theA,
		const G4int	theZ,
		const G4bool	pionPotential
	)

Create an INuclearPotential object.

This is the method that should be used to instantiate objects derived from INuclearPotential. It uses a caching mechanism to minimise thrashing and speed up the code.

Parameters

type	the type of the potential to be created
theA	mass number of the nucleus
theZ	charge number of the nucleus
pionPotential	whether pions should also feel the potential

Returns

a pointer to the nuclear potential

Definition at line 66 of file G4INCLINuclearPotential.cc.

References G4INCL::ConstantPotential, INCL_FATAL, G4INCL::IsospinEnergyPotential, G4INCL::IsospinEnergySmoothPotential, and G4INCL::IsospinPotential.

Referenced by G4INCL::Nucleus::Nucleus().

                                                                                                                                       {
      if(!nuclearPotentialCache)
        nuclearPotentialCache = new std::map<long,INuclearPotential const *>;

      const long nuclideID = (pionPotential ? 1 : -1) * (1000*theZ + theA + 1000000*type); // MCNP-style nuclide IDs
      const std::map<long,INuclearPotential const *>::const_iterator mapEntry = nuclearPotentialCache->find(nuclideID);
      if(mapEntry == nuclearPotentialCache->end()) {
        INuclearPotential const *thePotential = NULL;
        switch(type) {
          case IsospinEnergySmoothPotential:
            thePotential = new NuclearPotentialEnergyIsospinSmooth(theA, theZ, pionPotential);
            break;
          case IsospinEnergyPotential:
            thePotential = new NuclearPotentialEnergyIsospin(theA, theZ, pionPotential);
            break;
          case IsospinPotential:
            thePotential = new NuclearPotentialIsospin(theA, theZ, pionPotential);
            break;
          case ConstantPotential:
            thePotential = new NuclearPotentialConstant(theA, theZ, pionPotential);
            break;
          default:
            INCL_FATAL("Unrecognized potential type at Nucleus creation." << std::endl);
            break;
        }
        (*nuclearPotentialCache)[nuclideID] = thePotential;
        return thePotential;
      } else {
        return mapEntry->second;
      }
    }