G4CascadeCheckBalance Class Reference

#include <G4CascadeCheckBalance.hh>

Inheritance diagram for G4CascadeCheckBalance:

Public Member Functions
	G4CascadeCheckBalance (const G4String &owner="G4CascadeCheckBalance")
	G4CascadeCheckBalance (G4double relative, G4double absolute, const G4String &owner="G4CascadeCheckBalance")
virtual	~G4CascadeCheckBalance ()
void	setOwner (const G4String &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 G4String &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 G4String &name)

Additional Inherited Members
Protected Attributes inherited from G4VCascadeCollider
G4String	theName
G4int	verboseLevel

Detailed Description

Definition at line 64 of file G4CascadeCheckBalance.hh.

Constructor & Destructor Documentation

G4CascadeCheckBalance() [1/2]

G4CascadeCheckBalance::G4CascadeCheckBalance ( const G4String &  owner = "G4CascadeCheckBalance" )

explicit

Definition at line 71 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 G4String &	owner = "G4CascadeCheckBalance"
	)

Definition at line 77 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 72 of file G4CascadeCheckBalance.hh.

{};

Member Function Documentation

baryonOkay()

G4bool G4CascadeCheckBalance::baryonOkay

(

)

const

Definition at line 330 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 338 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/7]

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

Definition at line 156 of file G4CascadeCheckBalance.cc.

                                               {
  if (verboseLevel)
    G4cout << " >>> G4CascadeCheckBalance(" << theName << ")::collide(<FRAG>)"
       << G4endl;
 
  // Copy initial state directly from fragment (no bullet/target sums)
  initial = fragment.GetMomentum()/GeV;         // G4Fragment returns MeV
  initialCharge = fragment.GetZ_asInt();
  initialBaryon = fragment.GetA_asInt();
  initialStrange = 0;                           // No hypernuclei at present
 
  // Final state totals are computed for us
  final = output.getTotalOutputMomentum();
  
  // Remove electron masses when internal conversion occurs
  G4double elMass = 0.;
  std::vector<G4InuclElementaryParticle>& outParts = 
    output.getOutgoingParticles();
  G4int nelec = 0;
  for (G4int i = 0; i < output.numberOfOutgoingParticles(); i++) {
    if (outParts[i].getDefinition() == G4Electron::Electron() ) {
      elMass += outParts[i].getDefinition()->GetPDGMass();
      ++nelec;
    }
  }
  if(nelec > 0) {
    initial += G4LorentzVector(0.,0.,0.,elMass/GeV);
    initialCharge -= nelec;
  }
 
  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;
  }
}

collide() [2/7]

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

Definition at line 219 of file G4CascadeCheckBalance.cc.

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

collide() [3/7]

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

Definition at line 232 of file G4CascadeCheckBalance.cc.

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

collide() [4/7]

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

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

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

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

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

virtual

Implements G4VCascadeCollider.

Definition at line 88 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();
 
  G4int nelec = 0;
  G4double elMass = 0.;
  std::vector<G4InuclElementaryParticle>& outParts = 
    output.getOutgoingParticles();
  for (G4int i = 0; i < output.numberOfOutgoingParticles(); i++) {
    if (outParts[i].getDefinition() == G4Electron::Electron() ) {
      elMass += outParts[i].getDefinition()->GetPDGMass();
      ++nelec;
    }
  }
  if(nelec > 0) {
    initial += G4LorentzVector(0.,0.,0.,elMass/GeV);
    initialCharge -= nelec;
  }
 
  // 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(), G4CascadeDeexciteBase::validateOutput(), and G4CascadeColliderBase::validateOutput().

collide() [7/7]

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

Definition at line 259 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 147 of file G4CascadeCheckBalance.hh.

{ return (finalBaryon - initialBaryon); }

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

deltaE()

G4double G4CascadeCheckBalance::deltaE ( ) const

inline

Definition at line 126 of file G4CascadeCheckBalance.hh.

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

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

deltaKE()

G4double G4CascadeCheckBalance::deltaKE ( ) const

inline

Definition at line 132 of file G4CascadeCheckBalance.hh.

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

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

deltaLV()

G4LorentzVector G4CascadeCheckBalance::deltaLV ( ) const

inline

Definition at line 144 of file G4CascadeCheckBalance.hh.

{ return final - initial; }

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

deltaP()

G4double G4CascadeCheckBalance::deltaP ( ) const

inline

Definition at line 138 of file G4CascadeCheckBalance.hh.

{ return deltaLV().rho(); }

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

deltaQ()

G4int G4CascadeCheckBalance::deltaQ ( ) const

inline

Definition at line 148 of file G4CascadeCheckBalance.hh.

{ return (finalCharge - initialCharge); }

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

deltaS()

G4int G4CascadeCheckBalance::deltaS ( ) const

inline

Definition at line 149 of file G4CascadeCheckBalance.hh.

{ return (finalStrange- initialStrange); }

Referenced by strangeOkay().

ekin()

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

inlineprotected

Definition at line 153 of file G4CascadeCheckBalance.hh.

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

Referenced by deltaKE(), and relativeKE().

ekinOkay()

G4bool G4CascadeCheckBalance::ekinOkay

(

)

const

Definition at line 293 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 275 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 312 of file G4CascadeCheckBalance.cc.

                                                 {
  G4bool relokay = (std::abs(relativeP()) < 10.*relativeLimit);
  G4bool absokay = (std::abs(deltaP()) < 10.*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 121 of file G4CascadeCheckBalance.hh.

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

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

relativeE()

G4double G4CascadeCheckBalance::relativeE ( ) const

inline

Definition at line 127 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 133 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 139 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 82 of file G4CascadeCheckBalance.hh.

{ absoluteLimit = limit; }

Referenced by setLimits().

setLimits()

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

inline

Definition at line 76 of file G4CascadeCheckBalance.hh.

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

Referenced by G4CascadeInterface::G4CascadeInterface().

setOwner()

void G4CascadeCheckBalance::setOwner ( const G4String &  owner )

inline

Definition at line 74 of file G4CascadeCheckBalance.hh.

{ setName(owner); }

setRelativeLimit()

void G4CascadeCheckBalance::setRelativeLimit ( G4double  limit )

inline

Definition at line 81 of file G4CascadeCheckBalance.hh.

{ relativeLimit = limit; }

Referenced by setLimits().

strangeOkay()

G4bool G4CascadeCheckBalance::strangeOkay

(

)

const

Definition at line 347 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 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