G4TauMinus.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 "G4TauMinus.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 // ###                          TAUMINUS                             ###
00051 // ######################################################################
00052 G4TauMinus* G4TauMinus::theInstance = 0;
00053 
00054 G4TauMinus* G4TauMinus::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 
00081    // Bohr Magnetron
00082    G4double muB =  -0.5*eplus*hbar_Planck/(1.77699*GeV/c_squared) ;
00083    
00084    anInstance->SetPDGMagneticMoment( muB * 2. * 1.00116);
00085 
00086   //create Decay Table 
00087   G4DecayTable* table = new G4DecayTable();
00088 
00089   // create decay channels
00090   G4VDecayChannel* mode;
00091   // tau- -> mu- + anti_nu_mu + nu_tau
00092   mode = new G4TauLeptonicDecayChannel("tau-",0.1736,"mu-");
00093   table->Insert(mode);
00094   // tau- -> e- + anti_nu_e + nu_tau
00095   mode = new G4TauLeptonicDecayChannel("tau-",0.1784,"e-");
00096   table->Insert(mode);
00097   // tau- -> pi- + nu_tau
00098   mode = new G4PhaseSpaceDecayChannel("tau-",0.1106,2,"pi-","nu_tau");
00099   table->Insert(mode);
00100   // tau- -> pi0 + pi- + nu_tau
00101   mode = new G4PhaseSpaceDecayChannel("tau-",0.2541,3,"pi0","pi-","nu_tau");
00102   table->Insert(mode);
00103   // tau- -> pi0 + pi0 + pi- + nu_tau
00104   mode = new G4PhaseSpaceDecayChannel();
00105   mode->SetParent("tau-");
00106   mode->SetBR(0.0917);
00107   mode->SetNumberOfDaughters(4);
00108   mode->SetDaughter(0,"pi0");
00109   mode->SetDaughter(1,"pi0");
00110   mode->SetDaughter(2,"pi-");
00111   mode->SetDaughter(3,"nu_tau");
00112   table->Insert(mode);
00113   // tau- -> pi- + pi- + pi+ + nu_tau
00114   mode = new G4PhaseSpaceDecayChannel();
00115   mode->SetParent("tau-");
00116   mode->SetBR(0.0946);
00117   mode->SetNumberOfDaughters(4);
00118   mode->SetDaughter(0,"pi-");
00119   mode->SetDaughter(1,"pi-");
00120   mode->SetDaughter(2,"pi+");
00121   mode->SetDaughter(3,"nu_tau");
00122   table->Insert(mode);
00123 
00124   anInstance->SetDecayTable(table);
00125   }
00126   theInstance = reinterpret_cast<G4TauMinus*>(anInstance);
00127   return theInstance;
00128 }
00129 
00130 G4TauMinus*  G4TauMinus::TauMinusDefinition()
00131 {
00132   return Definition();
00133 }
00134 
00135 G4TauMinus*  G4TauMinus::TauMinus()
00136 {
00137   return Definition();
00138 }
00139 
00140 

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