221 {
222 G4DecayTable* decayTable;
223 G4VDecayChannel* mode;
224 G4ExcitedMesons* particle;
225
226
227
228
229
230
231
232
233
234
235
236 particle = new G4ExcitedMesons( "omega", 781.94 * MeV, 8.41 * MeV, +0.0 * eplus, 2, -1, -1,
237 0, +0, -1, "meson", 0, 0, 223, false, 0.0, NULL );
238 particle->SetAntiPDGEncoding( 223 );
239
240 particle->SetMultipletName( "omega" );
241
242 decayTable = new G4DecayTable();
243
244
245 mode = new G4PhaseSpaceDecayChannel( "omega", 0.888, 3, "pi+", "pi-", "pi0" );
246
247 decayTable->Insert( mode );
248 particle->SetDecayTable( decayTable );
249
250
251 particle = new G4ExcitedMesons( "phi", 1019.4 * MeV, 4.43 * MeV, +0.0 * eplus, 2, -1, -1, 0,
252 +0, -1, "meson", 0, 0, 333, false, 0.0, NULL );
253 particle->SetAntiPDGEncoding( 333 );
254
255 particle->SetMultipletName( "phi" );
256
257 decayTable = new G4DecayTable();
258
259
260 mode = new G4PhaseSpaceDecayChannel( "phi", 0.491, 2, "kaon+", "kaon-" );
261 decayTable->Insert( mode );
262
263
264 mode = new G4PhaseSpaceDecayChannel( "phi", 0.343, 2, "kaon0S", "kaon0L" );
265
266 decayTable->Insert( mode );
267
268
269 mode = new G4PhaseSpaceDecayChannel( "phi", 0.129, 2, "rho0", "pi0" );
270
271 decayTable->Insert( mode );
272 particle->SetDecayTable( decayTable );
273
274
275 particle = new G4ExcitedMesons( "rho+", 769.9 * MeV, 151.2 * MeV, +1.0 * eplus, 2, -1, -1, 2,
276 +2, +1, "meson", 0, 0, 213, false, 0.0, NULL );
277
278 particle->SetMultipletName( "rho" );
279
280 decayTable = new G4DecayTable();
281
282
283 mode = new G4PhaseSpaceDecayChannel( "rho+", 1.000, 2, "pi+", "pi0" );
284
285 decayTable->Insert( mode );
286 particle->SetDecayTable( decayTable );
287
288
289 particle = new G4ExcitedMesons( "rho-", 769.9 * MeV, 151.2 * MeV, -1.0 * eplus, 2, -1, -1, 2,
290 -2, +1, "meson", 0, 0, -213, false, 0.0, NULL );
291
292 particle->SetMultipletName( "rho" );
293
294 decayTable = new G4DecayTable();
295
296
297 mode = new G4PhaseSpaceDecayChannel( "rho-", 1.000, 2, "pi-", "pi0" );
298
299 decayTable->Insert( mode );
300 particle->SetDecayTable( decayTable );
301
302
303 particle = new G4ExcitedMesons( "rho0", 770.0 * MeV, 150.7 * MeV, 0.0, 2, -1, -1, 2, 0, +1,
304 "meson", 0, 0, 113,
false, 0.0 *
ns, NULL );
305 particle->SetAntiPDGEncoding( 113 );
306
307 particle->SetMultipletName( "rho" );
308
309 decayTable = new G4DecayTable();
310
311
312 mode = new G4PhaseSpaceDecayChannel( "rho0", 1.000, 2, "pi+", "pi-" );
313
314 decayTable->Insert( mode );
315 particle->SetDecayTable( decayTable );
316
317
318 particle = new G4ExcitedMesons( "a0(980)+", 984.8 * MeV, 60.0 * MeV, +1.0 * eplus, 0, +1, +1,
319 2, +2, -1, "meson", 0, 0, 9000211, false, 0.0, NULL );
320
321 particle->SetMultipletName( "a0(980)" );
322
323 decayTable = new G4DecayTable();
324
325
326 mode = new G4PhaseSpaceDecayChannel( "a0(980)+", 1.000, 2, "pi+", "eta" );
327
328 decayTable->Insert( mode );
329 particle->SetDecayTable( decayTable );
330
331
332 particle = new G4ExcitedMesons( "a0(980)-", 984.8 * MeV, 60.0 * MeV, -1.0 * eplus, 0, +1, +1,
333 2, -2, -1, "meson", 0, 0, -9000211, false, 0.0, NULL );
334
335 particle->SetMultipletName( "a0(980)" );
336
337 decayTable = new G4DecayTable();
338
339
340 mode = new G4PhaseSpaceDecayChannel( "a0(980)-", 1.000, 2, "pi-", "eta" );
341
342 decayTable->Insert( mode );
343 particle->SetDecayTable( decayTable );
344
345
346 particle = new G4ExcitedMesons( "a0(980)0", 984.8 * MeV, 60.0 * MeV, 0.0, 0, +1, +1, 2, 0,
347 -1, "meson", 0, 0, 9000111, false, 0.0, NULL );
348 particle->SetAntiPDGEncoding( 9000111 );
349
350 particle->SetMultipletName( "a0(980)" );
351
352 decayTable = new G4DecayTable();
353
354
355 mode = new G4PhaseSpaceDecayChannel( "a0(980)0", 1.000, 2, "pi0", "eta" );
356
357 decayTable->Insert( mode );
358 particle->SetDecayTable( decayTable );
359
360
361 particle = new G4ExcitedMesons( "f0(400-1200)", 470.0 * MeV, 400.0 * MeV, 0.0, 0, +1, +1, 0,
362 0, +1, "meson", 0, 0, 9000221, false, 0.0, NULL );
363 particle->SetAntiPDGEncoding( 9000221 );
364
365 particle->SetMultipletName( "f0(400-1200)" );
366
367 decayTable = new G4DecayTable();
368
369
370 mode = new G4PhaseSpaceDecayChannel( "f0(400-1200)", 1.000, 2, "pi+", "pi-" );
371
372 decayTable->Insert( mode );
373 particle->SetDecayTable( decayTable );
374
375
376 particle = new G4ExcitedMesons( "f0(980)", 980.0 * MeV, 40.0 * MeV, 0.0, 0, +1, +1, 0, 0, +1,
377 "meson", 0, 0, 9010221, false, 0.0, NULL );
378 particle->SetAntiPDGEncoding( 9010221 );
379
380 particle->SetMultipletName( "f0(980)" );
381
382 decayTable = new G4DecayTable();
383
384
385 mode = new G4PhaseSpaceDecayChannel( "f0(980)", 1.000, 2, "pi+", "pi-" );
386
387 decayTable->Insert( mode );
388 particle->SetDecayTable( decayTable );
389
390
391 particle = new G4ExcitedMesons( "k_star+", 891.6 * MeV, 50.8 * MeV, +1.0 * eplus, 2, -1, 0,
392 1, +1, 0, "meson", 0, 0, 323, false, 0.0, NULL );
393
394 particle->SetMultipletName( "k_star" );
395
396 decayTable = new G4DecayTable();
397
398
399 mode = new G4PhaseSpaceDecayChannel( "k_star+", 0.500, 2, "kaon+", "pi0" );
400
401 decayTable->Insert( mode );
402
403
404 mode = new G4PhaseSpaceDecayChannel( "k_star+", 0.500, 2, "kaon0", "pi+" );
405
406 decayTable->Insert( mode );
407 particle->SetDecayTable( decayTable );
408
409
410 particle = new G4ExcitedMesons( "k_star0", 896.1 * MeV, 50.5 * MeV, 0.0 * eplus, 2, -1, 0, 1,
411 -1, 0, "meson", 0, 0, 313, false, 0.0, NULL );
412
413 particle->SetMultipletName( "k_star" );
414
415 decayTable = new G4DecayTable();
416
417
418 mode = new G4PhaseSpaceDecayChannel( "k_star0", 0.500, 2, "kaon+", "pi-" );
419
420 decayTable->Insert( mode );
421
422
423 mode = new G4PhaseSpaceDecayChannel( "k_star0", 0.500, 2, "kaon0", "pi0" );
424
425 decayTable->Insert( mode );
426 particle->SetDecayTable( decayTable );
427
428
429 particle = new G4ExcitedMesons( "k_star-", 891.6 * MeV, 50.8 * MeV, -1.0 * eplus, 2, -1, 0,
430 1, +1, 0, "meson", 0, 0, -323, false, 0.0, NULL );
431
432 particle->SetMultipletName( "k_star" );
433
434 decayTable = new G4DecayTable();
435
436
437 mode = new G4PhaseSpaceDecayChannel( "k_star-", 0.500, 2, "kaon-", "pi0" );
438
439 decayTable->Insert( mode );
440
441
442 mode = new G4PhaseSpaceDecayChannel( "k_star-", 0.500, 2, "anti_kaon0", "pi+" );
443
444 decayTable->Insert( mode );
445 particle->SetDecayTable( decayTable );
446
447
448 particle = new G4ExcitedMesons( "anti_k_star0", 896.1 * MeV, 50.5 * MeV, 0.0 * eplus, 2, -1,
449 0, 1, -1, 0, "meson", 0, 0, -313, false, 0.0, NULL );
450
451 particle->SetMultipletName( "k_star" );
452
453 decayTable = new G4DecayTable();
454
455
456 mode = new G4PhaseSpaceDecayChannel( "anti_k_star0", 0.500, 2, "kaon-", "pi+" );
457
458 decayTable->Insert( mode );
459
460
461 mode = new G4PhaseSpaceDecayChannel( "anti_k_star0", 0.500, 2, "anti_kaon0", "pi0" );
462
463 decayTable->Insert( mode );
464 particle->SetDecayTable( decayTable );
465}