From Chris: Added possibility to cut on NEF
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Sep 2013 08:22:45 +0000 (08:22 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Sep 2013 08:22:45 +0000 (08:22 +0000)
PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C

index 9464096..8da06c0 100644 (file)
@@ -24,7 +24,9 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
   Double_t    maxCent            = -999,
   const char *taskname           = "AliJetResponseMaker",
   Bool_t      biggerMatrix       = kFALSE,
-  AliJetResponseMaker* address   = 0
+  AliJetResponseMaker* address   = 0,
+  Double_t    nefmincut          = -10,
+  Double_t    nefmaxcut          = 10
 )
 {  
   // Get the pointer to the existing analysis manager via the static access method.
@@ -55,7 +57,9 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
 
   if (ptHardBin != -999) 
     name += Form("_PtHard%d", ptHardBin);
-
+  if (nefmaxcut<1.0)
+    name += Form("_NEF%d", (Int_t)(100*nefmaxcut));
+    
   AliJetResponseMaker* jetTask = address;
   if (jetTask)
     new (jetTask) AliJetResponseMaker(name);
@@ -74,7 +78,8 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
   jetCont1->SetIsParticleLevel(kFALSE);
   jetCont1->ConnectParticleContainer(trackCont1);
   jetCont1->ConnectClusterContainer(clusCont1);
-
+  jetCont1->SetNEFCut(nefmincut,nefmaxcut);
+    
   AliParticleContainer *trackCont2 = jetTask->AddParticleContainer(ntracks2);
   AliClusterContainer *clusCont2 = jetTask->AddClusterContainer(nclusters2);
   AliJetContainer *jetCont2 = jetTask->AddJetContainer(njets2, cutType, jetradius2);
@@ -87,7 +92,7 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
   jetCont2->SetIsParticleLevel(kTRUE);
   jetCont2->ConnectParticleContainer(trackCont2);
   jetCont2->ConnectClusterContainer(clusCont2);
-
+    
   jetTask->SetMatching(matching, maxDistance1, maxDistance2);
   jetTask->SetVzRange(-10,10);
   jetTask->SetIsPythia(kTRUE);