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)
bool	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)
bool	operator== (const G4InuclParticle &right)
bool	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
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 (G4ParticleDefinition *pd, Model model=DefaultModel)
	G4InuclParticle (G4ParticleDefinition *pd, const G4LorentzVector &mom, Model model=DefaultModel)
	G4InuclParticle (G4ParticleDefinition *pd, G4double ekin, Model model=DefaultModel)
void	setDefinition (G4ParticleDefinition *pd)

Detailed Description

Definition at line 68 of file G4InuclNuclei.hh.

Constructor & Destructor Documentation

G4InuclNuclei() [1/8]

G4InuclNuclei::G4InuclNuclei ( )

inline

Definition at line 70 of file G4InuclNuclei.hh.

: G4InuclParticle() {}

G4InuclNuclei() [2/8]

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

inline

Definition at line 72 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 75 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 80 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 86 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 75 of file G4InuclNuclei.cc.

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

G4InuclNuclei() [7/8]

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

Definition at line 121 of file G4InuclNuclei.cc.

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

~G4InuclNuclei()

virtual G4InuclNuclei::~G4InuclNuclei ( )

inlinevirtual

Definition at line 96 of file G4InuclNuclei.hh.

{}

G4InuclNuclei() [8/8]

G4InuclNuclei::G4InuclNuclei ( const G4InuclNuclei &  right )

inline

Definition at line 99 of file G4InuclNuclei.hh.

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

Member Function Documentation

clear()

void G4InuclNuclei::clear

(

)

Definition at line 168 of file G4InuclNuclei.cc.

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

clearExitonConfiguration()

void G4InuclNuclei::clearExitonConfiguration ( )

inline

Definition at line 135 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 81 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 127 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;
    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 150 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 159 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 112 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 137 of file G4InuclNuclei.hh.

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

Referenced by G4PreCompoundDeexcitation::collide(), G4InuclCollider::collide(), G4BigBanger::collide(), G4CascadeCheckBalance::collide(), G4EquilibriumEvaporator::collide(), G4Fissioner::collide(), G4NonEquilibriumEvaporator::collide(), G4CascadeColliderBase::explosion(), G4CascadeColliderBase::inelasticInteractionPossible(), G4IntraNucleiCascader::initialize(), G4NucleiModel::initializeCascad(), makeG4Fragment(), print(), G4InteractionCase::set(), and G4IntraNucleiCascader::setupCascade().

getExitationEnergy()

G4double G4InuclNuclei::getExitationEnergy ( ) const

inline

Definition at line 144 of file G4InuclNuclei.hh.

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

Referenced by G4BigBanger::collide(), G4EquilibriumEvaporator::collide(), G4Fissioner::collide(), G4NonEquilibriumEvaporator::collide(), G4CascadeColliderBase::explosion(), getExitationEnergyInGeV(), and print().

getExitationEnergyInGeV()

G4double G4InuclNuclei::getExitationEnergyInGeV ( ) const

inline

Definition at line 148 of file G4InuclNuclei.hh.

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

getExitonConfiguration()

const G4ExitonConfiguration & G4InuclNuclei::getExitonConfiguration ( ) const

inline

Definition at line 150 of file G4InuclNuclei.hh.

                                                              {
    return theExitonConfiguration;
  }

Referenced by G4NonEquilibriumEvaporator::collide().

getNucleiMass() [1/2]

G4double G4InuclNuclei::getNucleiMass ( ) const

inline

Definition at line 140 of file G4InuclNuclei.hh.

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

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

getNucleiMass() [2/2]

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

static

Definition at line 260 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 138 of file G4InuclNuclei.hh.

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

Referenced by G4PreCompoundDeexcitation::collide(), G4InuclCollider::collide(), G4BigBanger::collide(), G4EquilibriumEvaporator::collide(), G4Fissioner::collide(), G4NonEquilibriumEvaporator::collide(), G4CascadeColliderBase::explosion(), 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 195 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->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 99 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 214 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 G4IonTable.hh) to maintain singletons
  // NOTE:  G4ParticleDefinitions don't need to be explicitly deleted
  //        (see comments in G4IonTable.cc::~G4IonTable)
 
  // If correct nucleus already created return it
  static std::map<G4int, G4ParticleDefinition*> fragmentList;
  if (fragmentList.find(code) != fragmentList.end()) return fragmentList[code];
 
  // 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);
 
  return (fragmentList[code] = fragPD);     // Store in table for next lookup
}

Referenced by makeDefinition().

operator G4Fragment()

G4InuclNuclei::operator G4Fragment

(

)

const

Definition at line 114 of file G4InuclNuclei.cc.

                                         {
  return makeG4Fragment();
}

operator=()

G4InuclNuclei & G4InuclNuclei::operator=

(

const G4InuclNuclei &

right

)

Definition at line 268 of file G4InuclNuclei.cc.

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

operator==()

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

inline

Definition at line 106 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 276 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 177 of file G4InuclNuclei.cc.

                                                 {
  G4double ekin = getKineticEnergy();       // Current kinetic energy
 
  G4double emass = getNucleiMass() + e*MeV/GeV; // From Bertini to G4 units
 
  // Directly compute new kinetic energy from old
  G4double ekin_new = 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 131 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