From: saiola <salvatore.aiola@cern.ch>
Date: Wed, 11 Jun 2014 15:41:43 +0000 (-0400)
Subject: random track rejection (tracking efficiency studies) only after all cuts have been... 
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=c04400125176f83a8af92dff26905b7832a3065d;p=u%2Fmrichter%2FAliRoot.git

random track rejection (tracking efficiency studies) only after all cuts have been applied
---

diff --git a/PWG/EMCAL/AliEmcalEsdTrackFilterTask.cxx b/PWG/EMCAL/AliEmcalEsdTrackFilterTask.cxx
index a0052ba9c3c..5b18ce68ec4 100644
--- a/PWG/EMCAL/AliEmcalEsdTrackFilterTask.cxx
+++ b/PWG/EMCAL/AliEmcalEsdTrackFilterTask.cxx
@@ -129,12 +129,6 @@ void AliEmcalEsdTrackFilterTask::UserExec(Option_t *)
         if (!etrack)
           continue;
 
-	if (fTrackEfficiency < 1) {
-	  Double_t r = gRandom->Rndm();
-	  if (fTrackEfficiency < r)
-	    continue;
-	}
-
         if (!fEsdTrackCuts->AcceptTrack(etrack))
           continue;
 
@@ -159,6 +153,13 @@ void AliEmcalEsdTrackFilterTask::UserExec(Option_t *)
           delete ntrack;
           continue;
         }
+
+	if (fTrackEfficiency < 1) {
+	  Double_t r = gRandom->Rndm();
+	  if (fTrackEfficiency < r)
+	    continue;
+	}
+
 	if (fDoPropagation) 	
 	  AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(ntrack,fDist);
         new ((*fTracks)[ntrnew++]) AliESDtrack(*ntrack);
@@ -170,13 +171,16 @@ void AliEmcalEsdTrackFilterTask::UserExec(Option_t *)
         AliESDtrack *etrack = fEsdEv->GetTrack(i);
         if (!etrack)
           continue;
+
+	if ((fEsdTrackCuts!=0) && !fEsdTrackCuts->AcceptTrack(etrack))
+          continue;
+	
 	if (fTrackEfficiency < 1) {
 	  Double_t r = gRandom->Rndm();
 	  if (fTrackEfficiency < r)
 	    continue;
 	}
-	if ((fEsdTrackCuts!=0) && !fEsdTrackCuts->AcceptTrack(etrack))
-          continue;
+
         AliESDtrack *ntrack = new ((*fTracks)[ntrnew++]) AliESDtrack(*etrack);
 	if (fDoPropagation) 	
 	  AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(ntrack,fDist);
@@ -192,12 +196,6 @@ void AliEmcalEsdTrackFilterTask::UserExec(Option_t *)
       if (!etrack) 
 	continue;
 
-      if (fTrackEfficiency < 1) {
-	Double_t r = gRandom->Rndm();
-	if (fTrackEfficiency < r)
-	  continue;
-      }
-
       if (fEsdTrackCuts->AcceptTrack(etrack)) {
         AliESDtrack *newTrack = new ((*fTracks)[ntrnew]) AliESDtrack(*etrack);
 	if (fDoPropagation) 
@@ -211,6 +209,13 @@ void AliEmcalEsdTrackFilterTask::UserExec(Option_t *)
 	UInt_t status = etrack->GetStatus();
 	if (!fIncludeNoITS && ((status&AliESDtrack::kITSrefit)==0))
 	  continue;
+
+	if (fTrackEfficiency < 1) {
+	  Double_t r = gRandom->Rndm();
+	  if (fTrackEfficiency < r)
+	    continue;
+	}
+
 	AliESDtrack *newTrack = new ((*fTracks)[ntrnew]) AliESDtrack(*etrack);
 	if (fDoPropagation) 	
 	  AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(newTrack,fDist);