]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PMD/AliPMDRecPoint.cxx
Better sorting algorithm
[u/mrichter/AliRoot.git] / PMD / AliPMDRecPoint.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 //_________________________________________________________________________
19 // Class for PMD reconstructed space points 
20 // usually coming from the clusterisation algorithms
21 // run on the digits
22 //
23 //////////////////////////////////////////////////////////////////////////////
24
25 // --- ROOT system ---
26
27 #include "TObjArray.h"
28
29 // --- Standard library ---
30
31 // --- AliRoot header files ---
32
33 #include "AliPMDRecPoint.h"
34 #include "AliGeometry.h"
35 #include "AliDigitNew.h"
36
37 ClassImp(AliPMDRecPoint)
38
39
40 //____________________________________________________________________________
41 //____________________________________________________________________________
42 void AliPMDRecPoint::AddDigit(AliDigitNew & digit)
43 {
44   // adds a digit to the digits list
45   // and accumulates the total amplitude and the multiplicity 
46   
47   
48   if ( fMulDigit >= fMaxDigit ) { // increase the size of the list 
49     int * tempo = new int[fMaxDigit*=2]; 
50     
51     Int_t index ; 
52     
53     for ( index = 0 ; index < fMulDigit ; index++ )
54       tempo[index] = fDigitsList[index] ; 
55     
56     delete fDigitsList ; 
57     fDigitsList = tempo ; 
58   }
59   
60   fDigitsList[fMulDigit] = digit.GetIndexInList()  ; 
61   fMulDigit++ ; 
62   fAmp += digit.GetAmp() ; 
63 }
64
65 //____________________________________________________________________________
66 void AliPMDRecPoint::Copy(AliPMDRecPoint& recp) const
67 {
68   //
69   // Copy *this onto pts
70   //
71   // Copy all first
72   if(this != &recp) {
73     ((TObject*) this)->Copy((TObject&)recp);
74     recp.fAmp = fAmp;
75     recp.fGeom = fGeom;
76     recp.fIndexInList = fIndexInList;
77     recp.fLocPos = fLocPos;
78     recp.fLocPosM = new TMatrix(*fLocPosM);
79     recp.fMaxDigit = fMaxDigit;
80     recp.fMulDigit = fMulDigit;
81     recp.fMaxTrack = fMaxTrack;
82     recp.fMulTrack = fMulTrack;
83     
84     // Duplicate pointed objects
85     recp.fDigitsList = new Int_t[fMulDigit];
86     memcpy(recp.fDigitsList,fDigitsList,fMulDigit*sizeof(Int_t));
87     recp.fTracksList = new Int_t[fMulTrack];
88     memcpy(recp.fTracksList,fTracksList,fMulTrack*sizeof(Int_t));
89   }
90 }
91
92 //____________________________________________________________________________
93 void AliPMDRecPoint::GetCovarianceMatrix(TMatrix & mat)
94 {
95   // returns the covariant matrix for the local position
96   
97   mat = *fLocPosM ; 
98
99 }
100
101 //____________________________________________________________________________
102 void AliPMDRecPoint::GetLocalPosition(TVector3 & pos) const
103 {
104   // returns the position of the cluster in the local reference system of the sub-detector
105
106   pos = fLocPos;
107
108  
109 }
110
111 //____________________________________________________________________________
112 AliPMDRecPoint & AliPMDRecPoint::operator= (const AliPMDRecPoint &recp)
113 {
114   recp.Copy(*this);
115   return (*this);
116 }
117
118
119 //____________________________________________________________________________
120 void AliPMDRecPoint::GetGlobalPosition(TVector3 & gpos, TMatrix & gmat) const
121 {
122   // returns the position of the cluster in the global reference system of ALICE
123   // and the uncertainty on this position
124   
125
126   fGeom->GetGlobal(this, gpos, gmat) ;
127  
128 }
129
130
131
132
133
134
135
136
137