]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliPDG.cxx
Updated task list
[u/mrichter/AliRoot.git] / STEER / AliPDG.cxx
index 2b63100459501a49e19c08093bbe22b27fb21baa..2395e4aaad0eb8d3e533ebcc29ca379c5544cb08 100644 (file)
  * 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);
 }