#include <G4TauPlus.hh>

Inheritance diagram for G4TauPlus:

Static Public Member Functions
static G4TauPlus *	Definition ()

static G4TauPlus *	TauPlusDefinition ()

static G4TauPlus *	TauPlus ()

Static Public Member Functions inherited from G4ParticleDefinition
static const G4PDefManager &	GetSubInstanceManager ()

static void	Clean ()

Additional Inherited Members
Public Member Functions inherited from G4ParticleDefinition
	G4ParticleDefinition (const G4String &aName, G4double mass, G4double width, G4double charge, G4int iSpin, G4int iParity, G4int iConjugation, G4int iIsospin, G4int iIsospinZ, G4int gParity, const G4String &pType, G4int lepton, G4int baryon, G4int encoding, G4bool stable, G4double lifetime, G4DecayTable *decaytable, G4bool shortlived=false, const G4String &subType="", G4int anti_encoding=0, G4double magneticMoment=0.0)

virtual	~G4ParticleDefinition ()

	G4ParticleDefinition (const G4ParticleDefinition &)=delete

G4ParticleDefinition &	operator= (const G4ParticleDefinition &)=delete

G4bool	operator== (const G4ParticleDefinition &right) const

G4bool	operator!= (const G4ParticleDefinition &right) const

const G4String &	GetParticleName () const

G4double	GetPDGMass () const

G4double	GetPDGWidth () const

G4double	GetPDGCharge () const

G4double	GetPDGSpin () const

G4int	GetPDGiSpin () const

G4int	GetPDGiParity () const

G4int	GetPDGiConjugation () const

G4double	GetPDGIsospin () const

G4double	GetPDGIsospin3 () const

G4int	GetPDGiIsospin () const

G4int	GetPDGiIsospin3 () const

G4int	GetPDGiGParity () const

G4double	GetPDGMagneticMoment () const

void	SetPDGMagneticMoment (G4double mageticMoment)

G4double	CalculateAnomaly () const

const G4String &	GetParticleType () const

const G4String &	GetParticleSubType () const

G4int	GetLeptonNumber () const

G4int	GetBaryonNumber () const

G4int	GetPDGEncoding () const

G4int	GetAntiPDGEncoding () const

void	SetAntiPDGEncoding (G4int aEncoding)

G4int	GetQuarkContent (G4int flavor) const

G4int	GetAntiQuarkContent (G4int flavor) const

G4bool	IsShortLived () const

G4bool	GetPDGStable () const

void	SetPDGStable (const G4bool aFlag)

G4double	GetPDGLifeTime () const

void	SetPDGLifeTime (G4double aLifeTime)

G4double	GetIonLifeTime () const

G4DecayTable *	GetDecayTable () const

void	SetDecayTable (G4DecayTable *aDecayTable)

G4ProcessManager *	GetProcessManager () const

void	SetProcessManager (G4ProcessManager *aProcessManager)

G4VTrackingManager *	GetTrackingManager () const

void	SetTrackingManager (G4VTrackingManager *aTrackingManager)

G4ParticleTable *	GetParticleTable () const

G4int	GetAtomicNumber () const

G4int	GetAtomicMass () const

void	DumpTable () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

void	SetApplyCutsFlag (G4bool)

G4bool	GetApplyCutsFlag () const

G4bool	IsGeneralIon () const

G4bool	IsMuonicAtom () const

G4ProcessManager *	GetMasterProcessManager () const

void	SetMasterProcessManager (G4ProcessManager *aNewPM)

G4int	GetInstanceID () const

void	SetParticleDefinitionID (G4int id=-1)

G4int	GetParticleDefinitionID () const

G4bool	IsHypernucleus () const

G4int	GetNumberOfLambdasInHypernucleus () const

G4bool	IsAntiHypernucleus () const

G4int	GetNumberOfAntiLambdasInAntiHypernucleus () const

Protected Types inherited from G4ParticleDefinition
enum	{ NumberOfQuarkFlavor = 6 }

