NeutParams.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Description:
//  Code for the NeutParams neutral track parameterization class
//
// Environment:
//      Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author(s): Justin Albert, Valery Miftahov (based on HelixParams.cc by
//                                            Dave Brown)
//
//------------------------------------------------------------------------
#include "TrkBase/NeutParams.h"
#include "TrkBase/TrkExchangePar.h"
#include <iostream>
using namespace std;
// construct from vector and covariance matrix
//------------------------------------------------------------------------
NeutParams::NeutParams( const HepVector& pvec, const HepSymMatrix& pcov )
    : TrkParams( pvec, pcov ) {
  //------------------------------------------------------------------------
 
  //  Make sure the dimensions of the input matrix and vector are correct
 
  if ( pvec.num_row() != _nneutprm || pcov.num_row() != _nneutprm )
  {
    cout << "ErrMsg(error)"
         << "NeutParams: incorrect constructor vector/matrix dimension" << endl;
    parameter()  = HepVector( _nneutprm, 0 );
    covariance() = HepSymMatrix( _nneutprm, 0 );
  }
}
 
//  Construct from the fit parameters directly
//------------------------------------------------------------------------
NeutParams::NeutParams( double n_d0, double n_phi0, double n_p, double n_z0, double n_tanDip,
                        double n_s0 )
    : //------------------------------------------------------------------------
    TrkParams( _nneutprm ) {
  d0()     = n_d0;
  phi0()   = n_phi0;
  p()      = n_p;
  z0()     = n_z0;
  s0()     = n_s0;
  tanDip() = n_tanDip;
}
 
//  Copy constructor
//------------------------------------------------------------------------
NeutParams::NeutParams( const NeutParams& old )
    : //------------------------------------------------------------------------
    TrkParams( old ) {}
//------------------------------------------------------------------------
double NeutParams::sinPhi0() const {
  //------------------------------------------------------------------------
  return sin( parameter()[_phi0] );
}
 
//------------------------------------------------------------------------
double NeutParams::cosPhi0() const {
  //------------------------------------------------------------------------
  return cos( parameter()[_phi0] );
}
 
//------------------------------------------------------------------------
double NeutParams::arcRatio() const {
  //------------------------------------------------------------------------
  return sqrt( 1. + parameter()[_tanDip] * parameter()[_tanDip] );
}
 
//------------------------------------------------------------------------
NeutParams::~NeutParams() {}
//------------------------------------------------------------------------