PWGJE/macros/AddTaskJetHadronCorrelation.C

   1 AliAnalysisTaskJetHadronCorrelation *AddTaskJetHadronCorrelation(Char_t *jf="ANTIKT",Float_t radius=0.4,UInt_t filter=256,Int_t backM=0,Float_t tPtcut=0.15,Int_t skipCone=0,Bool_t IsMC=true,Float_t JetEScale=1.,Float_t TrackEScale=1.)\r
   2 {\r
   3 \r
   4    // Creates a JetQA task, configures it and adds it to the analysis manager.\r
   5    // Get the pointer to the existing analysis manager via the static access method.\r
   6    //==============================================================================\r
   7    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
   8    if (!mgr) {\r
   9       ::Error("AddTaskJetHadronCorrelation", "No analysis manager to connect to.");\r
  10       return NULL;\r
  11    }\r
  12 \r
  13    // Check the analysis type using the event handlers connected to the analysis manager.\r
  14    //==============================================================================\r
  15    if (!mgr->GetInputEventHandler()) {\r
  16       ::Error("AddTaskJetHadronCorrelation", "This task requires an input event handler");\r
  17       return NULL;\r
  18    }\r
  19 \r
  20    // Create the task and configure it.\r
  21    //===========================================================================\r
  22 \r
  23    AliAnalysisTaskJetHadronCorrelation *jethadron = new AliAnalysisTaskJetHadronCorrelation("TaskJetHadronCorrelation");\r
  24    jethadron->SetDebugLevel(3);\r
  25          jethadron->SetAlgorithm(jf);\r
  26          jethadron->SetRadius(radius);\r
  27          jethadron->SetFilterMask(filter);\r
  28          jethadron->SetBackSubMode(backM);\r
  29          jethadron->SetTrackPtCut(tPtcut);\r
  30          jethadron->SetSkipCone(skipCone);\r
  31          jethadron->SetMC(IsMC);\r
  32          jethadron->SetJetEScale(JetEScale);\r
  33          jethadron->SetTrackEScale(TrackEScale);\r
  34          mgr->AddTask(jethadron); \r
  35 \r
  36          TString cAdd = "";\r
  37          cAdd += Form("%02d_"           ,(int)((radius+0.01)*10.));\r
  38          cAdd += Form("B%d"             ,(int)backM);\r
  39          cAdd += Form("_Filter%05d"     ,filter);\r
  40          cAdd += Form("_Cut%05d"        ,(int)(1000.*tPtcut));\r
  41          cAdd += Form("_Skip%02d"       ,skipCone);\r
  42          cAdd += Form("_JetEScale%03d"  ,(int)(JetEScale*100.));\r
  43          cAdd += Form("_TrackEScale%03d",(int)(TrackEScale*100.));\r
  44          TString Branch;\r
  45          if(IsMC)Branch = Form("MC_clustersAOD_%s%s",jf,cAdd.Data());\r
  46          else    Branch = Form("Data_clustersAOD_%s%s",jf,cAdd.Data());\r
  47 \r
  48    AliAnalysisDataContainer *cout_JHC = mgr->CreateContainer(Form("JHC_%s%s",jf,cAdd.Data()), TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWGJE_JHC_%s",AliAnalysisManager::GetCommonFileName(),Branch.Data()));\r
  49 \r
  50    // Create ONLY the output containers for the data produced by the task.\r
  51    // Get and connect other common input/output containers via the manager as below\r
  52    //==============================================================================\r
  53 \r
  54    mgr->ConnectInput (jethadron,0, mgr->GetCommonInputContainer());  \r
  55    mgr->ConnectOutput(jethadron,0, mgr->GetCommonOutputContainer());\r
  56    mgr->ConnectOutput(jethadron,1, cout_JHC);\r
  57 \r
  58 \r
  59    return jethadron;\r
  60 }\r