Add the efficiency estimation to common FF task.
[u/mrichter/AliRoot.git] / PWG4 / macros / AddTaskFragmentationFunction.C
index 4a3aa65..169ee4d 100644 (file)
@@ -130,18 +130,27 @@ AliAnalysisTaskFragmentationFunction *AddTaskFragmentationFunction(
    else if(typeJets.Contains("jetTypeUndef")) task->SetJetTypeGen(0); // undefined
    else Printf("jetType %s not found", typeJets.Data());
    
+   if(typeJets.Contains("AODMCb")) task->SetJetTypeRecEff(AliAnalysisTaskFragmentationFunction::kJetsGenAcceptance); // kJetsRecAcceptance
+   else if(typeJets.Contains("AODb")) task->SetJetTypeRecEff(AliAnalysisTaskFragmentationFunction::kJetsRecAcceptance); 
+   else task->SetJetTypeRecEff(0);
+
    task->SetFilterMask(filterMask);
   
-   // set default parameter 
+   // Set default parameters 
+   // Cut selection 
    task->SetTrackCuts();       // default : pt > 0.150 GeV, |eta|<0.9, full phi acc
    task->SetJetCuts();         // default: jet pt > 5 GeV, |eta|<0.5, full phi acc
-   task->SetDiJetCuts();       // default: to be defined
+   task->SetDiJetCuts();       // default: type of cut = 1 (cut in deltaPhi), deltaPhi = 0., cdf = 0.5, fraction of pt = 0.6
+   task->SetKindSlices();      // default: kindSlice = 1 (inv mass)
    task->SetFFRadius();        // default: R = 0.4
-   
    task->SetHighPtThreshold(); // default: pt > 5 Gev
+   // Define histo bins
    task->SetFFHistoBins();
    task->SetQAJetHistoBins();
    task->SetQATrackHistoBins();
+   task->SetIJHistoBins();
+   task->SetDiJetHistoBins();
+   task->SetQADiJetHistoBins();
 
    mgr->AddTask(task);