G4VDecayChannel Class Reference

#include <G4VDecayChannel.hh>

Inheritance diagram for G4VDecayChannel:

Public Member Functions
	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
virtual G4DecayProducts *	DecayIt (G4double parentMass=-1.0)=0
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)
virtual G4bool	IsOKWithParentMass (G4double parentMass)
void	SetPolarization (const G4ThreeVector &)
const G4ThreeVector &	GetPolarization () const

Protected Member Functions
	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
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
static const G4String	noName = " "

Detailed Description

Definition at line 49 of file G4VDecayChannel.hh.

Constructor & Destructor Documentation

G4VDecayChannel() [1/4]

G4VDecayChannel::G4VDecayChannel	(	const G4String &	aName,
		G4int	Verbose = 1 )

Definition at line 49 of file G4VDecayChannel.cc.

  : kinematics_name(aName), verboseLevel(verbose)
{
  // set pointer to G4ParticleTable (static and singleton object)
  particletable = G4ParticleTable::GetParticleTable();
}

G4VDecayChannel() [2/4]

G4VDecayChannel::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 = "" )

Definition at line 56 of file G4VDecayChannel.cc.

  : kinematics_name(aName), rbranch(theBR), numberOfDaughters(theNumberOfDaughters)
{
  // set pointer to G4ParticleTable (static and singleton object)
  particletable = G4ParticleTable::GetParticleTable();
 
  // parent name
  parent_name = new G4String(theParentName);
 
  // cleate array
  daughters_name = new G4String*[numberOfDaughters];
  for (G4int index = 0; index < numberOfDaughters; ++index) {
    daughters_name[index] = nullptr;
  }
 
  // daughters' name
  if (numberOfDaughters > 0) daughters_name[0] = new G4String(theDaughterName1);
  if (numberOfDaughters > 1) daughters_name[1] = new G4String(theDaughterName2);
  if (numberOfDaughters > 2) daughters_name[2] = new G4String(theDaughterName3);
  if (numberOfDaughters > 3) daughters_name[3] = new G4String(theDaughterName4);
  if (numberOfDaughters > 4) daughters_name[4] = new G4String(theDaughterName5);
 
  if (rbranch < 0.)
    rbranch = 0.0;
  else if (rbranch > 1.0)
    rbranch = 1.0;
}

~G4VDecayChannel()

G4VDecayChannel::~G4VDecayChannel ( )

virtual

Definition at line 154 of file G4VDecayChannel.cc.

{
  ClearDaughtersName();
  delete parent_name;
  parent_name = nullptr;
  delete[] G4MT_daughters_mass;
  G4MT_daughters_mass = nullptr;
  delete[] G4MT_daughters_width;
  G4MT_daughters_width = nullptr;
  G4MUTEXDESTROY(daughtersMutex);
  G4MUTEXDESTROY(parentMutex);
}

G4VDecayChannel() [3/4]

G4VDecayChannel::G4VDecayChannel ( )

protected

Definition at line 43 of file G4VDecayChannel.cc.

{
  // set pointer to G4ParticleTable (static and singleton object)
  particletable = G4ParticleTable::GetParticleTable();
}

G4VDecayChannel() [4/4]

G4VDecayChannel::G4VDecayChannel ( const G4VDecayChannel & right )

protected

Definition at line 88 of file G4VDecayChannel.cc.

{
  kinematics_name = right.kinematics_name;
  verboseLevel = right.verboseLevel;
  rbranch = right.rbranch;
  rangeMass = right.rangeMass;
 
  // copy parent name
  parent_name = new G4String(*right.parent_name);
 
  // create array
  numberOfDaughters = right.numberOfDaughters;
 
  daughters_name = nullptr;
  if (numberOfDaughters > 0) {
    daughters_name = new G4String*[numberOfDaughters];
    // copy daughters name
    for (G4int index = 0; index < numberOfDaughters; ++index) {
      daughters_name[index] = new G4String(*right.daughters_name[index]);
    }
  }
 
  // particle table
  particletable = G4ParticleTable::GetParticleTable();
 
  parent_polarization = right.parent_polarization;
 
  G4MUTEXINIT(daughtersMutex);
  G4MUTEXINIT(parentMutex);
}

Member Function Documentation

CheckAndFillDaughters()

void G4VDecayChannel::CheckAndFillDaughters ( )

inlineprotected

Definition at line 310 of file G4VDecayChannel.hh.

{
  G4AutoLock l(&daughtersMutex);
  if (G4MT_daughters == nullptr) {
    l.unlock();
    FillDaughters();
  }
}

Referenced by G4AlphaDecay::DecayIt(), G4DalitzDecayChannel::DecayIt(), G4ECDecay::DecayIt(), G4GeneralPhaseSpaceDecay::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4NeutronBetaDecayChannel::DecayIt(), G4NeutronDecay::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), G4PionRadiativeDecayChannel::DecayIt(), G4ProtonDecay::DecayIt(), G4TauLeptonicDecayChannel::DecayIt(), G4TritonDecay::DecayIt(), G4BetaMinusDecay::G4BetaMinusDecay(), G4BetaPlusDecay::G4BetaPlusDecay(), GetAngularMomentum(), GetDaughter(), G4PhaseSpaceDecayChannel::IsOKWithParentMass(), and IsOKWithParentMass().

