Coding conventions - correcting includes
[u/mrichter/AliRoot.git] / MUON / AliMUONTrack.cxx
CommitLineData
a9e2aefa 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
88cb7938 16/* $Id$ */
a9e2aefa 17
3831f268 18///////////////////////////////////////////////////
a9e2aefa 19//
3831f268 20// Reconstructed track
21// in
22// ALICE
23// dimuon
24// spectrometer
25//
26///////////////////////////////////////////////////
a9e2aefa 27
63ed9c6b 28#include "AliMUONTrack.h"
34f1bfa0 29
de2cd600 30#include "AliMUONTrackParam.h"
a9e2aefa 31#include "AliMUONHitForRec.h"
d837040f 32#include "AliMUONConstants.h"
208f139e 33#include "AliMUONTrackExtrap.h"
a9e2aefa 34
63ed9c6b 35#include "AliLog.h"
36
63ed9c6b 37#include <TMath.h>
4889d34c 38#include <Riostream.h>
39#include <TMatrixD.h>
63ed9c6b 40
7945aae7 41/// \cond CLASSIMP
63ed9c6b 42ClassImp(AliMUONTrack) // Class implementation in ROOT context
7945aae7 43/// \endcond
63ed9c6b 44
208f139e 45const Double_t AliMUONTrack::fgkMaxTrackingDistanceBending = 2.;
46const Double_t AliMUONTrack::fgkMaxTrackingDistanceNonBending = 2.;
47
63ed9c6b 48//__________________________________________________________________________
30178c30 49AliMUONTrack::AliMUONTrack()
54d7ba50 50 : TObject(),
54d7ba50 51 fTrackParamAtVertex(),
52 fTrackParamAtHit(0x0),
53 fHitForRecAtHit(0x0),
54d7ba50 54 fNTrackHits(0),
208f139e 55 fExtrapTrackParam(),
56 fFitWithVertex(kFALSE),
57 fVertex(0x0),
54d7ba50 58 fFitFMin(-1.),
59 fMatchTrigger(kFALSE),
c6ba19f7 60 floTrgNum(-1),
54d7ba50 61 fChi2MatchTrigger(0.),
62 fTrackID(0)
d837040f 63{
2457f726 64 /// Default constructor
d837040f 65}
66
a9e2aefa 67 //__________________________________________________________________________
208f139e 68AliMUONTrack::AliMUONTrack(AliMUONHitForRec* hitForRec1, AliMUONHitForRec* hitForRec2)
54d7ba50 69 : TObject(),
54d7ba50 70 fTrackParamAtVertex(),
71 fTrackParamAtHit(0x0),
72 fHitForRecAtHit(0x0),
54d7ba50 73 fNTrackHits(0),
208f139e 74 fExtrapTrackParam(),
75 fFitWithVertex(kFALSE),
76 fVertex(0x0),
54d7ba50 77 fFitFMin(-1.),
78 fMatchTrigger(kFALSE),
c6ba19f7 79 floTrgNum(-1),
54d7ba50 80 fChi2MatchTrigger(0.),
81 fTrackID(0)
a9e2aefa 82{
208f139e 83 /// Constructor from thw hitForRec's
54d7ba50 84
d837040f 85 fTrackParamAtHit = new TClonesArray("AliMUONTrackParam",10);
b8dc484b 86 fHitForRecAtHit = new TClonesArray("AliMUONHitForRec",10);
de2cd600 87
208f139e 88 if (!hitForRec1) return; //AZ
89
90 // Add hits to the track
91 AddTrackParamAtHit(0,hitForRec1);
92 AddTrackParamAtHit(0,hitForRec2);
93
94 // sort TrackParamAtHit according to increasing -Z
95 fTrackParamAtHit->Sort();
96
97 // Set track parameters at first track hit
98 AliMUONTrackParam* trackParamAtFirstHit = (AliMUONTrackParam*) fTrackParamAtHit->First();
99 AliMUONHitForRec* firstHit = trackParamAtFirstHit->GetHitForRecPtr();
100 AliMUONHitForRec* lastHit = ((AliMUONTrackParam*) fTrackParamAtHit->Last())->GetHitForRecPtr();
847cbaef 101 Double_t dZ = firstHit->GetZ() - lastHit->GetZ();
208f139e 102 // Non bending plane
103 Double_t nonBendingCoor = firstHit->GetNonBendingCoor();
104 trackParamAtFirstHit->SetNonBendingCoor(nonBendingCoor);
105 trackParamAtFirstHit->SetNonBendingSlope((nonBendingCoor - lastHit->GetNonBendingCoor()) / dZ);
106 // Bending plane
107 Double_t bendingCoor = firstHit->GetBendingCoor();
108 trackParamAtFirstHit->SetBendingCoor(bendingCoor);
109 Double_t bendingSlope = (bendingCoor - lastHit->GetBendingCoor()) / dZ;
110 trackParamAtFirstHit->SetBendingSlope(bendingSlope);
111 // Inverse bending momentum
847cbaef 112 Double_t bendingImpact = bendingCoor - firstHit->GetZ() * bendingSlope;
208f139e 113 Double_t inverseBendingMomentum = 1. / AliMUONTrackExtrap::GetBendingMomentumFromImpactParam(bendingImpact);
114 trackParamAtFirstHit->SetInverseBendingMomentum(inverseBendingMomentum);
115
116 // Evaluate covariances
117 TMatrixD *paramCov = trackParamAtFirstHit->GetCovariances();
118 (*paramCov) = 0;
119 // Non bending plane
120 (*paramCov)(0,0) = firstHit->GetNonBendingReso2();
121 (*paramCov)(0,1) = firstHit->GetNonBendingReso2() / dZ;
122 (*paramCov)(1,0) = (*paramCov)(0,1);
123 (*paramCov)(1,1) = ( firstHit->GetNonBendingReso2() + lastHit->GetNonBendingReso2() ) / dZ / dZ;
124 // Bending plane
125 (*paramCov)(2,2) = firstHit->GetBendingReso2();
126 (*paramCov)(2,3) = firstHit->GetBendingReso2() / dZ;
127 (*paramCov)(3,2) = (*paramCov)(2,3);
128 (*paramCov)(3,3) = ( firstHit->GetBendingReso2() + lastHit->GetBendingReso2() ) / dZ / dZ;
129 // Inverse bending momentum (50% error)
130 (*paramCov)(4,4) = 0.5*inverseBendingMomentum * 0.5*inverseBendingMomentum;
de2cd600 131
a9e2aefa 132}
133
8429a5e4 134 //__________________________________________________________________________
135AliMUONTrack::~AliMUONTrack()
136{
2457f726 137 /// Destructor
d837040f 138 if (fTrackParamAtHit) {
139 // delete the TClonesArray of pointers to TrackParam
140 delete fTrackParamAtHit;
208f139e 141 fTrackParamAtHit = 0x0;
d837040f 142 }
b8dc484b 143
144 if (fHitForRecAtHit) {
145 // delete the TClonesArray of pointers to HitForRec
146 delete fHitForRecAtHit;
208f139e 147 fHitForRecAtHit = 0x0;
148 }
149
150 if (fVertex) {
151 // delete the vertex used during the tracking procedure
152 delete fVertex;
153 fVertex = 0x0;
b8dc484b 154 }
8429a5e4 155}
156
956019b6 157 //__________________________________________________________________________
30178c30 158AliMUONTrack::AliMUONTrack (const AliMUONTrack& theMUONTrack)
54d7ba50 159 : TObject(theMUONTrack),
54d7ba50 160 fTrackParamAtVertex(theMUONTrack.fTrackParamAtVertex),
161 fTrackParamAtHit(0x0),
162 fHitForRecAtHit(0x0),
54d7ba50 163 fNTrackHits(theMUONTrack.fNTrackHits),
208f139e 164 fExtrapTrackParam(theMUONTrack.fExtrapTrackParam),
165 fFitWithVertex(theMUONTrack.fFitWithVertex),
166 fVertex(0x0),
54d7ba50 167 fFitFMin(theMUONTrack.fFitFMin),
168 fMatchTrigger(theMUONTrack.fMatchTrigger),
c6ba19f7 169 floTrgNum(theMUONTrack.floTrgNum),
54d7ba50 170 fChi2MatchTrigger(theMUONTrack.fChi2MatchTrigger),
171 fTrackID(theMUONTrack.fTrackID)
a9e2aefa 172{
2457f726 173 ///copy constructor
de2cd600 174 Int_t maxIndex = 0;
175
e516b01d 176 // necessary to make a copy of the objects and not only the pointers in TClonesArray.
a8865299 177 if (theMUONTrack.fTrackParamAtHit) {
de2cd600 178 maxIndex = (theMUONTrack.fTrackParamAtHit)->GetEntriesFast();
208f139e 179 fTrackParamAtHit = new TClonesArray("AliMUONTrackParam",maxIndex);
de2cd600 180 for (Int_t index = 0; index < maxIndex; index++) {
181 new ((*fTrackParamAtHit)[index]) AliMUONTrackParam(*(AliMUONTrackParam*)theMUONTrack.fTrackParamAtHit->At(index));
a8865299 182 }
208f139e 183 }
184
b8dc484b 185 // necessary to make a copy of the objects and not only the pointers in TClonesArray.
a8865299 186 if (theMUONTrack.fHitForRecAtHit) {
de2cd600 187 maxIndex = (theMUONTrack.fHitForRecAtHit)->GetEntriesFast();
208f139e 188 fHitForRecAtHit = new TClonesArray("AliMUONHitForRec",maxIndex);
de2cd600 189 for (Int_t index = 0; index < maxIndex; index++) {
190 new ((*fHitForRecAtHit)[index]) AliMUONHitForRec(*(AliMUONHitForRec*)theMUONTrack.fHitForRecAtHit->At(index));
a8865299 191 }
208f139e 192 }
193
194 // copy vertex used during the tracking procedure if any
195 if (theMUONTrack.fVertex) fVertex = new AliMUONHitForRec(*(theMUONTrack.fVertex));
196
a9e2aefa 197}
198
956019b6 199 //__________________________________________________________________________
30178c30 200AliMUONTrack & AliMUONTrack::operator=(const AliMUONTrack& theMUONTrack)
a9e2aefa 201{
2457f726 202 /// Asignment operator
30178c30 203 // check assignement to self
204 if (this == &theMUONTrack)
a9e2aefa 205 return *this;
61adb9bd 206
30178c30 207 // base class assignement
208 TObject::operator=(theMUONTrack);
209
de2cd600 210 fTrackParamAtVertex = theMUONTrack.fTrackParamAtVertex;
e516b01d 211
de2cd600 212 Int_t maxIndex = 0;
213
e516b01d 214 // necessary to make a copy of the objects and not only the pointers in TClonesArray.
a8865299 215 if (theMUONTrack.fTrackParamAtHit) {
208f139e 216 if (fTrackParamAtHit) fTrackParamAtHit->Clear();
217 else fTrackParamAtHit = new TClonesArray("AliMUONTrackParam",10);
de2cd600 218 maxIndex = (theMUONTrack.fTrackParamAtHit)->GetEntriesFast();
219 for (Int_t index = 0; index < maxIndex; index++) {
220 new ((*fTrackParamAtHit)[fTrackParamAtHit->GetEntriesFast()])
221 AliMUONTrackParam(*(AliMUONTrackParam*)(theMUONTrack.fTrackParamAtHit)->At(index));
a8865299 222 }
208f139e 223 } else if (fTrackParamAtHit) {
224 delete fTrackParamAtHit;
225 fTrackParamAtHit = 0x0;
226 }
e516b01d 227
b8dc484b 228 // necessary to make a copy of the objects and not only the pointers in TClonesArray.
de2cd600 229 if (theMUONTrack.fHitForRecAtHit) {
208f139e 230 if (fHitForRecAtHit) fHitForRecAtHit->Clear();
231 else fHitForRecAtHit = new TClonesArray("AliMUONHitForRec",10);
de2cd600 232 maxIndex = (theMUONTrack.fHitForRecAtHit)->GetEntriesFast();
233 for (Int_t index = 0; index < maxIndex; index++) {
234 new ((*fHitForRecAtHit)[fHitForRecAtHit->GetEntriesFast()])
235 AliMUONHitForRec(*(AliMUONHitForRec*)(theMUONTrack.fHitForRecAtHit)->At(index));
a8865299 236 }
208f139e 237 } else if (fHitForRecAtHit) {
238 delete fHitForRecAtHit;
239 fHitForRecAtHit = 0x0;
240 }
241
242 // copy vertex used during the tracking procedure if any.
243 if (theMUONTrack.fVertex) {
244 if (fVertex) *fVertex = *(theMUONTrack.fVertex);
245 else fVertex = new AliMUONHitForRec(*(theMUONTrack.fVertex));
246 } else if (fVertex) {
247 delete fVertex;
248 fVertex = 0x0;
249 }
250
251 fExtrapTrackParam = theMUONTrack.fExtrapTrackParam;
de2cd600 252
30178c30 253 fNTrackHits = theMUONTrack.fNTrackHits;
208f139e 254 fFitWithVertex = theMUONTrack.fFitWithVertex;
30178c30 255 fFitFMin = theMUONTrack.fFitFMin;
30178c30 256 fMatchTrigger = theMUONTrack.fMatchTrigger;
c6ba19f7 257 floTrgNum = theMUONTrack.floTrgNum;
30178c30 258 fChi2MatchTrigger = theMUONTrack.fChi2MatchTrigger;
b8dc484b 259 fTrackID = theMUONTrack.fTrackID;
30178c30 260
61adb9bd 261 return *this;
a9e2aefa 262}
263
8429a5e4 264 //__________________________________________________________________________
de2cd600 265void AliMUONTrack::AddTrackParamAtHit(AliMUONTrackParam *trackParam, AliMUONHitForRec *hitForRec)
8429a5e4 266{
847cbaef 267 /// Add TrackParamAtHit if "trackParam" != NULL
268 /// else create empty TrackParamAtHit and set the z position to the one of "hitForRec" if any
2457f726 269 /// Update link to HitForRec if "hitForRec" != NULL
de2cd600 270 if (!fTrackParamAtHit) {
271 fTrackParamAtHit = new TClonesArray("AliMUONTrackParam",10);
272 fNTrackHits = 0;
956019b6 273 }
2457f726 274 AliMUONTrackParam* trackParamAtHit;
847cbaef 275 if (trackParam) {
276 trackParamAtHit = new ((*fTrackParamAtHit)[fNTrackHits]) AliMUONTrackParam(*trackParam);
277 if (hitForRec) {
278 if (hitForRec->GetZ() != trackParam->GetZ())
279 AliWarning("Added track parameters at a different z position than the one of the attached hit");
280 }
281 } else {
282 trackParamAtHit = new ((*fTrackParamAtHit)[fNTrackHits]) AliMUONTrackParam();
283 if (hitForRec) ((AliMUONTrackParam*) fTrackParamAtHit->UncheckedAt(fNTrackHits))->SetZ(hitForRec->GetZ());
284 }
2457f726 285 if (hitForRec) trackParamAtHit->SetHitForRecPtr(hitForRec);
de2cd600 286 fNTrackHits++;
956019b6 287}
288
289 //__________________________________________________________________________
de2cd600 290void AliMUONTrack::AddHitForRecAtHit(const AliMUONHitForRec *hitForRec)
956019b6 291{
2457f726 292 /// Add hitForRec to the array of hitForRec at hit
de2cd600 293 if (!fHitForRecAtHit)
294 fHitForRecAtHit = new TClonesArray("AliMUONHitForRec",10);
295
296 if (!hitForRec)
297 AliFatal("AliMUONTrack::AddHitForRecAtHit: hitForRec == NULL");
298
299 new ((*fHitForRecAtHit)[fHitForRecAtHit->GetEntriesFast()]) AliMUONHitForRec(*hitForRec);
04b5ea16 300}
301
b8dc484b 302 //__________________________________________________________________________
208f139e 303void AliMUONTrack::SetVertex(AliMUONHitForRec* vertex)
304{
305 /// Set the vertex used during the tracking procedure
306 if (!fVertex) fVertex = new AliMUONHitForRec(*vertex);
307 else *fVertex = *vertex;
308}
309
310 //__________________________________________________________________________
311Int_t AliMUONTrack::HitsInCommon(AliMUONTrack* track) const
312{
313 /// Returns the number of hits in common between the current track ("this")
314 /// and the track pointed to by "track".
315 Int_t hitsInCommon = 0;
316 AliMUONTrackParam *trackParamAtHit1, *trackParamAtHit2;
317 // Loop over hits of first track
318 trackParamAtHit1 = (AliMUONTrackParam*) this->fTrackParamAtHit->First();
319 while (trackParamAtHit1) {
320 // Loop over hits of second track
321 trackParamAtHit2 = (AliMUONTrackParam*) track->fTrackParamAtHit->First();
322 while (trackParamAtHit2) {
323 // Increment "hitsInCommon" if both TrackParamAtHits point to the same HitForRec
324 if ((trackParamAtHit1->GetHitForRecPtr()) == (trackParamAtHit2->GetHitForRecPtr())) {
325 hitsInCommon++;
326 break;
327 }
328 trackParamAtHit2 = (AliMUONTrackParam*) track->fTrackParamAtHit->After(trackParamAtHit2);
329 } // trackParamAtHit2
330 trackParamAtHit1 = (AliMUONTrackParam*) this->fTrackParamAtHit->After(trackParamAtHit1);
331 } // trackParamAtHit1
332 return hitsInCommon;
333}
334
335 //__________________________________________________________________________
336Bool_t* AliMUONTrack::CompatibleTrack(AliMUONTrack * track, Double_t sigma2Cut) const
b8dc484b 337{
2457f726 338 /// Return kTRUE/kFALSE for each chamber if hit is compatible or not
b8dc484b 339 TClonesArray *hitArray, *thisHitArray;
340 AliMUONHitForRec *hit, *thisHit;
341 Int_t chamberNumber;
342 Float_t deltaZ;
343 Float_t deltaZMax = 1.; // 1 cm
344 Float_t chi2 = 0;
345 Bool_t *nCompHit = new Bool_t[AliMUONConstants::NTrackingCh()];
346
347 for ( Int_t ch = 0; ch < AliMUONConstants::NTrackingCh(); ch++) {
348 nCompHit[ch] = kFALSE;
349 }
350
351 thisHitArray = this->GetHitForRecAtHit();
352
208f139e 353 hitArray = track->GetHitForRecAtHit();
b8dc484b 354
355 for (Int_t iHthis = 0; iHthis < thisHitArray->GetEntriesFast(); iHthis++) {
356 thisHit = (AliMUONHitForRec*) thisHitArray->At(iHthis);
357 chamberNumber = thisHit->GetChamberNumber();
358 if (chamberNumber < 0 || chamberNumber > AliMUONConstants::NTrackingCh()) continue;
359 nCompHit[chamberNumber] = kFALSE;
360 for (Int_t iH = 0; iH < hitArray->GetEntriesFast(); iH++) {
361 hit = (AliMUONHitForRec*) hitArray->At(iH);
362 deltaZ = TMath::Abs(thisHit->GetZ() - hit->GetZ());
208f139e 363 chi2 = thisHit->NormalizedChi2WithHitForRec(hit,sigma2Cut); // set cut to 4 sigmas
b8dc484b 364 if (chi2 < 3. && deltaZ < deltaZMax) {
365 nCompHit[chamberNumber] = kTRUE;
366 break;
367 }
368 }
369 }
370
371 return nCompHit;
372}
a9e2aefa 373
8429a5e4 374 //__________________________________________________________________________
208f139e 375Double_t AliMUONTrack::TryOneHitForRec(AliMUONHitForRec* hitForRec)
8429a5e4 376{
208f139e 377/// Test the compatibility between the track and the hitForRec:
378/// return the corresponding Chi2
379
380 // Get track parameters and their covariances at the z position of hitForRec
381 AliMUONTrackParam extrapTrackParam(fExtrapTrackParam);
382 AliMUONTrackExtrap::ExtrapToZCov(&extrapTrackParam, hitForRec->GetZ());
383
384 // Set differences between trackParam and hitForRec in the bending and non bending directions
385 TMatrixD dPos(2,1);
386 dPos(0,0) = hitForRec->GetNonBendingCoor() - extrapTrackParam.GetNonBendingCoor();
387 dPos(1,0) = hitForRec->GetBendingCoor() - extrapTrackParam.GetBendingCoor();
388
389 // quick test of hitForRec compatibility within a wide road of x*y = 10*1 cm2 to save computing time
390 if (TMath::Abs(dPos(0,0)) > fgkMaxTrackingDistanceNonBending ||
391 TMath::Abs(dPos(1,0)) > fgkMaxTrackingDistanceBending) return 1.e10;
392
393 // Set the error matrix from trackParam covariances and hitForRec resolution
394 TMatrixD* paramCov = extrapTrackParam.GetCovariances();
395 TMatrixD error(2,2);
396 error(0,0) = (*paramCov)(0,0) + hitForRec->GetNonBendingReso2();
397 error(0,1) = (*paramCov)(0,2);
398 error(1,0) = (*paramCov)(2,0);
399 error(1,1) = (*paramCov)(2,2) + hitForRec->GetBendingReso2();
400
401 // Invert the error matrix for Chi2 calculation
402 if (error.Determinant() != 0) {
403 error.Invert();
404 } else {
405 AliWarning(" Determinant error=0");
406 return 1.e10;
407 }
408
409 // Compute the Chi2 value
410 TMatrixD tmp(error,TMatrixD::kMult,dPos);
411 TMatrixD result(dPos,TMatrixD::kTransposeMult,tmp);
412
413 return result(0,0);
414
415}
416
417 //__________________________________________________________________________
418Double_t AliMUONTrack::TryTwoHitForRec(AliMUONHitForRec* hitForRec1, AliMUONHitForRec* hitForRec2)
419{
420/// Test the compatibility between the track and the 2 hitForRec together:
421/// return the corresponding Chi2 accounting for covariances between the 2 hitForRec
422
423 // Get track parameters and their covariances at the z position of the first hitForRec
424 AliMUONTrackParam extrapTrackParam1(fExtrapTrackParam);
425 AliMUONTrackExtrap::ExtrapToZCov(&extrapTrackParam1, hitForRec1->GetZ());
426
427 // Get track parameters at second hitForRec
428 AliMUONTrackParam extrapTrackParam2(extrapTrackParam1);
429 AliMUONTrackExtrap::ExtrapToZ(&extrapTrackParam2, hitForRec2->GetZ());
430
431 // Set differences between track and the 2 hitForRec in the bending and non bending directions
432 TMatrixD dPos(4,1);
433 dPos(0,0) = hitForRec1->GetNonBendingCoor() - extrapTrackParam1.GetNonBendingCoor();
434 dPos(1,0) = hitForRec1->GetBendingCoor() - extrapTrackParam1.GetBendingCoor();
435 dPos(2,0) = hitForRec2->GetNonBendingCoor() - extrapTrackParam2.GetNonBendingCoor();
436 dPos(3,0) = hitForRec2->GetBendingCoor() - extrapTrackParam2.GetBendingCoor();
437
438 // quick tests of hitForRec compatibility within a wide road of x*y = 1*1 cm2 to save computing time
439 if (TMath::Abs(dPos(0,0)) > fgkMaxTrackingDistanceNonBending ||
440 TMath::Abs(dPos(1,0)) > fgkMaxTrackingDistanceBending ||
441 TMath::Abs(dPos(2,0)) > fgkMaxTrackingDistanceNonBending ||
442 TMath::Abs(dPos(3,0)) > fgkMaxTrackingDistanceBending) return 1.e10;
443
444 // Calculate the error matrix from the track parameter covariances at first hitForRec
445 TMatrixD error(4,4);
446 error = 0.;
447 if (extrapTrackParam1.CovariancesExist()) {
448 // Get the pointer to the parameter covariance matrix at first hitForRec
449 TMatrixD* paramCov = extrapTrackParam1.GetCovariances();
450
451 // Save track parameters at first hitForRec
452 AliMUONTrackParam extrapTrackParam1Save(extrapTrackParam1);
453 Double_t nonBendingCoor1 = extrapTrackParam1Save.GetNonBendingCoor();
454 Double_t nonBendingSlope1 = extrapTrackParam1Save.GetNonBendingSlope();
455 Double_t bendingCoor1 = extrapTrackParam1Save.GetBendingCoor();
456 Double_t bendingSlope1 = extrapTrackParam1Save.GetBendingSlope();
457 Double_t inverseBendingMomentum1 = extrapTrackParam1Save.GetInverseBendingMomentum();
458 Double_t z1 = extrapTrackParam1Save.GetZ();
459
460 // Save track coordinates at second hitForRec
461 Double_t nonBendingCoor2 = extrapTrackParam2.GetNonBendingCoor();
462 Double_t bendingCoor2 = extrapTrackParam2.GetBendingCoor();
463
464 // Calculate the jacobian related to the transformation between track parameters
465 // at first hitForRec and track coordinates at the 2 hitForRec z-position
466 TMatrixD jacob(4,5);
467 jacob = 0.;
468 // first derivative at the first hitForRec:
469 jacob(0,0) = 1.; // dx1/dx
470 jacob(1,2) = 1.; // dy1/dy
471 // first derivative at the second hitForRec:
472 Double_t dParam[5];
473 for (Int_t i=0; i<5; i++) {
474 // Skip jacobian calculation for parameters with no associated error
475 if ((*paramCov)(i,i) == 0.) continue;
476 // Small variation of parameter i only
477 for (Int_t j=0; j<5; j++) {
478 if (j==i) {
479 dParam[j] = TMath::Sqrt((*paramCov)(i,i));
480 if (j == 4) dParam[j] *= TMath::Sign(1.,-inverseBendingMomentum1); // variation always in the same direction
481 } else dParam[j] = 0.;
482 }
483 // Set new track parameters at first hitForRec
484 extrapTrackParam1Save.SetNonBendingCoor (nonBendingCoor1 + dParam[0]);
485 extrapTrackParam1Save.SetNonBendingSlope (nonBendingSlope1 + dParam[1]);
486 extrapTrackParam1Save.SetBendingCoor (bendingCoor1 + dParam[2]);
487 extrapTrackParam1Save.SetBendingSlope (bendingSlope1 + dParam[3]);
488 extrapTrackParam1Save.SetInverseBendingMomentum(inverseBendingMomentum1 + dParam[4]);
489 extrapTrackParam1Save.SetZ (z1);
490 // Extrapolate new track parameters to the z position of the second hitForRec
491 AliMUONTrackExtrap::ExtrapToZ(&extrapTrackParam1Save,hitForRec2->GetZ());
492 // Calculate the jacobian
493 jacob(2,i) = (extrapTrackParam1Save.GetNonBendingCoor() - nonBendingCoor2) / dParam[i]; // dx2/dParami
494 jacob(3,i) = (extrapTrackParam1Save.GetBendingCoor() - bendingCoor2 ) / dParam[i]; // dy2/dParami
495 }
496
497 // Calculate the error matrix
498 TMatrixD tmp((*paramCov),TMatrixD::kMultTranspose,jacob);
499 error = TMatrixD(jacob,TMatrixD::kMult,tmp);
500 }
501
502 // Add hitForRec resolution to the error matrix
503 error(0,0) += hitForRec1->GetNonBendingReso2();
504 error(1,1) += hitForRec1->GetBendingReso2();
505 error(2,2) += hitForRec2->GetNonBendingReso2();
506 error(3,3) += hitForRec2->GetBendingReso2();
507
508 // invert the error matrix for Chi2 calculation
509 if (error.Determinant() != 0) {
510 error.Invert();
511 } else {
512 AliWarning(" Determinant error=0");
513 return 1.e10;
514 }
515
516 // Compute the Chi2 value
517 TMatrixD tmp2(error,TMatrixD::kMult,dPos);
518 TMatrixD result(dPos,TMatrixD::kTransposeMult,tmp2);
519
520 return result(0,0);
521
8429a5e4 522}
523
d837040f 524 //__________________________________________________________________________
de2cd600 525void AliMUONTrack::RecursiveDump(void) const
a9e2aefa 526{
2457f726 527 /// Recursive dump of AliMUONTrack, i.e. with dump of TrackParamAtHit's and attached HitForRec's
de2cd600 528 AliMUONTrackParam *trackParamAtHit;
529 AliMUONHitForRec *hitForRec;
530 cout << "Recursive dump of Track: " << this << endl;
531 // Track
532 this->Dump();
533 for (Int_t trackHitIndex = 0; trackHitIndex < fNTrackHits; trackHitIndex++) {
534 trackParamAtHit = (AliMUONTrackParam*) ((*fTrackParamAtHit)[trackHitIndex]);
535 // TrackHit
536 cout << "TrackParamAtHit: " << trackParamAtHit << " (index: " << trackHitIndex << ")" << endl;
537 trackParamAtHit->Dump();
538 hitForRec = trackParamAtHit->GetHitForRecPtr();
539 // HitForRec
540 cout << "HitForRec: " << hitForRec << endl;
541 hitForRec->Dump();
a9e2aefa 542 }
de2cd600 543 return;
a9e2aefa 544}
04b5ea16 545
6464217e 546//_____________________________________________-
547void AliMUONTrack::Print(Option_t* opt) const
548{
2457f726 549 /// Printing Track information
550 /// "full" option for printing all the information about the track
6464217e 551 TString sopt(opt);
552 sopt.ToUpper();
553
554 if ( sopt.Contains("FULL") ) {
555 cout << "<AliMUONTrack> No.Clusters=" << setw(2) << GetNTrackHits() <<
556 // ", Bending P="<< setw(8) << setprecision(5) << 1./GetInverseBendingMomentum() <<
557 //", NonBendSlope=" << setw(8) << setprecision(5) << GetNonBendingSlope()*180./TMath::Pi() <<
558 //", BendSlope=" << setw(8) << setprecision(5) << GetBendingSlope()*180./TMath::Pi() <<
559 ", Match2Trig=" << setw(1) << GetMatchTrigger() <<
c6ba19f7 560 ", LoTrgNum=" << setw(3) << GetLoTrgNum() <<
6464217e 561 ", Chi2-tracking-trigger=" << setw(8) << setprecision(5) << GetChi2MatchTrigger() << endl ;
562 GetTrackParamAtHit()->First()->Print("full");
563 }
564 else {
565 cout << "<AliMUONTrack>";
566 GetTrackParamAtHit()->First()->Print("");
567
568 }
569
570}