#include <G4MuonDecayChannelWithSpin.hh>

Inheritance diagram for G4MuonDecayChannelWithSpin:

Public Member Functions
virtual G4DecayProducts *	DecayIt (G4double)

void	DumpInfo ()

	G4MuonDecayChannelWithSpin (const G4String &theParentName, G4double theBR)

G4int	GetAngularMomentum ()

G4double	GetBR () const

G4ParticleDefinition *	GetDaughter (G4int anIndex)

G4double	GetDaughterMass (G4int anIndex) const

const G4String &	GetDaughterName (G4int anIndex) const

const G4String &	GetKinematicsName () const

G4int	GetNumberOfDaughters () const

G4ParticleDefinition *	GetParent ()

G4double	GetParentMass () const

const G4String &	GetParentName () const

const G4ThreeVector &	GetPolarization () const

G4double	GetRangeMass () const

G4int	GetVerboseLevel () const

virtual G4bool	IsOKWithParentMass (G4double parentMass)

G4bool	operator!= (const G4VDecayChannel &r) const

G4bool	operator< (const G4VDecayChannel &right) const

G4bool	operator== (const G4VDecayChannel &r) const

void	SetBR (G4double value)

void	SetDaughter (G4int anIndex, const G4ParticleDefinition *particle_type)

void	SetDaughter (G4int anIndex, const G4String &particle_name)

void	SetNumberOfDaughters (G4int value)

void	SetParent (const G4ParticleDefinition *particle_type)

void	SetParent (const G4String &particle_name)

void	SetPolarization (const G4ThreeVector &)

void	SetRangeMass (G4double val)

void	SetVerboseLevel (G4int value)

virtual	~G4MuonDecayChannelWithSpin ()

Protected Member Functions
void	CheckAndFillDaughters ()

void	CheckAndFillParent ()

void	ClearDaughtersName ()

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

	G4MuonDecayChannelWithSpin (const G4MuonDecayChannelWithSpin &)

G4MuonDecayChannelWithSpin &	operator= (const G4MuonDecayChannelWithSpin &)

Protected Attributes
G4String **	daughters_name = nullptr

G4Mutex	daughtersMutex

G4ParticleDefinition **	G4MT_daughters = nullptr

G4double *	G4MT_daughters_mass = nullptr

G4double *	G4MT_daughters_width = nullptr

G4ParticleDefinition *	G4MT_parent = nullptr

G4double	G4MT_parent_mass = 0.0

G4String	kinematics_name = ""

G4int	numberOfDaughters = 0

G4String *	parent_name = nullptr

G4ThreeVector	parent_polarization

G4Mutex	parentMutex

G4ParticleTable *	particletable = nullptr

G4double	rangeMass = 2.5

G4double	rbranch = 0.0

G4int	verboseLevel = 1

Static Protected Attributes
static const G4String	noName = " "

Private Member Functions
G4double	F_c (G4double x, G4double x0, G4double omega)

G4double	F_theta (G4double x, G4double x0, G4double omega)

void	FillDaughters ()

void	FillParent ()

	G4MuonDecayChannelWithSpin ()

const G4String &	GetNoName () const

G4double	R_c (G4double x, G4double omega)

Detailed Description

Definition at line 51 of file G4MuonDecayChannelWithSpin.hh.

Constructor & Destructor Documentation

◆ G4MuonDecayChannelWithSpin() [1/3]

G4MuonDecayChannelWithSpin::G4MuonDecayChannelWithSpin	(	const G4String &	theParentName,
		G4double	theBR
	)

Definition at line 51 of file G4MuonDecayChannelWithSpin.cc.

  : G4MuonDecayChannel(theParentName,theBR)
{
}

◆ ~G4MuonDecayChannelWithSpin()

G4MuonDecayChannelWithSpin::~G4MuonDecayChannelWithSpin ( )

virtual

Definition at line 58 of file G4MuonDecayChannelWithSpin.cc.

59{

60}

◆ G4MuonDecayChannelWithSpin() [2/3]

G4MuonDecayChannelWithSpin::G4MuonDecayChannelWithSpin ( const G4MuonDecayChannelWithSpin & right )

protected

Definition at line 62 of file G4MuonDecayChannelWithSpin.cc.

  : G4MuonDecayChannel(right)
{
}

◆ G4MuonDecayChannelWithSpin() [3/3]

G4MuonDecayChannelWithSpin::G4MuonDecayChannelWithSpin ( )

private

Definition at line 46 of file G4MuonDecayChannelWithSpin.cc.

  : G4MuonDecayChannel()
{
}

Member Function Documentation

◆ CheckAndFillDaughters()

void G4VDecayChannel::CheckAndFillDaughters ( )

inlineprotectedinherited

Definition at line 337 of file G4VDecayChannel.hh.

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

References G4VDecayChannel::daughtersMutex, G4VDecayChannel::FillDaughters(), G4VDecayChannel::G4MT_daughters, and G4TemplateAutoLock< _Mutex_t >::unlock().

Referenced by G4GeneralPhaseSpaceDecay::DecayIt(), G4DalitzDecayChannel::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), 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(), G4VDecayChannel::GetAngularMomentum(), G4VDecayChannel::GetDaughter(), G4PhaseSpaceDecayChannel::IsOKWithParentMass(), and G4VDecayChannel::IsOKWithParentMass().

