G4Eta.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 //                              H.Kurashige   7 Jul 96
00036 // **********************************************************************
00037 //  New impelemenataion as an utility class  M.Asai, 26 July 2004
00038 // ----------------------------------------------------------------------
00039 //      Update mass (PDG2006)              Oct. 11 2006 H.Kurashige 
00040 //
00041 
00042 #include "G4Eta.hh"
00043 #include "G4SystemOfUnits.hh"
00044 #include "G4ParticleTable.hh"
00045 
00046 #include "G4PhaseSpaceDecayChannel.hh"
00047 #include "G4DecayTable.hh"
00048 
00049 // ######################################################################
00050 // ###                         ETA                                    ###
00051 // ######################################################################
00052 
00053 G4Eta* G4Eta::theInstance = 0;
00054 
00055 G4Eta* G4Eta::Definition()
00056 {
00057   if (theInstance !=0) return theInstance;
00058   const G4String name = "eta";
00059   // search in particle table]
00060   G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
00061   G4ParticleDefinition* anInstance = pTable->FindParticle(name);
00062   if (anInstance ==0)
00063   {
00064   // create particle
00065   //
00066   //    Arguments for constructor are as follows
00067   //               name             mass          width         charge
00068   //             2*spin           parity  C-conjugation
00069   //          2*Isospin       2*Isospin3       G-parity
00070   //               type    lepton number  baryon number   PDG encoding
00071   //             stable         lifetime    decay table
00072   //             shortlived      subType    anti_encoding
00073 
00074    anInstance = new G4ParticleDefinition(
00075                  name,    0.547853*GeV,      1.30*keV,         0.0,
00076                     0,              -1,            +1,
00077                     0,               0,            +1,
00078               "meson",               0,             0,         221,
00079                 false,          0.0*ns,          NULL,
00080                 false,           "eta",           221);
00081  //create Decay Table
00082   G4DecayTable* table = new G4DecayTable();
00083 
00084  // create decay channels
00085   G4VDecayChannel** mode = new G4VDecayChannel*[4];
00086   // eta -> gamma + gamma
00087   mode[0] = new G4PhaseSpaceDecayChannel("eta",0.3942,2,"gamma","gamma");
00088   // eta -> pi0 + pi0 + pi0
00089   mode[1] = new G4PhaseSpaceDecayChannel("eta",0.3256,3,"pi0","pi0","pi0");
00090   // eta -> pi0 + pi+ + pi-
00091   mode[2] = new G4PhaseSpaceDecayChannel("eta",0.226,3,"pi0","pi+","pi-");
00092   // eta -> gamma + pi+ + pi-
00093   mode[3] = new G4PhaseSpaceDecayChannel("eta",0.0468,3,"gamma","pi+","pi-");
00094 
00095   for (G4int index=0; index <4; index++ ) table->Insert(mode[index]);
00096   delete [] mode;
00097 
00098    anInstance->SetDecayTable(table);
00099   }
00100   theInstance = reinterpret_cast<G4Eta*>(anInstance);
00101   return theInstance;
00102 }
00103 
00104 G4Eta*  G4Eta::EtaDefinition()
00105 {
00106   return Definition();
00107 }
00108 
00109 G4Eta*  G4Eta::Eta()
00110 {
00111   return Definition();
00112 }
00113 

Generated on Mon May 27 17:48:12 2013 for Geant4 by  doxygen 1.4.7