#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 (G4InuclParticle bullet, G4InuclParticle target, const std::vector< G4InuclElementaryParticle > &particles)

void	collide (G4InuclParticle bullet, G4InuclParticle target, const std::vector< G4InuclNuclei > &fragments)

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

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

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	collide (G4InuclParticle bullet, G4InuclParticle target, G4CollisionOutput &output)=0

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 62 of file G4CascadeCheckBalance.hh.

Constructor & Destructor Documentation

◆ G4CascadeCheckBalance() [1/2]

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

explicit

Definition at line 67 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() [2/2]

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

Definition at line 73 of file G4CascadeCheckBalance.cc.

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

◆ ~G4CascadeCheckBalance()

virtual G4CascadeCheckBalance::~G4CascadeCheckBalance ( )

inlinevirtual

Definition at line 70 of file G4CascadeCheckBalance.hh.

70{};

Member Function Documentation

◆ baryonOkay()

G4bool G4CascadeCheckBalance::baryonOkay ( ) const

Definition at line 252 of file G4CascadeCheckBalance.cc.

                                               {
  G4bool bokay = (deltaB() == 0);   // Must be perfect!
 
  if (verboseLevel && !bokay)
    G4cerr << theName << ": Baryon number VIOLATED " << deltaB() << G4endl;
 
  return bokay;
}

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

◆ chargeOkay()

G4bool G4CascadeCheckBalance::chargeOkay ( ) const

Definition at line 261 of file G4CascadeCheckBalance.cc.

                                               {
  G4bool qokay = (deltaQ() == 0);   // Must be perfect!
 
  if (verboseLevel && !qokay)
    G4cerr << theName << ": Charge conservation VIOLATED " << deltaQ()
       << G4endl;
 
  return qokay;
}

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

◆ collide() [1/5]

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

Definition at line 168 of file G4CascadeCheckBalance.cc.

                                                          {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName
       << ")::collide(<cparticles>)" << G4endl;
 
  tempOutput.reset();           // Buffer for processing
  tempOutput.addOutgoingParticles(particles);
  collide(bullet, target, tempOutput);
}

◆ collide() [2/5]

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

Definition at line 138 of file G4CascadeCheckBalance.cc.

                                                           {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName << ")::collide(<vector>)"
       << G4endl;
 
  tempOutput.reset();           // Buffer for processing
  tempOutput.addOutgoingParticles(particles);
  collide(bullet, target, tempOutput);
}

◆ collide() [3/5]

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

Definition at line 153 of file G4CascadeCheckBalance.cc.

                                               {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName << ")::collide(<vector>)"
       << G4endl;
 
  tempOutput.reset();           // Buffer for processing
  tempOutput.addOutgoingNuclei(fragments);
  collide(bullet, target, tempOutput);
}

◆ collide() [4/5]

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

virtual

Implements G4VCascadeCollider.

Definition at line 84 of file G4CascadeCheckBalance.cc.

                                               {
  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;
  }
}

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

◆ collide() [5/5]

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

Definition at line 183 of file G4CascadeCheckBalance.cc.

                                                        {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName
       << ")::collide(<EP>,<CP>)" << G4endl;
 
  tempOutput.reset();           // Buffer for processing
  tempOutput.add(output);
  tempOutput.addOutgoingParticles(cparticles);
  collide(bullet, target, tempOutput);
}

◆ deltaB()

G4int G4CascadeCheckBalance::deltaB ( ) const

inline

Definition at line 138 of file G4CascadeCheckBalance.hh.

138{ return (finalBaryon - initialBaryon); }

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

◆ deltaE()

G4double G4CascadeCheckBalance::deltaE ( ) const

inline

Definition at line 117 of file G4CascadeCheckBalance.hh.

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

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

◆ deltaKE()

G4double G4CascadeCheckBalance::deltaKE ( ) const

inline

Definition at line 123 of file G4CascadeCheckBalance.hh.

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

G4CascadeCheckBalance::ekin

G4double ekin(const G4LorentzVector &p) const

Definition: G4CascadeCheckBalance.hh:144

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

◆ deltaLV()

G4LorentzVector G4CascadeCheckBalance::deltaLV ( ) const

inline

Definition at line 135 of file G4CascadeCheckBalance.hh.

135{ return final - initial; }

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

◆ deltaP()

G4double G4CascadeCheckBalance::deltaP ( ) const

inline

Definition at line 129 of file G4CascadeCheckBalance.hh.

129{ return deltaLV().rho(); }

CLHEP::HepLorentzVector::rho

double rho() const

G4CascadeCheckBalance::deltaLV

G4LorentzVector deltaLV() const

Definition: G4CascadeCheckBalance.hh:135

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

◆ deltaQ()

G4int G4CascadeCheckBalance::deltaQ ( ) const

inline

Definition at line 139 of file G4CascadeCheckBalance.hh.

139{ return (finalCharge - initialCharge); }

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

◆ deltaS()

