46 :
code(0),theParticleType(
""),
48 exotic(0),radial(0),multiplet(0),
49 quark1(0),quark2(0),quark3(0),spin(0)
114 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
115 G4cout <<
" meson has three quark ";
123 if (std::abs(tempPDGcode)%10000 == 3122) {
126 }
else if (std::abs(tempPDGcode)%10000 == 3124) {
129 }
else if (std::abs(tempPDGcode)%10000 == 3126) {
132 }
else if (std::abs(tempPDGcode)%10000 == 3128) {
135 }
else if (std::abs(tempPDGcode)%10000 == 4122) {
138 }
else if (std::abs(tempPDGcode)%10000 == 5122) {
141 }
else if (std::abs(tempPDGcode)%10000 == 4132) {
144 }
else if (std::abs(tempPDGcode)%10000 == 4232) {
147 }
else if (std::abs(tempPDGcode)%10000 == 5132) {
150 }
else if (std::abs(tempPDGcode)%10000 == 5232) {
153 }
else if (std::abs(tempPDGcode)%10000 == 2122) {
156 }
else if (std::abs(tempPDGcode)%10000 == 1212) {
159 }
else if (std::abs(tempPDGcode)%10000 == 2126) {
162 }
else if (std::abs(tempPDGcode)%10000 == 1216) {
165 }
else if (std::abs(tempPDGcode)%10000 == 2128) {
168 }
else if (std::abs(tempPDGcode)%10000 == 1218) {
171 }
else if (std::abs(tempPDGcode)%10000 == 2124) {
174 }
else if (std::abs(tempPDGcode)%10000 == 1214) {
183 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
184 G4cout <<
" illegal code for baryon ";
193 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
194 G4cout <<
" ??? unknown quark ";
203 if (tempPDGcode >0) {
222 if (tempPDGcode == 310)
spin = 0;
223 if (tempPDGcode == 130) {
233 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
234 G4cout <<
" meson has only quark and anti-quark pair";
243 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
244 G4cout <<
" illegal code for meson ";
255 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
256 G4cout <<
" ??? unknown quark ";
267 if (tempPDGcode >0) {
276 if (tempPDGcode >0) {
304 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
305 G4cout <<
" ??? unknown quark ";
333 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
334 G4cout <<
" ??? unknown quark ";
360 totalCharge += 2./3.*
eplus*theQuarkContent[flavor+1];
361 totalCharge += (-2./3.)*
eplus*theAntiQuarkContent[flavor+1];
364 if (std::fabs(totalCharge-thePDGCharge)>0.1*
eplus) {
367 G4cout <<
" G4PDGCodeChecker::CheckCharge : ";
368 G4cout <<
" illegal electric charge " << thePDGCharge/
eplus;
381 if (pcode < 1000000000) {
388 pcode -= 10000000*
LL;
395 if (A < 2 || Z > A-LL || LL>A ) {
398 G4cout <<
" G4PDGCodeChecker::CheckPDGCode : ";
399 G4cout <<
" ??? Illegal PDG encoding for nucleus ";
406 G4int n_up = 2*Z + (A-Z-
LL) + LL;
407 G4int n_down = Z + 2*(A-Z-
LL) + LL;
427 G4int temp = std::abs(PDGcode);
static const G4int LL[nN]
G4int CheckPDGCode(G4int code, G4String type)
G4bool CheckCharge(G4double charge) const
double A(double temperature)
G4int theAntiQuarkContent[NumberOfQuarkFlavor]
static constexpr double eplus
void GetDigits(G4int code)
G4int theQuarkContent[NumberOfQuarkFlavor]
G4GLOB_DLL std::ostream G4cout