candidate->SetIndexT(it);
candidate->SetIndexC(ic);
candidate->SetDistance(dist);
- candidate->SetPt(track->GetPt());
+ candidate->SetPt(track->GetSignedPt());
fMatches->Add(candidate);
count++;
}
// This function compares tracks according to the their curvature
//-----------------------------------------------------------------
AliHLTITStrack *t=(AliHLTITStrack*)o;
- Double_t co=TMath::Abs(t->Get1Pt());
- Double_t c =TMath::Abs(Get1Pt());
+ Double_t co=t->OneOverPt();
+ Double_t c =OneOverPt();
// Double_t co=t->GetSigmaY2()*t->GetSigmaZ2();
// Double_t c =GetSigmaY2()*GetSigmaZ2();
if (c>co) return 1;
/**/
// inter-tracks competition ???
if ((TMath::Abs(track->GetC() - p4) / TMath::Abs(p4) < 0.2) ||
- (TMath::Abs(track->GetPt()) > 0.8)) {
+ (track->Pt() > 0.8)) {
fHBackfit->Fill(4);
// Sign only gold tracks
if (track->GetChi2() / track->GetNumberOfClusters() < 4) {
if ((seed->GetKinkIndex(0) == 0) &&
- (TMath::Abs(track->GetPt()) < 1.5)) {
+ (track->Pt() < 1.5)) {
UseClusters(track);
}
}
tany *= -1.0;
}
- Double_t sy2 = ExpectedSigmaY2(x0,track->GetTgl(),track->GetPt());
+ Double_t sy2 = ExpectedSigmaY2(x0,track->GetTgl(),track->GetSignedPt());
Double_t sz2 = ExpectedSigmaZ2(x0,track->GetTgl());
Double_t road = 15.0 * TMath::Sqrt(track->GetSigmaY2() + sy2);
if (road > 6.0) {
fHFindCl[iCheckPoint]->Fill(checkPoint[iCheckPoint]);
if (checkPoint[3]) {
- if (track->GetPt() > 0) fHMinYPos->Fill(minY);
+ if (track->GetSignedPt() > 0) fHMinYPos->Fill(minY);
else fHMinYNeg->Fill(minY);
fHMinD->Fill(minD[0], minD[1]);
//------------------------------------------------------------------
//d=x0=0.;
- Double_t p2=(1.+ GetTgl()*GetTgl())/(Get1Pt()*Get1Pt());
+ Double_t p2=(1.+ GetTgl()*GetTgl())/(GetSigned1Pt()*GetSigned1Pt());
Double_t et = p2 + GetMass()*GetMass();
Double_t beta2=p2/et;
et = sqrt(et);
//Double_t theta2=1.0259e-6*14*14/28/(beta2*p2)*TMath::Abs(d)*9.36*2.33;
Cov(5) += theta2*(1.- GetSnp()*GetSnp())*(1. + GetTgl()*GetTgl());
Cov(9) += theta2*(1. + GetTgl()*GetTgl())*(1. + GetTgl()*GetTgl());
- Cov(13)+= theta2*GetTgl()*Get1Pt()*(1. + GetTgl()*GetTgl());
- Cov(14)+= theta2*GetTgl()*Get1Pt()*GetTgl()*Get1Pt();
+ Cov(13)+= theta2*GetTgl()*GetSigned1Pt()*(1. + GetTgl()*GetTgl());
+ Cov(14)+= theta2*GetTgl()*GetSigned1Pt()*GetTgl()*GetSigned1Pt();
}
//Energy losses************************
//}
Par(4)*=(1.- et/p2*dE);
- Double_t delta44 = (dE*Get1Pt()*et/p2);
+ Double_t delta44 = (dE*GetSigned1Pt()*et/p2);
delta44*=delta44;
Cov(14)+= delta44/400.;
}
// This function compares tracks according to the their curvature
//-----------------------------------------------------------------
AliITStrackV2 *t=(AliITStrackV2*)o;
- //Double_t co=TMath::Abs(t->Get1Pt());
- //Double_t c =TMath::Abs(Get1Pt());
+ //Double_t co=OneOverPt();
+ //Double_t c =OneOverPt();
Double_t co=t->GetSigmaY2()*t->GetSigmaZ2();
Double_t c =GetSigmaY2()*GetSigmaZ2();
if (c>co) return 1;
Double_t dy = Par(0) - yv, dz = Par(1) - zv;
Double_t r2=GetX()*GetX() + dy*dy;
- Double_t p2=(1.+ GetTgl()*GetTgl())/(Get1Pt()*Get1Pt());
+ Double_t p2=(1.+ GetTgl()*GetTgl())/(GetSigned1Pt()*GetSigned1Pt());
Double_t beta2=p2/(p2 + GetMass()*GetMass());
x0*=TMath::Sqrt((1.+ GetTgl()*GetTgl())/(1.- GetSnp()*GetSnp()));
Double_t theta2=14.1*14.1/(beta2*p2*1e6)*x0;
delete t;
continue;
}
- if (TMath::Abs(1/t->Get1Pt())<AliITSReconstructor::GetRecoParam()->GetMinPtForProlongation()) {
+ if (t->Pt()<AliITSReconstructor::GetRecoParam()->GetMinPtForProlongation()) {
delete t;
continue;
}
vec(1,0)=track1->GetZ() - track3.GetZ();
vec(2,0)=track1->GetSnp() - track3.GetSnp();
vec(3,0)=track1->GetTgl() - track3.GetTgl();
- vec(4,0)=track1->Get1Pt() - track3.Get1Pt();
+ vec(4,0)=track1->GetSigned1Pt() - track3.GetSigned1Pt();
//
TMatrixD cov(5,5);
cov(0,0) = track1->GetSigmaY2()+track3.GetSigmaY2();
npoints++;
}
if (npoints>1) return
- TMath::Max(TMath::Abs(0.3*forwardtrack->Get1Pt())-0.5,0.)+
+ TMath::Max(0.3*forwardtrack->OneOverPt()-0.5,0.)+
//2*forwardtrack->fNUsed+
res/TMath::Max(double(npoints-forwardtrack->GetNSkipped()),
1./(1.+forwardtrack->GetNSkipped()));
//
AliITStrackMI * track10=(AliITStrackMI*) arr1->UncheckedAt(0);
AliITStrackMI * track20=(AliITStrackMI*) arr2->UncheckedAt(0);
- if (TMath::Abs(1./track10->Get1Pt())>0.5+TMath::Abs(1/track20->Get1Pt())) return track10;
+ if (track10->Pt()>0.5+track20->Pt()) return track10;
for (Int_t itrack=0;itrack<entries1;itrack++){
AliITStrackMI * track=(AliITStrackMI*) arr1->UncheckedAt(itrack);
//
w1 = (d2/(d1+d2)+ 2*s2/(s1+s2)+
+s2/(s1+s2)*0.5*(chi22+2.)/(chi21+chi22+4.)
- +1.*TMath::Abs(1./track10->Get1Pt())/(TMath::Abs(1./track10->Get1Pt())+TMath::Abs(1./track20->Get1Pt()))
+ +1.*track10->Pt()/(track10->Pt()+track20->Pt())
);
w2 = (d1/(d1+d2)+ 2*s1/(s1+s2)+
s1/(s1+s2)*0.5*(chi21+2.)/(chi21+chi22+4.)
- +1.*TMath::Abs(1./track20->Get1Pt())/(TMath::Abs(1./track10->Get1Pt())+TMath::Abs(1./track20->Get1Pt()))
+ +1.*track20->Pt()/(track10->Pt()+track20->Pt())
);
Double_t sumw = w1+w2;
itstrack = (AliITStrackV2*)tArray.At(i);
// pt
- pt = 1./TMath::Abs(itstrack->Get1Pt());
+ pt = itstrack->Pt();
// primary vertex from other tracks in event
Bool_t goodVtx = kFALSE;
Double_t energy = TMath::Sqrt( mass*mass + fP[3]*fP[3] + fP[4]*fP[4] + fP[5]*fP[5]);
fP[6] = energy;
fP[7] = 0;
- fQ = (track.Get1Pt() >0 ) ?1 :-1;
+ fQ = (track.GetSigned1Pt() >0 ) ?1 :-1; // fQ = track->GetSign() would do the same thing
fNDF = 0;
fChi2 = 0;
fAtProductionVertex = 0;
if ( !( pTrack->GetStatus()&AliESDtrack::kITSrefit ) ) continue;
Int_t indi[12];
if( pTrack->GetITSclusters(indi) <5 ) continue;
- Int_t PDG = ( pTrack->Get1Pt() <0 ) ?321 :211;
+ Int_t PDG = ( pTrack->GetSigned1Pt() <0 ) ?321 :211;
//* take MC PDG
{
//if multifound find track with biggest pt - returning tracks - loos energy
//
if (fIndexRecTracks[fMCInfo->fLabel*4]+1){
- Double_t pt = TMath::Abs(1/fTPCTrack->Get1Pt());
+ Double_t pt = fTPCTrack-Pt();
for (Int_t i=1;i<4;i++){
if (fIndexRecTracks[fMCInfo->fLabel*4+i]>=0){
AliTPCtrack * track = (AliTPCtrack*)fTracks->UncheckedAt(fIndexRecTracks[fMCInfo->fLabel*4+i]);
p2->cd();
p2->Draw();
TH1F* his2 = new TH1F("Ptresolution","Ptresolution",40,-5.,5.);
- fTree->Draw("100.*(abs(1./fTPCTrack.Get1Pt())-MC.fTrackRef.Pt())/MC.fTrackRef.Pt()>>Ptresolution","MC.fRowsWithDigits>100&&RC.fTPCTrack.fN>50&&RC.fMultiple==1"+cteta1+cpos1+cprim);
+ fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt()>>Ptresolution","MC.fRowsWithDigits>100&&RC.fTPCTrack.fN>50&&RC.fMultiple==1"+cteta1+cpos1+cprim);
AliLabelAxes(his2, "#Delta p_{t} / p_{t} [%]", "entries");
his2->Fit("gaus");
his2->Draw();
Double_t maxRes = 10.;
TH2F* hRes2 = new TH2F("hRes2", "residuals", nBins, bins, nBinsRes, -maxRes, maxRes);
- fTree->Draw("100.*(abs(1./fTPCTrack.Get1Pt())-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
+ fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
TH1F* hMean=0;
TH1F* hRes = CreateResHisto(hRes2, &hMean);
Double_t maxRes = 10.;
TH2F* hRes2 = new TH2F("hRes2", "residuals", nBins, bins, nBinsRes, -maxRes, maxRes);
- fTree->Draw("100.*(1./fTPCTrack.Get1Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
+ fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
TH1F* hMean=0;
TH1F* hRes = CreateResHisto(hRes2, &hMean);
}
else {
Float_t f2 =1;
- f2 = 1-20*TMath::Sqrt(t->GetSigma1Pt2())/(TMath::Abs(t->Get1Pt())+0.0066);
+ f2 = 1-20*TMath::Sqrt(t->GetSigma1Pt2())/(t->OneOverPt()+0.0066);
if (t->fBConstrain) f2=1.2;
Float_t f1 =1;
- f1 = 1-20*TMath::Sqrt(GetSigma1Pt2())/(TMath::Abs(Get1Pt())+0.0066);
+ f1 = 1-20*TMath::Sqrt(GetSigma1Pt2())/(OneOverPt()+0.0066);
if (fBConstrain) f1=1.2;
// This function compares tracks according to the their curvature
//-----------------------------------------------------------------
AliTPCtrack *t=(AliTPCtrack*)o;
- //Double_t co=TMath::Abs(t->Get1Pt());
- //Double_t c =TMath::Abs(Get1Pt());
+ //Double_t co=t->OneOverPt();
+ //Double_t c = OneOverPt();
Double_t co=t->GetSigmaY2()*t->GetSigmaZ2();
Double_t c =GetSigmaY2()*GetSigmaZ2();
if (c>co) return 1;
UInt_t index=0;
Double_t maxchi2=kMaxCHI2;
const AliTPCRow &krow=fSectors[s][nr];
- Double_t pt=1./t.Get1Pt();
+ Double_t pt=t.GetSignedPt();
Double_t sy2=AliTPCcluster::SigmaY2(t.GetX(),t.GetTgl(),pt);
Double_t sz2=AliTPCcluster::SigmaZ2(t.GetX(),t.GetTgl());
Double_t road=4.*sqrt(t.GetSigmaY2() + sy2), y=t.GetY(), z=t.GetZ();
AliTPCcluster *cl=0;
Int_t index=0;
Double_t maxchi2=kMaxCHI2;
- Double_t pt=1./seed.Get1Pt();
+ Double_t pt=seed.GetSignedPt();
Double_t sy2=AliTPCcluster::SigmaY2(seed.GetX(),seed.GetTgl(),pt);
Double_t sz2=AliTPCcluster::SigmaZ2(seed.GetX(),seed.GetTgl());
Double_t road=4.*sqrt(seed.GetSigmaY2() + sy2), z=seed.GetZ();
cmptrk.eta = part->Eta();
cmptrk.r = TMath::Sqrt(part->Vx()*part->Vx()+part->Vy()*part->Vy());
- cmptrk.pt = 1/TMath::Abs(geatrack->Get1Pt());
+ cmptrk.pt = geatrack->Pt();
cmptrk.cosl = TMath::Cos(TMath::ATan(geatrack->GetTgl()));
cmptrk.p = cmptrk.pt/cmptrk.cosl;
cmptrk.dP2 = kaltrack->GetSnp()-geatrack->GetSnp();
cmptrk.dP3 = kaltrack->GetTgl()-geatrack->GetTgl();
cmptrk.dP4 = kaltrack->GetC()-geatrack->GetC();
- cmptrk.dpt = 1/kaltrack->Get1Pt()-1/geatrack->Get1Pt();
+ cmptrk.dpt = 1/kaltrack->GetSigned1Pt()-1/geatrack->GetSigned1Pt();
// get covariance matrix
// beware: lines 3 and 4 in the matrix are inverted!
// get P and Cosl from track
cosl = TMath::Cos(TMath::ATan(fTrack.GetTgl()));
- p = 1./TMath::Abs(fTrack.Get1Pt())/cosl;
+ p = fTrack.Pt()/cosl;
trkKine[0] = p;
if (track->GetNumberOfClusters() < kMinClusters) return kFALSE;
Float_t ratio = track->GetNumberOfClusters() / (track->GetNFoundable() + 1.);
if (ratio < kMinRatio) return kFALSE;
- Float_t mpt = track->Get1Pt();
+ Float_t mpt = track->GetSigned1Pt();
if (TMath::Abs(mpt) > kMax1pt) return kFALSE;
//if (TMath::Abs(track->GetZ())>240.) return kFALSE;
//if (TMath::Abs(track->GetZ())<10.) return kFALSE;
fPar1[time] = t->GetZ();
fPar2[time] = t->GetSnp();
fPar3[time] = t->GetTgl();
- fPar4[time] = t->Get1Pt();
+ fPar4[time] = t->GetSigned1Pt();
// Store the info bis of the tracklet
Int_t *rowcol = CalculateRowCol(cl);
{
// Energy losses************************
- Double_t p2 = (1.0 + GetTgl()*GetTgl()) / (Get1Pt()*Get1Pt());
+ Double_t p2 = (1.0 + GetTgl()*GetTgl()) / (GetSigned1Pt()*GetSigned1Pt());
Double_t beta2 = p2 / (p2 + GetMass()*GetMass());
if (beta2<1.0e-10 || (5940.0 * beta2/(1.0 - beta2 + 1.0e-10) - beta2) < 0.0) {
return kFALSE;
return kTRUE;
+
}
// //_____________________________________________________________________________
// fCtt-=k30*oldty+k31*oldtz;
// fCct-=k40*oldty+k41*oldtz;
// //
+
// fCcc-=k40*oldcy+k41*oldcz;
// //
/**/
// inter-tracks competition ???
if ((TMath::Abs(track->GetC() - p4) / TMath::Abs(p4) < 0.2) ||
- (TMath::Abs(track->GetPt()) > 0.8)) {
+ (track->Pt() > 0.8)) {
fHBackfit->Fill(4);
// Sign only gold tracks
if (track->GetChi2() / track->GetNumberOfClusters() < 4) {
if ((seed->GetKinkIndex(0) == 0) &&
- (TMath::Abs(track->GetPt()) < 1.5)) {
+ (track->Pt() < 1.5)) {
UseClusters(track);
}
}
tany *= -1.0;
}
- Double_t sy2 = ExpectedSigmaY2(x0,track->GetTgl(),track->GetPt());
+ Double_t sy2 = ExpectedSigmaY2(x0,track->GetTgl(),track->GetSignedPt());
Double_t sz2 = ExpectedSigmaZ2(x0,track->GetTgl());
Double_t road = 15.0 * TMath::Sqrt(track->GetSigmaY2() + sy2);
if (road > 6.0) {
fHFindCl[iCheckPoint]->Fill(checkPoint[iCheckPoint]);
if (checkPoint[3]) {
- if (track->GetPt() > 0) fHMinYPos->Fill(minY);
+ if (track->GetSignedPt() > 0) fHMinYPos->Fill(minY);
else fHMinYNeg->Fill(minY);
fHMinD->Fill(minD[0], minD[1]);
if (!(status & AliESDtrack::kTRDout)) continue;
if (!(status & AliESDtrack::kTRDrefit)) continue;
- if (TMath::Abs(esdTrack->GetOuterParam()->GetPt()) < 1.0) continue;
+ if (esdTrack->GetOuterParam()->Pt() < 1.0) continue;
int ch=0;
while(param->GetX() > fGeo->GetTime0(ch)+2) ch++;
AliExternalTrackParam out(*param);
out.PropagateTo(aRef->LocalX(),bz);
- double dp = aRef->Pt() + out.GetPt();
+ double dp = aRef->Pt() + out.GetSignedPt();
double dy = 10. * (aRef->LocalY() - out.GetY()); // in mm
fDeltaPt->Fill(100. * dp / aRef->Pt());
comp.fTree->SetAlias("dt","(abs(MC.fPdg)==2212)*dtp+(abs(MC.fPdg)==211)*dtpi+(abs(MC.fPdg)==321)*dtk"); //delta time
comp.fTree->SetAlias("Beta","sqrt(TR.P()**2/(TR.P()**2+MC.fMass**2))");
comp.fTree->SetAlias("dtsigma","10+(1-Beta)*280");
- comp.fTree->SetAlias("dptrel","(abs(TRD0.GetPt())-MC.fTOFReferences[0].Pt())/MC.fTOFReferences[0].Pt()");
+ comp.fTree->SetAlias("dptrel","(abs(TRD0.GetSignedPt())-MC.fTOFReferences[0].Pt())/MC.fTOFReferences[0].Pt()");
//
comp.fTree->SetAlias("dphi","atan2(TR.fY,TR..fX)-atan2(TR.fPy,TR..fPx)");
comp.fTree->SetAlias("dphi0","acos((TR.fX*TR.fPx+TR.fY*TR.fPy)/(TR.Pt()*TR.R()))");
hc->Fill(esd->GetTRDclusters(0));
- Float_t pt_1=TMath::Abs(out.Get1Pt());
+ Float_t pt_1=out.OneOverPt();
hpt->Fill((pt_1 - 1/ptg)/(1/ptg)*100.);
Float_t y=out.GetY();
hfound->Fill(p->Pt());
- Pt = TMath::Abs(rt->GetPt());
+ Pt = rt->Pt();
Double_t cc = TMath::Abs(rt->GetSigmaC2());
mct->GetPxPyPzXYZ(mcPx,mcPy,mcPz,x,y,z,-1);
rt->GetPxPyPz(Px,Py,Pz);
label, rt->GetSeedLabel());
return;
}
- Pt = TMath::Abs(rt->GetPt());
+ Pt = rt->GetPt();
Double_t cc = TMath::Abs(rt->GetSigmaC2());
mct->GetPxPyPzXYZ(mcPx,mcPy,mcPz,x,y,z,-1);
rt->GetPxPyPz(Px,Py,Pz);