G4InuclNuclei Class Reference

#include <G4InuclNuclei.hh>

Inheritance diagram for G4InuclNuclei:

Public Member Functions
	G4InuclNuclei ()
	G4InuclNuclei (const G4DynamicParticle &dynPart, Model model=DefaultModel)
	G4InuclNuclei (G4int a, G4int z, G4double exc=0., Model model=DefaultModel)
	G4InuclNuclei (const G4LorentzVector &mom, G4int a, G4int z, G4double exc=0., Model model=DefaultModel)
	G4InuclNuclei (G4double ekin, G4int a, G4int z, G4double exc, Model model=DefaultModel)
	G4InuclNuclei (const G4Fragment &aFragment, Model model=DefaultModel)
	G4InuclNuclei (G4V3DNucleus *a3DNucleus, Model model=DefaultModel)
virtual	~G4InuclNuclei ()
	G4InuclNuclei (const G4InuclNuclei &right)
G4InuclNuclei &	operator= (const G4InuclNuclei &right)
G4bool	operator== (const G4InuclNuclei &right)
void	fill (G4int a, G4int z, G4double exc=0., Model model=DefaultModel)
void	fill (const G4LorentzVector &mom, G4int a, G4int z, G4double exc=0., Model model=DefaultModel)
void	fill (G4double ekin, G4int a, G4int z, G4double exc, Model model=DefaultModel)
void	copy (const G4Fragment &aFragment, Model model=DefaultModel)
void	copy (G4V3DNucleus *a3DNucleus, Model model=DefaultModel)
void	clear ()
void	setExitationEnergy (G4double e)
void	setExitonConfiguration (const G4ExitonConfiguration &config)
void	clearExitonConfiguration ()
G4int	getA () const
G4int	getZ () const
G4double	getNucleiMass () const
G4double	getExitationEnergy () const
G4double	getExitationEnergyInGeV () const
const G4ExitonConfiguration &	getExitonConfiguration () const
virtual void	print (std::ostream &os) const
G4Fragment	makeG4Fragment () const
	operator G4Fragment () const
Public Member Functions inherited from G4InuclParticle
	G4InuclParticle ()
	G4InuclParticle (const G4DynamicParticle &dynPart, Model model=DefaultModel)
	G4InuclParticle (const G4LorentzVector &mom, Model model=DefaultModel)
virtual	~G4InuclParticle ()
	G4InuclParticle (const G4InuclParticle &right)
G4InuclParticle &	operator= (const G4InuclParticle &right)
G4bool	operator== (const G4InuclParticle &right)
G4bool	operator!= (const G4InuclParticle &right)
void	setEnergy ()
void	setMomentum (const G4LorentzVector &mom)
void	setKineticEnergy (G4double ekin)
void	setMass (G4double mass)
G4double	getMass () const
G4double	getCharge () const
G4double	getKineticEnergy () const
G4double	getEnergy () const
G4double	getMomModule () const
G4LorentzVector	getMomentum () const
virtual void	print (std::ostream &os) const
const G4ParticleDefinition *	getDefinition () const
const G4DynamicParticle &	getDynamicParticle () const
void	setModel (Model model)
Model	getModel () const

Static Public Member Functions
static G4double	getNucleiMass (G4int a, G4int z, G4double exc=0.)

Static Protected Member Functions
static G4ParticleDefinition *	makeDefinition (G4int a, G4int z)
static G4ParticleDefinition *	makeNuclearFragment (G4int a, G4int z)

Additional Inherited Members
Public Types inherited from G4InuclParticle
enum	Model {   DefaultModel , bullet , target , EPCollider ,   INCascader , NonEquilib , Equilib , Fissioner ,   BigBanger , PreCompound , Coalescence }
Protected Member Functions inherited from G4InuclParticle
	G4InuclParticle (const G4ParticleDefinition *pd, Model model=DefaultModel)
	G4InuclParticle (const G4ParticleDefinition *pd, const G4LorentzVector &mom, Model model=DefaultModel)
	G4InuclParticle (const G4ParticleDefinition *pd, G4double ekin, Model model=DefaultModel)
