]>
Commit | Line | Data |
---|---|---|
df9db588 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-2007, 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 | //------------------------------------------------------------------------- | |
a9255000 | 19 | // AOD cluster base class |
df9db588 | 20 | // Author: Markus Oldenburg, CERN |
21 | //------------------------------------------------------------------------- | |
22 | ||
a9255000 | 23 | #include "AliAODCluster.h" |
df9db588 | 24 | |
a9255000 | 25 | ClassImp(AliAODCluster) |
df9db588 | 26 | |
27 | //______________________________________________________________________________ | |
a9255000 | 28 | AliAODCluster::AliAODCluster() : |
df9db588 | 29 | AliVirtualParticle(), |
30 | fEnergy(0), | |
31 | fChi2(-999.), | |
32 | fID(-999), | |
33 | fLabel(-999), | |
34 | fCovMatrix(NULL), | |
35 | fProdVertex(0x0), | |
36 | fPrimTrack(NULL), | |
37 | fType(kUndef) | |
38 | { | |
39 | // default constructor | |
40 | ||
41 | SetPosition((Float_t*)NULL); | |
42 | SetPID((Float_t*)NULL); | |
43 | } | |
44 | ||
45 | //______________________________________________________________________________ | |
a9255000 | 46 | AliAODCluster::AliAODCluster(Int_t id, |
df9db588 | 47 | Int_t label, |
48 | Double_t energy, | |
49 | Double_t x[3], | |
01afc3fc | 50 | Double_t covMatrix[10], |
14d55e62 | 51 | Double_t pid[9], |
df9db588 | 52 | AliAODVertex *prodVertex, |
53 | AliAODTrack *primTrack, | |
54 | Char_t ttype) : | |
55 | AliVirtualParticle(), | |
56 | fEnergy(energy), | |
57 | fChi2(-999.), | |
58 | fID(id), | |
59 | fLabel(label), | |
60 | fCovMatrix(NULL), | |
61 | fProdVertex(prodVertex), | |
62 | fPrimTrack(primTrack), | |
63 | fType(ttype) | |
64 | { | |
65 | // constructor | |
66 | ||
01afc3fc | 67 | SetPosition(x); |
df9db588 | 68 | if(covMatrix) SetCovMatrix(covMatrix); |
69 | SetPID(pid); | |
70 | ||
71 | } | |
72 | ||
73 | //______________________________________________________________________________ | |
a9255000 | 74 | AliAODCluster::AliAODCluster(Int_t id, |
df9db588 | 75 | Int_t label, |
76 | Float_t energy, | |
77 | Float_t x[3], | |
01afc3fc | 78 | Float_t covMatrix[10], |
14d55e62 | 79 | Float_t pid[9], |
df9db588 | 80 | AliAODVertex *prodVertex, |
81 | AliAODTrack *primTrack, | |
82 | Char_t ttype) : | |
83 | AliVirtualParticle(), | |
84 | fEnergy(energy), | |
85 | fChi2(-999.), | |
86 | fID(id), | |
87 | fLabel(label), | |
88 | fCovMatrix(NULL), | |
89 | fProdVertex(prodVertex), | |
90 | fPrimTrack(primTrack), | |
91 | fType(ttype) | |
92 | { | |
93 | // constructor | |
94 | ||
01afc3fc | 95 | SetPosition(x); |
df9db588 | 96 | if(covMatrix) SetCovMatrix(covMatrix); |
97 | SetPID(pid); | |
98 | ||
99 | } | |
100 | ||
101 | ||
102 | //______________________________________________________________________________ | |
a9255000 | 103 | AliAODCluster::~AliAODCluster() |
df9db588 | 104 | { |
105 | // destructor | |
106 | delete fCovMatrix; | |
107 | } | |
108 | ||
109 | ||
110 | //______________________________________________________________________________ | |
14d55e62 | 111 | AliAODCluster::AliAODCluster(const AliAODCluster& clus) : |
112 | AliVirtualParticle(clus), | |
113 | fEnergy(clus.fEnergy), | |
114 | fChi2(clus.fChi2), | |
115 | fID(clus.fID), | |
116 | fLabel(clus.fLabel), | |
df9db588 | 117 | fCovMatrix(NULL), |
14d55e62 | 118 | fProdVertex(clus.fProdVertex), |
119 | fPrimTrack(clus.fPrimTrack), | |
120 | fType(clus.fType) | |
df9db588 | 121 | { |
122 | // Copy constructor | |
123 | ||
14d55e62 | 124 | clus.GetPosition(fPosition); |
125 | if(clus.fCovMatrix) fCovMatrix=new AliAODRedCov<4>(*clus.fCovMatrix); | |
126 | SetPID(clus.fPID); | |
df9db588 | 127 | |
128 | } | |
129 | ||
130 | //______________________________________________________________________________ | |
14d55e62 | 131 | AliAODCluster& AliAODCluster::operator=(const AliAODCluster& clus) |
df9db588 | 132 | { |
133 | // Assignment operator | |
14d55e62 | 134 | if(this!=&clus) { |
df9db588 | 135 | |
14d55e62 | 136 | AliVirtualParticle::operator=(clus); |
df9db588 | 137 | |
14d55e62 | 138 | clus.GetPosition(fPosition); |
139 | clus.GetPID(fPID); | |
df9db588 | 140 | |
14d55e62 | 141 | fEnergy = clus.fEnergy; |
142 | fChi2 = clus.fChi2; | |
df9db588 | 143 | |
14d55e62 | 144 | fID = clus.fID; |
145 | fLabel = clus.fLabel; | |
df9db588 | 146 | |
147 | delete fCovMatrix; | |
14d55e62 | 148 | if(clus.fCovMatrix) fCovMatrix=new AliAODRedCov<4>(*clus.fCovMatrix); |
df9db588 | 149 | else fCovMatrix=NULL; |
14d55e62 | 150 | fProdVertex = clus.fProdVertex; |
151 | fPrimTrack = clus.fPrimTrack; | |
df9db588 | 152 | |
14d55e62 | 153 | fType = clus.fType; |
df9db588 | 154 | } |
155 | ||
156 | return *this; | |
157 | } | |
158 | ||
159 | //______________________________________________________________________________ | |
a9255000 | 160 | template <class T> void AliAODCluster::SetPosition(const T *x) |
df9db588 | 161 | { |
162 | // set the position | |
163 | ||
164 | if (x) { | |
df9db588 | 165 | fPosition[0] = x[0]; |
166 | fPosition[1] = x[1]; | |
167 | fPosition[2] = x[2]; | |
df9db588 | 168 | } else { |
df9db588 | 169 | |
170 | fPosition[0] = -999.; | |
171 | fPosition[1] = -999.; | |
172 | fPosition[2] = -999.; | |
173 | } | |
174 | } | |
175 | ||
14d55e62 | 176 | //______________________________________________________________________________ |
177 | AliAODCluster::AODCluPID_t AliAODCluster::GetMostProbablePID() const | |
178 | { | |
179 | // Returns the most probable PID array element. | |
180 | ||
181 | Int_t nPID = 9; | |
182 | if (fPID) { | |
183 | AODCluPID_t loc = kUnknown; | |
184 | Double_t max = 0.; | |
185 | Bool_t allTheSame = kTRUE; | |
186 | ||
187 | for (Int_t iPID = 0; iPID < nPID; iPID++) { | |
188 | if (fPID[iPID] >= max) { | |
189 | if (fPID[iPID] > max) { | |
190 | allTheSame = kFALSE; | |
191 | max = fPID[iPID]; | |
192 | loc = (AODCluPID_t)iPID; | |
193 | } else { | |
194 | allTheSame = kTRUE; | |
195 | } | |
196 | } | |
197 | } | |
198 | ||
199 | return allTheSame ? kUnknown : loc; | |
200 | } else { | |
201 | return kUnknown; | |
202 | } | |
203 | } | |
204 | ||
df9db588 | 205 | //______________________________________________________________________________ |
a9255000 | 206 | void AliAODCluster::Print(Option_t* /* option */) const |
df9db588 | 207 | { |
a9255000 | 208 | // prints information about AliAODCluster |
df9db588 | 209 | |
a9255000 | 210 | printf("Object name: %s Cluster type: %s\n", GetName(), GetTitle()); |
df9db588 | 211 | printf(" energy = %f\n", E()); |
212 | printf(" chi2 = %f\n", Chi2()); | |
213 | printf(" PID object: %p\n", PID()); | |
214 | } | |
215 |