G4ReactionProduct.cc File Reference

#include "G4ReactionProduct.hh"

Go to the source code of this file.

Functions
G4Allocator< G4ReactionProduct > *&	aRPAllocator ()
G4ReactionProduct	operator+ (const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4ReactionProduct	operator- (const G4ReactionProduct &p1, const G4ReactionProduct &p2)

Function Documentation

aRPAllocator()

G4Allocator< G4ReactionProduct > *& aRPAllocator

(

)

Definition at line 33 of file G4ReactionProduct.cc.

{
    G4ThreadLocalStatic G4Allocator<G4ReactionProduct>* _instance = nullptr;
    return _instance;
}

References G4ThreadLocalStatic.

Referenced by G4ReactionProduct::operator delete(), and G4ReactionProduct::operator new().

operator+()

G4ReactionProduct operator+	(	const G4ReactionProduct &	p1,
		const G4ReactionProduct &	p2
	)

Definition at line 245 of file G4ReactionProduct.cc.

  {
    G4double totEnergy = p1.totalEnergy + p2.totalEnergy;
    G4double x = p1.momentum.x() + p2.momentum.x();
    G4double y = p1.momentum.y() + p2.momentum.y();
    G4double z = p1.momentum.z() + p2.momentum.z();
    G4double newMass = totEnergy*totEnergy - ( x*x + y*y + z*z );
    if( newMass < 0.0 )
      newMass = -1. * std::sqrt( -newMass );
    else
      newMass = std::sqrt( newMass );
    G4ReactionProduct result;
    result.SetMass( newMass );
    result.SetMomentum( x, y, z );
    result.SetTotalEnergy( totEnergy );
    result.SetPositionInNucleus( 0.0, 0.0, 0.0 );
    result.SetFormationTime(0.0);
    result.HasInitialStateParton(false);
    return result;
  }

operator-()

G4ReactionProduct operator-	(	const G4ReactionProduct &	p1,
		const G4ReactionProduct &	p2
	)

Definition at line 267 of file G4ReactionProduct.cc.

  {
    G4double totEnergy = p1.totalEnergy - p2.totalEnergy;
    G4double x = p1.momentum.x() - p2.momentum.x();
    G4double y = p1.momentum.y() - p2.momentum.y();
    G4double z = p1.momentum.z() - p2.momentum.z();
    G4double newMass = totEnergy*totEnergy - ( x*x + y*y + z*z );
    if( newMass < 0.0 )
      newMass = -1. * std::sqrt( -newMass );
    else
      newMass = std::sqrt( newMass );
    G4ReactionProduct result;
    result.SetMass( newMass );
    result.SetMomentum( x, y, z );
    result.SetTotalEnergy( totEnergy );
    result.SetPositionInNucleus( 0.0, 0.0, 0.0 );
    result.SetFormationTime(0.0);
    result.HasInitialStateParton(false);
    return result;
  }