CheckAndFillParent()

void G4VDecayChannel::CheckAndFillParent ( )

inlineprotected

Definition at line 319 of file G4VDecayChannel.hh.

{
  G4AutoLock l(&parentMutex);
  if (G4MT_parent == nullptr) {
    l.unlock();
    FillParent();
  }
}

Referenced by G4AlphaDecay::DecayIt(), G4DalitzDecayChannel::DecayIt(), G4ECDecay::DecayIt(), G4GeneralPhaseSpaceDecay::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4NeutronBetaDecayChannel::DecayIt(), G4NeutronDecay::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), G4PionRadiativeDecayChannel::DecayIt(), G4ProtonDecay::DecayIt(), G4SFDecay::DecayIt(), G4TauLeptonicDecayChannel::DecayIt(), G4TritonDecay::DecayIt(), G4BetaMinusDecay::G4BetaMinusDecay(), G4BetaPlusDecay::G4BetaPlusDecay(), GetParent(), G4PhaseSpaceDecayChannel::IsOKWithParentMass(), and IsOKWithParentMass().

ClearDaughtersName()

void G4VDecayChannel::ClearDaughtersName ( )

protected

Definition at line 167 of file G4VDecayChannel.cc.

{
  G4AutoLock l(&daughtersMutex);
  if (daughters_name != nullptr) {
    if (numberOfDaughters > 0) {
#ifdef G4VERBOSE
      if (verboseLevel > 1) {
        G4cout << "G4VDecayChannel::ClearDaughtersName() "
               << " for " << *parent_name << G4endl;
      }
#endif
      for (G4int index = 0; index < numberOfDaughters; ++index) {
        delete daughters_name[index];
      }
    }
    delete[] daughters_name;
    daughters_name = nullptr;
  }
 
  delete[] G4MT_daughters;
  delete[] G4MT_daughters_mass;
  delete[] G4MT_daughters_width;
  G4MT_daughters_width = nullptr;
  G4MT_daughters = nullptr;
  G4MT_daughters_mass = nullptr;
 
  numberOfDaughters = 0;
}

Referenced by G4DalitzDecayChannel::operator=(), G4KL3DecayChannel::operator=(), G4MuonDecayChannel::operator=(), G4MuonDecayChannelWithSpin::operator=(), G4MuonRadiativeDecayChannelWithSpin::operator=(), G4NeutronBetaDecayChannel::operator=(), G4PionRadiativeDecayChannel::operator=(), G4TauLeptonicDecayChannel::operator=(), operator=(), SetNumberOfDaughters(), and ~G4VDecayChannel().

DecayIt()

virtual G4DecayProducts * G4VDecayChannel::DecayIt ( G4double parentMass = -1.0 )

pure virtual

Implemented in G4AlphaDecay, G4BetaMinusDecay, G4BetaPlusDecay, G4DalitzDecayChannel, G4ECDecay, G4GeneralPhaseSpaceDecay, G4ITDecay, G4KL3DecayChannel, G4MuonDecayChannel, G4MuonDecayChannelWithSpin, G4MuonRadiativeDecayChannelWithSpin, G4NeutronBetaDecayChannel, G4NeutronDecay, G4PhaseSpaceDecayChannel, G4PionRadiativeDecayChannel, G4ProtonDecay, G4SFDecay, G4TauLeptonicDecayChannel, and G4TritonDecay.

Referenced by G4Decay::DecayIt(), G4Radioactivation::DecayIt(), G4IntraNucleiCascader::decayTrappedParticle(), and G4RadioactiveDecay::DoDecay().

DumpInfo()

void G4VDecayChannel::DumpInfo

(

)

Definition at line 484 of file G4VDecayChannel.cc.

{
  G4cout << " BR:  " << rbranch << "  [" << kinematics_name << "]";
  G4cout << "   :  ";
  for (G4int index = 0; index < numberOfDaughters; ++index) {
    if (daughters_name[index] != nullptr) {
      G4cout << " " << *(daughters_name[index]);
    }
    else {
      G4cout << " not defined ";
    }
  }
  G4cout << G4endl;
}

Referenced by G4GeneralPhaseSpaceDecay::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), and G4KL3DecayChannel::G4KL3DecayChannel().

DynamicalMass()

G4double G4VDecayChannel::DynamicalMass ( G4double massPDG, G4double width, G4double maxDev = 1.0 ) const

protected

Definition at line 504 of file G4VDecayChannel.cc.

{
  if (width <= 0.0) return massPDG;
  if (maxDev > rangeMass) maxDev = rangeMass;
  if (maxDev <= -1. * rangeMass) return massPDG;  // cannot calculate
 
  G4double x = G4UniformRand() * (maxDev + rangeMass) - rangeMass;
  G4double y = G4UniformRand();
  const std::size_t MAX_LOOP = 10000;
  for (std::size_t loop_counter = 0; loop_counter < MAX_LOOP; ++loop_counter) {
    if (y * (width * width * x * x + massPDG * massPDG * width * width)
        <= massPDG * massPDG * width * width)
      break;
    x = G4UniformRand() * (maxDev + rangeMass) - rangeMass;
    y = G4UniformRand();
  }
  G4double mass = massPDG + x * width;
  return mass;
}

