#include <stdio.h>
#include <string.h>
#include "Rtypes.h"
#include "AliESDtrack.h"
#include "AliTPCParam.h"
#include "AliTrackReference.h"
#include "AliTPCParamSR.h"
#include "AliESD.h"
#include "AliESDfriend.h"
#include "AliESDtrack.h"
#include "AliTPCseed.h"
#include "AliITStrackMI.h"
#include "AliHelix.h"
#include "AliESDVertex.h"
#include "AliExternalTrackParam.h"
#include "AliESDkink.h"
#include "AliESDv0.h"
#include "AliV0.h"
#include "AliTreeDraw.h"
#include "AliMCInfo.h"
#include "AliGenKinkInfo.h"
#include "AliGenV0Info.h"
#include "AliESDRecKinkInfo.h"
ClassImp(AliESDRecKinkInfo)
AliESDRecKinkInfo::AliESDRecKinkInfo():
fT1(),
fT2(),
fKink(),
fDist1(0),
fDist2(0),
fInvMass(0),
fRr(0),
fMinR(0),
fDistMinR(0),
fPointAngleFi(0),
fPointAngleTh(0),
fPointAngle(0),
fStatus(0),
fRecStatus(0),
fMultiple(0),
fKinkMultiple(0)
{
for (Int_t i=0; i<3; i++) {
fPdr[i] = 0;
fXr[i] = 0;
fPm[i] = 0;
fAngle[i] = 0;
}
for (Int_t i=0; i<2; i++) { fLab[i] = 0; }
}
void AliESDRecKinkInfo::Update()
{
if ( (fT1.fTPCOn)&& (fT2.fTPCOn)){
Float_t distance1,distance2;
Double_t xx[3],pp[3];
Double_t xd[3],pd[3],signd;
Double_t xm[3],pm[3],signm;
for (Int_t i=0;i<3;i++){
xd[i] = fT2.fTPCinR1[i];
pd[i] = fT2.fTPCinP1[i];
xm[i] = fT1.fTPCinR1[i];
pm[i] = fT1.fTPCinP1[i];
}
signd = fT2.fSign<0 ? -1:1;
signm = fT1.fSign<0 ? -1:1;
AliHelix dhelix1(xd,pd,signd);
dhelix1.GetMomentum(0,pp,0);
dhelix1.Evaluate(0,xx);
AliHelix mhelix(xm,pm,signm);
Double_t phase[2][2] = { {0,0}, {0,0}};
Double_t radius[2] = {0} ;
Int_t points = dhelix1.GetRPHIintersections(mhelix, phase, radius,200);
Double_t delta1=10000,delta2=10000;
if (points==1){
fMinR = TMath::Sqrt(radius[0]);
}
if (points==2){
fMinR =TMath::Min(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
}
if (points>0){
dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
}
if (points==2){
dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
}
if (points==1){
fMinR = TMath::Sqrt(radius[0]);
fDistMinR = delta1;
}
if (points==2){
if (radius[0]<radius[1]){
fMinR = TMath::Sqrt(radius[0]);
fDistMinR = delta1;
}
else{
fMinR = TMath::Sqrt(radius[1]);
fDistMinR = delta2;
}
}
distance1 = TMath::Min(delta1,delta2);
points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
delta1=10000,delta2=10000;
if (points>0){
dhelix1.ParabolicDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
}
if (points==2){
dhelix1.ParabolicDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
}
distance2 = TMath::Min(delta1,delta2);
if (delta1<delta2){
dhelix1.Evaluate(phase[0][0],fXr);
dhelix1.GetMomentum(phase[0][0],fPdr);
mhelix.GetMomentum(phase[0][1],fPm);
dhelix1.GetAngle(phase[0][0],mhelix,phase[0][1],fAngle);
fRr = TMath::Sqrt(radius[0]);
}
else{
dhelix1.Evaluate(phase[1][0],fXr);
dhelix1.GetMomentum(phase[1][0], fPdr);
mhelix.GetMomentum(phase[1][1], fPm);
dhelix1.GetAngle(phase[1][0],mhelix,phase[1][1],fAngle);
fRr = TMath::Sqrt(radius[1]);
}
fDist1 = TMath::Sqrt(distance1);
fDist2 = TMath::Sqrt(distance2);
if (fDist2<10.5){
Double_t x,alpha,param[5],cov[15];
fT1.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
fT1.GetESDtrack()->GetInnerExternalCovariance(cov);
AliExternalTrackParam paramm(x,alpha,param,cov);
fT2.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
fT2.GetESDtrack()->GetInnerExternalCovariance(cov);
AliExternalTrackParam paramd(x,alpha,param,cov);
}
}
}
AliESDRecKinkInfo.cxx:100 AliESDRecKinkInfo.cxx:101 AliESDRecKinkInfo.cxx:102 AliESDRecKinkInfo.cxx:103 AliESDRecKinkInfo.cxx:104 AliESDRecKinkInfo.cxx:105 AliESDRecKinkInfo.cxx:106 AliESDRecKinkInfo.cxx:107 AliESDRecKinkInfo.cxx:108 AliESDRecKinkInfo.cxx:109 AliESDRecKinkInfo.cxx:110 AliESDRecKinkInfo.cxx:111 AliESDRecKinkInfo.cxx:112 AliESDRecKinkInfo.cxx:113 AliESDRecKinkInfo.cxx:114 AliESDRecKinkInfo.cxx:115 AliESDRecKinkInfo.cxx:116 AliESDRecKinkInfo.cxx:117 AliESDRecKinkInfo.cxx:118 AliESDRecKinkInfo.cxx:119 AliESDRecKinkInfo.cxx:120 AliESDRecKinkInfo.cxx:121 AliESDRecKinkInfo.cxx:122 AliESDRecKinkInfo.cxx:123 AliESDRecKinkInfo.cxx:124 AliESDRecKinkInfo.cxx:125 AliESDRecKinkInfo.cxx:126 AliESDRecKinkInfo.cxx:127 AliESDRecKinkInfo.cxx:128 AliESDRecKinkInfo.cxx:129 AliESDRecKinkInfo.cxx:130 AliESDRecKinkInfo.cxx:131 AliESDRecKinkInfo.cxx:132 AliESDRecKinkInfo.cxx:133 AliESDRecKinkInfo.cxx:134 AliESDRecKinkInfo.cxx:135 AliESDRecKinkInfo.cxx:136 AliESDRecKinkInfo.cxx:137 AliESDRecKinkInfo.cxx:138 AliESDRecKinkInfo.cxx:139 AliESDRecKinkInfo.cxx:140 AliESDRecKinkInfo.cxx:141 AliESDRecKinkInfo.cxx:142 AliESDRecKinkInfo.cxx:143 AliESDRecKinkInfo.cxx:144 AliESDRecKinkInfo.cxx:145 AliESDRecKinkInfo.cxx:146 AliESDRecKinkInfo.cxx:147 AliESDRecKinkInfo.cxx:148 AliESDRecKinkInfo.cxx:149 AliESDRecKinkInfo.cxx:150 AliESDRecKinkInfo.cxx:151 AliESDRecKinkInfo.cxx:152 AliESDRecKinkInfo.cxx:153 AliESDRecKinkInfo.cxx:154 AliESDRecKinkInfo.cxx:155 AliESDRecKinkInfo.cxx:156 AliESDRecKinkInfo.cxx:157 AliESDRecKinkInfo.cxx:158 AliESDRecKinkInfo.cxx:159 AliESDRecKinkInfo.cxx:160 AliESDRecKinkInfo.cxx:161 AliESDRecKinkInfo.cxx:162 AliESDRecKinkInfo.cxx:163 AliESDRecKinkInfo.cxx:164 AliESDRecKinkInfo.cxx:165 AliESDRecKinkInfo.cxx:166 AliESDRecKinkInfo.cxx:167 AliESDRecKinkInfo.cxx:168 AliESDRecKinkInfo.cxx:169 AliESDRecKinkInfo.cxx:170 AliESDRecKinkInfo.cxx:171 AliESDRecKinkInfo.cxx:172 AliESDRecKinkInfo.cxx:173 AliESDRecKinkInfo.cxx:174 AliESDRecKinkInfo.cxx:175 AliESDRecKinkInfo.cxx:176 AliESDRecKinkInfo.cxx:177 AliESDRecKinkInfo.cxx:178 AliESDRecKinkInfo.cxx:179 AliESDRecKinkInfo.cxx:180 AliESDRecKinkInfo.cxx:181 AliESDRecKinkInfo.cxx:182 AliESDRecKinkInfo.cxx:183 AliESDRecKinkInfo.cxx:184 AliESDRecKinkInfo.cxx:185 AliESDRecKinkInfo.cxx:186 AliESDRecKinkInfo.cxx:187 AliESDRecKinkInfo.cxx:188 AliESDRecKinkInfo.cxx:189 AliESDRecKinkInfo.cxx:190 AliESDRecKinkInfo.cxx:191 AliESDRecKinkInfo.cxx:192 AliESDRecKinkInfo.cxx:193 AliESDRecKinkInfo.cxx:194 AliESDRecKinkInfo.cxx:195 AliESDRecKinkInfo.cxx:196 AliESDRecKinkInfo.cxx:197 AliESDRecKinkInfo.cxx:198 AliESDRecKinkInfo.cxx:199 AliESDRecKinkInfo.cxx:200 AliESDRecKinkInfo.cxx:201 AliESDRecKinkInfo.cxx:202 AliESDRecKinkInfo.cxx:203 AliESDRecKinkInfo.cxx:204 AliESDRecKinkInfo.cxx:205 AliESDRecKinkInfo.cxx:206 AliESDRecKinkInfo.cxx:207 AliESDRecKinkInfo.cxx:208 AliESDRecKinkInfo.cxx:209 AliESDRecKinkInfo.cxx:210 AliESDRecKinkInfo.cxx:211 AliESDRecKinkInfo.cxx:212 AliESDRecKinkInfo.cxx:213 AliESDRecKinkInfo.cxx:214 AliESDRecKinkInfo.cxx:215 AliESDRecKinkInfo.cxx:216 AliESDRecKinkInfo.cxx:217 AliESDRecKinkInfo.cxx:218 AliESDRecKinkInfo.cxx:219 AliESDRecKinkInfo.cxx:220 AliESDRecKinkInfo.cxx:221 AliESDRecKinkInfo.cxx:222 AliESDRecKinkInfo.cxx:223