◆ CheckAndFillParent()

void G4VDecayChannel::CheckAndFillParent ( )

inlineprotectedinherited

Definition at line 348 of file G4VDecayChannel.hh.

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

References G4VDecayChannel::FillParent(), G4VDecayChannel::G4MT_parent, G4VDecayChannel::parentMutex, and G4TemplateAutoLock< _Mutex_t >::unlock().

Referenced by G4GeneralPhaseSpaceDecay::DecayIt(), G4DalitzDecayChannel::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), 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(), G4VDecayChannel::FillDaughters(), G4VDecayChannel::GetParent(), G4PhaseSpaceDecayChannel::IsOKWithParentMass(), and G4VDecayChannel::IsOKWithParentMass().

◆ ClearDaughtersName()

void G4VDecayChannel::ClearDaughtersName ( )

protectedinherited

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;
}

References G4VDecayChannel::daughters_name, G4VDecayChannel::daughtersMutex, G4cout, G4endl, G4VDecayChannel::G4MT_daughters, G4VDecayChannel::G4MT_daughters_mass, G4VDecayChannel::G4MT_daughters_width, G4VDecayChannel::numberOfDaughters, G4VDecayChannel::parent_name, and G4VDecayChannel::verboseLevel.

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

◆ DecayIt()

G4DecayProducts * G4MuonDecayChannelWithSpin::DecayIt ( G4double )

virtual

Reimplemented from G4MuonDecayChannel.

Definition at line 99 of file G4MuonDecayChannelWithSpin.cc.

