STEER/AliTrackFitter.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 //-----------------------------------------------------------------
  17 //   Implementation of the base class for fast track fitters
  18 //
  19 //
  20 //-----------------------------------------------------------------
  21
  22 #include <TMatrixDSym.h>
  23 #include <TArrayI.h>
  24
  25 #include "AliTrackFitter.h"
  26 #include "AliTrackPointArray.h"
  27 #include "AliLog.h"
  28
  29 ClassImp(AliTrackFitter)
  30
  31 //_____________________________________________________________________________
  32 AliTrackFitter::AliTrackFitter() :
  33   TObject(),
  34   fCov(0),
  35   fPoints(0),
  36   fPVolId(0),
  37   fPTrack(0),
  38   fChi2(0),
  39   fNdf(0),
  40   fMinNPoints(0),
  41   fIsOwner(kFALSE)
  42 {
  43   // default constructor
  44   //
  45   for (Int_t i=0;i<6;i++) fParams[i] = 0;
  46 }
  47
  48 //_____________________________________________________________________________
  49 AliTrackFitter::AliTrackFitter(AliTrackPointArray *array, Bool_t owner) :
  50   TObject(),
  51   fCov(new TMatrixDSym(6)),
  52   fPoints(0),
  53   fPVolId(0),
  54   fPTrack(0),
  55   fChi2(0),
  56   fNdf(0),
  57   fMinNPoints(0),
  58   fIsOwner(kFALSE)
  59
  60 {
  61   // constructor from space points array
  62   //
  63   for (Int_t i=0;i<6;i++) fParams[i] = 0;
  64   SetTrackPointArray(array,owner);
  65 }
  66
  67 //_____________________________________________________________________________
  68 AliTrackFitter::AliTrackFitter(const AliTrackFitter &fitter):
  69   TObject(fitter),
  70   fCov(new TMatrixDSym(*fitter.fCov)),
  71   fPoints(0),
  72   fPVolId(0),
  73   fPTrack(0),
  74   fChi2(fitter.fChi2),
  75   fNdf(fitter.fNdf),
  76   fMinNPoints(fitter.fMinNPoints),
  77   fIsOwner(kFALSE)
  78 {
  79   // Copy constructor
  80   //
  81   SetTrackPointArray(fitter.fPoints,fitter.fIsOwner);
  82   for (Int_t i=0;i<6;i++) fParams[i] = fitter.fParams[i];
  83 }
  84
  85 //_____________________________________________________________________________
  86 AliTrackFitter &AliTrackFitter::operator =(const AliTrackFitter& fitter)
  87 {
  88   // assignment operator
  89   //
  90   if(this==&fitter) return *this;
  91
  92   SetTrackPointArray(fitter.fPoints);
  93   for (Int_t i=0;i<6;i++) fParams[i] = fitter.fParams[i];
  94   fCov = new TMatrixDSym(*fitter.fCov);
  95   fChi2 = fitter.fChi2;
  96   fNdf = fitter.fNdf;
  97   fMinNPoints = fitter.fMinNPoints;
  98   fIsOwner = kFALSE;
  99
 100   return *this;
 101 }
 102
 103 //_____________________________________________________________________________
 104 AliTrackFitter::~AliTrackFitter()
 105 {
 106   if (fIsOwner)
 107     delete fPoints;
 108   delete fCov;
 109 }
 110
 111 //_____________________________________________________________________________
 112 void AliTrackFitter::Reset()
 113 {
 114   for (Int_t i=0;i<6;i++) fParams[i] = 0;
 115   delete fCov;
 116   fCov = new TMatrixDSym(6);
 117   fPVolId = fPTrack = 0;
 118   fChi2 = 0;
 119   fNdf = 0;
 120 }
 121
 122 void AliTrackFitter::SetTrackPointArray(AliTrackPointArray *array, Bool_t owner)
 123 {
 124   // Load space points from array
 125   // By default we don't copy them but
 126   // just put the pointers to them
 127   if (!array) {
 128     AliWarning("Invalid pointer to the space-points array !");
 129     if (fIsOwner) delete fPoints;
 130     fPoints = NULL;
 131     return;
 132   }
 133
 134   Reset();
 135
 136   if (fIsOwner) delete fPoints;
 137
 138   if (owner) {
 139     fPoints = new AliTrackPointArray(*array);
 140     fIsOwner = kTRUE;
 141   }
 142   else {
 143     fPoints = array;
 144     fIsOwner = kFALSE;
 145   }
 146 }
 147
 148 Bool_t AliTrackFitter::FindVolId(const TArrayI *array, UShort_t volid) const
 149 {
 150   // The method is used to check whenever
 151   // the volume id (volid) is contained in
 152   // a array of integers
 153   Int_t nVolIds = array->GetSize();
 154   if (nVolIds == 0) return kFALSE;
 155
 156   Bool_t found = kFALSE;
 157   for (Int_t iVolId = 0; iVolId < nVolIds; iVolId++) {
 158     if ((*array)[iVolId] == volid) {
 159       found = kTRUE;
 160       break;
 161     }
 162   }
 163
 164   return found;
 165 }