Adding the AliAnalysisGUI class which is the main class that controls the GUI.
[u/mrichter/AliRoot.git] / ITS / AliITSPid.cxx
CommitLineData
8057c644 1//------------------------------------------------------------//
2// Class for identification of pions,kaons and protons in ITS //
3// Prior particles population (probabilities) are taken from //
4// Hijing event generator. //
5//------------------------------------------------------------//
6// #include <stdlib.h>
23efe5f1 7#include "AliITSPid.h"
7d62fb64 8//#include "TMath.h"
4ae5bbc4 9#include <Riostream.h>
e75c69cd 10#include <TClonesArray.h>
11#include <TVector.h>
5aa865de 12#include "AliKalmanTrack.h"
8057c644 13#include "AliITSIOTrack.h"
e75c69cd 14#include "AliITStrackV2.h"
15#include <TF1.h>
16
23efe5f1 17ClassImp(AliITSPid)
8057c644 18
19
8221b41b 20AliITSPid::AliITSPid(const AliITSPid &source) : TObject(source),
21fMxtrs(source.fMxtrs),
22fTrs(source.fTrs),
23fWpi(source.fWpi),
24fWk(source.fWk),
25fWp(source.fWp),
26fRpik(source.fRpik),
27fRppi(source.fRppi),
28fRpka(source.fRpka),
29fRp(source.fRp),
30fPcode(source.fPcode),
31fSigmin(source.fSigmin),
32fSilent(source.fSilent),
33fCutKa(source.fCutKa),
34fCutPr(source.fCutPr),
35fggpi(source.fggpi),
36fggka(source.fggka),
37fggpr(source.fggpr){
8057c644 38 // Copy constructor. This is a function which is not allowed to be
8057c644 39
8221b41b 40
8057c644 41}
8221b41b 42
8057c644 43//______________________________________________________________________
8221b41b 44AliITSPid& AliITSPid::operator=(const AliITSPid& source){
45 // Assignment operator. This is a function which is not allowed to be
46 this->~AliITSPid();
47 new(this) AliITSPid(source);
8057c644 48 return *this;
49}
50
51//
52 Float_t AliITSPid::Qtrm(Int_t track) {
53//
54// This calculates truncated mean signal for given track.
11b3ae13 55//
23efe5f1 56 TVector q(*( this->GetVec(track) ));
38f7c306 57 Int_t ml=(Int_t)q(0);
58 if(ml<1)return 0.;
59 if(ml>6)ml=6;
60 float vf[6];
61 Int_t nl=0; for(Int_t i=0;i<ml;i++){if(q(i)>fSigmin){vf[nl]=q(i+1);nl++;}}
62 if(nl==0)return 0.;
23efe5f1 63 switch(nl){
64 case 1:q(6)=q(1); break;
65 case 2:q(6)=(q(1)+q(2))/2.; break;
66 default:
38f7c306 67 for(int fi=0;fi<2;fi++){
68 Int_t swap;
69 do{ swap=0; float qmin=vf[fi];
70 for(int j=fi+1;j<nl;j++)
71 if(qmin>vf[j]){qmin=vf[j];vf[j]=vf[fi];vf[fi]=qmin;swap=1;};
72 }while(swap==1);
73 }
23efe5f1 74 q(6)= (vf[0]+vf[1])/2.;
75 break;
76 }
77 for(Int_t i=0;i<nl;i++){q(i+1)=vf[i];} this->SetVec(track,q);
78 return (q(6));
79}
38f7c306 80
8057c644 81Float_t AliITSPid::Qtrm(Float_t qarr[6],Int_t narr) const{
82//
83//This calculates truncated mean signal for given signal set.
84//
38f7c306 85 Float_t q[6],qm,qmin;
79b921e1 86 Int_t nl,ml;
38f7c306 87 if(narr>0&&narr<7){ml=narr;}else{return 0;};
88 nl=0; for(Int_t i=0;i<ml;i++){if(qarr[i]>fSigmin){q[nl]=qarr[i];nl++;}}
79b921e1 89 if(nl==0)return 0.;
79b921e1 90 switch(nl){
91 case 1:qm=q[0]; break;
92 case 2:qm=(q[0]+q[1])/2.; break;
93 default:
38f7c306 94 Int_t swap;
95 for(int fi=0;fi<2;fi++){
96 do{ swap=0; qmin=q[fi];
97 for(int j=fi+1;j<nl;j++)
98 if(qmin>q[j]){qmin=q[j];q[j]=q[fi];q[fi]=qmin;swap=1;};
99 }while(swap==1);
100 }
101 qm= (q[0]+q[1])/2.;
79b921e1 102 break;
103 }
104 return qm;
105}
11b3ae13 106
8057c644 107Int_t AliITSPid::Wpik(Float_t pm,Float_t q){
108 //Calcutates probabilityes of pions and kaons
109 //Returns particle code for dominant probability.
11b3ae13 110 Double_t par[6];
111 for(int i=0;i<6;i++){par[i]=fGGpi[i]->Eval(pm);}
112 fggpi->SetParameters(par);
23efe5f1 113
11b3ae13 114 for(int i=0;i<3;i++){par[i]=fGGka[i]->Eval(pm);}
115 fggka->SetParameters(par);
38f7c306 116
11b3ae13 117 Float_t ppi=fggpi->Eval(q);
118 Float_t pka=fggka->Eval(q);
119 Float_t p=ppi+pka;
120 /*
121 if(!fSilent){
122 fggka->Print();
123 fggpi->Print();
124 if(p>0)cout<<" ppi,pka="<<ppi/p<<" "<<pka/p<<endl;
125 }
126 */
38f7c306 127
11b3ae13 128 if(p>0){
129 ppi=ppi/p;
130 pka=pka/p;
131 fWp=0.; fWpi=ppi; fWk=pka;
132 if( pka>ppi){return fPcode=321;}else{return fPcode=211;}
133 }else{return 0;}
23efe5f1 134}
135//-----------------------------------------------------------
8057c644 136Int_t AliITSPid::Wpikp(Float_t pm,Float_t q){
137 //
138 //Calcutates probabilityes of pions,kaons and protons.
139 //Returns particle code for dominant probability.
11b3ae13 140 Double_t par[6];
141 for(int i=0;i<6;i++){par[i]=fGGpi[i]->Eval(pm);}
142 fggpi->SetParameters(par);
23efe5f1 143
11b3ae13 144 for(int i=0;i<3;i++){par[i]=fGGka[i]->Eval(pm);}
145 fggka->SetParameters(par);
23efe5f1 146
11b3ae13 147 for(int i=0;i<3;i++){par[i]=fGGpr[i]->Eval(pm);}
148 fggpr->SetParameters(par);
38f7c306 149
11b3ae13 150 Float_t p,ppi,pka,ppr;
151 if( q>(fggpr->GetParameter(1)+fggpr->GetParameter(2)) )
152 { p=1.0; ppr=1.0; ppi=pka=0.0;
153 }else{
154 ppi=fggpi->Eval(q);
155 pka=fggka->Eval(q);
156 ppr=fggpr->Eval(q);
157 p=ppi+pka+ppr;
158 }
159 if(p>0){
160 ppi=ppi/p;
161 pka=pka/p;
162 ppr=ppr/p;
163 fWp=ppr; fWpi=ppi; fWk=pka;
164 //if(!fSilent)cout<<" ppi,pka,ppr="<<ppi<<" "<<pka<<" "<<ppr<<endl;
38f7c306 165
11b3ae13 166 if( ppi>pka&&ppi>ppr )
167 {return fPcode=211;}
168 else{ if(pka>ppr){return fPcode=321;}else{return fPcode=2212;}
169 }
38f7c306 170
11b3ae13 171 }else{return 0;}
23efe5f1 172}
173//-----------------------------------------------------------
174Int_t AliITSPid::GetPcode(TClonesArray* rps,Float_t pm)
175{
8057c644 176 //Returns particle code
088e0b8d 177 Info("GetPcode","method not implemented - Inputs TClonesArray *%x , Float_t %f",rps,pm);
23efe5f1 178 return 0;
179}
180//-----------------------------------------------------------
5aa865de 181Int_t AliITSPid::GetPcode(AliKalmanTrack *track)
79b921e1 182{
8057c644 183 //Returns particle code for given track.
79b921e1 184 Double_t xk,par[5]; track->GetExternalParameters(xk,par);
185 Float_t phi=TMath::ASin(par[2]) + track->GetAlpha();
186 if (phi<-TMath::Pi()) phi+=2*TMath::Pi();
187 if (phi>=TMath::Pi()) phi-=2*TMath::Pi();
188 Float_t lam=TMath::ATan(par[3]);
e75c69cd 189 Float_t pt1=TMath::Abs(par[4]);
190 Float_t mom=1./(pt1*TMath::Cos(lam));
bbc6cd2c 191 Float_t dedx=track->GetPIDsignal();
79b921e1 192 Int_t pcode=GetPcode(dedx/40.,mom);
def162f4 193// cout<<"TPCtrack dedx,mom,pcode="<<dedx<<","<<mom<<","<<pcode<<endl;
79b921e1 194 return pcode?pcode:211;
195 }
196//------------------------------------------------------------
79b921e1 197Int_t AliITSPid::GetPcode(AliITSIOTrack *track)
198{
8057c644 199 //Returns particle code for given track(V1).
79b921e1 200 Double_t px,py,pz;
201 px=track->GetPx();
202 py=track->GetPy();
203 pz=track->GetPz();
204 Float_t mom=TMath::Sqrt(px*px+py*py+pz*pz);
205//???????????????????
def162f4 206 // Float_t dedx=1.0;
207 Float_t dedx=track->GetdEdx();
79b921e1 208//???????????????????
209 Int_t pcode=GetPcode(dedx,mom);
def162f4 210// cout<<"ITSV1 dedx,mom,pcode="<<dedx<<","<<mom<<","<<pcode<<endl;
79b921e1 211return pcode?pcode:211;
212}
79b921e1 213//-----------------------------------------------------------
214Int_t AliITSPid::GetPcode(AliITStrackV2 *track)
215{
8057c644 216//Returns particle code for given track(V2).
79b921e1 217 if(track==0)return 0;
38f7c306 218 // track->Propagate(track->GetAlpha(),3.,0.1/65.19*1.848,0.1*1.848);
79b921e1 219 track->PropagateTo(3.,0.0028,65.19);
6c94f330 220 //track->PropagateToVertex(); Not needed. (I.B.)
79b921e1 221 Double_t xk,par[5]; track->GetExternalParameters(xk,par);
222 Float_t lam=TMath::ATan(par[3]);
e75c69cd 223 Float_t pt1=TMath::Abs(par[4]);
79b921e1 224 Float_t mom=0.;
e75c69cd 225 if( (pt1*TMath::Cos(lam))!=0. ){ mom=1./(pt1*TMath::Cos(lam)); }else{mom=0.;};
79b921e1 226 Float_t dedx=track->GetdEdx();
e75c69cd 227// cout<<"lam,pt1,mom,dedx="<<lam<<","<<pt1<<","<<mom<<","<<dedx<<endl;
79b921e1 228 Int_t pcode=GetPcode(dedx,mom);
def162f4 229// cout<<"ITS V2 dedx,mom,pcode="<<dedx<<","<<mom<<","<<pcode<<endl;
79b921e1 230return pcode?pcode:211;
231}
232//-----------------------------------------------------------
23efe5f1 233Int_t AliITSPid::GetPcode(Float_t q,Float_t pm)
234{
8057c644 235 //Returns particle code for given signal and momentum.
38f7c306 236 fWpi=fWk=fWp=0.; fPcode=0;
23efe5f1 237
11b3ae13 238 if ( pm<=0.400 )
239 { if( q<fCutKa->Eval(pm) )
e75c69cd 240 {return Pion();}
11b3ae13 241 else{ if( q<fCutPr->Eval(pm) )
e75c69cd 242 {return Kaon();}
243 else{return Proton();}
11b3ae13 244 }
245 }
246 if ( pm<=0.750 )
247 if ( q>fCutPr->Eval(pm) )
e75c69cd 248 {return Proton();} else {return Wpik(pm,q);};
249 if( pm<=1.10 ){ return Wpikp(pm,q); }
23efe5f1 250 return fPcode;
251}
252//-----------------------------------------------------------
253void AliITSPid::SetCut(Int_t n,Float_t pm,Float_t pilo,Float_t pihi,
254 Float_t klo,Float_t khi,Float_t plo,Float_t phi)
255{
8057c644 256 // The cut-table initializer method.
e75c69cd 257 fCut[n][0]=pm;
258 fCut[n][1]=pilo;
259 fCut[n][2]=pihi;
260 fCut[n][3]=klo;
261 fCut[n][4]=khi;
262 fCut[n][5]=plo;
263 fCut[n][6]=phi;
23efe5f1 264 return ;
265}
38f7c306 266//------------------------------------------------------------
e8d02863 267void AliITSPid::SetVec(Int_t ntrack,const TVector& info) const
23efe5f1 268{
8057c644 269 //Store track info in tracls table
e75c69cd 270TClonesArray& arr=*fTrs;
23efe5f1 271 new( arr[ntrack] ) TVector(info);
272}
273//-----------------------------------------------------------
8057c644 274TVector* AliITSPid::GetVec(Int_t ntrack) const
23efe5f1 275{
8057c644 276 //Get given track from track table
e75c69cd 277TClonesArray& arr=*fTrs;
23efe5f1 278 return (TVector*)arr[ntrack];
279}
280//-----------------------------------------------------------
281void AliITSPid::SetEdep(Int_t track,Float_t Edep)
282{
8057c644 283 //Set dEdx for given track
23efe5f1 284 TVector xx(0,11);
e75c69cd 285 if( ((TVector*)fTrs->At(track))->IsValid() )
286 {TVector yy( *((TVector*)fTrs->At(track)) );xx=yy; }
23efe5f1 287 Int_t j=(Int_t)xx(0); if(j>4)return;
288 xx(++j)=Edep;xx(0)=j;
e75c69cd 289 TClonesArray &arr=*fTrs;
23efe5f1 290 new(arr[track])TVector(xx);
291}
292//-----------------------------------------------------------
293void AliITSPid::SetPmom(Int_t track,Float_t Pmom)
294{
8057c644 295 //Set momentum for given track
23efe5f1 296 TVector xx(0,11);
e75c69cd 297 if( ((TVector*)fTrs->At(track))->IsValid() )
298 {TVector yy( *((TVector*)fTrs->At(track)) );xx=yy; }
23efe5f1 299 xx(10)=Pmom;
e75c69cd 300 TClonesArray &arr=*fTrs;
23efe5f1 301 new(arr[track])TVector(xx);
302}
303//-----------------------------------------------------------
304void AliITSPid::SetPcod(Int_t track,Int_t partcode)
305{
8057c644 306 //Set particle code for given track
23efe5f1 307 TVector xx(0,11);
e75c69cd 308 if( ((TVector*)fTrs->At(track))->IsValid() )
309 {TVector yy( *((TVector*)fTrs->At(track)) );xx=yy; }
23efe5f1 310 if(xx(11)==0)
e75c69cd 311 {xx(11)=partcode; fMxtrs++;
312 TClonesArray &arr=*fTrs;
23efe5f1 313 new(arr[track])TVector(xx);
314 }
315}
316//-----------------------------------------------------------
317void AliITSPid::Print(Int_t track)
8057c644 318{
319 //Prints information for given track
320cout<<fMxtrs<<" tracks in AliITSPid obj."<<endl;
e75c69cd 321 if( ((TVector*)fTrs->At(track))->IsValid() )
322 {TVector xx( *((TVector*)fTrs->At(track)) );
23efe5f1 323 xx.Print();
324 }
325 else
326 {cout<<"No data for track "<<track<<endl;return;}
327}
328//-----------------------------------------------------------
329void AliITSPid::Tab(void)
330{
8057c644 331 //Make PID for tracks stored in tracks table
e75c69cd 332if(fTrs->GetEntries()==0){cout<<"No entries in TAB"<<endl;return;}
23efe5f1 333cout<<"------------------------------------------------------------------------"<<endl;
334cout<<"Nq"<<" q1 "<<" q2 "<<" q3 "<<" q4 "<<" q5 "<<
335 " Qtrm " <<" Wpi "<<" Wk "<<" Wp "<<"Pmom "<<endl;
336cout<<"------------------------------------------------------------------------"<<endl;
e75c69cd 337for(Int_t i=0;i<fTrs->GetEntries();i++)
23efe5f1 338{
e75c69cd 339 TVector xx( *((TVector*)fTrs->At(i)) );
23efe5f1 340 if( xx.IsValid() && xx(0)>0 )
341 {
e75c69cd 342 TVector xx( *((TVector*)fTrs->At(i)) );
23efe5f1 343 if(xx(0)>=2)
344 {
345// 1)Calculate Qtrm
e75c69cd 346 xx(6)=(this->Qtrm(i));
23efe5f1 347
348 }else{
349 xx(6)=xx(1);
350 }
351// 2)Calculate Wpi,Wk,Wp
352 this->GetPcode(xx(6),xx(10)/1000.);
353 xx(7)=GetWpi();
354 xx(8)=GetWk();
355 xx(9)=GetWp();
356// 3)Print table
357 if(xx(0)>0){
def162f4 358// cout<<xx(0)<<" ";
23efe5f1 359 for(Int_t j=1;j<11;j++){
38f7c306 360 if(i<7){ cout.width(7);cout.precision(4);cout<<xx(j);}
361 if(i>7){ cout.width(7);cout.precision(5);cout<<xx(j);}
23efe5f1 362 }
363 cout<<endl;
364 }
365// 4)Update data in TVector
e75c69cd 366 TClonesArray &arr=*fTrs;
23efe5f1 367 new(arr[i])TVector(xx);
368 }
369 else
370 {/*cout<<"No data for track "<<i<<endl;*/}
371}// End loop for tracks
372}
373void AliITSPid::Reset(void)
374{
8057c644 375 //Reset tracks table
e75c69cd 376 for(Int_t i=0;i<fTrs->GetEntries();i++){
23efe5f1 377 TVector xx(0,11);
e75c69cd 378 TClonesArray &arr=*fTrs;
23efe5f1 379 new(arr[i])TVector(xx);
380 }
381}
382//-----------------------------------------------------------
8221b41b 383AliITSPid::AliITSPid(Int_t ntrack):
384fMxtrs(0),
385fTrs(0),
386fWpi(0),
387fWk(0),
388fWp(0),
389fRpik(0),
390fRppi(0),
391fRpka(0),
392fRp(0),
393fPcode(0),
394fSigmin(0.01),
395fSilent(0),
396fCutKa(0),
397fCutPr(0),
398fggpi(0),
399fggka(0),
400fggpr(0){
8057c644 401 //Constructor for AliITSPid class
8221b41b 402
e75c69cd 403 fTrs = new TClonesArray("TVector",ntrack);
404 TClonesArray &arr=*fTrs;
23efe5f1 405 for(Int_t i=0;i<ntrack;i++)new(arr[i])TVector(0,11);
8221b41b 406 //
11b3ae13 407 fCutKa=new TF1("fcutka","pol4",0.05,0.4);
408 Double_t ka[5]={25.616, -161.59, 408.97, -462.17, 192.86};
409 fCutKa->SetParameters(ka);
410 //
411 fCutPr=new TF1("fcutpr","[0]/x/x+[1]",0.05,1.1);
412 Double_t pr[2]={0.70675,0.4455};
413 fCutPr->SetParameters(pr);
414 //
415 //---------- signal fit ----------
416{//Pions
417fGGpi[0]=new TF1("fp1pi","pol4",0.34,1.2);
e75c69cd 418 Double_t parpi0[10]={ -1.9096471071e+03, 4.5354331545e+04, -1.1860738840e+05,
11b3ae13 419 1.1405329025e+05, -3.8289694496e+04 };
e75c69cd 420 fGGpi[0]->SetParameters(parpi0);
11b3ae13 421fGGpi[1]=new TF1("fp2pi","[0]/x/x+[1]",0.34,1.2);
e75c69cd 422 Double_t parpi1[10]={ 1.0791668283e-02, 9.7347716496e-01 };
423 fGGpi[1]->SetParameters(parpi1);
11b3ae13 424fGGpi[2]=new TF1("fp3pi","[0]/x/x+[1]",0.34,1.2);
e75c69cd 425 Double_t parpi2[10]={ 5.8191602279e-04, 9.7285601334e-02 };
426 fGGpi[2]->SetParameters(parpi2);
11b3ae13 427fGGpi[3]=new TF1("fp4pi","pol4",0.34,1.2);
e75c69cd 428 Double_t parpi3[10]={ 6.6267353195e+02, 7.1595101104e+02, -5.3095111914e+03,
11b3ae13 429 6.2900977606e+03, -2.2935862292e+03 };
e75c69cd 430 fGGpi[3]->SetParameters(parpi3);
11b3ae13 431fGGpi[4]=new TF1("fp5pi","[0]/x/x+[1]",0.34,1.2);
e75c69cd 432 Double_t parpi4[10]={ 9.0419011783e-03, 1.1628922525e+00 };
433 fGGpi[4]->SetParameters(parpi4);
11b3ae13 434fGGpi[5]=new TF1("fp6pi","[0]/x/x+[1]",0.34,1.2);
e75c69cd 435 Double_t parpi5[10]={ 1.8324872519e-03, 2.1503968838e-01 };
436 fGGpi[5]->SetParameters(parpi5);
11b3ae13 437}//End Pions
438{//Kaons
439fGGka[0]=new TF1("fp1ka","pol4",0.24,1.2);
e75c69cd 440 Double_t parka0[20]={
11b3ae13 441 -1.1204243395e+02,4.6716191428e+01,2.2584059281e+03,
442 -3.7123338009e+03,1.6003647641e+03 };
e75c69cd 443 fGGka[0]->SetParameters(parka0);
11b3ae13 444fGGka[1]=new TF1("fp2ka","[0]/x/x+[1]",0.24,1.2);
e75c69cd 445 Double_t parka1[20]={
11b3ae13 446 2.5181172905e-01,8.7566001814e-01 };
e75c69cd 447 fGGka[1]->SetParameters(parka1);
11b3ae13 448fGGka[2]=new TF1("fp3ka","pol6",0.24,1.2);
e75c69cd 449 Double_t parka2[20]={
11b3ae13 450 8.6236021573e+00,-7.0970427531e+01,2.4846827669e+02,
451 -4.6094401290e+02,4.7546751408e+02,-2.5807112462e+02,
452 5.7545491696e+01 };
e75c69cd 453 fGGka[2]->SetParameters(parka2);
11b3ae13 454}//End Kaons
455{//Protons
456fGGpr[0]=new TF1("fp1pr","pol4",0.4,1.2);
e75c69cd 457 Double_t parpr0[10]={
11b3ae13 458 6.0150106543e+01,-8.8176206410e+02,3.1222644604e+03,
459 -3.5269200901e+03,1.2859128345e+03 };
e75c69cd 460 fGGpr[0]->SetParameters(parpr0);
11b3ae13 461fGGpr[1]=new TF1("fp2pr","[0]/x/x+[1]",0.4,1.2);
e75c69cd 462 Double_t parpr1[10]={
11b3ae13 463 9.4970837607e-01,7.3573504201e-01 };
e75c69cd 464 fGGpr[1]->SetParameters(parpr1);
11b3ae13 465fGGpr[2]=new TF1("fp3pr","[0]/x/x+[1]",0.4,1.2);
e75c69cd 466 Double_t parpr2[10]={
11b3ae13 467 1.2498403757e-01,2.7845072306e-02 };
e75c69cd 468 fGGpr[2]->SetParameters(parpr2);
11b3ae13 469}//End Protons
470 //----------- end fit -----------
471
472 fggpr=new TF1("ggpr","gaus",0.4,1.2);
473 fggpi=new TF1("ggpi","gaus+gaus(3)",0.4,1.2);
474 fggka=new TF1("ggka","gaus",0.4,1.2);
475
476 //-------------------------------------------------
e75c69cd 477const int kInf=10;
23efe5f1 478// Ncut Pmom pilo pihi klo khi plo phi
479// cut[j] [0] [1] [2] [3] [4] [5] [6]
480//----------------------------------------------------------------
e75c69cd 481 SetCut( 1, 0.12 , 0. , 0. , kInf , kInf , kInf , kInf );
482 SetCut( 2, 0.20 , 0. , 6.0 , 6.0 , kInf , kInf , kInf );
483 SetCut( 3, 0.30 , 0. , 3.5 , 3.5 , 9.0 , 9.0 , kInf );
484 SetCut( 4, 0.41 , 0. , 1.9 , 1.9 , 4.0 , 4.0 , kInf );
23efe5f1 485//----------------------------------------------------------------
e75c69cd 486 SetCut( 5, 0.47 , 0.935, 0.139, 1.738 , 0.498 , 3.5 , kInf ); //410-470
487 SetCut( 6, 0.53 , 0.914, 0.136, 1.493 , 0.436 , 3.0 , kInf ); //470-530
23efe5f1 488//----------------------------------------------------------------
e75c69cd 489 SetCut( 7, 0.59 , 0.895, 0.131, 1.384 , 0.290 , 2.7 , kInf ); //530-590
490 SetCut( 8, 0.65 , 0.887, 0.121, 1.167 , 0.287 , 2.5 , kInf ); //590-650
491 SetCut( 9, 0.73 , 0.879, 0.120, 1.153 , 0.257 , 2.0 , kInf ); //650-730
23efe5f1 492//----------------------------------------------------------------
38f7c306 493 SetCut( 10, 0.83 , 0.880, 0.126, 1.164 , 0.204 , 2.308 , 0.297 ); //730-830
494 SetCut( 11, 0.93 , 0.918, 0.145, 1.164 , 0.204 , 2.00 , 0.168 ); //830-930
495 SetCut( 12, 1.03 , 0.899, 0.128, 1.164 , 0.204 ,1.80 , 0.168);
496 //------------------------ pi,K ---------------------
e75c69cd 497 fAprob[0][0]=1212; fAprob[1][0]=33.; // fAprob[0][i] - const for pions,cut[i+5]
498 fAprob[0][1]=1022; fAprob[1][1]=46.2 ; // fAprob[1][i] - kaons
38f7c306 499 //---------------------------------------------------
e75c69cd 500 fAprob[0][2]= 889.7; fAprob[1][2]=66.58; fAprob[2][2]=14.53;
501 fAprob[0][3]= 686.; fAprob[1][3]=88.8; fAprob[2][3]=19.27;
502 fAprob[0][4]= 697.; fAprob[1][4]=125.6; fAprob[2][4]=28.67;
38f7c306 503 //------------------------ pi,K,p -------------------
e75c69cd 504 fAprob[0][5]= 633.7; fAprob[1][5]=100.1; fAprob[2][5]=37.99; // fAprob[2][i] - protons
505 fAprob[0][6]= 469.5; fAprob[1][6]=20.74; fAprob[2][6]=25.43;
506 fAprob[0][7]= 355.; fAprob[1][7]=
507 355.*(20.74/469.5); fAprob[2][7]=34.08;
23efe5f1 508}
11b3ae13 509//End AliITSPid.cxx