void	setDefinition (const G4ParticleDefinition *pd)

Detailed Description

Definition at line 67 of file G4InuclNuclei.hh.

Constructor & Destructor Documentation

G4InuclNuclei() [1/8]

G4InuclNuclei::G4InuclNuclei ( )

inline

Definition at line 69 of file G4InuclNuclei.hh.

: G4InuclParticle() {}

G4InuclNuclei() [2/8]

G4InuclNuclei::G4InuclNuclei ( const G4DynamicParticle &  dynPart, Model  model = DefaultModel  )

inline

Definition at line 71 of file G4InuclNuclei.hh.

    : G4InuclParticle(dynPart, model) {}

G4InuclNuclei() [3/8]

G4InuclNuclei::G4InuclNuclei ( G4int  a, G4int  z, G4double  exc = 0., Model  model = DefaultModel  )

inline

Definition at line 74 of file G4InuclNuclei.hh.

    : G4InuclParticle(makeDefinition(a,z), model) {
    setExitationEnergy(exc);
  }

G4InuclNuclei() [4/8]

G4InuclNuclei::G4InuclNuclei ( const G4LorentzVector &  mom, G4int  a, G4int  z, G4double  exc = 0., Model  model = DefaultModel  )

inline

Definition at line 79 of file G4InuclNuclei.hh.

    : G4InuclParticle(makeDefinition(a,z), mom, model) {
    setExitationEnergy(exc);
  }

G4InuclNuclei() [5/8]

G4InuclNuclei::G4InuclNuclei ( G4double  ekin, G4int  a, G4int  z, G4double  exc, Model  model = DefaultModel  )

inline

Definition at line 85 of file G4InuclNuclei.hh.

    : G4InuclParticle(makeDefinition(a,z), ekin, model) {
    setExitationEnergy(exc);
  }

G4InuclNuclei() [6/8]

G4InuclNuclei::G4InuclNuclei	(	const G4Fragment &	aFragment,
		G4InuclParticle::Model	model = DefaultModel
	)

Definition at line 81 of file G4InuclNuclei.cc.

  : G4InuclParticle() {
  copy(aFragment, model);
}

G4InuclNuclei() [7/8]

G4InuclNuclei::G4InuclNuclei	(	G4V3DNucleus *	a3DNucleus,
		G4InuclParticle::Model	model = DefaultModel
	)

Definition at line 127 of file G4InuclNuclei.cc.

  : G4InuclParticle() {
  copy(a3DNucleus, model);
}

~G4InuclNuclei()

virtual G4InuclNuclei::~G4InuclNuclei ( )

inlinevirtual

Definition at line 95 of file G4InuclNuclei.hh.

{}

G4InuclNuclei() [8/8]

G4InuclNuclei::G4InuclNuclei ( const G4InuclNuclei &  right )

inline

Definition at line 98 of file G4InuclNuclei.hh.

    : G4InuclParticle(right),
      theExitonConfiguration(right.theExitonConfiguration) {}

Member Function Documentation

clear()

void G4InuclNuclei::clear

(

)

Definition at line 176 of file G4InuclNuclei.cc.

                          {
  setDefinition(0);
  clearExitonConfiguration();
  setModel(G4InuclParticle::DefaultModel);
}

clearExitonConfiguration()

void G4InuclNuclei::clearExitonConfiguration ( )

inline

Definition at line 134 of file G4InuclNuclei.hh.

{ theExitonConfiguration.clear(); }

Referenced by clear(), and fill().

copy() [1/2]

void G4InuclNuclei::copy	(	const G4Fragment &	aFragment,
		Model	model = DefaultModel
	)

Definition at line 87 of file G4InuclNuclei.cc.

                                                                 {
  fill(aFragment.GetMomentum()/GeV, aFragment.GetA_asInt(),
       aFragment.GetZ_asInt(), aFragment.GetExcitationEnergy(), model);
 
  // Exciton configuration must be set by hand
  theExitonConfiguration.protonQuasiParticles = aFragment.GetNumberOfCharged();
 
  theExitonConfiguration.neutronQuasiParticles =
    aFragment.GetNumberOfParticles() - aFragment.GetNumberOfCharged();
 
  theExitonConfiguration.protonHoles = aFragment.GetNumberOfChargedHoles();
 
  theExitonConfiguration.neutronHoles =
    aFragment.GetNumberOfHoles() - theExitonConfiguration.protonHoles;
}

