G4SFDecay Class Reference

#include <G4SFDecay.hh>

Inheritance diagram for G4SFDecay:

Public Member Functions
	G4SFDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &Qvalue, const G4double &excitation, const G4Ions::G4FloatLevelBase &flb)
virtual	~G4SFDecay ()
virtual G4DecayProducts *	DecayIt (G4double)
virtual void	DumpNuclearInfo ()
Public Member Functions inherited from G4NuclearDecay
	G4NuclearDecay (const G4String &channelName, const G4RadioactiveDecayMode &mode, const G4double &excitation, const G4Ions::G4FloatLevelBase &floatingLevel)
	~G4NuclearDecay () override=default
G4bool	IsOKWithParentMass (G4double parentMass) override
G4RadioactiveDecayMode	GetDecayMode () const
G4double	GetDaughterExcitation () const
G4Ions::G4FloatLevelBase	GetFloatingLevel () const
G4ParticleDefinition *	GetDaughterNucleus ()
Public Member Functions inherited from G4VDecayChannel
	G4VDecayChannel (const G4String &aName, G4int Verbose=1)
	G4VDecayChannel (const G4String &aName, const G4String &theParentName, G4double theBR, G4int theNumberOfDaughters, const G4String &theDaughterName1, const G4String &theDaughterName2="", const G4String &theDaughterName3="", const G4String &theDaughterName4="", const G4String &theDaughterName5="")
virtual	~G4VDecayChannel ()
G4bool	operator== (const G4VDecayChannel &r) const
G4bool	operator!= (const G4VDecayChannel &r) const
G4bool	operator< (const G4VDecayChannel &right) const
const G4String &	GetKinematicsName () const
G4double	GetBR () const
G4int	GetNumberOfDaughters () const
G4ParticleDefinition *	GetParent ()
G4ParticleDefinition *	GetDaughter (G4int anIndex)
G4int	GetAngularMomentum ()
const G4String &	GetParentName () const
const G4String &	GetDaughterName (G4int anIndex) const
G4double	GetParentMass () const
G4double	GetDaughterMass (G4int anIndex) const
void	SetParent (const G4ParticleDefinition *particle_type)
void	SetParent (const G4String &particle_name)
void	SetBR (G4double value)
void	SetNumberOfDaughters (G4int value)
void	SetDaughter (G4int anIndex, const G4ParticleDefinition *particle_type)
void	SetDaughter (G4int anIndex, const G4String &particle_name)
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
void	DumpInfo ()
G4double	GetRangeMass () const
void	SetRangeMass (G4double val)
void	SetPolarization (const G4ThreeVector &)
const G4ThreeVector &	GetPolarization () const

Additional Inherited Members
Protected Member Functions inherited from G4VDecayChannel
	G4VDecayChannel ()
	G4VDecayChannel (const G4VDecayChannel &)
G4VDecayChannel &	operator= (const G4VDecayChannel &)
void	ClearDaughtersName ()
void	CheckAndFillDaughters ()
void	CheckAndFillParent ()
G4double	DynamicalMass (G4double massPDG, G4double width, G4double maxDev=1.0) const
Protected Attributes inherited from G4NuclearDecay
const G4RadioactiveDecayMode	theMode
Protected Attributes inherited from G4VDecayChannel
G4String	kinematics_name = ""
G4double	rbranch = 0.0
G4String *	parent_name = nullptr
G4String **	daughters_name = nullptr
G4double	rangeMass = 2.5
G4ThreeVector	parent_polarization
G4ParticleTable *	particletable = nullptr
G4ParticleDefinition *	G4MT_parent = nullptr
G4ParticleDefinition **	G4MT_daughters = nullptr
G4double	G4MT_parent_mass = 0.0
G4double *	G4MT_daughters_mass = nullptr
G4double *	G4MT_daughters_width = nullptr
G4Mutex	daughtersMutex
G4Mutex	parentMutex
G4int	numberOfDaughters = 0
G4int	verboseLevel = 1
Static Protected Attributes inherited from G4VDecayChannel
static const G4String	noName = " "

Detailed Description

Definition at line 44 of file G4SFDecay.hh.

Constructor & Destructor Documentation

G4SFDecay()

G4SFDecay::G4SFDecay	(	const G4ParticleDefinition *	theParentNucleus,
		const G4double &	theBR,
		const G4double &	Qvalue,
		const G4double &	excitation,
		const G4Ions::G4FloatLevelBase &	flb )