GetAngularMomentum()

G4int G4VDecayChannel::GetAngularMomentum

(

)

Definition at line 428 of file G4VDecayChannel.cc.

{
  // determine angular momentum
 
  // fill pointers to daughter particles if not yet set
  CheckAndFillDaughters();
 
  const G4int PiSpin = G4MT_parent->GetPDGiSpin();
  const G4int PParity = G4MT_parent->GetPDGiParity();
  if (2 == numberOfDaughters)  // up to now we can only handle two particle decays
  {
    const G4int D1iSpin = G4MT_daughters[0]->GetPDGiSpin();
    const G4int D1Parity = G4MT_daughters[0]->GetPDGiParity();
    const G4int D2iSpin = G4MT_daughters[1]->GetPDGiSpin();
    const G4int D2Parity = G4MT_daughters[1]->GetPDGiParity();
    const G4int MiniSpin = std::abs(D1iSpin - D2iSpin);
    const G4int MaxiSpin = D1iSpin + D2iSpin;
    const G4int lMax = (PiSpin + D1iSpin + D2iSpin) / 2;  // l is always int
    G4int lMin;
#ifdef G4VERBOSE
    if (verboseLevel > 1) {
      G4cout << "iSpin: " << PiSpin << " -> " << D1iSpin << " + " << D2iSpin << G4endl;
      G4cout << "2*jmin, 2*jmax, lmax " << MiniSpin << " " << MaxiSpin << " " << lMax << G4endl;
    }
#endif
    for (G4int j = MiniSpin; j <= MaxiSpin; j += 2)  // loop over all possible
    {  // spin couplings
      lMin = std::abs(PiSpin - j) / 2;
#ifdef G4VERBOSE
      if (verboseLevel > 1) G4cout << "-> checking 2*j=" << j << G4endl;
#endif
      for (G4int l = lMin; l <= lMax; ++l) {
#ifdef G4VERBOSE
        if (verboseLevel > 1) G4cout << " checking l=" << l << G4endl;
#endif
        if (l % 2 == 0) {
          if (PParity == D1Parity * D2Parity)  // check parity for this l
            return l;
        }
        else {
          if (PParity == -1 * D1Parity * D2Parity)  // check parity for this l
            return l;
        }
      }
    }
  }
  else {
    G4Exception("G4VDecayChannel::GetAngularMomentum()", "PART111", JustWarning,
                "Sorry, can't handle 3 particle decays (up to now)");
    return 0;
  }
  G4Exception("G4VDecayChannel::GetAngularMomentum()", "PART111", JustWarning,
              "Can't find angular momentum for this decay");
  return 0;
}

GetBR()

G4double G4VDecayChannel::GetBR ( ) const

inline

Definition at line 275 of file G4VDecayChannel.hh.

{
  return rbranch;
}

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4Decay::DecayIt(), G4Radioactivation::DecayIt(), G4AlphaDecay::DumpNuclearInfo(), G4BetaMinusDecay::DumpNuclearInfo(), G4BetaPlusDecay::DumpNuclearInfo(), G4ECDecay::DumpNuclearInfo(), G4NeutronDecay::DumpNuclearInfo(), G4ProtonDecay::DumpNuclearInfo(), G4SFDecay::DumpNuclearInfo(), G4TritonDecay::DumpNuclearInfo(), G4KineticTrack::G4KineticTrack(), G4TextPPReporter::GeneratePropertyTable(), G4DecayTableMessenger::GetCurrentValue(), G4SampleResonance::GetMinimumMass(), G4DecayTable::Insert(), and G4RadioactiveDecay::LoadDecayTable().

GetDaughter()

G4ParticleDefinition * G4VDecayChannel::GetDaughter ( G4int anIndex )

inline

Definition at line 199 of file G4VDecayChannel.hh.

{
  // pointers to daughter particles are filled, if they are not set yet
  CheckAndFillDaughters();
 
  // get the pointer to a daughter particle
  if ((anIndex >= 0) && (anIndex < numberOfDaughters)) {
    return G4MT_daughters[anIndex];
  }
 
  if (verboseLevel > 0)
    G4cout << "G4VDecayChannel::GetDaughter  index out of range " << anIndex << G4endl;
  return nullptr;
}

Referenced by G4IonTable::CreateIon(), G4KineticTrack::Decay(), G4KineticTrack::G4KineticTrack(), G4TextPPReporter::GeneratePropertyTable(), G4NuclearDecay::GetDaughterNucleus(), G4MTRunManagerKernel::SetUpDecayChannels(), and G4TaskRunManagerKernel::SetUpDecayChannels().

GetDaughterMass()

G4double G4VDecayChannel::GetDaughterMass ( G4int anIndex ) const

inline

Definition at line 227 of file G4VDecayChannel.hh.

