86 antinucleon = particle2;
90 antinucleon = particle1;
97 const std::vector<G4double> BFMM161 = {-6.434, 1.351, -5.185, 7.754, -1.692, 1.604};
99 const std::vector<G4double> BFMM169 = {3.696, -5.356, -0.053, 1.941, -0.432, 1.624};
101 const std::vector<G4double> BFMM201 = {-1.070, -0.636, -0.009, 2.335, -0.499, 1.624};
103 const std::vector<G4double> BFMM197 = {1.857, -21.213, -3.448, 0.827, -0.390, 1.616};
123 if(nucleon->getType()==
Proton){
135 antinucleon->setType(
Proton);
149 antinucleon->setType(
Proton);
192 antinucleon->setType(
Proton);
233 antinucleon->setType(
Proton);
250 antinucleon->setType(
Proton);
264 antinucleon->setType(
Proton);
320 antinucleon->setType(
Proton);
348 antinucleon->setType(
Proton);
355 list.push_back(nucleon);
356 list.push_back(antinucleon);
366 if (rdm < 1.0 / 3.0) {
367 pion1 =
new Particle(Pion1, zero, rcol);
368 pion2 =
new Particle(Pion2, zero, rcol);
369 pion3 =
new Particle(Pion3, zero, rcol);
370 }
else if (rdm < 2.0 / 3.0) {
371 pion1 =
new Particle(Pion1, zero, rcol);
372 pion2 =
new Particle(Pion3, zero, rcol);
373 pion3 =
new Particle(Pion2, zero, rcol);
375 pion1 =
new Particle(Pion2, zero, rcol);
376 pion2 =
new Particle(Pion1, zero, rcol);
377 pion3 =
new Particle(Pion3, zero, rcol);
380 list.push_back(pion1);
381 list.push_back(pion2);
382 list.push_back(pion3);