PWG4/macros/AddTaskUE.C

   1 \r
   2 void ConfigTaskUE(AliAnalysisTaskUE * ueana); // common config, extend with different cases\r
   3                   \r
   4 AliAnalysisTaskUE *AddTaskUE(char *jetBranch = "")\r
   5 {\r
   6 // Creates a jet fider task, configures it and adds it to the analysis manager.\r
   7 \r
   8    // Get the pointer to the existing analysis manager via the static access method.\r
   9    //==============================================================================\r
  10    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
  11    if (!mgr) {\r
  12       ::Error("AddTaskUE", "No analysis manager to connect to.");\r
  13       return NULL;\r
  14    }  \r
  15    \r
  16    // Check the analysis type using the event handlers connected to the analysis manager.\r
  17    //==============================================================================\r
  18    if (!mgr->GetInputEventHandler()) {\r
  19       ::Error("AddTaskUE", "This task requires an input event handler");\r
  20       return NULL;\r
  21    }\r
  22 \r
  23    // Create the task and configure it.\r
  24    //===========================================================================\r
  25    \r
  26    AliAnalysisTaskUE* ueana = new  AliAnalysisTaskUE("Underlying Event");\r
  27    ConfigTaskUE(ueana);\r
  28 \r
  29    TString jb(jetBranch);\r
  30 \r
  31    if(jb.Length()>0){\r
  32      ueana->ReadDeltaAOD(kTRUE);\r
  33      ueana->SelectDeltaAODBranch(jb.Data());\r
  34    }\r
  35 \r
  36    mgr->AddTask(ueana);\r
  37    \r
  38    // Create ONLY the output containers for the data produced by the task.\r
  39    // Get and connect other common input/output containers via the manager as below\r
  40    //==============================================================================\r
  41    AliAnalysisDataContainer *coutput1_UE = 0;\r
  42    if(jb.Length()==0)coutput1_UE = mgr->CreateContainer("histosUE", TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWG4_UE",AliAnalysisManager::GetCommonFileName()));\r
  43    else coutput1_UE = mgr->CreateContainer(Form("histosUE_%s",jb.Data()), TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWG4_UE_%s",AliAnalysisManager::GetCommonFileName(),jb.Data()));\r
  44    \r
  45    mgr->ConnectInput  (ueana, 0, mgr->GetCommonInputContainer());\r
  46    mgr->ConnectOutput (ueana,     0, coutput1_UE );\r
  47    \r
  48    return ueana;\r
  49 }\r
  50 \r
  51 void ConfigTaskUE(AliAnalysisTaskUE * ueana){\r
  52   // common config, extend with different cases\r
  53   Int_t anaType =1; \r
  54   Int_t regType =1;\r
  55   Double_t jetEtaCut=0.2;\r
  56   Double_t trackPtCut=0.5; \r
  57   Double_t trackEtaCut= 0.9; \r
  58   Double_t rad=0.7; \r
  59   Double_t deltaPhiCut = 2.616;\r
  60   \r
  61   ueana->SetDebugLevel(0); \r
  62   ueana->SetPtRangeInHist(25, 0., 250.);\r
  63   ueana->SetAnaTopology(anaType);      \r
  64   ueana->SetRegionType(regType);        \r
  65   ueana->SetJet1EtaCut(jetEtaCut);     \r
  66   ueana->SetTrackPtCut(trackPtCut); \r
  67   ueana->SetPtSumOrdering(2);\r
  68   ueana->SetConeRadius(rad);   \r
  69   ueana->SetTrackEtaCut(trackEtaCut);\r
  70   ueana->SetJet2DeltaPhiCut(deltaPhiCut);\r
  71 }\r