]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/AliITStrack.cxx.bck
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Source file of the ITS Kalman tracking
[u/mrichter/AliRoot.git] / ITS / AliITStrack.cxx.bck
CommitLineData
dce1757b 1#include <TMath.h>
2#include <TVector.h>
3#include <TMatrix.h>
4#include <TObjArray.h>
5
6#include "AliITS.h"
7#include "AliRun.h"
8#include "AliITStrack.h"
9
10
11ClassImp(AliITStrack)
12
13AliITStrack::AliITStrack() {
14//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
15// default constructor
16
17 fvTrack.ResizeTo(5);
18 fmCovariance = new TMatrix(5,5);
19 flistCluster = new TObjArray;
20 fNumClustInTrack =0;
21 fChi2=-1;
22 flabel =0;
23 fVertex.ResizeTo(3);
24 fErrorVertex.ResizeTo(3);
25 fLayer = -1;
26 ClusterInTrack = new TMatrix(6,9);
27 for(Int_t i=0; i<6; i++) (*ClusterInTrack)(i,6)=(*ClusterInTrack)(i,7)=
28 (*ClusterInTrack)(i,8)=-1.;
29 rtrack=0.;
30 //alphaprov=-50.; //provvisorio
31 d2.ResizeTo(6);
32 tgl2.ResizeTo(6);
33 dtgl.ResizeTo(6);
34}
35
36
37
38AliITStrack::AliITStrack(const AliITStrack &cobj) {
39//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
40
41 fvTrack.ResizeTo(5);
42 fmCovariance = new TMatrix(5,5);
43 ClusterInTrack = new TMatrix(6,9);
44 for(Int_t i=0; i<6; i++) (*ClusterInTrack)(i,6)=(*ClusterInTrack)(i,7)=
45 (*ClusterInTrack)(i,8)=-1.;
46 flistCluster = new TObjArray;
47 fVertex.ResizeTo(3);
48 fErrorVertex.ResizeTo(3);
49 fVertex = cobj.fVertex;
50 fErrorVertex = cobj.fErrorVertex;
51 flabel = cobj.flabel;
52 fLayer=cobj.fLayer;
53 fTPCtrack = cobj.fTPCtrack;
54 fNumClustInTrack = cobj.fNumClustInTrack;
55 fChi2= cobj.fChi2;
56 fvTrack = cobj.fvTrack;
57 rtrack=cobj.rtrack;
58 Dv=cobj.Dv;
59 Zv=cobj.Zv;
60 sigmaDv=cobj.sigmaDv;
61 sigmaZv=cobj.sigmaZv;
62 d2.ResizeTo(6);
63 tgl2.ResizeTo(6);
64 dtgl.ResizeTo(6);
65 d2=cobj.d2;
66 tgl2=cobj.tgl2;
67 dtgl=cobj.dtgl;
68 //alphaprov=cobj.alphaprov; //provvisorio
69
70 *fmCovariance = *cobj.fmCovariance;
71
72 *ClusterInTrack = *cobj.ClusterInTrack;
73
74 for(Int_t i=0; i<cobj.flistCluster->GetSize(); i++)
75 flistCluster->AddLast(cobj.flistCluster->At(i));
76
77}
78
79AliITStrack::AliITStrack(AliTPCtrack &obj)
80{
81//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
82
83 fTPCtrack = &obj;
84 fvTrack.ResizeTo(5);
85 fVertex.ResizeTo(3);
86 fErrorVertex.ResizeTo(3);
87
88 d2.ResizeTo(6);
89 tgl2.ResizeTo(6);
90 dtgl.ResizeTo(6);
91 AliGenerator *gener = gAlice->Generator();
92 Float_t Vxg,Vyg,Vzg;
93 gener->GetOrigin(Vxg,Vyg,Vzg);
94
95
96 fVertex(0)=(Double_t)Vxg;
97 fVertex(1)=(Double_t)Vyg;
98 fVertex(2)=(Double_t)Vzg;
99
100 fLayer = 7;
101 fmCovariance = new TMatrix(5,5);
102 ClusterInTrack = new TMatrix(6,9);
103 for(Int_t i=0; i<6; i++) (*ClusterInTrack)(i,6)=(*ClusterInTrack)(i,7)=
104 (*ClusterInTrack)(i,8)=-1.;
105 flistCluster = new TObjArray;
106 fNumClustInTrack = 0;
107 LmTPC();
108
109}
110
111AliITStrack::~AliITStrack() {
112
113 //destructor
114
115 if(fmCovariance) delete fmCovariance;
116 if(flistCluster) delete flistCluster;
117 if(ClusterInTrack) delete ClusterInTrack;
118
119}
120
121
122
123void AliITStrack::LmTPC() {
124
125// Transform the TPC state vector from TPC-local to master and build a new state vector ITS-type
126// The covariance matrix is also modified accordingly
127//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
128
129
130 TVector tpctrack(5);
131 Double_t Alpha = fTPCtrack->GetAlpha();
132
133 //printf("LmTPC: Alpha %f\n",Alpha);
134
135 tpctrack(0) = fTPCtrack->GetY();
136 tpctrack(1) = fTPCtrack->GetZ();
137 tpctrack(2) = fTPCtrack->GetC();
138 tpctrack(3) = fTPCtrack->GetEta();
139 tpctrack(4) = fTPCtrack->GetTgl();
140 //cout<< " tpctracks = "<<tpctrack(0)<<" "<<tpctrack(1)<<" "<<tpctrack(2)<<" "<<tpctrack(3)<<" "<<tpctrack(4)<<
141 //"\n ";
142
143 Double_t xm, ym, zm;
144 Double_t sina = TMath::Sin(Alpha);
145 Double_t cosa = TMath::Cos(Alpha);
146 Double_t xl= fTPCtrack->GetX();
147 xm = xl * cosa - tpctrack(0)*sina;
148 ym = xl * sina + tpctrack(0)*cosa;
149 zm = tpctrack(1);
150 //cout<<" xl e alpha = "<<xl<<" "<<Alpha<<"\n"; getchar();
151
152 Double_t x0m,y0m;
153
154 ///////////////////////////////////// determine yo //////////////////////////////////////////////////
155 Double_t Vxl=fVertex(0)*cosa+fVertex(1)*sina;
156 Double_t Vyl= -fVertex(0)*sina+fVertex(1)*cosa;
157 Double_t Xo,Yo, signy;
158 Double_t R = 1./tpctrack(2);
159 Xo = tpctrack(3) / tpctrack(2);
160 xoTPC=Xo;
161 Double_t Yo1, Yo2, diffsq1, diffsq2;
162 Yo1 = tpctrack(0) + TMath::Sqrt(R*R - (xl-Xo)*(xl-Xo));
163 Yo2 = tpctrack(0) - TMath::Sqrt(R*R - (xl-Xo)*(xl-Xo));
164 diffsq1=TMath::Abs((Yo1-Vyl)*(Yo1-Vyl)+(Xo-Vxl)*(Xo-Vxl)-R*R);
165 diffsq2=TMath::Abs((Yo2-Vyl)*(Yo2- Vyl)+(Xo-Vxl)*(Xo-Vxl)-R*R);
166 if(diffsq1<diffsq2) {Yo=Yo1; signy=1.;} else {Yo=Yo2; signy=-1.;};
167 ////////////////////////////////////////////////////////////////////////////////////////////////////
168
169 x0m = Xo * cosa - Yo * sina;
170 y0m = Xo * sina + Yo * cosa;
171
172 rtrack=TMath::Sqrt(xm*xm+ym*ym);
173 Double_t pigre=TMath::Pi();
174 Double_t phi;
175 if(ym == 0. || xm == 0.) {
176 if(ym == 0. && xm == 0.) {cout << " Error in AliITStrack::LmTPC x=0 and y=0 \n"; getchar();}
177 if(ym ==0. && xm>0.) phi=0.;
178 if(ym==0. && xm<0.) phi=pigre;
179 if(xm==0 && ym>0.) phi=pigre/2.;
180 if(xm==0 && ym<0.) phi=1.5*pigre;
181 }
182 else {
183 if (xm>0. && ym>0.) phi=TMath::ATan(ym/xm);
184 if (xm<0. && ym>0.) phi=pigre+TMath::ATan(ym/xm);
185 if (xm<0. && ym<0.) phi=pigre+TMath::ATan(ym/xm);
186 if (xm>0. && ym<0.) phi=(2*pigre)+TMath::ATan(ym/xm);
187 };
188 if(phi<0. || phi>(2*pigre)) {cout<<"attention error on phi in AliITStrack:LmTPC \n"; getchar();}
189
190 fvTrack(0)= phi;
191 fvTrack(1)=zm;
192 fvTrack(3)=tpctrack(4);
193 fvTrack(4)=tpctrack(2);
194 /*
195 //provvisorio
196 alphaprov=pigre-PhiDef(x0m,y0m)-PhiDef(fVertex(0),fVertex(1));
197 if(alphaprov<0.) alphaprov=alphaprov+2.*pigre;
198 //cout<<" PhiDef(x0m,y0m) ="<<PhiDef(x0m,y0m)<<"\n";
199 //cout<<" PhiDef(fVertex(0),fVertex(1)) = "<<PhiDef(fVertex(0),fVertex(1))<<"\n";
200 //cout<<"alphaprov = "<<alphaprov<<"\n"; getchar();
201 //cout<<" fvTrack(0) prima = "<<fvTrack(0)<<"\n";
202 fvTrack(0)-=alphaprov;
203 if(fvTrack(0)<0.) fvTrack(0)+= 2.*pigre;
204 //cout<<" fvTrack(0) dopo = "<<fvTrack(0)<<" alphaprov = "<<alphaprov<<"\n";
205 cout<<" fvertex = "<<fVertex(0)<<" "<<fVertex(1)<<" "<<fVertex(2)<<"\n";
206 fVertex(0)=0.;
207 fVertex(1)=0.;
208///////////////////////////////////////////////////////////////////////////////////////////////////
209*/
210
211 Double_t dd=TMath::Sqrt((x0m-fVertex(0))*(x0m-fVertex(0))+(y0m-fVertex(1))*(y0m-fVertex(1)));
212 Double_t signdd;
213 if (R>0) signdd=1.; else signdd=-1.;
214 fvTrack(2)=signdd*dd-R;
215 //cout<<" fvertex = "<<fVertex(0)<<" "<<fVertex(1)<<" "<<fVertex(2)<<"\n";
216 //cout<<" fvTrack(2) ="<<fvTrack(2)<<"\n"; getchar();
217
218 TMatrix localM(5,5);
219 Double_t cov[15];
220 fTPCtrack->GetCovariance(cov);
221 localM(0,0)=cov[0];
222 localM(1,1)=cov[2];
223 localM(2,2)=cov[5];
224 localM(3,3)=cov[9];
225 localM(4,4)=cov[14];
226 localM(1,0)=localM(0,1)=cov[1];
227 localM(2,0)=localM(0,2)=cov[3];
228 localM(2,1)=localM(1,2)=cov[4];
229 localM(3,0)=localM(0,3)=cov[6];
230 localM(3,1)=localM(1,3)=cov[7];
231 localM(3,2)=localM(2,3)=cov[8];
232 localM(4,0)=localM(0,4)=cov[10];
233 localM(4,1)=localM(1,4)=cov[11];
234 localM(4,2)=localM(2,4)=cov[12];
235 localM(4,3)=localM(3,4)=cov[13];
236
237 TMatrix F(5,5);
238
239 Double_t dfidy, dDdy, dDdC, dDdeta;
240
241 dfidy=(xm*cosa+ym*sina)/(rtrack*rtrack);
242 dDdy=signdd*((y0m-fVertex(1))*cosa-(x0m-fVertex(0))*sina)/dd;
243 Double_t dyodr=signy*(R+(xl-Xo)*tpctrack(3))/TMath::Sqrt(R*R-(xl-Xo)*(xl-Xo));
244 Double_t dyomdr=sina*tpctrack(3)+cosa*dyodr;
245 Double_t dxomdr=cosa*tpctrack(3)-sina*dyodr;
246 Double_t ddddR=((x0m-fVertex(0))*dxomdr+(y0m-fVertex(1))*dyomdr)/dd;
247 dDdC=-R*R*(signdd*ddddR-1.);
248 Double_t dyoldxol=signy*(xl-Xo)/TMath::Sqrt(R*R-(xl-Xo)*(xl-Xo));
249 Double_t dxomdeta=R*(cosa-sina*dyoldxol);
250 Double_t dyomdeta=R*(sina+cosa*dyoldxol);
251 dDdeta=signdd*((x0m-fVertex(0))*dxomdeta+(y0m-fVertex(1))*dyomdeta)/dd;
252 F(0,0)=dfidy;
253 F(0,1)=F(0,2)=F(0,3)=F(0,4)=0.;
254 F(1,1)=1.;
255 F(1,0)=F(1,2)=F(1,3)=F(1,4)=0.;
256 F(2,0)=dDdy;
257 F(2,2)=dDdC;
258 F(2,3)=dDdeta;
259 F(2,1)=F(2,4)=0.;
260 F(3,0)=F(3,1)=F(3,2)=F(3,3)=0.;
261 F(3,4)=1.;
262 F(4,0)=F(4,1)=F(4,3)=F(4,4)=0.;
263 F(4,2)=1.;
264// cout<<" Matrice F \n";
265// F.Print(); getchar();
266 TMatrix Ft(TMatrix::kTransposed,F);
267 TMatrix temp(localM);
268 temp*=Ft;
269 TMatrix B(F);
270 B*=temp; // B=F*localM*Ft
271// cout<<" Matrice C\n";
272// B.Print(); getchar();
273 *fmCovariance = B;
274}
275
276
277AliITStrack &AliITStrack::operator=(AliITStrack obj) {
278//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
279
280 delete fmCovariance;
281 delete flistCluster;
282 delete ClusterInTrack;
283 fmCovariance = new TMatrix(5,5);
284 ClusterInTrack = new TMatrix(6,9);
285 flistCluster = new TObjArray;
286 flabel = obj.flabel;
287 fTPCtrack = obj.fTPCtrack;
288 fvTrack.ResizeTo(5);
289 fNumClustInTrack = obj.fNumClustInTrack;
290 fChi2= obj.fChi2;
291 fVertex=obj.fVertex;
292 fErrorVertex=obj.fErrorVertex;
293 fvTrack = obj.fvTrack;
294 fLayer=obj.fLayer;
295 rtrack=obj.rtrack;
296 Dv=obj.Dv;
297 Zv=obj.Zv;
298 sigmaDv=obj.sigmaDv;
299 sigmaZv=obj.sigmaZv;
300 d2=obj.d2;
301 tgl2=obj.tgl2;
302 dtgl=obj.dtgl;
303 //alphaprov=obj.alphaprov; //proviisorio
304
305 *fmCovariance = *obj.fmCovariance;
306 *ClusterInTrack = *obj.ClusterInTrack;
307 for(Int_t i=0; i<obj.flistCluster->GetSize(); i++) flistCluster->AddLast(obj.flistCluster->At(i));
308
309 return *this;
310
311}
312
313void AliITStrack::PutCluster(Int_t layerc, TVector vecclust) {
314//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
315
316 (*ClusterInTrack)(layerc,0) = vecclust(0);
317 (*ClusterInTrack)(layerc,1) = vecclust(1);
318 (*ClusterInTrack)(layerc,2) = vecclust(2);
319 (*ClusterInTrack)(layerc,3) = vecclust(3);
320 (*ClusterInTrack)(layerc,4) = vecclust(4);
321 (*ClusterInTrack)(layerc,5) = vecclust(5);
322 (*ClusterInTrack)(layerc,6) = vecclust(6);
323 (*ClusterInTrack)(layerc,7) = vecclust(7);
324 (*ClusterInTrack)(layerc,8) = vecclust(8);
325
326}
327
328
329void AliITStrack::GetClusters() {
330
331 TMatrix A(*ClusterInTrack);
332 TMatrix B(6,3);
333 for(Int_t i=0;i<6; i++){
334 B(i,0)=A(i,6); B(i,1)=A(i,7); B(i,2)=A(i,8);
335 }
336 A.Print();
337 // B.Print();
338
339}
340
341
342TVector AliITStrack::GetLabTrack(Int_t lay) {
343 TVector VecLabel(3);
344 VecLabel(0)=( (Float_t) (*ClusterInTrack)(lay,6) );
345 VecLabel(1)=( (Float_t) (*ClusterInTrack)(lay,7) );
346 VecLabel(2)=( (Float_t) (*ClusterInTrack)(lay,8) );
347 return VecLabel;
348}
349
350void AliITStrack::Search(TVector VecTotLabref, Long_t &labref, Int_t &freq){
351//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
352// define label
353
354 Int_t vecfreq[18];
355
356 for(Int_t i=0; i<18; i++) vecfreq[i]=0;
357
358 for(Int_t i=0; i<18; i++) {
359 for(Int_t j=0; j<18; j++) {
360 if(VecTotLabref(i) == 0.) VecTotLabref(i)=-3.;
361 if( (VecTotLabref(i)>=0.) && (VecTotLabref(i)==VecTotLabref(j)) ) vecfreq[i]++;
362 }
363 }
364 Int_t imax=-1000;
365 Long_t labdefault= (Long_t)1000000.;
366 freq=0;
367 for(Int_t i=0; i<18; i++) {
368 if(vecfreq[i]>freq) {freq=vecfreq[i]; imax=i;}
369 }
370 if(imax<0) labref=labdefault; else labref=(Long_t) VecTotLabref(imax);
371}
372
373
374void AliITStrack::Propagation(Double_t rk) {
375//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
376//Propagation of track
377 Double_t duepi=2.*TMath::Pi();
378 Double_t rkm1=rtrack;
379//cout<<" rk e rkm1 dentro Propagation "<<rk<<" "<<rkm1<<"\n";
380
381 //
382 Double_t Ak=argA(rk), Akm1=argA(rkm1);
383 Double_t ak=arga(rk), akm1=arga(rkm1);
384 fvTrack(0)+=TMath::ASin(Ak)-TMath::ASin(Akm1);
385 if(fvTrack(0)>duepi) fvTrack(0)-=duepi;
386 if(fvTrack(0)<0.) fvTrack(0)+=duepi;
387
388 Double_t tgl=fvTrack(3);
389 Double_t C=fvTrack(4);
390 Double_t D=fvTrack(2);
391 Double_t Cy=C/2;
392 fvTrack(1)+=tgl/Cy*(TMath::ASin(ak)-TMath::ASin(akm1));
393 rtrack=rk;
394//cout<<"fvTrack(0) fvTrack(1) e rtrack dentro Propagation = "<<fvTrack(0)<<" "<<fvTrack(1)<<" "<<rtrack<<"\n";
395//getchar();
396
397 Double_t Bk=argB(rk), Bkm1=argB(rkm1);
398 Double_t Ck=argC(rk), Ckm1=argC(rkm1);
399 TMatrix F(5,5);
400 F(0,2)=Ck/TMath::Sqrt(1.-Ak*Ak) - Ckm1/TMath::Sqrt(1.-Akm1*Akm1);
401 F(0,4)=Bk/TMath::Sqrt(1.-Ak*Ak) - Bkm1/TMath::Sqrt(1.-Akm1*Akm1);
402 F(1,2)=tgl*D*(1./rk - 1./rkm1);
403 F(1,3) = rk - rkm1;
404 F(0,0)=F(1,1)=F(2,2)=F(3,3)=F(4,4)=1.;
405
406 TMatrix Ft(TMatrix::kTransposed,F);
407 TMatrix temp(*fmCovariance);
408//cout<<" C Matrix prima propagation =\n";
409//temp.Print(); getchar();
410 temp*=Ft;
411 TMatrix B(F);
412 B*=temp;
413//cout<<" C Matrix dopo propagation =\n"; // B=F*C*Ft
414//B.Print(); getchar();
415 *fmCovariance = B;
416
417}
418void AliITStrack::AddEL(Double_t signdE, Bool_t flagtot, Double_t mass=0.1396) {
419//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
420// add energy loss
421
422 TVector s(6);
423 s(0)=0.0026+0.00283; s(1)=0.018; s(2)=0.0094; s(3)=0.0095; s(4)=0.0091; s(5)=0.0087;
424//0.00277 is added in the first layer to take into account the energy loss in the beam pipe
425
426 Double_t cl=1.+fvTrack(3)*fvTrack(3); // cl=1/(cosl)**2 = 1 + (tgl)**2
427 Double_t sqcl=TMath::Sqrt(cl);
428 Double_t pt=GetPt();
429
430 Double_t p2=pt*pt*cl;
431 Double_t E=TMath::Sqrt(p2+mass*mass);
432 Double_t beta2=p2/(p2+mass*mass);
433
434 Double_t dE;
435 if(flagtot) {
436 Double_t stot=s(0)+s(1)+s(2)+s(3)+s(4)+s(5);
437 dE=0.153/beta2*(log(5940*beta2/(1-beta2)) - beta2)*stot*21.82*sqcl;
438 } else {
439 dE=0.153/beta2*(log(5940*beta2/(1-beta2)) - beta2)*s(fLayer-1)*21.82*sqcl;
440 }
441 dE=signdE*dE/1000.;
442
443 E+=dE;
444 Double_t p=TMath::Sqrt(E*E-mass*mass);
445 Double_t sign=1.;
446 if((fvTrack)(4) < 0.) sign=-1.;
447 pt=sign*p/sqcl;
448 Double_t CC=(0.3*0.2)/(pt*100.);
449 fvTrack(4)=CC;
450
451}
452
453void AliITStrack::Correct(Double_t rk) {
454//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
455// correct track to take into account real geometry detector
456
457 Double_t duepi=2.*TMath::Pi();
458 Double_t rkm1=rtrack;
459 Double_t Ak=argA(rk), Akm1=argA(rkm1);
460 Double_t ak=arga(rk), akm1=arga(rkm1);
461
462 fvTrack(0)+=TMath::ASin(Ak)-TMath::ASin(Akm1);
463 if(fvTrack(0)>duepi) fvTrack(0)-=duepi;
464 if(fvTrack(0)<0.) fvTrack(0)+=duepi;
465
466 Double_t tgl=fvTrack(3);
467 Double_t C=fvTrack(4);
468 Double_t Cy=C/2;
469 fvTrack(1)+=tgl/Cy*(TMath::ASin(ak)-TMath::ASin(akm1));
470 rtrack=rk;
471
472}
473
474void AliITStrack::AddMS() {
475
476////////// Modification of the covariance matrix to take into account multiple scattering ///////////
477//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
478
479 TVector s(6);
480
481 s(0)=0.0026+0.00283; s(1)=0.018; s(2)=0.0094; s(3)=0.0095; s(4)=0.0091; s(5)=0.0087;
482//0.00277 is added in the first layer to take into account the energy loss in the beam pipe
483 Double_t mass=0.1396;
484 Int_t layer=(Int_t)GetLayer();
485
486 Double_t tgl=fvTrack(3);
487 Double_t cosl=TMath::Cos(TMath::ATan(tgl));
488 Double_t D=fvTrack(2);
489 Double_t C=fvTrack(4);
490 Double_t Cy=C/2.;
491 Double_t Q20=1./(cosl*cosl);
492 Double_t Q30=C*tgl;
493
494 Double_t Q40=Cy*(rtrack*rtrack-D*D)/(1.+ 2.*Cy*D);
495 Double_t dd=D+Cy*D*D-Cy*rtrack*rtrack;
496 Double_t Q41=-1./cosl*TMath::Sqrt(rtrack*rtrack - dd*dd)/(1.+ 2.*Cy*D);
497
498 /*
499 Double_t xk=rtrack*TMath::Cos(fvTrack(0));
500 Double_t yk=rtrack*TMath::Sin(fvTrack(0));
501 Double_t rvertex=TMath::Sqrt((xk-fVertex(0))*(xk-fVertex(0)) + (yk-fVertex(1))*(yk-fVertex(1)));
502 Double_t Q40=Cy*(rvertex*rvertex-D*D)/(1.+ 2.*Cy*D);
503 Double_t dd=D+Cy*D*D-Cy*rvertex*rvertex;
504 Double_t Q41=-1./cosl*TMath::Sqrt(rvertex*rvertex - dd*dd)/(1.+ 2.*Cy*D);
505 */
506
507 TMatrix J(5,2);
508 J(0,0)=0.; J(0,1)=0.;
509 J(1,0)=0.; J(1,1)=0.;
510 J(2,0)=Q40;
511 J(2,1)=Q41;
512 J(3,0)=Q20; J(3,1)=0.;
513 J(4,0)=Q30; J(4,1)=0.;
514
515 Double_t p2=(GetPt()*GetPt())/(cosl*cosl);
516 Double_t beta2=p2/(p2+mass*mass);
517 Double_t theta2=14.1*14.1/(beta2*p2*1.e6)*(s(layer-1)/cosl);
518
519 TMatrix Jt(TMatrix::kTransposed,J);
520 TMatrix Q(J,TMatrix::kMult,Jt);
521 Q*=theta2;
522
523 (*fmCovariance)+=Q;
524
525}
526
527void AliITStrack::PrimaryTrack() {
528//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
529// calculation of part of covariance matrix for vertex constraint
530
531 Double_t Rlayer[6];
532
533 Rlayer[0]=4.; Rlayer[1]=7.; Rlayer[2]=14.9; Rlayer[3]=23.8;
534 Rlayer[4]=39.1; Rlayer[5]=43.6;
535
536 Double_t Cy=fvTrack(4)/2.;
537 Double_t tgl=(fvTrack(1)-Zv)*Cy/TMath::ASin(Cy*rtrack);
538 Double_t rtrack=1.;
539 fvTrack(0)=0.;
540 fvTrack(1)=rtrack*tgl;
541 fvTrack(2)=Dv;
542 fvTrack(3)=tgl;
543
544 TMatrix newC(5,5);
545 newC(4,4)=(*fmCovariance)(4,4);
546 (*fmCovariance)=newC;
547 AddEL(1.,1);
548 fLayer=0;
549 for (Int_t i=0; i<6; i++) {
550 Propagation(Rlayer[i]);
551 fLayer++;
552 d2(i)=(*fmCovariance)(2,2);
553 tgl2(i)=(*fmCovariance)(3,3);
554 dtgl(i)=(*fmCovariance)(2,3);
555 AddMS();
556 AddEL(-1,0);
557 }
558}
559
560
561Int_t AliITStrack::DoNotCross(Double_t rk) const{
562 Double_t C=fvTrack(4);
563 Double_t D=fvTrack(2);
564 Double_t Cy=C/2.;
565 return (TMath::Abs((Cy*rk+(1.+Cy*D)*D/rk)/(1.+2.*Cy*D))>=1.)?1:0;
566}
567
568
569Double_t AliITStrack::argA(Double_t rk) const {
570 Double_t C=fvTrack(4);
571 Double_t D=fvTrack(2);
572 Double_t Cy=C/2.;
573 Double_t arg=(Cy*rk + (1 + Cy*D)*D/rk)/(1.+ 2.*Cy*D);
574 if (TMath::Abs(arg) < 1.) return arg;
575 //cout<<"class AliITSTrack: argA out of range !\n";/* getchar();*/
576 return (arg>0) ? 1. : -1.;
577}
578
579Double_t AliITStrack::arga(Double_t rk) const {
580 Double_t C=fvTrack(4);
581 Double_t D=fvTrack(2);
582 Double_t Cy=C/2.;
583 Double_t arg=(rk*rk - D*D)/(1.+ 2.*Cy*D);
584 if (arg<0.) {/*cout<<"class AliITSTrack: arga out of range !\n";*/ arg=0.;}
585 return Cy*TMath::Sqrt(arg);
586}
587
588Double_t AliITStrack::argB(Double_t rk) const {
589 Double_t C=fvTrack(4);
590 Double_t D=fvTrack(2);
591 Double_t Cy=C/2.;
592 return (rk*rk - D*D)/(rk*(1.+ 2.*Cy*D)*(1.+ 2.*Cy*D));
593}
594
595Double_t AliITStrack::argC(Double_t rk) const {
596 Double_t C=fvTrack(4);
597 Double_t D=fvTrack(2);
598 Double_t Cy=C/2.;
599 return (1./rk - 2.*Cy*argA(rk)/(1.+ 2.*Cy*D));
600}
601/*
602Double_t AliITStrack::PhiDef(Double_t x, Double_t y){
603 Double_t pigre= TMath::Pi();
604 Double_t phi;
605 if(y == 0. || x == 0.) {
606 if(y == 0. && x == 0.) {
607 cout << " Error in AliITStracking::PhiDef x=0 and y=0 \n"; getchar();
608 }
609 if(y==0. && x>0.) phi=0.;
610 if(y==0. && x<0.) phi=pigre;
611 if(x==0 && y>0.) phi=pigre/2.;
612 if(x==0 && y<0.) phi=1.5*pigre;
613 }
614 else {
615 if (x>0. && y>0.) phi=TMath::ATan(y/x);
616 if (x<0. && y>0.) phi=pigre+TMath::ATan(y/x);
617 if (x<0. && y<0.) phi=pigre+TMath::ATan(y/x);
618 if (x>0. && y<0.) phi=(2.*pigre)+TMath::ATan(y/x);
619 }
620 if(phi<0. || phi>(2*pigre)) {
621 cout<<" Error on phi in AliITStracking::PhiDef \n"; getchar();
622 }
623 return phi;
624}
625*/