Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // INCL++ intra-nuclear cascade model 27 // Alain Boudard, CEA-Saclay, France 28 // Joseph Cugnon, University of Liege, Belgium 29 // Jean-Christophe David, CEA-Saclay, France 30 // Pekka Kaitaniemi, CEA-Saclay, France, and H 31 // Sylvie Leray, CEA-Saclay, France 32 // Davide Mancusi, CEA-Saclay, France 33 // 34 #define INCLXX_IN_GEANT4_MODE 1 35 36 #include "globals.hh" 37 38 #include "G4INCLNNOmegaToMultiPionsChannel.hh" 39 #include "G4INCLKinematicsUtils.hh" 40 #include "G4INCLBinaryCollisionAvatar.hh" 41 #include "G4INCLRandom.hh" 42 #include "G4INCLGlobals.hh" 43 #include "G4INCLLogger.hh" 44 #include <algorithm> 45 #include "G4INCLPhaseSpaceGenerator.hh" 46 47 namespace G4INCL { 48 49 const G4double NNOmegaToMultiPionsChannel::a 50 51 NNOmegaToMultiPionsChannel::NNOmegaToMultiPi 52 : npion(npi), 53 iso1(0), 54 iso2(0), 55 particle1(p1), 56 particle2(p2) 57 { 58 std::fill(isosp, isosp+4, 0); 59 } 60 61 NNOmegaToMultiPionsChannel::~NNOmegaToMultiP 62 63 } 64 65 void NNOmegaToMultiPionsChannel::fillFinalSt 66 // assert(npion > 0 && npion < 5); 67 68 iso1=ParticleTable::getIsospin(particle1 69 iso2=ParticleTable::getIsospin(particle2 70 71 ParticleList list; 72 list.push_back(particle1); 73 list.push_back(particle2); 74 fs->addModifiedParticle(particle1); 75 fs->addModifiedParticle(particle2); 76 77 isospinRepartition(); 78 79 const ParticleType tn1=ParticleTable::ge 80 particle1->setType(tn1); 81 const ParticleType tn2=ParticleTable::ge 82 particle2->setType(tn2); 83 const ThreeVector &rcolnucleon1 = partic 84 const ThreeVector &rcolnucleon2 = partic 85 const ThreeVector rcol = (rcolnucleon1+r 86 const ThreeVector zero; 87 for(G4int i=0; i<npion; ++i) { 88 const ParticleType pionType=ParticleTa 89 Particle *pion = new Particle(pionType 90 list.push_back(pion); 91 fs->addCreatedParticle(pion); 92 } 93 Particle *omega = new Particle(Omega,zer 94 list.push_back(omega); 95 fs->addCreatedParticle(omega); 96 97 const G4double sqrtS = KinematicsUtils:: 98 G4int biasIndex = ((Random::shoot()<0.5) 99 PhaseSpaceGenerator::generateBiased(sqrt 100 101 } 102 103 void NNOmegaToMultiPionsChannel::isospinRe 104 const G4double rjcd=Random::shoot(); 105 G4double p; 106 const G4int itot=iso1+iso2; 107 108 if (npion == 1) { 109 p=3.*rjcd; 110 if (p < 1.) pn_ppPim(); 111 else if (p < 2.) pn_pnPi0(); 112 else pn_nnPip(); 113 } 114 else if (npion == 2) { 115 if (itot == 2) { 116 p=20.*rjcd; 117 if (p >= 14.) pp_nnPipPip 118 else if (p >= 11.) pp_pnPipPi0 119 else if (p >= 7.) pp_ppPi0Pi0 120 else pp_ppPipPim( 121 } 122 else if (itot == -2) { 123 p=20.*rjcd; 124 if (p >= 14.) nn_ppPimPim 125 else if (p >= 11.) nn_pnPimPi0 126 else if (p >= 7.) nn_nnPi0Pi0 127 else nn_nnPipPim( 128 } 129 else { 130 G4double pp=Random::shoot(); 131 if (pp > 0.5) { 132 p=3.*rjcd; 133 if (p < 2.) { 134 pn_pnPipPim(); 135 } 136 else { 137 pn_pnPi0Pi0(); 138 } 139 } 140 else { 141 p=60.*rjcd; 142 if (p >= 51.) pn_nnPi 143 else if (p >= 33.) pn_pnPi 144 else if (p >= 9.) pn_pnPi 145 else pn_ppPim 146 } 147 } 148 } 149 else if (npion == 3) { 150 p=60.*rjcd; 151 if (itot == 2) { 152 if (p >= 42.) pp_nnPipPip 153 else if (p >= 39.) pp_pnPipPi0 154 else if (p >= 33.) pp_pnPipPip 155 else if (p >= 22.) pp_ppPi0Pi0 156 else pp_ppPipPimP 157 } 158 else if (itot == -2) { 159 if (p >= 42.) nn_ppPimPim 160 else if (p >= 39.) nn_pnPimPi0 161 else if (p >= 33.) nn_pnPipPim 162 else if (p >= 22.) nn_nnPi0Pi0 163 else nn_nnPipPimP 164 } 165 else { 166 if (p >= 57.) pn_nnPipPi0 167 else if (p >= 51.) pn_nnPipPip 168 else if (p >= 37.) pn_pnPi0Pi0 169 else if (p >= 9.) pn_pnPi0Pip 170 else if (p >= 6.) pn_ppPimPi0 171 else pn_ppPimPimP 172 173 } 174 } 175 else if (npion == 4) { 176 p=60.*rjcd; 177 if (itot == 2) { 178 if (p >= 48.) pp_nnPipPip 179 else if (p >= 42.) pp_nnPipPip 180 else if (p >= 36.) pp_pnPipPip 181 else if (p >= 33.) pp_pnPipPi0 182 else if (p >= 19.) pp_ppPipPip 183 else if (p >= 4.) pp_ppPipPi0 184 else pp_ppPi0Pi0P 185 } 186 else if (itot == -2) { 187 if (p >= 48.) nn_ppPipPim 188 else if (p >= 42.) nn_ppPi0Pi0 189 else if (p >= 36.) nn_pnPipPi0 190 else if (p >= 33.) nn_pnPi0Pi0 191 else if (p >= 19.) nn_nnPipPip 192 else if (p >= 4.) nn_nnPipPi0 193 else nn_nnPi0Pi0P 194 } 195 else { 196 G4double pp=Random::shoot(); 197 if (pp > 0.5) { 198 p=9.*rjcd; 199 if (p < 1.) pn_pnPi0P 200 else if (p < 5.) pn_pnPipP 201 else pn_pnPipPi 202 } 203 else { 204 if (p < 3.) pn_ppPi0 205 else if (p < 9.) pn_ppPip 206 else if (p < 15.) pn_pnPi0 207 else if (p < 35.) pn_pnPip 208 else if (p < 51.) pn_pnPip 209 else if (p < 54.) pn_nnPip 210 else pn_nnPipP 211 } 212 } 213 } 214 215 std::shuffle(isosp,isosp+npion,Random: 216 inter2Part(0.5); 217 } 218 219 220 void NNOmegaToMultiPionsChannel::pn_ppPim( 221 isosp[0]=-2; 222 iso1=1; 223 iso2=1; 224 } 225 void NNOmegaToMultiPionsChannel::pn_pnPi0( 226 isosp[0]=0; 227 } 228 void NNOmegaToMultiPionsChannel::pn_nnPip( 229 isosp[0]=2; 230 iso1=-1; 231 iso2=-1; 232 } 233 void NNOmegaToMultiPionsChannel::pp_nnPipP 234 isosp[0]=2; 235 isosp[1]=2; 236 iso1=-1; 237 iso2=-1; 238 } 239 void NNOmegaToMultiPionsChannel::nn_ppPimP 240 isosp[0]=-2; 241 isosp[1]=-2; 242 iso1=1; 243 iso2=1; 244 } 245 void NNOmegaToMultiPionsChannel::pn_pnPipP 246 isosp[0]=2; 247 isosp[1]=-2; 248 } 249 void NNOmegaToMultiPionsChannel::pn_pnPi0P 250 isosp[0]=0; 251 isosp[1]=0; 252 } 253 void NNOmegaToMultiPionsChannel::pp_ppPipP 254 isosp[0]=2; 255 isosp[1]=-2; 256 } 257 void NNOmegaToMultiPionsChannel::nn_nnPipP 258 isosp[0]=2; 259 isosp[1]=-2; 260 } 261 void NNOmegaToMultiPionsChannel::pp_ppPi0P 262 isosp[0]=0; 263 isosp[1]=0; 264 } 265 void NNOmegaToMultiPionsChannel::nn_nnPi0P 266 isosp[0]=0; 267 isosp[1]=0; 268 } 269 void NNOmegaToMultiPionsChannel::pp_pnPipP 270 isosp[0]=2; 271 isosp[1]=0; 272 iso1=1; 273 iso2=-1; 274 } 275 void NNOmegaToMultiPionsChannel::pn_ppPimP 276 isosp[0]=-2; 277 isosp[1]=0; 278 iso1=1; 279 iso2=1; 280 } 281 void NNOmegaToMultiPionsChannel::pn_nnPipP 282 isosp[0]=2; 283 isosp[1]=0; 284 iso1=-1; 285 iso2=-1; 286 } 287 void NNOmegaToMultiPionsChannel::nn_pnPimP 288 isosp[0]=-2; 289 isosp[1]=0; 290 iso1=1; 291 iso2=-1; 292 } 293 void NNOmegaToMultiPionsChannel::pp_pnPipP 294 isosp[0]=2; 295 isosp[1]=0; 296 isosp[2]=0; 297 iso1=1; 298 iso2=-1; 299 } 300 void NNOmegaToMultiPionsChannel::nn_pnPimP 301 isosp[0]=-2; 302 isosp[1]=0; 303 isosp[2]=0; 304 iso1=1; 305 iso2=-1; 306 } 307 void NNOmegaToMultiPionsChannel::pn_nnPipP 308 isosp[0]=2; 309 isosp[1]=0; 310 isosp[2]=0; 311 iso1=-1; 312 iso2=-1; 313 } 314 void NNOmegaToMultiPionsChannel::pp_ppPipP 315 isosp[0]=2; 316 isosp[1]=-2; 317 isosp[2]=0; 318 } 319 void NNOmegaToMultiPionsChannel::nn_nnPipP 320 isosp[0]=2; 321 isosp[1]=-2; 322 isosp[2]=0; 323 } 324 void NNOmegaToMultiPionsChannel::pp_ppPi0P 325 isosp[0]=0; 326 isosp[1]=0; 327 isosp[2]=0; 328 } 329 void NNOmegaToMultiPionsChannel::nn_nnPi0P 330 isosp[0]=0; 331 isosp[1]=0; 332 isosp[2]=0; 333 } 334 void NNOmegaToMultiPionsChannel::pp_pnPipP 335 isosp[0]=2; 336 isosp[1]=2; 337 isosp[2]=-2; 338 iso1=1; 339 iso2=-1; 340 } 341 void NNOmegaToMultiPionsChannel::pp_nnPipP 342 isosp[0]=2; 343 isosp[1]=2; 344 isosp[2]=0; 345 iso1=-1; 346 iso2=-1; 347 } 348 void NNOmegaToMultiPionsChannel::pn_ppPimP 349 isosp[0]=-2; 350 isosp[1]=0; 351 isosp[2]=0; 352 iso1=1; 353 iso2=1; 354 } 355 void NNOmegaToMultiPionsChannel::pn_ppPimP 356 isosp[0]=-2; 357 isosp[1]=-2; 358 isosp[2]=2; 359 iso1=1; 360 iso2=1; 361 } 362 void NNOmegaToMultiPionsChannel::pn_pnPi0P 363 isosp[0]=0; 364 isosp[1]=2; 365 isosp[2]=-2; 366 } 367 void NNOmegaToMultiPionsChannel::pn_pnPi0P 368 isosp[0]=0; 369 isosp[1]=0; 370 isosp[2]=0; 371 } 372 void NNOmegaToMultiPionsChannel::pn_nnPipP 373 isosp[0]=2; 374 isosp[1]=2; 375 isosp[2]=-2; 376 iso1=-1; 377 iso2=-1; 378 } 379 void NNOmegaToMultiPionsChannel::nn_pnPipP 380 isosp[0]=2; 381 isosp[1]=-2; 382 isosp[2]=-2; 383 iso1=1; 384 iso2=-1; 385 } 386 void NNOmegaToMultiPionsChannel::nn_ppPimP 387 isosp[0]=-2; 388 isosp[1]=-2; 389 isosp[2]=0; 390 iso1=1; 391 iso2=1; 392 } 393 void NNOmegaToMultiPionsChannel::pp_nnPipP 394 isosp[0]=2; 395 isosp[1]=2; 396 isosp[2]=0; 397 isosp[3]=0; 398 iso1=-1; 399 iso2=-1; 400 } 401 void NNOmegaToMultiPionsChannel::pp_nnPipP 402 isosp[0]=2; 403 isosp[1]=2; 404 isosp[2]=2; 405 isosp[3]=-2; 406 iso1=-1; 407 iso2=-1; 408 } 409 void NNOmegaToMultiPionsChannel::nn_ppPi0P 410 isosp[0]=0; 411 isosp[1]=0; 412 isosp[2]=-2; 413 isosp[3]=-2; 414 iso1=1; 415 iso2=1; 416 } 417 void NNOmegaToMultiPionsChannel::nn_ppPipP 418 isosp[0]=2; 419 isosp[1]=-2; 420 isosp[2]=-2; 421 isosp[3]=-2; 422 iso1=1; 423 iso2=1; 424 } 425 void NNOmegaToMultiPionsChannel::pp_ppPi0P 426 isosp[0]=0; 427 isosp[1]=0; 428 isosp[2]=0; 429 isosp[3]=0; 430 } 431 void NNOmegaToMultiPionsChannel::nn_nnPi0P 432 isosp[0]=0; 433 isosp[1]=0; 434 isosp[2]=0; 435 isosp[3]=0; 436 } 437 void NNOmegaToMultiPionsChannel::pn_pnPi0P 438 isosp[0]=0; 439 isosp[1]=0; 440 isosp[2]=0; 441 isosp[3]=0; 442 } 443 void NNOmegaToMultiPionsChannel::pp_ppPipP 444 isosp[0]=2; 445 isosp[1]=0; 446 isosp[2]=0; 447 isosp[3]=-2; 448 } 449 void NNOmegaToMultiPionsChannel::nn_nnPipP 450 isosp[0]=2; 451 isosp[1]=0; 452 isosp[2]=0; 453 isosp[3]=-2; 454 } 455 void NNOmegaToMultiPionsChannel::pn_pnPipP 456 isosp[0]=2; 457 isosp[1]=0; 458 isosp[2]=0; 459 isosp[3]=-2; 460 } 461 void NNOmegaToMultiPionsChannel::pp_ppPipP 462 isosp[0]=2; 463 isosp[1]=2; 464 isosp[2]=-2; 465 isosp[3]=-2; 466 } 467 void NNOmegaToMultiPionsChannel::nn_nnPipP 468 isosp[0]=2; 469 isosp[1]=2; 470 isosp[2]=-2; 471 isosp[3]=-2; 472 } 473 void NNOmegaToMultiPionsChannel::pn_pnPipP 474 isosp[0]=2; 475 isosp[1]=2; 476 isosp[2]=-2; 477 isosp[3]=-2; 478 } 479 void NNOmegaToMultiPionsChannel::pp_pnPipP 480 isosp[0]=2; 481 isosp[1]=0; 482 isosp[2]=0; 483 isosp[3]=0; 484 iso1=1; 485 iso2=-1; 486 } 487 void NNOmegaToMultiPionsChannel::pn_nnPipP 488 isosp[0]=2; 489 isosp[1]=0; 490 isosp[2]=0; 491 isosp[3]=0; 492 iso1=-1; 493 iso2=-1; 494 } 495 void NNOmegaToMultiPionsChannel::pp_nnPipP 496 isosp[0]=2; 497 isosp[1]=0; 498 isosp[2]=0; 499 isosp[3]=0; 500 iso1=-1; 501 iso2=-1; 502 } 503 void NNOmegaToMultiPionsChannel::pp_pnPipP 504 isosp[0]=2; 505 isosp[1]=2; 506 isosp[2]=0; 507 isosp[3]=-2; 508 iso1=1; 509 iso2=-1; 510 } 511 void NNOmegaToMultiPionsChannel::pn_nnPipP 512 isosp[0]=2; 513 isosp[1]=2; 514 isosp[2]=0; 515 isosp[3]=-2; 516 iso1=-1; 517 iso2=-1; 518 } 519 void NNOmegaToMultiPionsChannel::pp_nnPipP 520 isosp[0]=2; 521 isosp[1]=2; 522 isosp[2]=0; 523 isosp[3]=-2; 524 iso1=-1; 525 iso2=-1; 526 } 527 void NNOmegaToMultiPionsChannel::nn_pnPi0P 528 isosp[0]=0; 529 isosp[1]=0; 530 isosp[2]=0; 531 isosp[3]=-2; 532 iso1=1; 533 iso2=-1; 534 } 535 void NNOmegaToMultiPionsChannel::pn_ppPi0P 536 isosp[0]=0; 537 isosp[1]=0; 538 isosp[2]=0; 539 isosp[3]=-2; 540 iso1=1; 541 iso2=1; 542 } 543 void NNOmegaToMultiPionsChannel::nn_pnPipP 544 isosp[0]=2; 545 isosp[1]=0; 546 isosp[2]=-2; 547 isosp[3]=-2; 548 iso1=1; 549 iso2=-1; 550 } 551 void NNOmegaToMultiPionsChannel::pn_ppPipP 552 isosp[0]=2; 553 isosp[1]=0; 554 isosp[2]=-2; 555 isosp[3]=-2; 556 iso1=1; 557 iso2=1; 558 } 559 560 void NNOmegaToMultiPionsChannel::inter2Par 561 562 if (Random::shoot() < p) std::swap(iso 563 564 } 565 566 567 } 568