- Fixed TPC multiplicity estimator (nTotTracks -> nTotESDTracks); not an issue so...
authorbhess <bhess@cern.ch>
Mon, 30 Jun 2014 07:57:07 +0000 (09:57 +0200)
committermivanov <marian.ivanov@cern.ch>
Mon, 30 Jun 2014 20:14:53 +0000 (22:14 +0200)
PWGPP/TPC/AliTPCPIDEtaQA.cxx
PWGPP/TPC/AliTPCPIDEtaTree.cxx
PWGPP/TPC/AliTPCcalibResidualPID.cxx

index 3aa24ab..9bd5586 100644 (file)
@@ -220,7 +220,7 @@ void AliTPCPIDEtaQA::UserExec(Option_t *)
   // Fill V0 arrays with V0s
   FillV0PIDlist();
   
-  Int_t multiplicity = fESD->GetNumberOfTracks();
+  Int_t multiplicity = fESD->GetNumberOfESDTracks();
   
   // Track loop to fill a Train spectrum
   for (Int_t iTracks = 0; iTracks < fESD->GetNumberOfTracks(); iTracks++) {
index 194e945..269e993 100644 (file)
@@ -244,10 +244,10 @@ void AliTPCPIDEtaTree::UserExec(Option_t *)
   if(primaryVertex->GetNContributors() <= 0) 
     return;
   
-  fMultiplicity = fESD->GetNumberOfTracks(); 
+  fMultiplicity = fESD->GetNumberOfESDTracks(); 
     
   if (fDoAdditionalQA) {
-    Int_t nTotTracks = fESD->GetNumberOfTracks();
+    Int_t nTotTracks = fESD->GetNumberOfESDTracks();
     fhMultiplicityQA->Fill(primaryVertex->GetNContributors(), nTotTracks);
   }
   
index 13a1073..6ec06cb 100644 (file)
@@ -392,6 +392,8 @@ void AliTPCcalibResidualPID::Process(AliESDEvent *const esdEvent, AliMCEvent *co
     v0QAadded[i] = kFALSE;
   
   Int_t nTotTracks = esdEvent->GetNumberOfTracks();
+  const Int_t nTotESDTracks = esdEvent->GetNumberOfESDTracks();
+  
   for (Int_t iTracks = 0; iTracks < nTotTracks; iTracks++){//begin track loop 
     AliESDtrack *trackESD = esdEvent->GetTrack(iTracks);
     if(!trackESD) {
@@ -658,7 +660,7 @@ void AliTPCcalibResidualPID::Process(AliESDEvent *const esdEvent, AliMCEvent *co
         continue;
       
       Double_t vecHistQA[7] = {precin, processedTPCsignal[iPart], (Double_t)particleID, (Double_t)iPart, tpcQA[iPart], tofQA[iPart],
-                               (Double_t)nTotTracks};
+                               (Double_t)nTotESDTracks};
       if (fProduceGlobal) fHistPidQA->Fill(vecHistQA);
       vecHistQA[1] = signal[0]; vecHistQA[4] = ncl[0];
       if ((fProduceAllPadTypes || fProduceShortPads)) fHistPidQAshort->Fill(vecHistQA);