G4CascadeCheckBalance Class Reference

#include <G4CascadeCheckBalance.hh>

Inheritance diagram for G4CascadeCheckBalance:

Public Member Functions
	G4CascadeCheckBalance (const char *owner="G4CascadeCheckBalance")
	G4CascadeCheckBalance (G4double relative, G4double absolute, const char *owner="G4CascadeCheckBalance")
virtual	~G4CascadeCheckBalance ()
void	setOwner (const char *owner)
void	setLimits (G4double relative, G4double absolute)
void	setRelativeLimit (G4double limit)
void	setAbsoluteLimit (G4double limit)
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
void	collide (const G4Fragment &fragment, G4CollisionOutput &output)
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, const std::vector< G4InuclElementaryParticle > &particles)
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, const std::vector< G4CascadParticle > &particles)
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output, const std::vector< G4CascadParticle > &cparticles)
void	collide (const G4Fragment &target, const std::vector< G4InuclElementaryParticle > &particles)
void	collide (const G4Fragment &target, const std::vector< G4InuclNuclei > &fragments)
G4bool	energyOkay () const
G4bool	ekinOkay () const
G4bool	momentumOkay () const
G4bool	baryonOkay () const
G4bool	chargeOkay () const
G4bool	strangeOkay () const
G4bool	okay () const
G4double	deltaE () const
G4double	relativeE () const
G4double	deltaKE () const
G4double	relativeKE () const
G4double	deltaP () const
G4double	relativeP () const
G4LorentzVector	deltaLV () const
G4int	deltaB () const
G4int	deltaQ () const
G4int	deltaS () const
Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const char *name, G4int verbose=0)
virtual	~G4VCascadeCollider ()
virtual void	setVerboseLevel (G4int verbose=0)

Static Public Attributes
static const G4double	tolerance = 1e-6

Protected Member Functions
G4double	ekin (const G4LorentzVector &p) const
Protected Member Functions inherited from G4VCascadeCollider
virtual void	setName (const char *name)

Additional Inherited Members
Protected Attributes inherited from G4VCascadeCollider
const char *	theName
G4int	verboseLevel

Detailed Description

Definition at line 64 of file G4CascadeCheckBalance.hh.

Constructor & Destructor Documentation

G4CascadeCheckBalance::G4CascadeCheckBalance ( const char *  owner = "G4CascadeCheckBalance" )

explicit

Definition at line 68 of file G4CascadeCheckBalance.cc.

  : G4VCascadeCollider(owner), relativeLimit(G4CascadeCheckBalance::tolerance),
    absoluteLimit(G4CascadeCheckBalance::tolerance), initialBaryon(0),
    finalBaryon(0), initialCharge(0), finalCharge(0), initialStrange(0),
    finalStrange(0) {}

G4CascadeCheckBalance::G4CascadeCheckBalance	(	G4double	relative,
		G4double	absolute,
		const char *	owner = "G4CascadeCheckBalance"
	)

Definition at line 74 of file G4CascadeCheckBalance.cc.

  : G4VCascadeCollider(owner), relativeLimit(relative),
    absoluteLimit(absolute), initialBaryon(0), finalBaryon(0),
    initialCharge(0), finalCharge(0), initialStrange(0),
    finalStrange(0) {}

virtual G4CascadeCheckBalance::~G4CascadeCheckBalance ( )

inlinevirtual

Definition at line 72 of file G4CascadeCheckBalance.hh.

{};

Member Function Documentation

G4bool G4CascadeCheckBalance::baryonOkay

(

)

const

Definition at line 296 of file G4CascadeCheckBalance.cc.

References deltaB(), G4cerr, G4endl, G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::checkFinalResult(), okay(), and G4CascadeInterface::throwNonConservationFailure().

                                               {
  G4bool bokay = (deltaB() == 0);   // Must be perfect!

  if (verboseLevel && !bokay)
    G4cerr << theName << ": Baryon number VIOLATED " << deltaB() << G4endl;

  return bokay;
}

G4bool G4CascadeCheckBalance::chargeOkay

(

)

const

Definition at line 305 of file G4CascadeCheckBalance.cc.

References deltaQ(), G4cerr, G4endl, G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::checkFinalResult(), okay(), and G4CascadeInterface::throwNonConservationFailure().

                                               {
  G4bool qokay = (deltaQ() == 0);   // Must be perfect!

  if (verboseLevel && !qokay)
    G4cerr << theName << ": Charge conservation VIOLATED " << deltaQ()
       << G4endl;

  return qokay;
}

