G4PhaseSpaceDecayChannel Class Reference

#include <G4PhaseSpaceDecayChannel.hh>

Inheritance diagram for G4PhaseSpaceDecayChannel:

Public Types
enum	{ MAX_N_DAUGHTERS = 5 }

Public Member Functions
	G4PhaseSpaceDecayChannel (G4int Verbose=1)
	G4PhaseSpaceDecayChannel (const G4String &theParentName, G4double theBR, G4int theNumberOfDaughters, const G4String &theDaughterName1, const G4String &theDaughterName2="", const G4String &theDaughterName3="", const G4String &theDaughterName4="", const G4String &theDaughterName5="")
	~G4PhaseSpaceDecayChannel () override=default
G4bool	SetDaughterMasses (G4double masses[])
G4bool	SampleDaughterMasses ()
G4DecayProducts *	DecayIt (G4double) override
G4bool	IsOKWithParentMass (G4double parentMass) override
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

Static Public Member Functions
static G4double	Pmx (G4double e, G4double p1, G4double p2)

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 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 42 of file G4PhaseSpaceDecayChannel.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
MAX_N_DAUGHTERS

Definition at line 45 of file G4PhaseSpaceDecayChannel.hh.

    {
      MAX_N_DAUGHTERS = 5
    };

Constructor & Destructor Documentation

G4PhaseSpaceDecayChannel() [1/2]

G4PhaseSpaceDecayChannel::G4PhaseSpaceDecayChannel

(

G4int

Verbose = 1

)

Definition at line 42 of file G4PhaseSpaceDecayChannel.cc.

  : G4VDecayChannel("Phase Space", Verbose)
{}

G4PhaseSpaceDecayChannel() [2/2]

G4PhaseSpaceDecayChannel::G4PhaseSpaceDecayChannel	(	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 46 of file G4PhaseSpaceDecayChannel.cc.

  : G4VDecayChannel("Phase Space", theParentName, theBR, theNumberOfDaughters, theDaughterName1,
                    theDaughterName2, theDaughterName3, theDaughterName4, theDaughterName5)
{}

~G4PhaseSpaceDecayChannel()

G4PhaseSpaceDecayChannel::~G4PhaseSpaceDecayChannel ( )

overridedefault

Member Function Documentation

DecayIt()

G4DecayProducts * G4PhaseSpaceDecayChannel::DecayIt ( G4double parentMass )

overridevirtual

Implements G4VDecayChannel.

Definition at line 57 of file G4PhaseSpaceDecayChannel.cc.

{
#ifdef G4VERBOSE
  if (GetVerboseLevel() > 1) G4cout << "G4PhaseSpaceDecayChannel::DecayIt()" << G4endl;
#endif
 
  G4DecayProducts* products = nullptr;
 
  CheckAndFillParent();
  CheckAndFillDaughters();
 
  if (parentMass > 0.0)
    current_parent_mass.Put(parentMass);
  else
    current_parent_mass.Put(G4MT_parent_mass);
 
  switch (numberOfDaughters) {
    case 0:
#ifdef G4VERBOSE
      if (GetVerboseLevel() > 0) {
        G4cout << "G4PhaseSpaceDecayChannel::DecayIt() -";
        G4cout << " daughters not defined " << G4endl;
      }
#endif
      break;
    case 1:
      products = OneBodyDecayIt();
      break;
    case 2:
      products = TwoBodyDecayIt();
      break;
    case 3:
      products = ThreeBodyDecayIt();
      break;
    default:
      products = ManyBodyDecayIt();
      break;
  }
#ifdef G4VERBOSE
  if ((products == nullptr) && (GetVerboseLevel() > 0)) {
    G4cout << "G4PhaseSpaceDecayChannel::DecayIt() - ";
    G4cout << *parent_name << " cannot decay " << G4endl;
    DumpInfo();
  }
#endif
  return products;
}

IsOKWithParentMass()

G4bool G4PhaseSpaceDecayChannel::IsOKWithParentMass ( G4double parentMass )

overridevirtual

Reimplemented from G4VDecayChannel.

Definition at line 736 of file G4PhaseSpaceDecayChannel.cc.

{
  if (!useGivenDaughterMass) return G4VDecayChannel::IsOKWithParentMass(parentMass);
 
  CheckAndFillParent();
  CheckAndFillDaughters();
 
  G4double sumOfDaughterMassMin = 0.0;
  for (G4int index = 0; index < numberOfDaughters; ++index) {
    sumOfDaughterMassMin += givenDaughterMasses[index];
  }
  return (parentMass >= sumOfDaughterMassMin);
}

Pmx()

G4double G4PhaseSpaceDecayChannel::Pmx ( G4double e, G4double p1, G4double p2 )

static

Definition at line 750 of file G4PhaseSpaceDecayChannel.cc.

{
  // calcurate momentum of daughter particles in two-body decay
  G4double ppp = (e + p1 + p2) * (e + p1 - p2) * (e - p1 + p2) * (e - p1 - p2) / (4.0 * e * e);
  if (ppp > 0) return std::sqrt(ppp);
  return -1.;
}

Referenced by G4DalitzDecayChannel::DecayIt().

SampleDaughterMasses()

G4bool G4PhaseSpaceDecayChannel::SampleDaughterMasses

(

)

Definition at line 730 of file G4PhaseSpaceDecayChannel.cc.

{
  useGivenDaughterMass = false;
  return useGivenDaughterMass;
}

SetDaughterMasses()

G4bool G4PhaseSpaceDecayChannel::SetDaughterMasses

(

G4double

masses[]

)

Definition at line 721 of file G4PhaseSpaceDecayChannel.cc.

{
  for (G4int idx = 0; idx < numberOfDaughters; ++idx) {
    givenDaughterMasses[idx] = masses[idx];
  }
  useGivenDaughterMass = true;
  return useGivenDaughterMass;
}

---

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

• G4PhaseSpaceDecayChannel.hh
• G4PhaseSpaceDecayChannel.cc