from Redmer

[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / macros / AddTaskChargedJetsPA.C

diff --git a/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsPA.C b/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsPA.C
index 32ae938bdf0814c7222305566b14088f681e2825..2acf9cc79bd50d5c07cbdff21d75e347ea619b77 100644 (file)
--- a/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsPA.C
+++ b/PWGJE/EMCALJetTasks/macros/AddTaskChargedJetsPA.C
@@ -5,15 +5,16 @@ AliAnalysisTaskChargedJetsPA* AddTaskChargedJetsPA(
   Double_t            randomConeR             = 0.4,
   Double_t            trackBgrdConeR          = 0.6,
   const char*         usedTracks              = "PicoTracks",
-  const char*         usedClusters            = "CaloClustersCorr",
+  const char*         centralityType          = "V0A",
   Double_t            trackEtaWindow          = 0.9,
   Double_t            vertexWindow            = 10.0,
   Double_t            vertexMaxR              = 1.0,
-  Int_t               minVertexContributors   = 1,
   Double_t            minJetPt                = 5.0, // signal jet min pt
   Double_t            dijetLeadingMinPt       = 10.0,
   Double_t            dijetMaxAngleDev        = 10.0,
-  Int_t               numberOfPtHardBins      = 0
+  Int_t               numberOfPtHardBins      = 0,
+  const char*         fileEtaCorrectionFactors= "alien:///alice/cern.ch/user/r/rhaake/pA/EtaCorrectionFactors.root",
+  const char*         externalMacro           = NULL
 )
 {
   // #### Detect the demanded trigger with its readable name
@@ -24,6 +25,8 @@ AliAnalysisTaskChargedJetsPA* AddTaskChargedJetsPA(
     triggerName = "kAny";
   else if(trigger == AliVEvent::kINT7)
     triggerName = "kINT7";
+  else if(trigger == AliVEvent::kMB)
+    triggerName = "kMB";
   else if(trigger == AliVEvent::kEMC7)
     triggerName = "kEMC7";
   else if(trigger == AliVEvent::kEMCEJE)
@@ -32,7 +35,6 @@ AliAnalysisTaskChargedJetsPA* AddTaskChargedJetsPA(
     triggerName = "kEMCEGA";
 
   // #### Define manager and data container names
-  const char*         kFileName               = "ChargedJetsPA.root"; // hard coded
   AliAnalysisManager *manager = AliAnalysisManager::GetAnalysisManager();
   if (!manager) {
     ::Error("AddTaskChargedJetsPA", "No analysis manager to connect to.");
@@ -40,29 +42,30 @@ AliAnalysisTaskChargedJetsPA* AddTaskChargedJetsPA(
   }
   TString myContName("");
   if(isMC)
-    myContName = Form("MCChargedJets_pA_R0%2.0f_%s",jetRadius*100,triggerName.Data());
+    myContName = Form("AnalysisR0%2.0f_%s_MC",jetRadius*100,triggerName.Data());
   else
-    myContName = Form("ChargedJets_pA_R0%2.0f_%s",jetRadius*100,triggerName.Data());
+    myContName = Form("AnalysisR0%2.0f_%s",jetRadius*100,triggerName.Data());
 
   // #### Add necessary jet finder tasks
   gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
-  AliEmcalJetTask* jetFinderTask = AddTaskEmcalJet(usedTracks,"",1,jetRadius,1,0.150,0.300);// anti-kt
+  AliEmcalJetTask* jetFinderTask = AddTaskEmcalJet(usedTracks,"",1,jetRadius,1,0.150,0.300); // anti-kt
   AliEmcalJetTask* jetFinderTaskKT = AddTaskEmcalJet(usedTracks,"",0,jetRadius,1,0.150,0.300); // kt
 
   // #### Define analysis task
   AliAnalysisTaskChargedJetsPA *task = NULL;
-  contHistos = manager->CreateContainer(myContName.Data(), TList::Class(), AliAnalysisManager::kOutputContainer, kFileName);
-  task = new AliAnalysisTaskChargedJetsPA(Form("Analysis_pA_%s_%s", jetFinderTask->GetName(), triggerName.Data()), usedTracks, usedClusters, jetFinderTask->GetName(),jetFinderTaskKT->GetName());
+  contHistos = manager->CreateContainer(myContName.Data(), TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:ChargedJetsPA", AliAnalysisManager::GetCommonFileName()));
+  task = new AliAnalysisTaskChargedJetsPA(Form("AnalysisPA_%s_%s", jetFinderTask->GetName(), triggerName.Data()), usedTracks, jetFinderTask->GetName(),jetFinderTaskKT->GetName());
 
   // #### Task preferences
-  task->SetAcceptanceWindows(trackEtaWindow, vertexWindow, vertexMaxR, minVertexContributors, jetRadius, jetRadius);
+  task->SetAcceptanceWindows(trackEtaWindow, vertexWindow, vertexMaxR, jetRadius, jetRadius);
   task->SetSignalJetMinPt(minJetPt);
   task->SetSignalJetMinArea(0.6*jetRadius*jetRadius*TMath::Pi());
   task->SetDijetLeadingMinPt(dijetLeadingMinPt);
   task->SetDijetMaxAngleDeviation(dijetMaxAngleDev);
   task->SetRandConeRadius(randomConeR);
-  task->SetTrackBackgroundConeRadius(trackBgrdConeR);
+  task->SetTRBackgroundConeRadius(trackBgrdConeR);
   task->SelectCollisionCandidates(trigger);
+  task->SetCentralityType(centralityType);
   if(numberOfPtHardBins)
     task->SetNumberOfPtHardBins(numberOfPtHardBins);
 
@@ -70,5 +73,11 @@ AliAnalysisTaskChargedJetsPA* AddTaskChargedJetsPA(
   manager->AddTask(task);
   manager->ConnectInput(task, 0, manager->GetCommonInputContainer());
   manager->ConnectOutput(task, 1, contHistos);
+
+  // #### Do some nasty piggybacking on demand
+  if (externalMacro)
+    gROOT->Macro(externalMacro);
+
+
   return task;
 }