{
  // This version assumes V-A coupling with 1st order radiative correctons,
  //              the standard model Michel parameter values, but
  //              gives incorrect energy spectrum for neutrinos
 
#ifdef G4VERBOSE
  if (GetVerboseLevel()>1) G4cout << "G4MuonDecayChannelWithSpin::DecayIt ";
#endif
 
  CheckAndFillParent();
  CheckAndFillDaughters();
 
  // parent mass
  G4double parentmass = G4MT_parent->GetPDGMass();
 
  G4double EMMU = parentmass;
 
  //daughters'mass
  G4double daughtermass[3]; 
  //G4double sumofdaughtermass = 0.0;
  for (G4int index=0; index<3; ++index)
  {
    daughtermass[index] = G4MT_daughters[index]->GetPDGMass();
    //sumofdaughtermass += daughtermass[index];
  }
 
  G4double EMASS = daughtermass[0];
 
  // create parent G4DynamicParticle at rest
  G4ThreeVector dummy;
  G4DynamicParticle* parentparticle
    = new G4DynamicParticle( G4MT_parent, dummy, 0.0);
  // create G4Decayproducts
  G4DecayProducts *products = new G4DecayProducts(*parentparticle);
  delete parentparticle;
 
  // calculate electron energy
 
  G4double michel_rho   = 0.75; //Standard Model Michel rho
  G4double michel_delta = 0.75; //Standard Model Michel delta
  G4double michel_xsi   = 1.00; //Standard Model Michel xsi
  G4double michel_eta   = 0.00; //Standard Model eta
 
  G4double rndm, x, ctheta;
 
  G4double FG; 
  G4double FG_max = 2.00;
 
  G4double W_mue  = (EMMU*EMMU+EMASS*EMASS)/(2.*EMMU);
  G4double x0     =           EMASS/W_mue;
 
  G4double x0_squared = x0*x0;
  
  // ***************************************************
  //     x0 <= x <= 1.   and   -1 <= y <= 1
  //
  //     F(x,y) = f(x)*g(x,y);   g(x,y) = 1.+g(x)*y
  // ***************************************************
 
  // ***** sampling F(x,y) directly (brute force) *****
 
  const std::size_t MAX_LOOP=10000;
  for (std::size_t loop_count=0; loop_count<MAX_LOOP; ++loop_count)
  {
    // Sample the positron energy by sampling from F
 
    rndm = G4UniformRand();
 
    x = x0 + rndm*(1.-x0);
 
    G4double x_squared = x*x;
 
    G4double F_IS, F_AS, G_IS, G_AS;
 
    F_IS = 1./6.*(-2.*x_squared+3.*x-x0_squared);
    F_AS = 1./6.*std::sqrt(x_squared-x0_squared)
                *(2.*x-2.+std::sqrt(1.-x0_squared));
 
    G_IS = 2./9.*(michel_rho-0.75)*(4.*x_squared-3.*x-x0_squared);
    G_IS = G_IS + michel_eta*(1.-x)*x0;
 
    G_AS = 3.*(michel_xsi-1.)*(1.-x);
    G_AS = G_AS+2.*(michel_xsi*michel_delta-0.75)
                  *(4.*x-4.+std::sqrt(1.-x0_squared));
    G_AS = 1./9.*std::sqrt(x_squared-x0_squared)*G_AS;
 
    F_IS = F_IS + G_IS;
    F_AS = F_AS + G_AS;
 
    // *** Radiative Corrections ***
 
    const G4double omega =  std::log(EMMU/EMASS);
    G4double R_IS = F_c(x,x0,omega);
 
    G4double F = 6.*F_IS + R_IS/std::sqrt(x_squared-x0_squared);
 
    // *** Radiative Corrections ***
 
    G4double R_AS = F_theta(x,x0,omega);
 
    rndm = G4UniformRand();
 
    ctheta = 2.*rndm-1.;
 
    G4double G = 6.*F_AS - R_AS/std::sqrt(x_squared-x0_squared);
 
    FG = std::sqrt(x_squared-x0_squared)*F*(1.+(G/F)*ctheta);
 
    if(FG>FG_max)
    {
      G4Exception("G4MuonDecayChannelWithSpin::DecayIt()",
                  "PART113", JustWarning,
                  "Problem in Muon Decay: FG > FG_max");
      FG_max = FG;
    }
 
    rndm = G4UniformRand();
 
    if (FG >= rndm*FG_max) break;
  }
 
  G4double energy = x * W_mue;
 
  rndm = G4UniformRand();
 
  G4double phi = twopi * rndm;
 
  if(energy < EMASS) energy = EMASS;
 
  // Calculate daughter momentum
  G4double daughtermomentum[3];
 
  daughtermomentum[0] = std::sqrt(energy*energy - EMASS*EMASS);
 
  G4double stheta = std::sqrt(1.-ctheta*ctheta);
  G4double cphi = std::cos(phi);
  G4double sphi = std::sin(phi);
 
  // Coordinates of the decay positron with respect to the muon spin
  G4double px = stheta*cphi;
  G4double py = stheta*sphi;
  G4double pz = ctheta;
 
  G4ThreeVector direction0(px,py,pz);
 
  direction0.rotateUz(parent_polarization);
 
  G4DynamicParticle * daughterparticle0 
    = new G4DynamicParticle( G4MT_daughters[0], daughtermomentum[0]*direction0);
 
  products->PushProducts(daughterparticle0);
 
 
  // daughter 1 ,2 (neutrinos)
  // create neutrinos in the C.M frame of two neutrinos
  G4double energy2 = parentmass-energy; 
  G4double vmass = std::sqrt((energy2-daughtermomentum[0])
                 * (energy2+daughtermomentum[0]));
  G4double beta = -1.0*daughtermomentum[0]/energy2;
  G4double costhetan = 2.*G4UniformRand()-1.0;
  G4double sinthetan = std::sqrt((1.0-costhetan)*(1.0+costhetan));
  G4double phin  = twopi*G4UniformRand()*rad;
  G4double sinphin = std::sin(phin);
  G4double cosphin = std::cos(phin);
 
  G4ThreeVector direction1(sinthetan*cosphin,sinthetan*sinphin,costhetan);
  G4DynamicParticle * daughterparticle1 
    = new G4DynamicParticle( G4MT_daughters[1], direction1*(vmass/2.));
  G4DynamicParticle * daughterparticle2
    = new G4DynamicParticle( G4MT_daughters[2], direction1*(-1.0*vmass/2.));
 
  // boost to the muon rest frame
  G4LorentzVector p4;
  p4 = daughterparticle1->Get4Momentum();
  p4.boost( direction0.x()*beta, direction0.y()*beta, direction0.z()*beta);
  daughterparticle1->Set4Momentum(p4);
  p4 = daughterparticle2->Get4Momentum();
  p4.boost( direction0.x()*beta, direction0.y()*beta, direction0.z()*beta);
  daughterparticle2->Set4Momentum(p4);
  products->PushProducts(daughterparticle1);
  products->PushProducts(daughterparticle2);
  daughtermomentum[1] = daughterparticle1->GetTotalMomentum();
  daughtermomentum[2] = daughterparticle2->GetTotalMomentum();
 
  // output message
#ifdef G4VERBOSE
  if (GetVerboseLevel()>1)
  {
    G4cout << "G4MuonDecayChannelWithSpin::DecayIt ";
    G4cout << "  create decay products in rest frame " <<G4endl;
    G4double TT = daughterparticle0->GetTotalEnergy()
      + daughterparticle1->GetTotalEnergy()
      + daughterparticle2->GetTotalEnergy(); 
    G4cout << "e  " << daughterparticle0->GetTotalEnergy()/MeV << G4endl; 
    G4cout << "nu1" << daughterparticle1->GetTotalEnergy()/MeV << G4endl; 
    G4cout << "nu2" << daughterparticle2->GetTotalEnergy()/MeV << G4endl; 
    G4cout << "total" << (TT-parentmass)/keV << G4endl;
    if (GetVerboseLevel()>2) { products->DumpInfo(); }
  }
#endif
 
  return products;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CLHEP::HepLorentzVector::boost(), G4VDecayChannel::CheckAndFillDaughters(), G4VDecayChannel::CheckAndFillParent(), G4DecayProducts::DumpInfo(), G4INCL::KinematicsUtils::energy(), F_c(), F_theta(), G4cout, G4endl, G4Exception(), G4VDecayChannel::G4MT_daughters, G4VDecayChannel::G4MT_parent, G4UniformRand, G4DynamicParticle::Get4Momentum(), G4ParticleDefinition::GetPDGMass(), G4DynamicParticle::GetTotalEnergy(), G4DynamicParticle::GetTotalMomentum(), G4VDecayChannel::GetVerboseLevel(), JustWarning, keV, MeV, G4VDecayChannel::parent_polarization, G4DecayProducts::PushProducts(), rad, CLHEP::Hep3Vector::rotateUz(), G4DynamicParticle::Set4Momentum(), twopi, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

◆ DumpInfo()

void G4VDecayChannel::DumpInfo ( )

inherited

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;
}

