132 genericIonsManager->
GetIon(
"alpha+");
133 genericIonsManager->
GetIon(
"helium");
134 genericIonsManager->
GetIon(
"hydrogen");
147 G4int nreg = regnamesDNA.size();
155 G4cout <<
"### G4EmDNAPhysicsActivator::ConstructProcess for " << nreg
156 <<
" regions; DNA physics type " << typesDNA[0] <<
G4endl;
375 for(
G4int i = 0; i < nreg; ++i)
380 G4cout <<
"### DNA models type " << typesDNA[i]
381 <<
" are activated for G4Region " << reg <<
G4endl;
439 mod, reg, 0., elowest);
445 mod, reg, 0.0, elimel);
455 mod, reg, elim1, elimin);
459 mod, reg, 0.0, elim1);
464 mod, reg, 0.0, elim1);
468 mod, reg, elim1, elimin);
472 mod, reg, 0.0, elimvb);
476 mod, reg, 0.0, elimat);
505 em_config->
SetExtraEmModel(
"proton",
"CoulombScat", mod, reg, 0.0, emax);
510 mod, reg, 0.0, pminbb,
516 mod, reg, pminbb, emax,
521 mod, reg, 0.0, gmmax);
525 mod, reg, gmmax, pmax);
529 mod, reg, 0.0, gmmax);
533 mod, reg, gmmax, pmax);
537 mod, reg, 0.0, pmax);
541 mod, reg, 0.0, elimel);
546 mod, reg, 0.0, pmax);
550 mod, reg, 0.0, gmmax);
554 mod, reg, 0.0, pmax);
558 mod, reg, 0.0, elimel);
578 mod, reg, 0.0, pminbb,
584 mod, reg, pminbb, emax,
589 mod, reg, 0.0, iemax);
608 G4double pminbba = massRatio*pminbb;
610 G4cout <<
"AddHeliumModels0 for <" << reg <<
"> a1msc: " << a1msc <<
" a2msc: " << a2msc
611 <<
" elimel= " << elimel <<
" pminbba= " << pminbba <<
G4endl;
622 em_config->
SetExtraEmModel(
"alpha",
"CoulombScat", mod, reg, 0.0, emax);
629 mod, reg, 0.0, pminbba,
635 mod, reg, pminbba, emax,
640 mod, reg, 0.0, hemax);
644 mod, reg, 0.0, hemax);
648 mod, reg, 0.0, hemax);
652 mod, reg, 0.0, elimel);
664 em_config->
SetExtraEmModel(
"alpha+",
"CoulombScat", mod, reg, 0.0, emax);
671 mod, reg, 0.0, pminbba,
677 mod, reg, pminbba, emax,
682 mod, reg, 0.0, hemax);
686 mod, reg, 0.0, hemax);
690 mod, reg, 0.0, hemax);
694 mod, reg, 0.0, hemax);
698 mod, reg, 0.0, elimel);
704 mod, reg, 0.0, hemax);
708 mod, reg, 0.0, hemax);
712 mod, reg, 0.0, hemax);
716 mod, reg, 0.0, elimel);
726 for(
G4int i = 0; i < nproc; ++i) {
746 for(
G4int i = 0; i < nproc; ++i)