G4GHEKinematicsVector Class Reference

#include <G4GHEKinematicsVector.hh>

Collaboration diagram for G4GHEKinematicsVector:

Public Member Functions
	G4GHEKinematicsVector ()
	~G4GHEKinematicsVector ()
	G4GHEKinematicsVector (const G4GHEKinematicsVector &p)
G4GHEKinematicsVector &	operator= (const G4GHEKinematicsVector &p)
void	SetMomentum (G4ParticleMomentum mom)
void	SetMomentumAndUpdate (G4ParticleMomentum mom)
const G4ParticleMomentum	GetMomentum () const
void	SetMomentum (G4double x, G4double y, G4double z)
void	SetMomentumAndUpdate (G4double x, G4double y, G4double z)
void	SetMomentum (G4double x, G4double y)
void	SetMomentumAndUpdate (G4double x, G4double y)
void	SetMomentum (G4double z)
void	SetMomentumAndUpdate (G4double z)
void	SetEnergy (G4double e)
void	SetEnergyAndUpdate (G4double e)
void	SetKineticEnergy (G4double ekin)
void	SetKineticEnergyAndUpdate (G4double ekin)
G4double	GetEnergy ()
G4double	GetKineticEnergy ()
void	SetMass (G4double mas)
void	SetMassAndUpdate (G4double mas)
G4double	GetMass ()
void	SetCharge (G4double c)
G4double	GetCharge ()
void	SetTOF (G4double t)
G4double	GetTOF ()
void	SetSide (G4int sid)
G4int	GetSide ()
void	setFlag (G4bool f)
G4bool	getFlag ()
void	SetCode (G4int c)
void	SetParticleDef (G4ParticleDefinition *c)
G4int	GetCode ()
G4ParticleDefinition *	GetParticleDef ()
void	SetZero ()
void	Add (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
void	Sub (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
void	Lor (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
G4double	CosAng (const G4GHEKinematicsVector &p)
G4double	Ang (const G4GHEKinematicsVector &p)
G4double	Dot4 (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
G4double	Impu (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
void	Add3 (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
void	Sub3 (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
void	Cross (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
G4double	Dot (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
void	Smul (const G4GHEKinematicsVector &p, G4double h)
void	SmulAndUpdate (const G4GHEKinematicsVector &p, G4double h)
void	Norz (const G4GHEKinematicsVector &p)
G4double	Length ()
void	Exch (G4GHEKinematicsVector &p1)
void	Defs1 (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2)
void	Defs (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &p2, G4GHEKinematicsVector &my, G4GHEKinematicsVector &mz)
void	Trac (const G4GHEKinematicsVector &p1, const G4GHEKinematicsVector &mx, const G4GHEKinematicsVector &my, const G4GHEKinematicsVector &mz)
void	Print (G4int L)

Public Attributes
G4ParticleMomentum	momentum
G4double	energy
G4double	kineticEnergy
G4double	mass
G4double	charge
G4double	timeOfFlight
G4int	side
G4bool	flag
G4int	code
G4ParticleDefinition *	particleDef

Detailed Description

Definition at line 50 of file G4GHEKinematicsVector.hh.

Constructor & Destructor Documentation

G4GHEKinematicsVector::G4GHEKinematicsVector ( )

inline

Definition at line 54 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX(  0.0 );
     momentum.setY(  0.0 );
     momentum.setZ(  0.0 );
     energy        = 0.0;
     kineticEnergy = 0.0;
     mass          = 0.0;
     charge        = 0.0;
     timeOfFlight  = 0.0;
     side          = 0;
     flag          = false;
     code          = 0;
     particleDef   = NULL;
   }

Here is the call graph for this function:

G4GHEKinematicsVector::~G4GHEKinematicsVector ( )

inline

Definition at line 70 of file G4GHEKinematicsVector.hh.

{}

G4GHEKinematicsVector::G4GHEKinematicsVector ( const G4GHEKinematicsVector &  p )

inline

Definition at line 73 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( p.momentum.x() );
     momentum.setY( p.momentum.y() );
     momentum.setZ( p.momentum.z() );
     energy        = p.energy;
     kineticEnergy = p.kineticEnergy;
     mass          = p.mass;
     charge        = p.charge;
     timeOfFlight  = p.timeOfFlight;
     side          = p.side;
     flag          = p.flag;
     code          = p.code;
     particleDef   = p.particleDef;
   }

Here is the call graph for this function:

Member Function Documentation

void G4GHEKinematicsVector::Add ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 361 of file G4GHEKinematicsVector.hh.

   {
     momentum = p1.momentum + p2.momentum;
     energy = p1.energy + p2.energy;
     G4double b = energy*energy - momentum.mag2();
     if( b < 0 )
       mass = -1. * std::sqrt( -b );
     else
       mass = std::sqrt( b );
     kineticEnergy = std::max(0.,energy - mass);
     charge        = p1.charge + p2.charge;
     code          = p1.code   + p2.code;
     particleDef   = p1.particleDef;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Add3 ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 456 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( p1.momentum.x() + p2.momentum.x());
     momentum.setY( p1.momentum.y() + p2.momentum.y());
     momentum.setZ( p1.momentum.z() + p2.momentum.z());    
     return;
   }

Here is the call graph for this function:

G4double G4GHEKinematicsVector::Ang ( const G4GHEKinematicsVector &  p )

inline

Definition at line 424 of file G4GHEKinematicsVector.hh.

   {
     G4double a = std::sqrt( momentum.mag2() * p.momentum.mag2() );
     if( a != 0.0 ) 
       {
         a = (momentum.x()*p.momentum.x() +
              momentum.y()*p.momentum.y() +
              momentum.z()*p.momentum.z()) / a;
         if( std::fabs(a) > 1.0 ) a<0.0 ? a=-1.0 : a=1.0;
       }
     return std::acos(a);
   }      

Here is the call graph for this function:

G4double G4GHEKinematicsVector::CosAng ( const G4GHEKinematicsVector &  p )

inline

Definition at line 411 of file G4GHEKinematicsVector.hh.

   {
     G4double a = std::sqrt( momentum.mag2() * p.momentum.mag2() );
     if( a != 0.0 ) 
       {
         a = (momentum.x()*p.momentum.x() +
              momentum.y()*p.momentum.y() +
              momentum.z()*p.momentum.z()) / a;
         if( std::fabs(a) > 1.0 ) a<0.0 ? a=-1.0 : a=1.0;
       }
     return a;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Cross ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 474 of file G4GHEKinematicsVector.hh.

   {
     G4double px, py, pz;
     px = p1.momentum.y() * p2.momentum.z() - p1.momentum.z() * p2.momentum.y();
     py = p1.momentum.z() * p2.momentum.x() - p1.momentum.x() * p2.momentum.z();
     pz = p1.momentum.x() * p2.momentum.y() - p1.momentum.y() * p2.momentum.x();
     momentum.setX( px );
     momentum.setY( py );
     momentum.setZ( pz ); 
     return;
   }  

Here is the call graph for this function:

void G4GHEKinematicsVector::Defs ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2, G4GHEKinematicsVector &  my, G4GHEKinematicsVector &  mz  )

inline

Definition at line 608 of file G4GHEKinematicsVector.hh.

   {
     my = p1;
     mz = p2;
     momentum.setX(   my.momentum.y()*mz.momentum.z()
                    - my.momentum.z()*mz.momentum.y());
     momentum.setY(   my.momentum.z()*mz.momentum.x()
                    - my.momentum.x()*mz.momentum.z());
     momentum.setZ(   my.momentum.x()*mz.momentum.y()
                    - my.momentum.y()*mz.momentum.x());
     my.momentum.setX(   mz.momentum.y()*momentum.z()
                       - mz.momentum.z()*momentum.y());
     my.momentum.setY(   mz.momentum.z()*momentum.x()
                       - mz.momentum.x()*momentum.z());
     my.momentum.setZ(   mz.momentum.x()*momentum.y()
                       - mz.momentum.y()*momentum.x());
     G4double pp;
     pp = momentum.mag();
     if (pp > 0.)
        {
          pp = 1./pp; 
          momentum.setX( momentum.x()*pp );
          momentum.setY( momentum.y()*pp );
          momentum.setZ( momentum.z()*pp );
        }
     pp = my.momentum.mag();
     if (pp > 0.)
        {
          pp = 1./pp;
          my.momentum.setX( my.momentum.x()*pp );
          my.momentum.setY( my.momentum.y()*pp );
          my.momentum.setZ( my.momentum.z()*pp );
        }
     pp = mz.momentum.mag();
     if (pp > 0.)
        {
          pp = 1./pp;
          mz.momentum.setX( mz.momentum.x()*pp );
          mz.momentum.setY( mz.momentum.y()*pp );
          mz.momentum.setZ( mz.momentum.z()*pp );
        }
     return; 
   }  

Here is the call graph for this function:

void G4GHEKinematicsVector::Defs1 ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 579 of file G4GHEKinematicsVector.hh.

   {
     G4double pt2 = sqr(p1.momentum.x()) + sqr(p1.momentum.y());
     if (pt2 > 0.0)
        {
          G4double ph, px, py, pz;
          G4double cost  = p2.momentum.z()/p2.momentum.mag(); 
          G4double sint  = 0.5 * (  std::sqrt(std::fabs((1.-cost)*(1.+cost))) 
                                  + std::sqrt(pt2)/p2.momentum.mag());
          (p2.momentum.y() < 0.) ? ph = 1.5*CLHEP::pi : ph = CLHEP::halfpi;
          if( p2.momentum.x() != 0.0) 
             ph = std::atan2(p2.momentum.y(),p2.momentum.x());             
          px =   cost*std::cos(ph)*p1.momentum.x() - std::sin(ph)*p1.momentum.y()
               + sint*std::cos(ph)*p1.momentum.z();
          py =   cost*std::sin(ph)*p1.momentum.x() + std::cos(ph)*p1.momentum.y()
               + sint*std::sin(ph)*p1.momentum.z();
          pz = - sint        *p1.momentum.x() 
               + cost        *p1.momentum.z();
          momentum.setX( px ); 
          momentum.setY( py );
          momentum.setZ( pz );     
        }
     else
        {
          momentum = p1.momentum;             
    } 
   }

Here is the call graph for this function:

G4double G4GHEKinematicsVector::Dot ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 487 of file G4GHEKinematicsVector.hh.

   {
     return (   p1.momentum.x() * p2.momentum.x()
              + p1.momentum.y() * p2.momentum.y()
              + p1.momentum.z() * p2.momentum.z() );
   }

Here is the call graph for this function:

G4double G4GHEKinematicsVector::Dot4 ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 438 of file G4GHEKinematicsVector.hh.

   {
     return (   p1.energy       * p2.energy
              - p1.momentum.x() * p2.momentum.x()
              - p1.momentum.y() * p2.momentum.y()
              - p1.momentum.z() * p2.momentum.z() );
   } 

Here is the call graph for this function:

void G4GHEKinematicsVector::Exch ( G4GHEKinematicsVector &  p1 )

inline

Definition at line 548 of file G4GHEKinematicsVector.hh.

   {
     G4GHEKinematicsVector mx = *this;
//     mx.momentum.SetX( momentum.x());
//     mx.momentum.SetY( momentum.y());
//     mx.momentum.SetZ( momentum.z());
//     mx.energy        = energy;
//     mx.kineticEnergy = kineticEnergy;
//     mx.mass          = mass;
//     mx.charge        = charge;
//     mx.timeOfFlight  = timeOfFlight;
//     mx.side          = side;
//     mx.flag          = flag;
//     mx.code          = code; 
//     momentum.setX( p1.momentum.x());
//     momentum.setY( p1.momentum.y());
//     momentum.setZ( p1.momentum.z());
//     energy        = p1.energy;
//     kineticEnergy = p1.kineticEnergy;
//     mass          = p1.mass;
//     charge        = p1.charge;
//     timeOfFlight  = p1.timeOfFlight;
//     side          = p1.side
//     flag          = p1.flag;
//     code          = p1.code;
     *this = p1;
     p1 = mx;
     return; 
   }      

G4double G4GHEKinematicsVector::GetCharge ( )

inline

Definition at line 311 of file G4GHEKinematicsVector.hh.

{return charge; }

G4int G4GHEKinematicsVector::GetCode ( )

inline

Definition at line 338 of file G4GHEKinematicsVector.hh.

{ return code; } 

G4double G4GHEKinematicsVector::GetEnergy ( )

inline

Definition at line 260 of file G4GHEKinematicsVector.hh.

{return energy;}

G4bool G4GHEKinematicsVector::getFlag ( )

inline

Definition at line 329 of file G4GHEKinematicsVector.hh.

{ return flag; }

G4double G4GHEKinematicsVector::GetKineticEnergy ( )

inline

Definition at line 263 of file G4GHEKinematicsVector.hh.

{return kineticEnergy;}

G4double G4GHEKinematicsVector::GetMass ( )

inline

Definition at line 305 of file G4GHEKinematicsVector.hh.

{ return mass; }

const G4ParticleMomentum G4GHEKinematicsVector::GetMomentum ( void  ) const

inline

Definition at line 123 of file G4GHEKinematicsVector.hh.

{ return momentum; }

G4ParticleDefinition* G4GHEKinematicsVector::GetParticleDef ( )

inline

Definition at line 341 of file G4GHEKinematicsVector.hh.

{ return particleDef; } 

G4int G4GHEKinematicsVector::GetSide ( )

inline

Definition at line 323 of file G4GHEKinematicsVector.hh.

{ return side; }

G4double G4GHEKinematicsVector::GetTOF ( )

inline

Definition at line 317 of file G4GHEKinematicsVector.hh.

{ return timeOfFlight; }

Here is the caller graph for this function:

G4double G4GHEKinematicsVector::Impu ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 447 of file G4GHEKinematicsVector.hh.

   {
     return ( - sqr( p1.energy      - p2.energy)  
              + sqr(p1.momentum.x() - p2.momentum.x())
              + sqr(p1.momentum.y() - p2.momentum.y()) 
              + sqr(p1.momentum.z() - p2.momentum.z()) );
   }

Here is the call graph for this function:

G4double G4GHEKinematicsVector::Length ( )

inline

Definition at line 542 of file G4GHEKinematicsVector.hh.

   {
     return  momentum.mag() ;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Lor ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 393 of file G4GHEKinematicsVector.hh.

   {
     G4double a;
     a = ( p1.momentum.dot(p2.momentum)/(p2.energy+p2.mass) - p1.energy ) / p2.mass;
     momentum.setX( p1.momentum.x()+a*p2.momentum.x() );
     momentum.setY( p1.momentum.y()+a*p2.momentum.y() );
     momentum.setZ( p1.momentum.z()+a*p2.momentum.z() );
     energy = std::sqrt( sqr(p1.mass) + momentum.mag2() );
     mass = p1.mass;
     kineticEnergy = std::max(0.,energy - mass);
     timeOfFlight  = p1.timeOfFlight;
     side          = p1.side;
     flag          = p1.flag;
     code          = p1.code;
     particleDef   = p1.particleDef;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Norz ( const G4GHEKinematicsVector &  p )

inline

Definition at line 522 of file G4GHEKinematicsVector.hh.

   {
     G4double a =   p.momentum.mag2();
     if (a > 0.0) a = 1./std::sqrt(a);
     momentum.setX( a * p.momentum.x() );
     momentum.setY( a * p.momentum.y() );
     momentum.setZ( a * p.momentum.z() );
     mass          = p.mass;
     energy        = std::sqrt(momentum.mag2() + mass*mass);
     kineticEnergy = energy - mass;
     charge        = p.charge;
     timeOfFlight  = p.timeOfFlight;
     side          = p.side;
     flag          = p.flag;
     code          = p.code; 
     particleDef   = p.particleDef; 
     return;
   }

Here is the call graph for this function:

G4GHEKinematicsVector& G4GHEKinematicsVector::operator= ( const G4GHEKinematicsVector &  p )

inline

Definition at line 90 of file G4GHEKinematicsVector.hh.

   {
     if (this != &p)
     {
        momentum.setX( p.momentum.x() );
        momentum.setY( p.momentum.y() );
        momentum.setZ( p.momentum.z() );
        energy        = p.energy;
        kineticEnergy = p.kineticEnergy;
        mass          = p.mass;
        charge        = p.charge;
        timeOfFlight  = p.timeOfFlight;
        side          = p.side;
        flag          = p.flag;
        code          = p.code;
        particleDef   = p.particleDef;
     }
    return *this;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Print ( G4int  L )

inline

Definition at line 674 of file G4GHEKinematicsVector.hh.

   {
     G4cout << "G4GHEKinematicsVector: " 
          << L << " " << momentum.x() << " " <<  momentum.y() << " " <<  momentum.z() << " "
          << energy << " " << kineticEnergy << " " << mass << " " << charge << " " 
          << timeOfFlight << " " << side << " " << flag << " " << code << particleDef << G4endl;
     return;                         
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetCharge ( G4double  c )

inline

Definition at line 308 of file G4GHEKinematicsVector.hh.

{ charge = c; return; }

void G4GHEKinematicsVector::SetCode ( G4int  c )

inline

Definition at line 332 of file G4GHEKinematicsVector.hh.

{ code = c; return; }

void G4GHEKinematicsVector::SetEnergy ( G4double  e )

inline

Definition at line 180 of file G4GHEKinematicsVector.hh.

{ energy = e; return; }

void G4GHEKinematicsVector::SetEnergyAndUpdate ( G4double  e )

inline

Definition at line 183 of file G4GHEKinematicsVector.hh.

   {
     if (e <= mass)
       { 
         energy        = mass;
         kineticEnergy = 0.;
         momentum.setX(  0.);
         momentum.setY(  0.);
         momentum.setZ(  0.);
       }
     else
       {
         energy = e;
         kineticEnergy   = energy - mass;
         G4double momold = momentum.mag();
         G4double momnew = std::sqrt(energy*energy - mass*mass);
         if (momold == 0.)
           {
             G4double cost = 1.0- 2.0*G4UniformRand();
             G4double sint = std::sqrt(1. - cost*cost);
             G4double phi  = CLHEP::twopi* G4UniformRand();
             momentum.setX( momnew * sint * std::cos(phi));
             momentum.setY( momnew * sint * std::sin(phi));
             momentum.setZ( momnew * cost);
           }
         else
           {
             momnew /= momold;
             momentum.setX(momentum.x()*momnew);
             momentum.setY(momentum.y()*momnew);
             momentum.setZ(momentum.z()*momnew);
           }
       }    
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::setFlag ( G4bool  f )

inline

Definition at line 326 of file G4GHEKinematicsVector.hh.

{ flag = f; return; }

void G4GHEKinematicsVector::SetKineticEnergy ( G4double  ekin )

inline

Definition at line 220 of file G4GHEKinematicsVector.hh.

{ kineticEnergy = ekin; return; }

void G4GHEKinematicsVector::SetKineticEnergyAndUpdate ( G4double  ekin )

inline

Definition at line 223 of file G4GHEKinematicsVector.hh.

   {
     if (ekin <= 0.)
       { 
         energy        = mass;
         kineticEnergy = 0.;
         momentum.setX(  0.);
         momentum.setY(  0.);
         momentum.setZ(  0.);
       }
     else
       {
         energy = ekin + mass;
         kineticEnergy   = ekin;
         G4double momold = momentum.mag();
         G4double momnew = std::sqrt(energy*energy - mass*mass);
         if (momold == 0.)
           {
             G4double cost = 1.0-2.0*G4UniformRand();
             G4double sint = std::sqrt(1. - cost*cost);
             G4double phi  = CLHEP::twopi* G4UniformRand();
             momentum.setX( momnew * sint * std::cos(phi));
             momentum.setY( momnew * sint * std::sin(phi));
             momentum.setZ( momnew * cost);
           }
         else
           {
             momnew /= momold;
             momentum.setX(momentum.x()*momnew);
             momentum.setY(momentum.y()*momnew);
             momentum.setZ(momentum.z()*momnew);
           }
       }    
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMass ( G4double  mas )

inline

Definition at line 266 of file G4GHEKinematicsVector.hh.

{ mass = mas; return; }

void G4GHEKinematicsVector::SetMassAndUpdate ( G4double  mas )

inline

Definition at line 269 of file G4GHEKinematicsVector.hh.

  {
    kineticEnergy = std::max(0., energy - mas);
    mass = mas;
    energy = kineticEnergy + mass;
    G4double momnew = std::sqrt(std::max(0., energy*energy - mass*mass));
    if ( momnew == 0.0) 
       {
         momentum.setX( 0.0 );
         momentum.setY( 0.0 );
         momentum.setZ( 0.0 );
       }
    else
       {
         G4double momold = momentum.mag();
         if (momold == 0.)
            { 
              G4double cost = 1.-2.*G4UniformRand();
              G4double sint = std::sqrt(1.-cost*cost);
              G4double phi  = CLHEP::twopi*G4UniformRand();
              momentum.setX( momnew*sint*std::cos(phi));
              momentum.setY( momnew*sint*std::sin(phi));
              momentum.setZ( momnew*cost);
            }
         else
            {
              momnew /= momold;
              momentum.setX( momentum.x()*momnew );
              momentum.setY( momentum.y()*momnew );
              momentum.setZ( momentum.z()*momnew );
            }
       }     
    return;
  }    

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMomentum ( G4ParticleMomentum  mom )

inline

Definition at line 111 of file G4GHEKinematicsVector.hh.

{ momentum = mom; return; };

void G4GHEKinematicsVector::SetMomentum ( G4double  x, G4double  y, G4double  z  )

inline

Definition at line 126 of file G4GHEKinematicsVector.hh.

   { 
     momentum.setX( x );
     momentum.setY( y );
     momentum.setZ( z );
     return;
   } 

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMomentum ( G4double  x, G4double  y  )

inline

Definition at line 146 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( x );
     momentum.setY( y );
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMomentum ( G4double  z )

inline

Definition at line 164 of file G4GHEKinematicsVector.hh.

   {
     momentum.setZ( z );
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMomentumAndUpdate ( G4ParticleMomentum  mom )

inline

Definition at line 114 of file G4GHEKinematicsVector.hh.

   {
     momentum      = mom;
     energy        = std::sqrt(mass*mass + momentum.mag2());
     kineticEnergy = std::max(0.,energy - mass);
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMomentumAndUpdate ( G4double  x, G4double  y, G4double  z  )

inline

Definition at line 135 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( x );
     momentum.setY( y );
     momentum.setZ( z );
     energy        = std::sqrt(mass*mass + momentum.mag2());
     kineticEnergy = std::max(0.,energy-mass);
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMomentumAndUpdate ( G4double  x, G4double  y  )

inline

Definition at line 154 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( x );
     momentum.setY( y );
     energy = std::sqrt(mass*mass + momentum.mag2());
     kineticEnergy = std::max(0.,energy-mass);
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetMomentumAndUpdate ( G4double  z )

inline

Definition at line 171 of file G4GHEKinematicsVector.hh.

   {
     momentum.setZ( z );
     energy = std::sqrt(mass*mass + momentum.mag2());
     kineticEnergy = std::max(0.,energy-mass);
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::SetParticleDef ( G4ParticleDefinition *  c )

inline

Definition at line 335 of file G4GHEKinematicsVector.hh.

{ particleDef = c; return; }

void G4GHEKinematicsVector::SetSide ( G4int  sid )

inline

Definition at line 320 of file G4GHEKinematicsVector.hh.

{ side = sid; return; }

void G4GHEKinematicsVector::SetTOF ( G4double  t )

inline

Definition at line 314 of file G4GHEKinematicsVector.hh.

{ timeOfFlight = t; return; }

void G4GHEKinematicsVector::SetZero ( )

inline

Definition at line 344 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX(  0.0 );
     momentum.setY(  0.0 );
     momentum.setZ(  0.0 );
     energy        = 0.0;
     kineticEnergy = 0.0;
     mass          = 0.0;
     charge        = 0.0;
     timeOfFlight  = 0.0;
     side          = 0;
     flag          = false;
     code          = 0;
     particleDef   = NULL;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Smul ( const G4GHEKinematicsVector &  p, G4double  h  )

inline

Definition at line 495 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( h * p.momentum.x());
     momentum.setY( h * p.momentum.y());
     momentum.setZ( h * p.momentum.z());
     return;
   }   

Here is the call graph for this function:

void G4GHEKinematicsVector::SmulAndUpdate ( const G4GHEKinematicsVector &  p, G4double  h  )

inline

Definition at line 504 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( h * p.momentum.x());
     momentum.setY( h * p.momentum.y());
     momentum.setZ( h * p.momentum.z());
     mass          = p.mass;
     energy        = std::sqrt(momentum.mag2() + mass*mass);
     kineticEnergy = energy - mass;
     charge        = p.charge;
     timeOfFlight  = p.timeOfFlight;
     side          = p.side;
     flag          = p.flag;
     code          = p.code;
     particleDef   = p.particleDef;
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Sub ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 377 of file G4GHEKinematicsVector.hh.

   {
     momentum = p1.momentum - p2.momentum;
     energy = p1.energy - p2.energy;
     G4double b = energy*energy - momentum.mag2();
     if( b < 0 )
       mass = -1. * std::sqrt( -b );
     else
       mass = std::sqrt( b );
     kineticEnergy = std::max(0.,energy - mass);
     charge        = p1.charge - p2.charge;
     code          = p1.code   - p2.code;
     particleDef   = p1.particleDef;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Sub3 ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  p2  )

inline

Definition at line 465 of file G4GHEKinematicsVector.hh.

   {
     momentum.setX( p1.momentum.x() - p2.momentum.x());
     momentum.setY( p1.momentum.y() - p2.momentum.y());
     momentum.setZ( p1.momentum.z() - p2.momentum.z());
     return;
   }

Here is the call graph for this function:

void G4GHEKinematicsVector::Trac ( const G4GHEKinematicsVector &  p1, const G4GHEKinematicsVector &  mx, const G4GHEKinematicsVector &  my, const G4GHEKinematicsVector &  mz  )

inline

Definition at line 654 of file G4GHEKinematicsVector.hh.

   {
     double px, py, pz;
     px =   mx.momentum.x()*p1.momentum.x()
          + mx.momentum.y()*p1.momentum.y()
          + mx.momentum.z()*p1.momentum.z();
     py =   my.momentum.x()*p1.momentum.x()
          + my.momentum.y()*p1.momentum.y()
          + my.momentum.z()*p1.momentum.z();
     pz =   mz.momentum.x()*p1.momentum.x()
          + mz.momentum.y()*p1.momentum.y()
          + mz.momentum.z()*p1.momentum.z();
     momentum.setX( px );
     momentum.setY( py );
     momentum.setZ( pz );
     return;
   }                 

Here is the call graph for this function:

Member Data Documentation

G4double G4GHEKinematicsVector::charge

Definition at line 687 of file G4GHEKinematicsVector.hh.

G4int G4GHEKinematicsVector::code

Definition at line 691 of file G4GHEKinematicsVector.hh.

G4double G4GHEKinematicsVector::energy

Definition at line 684 of file G4GHEKinematicsVector.hh.

G4bool G4GHEKinematicsVector::flag

Definition at line 690 of file G4GHEKinematicsVector.hh.

G4double G4GHEKinematicsVector::kineticEnergy

Definition at line 685 of file G4GHEKinematicsVector.hh.

G4double G4GHEKinematicsVector::mass

Definition at line 686 of file G4GHEKinematicsVector.hh.

G4ParticleMomentum G4GHEKinematicsVector::momentum

Definition at line 683 of file G4GHEKinematicsVector.hh.

G4ParticleDefinition* G4GHEKinematicsVector::particleDef

Definition at line 692 of file G4GHEKinematicsVector.hh.

G4int G4GHEKinematicsVector::side

Definition at line 689 of file G4GHEKinematicsVector.hh.

G4double G4GHEKinematicsVector::timeOfFlight

Definition at line 688 of file G4GHEKinematicsVector.hh.

---

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

• geant4.10.03.p01/source/processes/hadronic/util/include/G4GHEKinematicsVector.hh