+ if (ptg<1e-33) continue; // for those not crossing 0 pad row
+
+ if (ptg<ptcutl) continue;
+ if (ptg>ptcuth) continue;
+
+ allselected++;
+
+ hgood->Fill(ptg);
+
+ AliESDtrack *track=0;
+ for (i=0; i<nentr; i++) {
+ track=event->GetTrack(i);
+ tlab=track->GetTPCLabel();
+ if (lab==TMath::Abs(tlab)) break;
+ }
+ if (i==nentr) {
+ track_notfound[itrack_notfound++]=lab;
+ continue;
+ }
+
+ //MI change - addition
+ Int_t micount=0;
+ Int_t mi;
+ AliESDtrack * mitrack;
+ for (mi=0; mi<nentr; mi++) {
+ mitrack=event->GetTrack(mi);
+ if (lab==TMath::Abs(mitrack->GetTPCLabel())) micount++;
+ }
+ if (micount>1) {
+ track_multifound[itrack_multifound]=lab;
+ track_multifound_n[itrack_multifound]=micount;
+ itrack_multifound++;
+ }
+ if ((track->GetStatus()&AliESDtrack::kTPCrefit)==0) continue;
+ if (lab==tlab) hfound->Fill(ptg);
+ else {
+ track_fake[itrack_fake++]=lab;
+ hfake->Fill(ptg);
+ }
+
+ Double_t pxpypz[3]; track->GetInnerPxPyPz(pxpypz);
+ Float_t phi=TMath::ATan2(pxpypz[1],pxpypz[0]);
+ if (phi<-TMath::Pi()) phi+=2*TMath::Pi();
+ if (phi>=TMath::Pi()) phi-=2*TMath::Pi();
+ Double_t pt=TMath::Sqrt(pxpypz[0]*pxpypz[0]+pxpypz[1]*pxpypz[1]);
+ Float_t lam=TMath::ATan2(pxpypz[2],pt);
+ Float_t pt_1=1/pt;
+
+ Int_t pdg=(Int_t)ref->GetLength(); //this is particle's PDG !
+
+ if (TMath::Abs(pdg)==11 && ptg>4.) {//high momentum electrons
+ hmpt->Fill((pt_1 - 1/ptg)/(1/ptg)*100.);
+ } else {
+ Float_t phig=TMath::ATan2(ref->Py(),ref->Px());
+ hp->Fill((phi - phig)*1000.);
+
+ Float_t lamg=TMath::ATan2(ref->Pz(),ptg);
+ hl->Fill((lam - lamg)*1000.);
+
+ hpt->Fill((pt_1 - 1/ptg)/(1/ptg)*100.);
+ }
+
+ Float_t mom=pt/TMath::Cos(lam);
+ Float_t dedx=track->GetTPCsignal();
+ hep->Fill(mom,dedx,1.);
+ if (TMath::Abs(pdg)==211) //pions
+ if (mom>0.4 && mom<0.5) {
+ he->Fill(dedx,1.);
+ }