STEER/STEER/AliRecPoint.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 /* $Id$ */
  17
  18 //-*-C++-*-
  19 //_________________________________________________________________________
  20 // Base Class for reconstructed space points
  21 // usually coming from the clusterisation algorithms
  22 // run on the digits
  23 //
  24 //*-- Author : Yves Schutz  SUBATECH
  25 //////////////////////////////////////////////////////////////////////////////
  26
  27
  28 // --- ROOT system ---
  29
  30 // --- Standard library ---
  31
  32 // --- AliRoot header files ---
  33
  34 #include "AliRecPoint.h"
  35 #include "AliGeometry.h"
  36 #include "AliDigitNew.h"
  37
  38 ClassImp(AliRecPoint)
  39
  40
  41 //_______________________________________________________________________
  42 AliRecPoint::AliRecPoint():
  43   fAmp(0),
  44   fGeom(0),
  45   fIndexInList(-1), // to be set when the point is already stored
  46   fLocPos(0,0,0),
  47   fLocPosM(0),
  48   fMaxDigit(100),
  49   fMulDigit(0),
  50   fMaxTrack(5),
  51   fMulTrack(0),
  52   fDigitsList(0),
  53   fTracksList(0)
  54 {
  55   //
  56   // default ctor
  57   //
  58 }
  59
  60 //_______________________________________________________________________
  61 AliRecPoint::AliRecPoint(const char * ):
  62   fAmp(0),
  63   fGeom(0),
  64   fIndexInList(-1), // to be set when the point is already stored
  65   fLocPos(0,0,0),
  66   fLocPosM(new TMatrixF(3,3)),
  67   fMaxDigit(100),
  68   fMulDigit(0),
  69   fMaxTrack(5),
  70   fMulTrack(0),
  71   fDigitsList(new int[fMaxDigit]),
  72   fTracksList(new int[fMaxTrack])
  73 {
  74   //
  75   // Standard ctor
  76   //
  77 }
  78
  79 //_______________________________________________________________________
  80 AliRecPoint::AliRecPoint(const AliRecPoint& recp):
  81   TObject(recp),
  82   fAmp(recp.fAmp),
  83   fGeom(recp.fGeom),
  84   fIndexInList(-1), // to be set when the point is already stored
  85   fLocPos(recp.fLocPos),
  86   fLocPosM(new TMatrixF(*(recp.fLocPosM))),
  87   fMaxDigit(recp.fMaxDigit),
  88   fMulDigit(recp.fMulDigit),
  89   fMaxTrack(recp.fMaxTrack),
  90   fMulTrack(recp.fMulTrack),
  91   fDigitsList(new Int_t[fMulDigit]),
  92   fTracksList(new Int_t[fMulTrack])
  93 {
  94   //
  95   // Copy constructor
  96   //
  97   memcpy(fDigitsList,recp.fDigitsList,sizeof(Int_t)*fMulDigit);
  98   memcpy(fTracksList,recp.fTracksList,sizeof(Int_t)*fMulTrack);
  99 }
 100
 101 //_______________________________________________________________________
 102 AliRecPoint& AliRecPoint::operator=(const AliRecPoint& recp)
 103 {
 104   //
 105   // Assignement constructor
 106   //
 107   if(this!=&recp) {
 108     TObject::operator=(recp);
 109     fAmp=recp.fAmp;
 110     fGeom=recp.fGeom;
 111     fIndexInList=-1; // to be set when the point is already stored
 112     fLocPos=recp.fLocPos;
 113     delete fLocPosM;
 114     fLocPosM=recp.fLocPosM;
 115     fMaxDigit=recp.fMaxDigit;
 116     fMulDigit=recp.fMulDigit;
 117     fMaxTrack=recp.fMaxTrack;
 118     fMulTrack=recp.fMulTrack;
 119     delete [] fDigitsList;
 120     fDigitsList=new Int_t[fMulDigit];
 121     delete [] fTracksList;
 122     fTracksList=new Int_t[fMulTrack];
 123
 124     memcpy(fDigitsList,recp.fDigitsList,sizeof(Int_t)*fMulDigit);
 125     memcpy(fTracksList,recp.fTracksList,sizeof(Int_t)*fMulTrack);
 126   }
 127   return *this;
 128 }
 129
 130 //_______________________________________________________________________
 131 AliRecPoint::~AliRecPoint()
 132 {
 133   // dtor
 134
 135   delete fLocPosM ;
 136   delete [] fDigitsList ;
 137   delete [] fTracksList ;
 138
 139 }
 140
 141 //_______________________________________________________________________
 142 void AliRecPoint::AddDigit(AliDigitNew & digit)
 143 {
 144   // adds a digit to the digits list
 145   // and accumulates the total amplitude and the multiplicity
 146
 147
 148   if ( fMulDigit >= fMaxDigit ) { // increase the size of the list
 149     int * tempo = new int[fMaxDigit*2];
 150
 151     Int_t index ;
 152
 153     for ( index = 0 ; index < fMulDigit ; index++ )
 154       tempo[index] = fDigitsList[index] ;
 155
 156     delete [] fDigitsList ;
 157     fDigitsList = tempo ;
 158   }
 159
 160   fDigitsList[fMulDigit] = digit.GetIndexInList()  ;
 161   fMulDigit++ ;
 162   fAmp += digit.GetAmp() ;
 163 }
 164
 165 //_______________________________________________________________________
 166 // void AliRecPoint::AddTrack(AliTrack & track)
 167 // {
 168 //   // adds a digit to the digits list
 169 //   // and accumulates the total amplitude and the multiplicity
 170
 171
 172 //   if ( fMulTrack >= fMaxTrack ) { // increase the size of the list
 173 //     int * tempo = new int[fMaxTrack*=2] ;
 174 //     Int_t index ;
 175 //     for ( index = 0 ; index < fMulTrack ; index++ )
 176 //       tempo[index] = fTracksList[index] ;
 177 //     delete fTracksList ;
 178 //     fTracksList = tempo ;
 179 //   }
 180
 181 //   fTracksList[fMulTrack++]=  (int) &Track  ;
 182 // }
 183
 184 //_______________________________________________________________________
 185 void AliRecPoint::Copy(TObject& recp) const
 186 {
 187   //
 188   // Copy *this onto pts
 189   //
 190   // Copy all first
 191   if((TObject*)this != &recp) {
 192     ((TObject*) this)->Copy(recp);
 193     (dynamic_cast<AliRecPoint&>(recp)).fAmp = fAmp;
 194     (dynamic_cast<AliRecPoint&>(recp)).fGeom = fGeom;
 195     (dynamic_cast<AliRecPoint&>(recp)).fIndexInList = fIndexInList;
 196     (dynamic_cast<AliRecPoint&>(recp)).fLocPos = fLocPos;
 197     (dynamic_cast<AliRecPoint&>(recp)).fLocPosM = new TMatrixF(*fLocPosM);
 198     (dynamic_cast<AliRecPoint&>(recp)).fMaxDigit = fMaxDigit;
 199     (dynamic_cast<AliRecPoint&>(recp)).fMulDigit = fMulDigit;
 200     (dynamic_cast<AliRecPoint&>(recp)).fMaxTrack = fMaxTrack;
 201     (dynamic_cast<AliRecPoint&>(recp)).fMulTrack = fMulTrack;
 202
 203     // Duplicate pointed objects
 204     (dynamic_cast<AliRecPoint&>(recp)).fDigitsList = new Int_t[fMulDigit];
 205     memcpy((dynamic_cast<AliRecPoint&>(recp)).fDigitsList,fDigitsList,fMulDigit*sizeof(Int_t));
 206     (dynamic_cast<AliRecPoint&>(recp)).fTracksList = new Int_t[fMulTrack];
 207     memcpy((dynamic_cast<AliRecPoint&>(recp)).fTracksList,fTracksList,fMulTrack*sizeof(Int_t));
 208   }
 209 }
 210
 211 //_______________________________________________________________________
 212 void AliRecPoint::GetCovarianceMatrix(TMatrixF & mat) const
 213 {
 214   // returns the covariant matrix for the local position
 215
 216   mat = *fLocPosM ;
 217
 218 }
 219
 220 //____________________________________________________________________________
 221 void AliRecPoint::GetLocalPosition(TVector3 & pos) const
 222 {
 223   // returns the position of the cluster in the local reference system of the sub-detector
 224
 225   pos = fLocPos;
 226
 227
 228 }
 229
 230 //____________________________________________________________________________
 231 void AliRecPoint::GetGlobalPosition(TVector3 & gpos, TMatrixF & gmat) const
 232 {
 233   // returns the position of the cluster in the global reference system of ALICE
 234   // and the uncertainty on this position
 235
 236
 237   fGeom->GetGlobal(this, gpos, gmat) ;
 238
 239 }
 240
 241