SetDecay.cxx Source File

00001 // -*- C++ -*-  $Id: SetDecay.cxx,v 1.2 2000/01/23 03:37:54 pfkeb Exp $
00002 //
00003 // This file is part of Gismo 2
00004 //
00005 
00006 #include "HepPDT.h"
00007 #include "gismo/PhaseSpace.h"
00008 #include "Sage.h"
00009 #include "gismo/Interactor.h"
00010 #include "gismo/PData.h"
00011 #ifdef _MSC_VER
00012 # pragma warning(disable: 4305)
00013 #endif
00014 
00015 HepPDT&
00016 HepPDT::setDecays()
00017 {
00018  //--------------------------------------------------------------------------
00019  //                     D e c a y s
00020  //--------------------------------------------------------------------------
00021 
00022  // Create the default phase-space decayers. Note that anti-particle decays
00023  // are defined automatically, and that the branching fractions do not have
00024  // to add up to one
00025 
00026  addDecayer(new PhaseSpace("PhaseSpace") );
00027  addDecayer(new SageDecayer("Sage"));
00028 
00029  //     parent  bf   method      ------- daughters -------------
00030  //     ------  --   ------
00031  addDecay("mu-",  1.0,  "PhaseSpace",     "e-",    "nu(e)",   "nu(mu)");
00032  addDecay("pi+",  1.0,  "PhaseSpace",  "mu+",   "nu(mu)", "");
00033 
00034  addDecay("pi0",  1.0,  "PhaseSpace",  "gamma", "gamma");
00035  addDecay("K(S)0", 1.0,  "PhaseSpace",  "pi+",   "pi-", "");
00036  addDecay("K(L)0", 1.0,  "PhaseSpace",     "pi+",   "pi-",    "pi0");
00037  addDecay("K+",   0.635,"PhaseSpace",  "mu+",   "nu(mu)");
00038  addDecay("K+",   0.212,"PhaseSpace",  "pi+",   "pi0");
00039  addDecay("rho(770)0", 1.0,  "PhaseSpace",  "pi+",   "pi-");
00040  addDecay("rho(770)+", 1.0,  "PhaseSpace",  "pi+",   "pi0");
00041  addDecay("omega(783)",1.0,  "PhaseSpace",     "pi+",   "pi0",    "pi-");
00042  addDecay("phi(1020)",  0.495,"PhaseSpace",  "K+",    "K-");
00043  addDecay("phi(1020)",  0.344,"PhaseSpace",  "K(L)0",  "K(S)0");
00044  addDecay("phi(1020)",  0.129/3.,"PhaseSpace", "rho(770)0",   "pi0");
00045  addDecay("phi(1020)",  0.129/3.,"PhaseSpace", "rho(770)+",   "pi-");
00046  addDecay("phi(1020)",  0.129/3.,"PhaseSpace", "rho(770)-",   "pi+");
00047 
00048  addDecay("D+",   0.026/2, "PhaseSpace", "K(S)0", "pi+");
00049  addDecay("D+",   0.026/2, "PhaseSpace", "K(L)0", "pi+");
00050  addDecay("D+",   0.080,   "PhaseSpace",    "K-",   "pi+",  "pi+");
00051 
00052  addDecay("D0",   0.021/2, "PhaseSpace", "K(S)0", "pi0");
00053  addDecay("D0",   0.021/2, "PhaseSpace", "K(L)0", "pi0");
00054  addDecay("D0",   0.054/2, "PhaseSpace",    "K(S)0",  "pi+",  "pi-");
00055  addDecay("D0",   0.054/2, "PhaseSpace",    "K(L)0",  "pi+",  "pi-");
00056 
00057  addDecay("D*(2010)+",  0.55,    "PhaseSpace", "D0",    "pi+");
00058  addDecay("D*(2010)+",  0.272,   "PhaseSpace", "D+",    "pi0");
00059  addDecay("D*(2010)+",  0.18,    "PhaseSpace", "D+",    "gamma");
00060 
00061  addDecay("D*(2010)0",  0.55,    "PhaseSpace", "D0",    "pi0");
00062  addDecay("D*(2010)0",  0.45,    "PhaseSpace", "D0",    "gamma");
00063 
00064  addDecay("B+",   0.0038,   "PhaseSpace",  "D~0", "pi+");
00065  addDecay("B+",   0.013,    "PhaseSpace",  "D~0", "rho(770)+");
00066  addDecay("B+",   0.005,    "Sage",     "D~0", "pi+", "pi+", "pi-");
00067  addDecay("B+",   0.010,    "PhaseSpace",  "D*(2010)0",    "rho(770)+");
00068  addDecay("B+",   0.018,    "Sage",     "D*(2010)-",    "pi+", "pi+", "pi0");
00069 
00070  addDecay("B0",   0.0032,   "PhaseSpace",  "D-", "pi+");
00071  addDecay("B0",   0.009,    "PhaseSpace",  "D-", "rho(770)+");
00072  addDecay("B0",   0.0032,   "PhaseSpace",  "D*(2010)-","pi+");
00073  addDecay("B0",   0.008,    "Sage",     "D-", "pi+", "pi+", "pi-");
00074 
00075  addDecay("Upsilon(4S)",0.5,"PhaseSpace",  "B0", "B~0");
00076  addDecay("Upsilon(4S)",0.5,"PhaseSpace",  "B+", "B-");
00077  return *this;
00078 }
00079