00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038 #include "G4ShortLivedConstructor.hh"
00039
00040 #include "G4SystemOfUnits.hh"
00041 #include "G4ParticleDefinition.hh"
00042 #include "G4ParticleTable.hh"
00043 #include "G4ShortLivedTable.hh"
00044 #include "G4PhaseSpaceDecayChannel.hh"
00045 #include "G4VDecayChannel.hh"
00046 #include "G4DecayTable.hh"
00047
00048 G4bool G4ShortLivedConstructor::isConstructed = false;
00049
00050 G4ShortLivedConstructor::G4ShortLivedConstructor()
00051 {
00052 }
00053
00054 G4ShortLivedConstructor::~G4ShortLivedConstructor()
00055 {
00056 }
00057
00058
00059 void G4ShortLivedConstructor::ConstructParticle()
00060 {
00061 if (!isConstructed){
00062 ConstructQuarks();
00063 ConstructResonances();
00064 isConstructed = true;
00065 }
00066 }
00067
00068 #include "G4DiQuarks.hh"
00069 #include "G4Quarks.hh"
00070 #include "G4Gluons.hh"
00071 void G4ShortLivedConstructor::ConstructQuarks()
00072 {
00073 G4ParticleDefinition* particle;
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084 particle = new G4Gluons(
00085 "gluon", 0.0*MeV, 0.0*MeV, 0.0*eplus,
00086 2, -1, 0,
00087 0, 0, 0,
00088 "gluons", 0, 0, 21,
00089 true, -1.0, NULL);
00090 particle->SetAntiPDGEncoding(21);
00091
00092 particle = new G4Quarks(
00093 "u_quark", 2.4*MeV, 0.0*MeV, 2./3.*eplus,
00094 1, +1, 0,
00095 1, +1, 0,
00096 "quarks", 0, 0, 2,
00097 true, -1.0, NULL);
00098
00099 particle = new G4Quarks(
00100 "d_quark", 4.9*MeV, 0.0*MeV, -1./3.*eplus,
00101 1, +1, 0,
00102 1, -1, 0,
00103 "quarks", 0, 0, 1,
00104 true, -1.0, NULL);
00105
00106 particle = new G4Quarks(
00107 "s_quark", 100.0*MeV, 0.0*MeV, -1./3.*eplus,
00108 1, +1, 0,
00109 0, 0, 0,
00110 "quarks", 0, 0, 3,
00111 true, -1.0, NULL);
00112
00113 particle = new G4Quarks(
00114 "c_quark", 1.29*GeV, 0.0*MeV, +2./3.*eplus,
00115 1, +1, 0,
00116 0, 0, 0,
00117 "quarks", 0, 0, 4,
00118 true, -1.0, NULL);
00119
00120 particle = new G4Quarks(
00121 "b_quark", 4.19*GeV, 0.0*MeV, -1./3.*eplus,
00122 1, +1, 0,
00123 0, 0, 0,
00124 "quarks", 0, 0, 5,
00125 true, -1.0, NULL);
00126
00127 particle = new G4Quarks(
00128 "t_quark", 172.9*GeV, 2.0*GeV, +2./3.*eplus,
00129 1, +1, 0,
00130 0, 0, 0,
00131 "quarks", 0, 0, 6,
00132 true, -1.0, NULL);
00133
00134 particle = new G4Quarks(
00135 "anti_u_quark", 2.4*MeV, 0.0*MeV, -2./3.*eplus,
00136 1, +1, 0,
00137 1, -1, 0,
00138 "quarks", 0, 0, -2,
00139 true, -1.0, NULL);
00140
00141 particle = new G4Quarks(
00142 "anti_d_quark", 4.9*MeV, 0.0*MeV, 1./3.*eplus,
00143 1, +1, 0,
00144 1, +1, 0,
00145 "quarks", 0, 0, -1,
00146 true, -1.0, NULL);
00147
00148 particle = new G4Quarks(
00149 "anti_s_quark", 100.0*MeV, 0.0*MeV, 1./3.*eplus,
00150 1, +1, 0,
00151 0, 0, 0,
00152 "quarks", 0, 0, -3,
00153 true, -1.0, NULL);
00154
00155 particle = new G4Quarks(
00156 "anti_c_quark", 1.29*GeV, 0.0*MeV, -2./3.*eplus,
00157 1, +1, 0,
00158 0, 0, 0,
00159 "quarks", 0, 0, -4,
00160 true, -1.0, NULL);
00161
00162 particle = new G4Quarks(
00163 "anti_b_quark", 4.19*GeV, 0.0*MeV, 1./3.*eplus,
00164 1, +1, 0,
00165 0, 0, 0,
00166 "quarks", 0, 0, -5,
00167 true, -1.0, NULL);
00168
00169 particle = new G4Quarks(
00170 "anti_t_quark", 172.9*GeV, 2.0*GeV, -2./3.*eplus,
00171 1, +1, 0,
00172 0, 0, 0,
00173 "quarks", 0, 0, -6,
00174 true, -1.0, NULL);
00175
00176
00177 particle = new G4DiQuarks(
00178 "uu1_diquark", 4.8*MeV, 0.0*MeV, 4./3.*eplus,
00179 2, +1, 0,
00180 2, +2, 0,
00181 "diquarks", 0, 0, 2203,
00182 true, -1.0, NULL);
00183
00184 particle = new G4DiQuarks(
00185 "ud1_diquark", 7.2*MeV, 0.0*MeV, 1./3.*eplus,
00186 2, +1, 0,
00187 2, +0, 0,
00188 "diquarks", 0, 0, 2103,
00189 true, -1.0, NULL);
00190
00191 particle = new G4DiQuarks(
00192 "dd1_diquark", 9.6*MeV, 0.0*MeV, -2./3.*eplus,
00193 2, +1, 0,
00194 2, -2, 0,
00195 "diquarks", 0, 0, 1103,
00196 true, -1.0, NULL);
00197
00198
00199 particle = new G4DiQuarks(
00200 "ud0_diquark", 7.3*MeV, 0.0*MeV, 1./3.*eplus,
00201 0, +1, 0,
00202 0, +0, 0,
00203 "diquarks", 0, 0, 2101,
00204 true, -1.0, NULL);
00205
00206
00207 particle = new G4DiQuarks(
00208 "sd1_diquark", 108.8*MeV, 0.0*MeV, -2./3.*eplus,
00209 2, +1, 0,
00210 1, -1, 0,
00211 "diquarks", 0, 0, 3103,
00212 true, -1.0, NULL);
00213
00214
00215 particle = new G4DiQuarks(
00216 "su1_diquark", 106.4*MeV, 0.0*MeV, 1./3.*eplus,
00217 2, +1, 0,
00218 1, +1, 0,
00219 "diquarks", 0, 0, 3203,
00220 true, -1.0, NULL);
00221
00222
00223 particle = new G4DiQuarks(
00224 "sd0_diquark", 108.0*MeV, 0.0*MeV, -2./3.*eplus,
00225 0, +1, 0,
00226 1, -1, 0,
00227 "diquarks", 0, 0, 3101,
00228 true, -1.0, NULL);
00229
00230
00231 particle = new G4DiQuarks(
00232 "su0_diquark", 106.4*MeV, 0.0*MeV, 1./3.*eplus,
00233 0, +1, 0,
00234 1, +1, 0,
00235 "diquarks", 0, 0, 3201,
00236 true, -1.0, NULL);
00237
00238
00239 particle = new G4DiQuarks(
00240 "anti_uu1_diquark", 4.8*MeV, 0.0*MeV, -4./3.*eplus,
00241 2, +1, 0,
00242 2, -2, 0,
00243 "diquarks", 0, 0, -2203,
00244 true, -1.0, NULL);
00245
00246 particle = new G4DiQuarks(
00247 "anti_ud1_diquark", 7.2*MeV, 0.0*MeV, -1./3.*eplus,
00248 2, +1, 0,
00249 2, +0, 0,
00250 "diquarks", 0, 0, -2103,
00251 true, -1.0, NULL);
00252
00253 particle = new G4DiQuarks(
00254 "anti_dd1_diquark", 9.6*MeV, 0.0*MeV, 2./3.*eplus,
00255 2, +1, 0,
00256 2, +2, 0,
00257 "diquarks", 0, 0, -1103,
00258 true, -1.0, NULL);
00259
00260
00261 particle = new G4DiQuarks(
00262 "anti_ud0_diquark", 7.3*MeV, 0.0*MeV, -1./3.*eplus,
00263 0, +1, 0,
00264 0, +0, 0,
00265 "diquarks", 0, 0, -2101,
00266 true, -1.0, NULL);
00267
00268
00269 particle = new G4DiQuarks(
00270 "anti_sd1_diquark", 108.8*MeV, 0.0*MeV, 2./3.*eplus,
00271 2, +1, 0,
00272 1, +1, 0,
00273 "diquarks", 0, 0, -3103,
00274 true, -1.0, NULL);
00275
00276
00277 particle = new G4DiQuarks(
00278 "anti_su1_diquark", 106.4*MeV, 0.0*MeV, -1./3.*eplus,
00279 2, +1, 0,
00280 1, -1, 0,
00281 "diquarks", 0, 0, -3203,
00282 true, -1.0, NULL);
00283
00284
00285 particle = new G4DiQuarks(
00286 "anti_sd0_diquark", 108.0*MeV, 0.0*MeV, 2./3.*eplus,
00287 0, +1, 0,
00288 1, +1, 0,
00289 "diquarks", 0, 0, -3101,
00290 true, -1.0, NULL);
00291
00292
00293 particle = new G4DiQuarks(
00294 "anti_su0_diquark", 106.4*MeV, 0.0*MeV, -1./3.*eplus,
00295 0, +1, 0,
00296 1, -1, 0,
00297 "diquarks", 0, 0, -3201,
00298 true, -1.0, NULL);
00299
00300 particle = new G4DiQuarks(
00301 "ss1_diquark", 208.0*MeV, 0.0*MeV, -2./3.*eplus,
00302 2, +1, 0,
00303 0, 0, 0,
00304 "diquarks", 0, 0, 3303,
00305 true, -1.0, NULL);
00306
00307
00308 particle = new G4DiQuarks(
00309 "anti_ss1_diquark", 208.0*MeV, 0.0*MeV, 2./3.*eplus,
00310 2, +1, 0,
00311 0, 0, 0,
00312 "diquarks", 0, 0, -3303,
00313 true, -1.0, NULL);
00314
00315 particle = NULL;
00316 }
00317
00318 #include "G4ExcitedNucleonConstructor.hh"
00319 #include "G4ExcitedDeltaConstructor.hh"
00320 #include "G4ExcitedLambdaConstructor.hh"
00321 #include "G4ExcitedSigmaConstructor.hh"
00322 #include "G4ExcitedXiConstructor.hh"
00323 #include "G4ExcitedMesonConstructor.hh"
00324 void G4ShortLivedConstructor::ConstructResonances()
00325 {
00326 ConstructBaryons();
00327 ConstructMesons();
00328
00329
00330 G4ExcitedNucleonConstructor nucleons;
00331 nucleons.Construct();
00332
00333
00334 G4ExcitedDeltaConstructor deltas;
00335 deltas.Construct();
00336
00337
00338 G4ExcitedLambdaConstructor lamdas;
00339 lamdas.Construct();
00340
00341
00342 G4ExcitedSigmaConstructor sigmas;
00343 sigmas.Construct();
00344
00345
00346 G4ExcitedXiConstructor xis;
00347 xis.Construct();
00348
00349
00350 G4ExcitedMesonConstructor mesons;
00351 mesons.Construct();
00352
00353 }
00354
00355
00356 #include "G4ExcitedBaryons.hh"
00357 void G4ShortLivedConstructor::ConstructBaryons()
00358 {
00359 G4DecayTable* decayTable;
00360 G4VDecayChannel* mode;
00361 G4ExcitedBaryons* particle;
00362
00363
00364
00365
00366
00367
00368
00369
00370
00371
00372
00373 particle = new G4ExcitedBaryons(
00374 "delta++", 1.232*GeV, 120.0*MeV, +2.0*eplus,
00375 3, +1, 0,
00376 3, +3, 0,
00377 "baryon", 0, +1, 2224,
00378 false, 0.0, NULL);
00379
00380 particle->SetMultipletName("delta");
00381
00382 decayTable = new G4DecayTable();
00383
00384
00385 mode = new G4PhaseSpaceDecayChannel("delta++",1.000, 2,
00386 "proton","pi+");
00387
00388 decayTable->Insert(mode);
00389 particle->SetDecayTable(decayTable);
00390
00391
00392 particle = new G4ExcitedBaryons(
00393 "delta+", 1.232*GeV, 120.0*MeV, +1.0*eplus,
00394 3, +1, 0,
00395 3, +1, 0,
00396 "baryon", 0, +1, 2214,
00397 false, 0.0, NULL);
00398
00399 particle->SetMultipletName("delta(1232)");
00400
00401 decayTable = new G4DecayTable();
00402
00403
00404 mode = new G4PhaseSpaceDecayChannel("delta+", 0.01, 2,
00405 "proton","gamma");
00406 decayTable->Insert(mode);
00407
00408
00409
00410
00411 mode = new G4PhaseSpaceDecayChannel("delta+", 0.495, 2,
00412 "proton","pi0");
00413 decayTable->Insert(mode);
00414
00415
00416 mode = new G4PhaseSpaceDecayChannel("delta+", 0.495, 2,
00417 "neutron","pi+");
00418 decayTable->Insert(mode);
00419 particle->SetDecayTable(decayTable);
00420
00421
00422 particle = new G4ExcitedBaryons(
00423 "delta0", 1.232*GeV, 120.0*MeV, +0.0*eplus,
00424 3, +1, 0,
00425 3, -1, 0,
00426 "baryon", 0, +1, 2114,
00427 false, 0.0, NULL);
00428
00429 particle->SetMultipletName("delta(1232)");
00430
00431 decayTable = new G4DecayTable();
00432
00433
00434 mode = new G4PhaseSpaceDecayChannel("delta0", 0.01, 2,
00435 "neutron","gamma");
00436 decayTable->Insert(mode);
00437
00438
00439 mode = new G4PhaseSpaceDecayChannel("delta0", 0.495, 2,
00440 "proton","pi-");
00441 decayTable->Insert(mode);
00442
00443
00444 mode = new G4PhaseSpaceDecayChannel("delta0", 0.495, 2,
00445 "neutron","pi0");
00446 decayTable->Insert(mode);
00447 particle->SetDecayTable(decayTable);
00448
00449
00450 particle = new G4ExcitedBaryons(
00451 "delta-", 1.232*GeV, 118.0*MeV, -1.0*eplus,
00452 3, +1, 0,
00453 3, -3, 0,
00454 "baryon", 0, +1, 1114,
00455 false, 0.0, NULL);
00456
00457 particle->SetMultipletName("delta(1232)");
00458
00459 decayTable = new G4DecayTable();
00460
00461
00462 mode = new G4PhaseSpaceDecayChannel("delta-", 1.000, 2,
00463 "neutron","pi-");
00464 decayTable->Insert(mode);
00465 particle->SetDecayTable(decayTable);
00466
00467
00469
00470
00471 particle = new G4ExcitedBaryons(
00472 "anti_delta++", 1.232*GeV, 120.0*MeV, -2.0*eplus,
00473 3, +1, 0,
00474 3, -3, 0,
00475 "baryon", 0, -1, -2224,
00476 false, 0.0, NULL);
00477
00478 particle->SetMultipletName("delta(1232)");
00479
00480 decayTable = new G4DecayTable();
00481
00482
00483 mode = new G4PhaseSpaceDecayChannel("anti_delta++",1.000, 2,
00484 "anti_proton","pi-");
00485
00486 decayTable->Insert(mode);
00487 particle->SetDecayTable(decayTable);
00488
00489
00490 particle = new G4ExcitedBaryons(
00491 "anti_delta+", 1.232*GeV, 120.0*MeV, -1.0*eplus,
00492 3, +1, 0,
00493 3, -1, 0,
00494 "baryon", 0, -1, -2214,
00495 false, 0.0, NULL);
00496
00497 particle->SetMultipletName("delta(1232)");
00498
00499 decayTable = new G4DecayTable();
00500
00501
00502 mode = new G4PhaseSpaceDecayChannel("anti_delta+", 0.500, 2,
00503 "anti_proton","pi0");
00504 decayTable->Insert(mode);
00505
00506
00507 mode = new G4PhaseSpaceDecayChannel("anti_delta+", 0.500, 2,
00508 "anti_neutron","pi-");
00509 decayTable->Insert(mode);
00510 particle->SetDecayTable(decayTable);
00511
00512
00513 particle = new G4ExcitedBaryons(
00514 "anti_delta0", 1.232*GeV, 120.0*MeV, +0.0*eplus,
00515 3, +1, 0,
00516 3, +1, 0,
00517 "baryon", 0, -1, -2114,
00518 false, 0.0, NULL);
00519
00520 particle->SetMultipletName("delta(1232)");
00521
00522 decayTable = new G4DecayTable();
00523
00524
00525 mode = new G4PhaseSpaceDecayChannel("anti_delta0", 0.500, 2,
00526 "anti_proton","pi+");
00527 decayTable->Insert(mode);
00528
00529
00530 mode = new G4PhaseSpaceDecayChannel("anti_delta0", 0.500, 2,
00531 "anti_neutron","pi0");
00532 decayTable->Insert(mode);
00533 particle->SetDecayTable(decayTable);
00534
00535
00536 particle = new G4ExcitedBaryons(
00537 "anti_delta-", 1.232*GeV, 118.0*MeV, +1.0*eplus,
00538 3, +1, 0,
00539 3, +3, 0,
00540 "baryon", 0, -1, -1114,
00541 false, 0.0, NULL);
00542
00543 particle->SetMultipletName("delta(1232)");
00544
00545 decayTable = new G4DecayTable();
00546
00547
00548 mode = new G4PhaseSpaceDecayChannel("anti_delta-", 1.000, 2,
00549 "anti_neutron","pi+");
00550 decayTable->Insert(mode);
00551 particle->SetDecayTable(decayTable);
00552
00553
00554
00555 }
00556 #include "G4ExcitedMesons.hh"
00557 void G4ShortLivedConstructor::ConstructMesons()
00558 {
00559 G4DecayTable* decayTable;
00560 G4VDecayChannel* mode;
00561 G4ExcitedMesons* particle;
00562
00563
00564
00565
00566
00567
00568
00569
00570
00571
00572
00573 particle = new G4ExcitedMesons(
00574 "omega", 782.65*MeV, 8.49*MeV, +0.0*eplus,
00575 2, -1, -1,
00576 0, +0, -1,
00577 "meson", 0, 0, 223,
00578 false, 0.0, NULL);
00579 particle->SetAntiPDGEncoding(223);
00580
00581 particle->SetMultipletName("omega");
00582
00583 decayTable = new G4DecayTable();
00584
00585
00586 mode = new G4PhaseSpaceDecayChannel("omega",0.891, 3,
00587 "pi+","pi-","pi0");
00588
00589 decayTable->Insert(mode);
00590
00591
00592
00593 mode = new G4PhaseSpaceDecayChannel("omega",0.0890, 2,
00594 "gamma","pi0");
00595
00596 decayTable->Insert(mode);
00597
00598
00599
00600 mode = new G4PhaseSpaceDecayChannel("omega",0.0170, 2,
00601 "pi+","pi-");
00602
00603 decayTable->Insert(mode);
00604 particle->SetDecayTable(decayTable);
00605
00606
00607 particle = new G4ExcitedMesons(
00608 "phi", 1019.46*MeV, 4.26*MeV, +0.0*eplus,
00609 2, -1, -1,
00610 0, +0, -1,
00611 "meson", 0, 0, 333,
00612 false, 0.0, NULL);
00613 particle->SetAntiPDGEncoding(333);
00614
00615 particle->SetMultipletName("phi");
00616
00617 decayTable = new G4DecayTable();
00618
00619
00620 mode = new G4PhaseSpaceDecayChannel("phi",0.492, 2,
00621 "kaon+","kaon-");
00622 decayTable->Insert(mode);
00623
00624
00625 mode = new G4PhaseSpaceDecayChannel("phi",0.340, 2,
00626 "kaon0S","kaon0L");
00627
00628 decayTable->Insert(mode);
00629
00630
00631 mode = new G4PhaseSpaceDecayChannel("phi",0.153, 2,
00632 "rho0","pi0");
00633
00634 decayTable->Insert(mode);
00635 particle->SetDecayTable(decayTable);
00636
00637
00638 particle = new G4ExcitedMesons(
00639 "rho+", 775.8*MeV, 150.3*MeV, +1.0*eplus,
00640 2, -1, -1,
00641 2, +2, +1,
00642 "meson", 0, 0, 213,
00643 false, 0.0, NULL);
00644
00645 particle->SetMultipletName("rho");
00646
00647 decayTable = new G4DecayTable();
00648
00649
00650 mode = new G4PhaseSpaceDecayChannel("rho+",1.000, 2,
00651 "pi+","pi0");
00652
00653 decayTable->Insert(mode);
00654 particle->SetDecayTable(decayTable);
00655
00656
00657 particle = new G4ExcitedMesons(
00658 "rho-", 775.8*MeV, 150.3*MeV, -1.0*eplus,
00659 2, -1, -1,
00660 2, -2, +1,
00661 "meson", 0, 0, -213,
00662 false, 0.0, NULL);
00663
00664 particle->SetMultipletName("rho");
00665
00666 decayTable = new G4DecayTable();
00667
00668
00669 mode = new G4PhaseSpaceDecayChannel("rho-",1.000, 2,
00670 "pi-","pi0");
00671
00672 decayTable->Insert(mode);
00673 particle->SetDecayTable(decayTable);
00674
00675
00676 particle = new G4ExcitedMesons(
00677 "rho0", 775.49*MeV, 149.1*MeV, 0.0,
00678 2, -1, -1,
00679 2, 0, +1,
00680 "meson", 0, 0, 113,
00681 false, 0.0*ns, NULL );
00682 particle->SetAntiPDGEncoding(113);
00683
00684 particle->SetMultipletName("rho");
00685
00686 decayTable = new G4DecayTable();
00687
00688
00689 mode = new G4PhaseSpaceDecayChannel("rho0",1.000, 2,
00690 "pi+","pi-");
00691
00692 decayTable->Insert(mode);
00693 particle->SetDecayTable(decayTable);
00694
00695
00696 particle = new G4ExcitedMesons(
00697 "a0(980)+", 980.0*MeV, 60.0*MeV, +1.0*eplus,
00698 0, +1, +1,
00699 2, +2, -1,
00700 "meson", 0, 0, 9000211,
00701 false, 0.0, NULL);
00702
00703 particle->SetMultipletName("a0(980)");
00704
00705 decayTable = new G4DecayTable();
00706
00707
00708 mode = new G4PhaseSpaceDecayChannel("a0(980)+",1.000, 2,
00709 "pi+","eta");
00710
00711 decayTable->Insert(mode);
00712 particle->SetDecayTable(decayTable);
00713
00714
00715 particle = new G4ExcitedMesons(
00716 "a0(980)-", 980.0*MeV, 60.0*MeV, -1.0*eplus,
00717 0, +1, +1,
00718 2, -2, -1,
00719 "meson", 0, 0, -9000211,
00720 false, 0.0, NULL);
00721
00722 particle->SetMultipletName("a0(980)");
00723
00724 decayTable = new G4DecayTable();
00725
00726
00727 mode = new G4PhaseSpaceDecayChannel("a0(980)-",1.000, 2,
00728 "pi-","eta");
00729
00730 decayTable->Insert(mode);
00731 particle->SetDecayTable(decayTable);
00732
00733
00734 particle = new G4ExcitedMesons(
00735 "a0(980)0", 980.0*MeV, 75.0*MeV, 0.0,
00736 0, +1, +1,
00737 2, 0, -1,
00738 "meson", 0, 0, 9000111,
00739 false, 0.0, NULL);
00740 particle->SetAntiPDGEncoding(9000111);
00741
00742 particle->SetMultipletName("a0(980)");
00743
00744 decayTable = new G4DecayTable();
00745
00746
00747 mode = new G4PhaseSpaceDecayChannel("a0(980)0",1.000, 2,
00748 "pi0","eta");
00749
00750 decayTable->Insert(mode);
00751 particle->SetDecayTable(decayTable);
00752
00753
00754 particle = new G4ExcitedMesons(
00755 "f0(600)", 800.0*MeV, 800.0*MeV, 0.0,
00756 0, +1, +1,
00757 0, 0, +1,
00758 "meson", 0, 0, 9000221,
00759 false, 0.0, NULL);
00760 particle->SetAntiPDGEncoding(9000221);
00761
00762 particle->SetMultipletName("f0(600)");
00763
00764 decayTable = new G4DecayTable();
00765
00766
00767 mode = new G4PhaseSpaceDecayChannel("f0(600)",1.000, 2,
00768 "pi+","pi-");
00769
00770 decayTable->Insert(mode);
00771 particle->SetDecayTable(decayTable);
00772
00773
00774
00775 particle = new G4ExcitedMesons(
00776 "f0(980)", 980.0*MeV, 70.0*MeV, 0.0,
00777 0, +1, +1,
00778 0, 0, +1,
00779 "meson", 0, 0, 9010221,
00780 false, 0.0, NULL);
00781 particle->SetAntiPDGEncoding(9010221);
00782
00783 particle->SetMultipletName("f0(980)");
00784
00785 decayTable = new G4DecayTable();
00786
00787
00788 mode = new G4PhaseSpaceDecayChannel("f0(980)",1.000, 2,
00789 "pi+","pi-");
00790
00791 decayTable->Insert(mode);
00792 particle->SetDecayTable(decayTable);
00793
00794
00795 particle = new G4ExcitedMesons(
00796 "eta(1405)", 1409.8*MeV, 51.1*MeV, 0.0,
00797 0, -1, +1,
00798 0, 0, +1,
00799 "meson", 0, 0, 9020221,
00800 false, 0.0, NULL);
00801 particle->SetAntiPDGEncoding(9020221);
00802
00803 particle->SetMultipletName("eta(1405)");
00804
00805 decayTable = new G4DecayTable();
00806
00807
00808 mode = new G4PhaseSpaceDecayChannel("eta(1405)",1.000, 2,
00809 "rho+","rho-");
00810
00811 decayTable->Insert(mode);
00812 particle->SetDecayTable(decayTable);
00813
00814 particle = new G4ExcitedMesons(
00815 "f0(1500)", 1505.0*MeV, 109.0*MeV, 0.0,
00816 0, +1, +1,
00817 0, 0, +1,
00818 "meson", 0, 0, 9030221,
00819 false, 0.0, NULL);
00820 particle->SetAntiPDGEncoding(9030221);
00821
00822 particle->SetMultipletName("f0(1500)");
00823
00824 decayTable = new G4DecayTable();
00825
00826
00827 mode = new G4PhaseSpaceDecayChannel("f0(1500)",1.000, 2,
00828 "eta","eta");
00829
00830 decayTable->Insert(mode);
00831 particle->SetDecayTable(decayTable);
00832
00833
00834 particle = new G4ExcitedMesons(
00835 "f0(1710)", 1720.0*MeV, 135.0*MeV, 0.0,
00836 0, +1, +1,
00837 0, 0, +1,
00838 "meson", 0, 0, 10331,
00839 false, 0.0, NULL);
00840 particle->SetAntiPDGEncoding(10331);
00841
00842 particle->SetMultipletName("f0(1710)");
00843
00844 decayTable = new G4DecayTable();
00845
00846
00847
00848 mode = new G4PhaseSpaceDecayChannel("f0(1710)",0.40, 2,
00849 "kaon0S","kaon0S");
00850
00851 decayTable->Insert(mode);
00852
00853
00854
00855 mode = new G4PhaseSpaceDecayChannel("f0(1710)",0.40, 2,
00856 "kaon+","kaon-");
00857
00858 decayTable->Insert(mode);
00859
00860
00861
00862 mode = new G4PhaseSpaceDecayChannel("f0(1710)",0.20, 2,
00863 "eta","eta");
00864
00865 decayTable->Insert(mode);
00866 particle->SetDecayTable(decayTable);
00867
00868
00869
00870 particle = new G4ExcitedMesons(
00871 "k_star+", 891.66*MeV, 50.8*MeV, +1.0*eplus,
00872 2, -1, 0,
00873 1, +1, 0,
00874 "meson", 0, 0, 323,
00875 false, 0.0, NULL);
00876
00877 particle->SetMultipletName("k_star");
00878
00879 decayTable = new G4DecayTable();
00880
00881
00882 mode = new G4PhaseSpaceDecayChannel("k_star+",0.500, 2,
00883 "kaon+","pi0");
00884
00885 decayTable->Insert(mode);
00886
00887
00888 mode = new G4PhaseSpaceDecayChannel("k_star+",0.500, 2,
00889 "kaon0","pi+");
00890
00891 decayTable->Insert(mode);
00892 particle->SetDecayTable(decayTable);
00893
00894
00895 particle = new G4ExcitedMesons(
00896 "k_star0", 895.94*MeV, 48.7*MeV, 0.0*eplus,
00897 2, -1, 0,
00898 1, -1, 0,
00899 "meson", 0, 0, 313,
00900 false, 0.0, NULL);
00901
00902 particle->SetMultipletName("k_star");
00903
00904 decayTable = new G4DecayTable();
00905
00906
00907 mode = new G4PhaseSpaceDecayChannel("k_star0",0.500, 2,
00908 "kaon+","pi-");
00909
00910 decayTable->Insert(mode);
00911
00912
00913 mode = new G4PhaseSpaceDecayChannel("k_star0",0.500, 2,
00914 "kaon0","pi0");
00915
00916 decayTable->Insert(mode);
00917 particle->SetDecayTable(decayTable);
00918
00919
00920 particle = new G4ExcitedMesons(
00921 "k_star-", 891.66*MeV, 50.8*MeV, -1.0*eplus,
00922 2, -1, 0,
00923 1, +1, 0,
00924 "meson", 0, 0, -323,
00925 false, 0.0, NULL);
00926
00927 particle->SetMultipletName("k_star");
00928
00929 decayTable = new G4DecayTable();
00930
00931
00932 mode = new G4PhaseSpaceDecayChannel("k_star-",0.500, 2,
00933 "kaon-","pi0");
00934
00935 decayTable->Insert(mode);
00936
00937
00938 mode = new G4PhaseSpaceDecayChannel("k_star-",0.500, 2,
00939 "anti_kaon0","pi-");
00940
00941 decayTable->Insert(mode);
00942 particle->SetDecayTable(decayTable);
00943
00944
00945
00946 particle = new G4ExcitedMesons(
00947 "anti_k_star0", 895.94*MeV, 48.7*MeV, 0.0*eplus,
00948 2, -1, 0,
00949 1, -1, 0,
00950 "meson", 0, 0, -313,
00951 false, 0.0, NULL);
00952
00953 particle->SetMultipletName("k_star");
00954
00955 decayTable = new G4DecayTable();
00956
00957
00958 mode = new G4PhaseSpaceDecayChannel("anti_k_star0",0.500, 2,
00959 "kaon-","pi+");
00960
00961 decayTable->Insert(mode);
00962
00963
00964 mode = new G4PhaseSpaceDecayChannel("anti_k_star0",0.500, 2,
00965 "anti_kaon0","pi0");
00966
00967 decayTable->Insert(mode);
00968 particle->SetDecayTable(decayTable);
00969
00970 }
00971
00972
00973
00974
00975
00976
00977
00978
00979
00980
00981
00982
00983
00984
00985
00986
00987
00988
00989