References G4VDecayChannel::daughters_name, G4cout, G4endl, G4VDecayChannel::kinematics_name, G4VDecayChannel::numberOfDaughters, and G4VDecayChannel::rbranch.

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

◆ DynamicalMass()

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

protectedinherited

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;
}

References G4UniformRand, and G4VDecayChannel::rangeMass.

Referenced by G4PhaseSpaceDecayChannel::ThreeBodyDecayIt(), and G4PhaseSpaceDecayChannel::TwoBodyDecayIt().

◆ F_c()

G4double G4MuonDecayChannelWithSpin::F_c	(	G4double	x,
		G4double	x0,
		G4double	omega
	)

inlineprivate

Definition at line 83 of file G4MuonDecayChannelWithSpin.hh.

{
  G4double f_c;
 
  f_c = (5.+17.*x-34.*x*x)*(omega+std::log(x))-22.*x+34.*x*x;
  f_c = (1.-x)/(3.*x*x)*f_c;
  f_c = (6.-4.*x)*R_c(x,omega)+(6.-6.*x)*std::log(x) + f_c;
  f_c = (CLHEP::fine_structure_const/CLHEP::twopi) * (x*x-x0*x0) * f_c;
 
  return f_c;
}

References CLHEP::fine_structure_const, R_c(), and CLHEP::twopi.

Referenced by DecayIt().

◆ F_theta()

G4double G4MuonDecayChannelWithSpin::F_theta	(	G4double	x,
		G4double	x0,
		G4double	omega
	)

inlineprivate

Definition at line 96 of file G4MuonDecayChannelWithSpin.hh.

{
  G4double f_theta;
 
  f_theta = (1.+x+34*x*x)*(omega+std::log(x))+3.-7.*x-32.*x*x;
  f_theta = f_theta + ((4.*(1.-x)*(1.-x))/x)*std::log(1.-x);
  f_theta = (1.-x)/(3.*x*x) * f_theta;
  f_theta = (2.-4.*x)*R_c(x,omega)+(2.-6.*x)*std::log(x)-f_theta;
  f_theta = (CLHEP::fine_structure_const/CLHEP::twopi) * (x*x-x0*x0) * f_theta;
 
  return f_theta;
}

References CLHEP::fine_structure_const, R_c(), and CLHEP::twopi.

Referenced by DecayIt().

◆ FillDaughters()

void G4VDecayChannel::FillDaughters ( )

privateinherited

Definition at line 308 of file G4VDecayChannel.cc.

{
  G4AutoLock lock(&daughtersMutex);
 
  // Double check, check again if another thread has already filled this, in
  // case do not need to do anything
  if ( G4MT_daughters != nullptr ) return;
 
  G4int index;
  
#ifdef G4VERBOSE
  if (verboseLevel>1) G4cout << "G4VDecayChannel::FillDaughters()" << G4endl;
#endif
  if (G4MT_daughters != nullptr)
  {
    delete [] G4MT_daughters;
    G4MT_daughters = nullptr;
  }
 
  // parent mass
  CheckAndFillParent();
  G4double parentmass = G4MT_parent->GetPDGMass();
 
  //
  G4double sumofdaughtermass = 0.0;
  G4double sumofdaughterwidthsq = 0.0;
 
  if ((numberOfDaughters <=0) || (daughters_name == nullptr) )
  {
#ifdef G4VERBOSE
    if (verboseLevel>0)
    {
      G4cout << "G4VDecayChannel::FillDaughters() - "
             << "[ " << G4MT_parent->GetParticleName() << " ]"
             << "numberOfDaughters is not defined yet";
    }
#endif
    G4MT_daughters = nullptr;
    G4Exception("G4VDecayChannel::FillDaughters()",
                "PART011", FatalException,
                "Cannot fill daughters: numberOfDaughters is not defined yet");
  } 
 
  // create and set the array of pointers to daughter particles
  G4MT_daughters = new G4ParticleDefinition*[numberOfDaughters];
  if (G4MT_daughters_mass != nullptr) delete [] G4MT_daughters_mass;
  if (G4MT_daughters_width != nullptr) delete [] G4MT_daughters_width;
  G4MT_daughters_mass = new G4double[numberOfDaughters];
  G4MT_daughters_width = new G4double[numberOfDaughters];
  // loop over all daughters
  for (index=0; index<numberOfDaughters; ++index)
  { 
    if (daughters_name[index] == nullptr)
    {
      // daughter name is not defined
#ifdef G4VERBOSE
      if (verboseLevel>0)
      {
        G4cout << "G4VDecayChannel::FillDaughters() - "
               << "[ " << G4MT_parent->GetParticleName() << " ]"
               << index << "-th daughter is not defined yet" << G4endl;
      }
#endif
      G4MT_daughters[index] = nullptr;
      G4Exception("G4VDecayChannel::FillDaughters()",
                  "PART011", FatalException,
                  "Cannot fill daughters: name of daughter is not defined yet");
    } 
    // search daughter particles in the particle table 
    G4MT_daughters[index] = particletable->FindParticle(*daughters_name[index]);
    if (G4MT_daughters[index] == nullptr )
    {
      // cannot find the daughter particle
#ifdef G4VERBOSE
      if (verboseLevel>0)
      {
        G4cout << "G4VDecayChannel::FillDaughters() - "
                << "[ " << G4MT_parent->GetParticleName() << " ]"
               << index << ":" << *daughters_name[index]
               << " is not defined !!" << G4endl;
        G4cout << " The BR of this decay mode is set to zero." << G4endl;
      }
#endif
      SetBR(0.0);
      return;
    }
#ifdef G4VERBOSE
    if (verboseLevel>1)
    {
      G4cout << index << ":" << *daughters_name[index];
      G4cout << ":" << G4MT_daughters[index] << G4endl;
    }
#endif
    G4MT_daughters_mass[index] = G4MT_daughters[index]->GetPDGMass();
    G4double d_width = G4MT_daughters[index]->GetPDGWidth();
    G4MT_daughters_width[index] = d_width;
    sumofdaughtermass += G4MT_daughters[index]->GetPDGMass();
    sumofdaughterwidthsq += d_width*d_width;
  }  // end loop over all daughters
 
  // check sum of daghter mass
  G4double widthMass = std::sqrt(G4MT_parent->GetPDGWidth()
                                *G4MT_parent->GetPDGWidth()
                                +sumofdaughterwidthsq);
  if ( (G4MT_parent->GetParticleType() != "nucleus")
    && (numberOfDaughters !=1)
    && (sumofdaughtermass > parentmass + rangeMass*widthMass) )
  {
    // !!! illegal mass  !!!
#ifdef G4VERBOSE
    if (GetVerboseLevel()>0)
    {
      G4cout << "G4VDecayChannel::FillDaughters() - "
             << "[ " << G4MT_parent->GetParticleName() << " ]"
             << "    Energy/Momentum conserevation breaks " << G4endl;
     if (GetVerboseLevel()>1)
     {
       G4cout << "    parent:" << *parent_name
              << " mass:" << parentmass/GeV << "[GeV/c/c]" << G4endl;
       for (index=0; index<numberOfDaughters; ++index)
       {
         G4cout << "     daughter " << index << ":" << *daughters_name[index]
                << " mass:" << G4MT_daughters[index]->GetPDGMass()/GeV
                << "[GeV/c/c]" << G4endl;
       }
     }
   }
#endif
 }
}