void G4CascadeCheckBalance::collide ( G4InuclParticle *  bullet, G4InuclParticle *  target, G4CollisionOutput &  output  )

virtual

Implements G4VCascadeCollider.

Definition at line 85 of file G4CascadeCheckBalance.cc.

References G4InuclElementaryParticle::baryon(), CLHEP::HepLorentzVector::e(), G4cout, G4endl, G4InuclNuclei::getA(), G4InuclParticle::getCharge(), G4InuclParticle::getMomentum(), G4InuclElementaryParticle::getStrangeness(), G4CollisionOutput::getTotalBaryonNumber(), G4CollisionOutput::getTotalCharge(), G4CollisionOutput::getTotalOutputMomentum(), G4CollisionOutput::getTotalStrangeness(), CLHEP::HepLorentzVector::px(), CLHEP::HepLorentzVector::py(), CLHEP::HepLorentzVector::pz(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::ApplyYourself(), G4CascadeInterface::checkFinalResult(), G4CascadeRecoilMaker::collide(), collide(), G4NucleiModel::generateParticleFate(), G4CascadeInterface::Propagate(), G4CascadeDeexciteBase::validateOutput(), and G4CascadeColliderBase::validateOutput().

                                               {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName << ")::collide"
       << G4endl;

  initial *= 0.;    // Fast reset; some colliders only have one pointer
  if (bullet) initial += bullet->getMomentum();
  if (target) initial += target->getMomentum();

  // Baryon number, charge and strangeness must be computed "by hand"
  initialCharge = 0;
  if (bullet) initialCharge += G4int(bullet->getCharge());
  if (target) initialCharge += G4int(target->getCharge());

  G4InuclElementaryParticle* pbullet =
    dynamic_cast<G4InuclElementaryParticle*>(bullet);
  G4InuclElementaryParticle* ptarget =
    dynamic_cast<G4InuclElementaryParticle*>(target);

  G4InuclNuclei* nbullet = dynamic_cast<G4InuclNuclei*>(bullet);
  G4InuclNuclei* ntarget = dynamic_cast<G4InuclNuclei*>(target);

  initialBaryon =
    ((pbullet ? pbullet->baryon() : nbullet ? nbullet->getA() : 0) +
     (ptarget ? ptarget->baryon() : ntarget ? ntarget->getA() : 0) );

  // NOTE:  Currently we ignore possibility of hypernucleus target
  initialStrange = 0;
  if (pbullet) initialStrange += pbullet->getStrangeness();
  if (ptarget) initialStrange += ptarget->getStrangeness();

  // Final state totals are computed for us
  final = output.getTotalOutputMomentum();
  finalBaryon = output.getTotalBaryonNumber();
  finalCharge = output.getTotalCharge();
  finalStrange = output.getTotalStrangeness();

  // Report results
  if (verboseLevel) {
    G4cout << " initial px " << initial.px() << " py " << initial.py()
       << " pz " << initial.pz() << " E " << initial.e()
       << " baryon " << initialBaryon << " charge " << initialCharge
       << " strange " << initialStrange << G4endl
       << "   final px " << final.px() << " py " << final.py()
       << " pz " << final.pz() << " E " << final.e()
       << " baryon " << finalBaryon << " charge " << finalCharge
       << " strange " << finalStrange << G4endl;
  }
}

void G4CascadeCheckBalance::collide	(	const G4Fragment &	fragment,
		G4CollisionOutput &	output
	)

Definition at line 138 of file G4CascadeCheckBalance.cc.

References CLHEP::HepLorentzVector::e(), G4cout, G4endl, G4Fragment::GetA(), G4Fragment::GetMomentum(), G4CollisionOutput::getTotalBaryonNumber(), G4CollisionOutput::getTotalCharge(), G4CollisionOutput::getTotalOutputMomentum(), G4CollisionOutput::getTotalStrangeness(), G4Fragment::GetZ(), CLHEP::HepLorentzVector::px(), CLHEP::HepLorentzVector::py(), CLHEP::HepLorentzVector::pz(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                               {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName << ")::collide(<FRAG>)"
       << G4endl;

  // Copy initial state directly from fragment (no bullet/target sums)
  initial = fragment.GetMomentum();
  initialCharge = G4int(fragment.GetZ());
  initialBaryon = G4int(fragment.GetA());
  initialStrange = 0;               // No hypernuclei at present

  // Final state totals are computed for us
  final = output.getTotalOutputMomentum();
  finalBaryon = output.getTotalBaryonNumber();
  finalCharge = output.getTotalCharge();
  finalStrange = output.getTotalStrangeness();

  // Report results
  if (verboseLevel) {
    G4cout << " initial px " << initial.px() << " py " << initial.py()
       << " pz " << initial.pz() << " E " << initial.e()
       << " baryon " << initialBaryon << " charge " << initialCharge
       << " strange " << initialStrange << G4endl
       << "   final px " << final.px() << " py " << final.py()
       << " pz " << final.pz() << " E " << final.e()
       << " baryon " << finalBaryon << " charge " << finalCharge
       << " strange " << finalStrange << G4endl;
  }
}

void G4CascadeCheckBalance::collide	(	G4InuclParticle *	bullet,
		G4InuclParticle *	target,
		const std::vector< G4InuclElementaryParticle > &	particles
	)

Definition at line 172 of file G4CascadeCheckBalance.cc.

References G4CollisionOutput::addOutgoingParticles(), collide(), G4cout, G4endl, G4CollisionOutput::reset(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                           {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName << ")::collide(<vector>)"
       << G4endl;

  tempOutput.reset();           // Buffer for processing
  tempOutput.addOutgoingParticles(particles);
  collide(bullet, target, tempOutput);
}

void G4CascadeCheckBalance::collide	(	G4InuclParticle *	bullet,
		G4InuclParticle *	target,
		const std::vector< G4CascadParticle > &	particles
	)

Definition at line 212 of file G4CascadeCheckBalance.cc.

References G4CollisionOutput::addOutgoingParticles(), collide(), G4cout, G4endl, G4CollisionOutput::reset(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                          {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName
       << ")::collide(<cparticles>)" << G4endl;

  tempOutput.reset();           // Buffer for processing
  tempOutput.addOutgoingParticles(particles);
  collide(bullet, target, tempOutput);
}

void G4CascadeCheckBalance::collide	(	G4InuclParticle *	bullet,
		G4InuclParticle *	target,
		G4CollisionOutput &	output,
		const std::vector< G4CascadParticle > &	cparticles
	)

Definition at line 227 of file G4CascadeCheckBalance.cc.

References G4CollisionOutput::add(), G4CollisionOutput::addOutgoingParticles(), collide(), G4cout, G4endl, G4CollisionOutput::reset(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                        {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName
       << ")::collide(<EP>,<CP>)" << G4endl;

  tempOutput.reset();           // Buffer for processing
  tempOutput.add(output);
  tempOutput.addOutgoingParticles(cparticles);
  collide(bullet, target, tempOutput);
}

void G4CascadeCheckBalance::collide	(	const G4Fragment &	target,
		const std::vector< G4InuclElementaryParticle > &	particles
	)

Definition at line 184 of file G4CascadeCheckBalance.cc.

References G4CollisionOutput::addOutgoingParticles(), collide(), G4cout, G4endl, G4CollisionOutput::reset(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                           {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName
       << ")::collide(<FRAG>,<vector>)" << G4endl;

  tempOutput.reset();           // Buffer for processing
  tempOutput.addOutgoingParticles(particles);
  collide(target, tempOutput);
}

void G4CascadeCheckBalance::collide	(	const G4Fragment &	target,
		const std::vector< G4InuclNuclei > &	fragments
	)

Definition at line 198 of file G4CascadeCheckBalance.cc.

References G4CollisionOutput::addOutgoingNuclei(), collide(), G4cout, G4endl, G4CollisionOutput::reset(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                               {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName << ")::collide(<vector>)"
       << G4endl;

  tempOutput.reset();           // Buffer for processing
  tempOutput.addOutgoingNuclei(fragments);
  collide(target, tempOutput);
}

G4int G4CascadeCheckBalance::deltaB ( ) const

inline

Definition at line 147 of file G4CascadeCheckBalance.hh.

Referenced by baryonOkay(), G4CascadeInterface::checkFinalResult(), G4CascadeRecoilMaker::fillRecoil(), and G4CascadeInterface::throwNonConservationFailure().

{ return (finalBaryon - initialBaryon); }

G4double G4CascadeCheckBalance::deltaE ( ) const

inline

Definition at line 126 of file G4CascadeCheckBalance.hh.

References CLHEP::HepLorentzVector::e().

Referenced by G4CascadeInterface::checkFinalResult(), energyOkay(), relativeE(), and G4CascadeInterface::throwNonConservationFailure().

{ return (final.e() - initial.e()); }

G4double G4CascadeCheckBalance::deltaKE ( ) const

inline

Definition at line 132 of file G4CascadeCheckBalance.hh.

References ekin().

Referenced by G4CascadeInterface::checkFinalResult(), ekinOkay(), and relativeKE().

{ return (ekin(final) - ekin(initial)); }

G4LorentzVector G4CascadeCheckBalance::deltaLV ( ) const

inline

Definition at line 144 of file G4CascadeCheckBalance.hh.

Referenced by deltaP(), and G4CascadeRecoilMaker::fillRecoil().

{ return final - initial; }

G4double G4CascadeCheckBalance::deltaP ( ) const

inline

Definition at line 138 of file G4CascadeCheckBalance.hh.

References deltaLV(), and CLHEP::HepLorentzVector::rho().

Referenced by momentumOkay(), relativeP(), and G4CascadeInterface::throwNonConservationFailure().

{ return deltaLV().rho(); }

G4int G4CascadeCheckBalance::deltaQ ( ) const

inline

Definition at line 148 of file G4CascadeCheckBalance.hh.

Referenced by chargeOkay(), G4CascadeInterface::checkFinalResult(), G4CascadeRecoilMaker::fillRecoil(), and G4CascadeInterface::throwNonConservationFailure().

{ return (finalCharge - initialCharge); }

G4int G4CascadeCheckBalance::deltaS ( ) const

inline

Definition at line 149 of file G4CascadeCheckBalance.hh.

Referenced by strangeOkay().

{ return (finalStrange- initialStrange); }

G4double G4CascadeCheckBalance::ekin ( const G4LorentzVector &  p ) const

inlineprotected

Definition at line 153 of file G4CascadeCheckBalance.hh.

References CLHEP::HepLorentzVector::e(), and CLHEP::HepLorentzVector::m().

Referenced by deltaKE(), dsihae(), elhain(), relativeKE(), scn4ba(), sihnab(), and tsamcs().

{ return (p.e() - p.m()); }

G4bool G4CascadeCheckBalance::ekinOkay

(

)

const

Definition at line 261 of file G4CascadeCheckBalance.cc.

References deltaKE(), G4cerr, G4cout, G4endl, relativeKE(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                             {
  G4bool relokay = (std::abs(relativeKE()) < relativeLimit);
  G4bool absokay = (std::abs(deltaKE()) < absoluteLimit);

  if (verboseLevel && !(relokay || absokay)) {
    G4cerr << theName << ": Kinetic energy balance: relative "
       << relativeKE() << (relokay ? " conserved" : " VIOLATED")
       << " absolute " << deltaKE()
       << (absokay ? " conserved" : " VIOLATED") << G4endl;
  } else if (verboseLevel > 1) {
    G4cout << theName << ": Kinetic energy balance: relative "
       << relativeKE() << " conserved absolute " << deltaKE()
       << " conserved" << G4endl;
  }

  return (relokay && absokay);
}

G4bool G4CascadeCheckBalance::energyOkay

(

)

const

Definition at line 244 of file G4CascadeCheckBalance.cc.

References deltaE(), G4cerr, G4cout, G4endl, relativeE(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by okay(), and G4CascadeInterface::throwNonConservationFailure().

                                               {
  G4bool relokay = (std::abs(relativeE()) < relativeLimit);
  G4bool absokay = (std::abs(deltaE()) < absoluteLimit);

  if (verboseLevel && !(relokay || absokay)) {
    G4cerr << theName << ": Energy conservation: relative " << relativeE()
       << (relokay ? " conserved" : " VIOLATED")
       << " absolute " << deltaE()
       << (absokay ? " conserved" : " VIOLATED") << G4endl;
  } else if (verboseLevel > 1) {
    G4cout << theName << ": Energy conservation: relative " << relativeE()
       << " conserved absolute " << deltaE() << " conserved" << G4endl;
  }

  return (relokay && absokay);
}

G4bool G4CascadeCheckBalance::momentumOkay

(

)

const

Definition at line 279 of file G4CascadeCheckBalance.cc.

References deltaP(), G4cerr, G4cout, G4endl, relativeP(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by okay(), and G4CascadeInterface::throwNonConservationFailure().

                                                 {
  G4bool relokay = (std::abs(relativeP()) < relativeLimit);
  G4bool absokay = (std::abs(deltaP()) < absoluteLimit);

  if (verboseLevel && !(relokay || absokay)) {
    G4cerr << theName << ": Momentum conservation: relative " << relativeP()
       << (relokay ? " conserved" : " VIOLATED")
       << " absolute " << deltaP()
       << (absokay ? " conserved" : " VIOLATED") << G4endl;
  } else if (verboseLevel > 1) {
    G4cout << theName << ": Momentum conservation: relative " << relativeP()
       << " conserved absolute " << deltaP() << " conserved" << G4endl;
  }

  return (relokay && absokay);
}

G4bool G4CascadeCheckBalance::okay ( ) const

inline

Definition at line 121 of file G4CascadeCheckBalance.hh.

References baryonOkay(), chargeOkay(), energyOkay(), and momentumOkay().

Referenced by G4CascadeInterface::ApplyYourself(), G4NucleiModel::generateParticleFate(), G4CascadeInterface::Propagate(), G4CascadeInterface::retryInelasticNucleus(), G4CascadeDeexciteBase::validateOutput(), and G4CascadeColliderBase::validateOutput().

                      { return (energyOkay() && momentumOkay() &&
                baryonOkay() && chargeOkay()); }

G4double G4CascadeCheckBalance::relativeE ( ) const

inline

Definition at line 127 of file G4CascadeCheckBalance.hh.

References deltaE(), CLHEP::HepLorentzVector::e(), and tolerance.

Referenced by energyOkay(), and G4CascadeInterface::throwNonConservationFailure().

                             {
    return ( (std::abs(deltaE())<tolerance) ? 0. : 
         (initial.e()<tolerance) ? 1. : deltaE()/initial.e() );
  }

G4double G4CascadeCheckBalance::relativeKE ( ) const

inline

Definition at line 133 of file G4CascadeCheckBalance.hh.

References deltaKE(), ekin(), and tolerance.

Referenced by ekinOkay().

                              {
    return ( (std::abs(deltaKE())<tolerance) ? 0. : 
         (ekin(initial)<tolerance) ? 1. : deltaKE()/ekin(initial) );
  }

G4double G4CascadeCheckBalance::relativeP ( ) const

inline

Definition at line 139 of file G4CascadeCheckBalance.hh.

References deltaP(), CLHEP::HepLorentzVector::rho(), and tolerance.

Referenced by momentumOkay(), and G4CascadeInterface::throwNonConservationFailure().

                             {
    return ( (std::abs(deltaP())<tolerance) ? 0. : 
         (initial.rho()<tolerance) ? 1. : deltaP()/initial.rho() );
  }

void G4CascadeCheckBalance::setAbsoluteLimit ( G4double  limit )

inline

Definition at line 82 of file G4CascadeCheckBalance.hh.

Referenced by setLimits().

{ absoluteLimit = limit; }

void G4CascadeCheckBalance::setLimits ( G4double  relative, G4double  absolute  )

inline

Definition at line 76 of file G4CascadeCheckBalance.hh.

References setAbsoluteLimit(), and setRelativeLimit().

Referenced by G4CascadeInterface::G4CascadeInterface().

                                                       {
    setRelativeLimit(relative);
    setAbsoluteLimit(absolute);
  }

void G4CascadeCheckBalance::setOwner ( const char *  owner )

inline

Definition at line 74 of file G4CascadeCheckBalance.hh.

References G4VCascadeCollider::setName().

{ setName(owner); }

void G4CascadeCheckBalance::setRelativeLimit ( G4double  limit )

inline

Definition at line 81 of file G4CascadeCheckBalance.hh.

Referenced by setLimits().

{ relativeLimit = limit; }

G4bool G4CascadeCheckBalance::strangeOkay

(

)

const

Definition at line 315 of file G4CascadeCheckBalance.cc.

References deltaS(), G4cerr, G4endl, G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

                                                {
  G4bool sokay = (deltaS() == 0);   // Must be perfect!

  if (verboseLevel && !sokay)
    G4cerr << theName << ": Strangeness conservation VIOLATED " << deltaS()
       << G4endl;

  return sokay;
}

Field Documentation

const G4double G4CascadeCheckBalance::tolerance = 1e-6

static

Definition at line 66 of file G4CascadeCheckBalance.hh.

Referenced by relativeE(), relativeKE(), and relativeP().

---

The documentation for this class was generated from the following files:

• G4CascadeCheckBalance.hh
• G4CascadeCheckBalance.cc