Referenced by G4InuclNuclei().

copy() [2/2]

void G4InuclNuclei::copy	(	G4V3DNucleus *	a3DNucleus,
		Model	model = DefaultModel
	)

Definition at line 133 of file G4InuclNuclei.cc.

                                                              {
  if (!a3DNucleus) return;      // Null pointer means no action
 
  fill(0., a3DNucleus->GetMassNumber(), a3DNucleus->GetCharge(), 0., model);
 
  // Convert every hit nucleon into an exciton hole
  if (a3DNucleus->StartLoop()) {
    G4Nucleon* nucl = 0;
 
    /* Loop checking 08.06.2015 MHK */
    while ((nucl = a3DNucleus->GetNextNucleon())) {
      if (nucl->AreYouHit()) {  // Found previously interacted nucleon
    if (nucl->GetParticleType() == G4Proton::Definition())
      theExitonConfiguration.protonHoles++;
 
    if (nucl->GetParticleType() == G4Neutron::Definition())
      theExitonConfiguration.neutronHoles++;
      }
    }
  }
}

fill() [1/3]

void G4InuclNuclei::fill	(	const G4LorentzVector &	mom,
		G4int	a,
		G4int	z,
		G4double	exc = 0.,
		G4InuclParticle::Model	model = DefaultModel
	)

Definition at line 158 of file G4InuclNuclei.cc.

                                                       {
  setDefinition(makeDefinition(a,z));
  setMomentum(mom);
  setExitationEnergy(exc);
  clearExitonConfiguration();
  setModel(model);
}

fill() [2/3]

void G4InuclNuclei::fill	(	G4double	ekin,
		G4int	a,
		G4int	z,
		G4double	exc,
		G4InuclParticle::Model	model = DefaultModel
	)

Definition at line 167 of file G4InuclNuclei.cc.

                                         {
  setDefinition(makeDefinition(a,z));
  setKineticEnergy(ekin);
  setExitationEnergy(exc);
  clearExitonConfiguration();
  setModel(model);
}

fill() [3/3]

void G4InuclNuclei::fill ( G4int  a, G4int  z, G4double  exc = 0., Model  model = DefaultModel  )

inline

Definition at line 111 of file G4InuclNuclei.hh.

                                                                         {
    fill(0., a, z, exc, model);
  }

Referenced by copy(), G4CascadeInterface::createBullet(), G4CascadeInterface::createTarget(), G4IntraNucleiCascader::createTarget(), fill(), G4CascadeRecoilMaker::makeRecoilNuclei(), and G4IntraNucleiCascader::releaseSecondary().

getA()

G4int G4InuclNuclei::getA ( ) const

inline

Definition at line 136 of file G4InuclNuclei.hh.

{ return getDefinition()->GetAtomicMass(); }

Referenced by G4InuclCollider::collide(), G4CascadeCheckBalance::collide(), G4CascadeColliderBase::inelasticInteractionPossible(), G4IntraNucleiCascader::initialize(), G4NucleiModel::initializeCascad(), makeG4Fragment(), print(), G4InteractionCase::set(), and G4IntraNucleiCascader::setupCascade().

getExitationEnergy()

G4double G4InuclNuclei::getExitationEnergy ( ) const

inline

Definition at line 143 of file G4InuclNuclei.hh.

                                      {
    return (getMass()-getNucleiMass())*CLHEP::GeV/CLHEP::MeV; // Always in MeV
  }

Referenced by getExitationEnergyInGeV(), and print().

getExitationEnergyInGeV()

G4double G4InuclNuclei::getExitationEnergyInGeV ( ) const

inline

Definition at line 147 of file G4InuclNuclei.hh.

{ return getExitationEnergy()/CLHEP::GeV; }

getExitonConfiguration()

