83 std::pair<G4double,G4double> ratios;
89 G4cout <<
"G4QuasiElasticChannel::ratios " << ratios.first <<
" x " <<ratios.second
90 <<
" = " << ratios.first*ratios.second <<
G4endl;
93 return ratios.first*ratios.second;
109 }
while (index < 0 || index >= static_cast<G4int>(nucleons.size()));
127 G4cout<<
"G4QElChan::Scatter: neutron - proton? A ="<<A<<
", Z="<<Z<<
", projName="
132 G4double residualNucleusEnergy=std::sqrt(
sqr(residualNucleusMass) +
134 pNucleon.
setE(targetNucleusMass-residualNucleusEnergy);
137 std::pair<G4LorentzVector,G4LorentzVector> result;
143 if (result.first.e() > 0.)
144 scatteredHadron4Mom=result.second;
155 - result.first - result.second;
159 G4cout <<
"Warning - G4QuasiElasticChannel::Scatter E-p non conservation : "
160 << EpConservation <<
G4endl;
167 ktv->push_back(sPrim);
168 if (result.first.e() > 0.)
171 ktv->push_back(sNuc);
177 ktv->push_back(rNuc);
181 residualNucleus4Mom/=resA;
186 ktv->push_back(rNuc);
190 G4cout<<
"G4QElC::Scat: Nucleon: "<<result.first <<
" mass "<<result.first.mag() <<
G4endl;
191 G4cout<<
"G4QElC::Scat: Project: "<<result.second<<
" mass "<<result.second.mag()<<
G4endl;