45 :
code(0),theParticleType(
""),
47 exotic(0),radial(0),multiplet(0),
48 quark1(0),quark2(0),quark3(0),spin(0)
113 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
114 G4cout <<
" meson has three quark ";
122 if (
std::abs(tempPDGcode)%10000 == 3122) {
125 }
else if (
std::abs(tempPDGcode)%10000 == 3124) {
128 }
else if (
std::abs(tempPDGcode)%10000 == 3126) {
131 }
else if (
std::abs(tempPDGcode)%10000 == 3128) {
134 }
else if (
std::abs(tempPDGcode)%10000 == 4122) {
137 }
else if (
std::abs(tempPDGcode)%10000 == 5122) {
140 }
else if (
std::abs(tempPDGcode)%10000 == 4132) {
143 }
else if (
std::abs(tempPDGcode)%10000 == 4232) {
146 }
else if (
std::abs(tempPDGcode)%10000 == 5132) {
149 }
else if (
std::abs(tempPDGcode)%10000 == 5232) {
152 }
else if (
std::abs(tempPDGcode)%10000 == 2122) {
155 }
else if (
std::abs(tempPDGcode)%10000 == 1212) {
158 }
else if (
std::abs(tempPDGcode)%10000 == 2126) {
161 }
else if (
std::abs(tempPDGcode)%10000 == 1216) {
164 }
else if (
std::abs(tempPDGcode)%10000 == 2128) {
167 }
else if (
std::abs(tempPDGcode)%10000 == 1218) {
170 }
else if (
std::abs(tempPDGcode)%10000 == 2124) {
173 }
else if (
std::abs(tempPDGcode)%10000 == 1214) {
182 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
183 G4cout <<
" illegal code for baryon ";
192 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
193 G4cout <<
" ??? unknown quark ";
202 if (tempPDGcode >0) {
221 if (tempPDGcode == 310)
spin = 0;
222 if (tempPDGcode == 130) {
232 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
233 G4cout <<
" meson has only quark and anti-quark pair";
242 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
243 G4cout <<
" illegal code for meson ";
254 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
255 G4cout <<
" ??? unknown quark ";
266 if (tempPDGcode >0) {
275 if (tempPDGcode >0) {
303 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
304 G4cout <<
" ??? unknown quark ";
332 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
333 G4cout <<
" ??? unknown quark ";
359 totalCharge += 2./3.*
eplus*theQuarkContent[flavor+1];
360 totalCharge += (-2./3.)*
eplus*theAntiQuarkContent[flavor+1];
363 if (std::fabs(totalCharge-thePDGCharge)>0.1*
eplus) {
366 G4cout <<
" G4PDGCodeChecker::CheckCharge : ";
367 G4cout <<
" illegal electric charge " << thePDGCharge/
eplus;
380 if (pcode < 1000000000) {
387 pcode -= 10000000*
LL;
394 if (A < 2 || Z > A-LL || LL>A ) {
397 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
398 G4cout <<
" ??? Illegal PDG encoding for nucleus ";
405 G4int n_up = 2*Z + (A-Z-
LL) + LL;
406 G4int n_down = Z + 2*(A-Z-
LL) + LL;