Gabriele Cosmo
375 lines
12 KiB
C++
375 lines
12 KiB
C++
// -*- C++ -*-
|
|
// ---------------------------------------------------------------------------
|
|
//
|
|
// This file is a part of the CLHEP - a Class Library for High Energy Physics.
|
|
//
|
|
// This is the definitions of the inline member functions of the
|
|
// HepLorentzRotation class
|
|
//
|
|
|
|
namespace CLHEP {
|
|
|
|
// ---------- Constructors and Assignment:
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation() :
|
|
mxx(1.0), mxy(0.0), mxz(0.0), mxt(0.0),
|
|
myx(0.0), myy(1.0), myz(0.0), myt(0.0),
|
|
mzx(0.0), mzy(0.0), mzz(1.0), mzt(0.0),
|
|
mtx(0.0), mty(0.0), mtz(0.0), mtt(1.0) {}
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepLorentzRotation & r) :
|
|
mxx(r.mxx), mxy(r.mxy), mxz(r.mxz), mxt(r.mxt),
|
|
myx(r.myx), myy(r.myy), myz(r.myz), myt(r.myt),
|
|
mzx(r.mzx), mzy(r.mzy), mzz(r.mzz), mzt(r.mzt),
|
|
mtx(r.mtx), mty(r.mty), mtz(r.mtz), mtt(r.mtt) {}
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepRotation & r) {
|
|
set (r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepRotationX & r) {
|
|
set (r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepRotationY & r) {
|
|
set (r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepRotationZ & r) {
|
|
set (r.rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepBoost & b) {
|
|
set (b.rep4x4());
|
|
}
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepBoostX & b) {
|
|
set (b.rep4x4());
|
|
}
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepBoostY & b) {
|
|
set (b.rep4x4());
|
|
}
|
|
inline HepLorentzRotation::HepLorentzRotation(const HepBoostZ & b) {
|
|
set (b.rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::operator = (const HepLorentzRotation & r) {
|
|
mxx = r.mxx; mxy = r.mxy; mxz = r.mxz; mxt = r.mxt;
|
|
myx = r.myx; myy = r.myy; myz = r.myz; myt = r.myt;
|
|
mzx = r.mzx; mzy = r.mzy; mzz = r.mzz; mzt = r.mzt;
|
|
mtx = r.mtx; mty = r.mty; mtz = r.mtz; mtt = r.mtt;
|
|
return *this;
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::operator = (const HepRotation & m1) {
|
|
return set (m1.rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::operator = (const HepBoost & m1) {
|
|
return set (m1.rep4x4());
|
|
}
|
|
|
|
HepLorentzRotation & HepLorentzRotation::set (const Hep3Vector & p) {
|
|
return set (p.x(), p.y(), p.z());
|
|
}
|
|
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepRotation & r) {
|
|
return set (r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepRotationX & r) {
|
|
return set (r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepRotationY & r) {
|
|
return set (r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepRotationZ & r) {
|
|
return set (r.rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepBoost & bboost) {
|
|
return set (bboost.rep4x4());
|
|
}
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepBoostX & bboost) {
|
|
return set (bboost.rep4x4());
|
|
}
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepBoostY & bboost) {
|
|
return set (bboost.rep4x4());
|
|
}
|
|
inline HepLorentzRotation & HepLorentzRotation::set (const HepBoostZ & bboost) {
|
|
return set (bboost.rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(double bx,
|
|
double by,
|
|
double bz)
|
|
{
|
|
set(bx, by, bz);
|
|
}
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(const Hep3Vector & p)
|
|
{
|
|
set(p.x(), p.y(), p.z());
|
|
}
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(
|
|
const HepBoost & B, const HepRotation & R)
|
|
{
|
|
set(B, R);
|
|
}
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(
|
|
const HepRotation & R, const HepBoost & B)
|
|
{
|
|
set(R, B);
|
|
}
|
|
|
|
inline HepLorentzRotation & HepLorentzRotation::set( const HepRep4x4 & rep ) {
|
|
mxx=rep.xx_; mxy=rep.xy_; mxz=rep.xz_; mxt=rep.xt_;
|
|
myx=rep.yx_; myy=rep.yy_; myz=rep.yz_; myt=rep.yt_;
|
|
mzx=rep.zx_; mzy=rep.zy_; mzz=rep.zz_; mzt=rep.zt_;
|
|
mtx=rep.tx_; mty=rep.ty_; mtz=rep.tz_; mtt=rep.tt_;
|
|
return *this;
|
|
}
|
|
|
|
inline HepLorentzRotation ::HepLorentzRotation ( const HepRep4x4 & rep ) :
|
|
mxx(rep.xx_), mxy(rep.xy_), mxz(rep.xz_), mxt(rep.xt_),
|
|
myx(rep.yx_), myy(rep.yy_), myz(rep.yz_), myt(rep.yt_),
|
|
mzx(rep.zx_), mzy(rep.zy_), mzz(rep.zz_), mzt(rep.zt_),
|
|
mtx(rep.tx_), mty(rep.ty_), mtz(rep.tz_), mtt(rep.tt_) {}
|
|
|
|
// - Protected methods
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation(
|
|
double rxx, double rxy, double rxz, double rxt,
|
|
double ryx, double ryy, double ryz, double ryt,
|
|
double rzx, double rzy, double rzz, double rzt,
|
|
double rtx, double rty, double rtz, double rtt) :
|
|
mxx(rxx), mxy(rxy), mxz(rxz), mxt(rxt),
|
|
myx(ryx), myy(ryy), myz(ryz), myt(ryt),
|
|
mzx(rzx), mzy(rzy), mzz(rzz), mzt(rzt),
|
|
mtx(rtx), mty(rty), mtz(rtz), mtt(rtt) {}
|
|
|
|
inline void HepLorentzRotation::setBoost
|
|
(double bx, double by, double bz) {
|
|
set(bx, by, bz);
|
|
}
|
|
|
|
// ---------- Accessors:
|
|
|
|
inline double HepLorentzRotation::xx() const { return mxx; }
|
|
inline double HepLorentzRotation::xy() const { return mxy; }
|
|
inline double HepLorentzRotation::xz() const { return mxz; }
|
|
inline double HepLorentzRotation::xt() const { return mxt; }
|
|
inline double HepLorentzRotation::yx() const { return myx; }
|
|
inline double HepLorentzRotation::yy() const { return myy; }
|
|
inline double HepLorentzRotation::yz() const { return myz; }
|
|
inline double HepLorentzRotation::yt() const { return myt; }
|
|
inline double HepLorentzRotation::zx() const { return mzx; }
|
|
inline double HepLorentzRotation::zy() const { return mzy; }
|
|
inline double HepLorentzRotation::zz() const { return mzz; }
|
|
inline double HepLorentzRotation::zt() const { return mzt; }
|
|
inline double HepLorentzRotation::tx() const { return mtx; }
|
|
inline double HepLorentzRotation::ty() const { return mty; }
|
|
inline double HepLorentzRotation::tz() const { return mtz; }
|
|
inline double HepLorentzRotation::tt() const { return mtt; }
|
|
|
|
inline HepLorentzVector HepLorentzRotation::col1() const {
|
|
return HepLorentzVector ( mxx, myx, mzx, mtx );
|
|
}
|
|
inline HepLorentzVector HepLorentzRotation::col2() const {
|
|
return HepLorentzVector ( mxy, myy, mzy, mty );
|
|
}
|
|
inline HepLorentzVector HepLorentzRotation::col3() const {
|
|
return HepLorentzVector ( mxz, myz, mzz, mtz );
|
|
}
|
|
inline HepLorentzVector HepLorentzRotation::col4() const {
|
|
return HepLorentzVector ( mxt, myt, mzt, mtt );
|
|
}
|
|
|
|
inline HepLorentzVector HepLorentzRotation::row1() const {
|
|
return HepLorentzVector ( mxx, mxy, mxz, mxt );
|
|
}
|
|
inline HepLorentzVector HepLorentzRotation::row2() const {
|
|
return HepLorentzVector ( myx, myy, myz, myt );
|
|
}
|
|
inline HepLorentzVector HepLorentzRotation::row3() const {
|
|
return HepLorentzVector ( mzx, mzy, mzz, mzt );
|
|
}
|
|
inline HepLorentzVector HepLorentzRotation::row4() const {
|
|
return HepLorentzVector ( mtx, mty, mtz, mtt );
|
|
}
|
|
|
|
inline HepRep4x4 HepLorentzRotation::rep4x4() const {
|
|
return HepRep4x4( mxx, mxy, mxz, mxt,
|
|
myx, myy, myz, myt,
|
|
mzx, mzy, mzz, mzt,
|
|
mtx, mty, mtz, mtt );
|
|
}
|
|
|
|
|
|
// ------------ Subscripting:
|
|
|
|
inline HepLorentzRotation::HepLorentzRotation_row::HepLorentzRotation_row
|
|
(const HepLorentzRotation & r, int i) : rr(r), ii(i) {}
|
|
|
|
inline double
|
|
HepLorentzRotation::HepLorentzRotation_row::operator [] (int jj) const {
|
|
return rr(ii,jj);
|
|
}
|
|
|
|
inline const HepLorentzRotation::HepLorentzRotation_row
|
|
HepLorentzRotation::operator [] (int i) const {
|
|
return HepLorentzRotation_row(*this, i);
|
|
}
|
|
|
|
// ---------- Comparisons:
|
|
|
|
inline bool
|
|
HepLorentzRotation::operator == (const HepLorentzRotation & r) const {
|
|
return (mxx == r.xx() && mxy == r.xy() && mxz == r.xz() && mxt == r.xt() &&
|
|
myx == r.yx() && myy == r.yy() && myz == r.yz() && myt == r.yt() &&
|
|
mzx == r.zx() && mzy == r.zy() && mzz == r.zz() && mzt == r.zt() &&
|
|
mtx == r.tx() && mty == r.ty() && mtz == r.tz() && mtt == r.tt());
|
|
}
|
|
|
|
inline bool
|
|
HepLorentzRotation::operator != (const HepLorentzRotation & r) const {
|
|
return ! operator==(r);
|
|
}
|
|
|
|
inline bool
|
|
HepLorentzRotation::operator < ( const HepLorentzRotation & r ) const
|
|
{ return compare(r)< 0; }
|
|
inline bool
|
|
HepLorentzRotation::operator <= ( const HepLorentzRotation & r ) const
|
|
{ return compare(r)<=0; }
|
|
|
|
inline bool
|
|
HepLorentzRotation::operator >= ( const HepLorentzRotation & r ) const
|
|
{ return compare(r)>=0; }
|
|
inline bool
|
|
HepLorentzRotation::operator > ( const HepLorentzRotation & r ) const
|
|
{ return compare(r)> 0; }
|
|
|
|
inline bool HepLorentzRotation::isIdentity() const {
|
|
return (mxx == 1.0 && mxy == 0.0 && mxz == 0.0 && mxt == 0.0 &&
|
|
myx == 0.0 && myy == 1.0 && myz == 0.0 && myt == 0.0 &&
|
|
mzx == 0.0 && mzy == 0.0 && mzz == 1.0 && mzt == 0.0 &&
|
|
mtx == 0.0 && mty == 0.0 && mtz == 0.0 && mtt == 1.0);
|
|
}
|
|
|
|
// ---------- Properties:
|
|
|
|
// ---------- Application:
|
|
|
|
inline HepLorentzVector
|
|
HepLorentzRotation::vectorMultiplication(const HepLorentzVector & p) const {
|
|
double x(p.x());
|
|
double y(p.y());
|
|
double z(p.z());
|
|
double t(p.t());
|
|
return HepLorentzVector(mxx*x + mxy*y + mxz*z + mxt*t,
|
|
myx*x + myy*y + myz*z + myt*t,
|
|
mzx*x + mzy*y + mzz*z + mzt*t,
|
|
mtx*x + mty*y + mtz*z + mtt*t);
|
|
}
|
|
|
|
inline HepLorentzVector
|
|
HepLorentzRotation::operator() (const HepLorentzVector & w) const {
|
|
return vectorMultiplication(w);
|
|
}
|
|
|
|
inline HepLorentzVector
|
|
HepLorentzRotation::operator * (const HepLorentzVector & p) const {
|
|
return vectorMultiplication(p);
|
|
}
|
|
|
|
// ---------- Operations in the group of 4-Rotations
|
|
|
|
inline HepLorentzRotation
|
|
HepLorentzRotation::operator * (const HepBoost & b) const {
|
|
return matrixMultiplication(b.rep4x4());
|
|
}
|
|
inline HepLorentzRotation
|
|
HepLorentzRotation::operator * (const HepRotation & r) const {
|
|
return matrixMultiplication(r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation
|
|
HepLorentzRotation::operator * (const HepLorentzRotation & lt) const {
|
|
return matrixMultiplication(lt.rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::operator *= (const HepBoost & b) {
|
|
return *this = matrixMultiplication(b.rep4x4());
|
|
}
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::operator *= (const HepRotation & r) {
|
|
return *this = matrixMultiplication(r.rep4x4());
|
|
}
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::operator *= (const HepLorentzRotation & lt) {
|
|
return *this = matrixMultiplication(lt.rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::transform (const HepBoost & b) {
|
|
return *this = HepLorentzRotation(b).matrixMultiplication(rep4x4());
|
|
}
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::transform (const HepRotation & r) {
|
|
return *this = HepLorentzRotation(r).matrixMultiplication(rep4x4());
|
|
}
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::transform (const HepLorentzRotation & lt) {
|
|
return *this = lt.matrixMultiplication(rep4x4());
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::rotate(double angle, const Hep3Vector & axis) {
|
|
return transform(HepRotation().rotate(angle, axis));
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::rotate(double angle, const Hep3Vector * axis) {
|
|
return transform(HepRotation().rotate(angle, axis));
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::boost(double bx, double by, double bz) {
|
|
return transform(HepLorentzRotation(bx, by, bz));
|
|
}
|
|
|
|
inline HepLorentzRotation &
|
|
HepLorentzRotation::boost(const Hep3Vector & b) {
|
|
return transform(HepLorentzRotation(b));
|
|
}
|
|
|
|
inline HepLorentzRotation HepLorentzRotation::inverse() const {
|
|
return HepLorentzRotation( mxx, myx, mzx, -mtx,
|
|
mxy, myy, mzy, -mty,
|
|
mxz, myz, mzz, -mtz,
|
|
-mxt, -myt, -mzt, mtt );
|
|
}
|
|
|
|
inline HepLorentzRotation & HepLorentzRotation::invert() {
|
|
return *this = inverse();
|
|
}
|
|
|
|
inline HepLorentzRotation inverseOf ( const HepLorentzRotation & lt ) {
|
|
return HepLorentzRotation(
|
|
HepRep4x4(
|
|
lt.mxx, lt.myx, lt.mzx, -lt.mtx,
|
|
lt.mxy, lt.myy, lt.mzy, -lt.mty,
|
|
lt.mxz, lt.myz, lt.mzz, -lt.mtz,
|
|
-lt.mxt, -lt.myt, -lt.mzt, lt.mtt ) );
|
|
}
|
|
|
|
inline double HepLorentzRotation::getTolerance() {
|
|
return Hep4RotationInterface::tolerance;
|
|
}
|
|
inline double HepLorentzRotation::setTolerance(double tol) {
|
|
return Hep4RotationInterface::setTolerance(tol);
|
|
}
|
|
|
|
} // namespace CLHEP
|