References G4VDecayChannel::CheckAndFillParent(), G4VDecayChannel::daughters_name, G4VDecayChannel::daughtersMutex, FatalException, G4ParticleTable::FindParticle(), G4cout, G4endl, G4Exception(), G4VDecayChannel::G4MT_daughters, G4VDecayChannel::G4MT_daughters_mass, G4VDecayChannel::G4MT_daughters_width, G4VDecayChannel::G4MT_parent, G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetParticleType(), G4ParticleDefinition::GetPDGMass(), G4ParticleDefinition::GetPDGWidth(), G4VDecayChannel::GetVerboseLevel(), GeV, G4VDecayChannel::numberOfDaughters, G4VDecayChannel::parent_name, G4VDecayChannel::particletable, G4VDecayChannel::rangeMass, G4VDecayChannel::SetBR(), and G4VDecayChannel::verboseLevel.

Referenced by G4VDecayChannel::CheckAndFillDaughters().

◆ FillParent()

void G4VDecayChannel::FillParent ( )

privateinherited

Definition at line 440 of file G4VDecayChannel.cc.

{
  G4AutoLock lock(&parentMutex);
 
  // Double check, check again if another thread has already filled this, in
  // case do not need to do anything
  if ( G4MT_parent != nullptr ) return;
 
  if (parent_name == nullptr)
  {
    // parent name is not defined
#ifdef G4VERBOSE
    if (verboseLevel>0)
    {
      G4cout << "G4VDecayChannel::FillParent() - "
             << "parent name is not defined !!" << G4endl;
    }
#endif
    G4MT_parent = nullptr;
    G4Exception("G4VDecayChannel::FillParent()",
                "PART012", FatalException,
                "Cannot fill parent: parent name is not defined yet");    
    return;
  }
 
  // search parent particle in the particle table
  G4MT_parent = particletable->FindParticle(*parent_name);
  if (G4MT_parent == nullptr)
  {
    // parent particle does not exist
#ifdef G4VERBOSE
    if (verboseLevel>0)
    {
      G4cout << "G4VDecayChannel::FillParent() - "
             << *parent_name << " does not exist !!" << G4endl;
    }
#endif
    G4Exception("G4VDecayChannel::FillParent()",
                "PART012", FatalException,
                "Cannot fill parent: parent does not exist");
    return;
  }
  G4MT_parent_mass = G4MT_parent->GetPDGMass();
}

References FatalException, G4ParticleTable::FindParticle(), G4cout, G4endl, G4Exception(), G4VDecayChannel::G4MT_parent, G4VDecayChannel::G4MT_parent_mass, G4ParticleDefinition::GetPDGMass(), G4VDecayChannel::parent_name, G4VDecayChannel::parentMutex, G4VDecayChannel::particletable, and G4VDecayChannel::verboseLevel.

Referenced by G4VDecayChannel::CheckAndFillParent().

◆ GetAngularMomentum()

G4int G4VDecayChannel::GetAngularMomentum ( )

inherited

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;
}