Definition at line 47 of file G4SFDecay.cc.

 : G4NuclearDecay("SF decay", SpFission, excitationE, flb), transitionQ(Qvalue)
{
  SetParent(theParentNucleus);  // Store name of parent nucleus, delete G4MT_parent 
  SetBR(branch);
 
  parentZ = theParentNucleus->GetAtomicNumber();
  parentA = theParentNucleus->GetAtomicMass(); 
 
//  SetNumberOfDaughters(0);  doesn't seem to work
//  Set to 1 for now
//  When we later install proper fragment generation, set this to 2
//  and modify Z and A in DecayIt() 
  SetNumberOfDaughters(1);
  G4IonTable* theIonTable =
    (G4IonTable*)(G4ParticleTable::GetParticleTable()->GetIonTable());
  SetDaughter(0, theIonTable->GetIon(parentZ, parentA, excitationE, flb) );
}

~G4SFDecay()

G4SFDecay::~G4SFDecay ( )

virtual

Definition at line 70 of file G4SFDecay.cc.

{}

Member Function Documentation

DecayIt()

G4DecayProducts * G4SFDecay::DecayIt ( G4double )

virtual

Implements G4VDecayChannel.

Definition at line 74 of file G4SFDecay.cc.

{
  // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)  
  CheckAndFillParent();
 
  // Set up final state
  // parentParticle is set at rest here because boost with correct momentum 
  // is done later
  G4LorentzVector atRest(G4MT_parent->GetPDGMass(),
                         G4ThreeVector(0.,0.,0.) );
  G4DynamicParticle parentParticle(G4MT_parent, atRest);
  G4DecayProducts* products = new G4DecayProducts(parentParticle);
 
  // Generate LLNL code for parent nucleus
  G4int code = 1000*G4MT_parent->GetAtomicNumber() +
                    G4MT_parent->GetAtomicMass();
  
  // Let G4fissionEvent do the decay
  // Argument 2 : time, passed through to set time of secondary
  // Argument 3 : nubar, when set to -1, selects spontaneous fission
  // Argument 4 : energy used in Watt spectrum 
  G4fissionEvent aFission(code, 10.0, -1, 0.);
  G4int nNeut = aFission.getNeutronNu();
  G4int nGam = aFission.getPhotonNu();
 
  G4DynamicParticle* dynPart = 0; 
  G4ThreeVector direction(0.,0.,0.);
  G4double KE = 0.;   // energy from generator is in MeV 
 
  if (nNeut > 0) {
    // Put neutrons on the stack
    for (G4int i = 0; i < nNeut; i++) {
      KE = aFission.getNeutronEnergy(i);
      direction.setX(aFission.getNeutronDircosu(i) );
      direction.setY(aFission.getNeutronDircosv(i) );
      direction.setZ(aFission.getNeutronDircosw(i) );
 
      dynPart = new G4DynamicParticle(G4Neutron::Neutron(), direction, KE);
    
      // Not needed - time comes from parent track
      // dynPart->SetProperTime(aFission.getNeutronAge(i) );
  
      products->PushProducts(dynPart);
    }
 
    // Put gammas on the stack
    for (G4int i = 0; i < nGam; i++) {
      KE = aFission.getPhotonEnergy(i);
      direction.setX(aFission.getPhotonDircosu(i) );
      direction.setY(aFission.getPhotonDircosv(i) );
      direction.setZ(aFission.getPhotonDircosw(i) );
 
      dynPart = new G4DynamicParticle(G4Gamma::Gamma(), direction, KE);
 
      // Not needed - time comes from parent track
      // dynPart->SetProperTime(aFission.getPhotonAge(i) );
 
      products->PushProducts(dynPart);
      // No residual nucleus in this model
    }
 
  } else {
    // No data for this isotope, return parent nucleus for now
//    G4cout << " No fission data for this isotope " << G4endl;
    G4DynamicParticle* parent =
      new G4DynamicParticle(G4MT_parent, G4ThreeVector(0.,0.,0.) );
    products->PushProducts(parent);
  } 
 
  // Energy conservation check not valid in this model
 
  return products;
}

DumpNuclearInfo()

void G4SFDecay::DumpNuclearInfo ( )

virtual

Implements G4NuclearDecay.

Definition at line 149 of file G4SFDecay.cc.

{
  G4cout << " G4SFDecay for parent nucleus " << GetParentName() << G4endl;
  G4cout << " decays to neutrons and gammas, with branching ratio " << GetBR()
         << "% and Q value " << transitionQ << G4endl;
}

---

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

• G4SFDecay.hh
• G4SFDecay.cc