#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
void	ClearDaughtersName ()

void	CheckAndFillDaughters ()

void	CheckAndFillParent ()

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

	G4VDecayChannel ()

	G4VDecayChannel (const G4VDecayChannel &)

G4VDecayChannel &	operator= (const G4VDecayChannel &)

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 51 of file G4VDecayChannel.cc.

  : kinematics_name(aName), parent_polarization(), 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 59 of file G4VDecayChannel.cc.

  : kinematics_name(aName), rbranch(theBR), parent_polarization(),
    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 169 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.

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

◆ G4VDecayChannel() [4/4]

G4VDecayChannel::G4VDecayChannel ( const G4VDecayChannel & right )

protected

Definition at line 96 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 337 of file G4VDecayChannel.hh.

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

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

◆ CheckAndFillParent()

void G4VDecayChannel::CheckAndFillParent ( )

inlineprotected

Definition at line 348 of file G4VDecayChannel.hh.

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

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

◆ ClearDaughtersName()

void G4VDecayChannel::ClearDaughtersName ( )

protected

Definition at line 183 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

◆ DumpInfo()

void G4VDecayChannel::DumpInfo ( )

Definition at line 560 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 585 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 492 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

◆ GetDaughter()

G4ParticleDefinition * G4VDecayChannel::GetDaughter ( G4int anIndex )

inline

Definition at line 201 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];
  }
  else
  {
    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 239 of file G4VDecayChannel.hh.

{
  if ( (anIndex>=0) && (anIndex<numberOfDaughters) )
  {
    return G4MT_daughters_mass[anIndex];
  }
  else
  {
    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 221 of file G4VDecayChannel.hh.

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

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

◆ GetKinematicsName()

const G4String & G4VDecayChannel::GetKinematicsName ( ) const

inline

Definition at line 292 of file G4VDecayChannel.hh.

{
  return kinematics_name;
}

Referenced by G4TextPPReporter::GeneratePropertyTable().

◆ GetNumberOfDaughters()

G4int G4VDecayChannel::GetNumberOfDaughters ( ) const

inline

Definition at line 286 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 257 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 272 of file G4VDecayChannel.hh.

{
  return G4MT_parent_mass;
}

◆ GetParentName()

const G4String & G4VDecayChannel::GetParentName ( ) const

inline

Definition at line 266 of file G4VDecayChannel.hh.

{
  return *parent_name;
}

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

◆ GetPolarization()

const G4ThreeVector & G4VDecayChannel::GetPolarization ( ) const

inline

Definition at line 331 of file G4VDecayChannel.hh.

{
  return parent_polarization;
}

◆ GetRangeMass()

G4double G4VDecayChannel::GetRangeMass ( ) const

inline

Definition at line 316 of file G4VDecayChannel.hh.

{
  return rangeMass;
}

◆ GetVerboseLevel()

G4int G4VDecayChannel::GetVerboseLevel ( ) const

inline

Definition at line 310 of file G4VDecayChannel.hh.

{
  return verboseLevel;
}

Referenced by G4KL3DecayChannel::DalitzDensity(), G4Decay::DecayIt(), G4GeneralPhaseSpaceDecay::DecayIt(), G4DalitzDecayChannel::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4NeutronBetaDecayChannel::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), G4PionRadiativeDecayChannel::DecayIt(), G4TauLeptonicDecayChannel::DecayIt(), 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 G4PhaseSpaceDecayChannel.

Definition at line 607 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 69 of file G4VDecayChannel.hh.

69{ return (this != &r); }

◆ operator<()

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

inline

Definition at line 195 of file G4VDecayChannel.hh.

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

◆ operator=()

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

protected

Definition at line 130 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 68 of file G4VDecayChannel.hh.

68{ return (this == &r); }

◆ SetBR()

void G4VDecayChannel::SetBR ( G4double value )

Definition at line 624 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 300 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 234 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 217 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 486 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 278 of file G4VDecayChannel.hh.

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

◆ SetPolarization()

void G4VDecayChannel::SetPolarization ( const G4ThreeVector & polar )

inline

Definition at line 325 of file G4VDecayChannel.hh.

{
  parent_polarization = polar;
}

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

◆ SetRangeMass()

void G4VDecayChannel::SetRangeMass ( G4double val )

inline

Definition at line 322 of file G4VDecayChannel.hh.

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

◆ SetVerboseLevel()

void G4VDecayChannel::SetVerboseLevel ( G4int value )

inline

Definition at line 304 of file G4VDecayChannel.hh.

{
  verboseLevel = value;
}

Referenced by G4Decay::DecayIt().

Member Data Documentation

◆ daughters_name

G4String** G4VDecayChannel::daughters_name = nullptr

protected

Definition at line 158 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), G4KL3DecayChannel::DecayIt(), DumpInfo(), 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 177 of file G4VDecayChannel.hh.

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

◆ G4MT_daughters

G4ParticleDefinition** G4VDecayChannel::G4MT_daughters = nullptr

protected

◆ G4MT_daughters_mass

G4double* G4VDecayChannel::G4MT_daughters_mass = nullptr

protected

Definition at line 175 of file G4VDecayChannel.hh.

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

◆ G4MT_daughters_width

G4double* G4VDecayChannel::G4MT_daughters_width = nullptr

protected

Definition at line 176 of file G4VDecayChannel.hh.

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

◆ G4MT_parent

G4ParticleDefinition* G4VDecayChannel::G4MT_parent = nullptr

protected

Definition at line 172 of file G4VDecayChannel.hh.

Referenced by CheckAndFillParent(), G4DalitzDecayChannel::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), G4MuonDecayChannelWithSpin::DecayIt(), G4MuonRadiativeDecayChannelWithSpin::DecayIt(), G4NeutronBetaDecayChannel::DecayIt(), G4PionRadiativeDecayChannel::DecayIt(), G4TauLeptonicDecayChannel::DecayIt(), G4AlphaDecay::DecayIt(), G4BetaMinusDecay::DecayIt(), G4BetaPlusDecay::DecayIt(), G4ECDecay::DecayIt(), G4ITDecay::DecayIt(), G4NeutronDecay::DecayIt(), G4ProtonDecay::DecayIt(), G4SFDecay::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 174 of file G4VDecayChannel.hh.

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

