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