Fix bug in total probability
authorcblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 7 Mar 2007 14:59:23 +0000 (14:59 +0000)
committercblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 7 Mar 2007 14:59:23 +0000 (14:59 +0000)
TRD/AliTRDpidESD.cxx

index 55a3768..c16ef19 100644 (file)
@@ -41,7 +41,7 @@ ClassImp(AliTRDpidESD)
 
   Bool_t AliTRDpidESD::fCheckTrackStatus = kTRUE;
   Bool_t AliTRDpidESD::fCheckKinkStatus  = kFALSE;
-  Int_t AliTRDpidESD::fMinPlane         = 0;
+  Int_t AliTRDpidESD::fMinPlane          = 0;
 
 //_____________________________________________________________________________
 AliTRDpidESD::AliTRDpidESD():TObject()
@@ -50,7 +50,6 @@ AliTRDpidESD::AliTRDpidESD():TObject()
   // Default constructor
   //
 
-
 }
 
 //_____________________________________________________________________________
@@ -173,13 +172,15 @@ Int_t AliTRDpidESD::MakePID(AliESD *event)
          p[iSpecies] *= pd->GetProbabilityT(iSpecies,mom,timebin);
          p[iSpecies] *= 100.0; // ??????????????
 
-          probTotal   += p[iSpecies];
-
        }
 
       }
 
-    } 
+    }
+
+    for (Int_t iSpecies = 0; iSpecies < nSpecies; iSpecies++) {
+      probTotal += p[iSpecies];
+    }
 
     for (Int_t iSpecies = 0; iSpecies < nSpecies; iSpecies++) {
       if ((probTotal >       0.0) &&