{
  if ((anIndex >= 0) && (anIndex < numberOfDaughters)) {
    return G4MT_daughters_mass[anIndex];
  }
 
  if (verboseLevel > 0) {
    G4cout << "G4VDecayChannel::GetDaughterMass ";
    G4cout << "index out of range " << anIndex << G4endl;
  }
  return 0.0;
}

GetDaughterName()

const G4String & G4VDecayChannel::GetDaughterName ( G4int anIndex ) const

inline

Definition at line 214 of file G4VDecayChannel.hh.

{
  if ((anIndex >= 0) && (anIndex < numberOfDaughters)) {
    return *daughters_name[anIndex];
  }
 
  if (verboseLevel > 0) {
    G4cout << "G4VDecayChannel::GetDaughterName ";
    G4cout << "index out of range " << anIndex << G4endl;
  }
  return GetNoName();
}

Referenced by G4KineticTrack::Decay(), G4Decay::DecayIt(), G4AlphaDecay::DumpNuclearInfo(), G4ECDecay::DumpNuclearInfo(), G4NeutronDecay::DumpNuclearInfo(), G4ProtonDecay::DumpNuclearInfo(), and G4TritonDecay::DumpNuclearInfo().

GetKinematicsName()

const G4String & G4VDecayChannel::GetKinematicsName ( ) const

inline

Definition at line 270 of file G4VDecayChannel.hh.

{
  return kinematics_name;
}

Referenced by G4TextPPReporter::GeneratePropertyTable().

GetNumberOfDaughters()

G4int G4VDecayChannel::GetNumberOfDaughters ( ) const

inline

Definition at line 265 of file G4VDecayChannel.hh.

{
  return numberOfDaughters;
}

Referenced by G4KineticTrack::Decay(), G4Decay::DecayIt(), G4KineticTrack::G4KineticTrack(), G4TextPPReporter::GeneratePropertyTable(), G4SampleResonance::GetMinimumMass(), and G4TextPPRetriever::ModifyPropertyTable().

GetParent()

G4ParticleDefinition * G4VDecayChannel::GetParent ( )

inline

Definition at line 240 of file G4VDecayChannel.hh.

{
  // the pointer to the parent particle is filled, if it is not set yet
  CheckAndFillParent();
  // get the pointer to the parent particle
  return G4MT_parent;
}

Referenced by G4DecayTable::Insert().

GetParentMass()

G4double G4VDecayChannel::GetParentMass ( ) const

inline

Definition at line 253 of file G4VDecayChannel.hh.

{
  return G4MT_parent_mass;
}

GetParentName()

const G4String & G4VDecayChannel::GetParentName ( ) const

inline

Definition at line 248 of file G4VDecayChannel.hh.

{
  return *parent_name;
}

Referenced by G4KineticTrack::Decay(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4AlphaDecay::DumpNuclearInfo(), G4ECDecay::DumpNuclearInfo(), G4NeutronDecay::DumpNuclearInfo(), G4ProtonDecay::DumpNuclearInfo(), G4SFDecay::DumpNuclearInfo(), and G4TritonDecay::DumpNuclearInfo().

GetPolarization()

const G4ThreeVector & G4VDecayChannel::GetPolarization ( ) const

inline

Definition at line 305 of file G4VDecayChannel.hh.

{
  return parent_polarization;
}

GetRangeMass()

G4double G4VDecayChannel::GetRangeMass ( ) const

inline

Definition at line 290 of file G4VDecayChannel.hh.

{
  return rangeMass;
}

GetVerboseLevel()

G4int G4VDecayChannel::GetVerboseLevel ( ) const

inline

Definition at line 285 of file G4VDecayChannel.hh.

{
  return verboseLevel;
}

Referenced by G4KL3DecayChannel::DalitzDensity(), G4DalitzDecayChannel::DecayIt(), G4Decay::DecayIt(), G4GeneralPhaseSpaceDecay::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4NeutronBetaDecayChannel::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), G4PionRadiativeDecayChannel::DecayIt(), G4TauLeptonicDecayChannel::DecayIt(), G4GeneralPhaseSpaceDecay::G4GeneralPhaseSpaceDecay(), G4GeneralPhaseSpaceDecay::G4GeneralPhaseSpaceDecay(), G4GeneralPhaseSpaceDecay::G4GeneralPhaseSpaceDecay(), G4GeneralPhaseSpaceDecay::G4GeneralPhaseSpaceDecay(), G4GeneralPhaseSpaceDecay::G4GeneralPhaseSpaceDecay(), G4KL3DecayChannel::G4KL3DecayChannel(), G4MuonDecayChannel::G4MuonDecayChannel(), G4MuonRadiativeDecayChannelWithSpin::G4MuonRadiativeDecayChannelWithSpin(), G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel(), G4PionRadiativeDecayChannel::G4PionRadiativeDecayChannel(), G4TauLeptonicDecayChannel::G4TauLeptonicDecayChannel(), G4GeneralPhaseSpaceDecay::ManyBodyDecayIt(), G4GeneralPhaseSpaceDecay::OneBodyDecayIt(), G4KL3DecayChannel::PhaseSpace(), G4GeneralPhaseSpaceDecay::ThreeBodyDecayIt(), and G4GeneralPhaseSpaceDecay::TwoBodyDecayIt().

