85 G4bool IsScatProjToProjCase,
104 adjointPrimKinEnergy,
106 IsScatProjToProjCase);
115 G4double projectileTotalEnergy = projectileM0+projectileKinEnergy;
116 G4double projectileP2 = projectileTotalEnergy*projectileTotalEnergy - projectileM0*projectileM0;
123 if (IsScatProjToProjCase) {
126 G4double companionTotalEnergy =companionM0+ projectileKinEnergy-adjointPrimKinEnergy;
127 G4double companionP2 = companionTotalEnergy*companionTotalEnergy - companionM0*companionM0;
132 G4double P_parallel = (adjointPrimP*adjointPrimP + projectileP2 - companionP2)/(2.*adjointPrimP);
133 G4double P_perp = std::sqrt( projectileP2 - P_parallel*P_parallel);
137 projectileMomentum.
rotateUz(dir_parallel);
141 if (!IsScatProjToProjCase ){
174 if (kinEnergyProj>Emin_proj && kinEnergyProj<=Emax_proj){
186 dSigmadEprod=(sigma1-sigma2)/dE;
192 if (dSigmadEprod>1.) {
195 G4cout<<
"dsigma "<<kinEnergyProj/
MeV<<
'\t'<<kinEnergyProd/
MeV<<
'\t'<<dSigmadEprod<<
G4endl;
212 G4double deltaKinEnergy = kinEnergyProd;
226 G4double etot2 = totEnergy*totEnergy;
227 G4double beta2 = kinEnergyProj*(kinEnergyProj + 2.0*
mass)/etot2;
230 f = 1.0 - beta2*deltaKinEnergy/Tmax;
232 f1 = 0.5*deltaKinEnergy*deltaKinEnergy/etot2;
242 G4cout <<
"### G4BetheBlochModel in Adjoint Sim WARNING: gg= " << gg
285 if (UsedFwdCS >0) new_weight*= CorrectFwdCS/UsedFwdCS;
294 new_weight*=projectileKinEnergy/adjointPrimKinEnergy;
310 pname !=
"deuteron" && pname !=
"triton") {
351 {
return PrimAdjEnergy+Tcut;