55 applyARM(
true), photonEvaporation(aPhotoEvap)
112 if (shellIndex > -1) {
117 if (shellIndex >= nShells) shellIndex = nShells;
120 std::vector<G4DynamicParticle*> armProducts;
131 for (
G4int i = 0; i <
G4int(armProducts.size()); i++)
132 productEnergy += armProducts[i]->GetKineticEnergy();
138 G4double sinTh = std::sqrt(1.- cosTh*cosTh);
142 sinTh*std::cos(phi), cosTh);
146 armProducts.push_back(extra);
149 G4int nArm = armProducts.size();
152 for (
G4int i = 0; i < nArm; ++i) {
186 <<
" + gammas (or electrons), with branching ratio " <<
GetBR()