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 {
  34   // default constructor
  35   //
  36   for (Int_t i=0;i<6;i++) fParams[i] = 0;
  37   fCov = 0;
  38   fPoints = 0;
  39   fPVolId = fPTrack = 0;
  40   fChi2 = 0;
  41   fNdf = 0;
  42   fMinNPoints = 0;
  43   fIsOwner = kFALSE;
  44 }
  45
  46 //_____________________________________________________________________________
  47 AliTrackFitter::AliTrackFitter(AliTrackPointArray *array, Bool_t owner)
  48 {
  49   // constructor from space points array
  50   //
  51   for (Int_t i=0;i<6;i++) fParams[i] = 0;
  52   fCov = new TMatrixDSym(6);
  53   fPVolId = fPTrack = 0;
  54   fChi2 = 0;
  55   fNdf = 0;
  56   fMinNPoints = 0;
  57   fIsOwner = kFALSE;
  58   SetTrackPointArray(array,owner);
  59 }
  60
  61 //_____________________________________________________________________________
  62 AliTrackFitter::AliTrackFitter(const AliTrackFitter &fitter):
  63   TObject(fitter)
  64 {
  65   // Copy constructor
  66   //
  67   SetTrackPointArray(fitter.fPoints,fitter.fIsOwner);
  68   for (Int_t i=0;i<6;i++) fParams[i] = fitter.fParams[i];
  69   fCov = new TMatrixDSym(*fitter.fCov);
  70   fChi2 = fitter.fChi2;
  71   fNdf = fitter.fNdf;
  72   fMinNPoints = fitter.fMinNPoints;
  73   fIsOwner = kFALSE;
  74 }
  75
  76 //_____________________________________________________________________________
  77 AliTrackFitter &AliTrackFitter::operator =(const AliTrackFitter& fitter)
  78 {
  79   // assignment operator
  80   //
  81   if(this==&fitter) return *this;
  82
  83   SetTrackPointArray(fitter.fPoints);
  84   for (Int_t i=0;i<6;i++) fParams[i] = fitter.fParams[i];
  85   fCov = new TMatrixDSym(*fitter.fCov);
  86   fChi2 = fitter.fChi2;
  87   fNdf = fitter.fNdf;
  88   fMinNPoints = fitter.fMinNPoints;
  89   fIsOwner = kFALSE;
  90
  91   return *this;
  92 }
  93
  94 //_____________________________________________________________________________
  95 AliTrackFitter::~AliTrackFitter()
  96 {
  97   if (fIsOwner)
  98     delete fPoints;
  99   delete fCov;
 100 }
 101
 102 //_____________________________________________________________________________
 103 void AliTrackFitter::Reset()
 104 {
 105   for (Int_t i=0;i<6;i++) fParams[i] = 0;
 106   delete fCov;
 107   fCov = new TMatrixDSym(6);
 108   fPVolId = fPTrack = 0;
 109   fChi2 = 0;
 110   fNdf = 0;
 111 }
 112
 113 void AliTrackFitter::SetTrackPointArray(AliTrackPointArray *array, Bool_t owner)
 114 {
 115   // Load space points from array
 116   // By default we don't copy them but
 117   // just put the pointers to them
 118   if (!array) {
 119     AliWarning("Invalid pointer to the space-points array !");
 120     if (fIsOwner) delete fPoints;
 121     fPoints = NULL;
 122     return;
 123   }
 124
 125   Reset();
 126
 127   if (fIsOwner) delete fPoints;
 128
 129   if (owner) {
 130     fPoints = new AliTrackPointArray(*array);
 131     fIsOwner = kTRUE;
 132   }
 133   else {
 134     fPoints = array;
 135     fIsOwner = kFALSE;
 136   }
 137 }
 138
 139 Bool_t AliTrackFitter::FindVolId(const TArrayI *array, UShort_t volid) const
 140 {
 141   // The method is used to check whenever
 142   // the volume id (volid) is contained in
 143   // a array of integers
 144   Int_t nVolIds = array->GetSize();
 145   if (nVolIds == 0) return kFALSE;
 146
 147   Bool_t found = kFALSE;
 148   for (Int_t iVolId = 0; iVolId < nVolIds; iVolId++) {
 149     if ((*array)[iVolId] == volid) {
 150       found = kTRUE;
 151       break;
 152     }
 153   }
 154
 155   return found;
 156 }