References G4VDecayChannel::CheckAndFillDaughters(), G4cout, G4endl, G4Exception(), G4VDecayChannel::G4MT_daughters, G4VDecayChannel::G4MT_parent, G4ParticleDefinition::GetPDGiParity(), G4ParticleDefinition::GetPDGiSpin(), JustWarning, G4VDecayChannel::numberOfDaughters, and G4VDecayChannel::verboseLevel.

◆ GetBR()

G4double G4VDecayChannel::GetBR ( ) const

inlineinherited

◆ GetDaughter()

G4ParticleDefinition * G4VDecayChannel::GetDaughter ( G4int anIndex )

inlineinherited

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;
  }
}

References G4VDecayChannel::CheckAndFillDaughters(), G4cout, G4endl, G4VDecayChannel::G4MT_daughters, G4VDecayChannel::numberOfDaughters, and G4VDecayChannel::verboseLevel.

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

◆ GetDaughterMass()

G4double G4VDecayChannel::GetDaughterMass ( G4int anIndex ) const

inlineinherited

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;
  }
}

References G4cout, G4endl, G4VDecayChannel::G4MT_daughters_mass, G4VDecayChannel::numberOfDaughters, and G4VDecayChannel::verboseLevel.

◆ GetDaughterName()

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

inlineinherited

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();
  }
}

References G4VDecayChannel::daughters_name, G4cout, G4endl, G4VDecayChannel::GetNoName(), G4VDecayChannel::numberOfDaughters, and G4VDecayChannel::verboseLevel.

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

◆ GetKinematicsName()

const G4String & G4VDecayChannel::GetKinematicsName ( ) const

inlineinherited

Definition at line 292 of file G4VDecayChannel.hh.

{
  return kinematics_name;
}

References G4VDecayChannel::kinematics_name.

Referenced by G4HtmlPPReporter::GeneratePropertyTable(), and G4TextPPReporter::GeneratePropertyTable().

◆ GetNoName()

const G4String & G4VDecayChannel::GetNoName ( ) const

privateinherited

Definition at line 579 of file G4VDecayChannel.cc.

{
  return noName;
}

References G4VDecayChannel::noName.

Referenced by G4VDecayChannel::GetDaughterName().

◆ GetNumberOfDaughters()

G4int G4VDecayChannel::GetNumberOfDaughters ( ) const

inlineinherited

Definition at line 286 of file G4VDecayChannel.hh.

{ 
  return  numberOfDaughters;
}

References G4VDecayChannel::numberOfDaughters.

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

◆ GetParent()

G4ParticleDefinition * G4VDecayChannel::GetParent ( )

inlineinherited

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;
}

References G4VDecayChannel::CheckAndFillParent(), and G4VDecayChannel::G4MT_parent.

Referenced by G4DecayTable::Insert().

◆ GetParentMass()

G4double G4VDecayChannel::GetParentMass ( ) const

inlineinherited

Definition at line 272 of file G4VDecayChannel.hh.

{
  return G4MT_parent_mass;
}

References G4VDecayChannel::G4MT_parent_mass.

◆ GetParentName()

const G4String & G4VDecayChannel::GetParentName ( ) const

inlineinherited

Definition at line 266 of file G4VDecayChannel.hh.

{
  return *parent_name;
}

References G4VDecayChannel::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

inlineinherited

Definition at line 331 of file G4VDecayChannel.hh.

{
  return parent_polarization;
}

References G4VDecayChannel::parent_polarization.

◆ GetRangeMass()

G4double G4VDecayChannel::GetRangeMass ( ) const

inlineinherited

Definition at line 316 of file G4VDecayChannel.hh.

{
  return rangeMass;
}

References G4VDecayChannel::rangeMass.

◆ GetVerboseLevel()

G4int G4VDecayChannel::GetVerboseLevel ( ) const

inlineinherited

Definition at line 310 of file G4VDecayChannel.hh.

{
  return verboseLevel;
}

References G4VDecayChannel::verboseLevel.

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

◆ IsOKWithParentMass()

G4bool G4VDecayChannel::IsOKWithParentMass ( G4double parentMass )

virtualinherited

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); 
}

References G4VDecayChannel::CheckAndFillDaughters(), G4VDecayChannel::CheckAndFillParent(), G4VDecayChannel::G4MT_daughters_mass, G4VDecayChannel::G4MT_daughters_width, G4VDecayChannel::numberOfDaughters, and G4VDecayChannel::rangeMass.

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

◆ operator!=()

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

inlineinherited

Definition at line 69 of file G4VDecayChannel.hh.

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

◆ operator<()

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

inlineinherited

Definition at line 195 of file G4VDecayChannel.hh.

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

References G4VDecayChannel::rbranch.

◆ operator=()

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

protected

Definition at line 68 of file G4MuonDecayChannelWithSpin.cc.

{
  if (this != &right)
  { 
    kinematics_name = right.kinematics_name;
    verboseLevel = right.verboseLevel;
    rbranch = right.rbranch;
 
    // 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]);
      }
    }
  }
  return *this;
}

References G4VDecayChannel::ClearDaughtersName(), G4VDecayChannel::daughters_name, G4VDecayChannel::kinematics_name, G4VDecayChannel::numberOfDaughters, G4VDecayChannel::parent_name, G4VDecayChannel::rbranch, and G4VDecayChannel::verboseLevel.

◆ operator==()

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

inlineinherited

Definition at line 68 of file G4VDecayChannel.hh.

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