IsOKWithParentMass()

G4bool G4VDecayChannel::IsOKWithParentMass ( G4double parentMass )

virtual

Reimplemented in G4NuclearDecay, and G4PhaseSpaceDecayChannel.

Definition at line 524 of file G4VDecayChannel.cc.

{
  G4double sumOfDaughterMassMin = 0.0;
  CheckAndFillParent();
  CheckAndFillDaughters();
  // skip one body decay
  if (numberOfDaughters == 1) return true;
 
  for (G4int index = 0; index < numberOfDaughters; ++index) {
    sumOfDaughterMassMin += G4MT_daughters_mass[index] - rangeMass * G4MT_daughters_width[index];
  }
  return (parentMass >= sumOfDaughterMassMin);
}

Referenced by G4Decay::DecayIt(), and G4PhaseSpaceDecayChannel::IsOKWithParentMass().

operator!=()

G4bool G4VDecayChannel::operator!= ( const G4VDecayChannel & r ) const

inline

Definition at line 64 of file G4VDecayChannel.hh.

{ return (this != &r); }

operator<()

G4bool G4VDecayChannel::operator< ( const G4VDecayChannel & right ) const

inline

Definition at line 194 of file G4VDecayChannel.hh.

{
  return (this->rbranch < right.rbranch);
}

operator=()

G4VDecayChannel & G4VDecayChannel::operator= ( const G4VDecayChannel & right )

protected

Definition at line 119 of file G4VDecayChannel.cc.

{
  if (this != &right) {
    kinematics_name = right.kinematics_name;
    verboseLevel = right.verboseLevel;
    rbranch = right.rbranch;
    rangeMass = right.rangeMass;
    parent_polarization = right.parent_polarization;
    // copy parent name
    delete parent_name;
    parent_name = new G4String(*right.parent_name);
 
    // clear daughters_name array
    ClearDaughtersName();
 
    // recreate array
    numberOfDaughters = right.numberOfDaughters;
    if (numberOfDaughters > 0) {
      daughters_name = new G4String*[numberOfDaughters];
      // copy daughters name
      for (G4int index = 0; index < numberOfDaughters; ++index) {
        daughters_name[index] = new G4String(*right.daughters_name[index]);
      }
    }
  }
 
  // particle table
  particletable = G4ParticleTable::GetParticleTable();
 
  G4MUTEXINIT(daughtersMutex);
  G4MUTEXINIT(parentMutex);
 
  return *this;
}

operator==()

G4bool G4VDecayChannel::operator== ( const G4VDecayChannel & r ) const

inline

Definition at line 63 of file G4VDecayChannel.hh.

{ return (this == &r); }

SetBR()

void G4VDecayChannel::SetBR

(

G4double

value

)

Definition at line 538 of file G4VDecayChannel.cc.

{
  rbranch = value;
  if (rbranch < 0.)
    rbranch = 0.0;
  else if (rbranch > 1.0)
    rbranch = 1.0;
}

Referenced by G4TauMinus::Definition(), G4TauPlus::Definition(), G4AlphaDecay::G4AlphaDecay(), G4BetaMinusDecay::G4BetaMinusDecay(), G4BetaPlusDecay::G4BetaPlusDecay(), G4DalitzDecayChannel::G4DalitzDecayChannel(), G4ECDecay::G4ECDecay(), G4ITDecay::G4ITDecay(), G4MuonDecayChannel::G4MuonDecayChannel(), G4MuonRadiativeDecayChannelWithSpin::G4MuonRadiativeDecayChannelWithSpin(), G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel(), G4NeutronDecay::G4NeutronDecay(), G4PionRadiativeDecayChannel::G4PionRadiativeDecayChannel(), G4ProtonDecay::G4ProtonDecay(), G4SFDecay::G4SFDecay(), G4TauLeptonicDecayChannel::G4TauLeptonicDecayChannel(), G4TritonDecay::G4TritonDecay(), G4RadioactiveDecay::LoadDecayTable(), G4TextPPRetriever::ModifyPropertyTable(), and G4DecayTableMessenger::SetNewValue().

SetDaughter() [1/2]

void G4VDecayChannel::SetDaughter	(	G4int	anIndex,
		const G4ParticleDefinition *	particle_type )

Definition at line 266 of file G4VDecayChannel.cc.

{
  if (parent_type != nullptr) SetDaughter(anIndex, parent_type->GetParticleName());
}

Referenced by G4TauMinus::Definition(), G4TauPlus::Definition(), G4AlphaDecay::G4AlphaDecay(), G4BetaMinusDecay::G4BetaMinusDecay(), G4BetaPlusDecay::G4BetaPlusDecay(), G4DalitzDecayChannel::G4DalitzDecayChannel(), G4ECDecay::G4ECDecay(), G4ITDecay::G4ITDecay(), G4MuonDecayChannel::G4MuonDecayChannel(), G4MuonRadiativeDecayChannelWithSpin::G4MuonRadiativeDecayChannelWithSpin(), G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel(), G4NeutronDecay::G4NeutronDecay(), G4PionRadiativeDecayChannel::G4PionRadiativeDecayChannel(), G4ProtonDecay::G4ProtonDecay(), G4SFDecay::G4SFDecay(), G4TauLeptonicDecayChannel::G4TauLeptonicDecayChannel(), G4TritonDecay::G4TritonDecay(), and SetDaughter().

