G4TauPlus.cc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 //
00027 // $Id$
00028 //
00029 // 
00030 // ----------------------------------------------------------------------
00031 //      GEANT 4 class implementation file
00032 //
00033 //      History: first implementation, based on object model of
00034 //      4th April 1996, G.Cosmo
00035 //      7 July 1996                   H.Kurashige 
00036 // **********************************************************************
00037 //  New impelemenataion as an utility class  M.Asai, 26 July 2004
00038 // ----------------------------------------------------------------------
00039 
00040 #include "G4TauPlus.hh"
00041 #include "G4PhysicalConstants.hh"
00042 #include "G4SystemOfUnits.hh"
00043 #include "G4ParticleTable.hh"
00044 
00045 #include "G4PhaseSpaceDecayChannel.hh"
00046 #include "G4TauLeptonicDecayChannel.hh"
00047 #include "G4DecayTable.hh"
00048 
00049 // ######################################################################
00050 // ###                          TAPLUS                                ###
00051 // ######################################################################
00052 G4TauPlus* G4TauPlus::theInstance = 0;
00053 
00054 G4TauPlus* G4TauPlus::Definition()
00055 {
00056   if (theInstance !=0) return theInstance;
00057   const G4String name = "tau+";
00058   // search in particle table]
00059   G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
00060   G4ParticleDefinition* anInstance = pTable->FindParticle(name);
00061   if (anInstance ==0)
00062   {
00063   // create particle
00064   //
00065   //    Arguments for constructor are as follows
00066   //               name             mass          width         charge
00067   //             2*spin           parity  C-conjugation
00068   //          2*Isospin       2*Isospin3       G-parity
00069   //               type    lepton number  baryon number   PDG encoding
00070   //             stable         lifetime    decay table
00071   //             shortlived      subType    anti_encoding
00072    anInstance = new G4ParticleDefinition(
00073                  name,     1.77682*GeV,  2.265e-9*MeV,     1.*eplus, 
00074                     1,               0,             0,          
00075                     0,               0,             0,             
00076              "lepton",              -1,             0,         -15,
00077                 false,     290.6e-6*ns,          NULL,
00078              false,           "tau"
00079               );
00080    // Bohr Magnetron
00081    G4double muB =  0.5*eplus*hbar_Planck/(1.77699*GeV/c_squared) ;
00082    
00083    anInstance->SetPDGMagneticMoment( muB * 2.* 1.00116);
00084 
00085   //create Decay Table 
00086   G4DecayTable* table = new G4DecayTable();
00087 
00088   // create decay channels
00089   G4VDecayChannel* mode;
00090   // tau+ -> mu+ + nu_mu + anti_nu_tau
00091   mode = new G4TauLeptonicDecayChannel("tau+",0.1736,"mu+");
00092   table->Insert(mode);
00093   // tau+ -> e+ + nu_e + anti_nu_tau
00094   mode = new G4TauLeptonicDecayChannel("tau+",0.1784,"e+");
00095   table->Insert(mode);
00096   // tau+ -> pi+ + anti_nu_tau
00097   mode = new G4PhaseSpaceDecayChannel("tau+",0.1106,2,"pi+","anti_nu_tau");
00098   table->Insert(mode);
00099   // tau+ -> pi0 + pi0 + pi+ + anti_nu_tau
00100   mode = new G4PhaseSpaceDecayChannel("tau+",0.2541,3,"pi0","pi+","anti_nu_tau");
00101   table->Insert(mode);
00102   // tau+ -> pi0 + pi0 + pi+ + anti_nu_tau
00103   mode = new G4PhaseSpaceDecayChannel();
00104   mode->SetParent("tau+");
00105   mode->SetBR(0.0917);
00106   mode->SetNumberOfDaughters(4);
00107   mode->SetDaughter(0,"pi0");
00108   mode->SetDaughter(1,"pi0");
00109   mode->SetDaughter(2,"pi+");
00110   mode->SetDaughter(3,"anti_nu_tau");
00111   table->Insert(mode);
00112   // tau+ -> pi+ + pi+ + pi- + anti_nu_tau
00113   mode = new G4PhaseSpaceDecayChannel();
00114   mode->SetParent("tau+");
00115   mode->SetBR(0.0946);
00116   mode->SetNumberOfDaughters(4);
00117   mode->SetDaughter(0,"pi+");
00118   mode->SetDaughter(1,"pi+");
00119   mode->SetDaughter(2,"pi-");
00120   mode->SetDaughter(3,"anti_nu_tau");
00121   table->Insert(mode);
00122 
00123    anInstance->SetDecayTable(table);
00124   }
00125   theInstance = reinterpret_cast<G4TauPlus*>(anInstance);
00126   return theInstance;
00127 }
00128 
00129 G4TauPlus*  G4TauPlus::TauPlusDefinition()
00130 {
00131   return Definition();
00132 }
00133 
00134 G4TauPlus*  G4TauPlus::TauPlus()
00135 {
00136   return Definition();
00137 }
00138 

Generated on Mon May 27 17:49:56 2013 for Geant4 by  doxygen 1.4.7