◆ R_c()

G4double G4MuonDecayChannelWithSpin::R_c	(	G4double	x,
		G4double	omega
	)

private

Definition at line 304 of file G4MuonDecayChannelWithSpin.cc.

{
  G4int n_max = (G4int)(100.*x);
 
  if(n_max<10)n_max=10;
 
  G4double L2 = 0.0;
 
  for(G4int n=1; n<=n_max; ++n)
  {
    L2 += std::pow(x,n)/(n*n);
  }
 
  G4double r_c;
 
  r_c = 2.*L2-(pi*pi/3.)-2.;
  r_c = r_c + omega * (1.5+2.*std::log((1.-x)/x));
  r_c = r_c - std::log(x)*(2.*std::log(x)-1.);
  r_c = r_c + (3.*std::log(x)-1.-1./x)*std::log(1.-x);
 
  return r_c;
}

References CLHEP::detail::n, and pi.

Referenced by F_c(), and F_theta().

◆ SetBR()

void G4VDecayChannel::SetBR ( G4double value )

inherited

Definition at line 624 of file G4VDecayChannel.cc.

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

References G4VDecayChannel::rbranch.

Referenced by G4TauMinus::Definition(), G4TauPlus::Definition(), G4VDecayChannel::FillDaughters(), 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
	)

inherited

Definition at line 300 of file G4VDecayChannel.cc.

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

References G4ParticleDefinition::GetParticleName(), and G4VDecayChannel::SetDaughter().

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 G4VDecayChannel::SetDaughter().

◆ SetDaughter() [2/2]

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

inherited

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
  }
}

References G4VDecayChannel::daughters_name, FatalException, G4cout, G4endl, G4Exception(), G4VDecayChannel::G4MT_daughters, G4VDecayChannel::numberOfDaughters, and G4VDecayChannel::verboseLevel.

◆ SetNumberOfDaughters()

void G4VDecayChannel::SetNumberOfDaughters ( G4int value )

inherited

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;
  }
}

References G4VDecayChannel::ClearDaughtersName(), G4VDecayChannel::daughters_name, and G4VDecayChannel::numberOfDaughters.

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 )

inherited

Definition at line 486 of file G4VDecayChannel.cc.

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

References G4ParticleDefinition::GetParticleName(), and G4VDecayChannel::SetParent().

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 G4VDecayChannel::SetParent().

◆ SetParent() [2/2]

void G4VDecayChannel::SetParent ( const G4String & particle_name )

inlineinherited

Definition at line 278 of file G4VDecayChannel.hh.

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

References G4VDecayChannel::G4MT_parent, and G4VDecayChannel::parent_name.

◆ SetPolarization()

void G4VDecayChannel::SetPolarization ( const G4ThreeVector & polar )

inlineinherited

Definition at line 325 of file G4VDecayChannel.hh.

{
  parent_polarization = polar;
}

References G4VDecayChannel::parent_polarization.

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

◆ SetRangeMass()

void G4VDecayChannel::SetRangeMass ( G4double val )

inlineinherited

Definition at line 322 of file G4VDecayChannel.hh.

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

References G4VDecayChannel::rangeMass.

◆ SetVerboseLevel()

void G4VDecayChannel::SetVerboseLevel ( G4int value )

inlineinherited

Definition at line 304 of file G4VDecayChannel.hh.

{
  verboseLevel = value;
}

References G4VDecayChannel::verboseLevel.

Referenced by G4Decay::DecayIt(), and G4MuonicAtomDecay::DecayIt().

Field Documentation

◆ daughters_name

G4String** G4VDecayChannel::daughters_name = nullptr

protectedinherited

◆ daughtersMutex

G4Mutex G4VDecayChannel::daughtersMutex

protectedinherited

Definition at line 177 of file G4VDecayChannel.hh.

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

◆ G4MT_daughters

G4ParticleDefinition** G4VDecayChannel::G4MT_daughters = nullptr

protectedinherited

Definition at line 173 of file G4VDecayChannel.hh.

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

◆ G4MT_daughters_mass

G4double* G4VDecayChannel::G4MT_daughters_mass = nullptr

protectedinherited

Definition at line 175 of file G4VDecayChannel.hh.

Referenced by G4VDecayChannel::ClearDaughtersName(), G4VDecayChannel::FillDaughters(), G4VDecayChannel::GetDaughterMass(), G4VDecayChannel::IsOKWithParentMass(), G4PhaseSpaceDecayChannel::ManyBodyDecayIt(), G4PhaseSpaceDecayChannel::ThreeBodyDecayIt(), G4PhaseSpaceDecayChannel::TwoBodyDecayIt(), and G4VDecayChannel::~G4VDecayChannel().

◆ G4MT_daughters_width

G4double* G4VDecayChannel::G4MT_daughters_width = nullptr

protectedinherited

Definition at line 176 of file G4VDecayChannel.hh.

Referenced by G4VDecayChannel::ClearDaughtersName(), G4VDecayChannel::FillDaughters(), G4VDecayChannel::IsOKWithParentMass(), G4PhaseSpaceDecayChannel::ThreeBodyDecayIt(), G4PhaseSpaceDecayChannel::TwoBodyDecayIt(), and G4VDecayChannel::~G4VDecayChannel().

◆ G4MT_parent

