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"
46 #include "AliTRDtrack.h"
48 #include "AliESDVertex.h"
49 #include "AliExternalTrackParam.h"
50 #include "AliESDkink.h"
54 #include "AliTreeDraw.h"
55 #include "AliGenInfo.h"
56 #include "AliRecInfo.h"
60 ClassImp(AliESDRecInfo)
61 ClassImp(AliESDRecV0Info)
62 ClassImp(AliESDRecKinkInfo)
67 AliTPCParam * GetTPCParam(){
68 AliTPCParamSR * par = new AliTPCParamSR;
76 AliESDRecInfo::AliESDRecInfo():
77 fITSOn(0), // ITS refitted inward
78 fTRDOn(0), // ITS refitted inward
79 fDeltaP(0), //delta of momenta
81 fReconstructed(0), //flag if track was reconstructed
82 fFake(0), // fake track
83 fMultiple(0), // number of reconstructions
84 fTPCOn(0), // TPC refitted inward
85 fBestTOFmatch(0), //best matching between times
86 fESDtrack(0), // esd track
87 fTrackF(0), // friend track
88 fTPCtrack(0), // tpc track
89 fITStrack(0), // its track
90 fTRDtrack(0) // trd track
93 // default constructor
98 AliESDRecInfo::AliESDRecInfo(const AliESDRecInfo& recinfo):
104 memcpy(this,&recinfo, sizeof(recinfo));
105 fESDtrack=0; fTrackF=0; fTPCtrack=0;fITStrack=0;fTRDtrack=0;
106 SetESDtrack(recinfo.GetESDtrack());
110 AliESDRecInfo::~AliESDRecInfo()
116 if (fESDtrack) { delete fESDtrack; fESDtrack=0;}
117 if (fTrackF) { delete fTrackF; fTrackF=0;}
118 if (fTPCtrack) { delete fTPCtrack; fTPCtrack=0;}
119 if (fITStrack) { delete fITStrack; fITStrack=0;}
120 if (fTRDtrack) { delete fTRDtrack; fTRDtrack=0;}
126 void AliESDRecInfo::Reset()
134 if (fESDtrack) { delete fESDtrack; fESDtrack=0;}
135 if (fTrackF) { delete fTrackF; fTrackF=0;}
136 if (fTPCtrack) { delete fTPCtrack; fTPCtrack=0;}
137 if (fITStrack) { delete fITStrack; fITStrack=0;}
138 if (fTRDtrack) { delete fTRDtrack; fTRDtrack=0;}
141 void AliESDRecInfo::SetESDtrack(const AliESDtrack *track){
145 if (fESDtrack) delete fESDtrack;
146 fESDtrack = (AliESDtrack*)track->Clone();
147 if (0 &&track->GetFriendTrack()){
148 if (fTrackF) delete fTrackF;
149 fTrackF = (AliESDfriendTrack*)track->GetFriendTrack()->Clone();
150 if (fTrackF->GetCalibObject(0)){
151 if (fTPCtrack) delete fTPCtrack;
152 fTPCtrack = (AliTPCseed*)fTrackF->GetCalibObject(0)->Clone();
158 void AliESDRecInfo::UpdatePoints(AliESDtrack*track)
162 Int_t iclusters[200];
163 Float_t density[160];
164 for (Int_t i=0;i<160;i++) density[i]=-1.;
168 if (fTPCPoints[0]<fTPCPoints[1]) return;
169 // Int_t nclusters=track->GetTPCclusters(iclusters);
175 for (Int_t i=0;i<160;i++){
176 Int_t last = i-range;
177 if (nall<range) nall++;
179 if (iclusters[last]>0&& (iclusters[last]&0x8000)==0) ngood--;
180 if (iclusters[last]==-1) undeff--;
182 if (iclusters[i]>0&& (iclusters[i]&0x8000)==0) ngood++;
183 if (iclusters[i]==-1) undeff++;
184 if (nall==range &&undeff<range/2) density[i-range/2] = Float_t(ngood)/Float_t(nall-undeff);
188 for (Int_t i=0;i<160;i++){
189 if (density[i]<0) continue;
190 if (density[i]>maxdens){
197 fTPCPoints[3] = maxdens;
198 fTPCPoints[1] = indexmax;
201 for (Int_t i=indexmax;i<160;i++){
202 if (density[i]<0) continue;
203 if (density[i]<maxdens/2.) {
210 for (Int_t i=indexmax;i>0;i--){
211 if (density[i]<0) continue;
212 if (density[i]<maxdens/2.) {
218 // Density at the last 30 padrows
223 for (Int_t i=159;i>0;i--){
224 if (iclusters[i]==-1) continue; //dead zone
226 if (iclusters[i]>0) ngood++;
229 fTPCPoints[4] = Float_t(ngood)/Float_t(nall);
231 if ((track->GetStatus()&AliESDtrack::kITSrefit)>0) fTPCPoints[0]=-1;
238 void AliESDRecInfo::Update(AliMCInfo* info,AliTPCParam * /*par*/, Bool_t reconstructed)
242 //calculates derived variables
245 UpdatePoints(fESDtrack);
247 AliTrackReference * ref = &(info->fTrackRef);
248 fTPCinR0[0] = info->fTrackRef.X();
249 fTPCinR0[1] = info->fTrackRef.Y();
250 fTPCinR0[2] = info->fTrackRef.Z();
251 fTPCinR0[3] = TMath::Sqrt(fTPCinR0[0]*fTPCinR0[0]+fTPCinR0[1]*fTPCinR0[1]);
252 fTPCinR0[4] = TMath::ATan2(fTPCinR0[1],fTPCinR0[0]);
254 fTPCinP0[0] = ref->Px();
255 fTPCinP0[1] = ref->Py();
256 fTPCinP0[2] = ref->Pz();
257 fTPCinP0[3] = ref->Pt();
258 fTPCinP0[4] = ref->P();
259 fDeltaP = (ref->P()-info->fParticle.P())/info->fParticle.P();
262 if (fTPCinP0[3]>0.0000001){
264 fTPCAngle0[0] = TMath::ATan2(fTPCinP0[1],fTPCinP0[0]);
265 fTPCAngle0[1] = TMath::ATan(fTPCinP0[2]/fTPCinP0[3]);
269 fITSinP0[0]=info->fParticle.Px();
270 fITSinP0[1]=info->fParticle.Py();
271 fITSinP0[2]=info->fParticle.Pz();
272 fITSinP0[3]=info->fParticle.Pt();
274 fITSinR0[0]=info->fParticle.Vx();
275 fITSinR0[1]=info->fParticle.Vy();
276 fITSinR0[2]=info->fParticle.Vz();
277 fITSinR0[3] = TMath::Sqrt(fITSinR0[0]*fITSinR0[0]+fITSinR0[1]*fITSinR0[1]);
278 fITSinR0[4] = TMath::ATan2(fITSinR0[1],fITSinR0[0]);
281 if (fITSinP0[3]>0.0000001){
282 fITSAngle0[0] = TMath::ATan2(fITSinP0[1],fITSinP0[0]);
283 fITSAngle0[1] = TMath::ATan(fITSinP0[2]/fITSinP0[3]);
286 for (Int_t i=0;i<4;i++) fStatus[i] =0;
287 fReconstructed = kFALSE;
291 if (reconstructed==kFALSE) return;
293 fLabels[0] = info->fLabel;
294 fLabels[1] = info->fPrimPart;
295 fReconstructed = kTRUE;
296 fTPCOn = ((fESDtrack->GetStatus()&AliESDtrack::kTPCrefit)>0) ? kTRUE : kFALSE;
297 fITSOn = ((fESDtrack->GetStatus()&AliESDtrack::kITSrefit)>0) ? kTRUE : kFALSE;
298 fTRDOn = ((fESDtrack->GetStatus()&AliESDtrack::kTRDrefit)>0) ? kTRUE : kFALSE;
301 if ((fESDtrack->GetStatus()&AliESDtrack::kTPCrefit)>0){
305 if ((fESDtrack->GetStatus()&AliESDtrack::kTPCout)>0){
309 if ((fESDtrack->GetStatus()&AliESDtrack::kTPCin)>0)
314 if ((fESDtrack->GetStatus()&AliESDtrack::kITSout)>0){
318 if ((fESDtrack->GetStatus()&AliESDtrack::kITSrefit)>0){
328 if ((fESDtrack->GetStatus()&AliESDtrack::kTRDrefit)>0){
332 if ((fESDtrack->GetStatus()&AliESDtrack::kTRDout)>0){
336 if ((fESDtrack->GetStatus()&AliESDtrack::kTRDStop)>0){
343 if (((fESDtrack->GetStatus()&AliESDtrack::kTOFout)>0)){
347 fESDtrack->GetIntegratedTimes(times);
348 for (Int_t i=0;i<5;i++){
349 if ( TMath::Abs(fESDtrack->GetTOFsignal()-times[i]) <TMath::Abs(fBestTOFmatch) ){
350 fBestTOFmatch = fESDtrack->GetTOFsignal()-times[i];
354 fESDtrack->GetTOFLabel(toflabel);
355 Bool_t toffake=kTRUE;
356 Bool_t tofdaughter=kFALSE;
357 for (Int_t i=0;i<3;i++){
358 if (toflabel[i]<0) continue;
359 if (toflabel[i]== TMath::Abs(fESDtrack->GetLabel())) toffake=kFALSE;
360 if (toflabel[i]==info->fParticle.GetDaughter(0) || (toflabel[i]==info->fParticle.GetDaughter(1))) tofdaughter=kTRUE; // decay product of original particle
363 if (toffake) fStatus[3] =3; //total fake
364 if (tofdaughter) fStatus[3]=2; //fake because of decay
370 if (fStatus[1]>0 &&info->fNTPCRef>0&&TMath::Abs(fTPCinP0[3])>0.0001){
372 fESDtrack->GetInnerXYZ(fTPCinR1);
373 fTPCinR1[3] = TMath::Sqrt(fTPCinR1[0]*fTPCinR1[0]+fTPCinR1[1]*fTPCinR1[1]);
374 fTPCinR1[4] = TMath::ATan2(fTPCinR1[1],fTPCinR1[0]);
375 fESDtrack->GetInnerPxPyPz(fTPCinP1);
376 fTPCinP1[3] = TMath::Sqrt(fTPCinP1[0]*fTPCinP1[0]+fTPCinP1[1]*fTPCinP1[1]);
377 fTPCinP1[4] = TMath::Sqrt(fTPCinP1[3]*fTPCinP1[3]+fTPCinP1[2]*fTPCinP1[2]);
380 if (fTPCinP1[3]>0.000000000000001){
381 fTPCAngle1[0] = TMath::ATan2(fTPCinP1[1],fTPCinP1[0]);
382 fTPCAngle1[1] = TMath::ATan(fTPCinP1[2]/fTPCinP1[3]);
384 Double_t cov[15], param[5],x, alpha;
385 fESDtrack->GetInnerExternalCovariance(cov);
386 fESDtrack->GetInnerExternalParameters(alpha, x,param);
389 fTPCDelta[0] = (fTPCinR0[4]-fTPCinR1[4])*fTPCinR1[3]; //delta rfi
390 fTPCPools[0] = fTPCDelta[0]/TMath::Sqrt(cov[0]);
391 fTPCDelta[1] = (fTPCinR0[2]-fTPCinR1[2]); //delta z
392 fTPCPools[1] = fTPCDelta[1]/TMath::Sqrt(cov[2]);
393 fTPCDelta[2] = (fTPCAngle0[0]-fTPCAngle1[0]);
394 fTPCPools[2] = fTPCDelta[2]/TMath::Sqrt(cov[5]);
395 fTPCDelta[3] = (TMath::Tan(fTPCAngle0[1])-TMath::Tan(fTPCAngle1[1]));
396 fTPCPools[3] = fTPCDelta[3]/TMath::Sqrt(cov[9]);
397 fTPCDelta[4] = (fTPCinP0[3]-fTPCinP1[3]);
398 Double_t sign = (param[4]>0)? 1.:-1;
400 fTPCPools[4] = sign*(1./fTPCinP0[3]-1./fTPCinP1[3])/TMath::Sqrt(TMath::Abs(cov[14]));
405 fESDtrack->GetExternalParameters(x,param);
406 // fESDtrack->GetConstrainedExternalParameters(x,param);
408 fESDtrack->GetExternalCovariance(cov);
409 //fESDtrack->GetConstrainedExternalCovariance(cov);
410 if (TMath::Abs(param[4])<0.0000000001) return;
412 fESDtrack->GetXYZ(fITSinR1);
413 fESDtrack->GetPxPyPz(fITSinP1);
414 fITSinP1[3] = TMath::Sqrt(fITSinP1[0]*fITSinP1[0]+fITSinP1[1]*fITSinP1[1]);
416 fITSinR1[3] = TMath::Sqrt(fITSinR1[0]*fITSinR1[0]+fITSinR1[1]*fITSinR1[1]);
417 fITSinR1[4] = TMath::ATan2(fITSinR1[1],fITSinR1[0]);
420 if (fITSinP1[3]>0.0000001){
421 fITSAngle1[0] = TMath::ATan2(fITSinP1[1],fITSinP1[0]);
422 fITSAngle1[1] = TMath::ATan(fITSinP1[2]/fITSinP1[3]);
426 fITSDelta[0] = (fITSinR0[4]-fITSinR1[4])*fITSinR1[3]; //delta rfi
427 fITSPools[0] = fITSDelta[0]/TMath::Sqrt(cov[0]);
428 fITSDelta[1] = (fITSinR0[2]-fITSinR1[2]); //delta z
429 fITSPools[1] = fITSDelta[1]/TMath::Sqrt(cov[2]);
430 fITSDelta[2] = (fITSAngle0[0]-fITSAngle1[0]);
431 fITSPools[2] = fITSDelta[2]/TMath::Sqrt(cov[5]);
432 fITSDelta[3] = (TMath::Tan(fITSAngle0[1])-TMath::Tan(fITSAngle1[1]));
433 fITSPools[3] = fITSDelta[3]/TMath::Sqrt(cov[9]);
434 fITSDelta[4] = (fITSinP0[3]-fITSinP1[3]);
435 Double_t sign = (param[4]>0) ? 1:-1;
437 fITSPools[4] = sign*(1./fITSinP0[3]-1./fITSinP1[3])/TMath::Sqrt(cov[14]);
443 void AliESDRecV0Info::Update(Float_t vertex[3])
446 if ( (fT1.fStatus[1]>0)&& (fT2.fStatus[1]>0)){
447 Float_t distance1,distance2;
448 Double_t xx[3],pp[3];
450 Double_t xd[3],pd[3],signd;
451 Double_t xm[3],pm[3],signm;
454 if (fT1.fITSOn&&fT2.fITSOn){
455 for (Int_t i=0;i<3;i++){
456 xd[i] = fT2.fITSinR1[i];
457 pd[i] = fT2.fITSinP1[i];
458 xm[i] = fT1.fITSinR1[i];
459 pm[i] = fT1.fITSinP1[i];
464 for (Int_t i=0;i<3;i++){
465 xd[i] = fT2.fTPCinR1[i];
466 pd[i] = fT2.fTPCinP1[i];
467 xm[i] = fT1.fTPCinR1[i];
468 pm[i] = fT1.fTPCinP1[i];
473 signd = fT2.fSign<0 ? -1:1;
474 signm = fT1.fSign<0 ? -1:1;
476 AliHelix dhelix1(xd,pd,signd);
477 dhelix1.GetMomentum(0,pp,0);
478 dhelix1.Evaluate(0,xx);
480 // Double_t x2[3],p2[3];
482 AliHelix mhelix(xm,pm,signm);
484 //find intersection linear
486 Double_t phase[2][2],radius[2];
487 Int_t points = dhelix1.GetRPHIintersections(mhelix, phase, radius,200);
488 Double_t delta1=10000,delta2=10000;
491 fRs[0] = TMath::Sqrt(radius[0]);
492 fRs[1] = TMath::Sqrt(radius[0]);
495 fRs[0] =TMath::Min(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
496 fRs[1] =TMath::Max(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
500 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
501 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
502 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
505 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
506 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
507 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
510 fRs[0] = TMath::Sqrt(radius[0]);
511 fRs[1] = TMath::Sqrt(radius[0]);
515 if (radius[0]<radius[1]){
516 fRs[0] = TMath::Sqrt(radius[0]);
517 fRs[1] = TMath::Sqrt(radius[1]);
521 fRs[0] = TMath::Sqrt(radius[1]);
522 fRs[1] = TMath::Sqrt(radius[0]);
528 distance1 = TMath::Min(delta1,delta2);
530 //find intersection parabolic
532 points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
533 delta1=10000,delta2=10000;
536 dhelix1.ParabolicDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
539 dhelix1.ParabolicDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
542 distance2 = TMath::Min(delta1,delta2);
543 if (distance2>100) fDist2 =100;
547 dhelix1.Evaluate(phase[0][0],fXr);
548 dhelix1.GetMomentum(phase[0][0],fPdr);
549 mhelix.GetMomentum(phase[0][1],fPm);
550 dhelix1.GetAngle(phase[0][0],mhelix,phase[0][1],fAngle);
551 fRr = TMath::Sqrt(radius[0]);
554 dhelix1.Evaluate(phase[1][0],fXr);
555 dhelix1.GetMomentum(phase[1][0], fPdr);
556 mhelix.GetMomentum(phase[1][1], fPm);
557 dhelix1.GetAngle(phase[1][0],mhelix,phase[1][1],fAngle);
558 fRr = TMath::Sqrt(radius[1]);
560 fDist1 = TMath::Sqrt(distance1);
561 fDist2 = TMath::Sqrt(distance2);
564 Double_t x,alpha,param[5],cov[15];
566 fT1.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
567 fT1.GetESDtrack()->GetInnerExternalCovariance(cov);
568 AliExternalTrackParam paramm(x,alpha,param,cov);
570 fT2.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
571 fT2.GetESDtrack()->GetInnerExternalCovariance(cov);
572 AliExternalTrackParam paramd(x,alpha,param,cov);
577 Float_t v[3] = {fXr[0]-vertex[0],fXr[1]-vertex[1],fXr[2]-vertex[2]};
578 Float_t p[3] = {fPdr[0]+fPm[0], fPdr[1]+fPm[1],fPdr[2]+fPm[2]};
580 Float_t vnorm2 = v[0]*v[0]+v[1]*v[1];
581 Float_t vnorm3 = TMath::Sqrt(v[2]*v[2]+vnorm2);
582 vnorm2 = TMath::Sqrt(vnorm2);
583 Float_t pnorm2 = p[0]*p[0]+p[1]*p[1];
584 Float_t pnorm3 = TMath::Sqrt(p[2]*p[2]+pnorm2);
585 pnorm2 = TMath::Sqrt(pnorm2);
587 fPointAngleFi = (v[0]*p[0]+v[1]*p[1])/(vnorm2*pnorm2);
588 fPointAngleTh = (v[2]*p[2]+vnorm2*pnorm2)/(vnorm3*pnorm3);
589 fPointAngle = (v[0]*p[0]+v[1]*p[1]+v[2]*p[2])/(vnorm3*pnorm3);
594 void AliESDRecKinkInfo::Update()
597 if ( (fT1.fTPCOn)&& (fT2.fTPCOn)){
599 // IF BOTH RECONSTRUCTED
600 Float_t distance1,distance2;
601 Double_t xx[3],pp[3];
603 Double_t xd[3],pd[3],signd;
604 Double_t xm[3],pm[3],signm;
605 for (Int_t i=0;i<3;i++){
606 xd[i] = fT2.fTPCinR1[i];
607 pd[i] = fT2.fTPCinP1[i];
608 xm[i] = fT1.fTPCinR1[i];
609 pm[i] = fT1.fTPCinP1[i];
611 signd = fT2.fSign<0 ? -1:1;
612 signm = fT1.fSign<0 ? -1:1;
614 AliHelix dhelix1(xd,pd,signd);
615 dhelix1.GetMomentum(0,pp,0);
616 dhelix1.Evaluate(0,xx);
618 // Double_t x2[3],p2[3];
620 AliHelix mhelix(xm,pm,signm);
622 //find intersection linear
624 Double_t phase[2][2],radius[2];
625 Int_t points = dhelix1.GetRPHIintersections(mhelix, phase, radius,200);
626 Double_t delta1=10000,delta2=10000;
629 fMinR = TMath::Sqrt(radius[0]);
632 fMinR =TMath::Min(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
636 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
637 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
638 dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
641 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
642 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
643 dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
646 fMinR = TMath::Sqrt(radius[0]);
650 if (radius[0]<radius[1]){
651 fMinR = TMath::Sqrt(radius[0]);
655 fMinR = TMath::Sqrt(radius[1]);
661 distance1 = TMath::Min(delta1,delta2);
663 //find intersection parabolic
665 points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
666 delta1=10000,delta2=10000;
669 dhelix1.ParabolicDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
672 dhelix1.ParabolicDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
675 distance2 = TMath::Min(delta1,delta2);
678 dhelix1.Evaluate(phase[0][0],fXr);
679 dhelix1.GetMomentum(phase[0][0],fPdr);
680 mhelix.GetMomentum(phase[0][1],fPm);
681 dhelix1.GetAngle(phase[0][0],mhelix,phase[0][1],fAngle);
682 fRr = TMath::Sqrt(radius[0]);
685 dhelix1.Evaluate(phase[1][0],fXr);
686 dhelix1.GetMomentum(phase[1][0], fPdr);
687 mhelix.GetMomentum(phase[1][1], fPm);
688 dhelix1.GetAngle(phase[1][0],mhelix,phase[1][1],fAngle);
689 fRr = TMath::Sqrt(radius[1]);
691 fDist1 = TMath::Sqrt(distance1);
692 fDist2 = TMath::Sqrt(distance2);
695 Double_t x,alpha,param[5],cov[15];
697 fT1.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
698 fT1.GetESDtrack()->GetInnerExternalCovariance(cov);
699 AliExternalTrackParam paramm(x,alpha,param,cov);
701 fT2.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
702 fT2.GetESDtrack()->GetInnerExternalCovariance(cov);
703 AliExternalTrackParam paramd(x,alpha,param,cov);
706 kink.Update(¶mm,¶md);
708 Double_t diff = kink.fRr-fRr;
709 Double_t diff2 = kink.fDist2-fDist2;
710 printf("Diff\t%f\t%f\n",diff,diff2);