G4CascadeCheckBalance Class Reference

#include <G4CascadeCheckBalance.hh>

Inheritance diagram for G4CascadeCheckBalance:

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

Constructor & Destructor Documentation

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

explicit

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

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

{};

Member Function Documentation

G4bool G4CascadeCheckBalance::baryonOkay

(

)

const

Definition at line 298 of file G4CascadeCheckBalance.cc.

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

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

  return bokay;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4CascadeCheckBalance::chargeOkay

(

)

const

Definition at line 307 of file G4CascadeCheckBalance.cc.

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

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

  return qokay;
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Implements G4VCascadeCollider.

Definition at line 87 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;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 140 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();
  initialCharge = fragment.GetZ_asInt();
  initialBaryon = fragment.GetA_asInt();
  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;
  }
}

Here is the call graph for this function:

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

Definition at line 174 of file G4CascadeCheckBalance.cc.

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

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

Here is the call graph for this function:

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

Definition at line 214 of file G4CascadeCheckBalance.cc.

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

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

Here is the call graph for this function:

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

Definition at line 229 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);
}

Here is the call graph for this function:

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

Definition at line 186 of file G4CascadeCheckBalance.cc.

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

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

Here is the call graph for this function:

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

Definition at line 200 of file G4CascadeCheckBalance.cc.

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

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

Here is the call graph for this function:

G4int G4CascadeCheckBalance::deltaB ( ) const

inline

Definition at line 148 of file G4CascadeCheckBalance.hh.

{ return (finalBaryon - initialBaryon); }

Here is the caller graph for this function:

G4double G4CascadeCheckBalance::deltaE ( ) const

inline

Definition at line 127 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4CascadeCheckBalance::deltaKE ( ) const

inline

Definition at line 133 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4LorentzVector G4CascadeCheckBalance::deltaLV ( ) const

inline

Definition at line 145 of file G4CascadeCheckBalance.hh.

{ return final - initial; }

Here is the caller graph for this function:

G4double G4CascadeCheckBalance::deltaP ( ) const

inline

Definition at line 139 of file G4CascadeCheckBalance.hh.

{ return deltaLV().rho(); }

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4CascadeCheckBalance::deltaQ ( ) const

inline

Definition at line 149 of file G4CascadeCheckBalance.hh.

{ return (finalCharge - initialCharge); }

Here is the caller graph for this function:

G4int G4CascadeCheckBalance::deltaS ( ) const

inline

Definition at line 150 of file G4CascadeCheckBalance.hh.

{ return (finalStrange- initialStrange); }

Here is the caller graph for this function:

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

inlineprotected

Definition at line 154 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4CascadeCheckBalance::ekinOkay

(

)

const

Definition at line 263 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);
}

Here is the call graph for this function:

G4bool G4CascadeCheckBalance::energyOkay

(

)

const

Definition at line 246 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4CascadeCheckBalance::momentumOkay

(

)

const

Definition at line 281 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4CascadeCheckBalance::okay ( ) const

inline

Definition at line 122 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4CascadeCheckBalance::relativeE ( ) const

inline

Definition at line 128 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4CascadeCheckBalance::relativeKE ( ) const

inline

Definition at line 134 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4CascadeCheckBalance::relativeP ( ) const

inline

Definition at line 140 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

void G4CascadeCheckBalance::setAbsoluteLimit ( G4double  limit )

inline

Definition at line 83 of file G4CascadeCheckBalance.hh.

{ absoluteLimit = limit; }

Here is the caller graph for this function:

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

inline

Definition at line 77 of file G4CascadeCheckBalance.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

void G4CascadeCheckBalance::setOwner ( const G4String &  owner )

inline

Definition at line 75 of file G4CascadeCheckBalance.hh.

{ setName(owner); }

Here is the call graph for this function:

void G4CascadeCheckBalance::setRelativeLimit ( G4double  limit )

inline

Definition at line 82 of file G4CascadeCheckBalance.hh.

{ relativeLimit = limit; }

Here is the caller graph for this function:

G4bool G4CascadeCheckBalance::strangeOkay

(

)

const

Definition at line 317 of file G4CascadeCheckBalance.cc.

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

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

  return sokay;
}

Here is the call graph for this function:

Member Data Documentation

const G4double G4CascadeCheckBalance::tolerance = 1e-6

static

Definition at line 67 of file G4CascadeCheckBalance.hh.

---

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

• source/geant4.10.03.p02/source/processes/hadronic/models/cascade/cascade/include/G4CascadeCheckBalance.hh
• source/geant4.10.03.p02/source/processes/hadronic/models/cascade/cascade/src/G4CascadeCheckBalance.cc