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 **************************************************************************/
17 ///////////////////////////////////////////////////////////////////////////////
19 // Time Projection Chamber //
20 // Comparison macro for reconstructed tracks - ESDs V0s //
22 // marian.ivanov@cern.ch //
37 #include "AliESDtrack.h"
38 #include "AliTPCParam.h"
39 #include "AliTrackReference.h"
40 #include "AliTPCParamSR.h"
42 #include "AliESDfriend.h"
43 #include "AliESDtrack.h"
44 #include "AliTPCseed.h"
45 #include "AliITStrackMI.h"
47 #include "AliESDVertex.h"
48 #include "AliExternalTrackParam.h"
49 #include "AliESDkink.h"
53 #include "AliTreeDraw.h"
54 #include "AliMCInfo.h"
55 #include "AliGenKinkInfo.h"
56 #include "AliGenV0Info.h"
57 #include "AliESDRecKinkInfo.h"
61 ClassImp(AliESDRecKinkInfo)
65 AliESDRecKinkInfo::AliESDRecKinkInfo():
69 fDist1(0), //info about closest distance according closest MC - linear DCA
70 fDist2(0), //info about closest distance parabolic DCA
71 fInvMass(0), //reconstructed invariant mass -
73 fRr(0), // rec position of the vertex
74 fMinR(0), // minimum radius in rphi intersection
75 fDistMinR(0), // distance at minimal radius
76 fPointAngleFi(0), //point angle fi
77 fPointAngleTh(0), //point angle theta
78 fPointAngle(0), //point angle full
79 fStatus(0), //status -tracks
80 fRecStatus(0), //kink -status- 0 - not found 1-good - fake
81 fMultiple(0), // how many times was kink reconstructed
82 fKinkMultiple(0) // how many times was kink reconstructed
85 // Default constructor
87 for (Int_t i=0; i<3; i++) {
93 for (Int_t i=0; i<2; i++) { fLab[i] = 0; }
97 void AliESDRecKinkInfo::Update()
100 if ( (fT1.fTPCOn)&& (fT2.fTPCOn)){
102 // IF BOTH RECONSTRUCTED
103 Float_t distance1,distance2;
104 Double_t xx[3],pp[3];
106 Double_t xd[3],pd[3],signd;
107 Double_t xm[3],pm[3],signm;
108 for (Int_t i=0;i<3;i++){
109 xd[i] = fT2.fTPCinR1[i];
110 pd[i] = fT2.fTPCinP1[i];
111 xm[i] = fT1.fTPCinR1[i];
112 pm[i] = fT1.fTPCinP1[i];
114 signd = fT2.fSign<0 ? -1:1;
115 signm = fT1.fSign<0 ? -1:1;
117 AliHelix dhelix1(xd,pd,signd);
118 dhelix1.GetMomentum(0,pp,0);
119 dhelix1.Evaluate(0,xx);
121 // Double_t x2[3],p2[3];
123 AliHelix mhelix(xm,pm,signm);
125 //find intersection linear
127 Double_t phase[2][2] = { {0,0}, {0,0}};
128 Double_t radius[2] = {0} ;
129 Int_t points = dhelix1.GetRPHIintersections(mhelix, phase, radius,200);
130 Double_t delta1=10000,delta2=10000;
133 fMinR = TMath::Sqrt(radius[0]);
136 fMinR =TMath::Min(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
140 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
141 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
142 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
145 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
146 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
147 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
150 fMinR = TMath::Sqrt(radius[0]);
154 if (radius[0]<radius[1]){
155 fMinR = TMath::Sqrt(radius[0]);
159 fMinR = TMath::Sqrt(radius[1]);
165 distance1 = TMath::Min(delta1,delta2);
167 //find intersection parabolic
169 points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
170 delta1=10000,delta2=10000;
173 dhelix1.ParabolicDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
176 dhelix1.ParabolicDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
179 distance2 = TMath::Min(delta1,delta2);
182 dhelix1.Evaluate(phase[0][0],fXr);
183 dhelix1.GetMomentum(phase[0][0],fPdr);
184 mhelix.GetMomentum(phase[0][1],fPm);
185 dhelix1.GetAngle(phase[0][0],mhelix,phase[0][1],fAngle);
186 fRr = TMath::Sqrt(radius[0]);
189 dhelix1.Evaluate(phase[1][0],fXr);
190 dhelix1.GetMomentum(phase[1][0], fPdr);
191 mhelix.GetMomentum(phase[1][1], fPm);
192 dhelix1.GetAngle(phase[1][0],mhelix,phase[1][1],fAngle);
193 fRr = TMath::Sqrt(radius[1]);
195 fDist1 = TMath::Sqrt(distance1);
196 fDist2 = TMath::Sqrt(distance2);
199 Double_t x,alpha,param[5],cov[15];
201 fT1.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
202 fT1.GetESDtrack()->GetInnerExternalCovariance(cov);
203 AliExternalTrackParam paramm(x,alpha,param,cov);
205 fT2.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
206 fT2.GetESDtrack()->GetInnerExternalCovariance(cov);
207 AliExternalTrackParam paramd(x,alpha,param,cov);
210 kink.Update(¶mm,¶md);
212 Double_t diff = kink.fRr-fRr;
213 Double_t diff2 = kink.fDist2-fDist2;
214 printf("Diff\t%f\t%f\n",diff,diff2);