Protected Member Functions inherited from G4ParticleDefinition
	G4ParticleDefinition ()

G4int	FillQuarkContents ()

void	SetParticleSubType (const G4String &subtype)

void	SetAtomicNumber (G4int)

void	SetAtomicMass (G4int)

Protected Attributes inherited from G4ParticleDefinition
G4int	theQuarkContent [NumberOfQuarkFlavor]

G4int	theAntiQuarkContent [NumberOfQuarkFlavor]

G4bool	isGeneralIon = false

G4bool	isMuonicAtom = false

Detailed Description

Definition at line 40 of file G4TauPlus.hh.

Member Function Documentation

◆ Definition()

G4TauPlus * G4TauPlus::Definition ( )

static

Definition at line 50 of file G4TauPlus.cc.

{
  if (theInstance != nullptr) return theInstance;
  const G4String name = "tau+";
  // search in particle table]
  G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
  G4ParticleDefinition* anInstance = pTable->FindParticle(name);
  if (anInstance == nullptr) {
    // create particle
    //
    //    Arguments for constructor are as follows
    //               name             mass          width         charge
    //             2*spin           parity  C-conjugation
    //          2*Isospin       2*Isospin3       G-parity
    //               type    lepton number  baryon number   PDG encoding
    //             stable         lifetime    decay table
    //             shortlived      subType    anti_encoding
 
    // clang-format off
   anInstance = new G4ParticleDefinition(
                 name,     1.77686*GeV,  2.267e-9*MeV,     1.*eplus, 
            1,               0,             0,          
            0,               0,             0,             
         "lepton",              -1,             0,         -15,
        false,     290.3e-6*ns,          nullptr,
             false,           "tau"
              );
    // clang-format on
 
    // Bohr Magnetron
    G4double muB = 0.5 * eplus * hbar_Planck / (anInstance->GetPDGMass() / c_squared);
 
    anInstance->SetPDGMagneticMoment(muB * 1.00118);
 
    // create Decay Table
    auto table = new G4DecayTable();
 
    // create decay channels
    G4VDecayChannel* mode;
    // tau+ -> mu+ + nu_mu + anti_nu_tau
    mode = new G4TauLeptonicDecayChannel("tau+", 0.1736, "mu+");
    table->Insert(mode);
    // tau+ -> e+ + nu_e + anti_nu_tau
    mode = new G4TauLeptonicDecayChannel("tau+", 0.1784, "e+");
    table->Insert(mode);
    // tau+ -> pi+ + anti_nu_tau
    mode = new G4PhaseSpaceDecayChannel("tau+", 0.1106, 2, "pi+", "anti_nu_tau");
    table->Insert(mode);
    // tau+ -> pi0 + pi0 + pi+ + anti_nu_tau
    mode = new G4PhaseSpaceDecayChannel("tau+", 0.2541, 3, "pi0", "pi+", "anti_nu_tau");
    table->Insert(mode);
    // tau+ -> pi0 + pi0 + pi+ + anti_nu_tau
    mode = new G4PhaseSpaceDecayChannel();
    mode->SetParent("tau+");
    mode->SetBR(0.0917);
    mode->SetNumberOfDaughters(4);
    mode->SetDaughter(0, "pi0");
    mode->SetDaughter(1, "pi0");
    mode->SetDaughter(2, "pi+");
    mode->SetDaughter(3, "anti_nu_tau");
    table->Insert(mode);
    // tau+ -> pi+ + pi+ + pi- + anti_nu_tau
    mode = new G4PhaseSpaceDecayChannel();
    mode->SetParent("tau+");
    mode->SetBR(0.0946);
    mode->SetNumberOfDaughters(4);
    mode->SetDaughter(0, "pi+");
    mode->SetDaughter(1, "pi+");
    mode->SetDaughter(2, "pi-");
    mode->SetDaughter(3, "anti_nu_tau");
    table->Insert(mode);
 
    anInstance->SetDecayTable(table);
  }
  theInstance = static_cast<G4TauPlus*>(anInstance);
  return theInstance;
}