73 finalBaryon(0), initialCharge(0), finalCharge(0), initialStrange(0),
80 absoluteLimit(absolute), initialBaryon(0), finalBaryon(0),
81 initialCharge(0), finalCharge(0), initialStrange(0),
91 G4cout <<
" >>> G4CascadeCheckBalance(" <<
theName <<
")::collide"
112 ((pbullet ? pbullet->baryon() : nbullet ? nbullet->
getA() : 0) +
113 (ptarget ? ptarget->
baryon() : ntarget ? ntarget->
getA() : 0) );
132 <<
" final px " <<
final.px() <<
" py " <<
final.py()
133 <<
" pz " <<
final.pz() <<
" E " <<
final.e()
143 G4cout <<
" >>> G4CascadeCheckBalance(" <<
theName <<
")::collide(<FRAG>)"
164 <<
" final px " <<
final.px() <<
" py " <<
final.py()
165 <<
" pz " <<
final.pz() <<
" E " <<
final.e()
176 const std::vector<G4InuclElementaryParticle>& particles) {
178 G4cout <<
" >>> G4CascadeCheckBalance(" <<
theName <<
")::collide(<vector>)"
187 const std::vector<G4InuclElementaryParticle>& particles) {
190 <<
")::collide(<FRAG>,<vector>)" <<
G4endl;
201 const std::vector<G4InuclNuclei>& fragments) {
203 G4cout <<
" >>> G4CascadeCheckBalance(" <<
theName <<
")::collide(<vector>)"
216 const std::vector<G4CascadParticle>& particles) {
219 <<
")::collide(<cparticles>)" <<
G4endl;
232 const std::vector<G4CascadParticle>& cparticles) {
235 <<
")::collide(<EP>,<CP>)" <<
G4endl;
252 << (relokay ?
" conserved" :
" VIOLATED")
253 <<
" absolute " <<
deltaE()
254 << (absokay ?
" conserved" :
" VIOLATED") <<
G4endl;
257 <<
" conserved absolute " <<
deltaE() <<
" conserved" <<
G4endl;
260 return (relokay && absokay);
269 <<
relativeKE() << (relokay ?
" conserved" :
" VIOLATED")
271 << (absokay ?
" conserved" :
" VIOLATED") <<
G4endl;
275 <<
" conserved" <<
G4endl;
278 return (relokay && absokay);
287 << (relokay ?
" conserved" :
" VIOLATED")
288 <<
" absolute " <<
deltaP()
289 << (absokay ?
" conserved" :
" VIOLATED") <<
G4endl;
292 <<
" conserved absolute " <<
deltaP() <<
" conserved" <<
G4endl;
295 return (relokay && absokay);
const G4LorentzVector & GetMomentum() const
void collide(G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
G4bool momentumOkay() const
G4int getTotalStrangeness() const
void add(const G4CollisionOutput &right)
G4bool chargeOkay() const
void addOutgoingParticles(const std::vector< G4InuclElementaryParticle > &particles)
G4bool strangeOkay() const
G4bool baryonOkay() const
G4double relativeKE() const
void addOutgoingNuclei(const std::vector< G4InuclNuclei > &nuclea)
G4CollisionOutput tempOutput
G4double relativeE() const
G4int getStrangeness() const
static const G4double tolerance
G4CascadeCheckBalance(const G4String &owner="G4CascadeCheckBalance")
G4int getTotalCharge() const
G4GLOB_DLL std::ostream G4cerr
G4double relativeP() const
G4double getCharge() const
G4LorentzVector getMomentum() const
G4GLOB_DLL std::ostream G4cout
G4int getTotalBaryonNumber() const
G4bool energyOkay() const
G4LorentzVector getTotalOutputMomentum() const