#include <G4TauMinus.hh>

Inheritance diagram for G4TauMinus:

Collaboration diagram for G4TauMinus:

Static Public Member Functions
static G4TauMinus *	Definition ()

static G4TauMinus *	TauMinusDefinition ()

static G4TauMinus *	TauMinus ()

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

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

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)

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

G4int	operator== (const G4ParticleDefinition &right) const

G4int	operator!= (const G4ParticleDefinition &right) const

G4ProcessManager *	GetMasterProcessManager () const

void	SetMasterProcessManager (G4ProcessManager *aNewPM)

G4int	GetInstanceID () const

void	SetParticleDefinitionID (G4int id=-1)

G4int	GetParticleDefinitionID () const

Protected Types inherited from G4ParticleDefinition
enum	{ NumberOfQuarkFlavor = 6 }

Protected Member Functions inherited from G4ParticleDefinition
G4int	FillQuarkContents ()

void	SetParticleSubType (const G4String &subtype)

void	SetAtomicNumber (G4int)

void	SetAtomicMass (G4int)

	G4ParticleDefinition (const G4ParticleDefinition &right)

	G4ParticleDefinition ()

Protected Attributes inherited from G4ParticleDefinition
G4int	theQuarkContent [NumberOfQuarkFlavor]

G4int	theAntiQuarkContent [NumberOfQuarkFlavor]

G4bool	isGeneralIon

Detailed Description

Definition at line 50 of file G4TauMinus.hh.

Member Function Documentation

G4TauMinus * G4TauMinus::Definition ( )

static

Definition at line 54 of file G4TauMinus.cc.

 {
   if (theInstance !=0) return theInstance;
   const G4String name = "tau-";
   // search in particle table]
   G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
   G4ParticleDefinition* anInstance = pTable->FindParticle(name);
   if (anInstance ==0)
   {
   // 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
    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,          NULL,
              false,           "tau"
               );
 
    // Bohr Magnetron
    G4double muB =  -0.5*eplus*hbar_Planck/(anInstance->GetPDGMass()/c_squared) ;
    
    anInstance->SetPDGMagneticMoment( muB * 1.00118);
 
   //create Decay Table 
   G4DecayTable* table = new G4DecayTable();
 
   // create decay channels
   G4VDecayChannel* mode;
   // tau- -> mu- + anti_nu_mu + nu_tau
   mode = new G4TauLeptonicDecayChannel("tau-",0.1736,"mu-");
   table->Insert(mode);
   // tau- -> e- + anti_nu_e + nu_tau
   mode = new G4TauLeptonicDecayChannel("tau-",0.1784,"e-");
   table->Insert(mode);
   // tau- -> pi- + nu_tau
   mode = new G4PhaseSpaceDecayChannel("tau-",0.1106,2,"pi-","nu_tau");
   table->Insert(mode);
   // tau- -> pi0 + pi- + nu_tau
   mode = new G4PhaseSpaceDecayChannel("tau-",0.2541,3,"pi0","pi-","nu_tau");
   table->Insert(mode);
   // tau- -> pi0 + pi0 + pi- + 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,"nu_tau");
   table->Insert(mode);
   // tau- -> pi- + pi- + pi+ + 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,"nu_tau");
   table->Insert(mode);
 
   anInstance->SetDecayTable(table);
   }
   theInstance = reinterpret_cast<G4TauMinus*>(anInstance);
   return theInstance;
 }