G4int G4CascadeCheckBalance::deltaS ( ) const

inline

Definition at line 140 of file G4CascadeCheckBalance.hh.

140{ return (finalStrange- initialStrange); }

Referenced by strangeOkay().

◆ ekin()

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

inlineprotected

Definition at line 144 of file G4CascadeCheckBalance.hh.

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

CLHEP::HepLorentzVector::m

double m() const

Referenced by deltaKE(), and relativeKE().

◆ ekinOkay()

G4bool G4CascadeCheckBalance::ekinOkay ( ) const

Definition at line 217 of file G4CascadeCheckBalance.cc.

                                             {
  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);
}

◆ energyOkay()

G4bool G4CascadeCheckBalance::energyOkay ( ) const

Definition at line 200 of file G4CascadeCheckBalance.cc.

                                               {
  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);
}

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

◆ momentumOkay()

G4bool G4CascadeCheckBalance::momentumOkay ( ) const

Definition at line 235 of file G4CascadeCheckBalance.cc.

                                                 {
  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);
}

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

◆ okay()

G4bool G4CascadeCheckBalance::okay ( ) const

inline

Definition at line 112 of file G4CascadeCheckBalance.hh.

112 { return (energyOkay() && momentumOkay() &&

113 baryonOkay() && chargeOkay()); }

G4CascadeCheckBalance::baryonOkay

G4bool baryonOkay() const

Definition: G4CascadeCheckBalance.cc:252

G4CascadeCheckBalance::momentumOkay

G4bool momentumOkay() const

Definition: G4CascadeCheckBalance.cc:235

G4CascadeCheckBalance::chargeOkay

G4bool chargeOkay() const

Definition: G4CascadeCheckBalance.cc:261

G4CascadeCheckBalance::energyOkay

G4bool energyOkay() const

Definition: G4CascadeCheckBalance.cc:200

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

◆ relativeE()

G4double G4CascadeCheckBalance::relativeE ( ) const

inline

Definition at line 118 of file G4CascadeCheckBalance.hh.

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

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

◆ relativeKE()

G4double G4CascadeCheckBalance::relativeKE ( ) const

inline

Definition at line 124 of file G4CascadeCheckBalance.hh.

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

Referenced by ekinOkay().

◆ relativeP()

G4double G4CascadeCheckBalance::relativeP ( ) const

inline

Definition at line 130 of file G4CascadeCheckBalance.hh.

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

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

◆ setAbsoluteLimit()

void G4CascadeCheckBalance::setAbsoluteLimit ( G4double limit )

inline

Definition at line 80 of file G4CascadeCheckBalance.hh.

80{ absoluteLimit = limit; }

Referenced by setLimits().

◆ setLimits()

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

inline

Definition at line 74 of file G4CascadeCheckBalance.hh.

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

Referenced by G4CascadeInterface::G4CascadeInterface().

◆ setOwner()

void G4CascadeCheckBalance::setOwner ( const char * owner )

inline

Definition at line 72 of file G4CascadeCheckBalance.hh.

72{ setName(owner); }

G4VCascadeCollider::setName

virtual void setName(const char *name)

Definition: G4VCascadeCollider.hh:54

◆ setRelativeLimit()

void G4CascadeCheckBalance::setRelativeLimit ( G4double limit )

inline

Definition at line 79 of file G4CascadeCheckBalance.hh.

79{ relativeLimit = limit; }

Referenced by setLimits().

◆ strangeOkay()

G4bool G4CascadeCheckBalance::strangeOkay ( ) const

Definition at line 271 of file G4CascadeCheckBalance.cc.

                                                {
  G4bool sokay = (deltaS() == 0);   // Must be perfect!
 
  if (verboseLevel && !sokay)
    G4cerr << theName << ": Strangeness conservation VIOLATED " << deltaS()
       << G4endl;
 
  return sokay;
}

Member Data Documentation

◆ tolerance

const G4double G4CascadeCheckBalance::tolerance = 1e-6

static

Definition at line 64 of file G4CascadeCheckBalance.hh.

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

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

Public Member Functions

Static Public Attributes

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4CascadeCheckBalance() [1/2]

◆ G4CascadeCheckBalance() [2/2]

◆ ~G4CascadeCheckBalance()

Member Function Documentation

◆ baryonOkay()

◆ chargeOkay()

◆ collide() [1/5]

◆ collide() [2/5]

◆ collide() [3/5]

◆ collide() [4/5]

◆ collide() [5/5]

◆ deltaB()

◆ deltaE()

◆ deltaKE()

◆ deltaLV()

◆ deltaP()

◆ deltaQ()

◆ deltaS()

◆ ekin()

◆ ekinOkay()

◆ energyOkay()

◆ momentumOkay()

◆ okay()

◆ relativeE()

◆ relativeKE()

◆ relativeP()

◆ setAbsoluteLimit()

◆ setLimits()

◆ setOwner()

◆ setRelativeLimit()

◆ strangeOkay()

Member Data Documentation

◆ tolerance