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;
-/*
- const Int_t kion=10000000;
+ TDatabasePDG *pdgDB = TDatabasePDG::Instance();
+ const Int_t kspe=50000000;
+
+ // 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 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
0,0,"QCD diffr. state",9902110);
pdgDB->AddParticle("p_diffr+","p_diffr+",0,kTRUE,
0,1,"QCD diffr. state",9902210);
-// Done by default now from Pythia6 table!
-//
-//
-// Ions
-//
-/*
- 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);
-*/
+
+// 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);
+ }
+
+ 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);
+ }
+
+
// 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);
}