DifVector.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information: 
//  $Id: DifVector.cxx,v 1.2 2009/12/23 02:59:56 zhangy Exp $
//
// Description:
//  Class Implementation for |DifVector| 
//      What do i do ?
// Environment:
//  Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author List:
//  A. Snyder
//
// Copyright Information:
//  Copyright (C) 1996  SLAC
// 
// History:
//  Migration for BESIII MDC
//
//------------------------------------------------------------------------
 
#include "MdcRecoUtil/DifVector.h"
#include "MdcRecoUtil/DifRotation.h"
#include "CLHEP/Vector/ThreeVector.h"
using std::cout;
 
extern const DifVector xhat(1,0,0);
extern const DifVector yhat(0,1,0);
extern const DifVector zhat(0,0,1);
extern const DifVector nullVec(0,0,0);
 
DifVector::DifVector()
  :x(0.0),y(0.0),z(0.0)  
{} 
 
DifVector::DifVector
(const DifNumber& X,const DifNumber& Y,const DifNumber& Z)
  :x(X),y(Y),z(Z)
{}
 
DifVector::DifVector
(double X, double Y,double Z)
  :x(X),y(Y),z(Z)
{}
 
DifVector::DifVector(const Hep3Vector& v)
  :x(v.x()),y(v.y()),z(v.z())
{}
 
DifVector::DifVector(const DifVector& v)
  :x(v.x),y(v.y),z(v.z)
{}
 
HepSymMatrix DifVector::errorMatrix(const HepSymMatrix& e)const {
  HepSymMatrix temp(3);
  temp(1,1)=correlation(x,x,e);  
  temp(1,2)=correlation(x,y,e);
  temp(1,3)=correlation(x,z,e);
  temp(2,2)=correlation(y,y,e);
  temp(2,3)=correlation(y,z,e);
  temp(3,3)=correlation(z,z,e);
  return temp;
}
 
HepMatrix DifVector::jacobian()const{
  int npar=x.nPar();
  HepMatrix temp(3,npar);
  for(int i=1; i<=npar; i++){
    temp(1,i)=x.derivative(i);
    temp(2,i)=y.derivative(i);
    temp(3,i)=z.derivative(i);
  } // (int i=1; i<=npar; i++)
  return temp;
}
 
DifVector& DifVector::rotate(const DifRotation& r) {
  r.rotate(*this); return *this;
}
 
DifVector& DifVector::rotate(const DifNumber& a,const DifNumber& b,const DifNumber& g) {
  DifRotation(a,b,g).rotate(*this);
  return *this;
}
 
void DifVector::print()const {
  /*SKIP
  cout << "x:" << x;
  cout << "y:" << y;
  cout << "z:" << z;
  */
  //yzhang TEMP
  cout << "Yzhang::DifVector::print x,y,z SKIPED!\n";
}