Method added.

[u/mrichter/AliRoot.git] / TPC / AliTPCcalibTimeGain.cxx

diff --git a/TPC/AliTPCcalibTimeGain.cxx b/TPC/AliTPCcalibTimeGain.cxx
index 684bd4c6da7382aa5d4a3cc0f015b6aa79627d09..e03c9c5a67d92ac6222b9c4c00e8c27e8be1d2ed 100644 (file)
--- a/TPC/AliTPCcalibTimeGain.cxx
+++ b/TPC/AliTPCcalibTimeGain.cxx
@@ -314,9 +314,10 @@ void AliTPCcalibTimeGain::ProcessCosmicEvent(AliESDEvent *event) {
 
     AliESDtrack *track = event->GetTrack(i);
     if (!track) continue;
-        
+    AliESDfriendTrack *friendTrack = esdFriend->GetTrack(i);
+    if (!friendTrack) continue;        
     const AliExternalTrackParam * trackIn = track->GetInnerParam();
-    const AliExternalTrackParam * trackOut = track->GetOuterParam();
+    const AliExternalTrackParam * trackOut = friendTrack->GetTPCOut();
     if (!trackIn) continue;
     if (!trackOut) continue;
 
@@ -331,8 +332,6 @@ void AliTPCcalibTimeGain::ProcessCosmicEvent(AliESDEvent *event) {
     if (TMath::Abs(trackIn->GetSnp()) > 0.6) continue;
     
     // Get seeds
-    AliESDfriendTrack *friendTrack = esdFriend->GetTrack(i);
-    if (!friendTrack) continue;
     TObject *calibObject;
     AliTPCseed *seed = 0;
     for (Int_t l=0;(calibObject=friendTrack->GetCalibObject(l));++l) {
@@ -378,9 +377,11 @@ void AliTPCcalibTimeGain::ProcessBeamEvent(AliESDEvent *event) {
 
     AliESDtrack *track = event->GetTrack(i);
     if (!track) continue;
+    AliESDfriendTrack *friendTrack = esdFriend->GetTrack(i);
+    if (!friendTrack) continue;
         
     const AliExternalTrackParam * trackIn = track->GetInnerParam();
-    const AliExternalTrackParam * trackOut = track->GetOuterParam();
+    const AliExternalTrackParam * trackOut = friendTrack->GetTPCOut();
     if (!trackIn) continue;
     if (!trackOut) continue;
 
@@ -395,8 +396,6 @@ void AliTPCcalibTimeGain::ProcessBeamEvent(AliESDEvent *event) {
     if (TMath::Abs(trackIn->GetSnp()) > 0.6) continue;
     
     // Get seeds
-    AliESDfriendTrack *friendTrack = esdFriend->GetTrack(i);
-    if (!friendTrack) continue;
     TObject *calibObject;
     AliTPCseed *seed = 0;
     for (Int_t l=0;(calibObject=friendTrack->GetCalibObject(l));++l) {
@@ -404,13 +403,13 @@ void AliTPCcalibTimeGain::ProcessBeamEvent(AliESDEvent *event) {
     }    
 
     if (seed) { 
-      Double_t tpcSignal = GetTPCdEdx(seed);
-      fHistDeDxTotal->Fill(meanP, tpcSignal);
-      //
       if (fLowMemoryConsumption) {
        if (meanP > 0.5 || meanP < 0.4) continue;
        meanP = 0.45; // set all momenta to one in order to save memory
       }
+      Double_t tpcSignal = GetTPCdEdx(seed);
+      fHistDeDxTotal->Fill(meanP, tpcSignal);
+      //
       //dE/dx, time, type (1-muon cosmic,2-pion beam data), momenta
       Double_t vec[6] = {tpcSignal,time,2,meanDrift,meanP,runNumber};
       fHistGainTime->Fill(vec);