Call FindFixBin instead of FindBin for eta correction maps (NO impact
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 28 Nov 2013 15:34:03 +0000 (15:34 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 28 Nov 2013 15:34:03 +0000 (15:34 +0000)
on any analysis, since bit TAxis::kCanExtend is kFALSE for the maps; but
FindFixBin faster and safer (if bit is changed in future))

Benjamin Hess

STEER/STEERBase/AliTPCPIDResponse.cxx

index 779be05..ff8be78 100644 (file)
@@ -785,8 +785,8 @@ Double_t AliTPCPIDResponse::GetEtaCorrectionFast(const AliVTrack *track, Double_
     return 1.;
   
   Double_t tanTheta = GetTrackTanTheta(track); 
-  Int_t binX = fhEtaCorr->GetXaxis()->FindBin(tanTheta);
-  Int_t binY = fhEtaCorr->GetYaxis()->FindBin(1. / tpcSignal);
+  Int_t binX = fhEtaCorr->GetXaxis()->FindFixBin(tanTheta);
+  Int_t binY = fhEtaCorr->GetYaxis()->FindFixBin(1. / tpcSignal);
   
   if (binX == 0) 
     binX = 1;
@@ -907,8 +907,8 @@ Double_t AliTPCPIDResponse::GetSigmaPar1Fast(const AliVTrack *track, AliPID::EPa
     return 999;
   
   Double_t tanTheta = GetTrackTanTheta(track);
-  Int_t binX = fhEtaSigmaPar1->GetXaxis()->FindBin(tanTheta);
-  Int_t binY = fhEtaSigmaPar1->GetYaxis()->FindBin(1. / dEdxExpected);
+  Int_t binX = fhEtaSigmaPar1->GetXaxis()->FindFixBin(tanTheta);
+  Int_t binY = fhEtaSigmaPar1->GetYaxis()->FindFixBin(1. / dEdxExpected);
     
   if (binX == 0) 
     binX = 1;