1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 //-----------------------------------------------------------------------------
19 // Class AliMUONTriggerTrack
20 //---------------------------
21 // Reconstructed Trigger track in ALICE dimuon spectrometer
22 // Note: equivalent to AliMUONTriggerTrack for tracking,
23 // No need for a AliMUONTriggerTrackParam
24 // Author: Philippe Crochet
25 //-----------------------------------------------------------------------------
27 #include "AliMUONTriggerTrack.h"
28 #include "AliMUONTrackReconstructor.h"
30 #include <Riostream.h>
34 ClassImp(AliMUONTriggerTrack)
37 //__________________________________________________________________________
38 AliMUONTriggerTrack::AliMUONTriggerTrack()
49 fHitsPatternInTrigCh(0),
53 AliDebug(5,Form("this=%p",this));
55 //__________________________________________________________________________
56 AliMUONTriggerTrack::AliMUONTriggerTrack(Float_t x11, Float_t y11, Float_t z11, Float_t z21, Float_t slopeX, Float_t slopeY, Int_t loTrgNum, Long_t theGTPattern, Int_t ptCutLevel)
65 fGTPattern(theGTPattern),
66 fPtCutLevel(ptCutLevel),
67 fHitsPatternInTrigCh(0),
70 /// ctor from local trigger output
71 AliDebug(5,Form("this=%p x11=%f y11=%f z11=%f z21=%f slopeX=%f slopeY=%f loTrgNum=%d GTPattern=%ld HitsPatternInTrigCh %i",
72 this,x11,y11,z11,z21,slopeX,slopeY,loTrgNum,theGTPattern,fHitsPatternInTrigCh));
76 //__________________________________________________________________________
77 AliMUONTriggerTrack::~AliMUONTriggerTrack()
80 AliDebug(5,Form("this=%p",this));
87 //__________________________________________________________________________
88 AliMUONTriggerTrack::AliMUONTriggerTrack (const AliMUONTriggerTrack& theMUONTriggerTrack)
89 : TObject(theMUONTriggerTrack),
90 fx11(theMUONTriggerTrack.fx11),
91 fy11(theMUONTriggerTrack.fy11),
92 fz11(theMUONTriggerTrack.fz11),
93 fz21(theMUONTriggerTrack.fz21),
94 fSlopeX(theMUONTriggerTrack.fSlopeX),
95 fSlopeY(theMUONTriggerTrack.fSlopeY),
96 floTrgNum(theMUONTriggerTrack.floTrgNum),
97 fGTPattern(theMUONTriggerTrack.fGTPattern),
98 fPtCutLevel(theMUONTriggerTrack.fPtCutLevel),
99 fHitsPatternInTrigCh(theMUONTriggerTrack.fHitsPatternInTrigCh),
105 if (theMUONTriggerTrack.fCovariances) fCovariances = new TMatrixD(*(theMUONTriggerTrack.fCovariances));
106 AliDebug(5,Form("this=%p copy ctor",this));
110 //__________________________________________________________________________
111 AliMUONTriggerTrack & AliMUONTriggerTrack::operator=(const AliMUONTriggerTrack&
114 /// Assignment operator
116 // check assignement to self
117 if (this == &theMUONTriggerTrack)
120 /// base class assignement
121 TObject::operator=(theMUONTriggerTrack);
123 fx11 = theMUONTriggerTrack.fx11;
124 fy11 = theMUONTriggerTrack.fy11;
125 fz11 = theMUONTriggerTrack.fz11;
126 fz21 = theMUONTriggerTrack.fz21;
127 fSlopeX = theMUONTriggerTrack.fSlopeX;
128 fSlopeY = theMUONTriggerTrack.fSlopeY;
129 floTrgNum = theMUONTriggerTrack.floTrgNum;
130 fGTPattern = theMUONTriggerTrack.fGTPattern;
131 fHitsPatternInTrigCh = theMUONTriggerTrack.fHitsPatternInTrigCh;
133 if (theMUONTriggerTrack.fCovariances) {
134 if (fCovariances) *fCovariances = *(theMUONTriggerTrack.fCovariances);
135 else fCovariances = new TMatrixD(*(theMUONTriggerTrack.fCovariances));
144 //__________________________________________________________________________
145 void AliMUONTriggerTrack::Clear(Option_t* /*opt*/)
152 //__________________________________________________________________________
154 AliMUONTriggerTrack::Print(Option_t* opt) const
157 TString optString(opt);
159 if ( optString.Contains("FULL") ) optString = "PARAM COV";
161 if ( optString.Contains("PARAM"))
162 cout << Form("(X,Y,Z)11=(%7.2f,%7.2f,%7.2f) Z21=%7.2f Slope(X,Y)=(%7.2f,%7.2f) LocalBoard #%3d GlobalTriggerPattern %x HitsPatternInTrigCh %x",
163 fx11,fy11,fz11,fz21,fSlopeX,fSlopeY,floTrgNum,fGTPattern,fHitsPatternInTrigCh) << endl;
165 if ( optString.Contains("COV") ){
166 if ( ! fCovariances ) cout << "Covariances not initialized " << endl;
167 else fCovariances->Print();
171 //__________________________________________________________________________
172 void AliMUONTriggerTrack::SetCovariances(const TMatrixD& covariances)
174 /// Set the covariance matrix
175 if (fCovariances) *fCovariances = covariances;
176 else fCovariances = new TMatrixD(covariances);
179 //__________________________________________________________________________
180 void AliMUONTriggerTrack::SetCovariances(const Double_t matrix[3][3])
182 /// Set the covariance matrix
183 if (fCovariances) fCovariances->SetMatrixArray(&(matrix[0][0]));
184 else fCovariances = new TMatrixD(3,3,&(matrix[0][0]));
187 //__________________________________________________________________________
188 const TMatrixD& AliMUONTriggerTrack::GetCovariances() const
190 /// Return the covariance matrix (create it before if needed)
192 fCovariances = new TMatrixD(3,3);
193 fCovariances->Zero();
195 return *fCovariances;
198 //__________________________________________________________________________
199 Bool_t AliMUONTriggerTrack::Match(AliMUONTriggerTrack &track,
200 Double_t sigmaCut) const
202 /// Try to match this track with the given track. Matching conditions:
203 /// - x, y position and y slope within sigmaCut
205 // Find the track with the covariances correctly set
206 // Extrapolate to the z of the other track
207 Bool_t hasCov1 = ( GetCovariances().NonZeros() != 0 );
208 Bool_t hasCov2 = ( track.GetCovariances().NonZeros() != 0 );
210 const AliMUONTriggerTrack* trackToExtrap = ( hasCov2 ) ? &track : this;
211 const AliMUONTriggerTrack* fixedTrack = ( hasCov2 ) ? this : &track;
213 TMatrixD paramDiff(3,1);
214 Double_t deltaZ = fixedTrack->GetZ11() - trackToExtrap->GetZ11();
215 paramDiff(0,0) = fixedTrack->GetX11() - trackToExtrap->GetX11();
216 paramDiff(1,0) = fixedTrack->GetY11() - ( trackToExtrap->GetY11() + trackToExtrap->GetSlopeY() * deltaZ );
217 paramDiff(2,0) = fixedTrack->GetSlopeY() - trackToExtrap->GetSlopeY();
220 TMatrixD cov1(fixedTrack->GetCovariances());
221 TMatrixD cov2(trackToExtrap->GetCovariances());
223 // Extrapolate covariances to z
225 if ( hasCov1 || hasCov2 ){
229 TMatrixD tmp(trackToExtrap->GetCovariances(),TMatrixD::kMultTranspose,jacob);
230 TMatrixD tmp2(jacob,TMatrixD::kMult,tmp);
235 AliDebug(3, Form("track1 Y11 %f track2 Y11: %f (Z11 %f) -> %f (Z11 %f)", fixedTrack->GetY11(), trackToExtrap->GetY11(), trackToExtrap->GetZ11(), trackToExtrap->GetY11() + trackToExtrap->GetSlopeY() * deltaZ, fixedTrack->GetZ11()));
237 TMatrixD sumCov(cov1,TMatrixD::kPlus,cov2);
238 if (sumCov.Determinant() != 0) {
240 TMatrixD tmp(sumCov,TMatrixD::kMult,paramDiff);
241 TMatrixD chi2M(paramDiff,TMatrixD::kTransposeMult,tmp);
244 AliWarning(" Determinant = 0");
245 Double_t sigma2 = 0.;
246 for (Int_t iVar = 0; iVar < 3; iVar++) {
247 sigma2 = cov1(iVar,iVar) + cov2(iVar,iVar);
248 chi2 += paramDiff(iVar,0) * paramDiff(iVar,0) / sigma2;
252 if ( chi2/3 > sigmaCut * sigmaCut )