use extrapolation fct of AliEMCALRecoUtils
authormverweij <marta.verweij@cern.ch>
Wed, 11 Dec 2013 08:49:26 +0000 (09:49 +0100)
committerPeter Hristov <phristov@pb-d-128-141-110-112.cern.ch>
Wed, 11 Dec 2013 10:22:46 +0000 (11:22 +0100)
PWG/EMCAL/AliEmcalEsdTrackFilterTask.cxx
PWG/EMCAL/AliEmcalPicoTrackMaker.cxx
PWG/EMCAL/macros/AddTaskEmcalEsdTrackFilter.C [moved from PWG/EMCAL/macros/AddTaskEmcalEsdTpcTrack.C with 95% similarity]

index 3bea7e0..810402a 100644 (file)
@@ -174,7 +174,7 @@ void AliEmcalEsdTrackFilterTask::UserExec(Option_t *)
         AliESDtrack *newTrack = static_cast<AliESDtrack*>(fTracks->At(ntrnew));
        if (fDoPropagation) 
          AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(newTrack,fDist);
-        //newTrack->SetTRDNchamberdEdx(0); todo
+        newTrack->SetTRDNchamberdEdx(0); //todo
         ++ntrnew;
       } else if (fHybridTrackCuts->AcceptTrack(etrack)) {
 
@@ -187,10 +187,10 @@ void AliEmcalEsdTrackFilterTask::UserExec(Option_t *)
                         constrainParam->GetAlpha(),
                         constrainParam->GetParameter(),
                         constrainParam->GetCovariance());
-         //          if ((status&AliESDtrack::kITSrefit)==0) // todo
-         //   newTrack->SetTRDNchamberdEdx(2);
-         // else
-          //  newTrack->SetTRDNchamberdEdx(1);
+         if ((status&AliESDtrack::kITSrefit)==0) // todo
+           newTrack->SetTRDNchamberdEdx(2);
+         else
+            newTrack->SetTRDNchamberdEdx(1);
          if (fDoPropagation)   
            AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(newTrack,fDist);
           ++ntrnew;
index 91be15f..0b01816 100644 (file)
@@ -146,7 +146,11 @@ void AliEmcalPicoTrackMaker::UserExec(Option_t *)
       type = esdtrack->GetTRDNchamberdEdx();
       if (!fIncludeNoITS && (type==2))
        continue;
-      isEmc = track->IsEMCAL();
+      if (TMath::Abs(esdtrack->GetTrackEtaOnEMCal()) < 0.75 && 
+         esdtrack->GetTrackPhiOnEMCal() > 70 * TMath::DegToRad() &&
+         esdtrack->GetTrackPhiOnEMCal() < 190 * TMath::DegToRad())
+       isEmc = kTRUE;
+      //   isEmc = track->IsEMCAL();
     } else {
       AliAODTrack *aodtrack = static_cast<AliAODTrack*>(track);
       if (fAODfilterBits[0] < 0) {
similarity index 95%
rename from PWG/EMCAL/macros/AddTaskEmcalEsdTpcTrack.C
rename to PWG/EMCAL/macros/AddTaskEmcalEsdTrackFilter.C
index 0412c31..5aa6211 100644 (file)
@@ -1,7 +1,7 @@
 // $Id$
 
-AliEmcalEsdTpcTrackTask* AddTaskEmcalEsdTpcTrack(
-  const char *name              = "TpcSpdVertexConstrainedTracks",
+AliEmcalEsdTrackFilterTask* AddTaskEmcalEsdTrackFilter(
+  const char *name              = "AliEmcalEsdTrackFilter",
   const char *trackCuts         = "Hybrid_LHC11h"
 )
 { 
@@ -119,7 +119,7 @@ AliEmcalEsdTpcTrackTask* AddTaskEmcalEsdTpcTrack(
 
   gROOT->LoadMacro("$ALICE_ROOT/PWGJE/macros/CreateTrackCutsPWGJE.C");
   Bool_t includeNoITS = kFALSE;
-  AliEmcalEsdTpcTrackTask *eTask = new AliEmcalEsdTpcTrackTask(); // default is TPC only tracks constrained to the vertex
+  AliEmcalEsdTrackFilterTask *eTask = new AliEmcalEsdTrackFilterTask(); // default is TPC only tracks constrained to the vertex
 
   if ((dataSet == kLHC11c && cutsType == kHybrid) ||
       (dataSet == kLHC11d && cutsType == kHybrid) ||