SetDaughter() [2/2]

void G4VDecayChannel::SetDaughter	(	G4int	anIndex,
		const G4String &	particle_name )

Definition at line 210 of file G4VDecayChannel.cc.

{
  // check numberOfDaughters is positive
  if (numberOfDaughters <= 0) {
#ifdef G4VERBOSE
    if (verboseLevel > 0) {
      G4cout << "G4VDecayChannel::SetDaughter() - "
             << "Number of daughters is not defined" << G4endl;
    }
#endif
    return;
  }
 
  // An analysis of this code, shows that this method is called
  // only in the constructor of derived classes.
  // The general idea of this method is probably to support
  // the possibility to re-define daughters on the fly, however
  // this design is extremely problematic for MT mode, we thus
  // require (as practically happens) that the method is called only
  // at construction, i.e. when G4MT_daughters == 0
  // moreover this method can be called only after SetNumberOfDaugthers()
  // has been called (see previous if), in such a case daughters_name != nullptr
  //
  if (daughters_name == nullptr) {
    G4Exception("G4VDecayChannel::SetDaughter()", "PART112", FatalException,
                "Trying to add a daughter without specifying number of secondaries!");
    return;
  }
  if (G4MT_daughters != nullptr) {
    G4Exception("G4VDecayChannel::SetDaughter()", "PART111", FatalException,
                "Trying to modify a daughter of a decay channel, \
                 but decay channel already has daughters.");
    return;
  }
 
  // check an index
  if ((anIndex < 0) || (anIndex >= numberOfDaughters)) {
#ifdef G4VERBOSE
    if (verboseLevel > 0) {
      G4cout << "G4VDecayChannel::SetDaughter() - "
             << "index out of range " << anIndex << G4endl;
    }
#endif
  }
  else {
    // fill the name
    daughters_name[anIndex] = new G4String(particle_name);
#ifdef G4VERBOSE
    if (verboseLevel > 1) {
      G4cout << "G4VDecayChannel::SetDaughter[" << anIndex << "] :";
      G4cout << daughters_name[anIndex] << ":" << *daughters_name[anIndex] << G4endl;
    }
#endif
  }
}

SetNumberOfDaughters()

void G4VDecayChannel::SetNumberOfDaughters

(

G4int

value

)

Definition at line 196 of file G4VDecayChannel.cc.

{
  if (size > 0) {
    // remove old contents
    ClearDaughtersName();
    // cleate array
    daughters_name = new G4String*[size];
    for (G4int index = 0; index < size; ++index) {
      daughters_name[index] = nullptr;
    }
    numberOfDaughters = size;
  }
}

Referenced by G4TauMinus::Definition(), G4TauPlus::Definition(), G4AlphaDecay::G4AlphaDecay(), G4BetaMinusDecay::G4BetaMinusDecay(), G4BetaPlusDecay::G4BetaPlusDecay(), G4DalitzDecayChannel::G4DalitzDecayChannel(), G4ECDecay::G4ECDecay(), G4ITDecay::G4ITDecay(), G4MuonDecayChannel::G4MuonDecayChannel(), G4MuonRadiativeDecayChannelWithSpin::G4MuonRadiativeDecayChannelWithSpin(), G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel(), G4NeutronDecay::G4NeutronDecay(), G4PionRadiativeDecayChannel::G4PionRadiativeDecayChannel(), G4ProtonDecay::G4ProtonDecay(), G4SFDecay::G4SFDecay(), G4TauLeptonicDecayChannel::G4TauLeptonicDecayChannel(), and G4TritonDecay::G4TritonDecay().

SetParent() [1/2]

void G4VDecayChannel::SetParent

(

const G4ParticleDefinition *

particle_type

)

Definition at line 423 of file G4VDecayChannel.cc.

{
  if (parent_type != nullptr) SetParent(parent_type->GetParticleName());
}

Referenced by G4TauMinus::Definition(), G4TauPlus::Definition(), G4AlphaDecay::G4AlphaDecay(), G4BetaMinusDecay::G4BetaMinusDecay(), G4BetaPlusDecay::G4BetaPlusDecay(), G4DalitzDecayChannel::G4DalitzDecayChannel(), G4ECDecay::G4ECDecay(), G4ITDecay::G4ITDecay(), G4MuonDecayChannel::G4MuonDecayChannel(), G4MuonRadiativeDecayChannelWithSpin::G4MuonRadiativeDecayChannelWithSpin(), G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel(), G4NeutronDecay::G4NeutronDecay(), G4PionRadiativeDecayChannel::G4PionRadiativeDecayChannel(), G4ProtonDecay::G4ProtonDecay(), G4SFDecay::G4SFDecay(), G4TauLeptonicDecayChannel::G4TauLeptonicDecayChannel(), G4TritonDecay::G4TritonDecay(), and SetParent().

SetParent() [2/2]

void G4VDecayChannel::SetParent ( const G4String & particle_name )

inline