◆ kinematics_name

G4String G4VDecayChannel::kinematics_name = ""

protected

Definition at line 152 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 170 of file G4VDecayChannel.hh.

◆ numberOfDaughters

G4int G4VDecayChannel::numberOfDaughters = 0

protected

◆ parent_name

G4String* G4VDecayChannel::parent_name = nullptr

protected

Definition at line 156 of file G4VDecayChannel.hh.

Referenced by ClearDaughtersName(), G4GeneralPhaseSpaceDecay::DecayIt(), G4PhaseSpaceDecayChannel::DecayIt(), 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 164 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 178 of file G4VDecayChannel.hh.

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

◆ particletable

G4ParticleTable* G4VDecayChannel::particletable = nullptr

protected

Definition at line 167 of file G4VDecayChannel.hh.

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

◆ rangeMass

G4double G4VDecayChannel::rangeMass = 2.5

protected

Definition at line 161 of file G4VDecayChannel.hh.

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

◆ rbranch

G4double G4VDecayChannel::rbranch = 0.0

protected

Definition at line 154 of file G4VDecayChannel.hh.

Referenced by DumpInfo(), 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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4VDecayChannel() [1/4]

◆ G4VDecayChannel() [2/4]

◆ ~G4VDecayChannel()

◆ G4VDecayChannel() [3/4]

◆ G4VDecayChannel() [4/4]

Member Function Documentation

◆ CheckAndFillDaughters()

◆ CheckAndFillParent()

◆ ClearDaughtersName()

◆ DecayIt()

◆ DumpInfo()

◆ DynamicalMass()

◆ GetAngularMomentum()

◆ GetBR()

◆ GetDaughter()

◆ GetDaughterMass()

◆ GetDaughterName()

◆ GetKinematicsName()

◆ GetNumberOfDaughters()

◆ GetParent()

◆ GetParentMass()

◆ GetParentName()

◆ GetPolarization()

◆ GetRangeMass()

◆ GetVerboseLevel()

◆ IsOKWithParentMass()

◆ operator!=()

◆ operator<()

◆ operator=()

◆ operator==()

◆ SetBR()

◆ SetDaughter() [1/2]

◆ SetDaughter() [2/2]

◆ SetNumberOfDaughters()

◆ SetParent() [1/2]

◆ SetParent() [2/2]

◆ SetPolarization()

◆ SetRangeMass()

◆ SetVerboseLevel()

Member Data Documentation

◆ daughters_name

◆ daughtersMutex

◆ G4MT_daughters

◆ G4MT_daughters_mass

◆ G4MT_daughters_width

◆ G4MT_parent

◆ G4MT_parent_mass

◆ kinematics_name

◆ noName

◆ numberOfDaughters

◆ parent_name

◆ parent_polarization

◆ parentMutex

◆ particletable

◆ rangeMass

◆ rbranch

◆ verboseLevel