G4INCL::ClusterDecay Namespace Reference

Namespace for functions that handle decay of unstable clusters. More...

Enumerations
enum	ClusterDecayType {   StableCluster, NeutronDecay, ProtonDecay, AlphaDecay,   TwoProtonDecay, TwoNeutronDecay, ProtonUnbound, NeutronUnbound }

Functions
G4bool	isStable (Cluster const *const c)
	True if the cluster is stable. More...
ParticleList	decay (Cluster *const c)
	Carries out a cluster decay. More...

Variables
G4ThreadLocal ClusterDecayType	clusterDecayMode [ParticleTable::clusterTableZSize][ParticleTable::clusterTableASize]
	Table for cluster decays. More...

Detailed Description

Namespace for functions that handle decay of unstable clusters.

Enumeration Type Documentation

enum G4INCL::ClusterDecay::ClusterDecayType

Enumerator
StableCluster
NeutronDecay
ProtonDecay
AlphaDecay
TwoProtonDecay
TwoNeutronDecay
ProtonUnbound
NeutronUnbound

Definition at line 65 of file G4INCLClusterDecay.hh.

                            {
        StableCluster,
        NeutronDecay,
        ProtonDecay,
        AlphaDecay,
        TwoProtonDecay,
        TwoNeutronDecay,
        ProtonUnbound,
        NeutronUnbound
      };

Function Documentation

ParticleList G4INCL::ClusterDecay::decay

(

Cluster *const

c

)

Carries out a cluster decay.

Parameters

c	cluster that should decay

Returns

list of decay products

Definition at line 551 of file G4INCLClusterDecay.cc.

References G4INCL::Particle::getA(), G4INCL::Particle::getZ(), G4INCL::Neutron, G4INCL::Proton, G4INCL::Particle::setRealMass(), and G4INCL::Particle::setType().

Referenced by G4MolecularDecayTable::CheckDataConsistency(), F05PhysicsList::ConstructProcess(), P6DExtDecayerPhysics::ConstructProcess(), F04PhysicsList::ConstructProcess(), WLSPhysicsList::ConstructProcess(), G4INCL::Nucleus::decayInsideDeltas(), G4INCL::Nucleus::decayMe(), G4INCL::Nucleus::decayOutgoingClusters(), G4INCL::Nucleus::decayOutgoingDeltas(), G4DecayStrongResonances::Propagate(), G4GeneratorPrecompoundInterface::Propagate(), G4BinaryCascade::Propagate(), G4CascadeInterface::Propagate(), and G4GeneratorPrecompoundInterface::PropagateNuclNucl().

                                          {
      ParticleList decayProducts;
      recursiveDecay(c, &decayProducts);

      // Correctly update the particle type
      if(c->getA()==1) {
// assert(c->getZ()==1 || c->getZ()==0);
        if(c->getZ()==1)
          c->setType(Proton);
        else
          c->setType(Neutron);
        c->setRealMass();
      }

      return decayProducts;
    }

G4bool G4INCL::ClusterDecay::isStable

(

Cluster const *const

c

)

True if the cluster is stable.

Definition at line 521 of file G4INCLClusterDecay.cc.

References clusterDecayMode, G4INCL::Particle::getA(), G4INCL::Particle::getZ(), and StableCluster.

                                             {
      const G4int Z = c->getZ();
      const G4int A = c->getA();
      return (clusterDecayMode[Z][A]==StableCluster);
    }

Variable Documentation

G4ThreadLocal ClusterDecayType G4INCL::ClusterDecay::clusterDecayMode

Initial value:

=
    {
      
          {StableCluster, StableCluster,   NeutronDecay, NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound, NeutronUnbound},
          {StableCluster, StableCluster,  StableCluster,  StableCluster,   NeutronDecay, TwoNeutronDecay,   NeutronDecay, TwoNeutronDecay, NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound, NeutronUnbound},
          {StableCluster, StableCluster,    ProtonDecay,  StableCluster,  StableCluster,    NeutronDecay,  StableCluster,    NeutronDecay,  StableCluster,   NeutronDecay, TwoNeutronDecay, NeutronUnbound, NeutronUnbound},
          {StableCluster, StableCluster,  StableCluster,  ProtonUnbound,    ProtonDecay,     ProtonDecay,  StableCluster,   StableCluster,  StableCluster,  StableCluster,    NeutronDecay,  StableCluster,   NeutronDecay},
          {StableCluster, StableCluster,  StableCluster,  StableCluster,  ProtonUnbound,     ProtonDecay, TwoProtonDecay,   StableCluster,     AlphaDecay,  StableCluster,   StableCluster,  StableCluster,  StableCluster},
          {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   ProtonUnbound, TwoProtonDecay,     ProtonDecay,  StableCluster,    ProtonDecay,   StableCluster,  StableCluster,  StableCluster},
          {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   StableCluster,  ProtonUnbound,   ProtonUnbound, TwoProtonDecay,  StableCluster,   StableCluster,  StableCluster,  StableCluster},
          {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   StableCluster,  StableCluster,   ProtonUnbound,  ProtonUnbound,  ProtonUnbound,     ProtonDecay,    ProtonDecay,  StableCluster},
          {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   StableCluster,  StableCluster,   StableCluster,  ProtonUnbound,  ProtonUnbound,   ProtonUnbound,  ProtonUnbound,    ProtonDecay}
    }

Table for cluster decays.

Definition of "Stable": halflife > 1 ms

These table includes decay data for clusters that INCL presently does not produce. It can't hurt.

Unphysical nuclides (A<Z) are marked as stable, but should never be produced by INCL. If you find them in the output, something is fishy.

Definition at line 537 of file G4INCLClusterDecay.cc.

Referenced by isStable().