G4ParticleDefinition* G4VDecayChannel::G4MT_parent = nullptr

protectedinherited

Definition at line 172 of file G4VDecayChannel.hh.

Referenced by G4VDecayChannel::CheckAndFillParent(), G4DalitzDecayChannel::DecayIt(), G4KL3DecayChannel::DecayIt(), G4MuonDecayChannel::DecayIt(), 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(), G4VDecayChannel::FillDaughters(), G4VDecayChannel::FillParent(), G4GeneralPhaseSpaceDecay::G4GeneralPhaseSpaceDecay(), G4VDecayChannel::GetAngularMomentum(), G4VDecayChannel::GetParent(), G4PhaseSpaceDecayChannel::ManyBodyDecayIt(), G4GeneralPhaseSpaceDecay::ManyBodyDecayIt(), G4PhaseSpaceDecayChannel::OneBodyDecayIt(), G4GeneralPhaseSpaceDecay::OneBodyDecayIt(), G4VDecayChannel::SetParent(), G4PhaseSpaceDecayChannel::ThreeBodyDecayIt(), G4GeneralPhaseSpaceDecay::ThreeBodyDecayIt(), G4PhaseSpaceDecayChannel::TwoBodyDecayIt(), and G4GeneralPhaseSpaceDecay::TwoBodyDecayIt().

◆ G4MT_parent_mass

G4double G4VDecayChannel::G4MT_parent_mass = 0.0

protectedinherited

Definition at line 174 of file G4VDecayChannel.hh.

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

◆ kinematics_name

G4String G4VDecayChannel::kinematics_name = ""

protectedinherited

Definition at line 152 of file G4VDecayChannel.hh.

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

◆ noName

const G4String G4VDecayChannel::noName = " "

staticprotectedinherited

Definition at line 170 of file G4VDecayChannel.hh.

Referenced by G4VDecayChannel::GetNoName().

◆ numberOfDaughters

G4int G4VDecayChannel::numberOfDaughters = 0

protectedinherited

Definition at line 180 of file G4VDecayChannel.hh.

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

◆ parent_name

G4String* G4VDecayChannel::parent_name = nullptr

protectedinherited

◆ parent_polarization

G4ThreeVector G4VDecayChannel::parent_polarization

protectedinherited

Definition at line 164 of file G4VDecayChannel.hh.

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

◆ parentMutex

G4Mutex G4VDecayChannel::parentMutex

protectedinherited

Definition at line 178 of file G4VDecayChannel.hh.

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

◆ particletable

G4ParticleTable* G4VDecayChannel::particletable = nullptr

protectedinherited

Definition at line 167 of file G4VDecayChannel.hh.

Referenced by G4VDecayChannel::FillDaughters(), G4VDecayChannel::FillParent(), G4VDecayChannel::G4VDecayChannel(), and G4VDecayChannel::operator=().

◆ rangeMass

G4double G4VDecayChannel::rangeMass = 2.5

protectedinherited

Definition at line 161 of file G4VDecayChannel.hh.

Referenced by G4VDecayChannel::DynamicalMass(), G4VDecayChannel::FillDaughters(), G4VDecayChannel::G4VDecayChannel(), G4VDecayChannel::GetRangeMass(), G4VDecayChannel::IsOKWithParentMass(), G4VDecayChannel::operator=(), G4VDecayChannel::SetRangeMass(), G4PhaseSpaceDecayChannel::ThreeBodyDecayIt(), and G4PhaseSpaceDecayChannel::TwoBodyDecayIt().

◆ rbranch

G4double G4VDecayChannel::rbranch = 0.0

protectedinherited

Definition at line 154 of file G4VDecayChannel.hh.

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

◆ verboseLevel

G4int G4VDecayChannel::verboseLevel = 1

protectedinherited

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

source/particles/management/include/G4MuonDecayChannelWithSpin.hh
source/particles/management/src/G4MuonDecayChannelWithSpin.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ G4MuonDecayChannelWithSpin() [1/3]

◆ ~G4MuonDecayChannelWithSpin()

◆ G4MuonDecayChannelWithSpin() [2/3]

◆ G4MuonDecayChannelWithSpin() [3/3]

Member Function Documentation

◆ CheckAndFillDaughters()

◆ CheckAndFillParent()

◆ ClearDaughtersName()

◆ DecayIt()

◆ DumpInfo()

◆ DynamicalMass()

◆ F_c()

◆ F_theta()

◆ FillDaughters()

◆ FillParent()

◆ GetAngularMomentum()

◆ GetBR()

◆ GetDaughter()

◆ GetDaughterMass()

◆ GetDaughterName()

◆ GetKinematicsName()

◆ GetNoName()

◆ GetNumberOfDaughters()

◆ GetParent()

◆ GetParentMass()

◆ GetParentName()

◆ GetPolarization()

◆ GetRangeMass()

◆ GetVerboseLevel()

◆ IsOKWithParentMass()

◆ operator!=()

◆ operator<()

◆ operator=()

◆ operator==()

◆ R_c()

◆ SetBR()

◆ SetDaughter() [1/2]

◆ SetDaughter() [2/2]

◆ SetNumberOfDaughters()

◆ SetParent() [1/2]

◆ SetParent() [2/2]

◆ SetPolarization()

◆ SetRangeMass()

◆ SetVerboseLevel()

Field 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