+ fPDGCode[fNPDGCodes++]=kspe+51; // 51 = FeedbackPhoton
+
+ Gspart(52, "Lambda_c+", 4, 2.2849, +1., 2.06e-13);
+ fPDGCode[fNPDGCodes++]=4122; //52 = Lambda_c+
+
+ Gspart(53, "Lambda_c-", 4, 2.2849, -1., 2.06e-13);
+ fPDGCode[fNPDGCodes++]=-4122; //53 = Lambda_c-
+
+ Gspart(54, "D_s+", 4, 1.9685, +1., 4.67e-13);
+ fPDGCode[fNPDGCodes++]=431; //54 = D_s+
+
+ Gspart(55, "D_s-", 4, 1.9685, -1., 4.67e-13);
+ fPDGCode[fNPDGCodes++]=-431; //55 = D_s-
+
+ Gspart(56, "Tau+", 5, 1.77705, +1., 2.9e-13);
+ fPDGCode[fNPDGCodes++]=15; //56 = Tau+
+
+ Gspart(57, "Tau-", 5, 1.77705, -1., 2.9e-13);
+ fPDGCode[fNPDGCodes++]=-15; //57 = Tau-
+
+ Gspart(58, "B0", 3, 5.2792, +0., 1.56e-12);
+ fPDGCode[fNPDGCodes++]=511; //58 = B0
+
+ Gspart(59, "B0 bar", 3, 5.2792, -0., 1.56e-12);
+ fPDGCode[fNPDGCodes++]=-511; //58 = B0bar
+
+ Gspart(60, "B+", 4, 5.2789, +1., 1.65e-12);
+ fPDGCode[fNPDGCodes++]=521; //60 = B+
+
+ Gspart(61, "B-", 4, 5.2789, -1., 1.65e-12);
+ fPDGCode[fNPDGCodes++]=-521; //61 = B-
+
+ Gspart(62, "Bs", 3, 5.3693, +0., 1.54e-12);
+ fPDGCode[fNPDGCodes++]=521; //62 = B_s
+
+ Gspart(63, "Bs bar", 3, 5.3693, -0., 1.54e-12);
+ fPDGCode[fNPDGCodes++]=-521; //63 = B_s bar
+
+ Gspart(64, "Lambda_b", 3, 5.624, +0., 1.24e-12);
+ fPDGCode[fNPDGCodes++]=5122; //64 = Lambda_b
+
+ Gspart(65, "Lambda_b bar", 3, 5.624, -0., 1.24e-12);
+ fPDGCode[fNPDGCodes++]=-5122; //65 = Lambda_b bar
+
+ Gspart(66, "J/Psi", 3.09688, 3, 0., 0.);
+ fPDGCode[fNPDGCodes++]=443; // 66 = J/Psi
+
+ Gspart(67, "Psi Prime", 3, 3.686, 0., 0.);
+ fPDGCode[fNPDGCodes++]=20443; // 67 = Psi prime
+
+ Gspart(68, "Upsilon(1S)", 9.46037, 3, 0., 0.);
+ fPDGCode[fNPDGCodes++]=553; // 68 = Upsilon(1S)
+
+ Gspart(69, "Upsilon(2S)", 10.0233, 3, 0., 0.);
+ fPDGCode[fNPDGCodes++]=20553; // 69 = Upsilon(2S)
+
+ Gspart(70, "Upsilon(3S)", 10.3553, 3, 0., 0.);
+ fPDGCode[fNPDGCodes++]=30553; // 70 = Upsilon(3S)