Bug corrected.
[u/mrichter/AliRoot.git] / STEER / AliPDG.cxx
CommitLineData
56050439 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
acd84897 16/* $Id$ */
56050439 17
18// Class to encapsulate the ALICE updates to TDatabasePDG.h
19// Can be used by TGeant3 and TGeant4
20// It contains also the constants for the PDG particle IDs.
21// Should evolve towards dynamical loading from external data base.
22// Comments to: andreas.morsch@cern.ch
23
24#include "AliPDG.h"
25#include "TDatabasePDG.h"
26
27ClassImp(AliPDG)
28
f3069e25 29
30
31void AliPDG::AddParticlesToPdgDataBase()
56050439 32{
33
34//
35// Add particles to the PDG data base
f3069e25 36//
37
38 static Bool_t bAdded = kFALSE;
39 // Check if already called
40 if(bAdded)return;
41 bAdded = true;
42
43 TDatabasePDG *pdgDB = TDatabasePDG::Instance();
44 const Int_t kspe=50000000;
56050439 45
f3069e25 46 // PDG nuclear states are 10-digit numbers
47 // 10LZZZAAAI e.g. deuteron is
48 // 1000010020
49 const Int_t kion=1000000000;
50
b5076441 51/*
f3069e25 52 const Double_t kAu2Gev=0.9314943228;
b5076441 53*/
54
f3069e25 55 const Double_t khSlash = 1.0545726663e-27;
56 const Double_t kErg2Gev = 1/1.6021773349e-3;
57 const Double_t khShGev = khSlash*kErg2Gev;
58 const Double_t kYear2Sec = 3600*24*365.25;
56050439 59
56050439 60//
61// Bottom mesons
62// mass and life-time from PDG
56050439 63//
5c618045 64 pdgDB->AddParticle("Upsilon(3S)","Upsilon(3S)",10.3552,kTRUE,
65 0,1,"Bottonium",200553);
66
f51912a2 67// QCD diffractive states
68 pdgDB->AddParticle("rho_diff0","rho_diff0",0,kTRUE,
69 0,0,"QCD diffr. state",9900110);
70 pdgDB->AddParticle("pi_diffr+","pi_diffr+",0,kTRUE,
71 0,1,"QCD diffr. state",9900210);
72 pdgDB->AddParticle("omega_di","omega_di",0,kTRUE,
73 0,0,"QCD diffr. state",9900220);
74 pdgDB->AddParticle("phi_diff","phi_diff",0,kTRUE,
75 0,0,"QCD diffr. state",9900330);
76 pdgDB->AddParticle("J/psi_di","J/psi_di",0,kTRUE,
77 0,0,"QCD diffr. state",9900440);
78 pdgDB->AddParticle("n_diffr0","n_diffr0",0,kTRUE,
79 0,0,"QCD diffr. state",9902110);
80 pdgDB->AddParticle("p_diffr+","p_diffr+",0,kTRUE,
81 0,1,"QCD diffr. state",9902210);
aac346c3 82
169745d2 83// From Herwig
84 pdgDB->AddParticle("PSID ", " ", 3.7699, kFALSE, 0.0, 0, "meson", 30443);
85
65327f03 86 pdgDB->AddParticle("A_00 ", " ", 0.9960, kFALSE, 0.0, 0, "meson", 9000111);
b5076441 87 pdgDB->AddParticle("A_0+ ", " ", 0.9960, kFALSE, 0.0,+3, "meson", 9000211);
88 pdgDB->AddParticle("A_0- ", " ", 0.9960, kFALSE, 0.0,-3, "meson", -9000211);
65327f03 89
90 pdgDB->AddParticle("F0P0 ", " ", 0.9960, kFALSE, 0.0, 0, "meson", 9010221);
169745d2 91
b5076441 92 pdgDB->AddParticle("KDL_2+ ", " ", 1.773, kFALSE, 0.0,+3, "meson", 10325);
93 pdgDB->AddParticle("KDL_2- ", " ", 1.773, kFALSE, 0.0,-3, "meson", -10325);
65327f03 94
169745d2 95 pdgDB->AddParticle("KDL_20 ", " ", 1.773, kFALSE, 0.0, 0, "meson", 10315);
65327f03 96 pdgDB->AddParticle("KDL_2BR0", " ", 1.773, kFALSE, 0.0, 0, "meson", -10315);
97
b5076441 98 pdgDB->AddParticle("PI_2+ ", " ", 1.670, kFALSE, 0.0,+3, "meson", 10215);
99 pdgDB->AddParticle("PI_2- ", " ", 1.670, kFALSE, 0.0,-3, "meson", -10215);
169745d2 100 pdgDB->AddParticle("PI_20 ", " ", 1.670, kFALSE, 0.0, 0, "meson", 10115);
101
102
b5076441 103 pdgDB->AddParticle("KD*+ ", " ", 1.717, kFALSE, 0.0,+3, "meson", 30323);
104 pdgDB->AddParticle("KD*- ", " ", 1.717, kFALSE, 0.0,-3, "meson", -30323);
65327f03 105
c0797b51 106 pdgDB->AddParticle("KD*0 ", " ", 1.717, kFALSE, 0.0, 0, "meson", 30313);
65327f03 107 pdgDB->AddParticle("KDBR*0 ", " ", 1.717, kFALSE, 0.0, 0, "meson", -30313);
108
b5076441 109 pdgDB->AddParticle("RHOD+ ", " ", 1.700, kFALSE, 0.0,+3, "meson", 30213);
110 pdgDB->AddParticle("RHOD- ", " ", 1.700, kFALSE, 0.0,-3, "meson", -30213);
169745d2 111 pdgDB->AddParticle("RHOD0 ", " ", 1.700, kFALSE, 0.0, 0, "meson", 30113);
112
113 pdgDB->AddParticle("ETA_2(L)", " ", 1.632, kFALSE, 0.0, 0, "meson", 10225);
114 pdgDB->AddParticle("ETA_2(H)", " ", 1.854, kFALSE, 0.0, 0, "meson", 10335);
115 pdgDB->AddParticle("OMEGA(H)", " ", 1.649, kFALSE, 0.0, 0, "meson", 30223);
116
117
b5076441 118 pdgDB->AddParticle("KDH_2+ ", " ", 1.816, kFALSE, 0.0,+3, "meson", 20325);
119 pdgDB->AddParticle("KDH_2- ", " ", 1.816, kFALSE, 0.0,-3, "meson", -20325);
65327f03 120
c0797b51 121 pdgDB->AddParticle("KDH_20 ", " ", 1.816, kFALSE, 0.0, 0, "meson", 20315);
65327f03 122 pdgDB->AddParticle("KDH_2BR0", " ", 1.816, kFALSE, 0.0, 0, "meson", -20315);
169745d2 123
124
b5076441 125 pdgDB->AddParticle("KD_3+ ", " ", 1.773, kFALSE, 0.0,+3, "meson", 327);
126 pdgDB->AddParticle("KD_3- ", " ", 1.773, kFALSE, 0.0,-3, "meson", -327);
65327f03 127
169745d2 128 pdgDB->AddParticle("KD_30 ", " ", 1.773, kFALSE, 0.0, 0, "meson", 317);
65327f03 129 pdgDB->AddParticle("KD_3BR0 ", " ", 1.773, kFALSE, 0.0, 0, "meson", -317);
169745d2 130
b5076441 131 pdgDB->AddParticle("RHO_3+ ", " ", 1.691, kFALSE, 0.0,+3, "meson", 217);
132 pdgDB->AddParticle("RHO_3- ", " ", 1.691, kFALSE, 0.0,-3, "meson", -217);
169745d2 133 pdgDB->AddParticle("RHO_30 ", " ", 1.691, kFALSE, 0.0, 0, "meson", 117);
134 pdgDB->AddParticle("OMEGA_3 ", " ", 1.667, kFALSE, 0.0, 0, "meson", 227);
135 pdgDB->AddParticle("PHI_3 ", " ", 1.854, kFALSE, 0.0, 0, "meson", 337);
136
137 pdgDB->AddParticle("CHI2P_B0", " ", 10.232, kFALSE, 0.0, 0, "meson", 110551);
138 pdgDB->AddParticle("CHI2P_B1", " ", 10.255, kFALSE, 0.0, 0, "meson", 120553);
139 pdgDB->AddParticle("CHI2P_B2", " ", 10.269, kFALSE, 0.0, 0, "meson", 100555);
140 pdgDB->AddParticle("UPSLON4S", " ", 10.580, kFALSE, 0.0, 0, "meson", 300553);
141
aac346c3 142
f3069e25 143 // IONS
144 //
145 // Done by default now from Pythia6 table
146 // Needed for other generators
147 // So check if already defined
148
149
150 Int_t ionCode = kion+10020;
151 if(!pdgDB->GetParticle(ionCode)){
b5076441 152 pdgDB->AddParticle("Deuteron","Deuteron", 1.875613, kTRUE,
153 0,3,"Ion",ionCode);
f3069e25 154 }
b5076441 155 pdgDB->AddAntiParticle("AntiDeuteron", - ionCode);
5b9463ce 156
f3069e25 157 ionCode = kion+10030;
158 if(!pdgDB->GetParticle(ionCode)){
b5076441 159 pdgDB->AddParticle("Triton","Triton", 2.80925, kFALSE,
160 khShGev/(12.33*kYear2Sec),3,"Ion",ionCode);
f3069e25 161 }
b5076441 162 pdgDB->AddAntiParticle("AntiTriton", - ionCode);
f3069e25 163
164 ionCode = kion+20030;
165 if(!pdgDB->GetParticle(ionCode)){
b5076441 166 pdgDB->AddParticle("HE3","HE3", 2.80923,kFALSE,
167 0,6,"Ion",ionCode);
f3069e25 168 }
b5076441 169 pdgDB->AddAntiParticle("AntiHE3", - ionCode);
f3069e25 170
171 ionCode = kion+20040;
172 if(!pdgDB->GetParticle(ionCode)){
b5076441 173 pdgDB->AddParticle("Alpha","Alpha", 3.727417, kTRUE,
174 khShGev/(12.33*kYear2Sec), 6, "Ion", ionCode);
f3069e25 175 }
b5076441 176 pdgDB->AddAntiParticle("AntiAlpha", - ionCode);
f3069e25 177
41fc1188 178// Special particles
56050439 179//
180 pdgDB->AddParticle("Cherenkov","Cherenkov",0,kFALSE,
181 0,0,"Special",kspe+50);
182 pdgDB->AddParticle("FeedbackPhoton","FeedbackPhoton",0,kFALSE,
183 0,0,"Special",kspe+51);
13313168 184 pdgDB->AddParticle("Lambda1520","Lambda1520",1.5195,kFALSE,
185 0.0156,0,"Resonance",3124);
186 pdgDB->AddAntiParticle("Lambda1520bar",-3124);
56050439 187}
41fc1188 188
189