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 66 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 553 of file G4INCLClusterDecay.cc.

                                          {
      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;
    }

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4INCL::ClusterDecay::isStable

(

Cluster const *const

c

)

True if the cluster is stable.

Definition at line 523 of file G4INCLClusterDecay.cc.

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

Here is the call graph for this function:

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 539 of file G4INCLClusterDecay.cc.