/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/MdcPatRec/MdcRecoUtil/MdcRecoUtil-00-01-08/src/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];
  }

}

---