const G4ExitonConfiguration & G4InuclNuclei::getExitonConfiguration ( ) const

inline

Definition at line 149 of file G4InuclNuclei.hh.

                                                              {
    return theExitonConfiguration;
  }

getNucleiMass() [1/2]

G4double G4InuclNuclei::getNucleiMass ( ) const

inline

Definition at line 139 of file G4InuclNuclei.hh.

                                 {
    return getDefinition()->GetPDGMass()*CLHEP::MeV/CLHEP::GeV; // From G4 to Bertini
  }

Referenced by G4EquilibriumEvaporator::deExcite(), G4Fissioner::deExcite(), G4NonEquilibriumEvaporator::deExcite(), G4CascadeRecoilMaker::deltaM(), getExitationEnergy(), G4CascadeDeexciteBase::makeFragment(), makeNuclearFragment(), G4CascadeRecoilMaker::makeRecoilFragment(), and setExitationEnergy().

getNucleiMass() [2/2]

G4double G4InuclNuclei::getNucleiMass ( G4int  a, G4int  z, G4double  exc = 0.  )

static

Definition at line 276 of file G4InuclNuclei.cc.

                                                                    {
  // Simple minded mass calculation use constants in CLHEP (all in MeV)
  G4double mass = G4NucleiProperties::GetNuclearMass(a,z) + exc;
 
  return mass*MeV/GeV;      // Convert from GEANT4 to Bertini units
}

getZ()

G4int G4InuclNuclei::getZ ( ) const

inline

Definition at line 137 of file G4InuclNuclei.hh.

{ return getDefinition()->GetAtomicNumber(); }

Referenced by G4InuclCollider::collide(), G4IntraNucleiCascader::generateCascade(), G4CascadeColliderBase::inelasticInteractionPossible(), G4IntraNucleiCascader::initialize(), G4NucleiModel::initializeCascad(), makeG4Fragment(), print(), and G4IntraNucleiCascader::setupCascade().

makeDefinition()

G4ParticleDefinition * G4InuclNuclei::makeDefinition ( G4int  a, G4int  z  )

staticprotected

Definition at line 204 of file G4InuclNuclei.cc.

                                                                    {
  // SPECIAL CASE:  (0,0) means create dummy without definition
  if (0 == a && 0 == z) return 0;
 
  G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
  G4ParticleDefinition *pd = pTable->GetIonTable()->GetIon(z, a, 0);
 
  // SPECIAL CASE:  Non-physical nuclear fragment, for final-state return
  if (!pd) pd = makeNuclearFragment(a,z);
 
  return pd;        // This could return a null pointer if above fails
}

Referenced by fill().

makeG4Fragment()

G4Fragment G4InuclNuclei::makeG4Fragment

(

)

const

Definition at line 105 of file G4InuclNuclei.cc.

                                               {
  G4Fragment frag(getA(), getZ(), getMomentum()*GeV);   // From Bertini units
 
  // Note:  exciton configuration has to be set piece by piece
  frag.SetNumberOfHoles(theExitonConfiguration.protonHoles
            + theExitonConfiguration.neutronHoles,
            theExitonConfiguration.protonHoles);
 
  frag.SetNumberOfExcitedParticle(theExitonConfiguration.protonQuasiParticles 
          + theExitonConfiguration.neutronQuasiParticles,
          theExitonConfiguration.protonQuasiParticles);
 
  return frag;
}

makeNuclearFragment()

G4ParticleDefinition * G4InuclNuclei::makeNuclearFragment ( G4int  a, G4int  z  )

staticprotected

