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="")
virtual	~G4PhaseSpaceDecayChannel ()
G4bool	SetDaughterMasses (G4double masses[])
G4bool	SampleDaughterMasses ()
virtual G4DecayProducts *	DecayIt (G4double)
G4bool	IsOKWithParentMass (G4double parentMass)
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
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

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

Additional Inherited Members
Protected Member Functions inherited from G4VDecayChannel
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 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 46 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 47 of file G4PhaseSpaceDecayChannel.cc.

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

~G4PhaseSpaceDecayChannel()

G4PhaseSpaceDecayChannel::~G4PhaseSpaceDecayChannel ( )

virtual

Definition at line 63 of file G4PhaseSpaceDecayChannel.cc.

{
}

Member Function Documentation

DecayIt()

G4DecayProducts * G4PhaseSpaceDecayChannel::DecayIt ( G4double  parentMass )

virtual

Implements G4VDecayChannel.

Definition at line 68 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 )

virtual

Reimplemented from G4VDecayChannel.

Definition at line 853 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 870 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);
   else       return -1.;
}

Referenced by G4DalitzDecayChannel::DecayIt().

SampleDaughterMasses()

G4bool G4PhaseSpaceDecayChannel::SampleDaughterMasses

(

)

Definition at line 846 of file G4PhaseSpaceDecayChannel.cc.

{
  useGivenDaughterMass = false;
  return useGivenDaughterMass; 
}

SetDaughterMasses()

G4bool G4PhaseSpaceDecayChannel::SetDaughterMasses

(

G4double

masses[]

)

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