Definition at line 258 of file G4VDecayChannel.hh.

{
  delete parent_name;
  parent_name = new G4String(particle_name);
  G4MT_parent = nullptr;
}

SetPolarization()

void G4VDecayChannel::SetPolarization ( const G4ThreeVector & polar )

inline

Definition at line 300 of file G4VDecayChannel.hh.

{
  parent_polarization = polar;
}

Referenced by G4DecayWithSpin::AtRestDoIt(), and G4DecayWithSpin::PostStepDoIt().

SetRangeMass()

void G4VDecayChannel::SetRangeMass ( G4double val )

inline

Definition at line 295 of file G4VDecayChannel.hh.

{
  if (val >= 0.) rangeMass = val;
}

SetVerboseLevel()

void G4VDecayChannel::SetVerboseLevel ( G4int value )

inline

Definition at line 280 of file G4VDecayChannel.hh.

{
  verboseLevel = value;
}

Referenced by G4Decay::DecayIt().

Member Data Documentation

daughters_name

G4String** G4VDecayChannel::daughters_name = nullptr

protected

Definition at line 150 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), G4KL3DecayChannel::DecayIt(), DumpInfo(), G4VDecayChannel(), G4VDecayChannel(), GetDaughterName(), G4GeneralPhaseSpaceDecay::ManyBodyDecayIt(), G4DalitzDecayChannel::operator=(), G4KL3DecayChannel::operator=(), G4MuonDecayChannel::operator=(), G4MuonDecayChannelWithSpin::operator=(), G4MuonRadiativeDecayChannelWithSpin::operator=(), G4NeutronBetaDecayChannel::operator=(), G4PionRadiativeDecayChannel::operator=(), G4TauLeptonicDecayChannel::operator=(), operator=(), SetDaughter(), and SetNumberOfDaughters().

daughtersMutex

G4Mutex G4VDecayChannel::daughtersMutex

protected

Definition at line 168 of file G4VDecayChannel.hh.

Referenced by CheckAndFillDaughters(), ClearDaughtersName(), G4VDecayChannel(), operator=(), and ~G4VDecayChannel().

G4MT_daughters

G4ParticleDefinition** G4VDecayChannel::G4MT_daughters = nullptr

protected

Definition at line 164 of file G4VDecayChannel.hh.

Referenced by CheckAndFillDaughters(), ClearDaughtersName(), G4AlphaDecay::DecayIt(), G4DalitzDecayChannel::DecayIt(), G4ECDecay::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4NeutronBetaDecayChannel::DecayIt(), G4NeutronDecay::DecayIt(), G4PionRadiativeDecayChannel::DecayIt(), G4ProtonDecay::DecayIt(), G4TauLeptonicDecayChannel::DecayIt(), G4TritonDecay::DecayIt(), GetAngularMomentum(), GetDaughter(), G4GeneralPhaseSpaceDecay::ManyBodyDecayIt(), G4GeneralPhaseSpaceDecay::OneBodyDecayIt(), SetDaughter(), G4GeneralPhaseSpaceDecay::ThreeBodyDecayIt(), and G4GeneralPhaseSpaceDecay::TwoBodyDecayIt().

G4MT_daughters_mass

G4double* G4VDecayChannel::G4MT_daughters_mass = nullptr

protected

Definition at line 166 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), GetDaughterMass(), IsOKWithParentMass(), and ~G4VDecayChannel().

G4MT_daughters_width

G4double* G4VDecayChannel::G4MT_daughters_width = nullptr

protected

Definition at line 167 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), IsOKWithParentMass(), and ~G4VDecayChannel().

G4MT_parent

G4ParticleDefinition* G4VDecayChannel::G4MT_parent = nullptr

protected

Definition at line 163 of file G4VDecayChannel.hh.

Referenced by CheckAndFillParent(), G4AlphaDecay::DecayIt(), G4DalitzDecayChannel::DecayIt(), G4ECDecay::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4NeutronBetaDecayChannel::DecayIt(), G4NeutronDecay::DecayIt(), G4PionRadiativeDecayChannel::DecayIt(), G4ProtonDecay::DecayIt(), G4SFDecay::DecayIt(), G4TauLeptonicDecayChannel::DecayIt(), G4TritonDecay::DecayIt(), G4GeneralPhaseSpaceDecay::G4GeneralPhaseSpaceDecay(), GetAngularMomentum(), GetParent(), G4GeneralPhaseSpaceDecay::ManyBodyDecayIt(), G4GeneralPhaseSpaceDecay::OneBodyDecayIt(), SetParent(), G4GeneralPhaseSpaceDecay::ThreeBodyDecayIt(), and G4GeneralPhaseSpaceDecay::TwoBodyDecayIt().

G4MT_parent_mass

G4double G4VDecayChannel::G4MT_parent_mass = 0.0

protected

Definition at line 165 of file G4VDecayChannel.hh.

Referenced by G4PhaseSpaceDecayChannel::DecayIt(), and GetParentMass().

kinematics_name

G4String G4VDecayChannel::kinematics_name = ""

protected

Definition at line 144 of file G4VDecayChannel.hh.

