X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STEER%2FAliPDG.cxx;h=2395e4aaad0eb8d3e533ebcc29ca379c5544cb08;hb=f3be8296f226d7770da6be055300e20b3a4fddcc;hp=2b63100459501a49e19c08093bbe22b27fb21baa;hpb=5c618045bcb3d199a87008bc83e44ef7817b9c9e;p=u%2Fmrichter%2FAliRoot.git diff --git a/STEER/AliPDG.cxx b/STEER/AliPDG.cxx index 2b631004595..2395e4aaad0 100644 --- a/STEER/AliPDG.cxx +++ b/STEER/AliPDG.cxx @@ -13,18 +13,7 @@ * provided "as is" without express or implied warranty. * **************************************************************************/ -/* -$Log$ -Revision 1.3 2001/03/16 09:48:35 morsch -Exclude pdg particle definitions. Done by root now by default from Pythia6 table. - -Revision 1.2 2001/01/31 14:32:42 morsch -Some B mesons added - -Revision 1.1 2000/12/21 16:48:39 morsch -AliPDG class, first commit. - -*/ +/* $Id$ */ // Class to encapsulate the ALICE updates to TDatabasePDG.h // Can be used by TGeant3 and TGeant4 @@ -37,22 +26,35 @@ AliPDG class, first commit. ClassImp(AliPDG) - void AliPDG::AddParticlesToPdgDataBase() + + +void AliPDG::AddParticlesToPdgDataBase() { // // Add particles to the PDG data base +// + + static Bool_t bAdded = kFALSE; + // Check if already called + if(bAdded)return; + bAdded = true; + + TDatabasePDG *pdgDB = TDatabasePDG::Instance(); + const Int_t kspe=50000000; - TDatabasePDG *pdgDB = TDatabasePDG::Instance(); + // PDG nuclear states are 10-digit numbers + // 10LZZZAAAI e.g. deuteron is + // 1000010020 + const Int_t kion=1000000000; + + const Double_t kAu2Gev=0.9314943228; + const Double_t khSlash = 1.0545726663e-27; + const Double_t kErg2Gev = 1/1.6021773349e-3; + const Double_t khShGev = khSlash*kErg2Gev; + const Double_t kYear2Sec = 3600*24*365.25; - const Int_t kion=10000000; - const Int_t kspe=50000000; - const Double_t kAu2Gev=0.9314943228; - const Double_t khSlash = 1.0545726663e-27; - const Double_t kErg2Gev = 1/1.6021773349e-3; - const Double_t khShGev = khSlash*kErg2Gev; - const Double_t kYear2Sec = 3600*24*365.25; // // Bottom mesons // mass and life-time from PDG @@ -60,27 +62,125 @@ ClassImp(AliPDG) pdgDB->AddParticle("Upsilon(3S)","Upsilon(3S)",10.3552,kTRUE, 0,1,"Bottonium",200553); -// Done by default now from Pythia6 table! -// -// -// Ions -// +// QCD diffractive states + pdgDB->AddParticle("rho_diff0","rho_diff0",0,kTRUE, + 0,0,"QCD diffr. state",9900110); + pdgDB->AddParticle("pi_diffr+","pi_diffr+",0,kTRUE, + 0,1,"QCD diffr. state",9900210); + pdgDB->AddParticle("omega_di","omega_di",0,kTRUE, + 0,0,"QCD diffr. state",9900220); + pdgDB->AddParticle("phi_diff","phi_diff",0,kTRUE, + 0,0,"QCD diffr. state",9900330); + pdgDB->AddParticle("J/psi_di","J/psi_di",0,kTRUE, + 0,0,"QCD diffr. state",9900440); + pdgDB->AddParticle("n_diffr0","n_diffr0",0,kTRUE, + 0,0,"QCD diffr. state",9902110); + pdgDB->AddParticle("p_diffr+","p_diffr+",0,kTRUE, + 0,1,"QCD diffr. state",9902210); + +// From Herwig + pdgDB->AddParticle("PSID ", " ", 3.7699, kFALSE, 0.0, 0, "meson", 30443); + + pdgDB->AddParticle("A_00 ", " ", 0.9960, kFALSE, 0.0, 0, "meson", 9000111); + pdgDB->AddParticle("A_0+ ", " ", 0.9960, kFALSE, 0.0,+1, "meson", 9000211); + pdgDB->AddParticle("A_0- ", " ", 0.9960, kFALSE, 0.0,-1, "meson", -9000211); + + pdgDB->AddParticle("F0P0 ", " ", 0.9960, kFALSE, 0.0, 0, "meson", 9010221); + + pdgDB->AddParticle("KDL_2+ ", " ", 1.773, kFALSE, 0.0,+1, "meson", 10325); + pdgDB->AddParticle("KDL_2- ", " ", 1.773, kFALSE, 0.0,-1, "meson", -10325); + + pdgDB->AddParticle("KDL_20 ", " ", 1.773, kFALSE, 0.0, 0, "meson", 10315); + pdgDB->AddParticle("KDL_2BR0", " ", 1.773, kFALSE, 0.0, 0, "meson", -10315); + + pdgDB->AddParticle("PI_2+ ", " ", 1.670, kFALSE, 0.0,+1, "meson", 10215); + pdgDB->AddParticle("PI_2- ", " ", 1.670, kFALSE, 0.0,-1, "meson", -10215); + pdgDB->AddParticle("PI_20 ", " ", 1.670, kFALSE, 0.0, 0, "meson", 10115); + + + pdgDB->AddParticle("KD*+ ", " ", 1.717, kFALSE, 0.0,+1, "meson", 30323); + pdgDB->AddParticle("KD*- ", " ", 1.717, kFALSE, 0.0,-1, "meson", -30323); + + pdgDB->AddParticle("KD*0 ", " ", 1.717, kFALSE, 0.0, 0, "meson", 30313); + pdgDB->AddParticle("KDBR*0 ", " ", 1.717, kFALSE, 0.0, 0, "meson", -30313); + + pdgDB->AddParticle("RHOD+ ", " ", 1.700, kFALSE, 0.0,+1, "meson", 30213); + pdgDB->AddParticle("RHOD- ", " ", 1.700, kFALSE, 0.0,-1, "meson", -30213); + pdgDB->AddParticle("RHOD0 ", " ", 1.700, kFALSE, 0.0, 0, "meson", 30113); + + pdgDB->AddParticle("ETA_2(L)", " ", 1.632, kFALSE, 0.0, 0, "meson", 10225); + pdgDB->AddParticle("ETA_2(H)", " ", 1.854, kFALSE, 0.0, 0, "meson", 10335); + pdgDB->AddParticle("OMEGA(H)", " ", 1.649, kFALSE, 0.0, 0, "meson", 30223); + + + pdgDB->AddParticle("KDH_2+ ", " ", 1.816, kFALSE, 0.0,+1, "meson", 20325); + pdgDB->AddParticle("KDH_2- ", " ", 1.816, kFALSE, 0.0,-1, "meson", -20325); + + pdgDB->AddParticle("KDH_20 ", " ", 1.816, kFALSE, 0.0, 0, "meson", 20315); + pdgDB->AddParticle("KDH_2BR0", " ", 1.816, kFALSE, 0.0, 0, "meson", -20315); + + + pdgDB->AddParticle("KD_3+ ", " ", 1.773, kFALSE, 0.0,+1, "meson", 327); + pdgDB->AddParticle("KD_3- ", " ", 1.773, kFALSE, 0.0,-1, "meson", -327); + + pdgDB->AddParticle("KD_30 ", " ", 1.773, kFALSE, 0.0, 0, "meson", 317); + pdgDB->AddParticle("KD_3BR0 ", " ", 1.773, kFALSE, 0.0, 0, "meson", -317); + + pdgDB->AddParticle("RHO_3+ ", " ", 1.691, kFALSE, 0.0,+1, "meson", 217); + pdgDB->AddParticle("RHO_3- ", " ", 1.691, kFALSE, 0.0,-1, "meson", -217); + pdgDB->AddParticle("RHO_30 ", " ", 1.691, kFALSE, 0.0, 0, "meson", 117); + pdgDB->AddParticle("OMEGA_3 ", " ", 1.667, kFALSE, 0.0, 0, "meson", 227); + pdgDB->AddParticle("PHI_3 ", " ", 1.854, kFALSE, 0.0, 0, "meson", 337); + + pdgDB->AddParticle("CHI2P_B0", " ", 10.232, kFALSE, 0.0, 0, "meson", 110551); + pdgDB->AddParticle("CHI2P_B1", " ", 10.255, kFALSE, 0.0, 0, "meson", 120553); + pdgDB->AddParticle("CHI2P_B2", " ", 10.269, kFALSE, 0.0, 0, "meson", 100555); + pdgDB->AddParticle("UPSLON4S", " ", 10.580, kFALSE, 0.0, 0, "meson", 300553); + + + // IONS + // + // Done by default now from Pythia6 table + // Needed for other generators + // So check if already defined + + + Int_t ionCode = kion+10020; + if(!pdgDB->GetParticle(ionCode)){ + pdgDB->AddParticle("Deuteron","Deuteron",2*kAu2Gev+8.071e-3,kTRUE, + 0,1,"Ion",ionCode); + } + + pdgDB->AddAntiParticle("Anti Deuteron", - ionCode); + + ionCode = kion+10030; + if(!pdgDB->GetParticle(ionCode)){ + pdgDB->AddParticle("Triton","Triton",3*kAu2Gev+14.931e-3,kFALSE, + khShGev/(12.33*kYear2Sec),1,"Ion",ionCode); + } + + ionCode = kion+20030; + if(!pdgDB->GetParticle(ionCode)){ + pdgDB->AddParticle("HE3","HE3",3*kAu2Gev+14.931e-3,kFALSE, + 0,2,"Ion",ionCode); + } + + ionCode = kion+20040; + if(!pdgDB->GetParticle(ionCode)){ + pdgDB->AddParticle("Alpha","Alpha",4*kAu2Gev+2.424e-3,kTRUE, + khShGev/(12.33*kYear2Sec),2,"Ion",ionCode); + } + - pdgDB->AddParticle("Deuteron","Deuteron",2*kAu2Gev+8.071e-3,kTRUE, - 0,1,"Ion",kion+10020); - pdgDB->AddParticle("Triton","Triton",3*kAu2Gev+14.931e-3,kFALSE, - khShGev/(12.33*kYear2Sec),1,"Ion",kion+10030); - pdgDB->AddParticle("Alpha","Alpha",4*kAu2Gev+2.424e-3,kTRUE, - khShGev/(12.33*kYear2Sec),2,"Ion",kion+20040); - pdgDB->AddParticle("HE3","HE3",3*kAu2Gev+14.931e-3,kFALSE, - 0,2,"Ion",kion+20030); // Special particles // pdgDB->AddParticle("Cherenkov","Cherenkov",0,kFALSE, 0,0,"Special",kspe+50); pdgDB->AddParticle("FeedbackPhoton","FeedbackPhoton",0,kFALSE, 0,0,"Special",kspe+51); - + pdgDB->AddParticle("Lambda1520","Lambda1520",1.5195,kFALSE, + 0.0156,0,"Resonance",3124); + pdgDB->AddAntiParticle("Lambda1520bar",-3124); }