DifArray.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//  $Id: DifArray.cxx,v 1.3 2009/12/23 02:59:56 zhangy Exp $
//
// Description:
//  Class Implementation for |DifArray|
//      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/DifArray.h"
#include "MdcRecoUtil/DifNumber.h"
#include <iostream>
using std::endl;
// CHAGNE using std::ostream;
using std::cout;
 
DifArray::DifArray( int n, int npar ) : _nElem( n ), _pointer( new DifNumber[n] ) {
  zero( npar );
}
 
DifArray::DifArray( const HepVector& a, int npar )
    : _nElem( a.num_row() ), _pointer( new DifNumber[nElem()] ) {
  copy( a, npar );
}
 
DifArray::DifArray( const DifArray& a )
    : _nElem( a.nElem() ), _pointer( new DifNumber[nElem()] ) {
  copy( a );
}
 
DifArray::~DifArray() { delete[] _pointer; }
 
DifNumber& DifArray::operator[]( int i ) {
  assert( i >= 0 );
  assert( i < nElem() );
  return _pointer[i];
}
 
DifNumber& DifArray::operator()( int i ) {
  i = i - 1;
  assert( i >= 0 );
  assert( i < nElem() );
  return _pointer[i];
}
 
DifNumber DifArray::fetch( int i ) const {
  i = i - 1;
  assert( i >= 0 );
  assert( i < nElem() );
  return _pointer[i];
}
 
DifArray& DifArray::operator=( const DifArray& rhs ) {
  copy( rhs );
  return *this;
}
 
HepMatrix DifArray::jacobian() const {
  int       npar = _pointer[0].nPar();
  HepMatrix temp( nElem(), npar, 0 );
  for ( int i = 1; i <= nElem(); i++ )
  {
    for ( int j = 1; j <= npar; j++ ) { temp( i, j ) = _pointer[i - 1].derivative( j ); }
  } //(int i=1; i<=nElem(); i++)
  return temp;
}
 
void DifArray::copy( const HepVector& a, int npar ) {
  assert( nElem() == a.num_row() );
  for ( int i = 0; i < nElem(); i++ )
  {
    _pointer[i].setNumber( a( i + 1 ) );
    _pointer[i].setNPar( npar );
  }
}
 
void DifArray::copy( const DifArray& a ) {
  assert( nElem() == a.nElem() );
  for ( int i = 0; i < nElem(); i++ ) { _pointer[i] = a._pointer[i]; }
}
 
void DifArray::zero( int npar ) {
  for ( int i = 0; i < nElem(); i++ )
  {
    _pointer[i] = 0.0;
    _pointer[i].setNPar( npar );
  }
}
 
void DifArray::print( /*ostream& o*/ ) const {
  cout << "nElem=" << nElem() << endl;
  for ( int i = 1; i <= nElem(); i++ )
  {
    // SKIP cout << "element(" << i << ")=" << _pointer[i-1];
  }
}