Definition at line 229 of file G4InuclNuclei.cc.

                                                   {
  if (a<=0 || z<0 || a<z) {
    G4cerr << " >>> G4InuclNuclei::makeNuclearFragment() called with"
       << " impossible arguments A=" << a << " Z=" << z << G4endl;
    throw G4HadronicException(__FILE__, __LINE__,
                  "G4InuclNuclei impossible A/Z arguments");
  }
 
  G4int code = G4IonTable::GetNucleusEncoding(z, a);
 
  // Use local lookup table (see above) to maintain singletons
  // NOTE:  G4ParticleDefinitions don't need to be explicitly deleted
  //        (see comments in G4IonTable.cc::~G4IonTable)
 
  G4AutoLock fragListLock(&fragListMutex);
  if (fragmentList.find(code) != fragmentList.end()) return fragmentList[code];
  fragListLock.unlock();
 
  // Name string follows format in G4IonTable.cc::GetIonName(Z,A,E)
  std::stringstream zstr, astr;
  zstr << z;
  astr << a;
  
  G4String name = "Z" + zstr.str() + "A" + astr.str();
  
  G4double mass = getNucleiMass(a,z) *GeV/MeV;  // From Bertini to GEANT4 units
  
  //    Arguments for constructor are as follows
  //               name             mass          width         charge
  //             2*spin           parity  C-conjugation
  //          2*Isospin       2*Isospin3       G-parity
  //               type    lepton number  baryon number   PDG encoding
  //             stable         lifetime    decay table
  //             shortlived      subType    anti_encoding Excitation-energy
  
  G4Ions* fragPD = new G4Ions(name,       mass, 0., z*eplus,
                  0,          +1,   0,
                  0,          0,    0,
                  "nucleus",  0,    a, code,
                  true,   0.,   0,
                  true, "generic",  0,  0.);
  fragPD->SetAntiPDGEncoding(0);
 
  fragListLock.lock();          // Protect before saving new fragment
  return (fragmentList[code] = fragPD);     // Store in table for next lookup
}

Referenced by makeDefinition().

operator G4Fragment()

G4InuclNuclei::operator G4Fragment

(

)

const

Definition at line 120 of file G4InuclNuclei.cc.

                                         {
  return makeG4Fragment();
}

operator=()

G4InuclNuclei & G4InuclNuclei::operator=

(

const G4InuclNuclei &

right

)

Definition at line 284 of file G4InuclNuclei.cc.

                                                                  {
  if (this != &right) {
    theExitonConfiguration = right.theExitonConfiguration;
    G4InuclParticle::operator=(right);
  }
  return *this;
}

operator==()

G4bool G4InuclNuclei::operator== ( const G4InuclNuclei &  right )

inline

Definition at line 105 of file G4InuclNuclei.hh.

                                                {
    return ( G4InuclParticle::operator==(right) && 
         theExitonConfiguration == right.theExitonConfiguration );
  }

print()

void G4InuclNuclei::print ( std::ostream &  os ) const

virtual

Reimplemented from G4InuclParticle.

Definition at line 294 of file G4InuclNuclei.cc.

                                              {
  G4InuclParticle::print(os);
  os << G4endl << " Nucleus: " << getDefinition()->GetParticleName() 
     << " A " << getA() << " Z " << getZ() << " mass " << getMass()
     << " Eex (MeV) " << getExitationEnergy();
 
  if (!theExitonConfiguration.empty())
    os << G4endl << "         " << theExitonConfiguration;
}

setExitationEnergy()

void G4InuclNuclei::setExitationEnergy

(

G4double

e

)

Definition at line 185 of file G4InuclNuclei.cc.

                                                 {
  G4double ekin = getKineticEnergy();       // Current kinetic energy
 
  G4double emass = getNucleiMass() + e*MeV/GeV; // From Bertini to G4 units
 
  // Safety check -- if zero energy, don't do computation
  G4double ekin_new = (ekin == 0.) ? 0.
    : std::sqrt(emass*emass + ekin*(2.*getMass()+ekin)) - emass;
 
  setMass(emass);          // Momentum is computed from mass and Ekin
  setKineticEnergy(ekin_new);
}

Referenced by G4InuclEvaporation::BreakItUp(), fill(), and G4InuclNuclei().

setExitonConfiguration()

void G4InuclNuclei::setExitonConfiguration ( const G4ExitonConfiguration &  config )

inline

Definition at line 130 of file G4InuclNuclei.hh.

                                                                   { 
    theExitonConfiguration = config;
  }

Referenced by G4CascadeRecoilMaker::makeRecoilNuclei().

---

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

• G4InuclNuclei.hh
• G4InuclNuclei.cc