PWG4/macros/AddTaskJetChem.C

   1 AliAnalysisTask *AddTaskJetChem(){
   2
   3
   4
   5   //cout<<" OB : add JetTasks inlcude path ! "<<endl;
   6   //gSystem->AddIncludePath("-I$ALICE_ROOT/PWG4/JetTasks");
   7
   8   //get the current analysis manager
   9   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  10   if (!mgr) {
  11     Error("AddTask_obusch_jets", "No analysis manager found.");
  12     return 0;
  13   }
  14
  15   // physics event selection task - not for AOD analysis
  16
  17   //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
  18   //AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
  19   //AliPhysicsSelection* physSel = physSelTask->GetPhysicsSelection();
  20   //physSel->AddBackgroundIdentification(new AliBackgroundSelection());
  21   //physSel->SetAnalyzeMC();
  22
  23   /*
  24   //============= Set Task Name ===================
  25   TString taskName=("AliAnalysisTaskJetChem.cxx+");
  26   //===============================================
  27   //            Load the task
  28   gROOT->LoadMacro(taskName.Data());
  29   if (gProof){
  30     TString taskSO=gSystem->pwd();
  31     taskSO+="/";
  32     taskSO+=taskName(0,taskName.First('.'))+"_cxx.so";
  33     gProof->Exec(Form("gSystem->Load(\"%s\")",taskSO.Data()),kTRUE);
  34   }
  35   */
  36   //========= Add task to the ANALYSIS manager =====
  37   AliAnalysisTaskJetChem *task = new AliAnalysisTaskJetChem;
  38
  39   // configure task
  40
  41
  42   task->SetUseLOConeJets();
  43   task->SetUseLOConeMCJets();
  44
  45   //task->SetUsePythiaJets();
  46   task->SetConeRadius(0.4);
  47   task->SetTrackPtCutJF(0.150); //
  48   task->SetFilterBitJF(0x01);   // official PWG4 high pt filter bit 0x10, but not all AliAnalysisTaskESDFilter had configured this ESDTrackCut
  49   task->SetRequireITSRefitJF(); // 0x01 + ITS refit = 0x10
  50   task->SetRejectK0TracksJF();  // uncomment for K0 analysis running jet finder in task - modifies jet spectrum
  51
  52   task->SetJetPtCut(2.0);
  53   //task->SetJetPtCut(0.150); // lower pt cut: for plot of diffractive contribution to jet spectrum (goes up to 2 GeV ...)
  54   task->SetJetEtaCut(0.5);
  55
  56   task->SetFilterBit((UInt_t) 0X01); // std AOD track cuts
  57   task->SetTrackPtCut(0.150);
  58   task->SetTrackEtaCut(0.9);
  59
  60   task->SetUseOnFlyV0s();
  61   task->SetCutnSigdEdx(2);
  62
  63   //task->ReadDeltaAOD(); // uncomment for DeltaAODs
  64   //task->SelectDeltaAODBranch("bla");
  65   //task->SelectAODBranch("jetsAOD_FASTKT04");
  66   //task->SelectAODBranch("jetsAOD_UA107");
  67   //task->SelectAODBranch("jets");
  68
  69   task->SelectCollisionCandidates(); // either here or in userExec of task - but not for AODs ...
  70
  71   //  AliLog::SetGlobalLogLevel(AliLog::kInfo);    // kInfo // kDebug // kFatal
  72   //  task->SetDebugLevel(10);
  73   //  mgr->SetDebugLevel(10);
  74
  75   mgr->AddTask(task);
  76
  77
  78   //================================================
  79   //              data containers
  80   //================================================
  81   //            find input container
  82   //below the trunk version
  83   AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();
  84   AliAnalysisDataContainer *coutput = mgr->CreateContainer("PWG4_JetChem",
  85                                                            TList::Class(), AliAnalysisManager::kOutputContainer,
  86                                                            Form("%s:PWG4_JetChem",AliAnalysisManager::GetCommonFileName()));
  87
  88   mgr->ConnectInput(task,0,cinput );
  89   mgr->ConnectOutput(task,1,coutput);
  90
  91   return task;
  92 }