TrkExchangePar.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id: TrkExchangePar.cxx,v 1.2 2005/07/18 02:56:08 zhangy Exp $
//
// Description:
//     
//
// Environment:
//      Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Authors: Steve Schaffner
//------------------------------------------------------------------------
#include "TrkBase/TrkExchangePar.h"
#include "MdcGeom/BesAngle.h"
#include <assert.h>
#include <iostream>
using std::endl;
using std::ostream;
 
//----------------------------------------------------------------------
TrkExchangePar::TrkExchangePar(const HepVector& inV, const HepMatrix& inErr) 
  : paramVec(inV) { 
//----------------------------------------------------------------------
    assert(paramVec.num_row() == nParam);
    assert(inErr.num_row() == nParam);
    assert(inErr.num_col() == nParam);
    paramErr.assign(inErr);
    paramVec[ex_phi0] = BesAngle(paramVec[ex_phi0]).rad();
}
 
//----------------------------------------------------------------------
TrkExchangePar::TrkExchangePar(const HepVector& inV, const HepSymMatrix& err) 
  :  paramVec(inV), paramErr(err) { 
//----------------------------------------------------------------------
    assert(paramVec.num_row() == nParam);
    assert(paramErr.num_row() == nParam);
    paramVec[ex_phi0] = BesAngle(paramVec[ex_phi0]).rad();
}
 
//----------------------------------------------------------------------
TrkExchangePar::TrkExchangePar(const HepVector& inV) 
  : paramVec(inV), paramErr(nParam, 1) {
//----------------------------------------------------------------------
    assert(paramVec.num_row() == nParam);
    paramVec[ex_phi0] = BesAngle(paramVec[ex_phi0]).rad();
}
 
//----------------------------------------------------------------------
TrkExchangePar::TrkExchangePar(double d0In, double phi0In, double omegaIn, 
                   double z0In, double tanDipIn) 
  : paramVec(nParam), paramErr(nParam, 1) { 
//----------------------------------------------------------------------
    paramVec[ex_d0]     = d0In;
    paramVec[ex_phi0]   = phi0In;
    paramVec[ex_omega]  = omegaIn;
    paramVec[ex_z0]     = z0In;
    paramVec[ex_tanDip] = tanDipIn;
    paramVec[ex_phi0] = BesAngle(paramVec[ex_phi0]).rad();
}
 
//----------------------------------------------------------------------
TrkExchangePar::~TrkExchangePar(){}
 
//----------------------------------------------------------------------
// Output functions: simple, detailed, and stream
 
void TrkExchangePar::print(ostream& o) const {
  o << "d0="  << d0() << " phi0="   << phi0() << " omega=" << omega()
    << " z0=" << z0() << " tanDip=" << tanDip();
}
 
void TrkExchangePar::printAll(ostream& o) const {
  print(o);
  o << endl << covariance() << endl;
}
 
 
 
ostream& operator<<(ostream& o, const TrkExchangePar& helix) {
  helix.print(o);
  return o;
}