Referenced by DumpInfo(), G4VDecayChannel(), GetKinematicsName(), G4DalitzDecayChannel::operator=(), G4KL3DecayChannel::operator=(), G4MuonDecayChannel::operator=(), G4MuonDecayChannelWithSpin::operator=(), G4MuonRadiativeDecayChannelWithSpin::operator=(), G4NeutronBetaDecayChannel::operator=(), G4PionRadiativeDecayChannel::operator=(), G4TauLeptonicDecayChannel::operator=(), and operator=().

noName

const G4String G4VDecayChannel::noName = " "

staticprotected

Definition at line 161 of file G4VDecayChannel.hh.

numberOfDaughters

G4int G4VDecayChannel::numberOfDaughters = 0

protected

Definition at line 172 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), G4GeneralPhaseSpaceDecay::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), DumpInfo(), G4VDecayChannel(), G4VDecayChannel(), GetAngularMomentum(), GetDaughter(), GetDaughterMass(), GetDaughterName(), GetNumberOfDaughters(), G4PhaseSpaceDecayChannel::IsOKWithParentMass(), IsOKWithParentMass(), G4GeneralPhaseSpaceDecay::ManyBodyDecayIt(), G4DalitzDecayChannel::operator=(), G4KL3DecayChannel::operator=(), G4MuonDecayChannel::operator=(), G4MuonDecayChannelWithSpin::operator=(), G4MuonRadiativeDecayChannelWithSpin::operator=(), G4NeutronBetaDecayChannel::operator=(), G4PionRadiativeDecayChannel::operator=(), G4TauLeptonicDecayChannel::operator=(), operator=(), SetDaughter(), G4PhaseSpaceDecayChannel::SetDaughterMasses(), and SetNumberOfDaughters().

parent_name

G4String* G4VDecayChannel::parent_name = nullptr

protected

Definition at line 148 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), G4GeneralPhaseSpaceDecay::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), G4VDecayChannel(), G4VDecayChannel(), GetParentName(), G4GeneralPhaseSpaceDecay::ManyBodyDecayIt(), G4DalitzDecayChannel::operator=(), G4KL3DecayChannel::operator=(), G4MuonDecayChannel::operator=(), G4MuonDecayChannelWithSpin::operator=(), G4MuonRadiativeDecayChannelWithSpin::operator=(), G4NeutronBetaDecayChannel::operator=(), G4PionRadiativeDecayChannel::operator=(), G4TauLeptonicDecayChannel::operator=(), operator=(), SetParent(), and ~G4VDecayChannel().

parent_polarization

G4ThreeVector G4VDecayChannel::parent_polarization

protected

Definition at line 156 of file G4VDecayChannel.hh.

Referenced by G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4VDecayChannel(), GetPolarization(), G4MuonRadiativeDecayChannelWithSpin::operator=(), operator=(), and SetPolarization().

parentMutex

G4Mutex G4VDecayChannel::parentMutex

protected

Definition at line 169 of file G4VDecayChannel.hh.

Referenced by CheckAndFillParent(), G4VDecayChannel(), operator=(), and ~G4VDecayChannel().

particletable

G4ParticleTable* G4VDecayChannel::particletable = nullptr

protected

Definition at line 159 of file G4VDecayChannel.hh.

Referenced by G4VDecayChannel(), G4VDecayChannel(), G4VDecayChannel(), G4VDecayChannel(), and operator=().

rangeMass

G4double G4VDecayChannel::rangeMass = 2.5

protected

Definition at line 153 of file G4VDecayChannel.hh.

Referenced by DynamicalMass(), G4VDecayChannel(), GetRangeMass(), IsOKWithParentMass(), operator=(), and SetRangeMass().

rbranch

G4double G4VDecayChannel::rbranch = 0.0

protected

Definition at line 146 of file G4VDecayChannel.hh.

Referenced by DumpInfo(), G4VDecayChannel(), G4VDecayChannel(), GetBR(), operator<(), G4DalitzDecayChannel::operator=(), G4KL3DecayChannel::operator=(), G4MuonDecayChannel::operator=(), G4MuonDecayChannelWithSpin::operator=(), G4MuonRadiativeDecayChannelWithSpin::operator=(), G4NeutronBetaDecayChannel::operator=(), G4PionRadiativeDecayChannel::operator=(), G4TauLeptonicDecayChannel::operator=(), operator=(), and SetBR().

verboseLevel

G4int G4VDecayChannel::verboseLevel = 1

protected

Definition at line 178 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), G4VDecayChannel(), GetAngularMomentum(), GetDaughter(), GetDaughterMass(), GetDaughterName(), GetVerboseLevel(), G4DalitzDecayChannel::operator=(), G4KL3DecayChannel::operator=(), G4MuonDecayChannel::operator=(), G4MuonDecayChannelWithSpin::operator=(), G4MuonRadiativeDecayChannelWithSpin::operator=(), G4NeutronBetaDecayChannel::operator=(), G4PionRadiativeDecayChannel::operator=(), G4TauLeptonicDecayChannel::operator=(), operator=(), SetDaughter(), and SetVerboseLevel().

---

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

• G4VDecayChannel.hh
• G4VDecayChannel.cc