GeoXPEngine Class Reference

#include <GeoXPEngine.h>

Collaboration diagram for GeoXPEngine:

Public Member Functions
	GeoXPEngine (const Genfun::AbsFunction &Bx, const Genfun::AbsFunction &By, const Genfun::AbsFunction &Bz, const HepGeom::Point3D< double > &x0, const CLHEP::HepLorentzVector &p0, double q)
	~GeoXPEngine ()
const Genfun::AbsFunction &	x () const
const Genfun::AbsFunction &	y () const
const Genfun::AbsFunction &	z () const
const Genfun::AbsFunction &	px () const
const Genfun::AbsFunction &	py () const
const Genfun::AbsFunction &	pz () const

Private Attributes
std::unique_ptr< const Genfun::AbsFunction >	m_x
std::unique_ptr< const Genfun::AbsFunction >	m_y
std::unique_ptr< const Genfun::AbsFunction >	m_z
std::unique_ptr< const Genfun::AbsFunction >	m_px
std::unique_ptr< const Genfun::AbsFunction >	m_py
std::unique_ptr< const Genfun::AbsFunction >	m_pz

Detailed Description

Definition at line 23 of file GeoXPEngine.h.

Constructor & Destructor Documentation

GeoXPEngine()

GeoXPEngine::GeoXPEngine	(	const Genfun::AbsFunction &	Bx,
		const Genfun::AbsFunction &	By,
		const Genfun::AbsFunction &	Bz,
		const HepGeom::Point3D< double > &	x0,
		const CLHEP::HepLorentzVector &	p0,
		double	q )

Definition at line 14 of file GeoXPEngine.cxx.

{
 
 
  // Peform a Runge-Kutte Integration of the trajectory of the charged
  // particle in a magnetic field.  Physics input is the Lorentz force
  // law; note, relativistic expression applies.
 
  double m = p0.m();
  double E = p0.e();
  Genfun::Variable t;
 
 
  if (m<100*CLHEP::keV) { 
    m_x.reset ( (Genfun::FixedConstant(x0.x()) + (CLHEP::c_light*p0.x()/p0.vect().mag())*t).clone() );
    m_y.reset ( (Genfun::FixedConstant(x0.y()) + (CLHEP::c_light*p0.y()/p0.vect().mag())*t).clone() );
    m_z.reset ( (Genfun::FixedConstant(x0.z()) + (CLHEP::c_light*p0.z()/p0.vect().mag())*t).clone() );
    m_px.reset ( Genfun::FixedConstant(p0.x()).clone() );
    m_py.reset ( Genfun::FixedConstant(p0.y()).clone() );
    m_pz.reset ( Genfun::FixedConstant(p0.z()).clone() );
    
 
  }
  else {
 
    Genfun::RKIntegrator rkIntegrator;
    Genfun::Variable Px(0,6),Py(1,6), Pz(2,6), X(3,6), Y(4,6), Z(5,6);
    Genfun::FixedConstant I(1.0);
    
    const double inv_E = 1. / E;
    Genfun::GENFUNCTION DPxDt  = (q*inv_E)*(Py*(I%I%I%Bz) -Pz*(I%I%I%By))*CLHEP::c_light;
    Genfun::GENFUNCTION DPyDt  = (q*inv_E)*(Pz*(I%I%I%Bx) -Px*(I%I%I%Bz))*CLHEP::c_light;
    Genfun::GENFUNCTION DPzDt  = (q*inv_E)*(Px*(I%I%I%By) -Py*(I%I%I%Bx))*CLHEP::c_light;
    Genfun::GENFUNCTION DxDt   = Px*inv_E;
    Genfun::GENFUNCTION DyDt   = Py*inv_E;
    Genfun::GENFUNCTION DzDt   = Pz*inv_E;
    
    rkIntegrator.addDiffEquation(&DPxDt, "Px",p0.x(), p0.x(), p0.x());
    rkIntegrator.addDiffEquation(&DPyDt, "Py",p0.y(), p0.y(), p0.y());
    rkIntegrator.addDiffEquation(&DPzDt, "Pz",p0.z(), p0.z(), p0.z());
    rkIntegrator.addDiffEquation(&DxDt,  "x", x0.x(), x0.x(), x0.x());
    rkIntegrator.addDiffEquation(&DyDt,  "y", x0.y(), x0.y(), x0.y());
    rkIntegrator.addDiffEquation(&DzDt,  "z", x0.z(), x0.z(), x0.z());
    
   
 
    m_px.reset ( (*rkIntegrator.getFunction(0))(CLHEP::c_light*t).clone() );
    m_py.reset ( (*rkIntegrator.getFunction(1))(CLHEP::c_light*t).clone() );
    m_pz.reset ( (*rkIntegrator.getFunction(2))(CLHEP::c_light*t).clone() );
    m_x.reset ( (*rkIntegrator.getFunction(3))(CLHEP::c_light*t).clone() );
    m_y.reset ( (*rkIntegrator.getFunction(4))(CLHEP::c_light*t).clone() );
    m_z.reset ( (*rkIntegrator.getFunction(5))(CLHEP::c_light*t).clone() );
  }
}

~GeoXPEngine()

GeoXPEngine::~GeoXPEngine

(

)

Definition at line 74 of file GeoXPEngine.cxx.

                          {
}

Member Function Documentation

px()

const Genfun::AbsFunction & GeoXPEngine::px

(

)

const

Definition at line 86 of file GeoXPEngine.cxx.

                                                 {
  return *m_px;
}

py()

const Genfun::AbsFunction & GeoXPEngine::py

(

)

const

Definition at line 89 of file GeoXPEngine.cxx.

                                                 {
  return *m_py;
}

pz()

const Genfun::AbsFunction & GeoXPEngine::pz

(

)

const

Definition at line 92 of file GeoXPEngine.cxx.

                                                 {
  return *m_pz;
}

x()

const Genfun::AbsFunction & GeoXPEngine::x

(

)

const

Definition at line 77 of file GeoXPEngine.cxx.

                                                {
  return *m_x;
}

y()

const Genfun::AbsFunction & GeoXPEngine::y

(

)

const

Definition at line 80 of file GeoXPEngine.cxx.

                                                {
  return *m_y;
}

z()

const Genfun::AbsFunction & GeoXPEngine::z

(

)

const

Definition at line 83 of file GeoXPEngine.cxx.

                                                {
  return *m_z;
}

Member Data Documentation

m_px

std::unique_ptr<const Genfun::AbsFunction> GeoXPEngine::m_px

private

Definition at line 59 of file GeoXPEngine.h.

m_py

std::unique_ptr<const Genfun::AbsFunction> GeoXPEngine::m_py

private

Definition at line 60 of file GeoXPEngine.h.

m_pz

std::unique_ptr<const Genfun::AbsFunction> GeoXPEngine::m_pz

private

Definition at line 61 of file GeoXPEngine.h.

m_x

std::unique_ptr<const Genfun::AbsFunction> GeoXPEngine::m_x

private

Definition at line 56 of file GeoXPEngine.h.

m_y

std::unique_ptr<const Genfun::AbsFunction> GeoXPEngine::m_y

private

Definition at line 57 of file GeoXPEngine.h.

m_z

std::unique_ptr<const Genfun::AbsFunction> GeoXPEngine::m_z

private

Definition at line 58 of file GeoXPEngine.h.

---

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

• GeoXPEngine.h
• GeoXPEngine.cxx