PWGHF/hfe/macros/AddTaskHFECal.C

   1 AliAnalysisTask *AddTaskHFECal()
   2 {
   3   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   4   if (!mgr) {
   5     Error("AddTaskHFEECal", "No analysis manager found.");
   6     return NULL;
   7   }
   8
   9   if (!mgr->GetInputEventHandler()) {
  10     ::Error("AddTaskHFECal", "This task requires an input event handler");
  11     return NULL;
  12   }
  13   TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
  14   if (type=="AOD"){
  15     ::Error("AddTaskHFECal", "The tasks exits because AODs are in input");
  16     return NULL;
  17   }
  18   Bool_t MCthere=kFALSE;
  19   AliMCEventHandler *mcH = dynamic_cast<AliMCEventHandler*>(mgr->GetMCtruthEventHandler());
  20   if(!mcH){
  21     MCthere=kFALSE;
  22   }else{
  23     MCthere=kTRUE;
  24   }
  25
  26
  27   //analysis task
  28   gROOT->LoadMacro("$ALICE_ROOT/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C");
  29
  30   AliAnalysisTaskHFECal *hfetaskCent = ConfigHFECal(MCthere);
  31   AliAnalysisTaskHFECal *hfetaskTrig= ConfigHFECal(MCthere);
  32   AliAnalysisTaskHFECal *hfetaskTrig2= ConfigHFECal(MCthere);
  33
  34   mgr->AddTask(hfetaskCent);
  35   mgr->AddTask(hfetaskTrig);
  36   mgr->AddTask(hfetaskTrig2);
  37
  38   // central trigger
  39   hfetaskCent->SelectCollisionCandidates(AliVEvent::kCentral);
  40
  41   TString containerName = mgr->GetCommonFileName();
  42   containerName += ":PWGHF_hfeCalCentral";
  43
  44   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  45   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_EMCalCentral", TList::Class(),AliAnalysisManager::kOutputContainer, containerName.Data());
  46   mgr->ConnectInput(hfetaskCent, 0, cinput);
  47   mgr->ConnectOutput(hfetaskCent, 1, coutput1);
  48
  49   //L1 gamma trigger
  50   hfetaskTrig->SelectCollisionCandidates(AliVEvent::kEMCEGA);
  51
  52   TString containerName2 = mgr->GetCommonFileName();
  53   containerName2 += ":PWGHF_hfeCalTrigEGA";
  54
  55   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  56   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_EMCalTrigEGA", TList::Class(),AliAnalysisManager::kOutputContainer, containerName2.Data());
  57   mgr->ConnectInput(hfetaskTrig, 0, cinput);
  58   mgr->ConnectOutput(hfetaskTrig, 1, coutput1);
  59
  60
  61   //Jet trigger
  62   // semi-central
  63   hfetaskTrig2->SelectCollisionCandidates(AliVEvent::kSemiCentral);
  64
  65   TString containerName3 = mgr->GetCommonFileName();
  66   containerName3 += ":PWGHF_hfeCalSemiCentral";
  67
  68   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  69   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_SemiCentral", TList::Class(),AliAnalysisManager::kOutputContainer, containerName3.Data());
  70   mgr->ConnectInput(hfetaskTrig2, 0, cinput);
  71   mgr->ConnectOutput(hfetaskTrig2, 1, coutput1);
  72
  73
  74   //if(MCthere)
  75     //{
  76     //MB trigger
  77     AliAnalysisTaskHFECal *hfetaskMB = ConfigHFECal(MCthere);
  78     mgr->AddTask(hfetaskMB);
  79     hfetaskMB->SelectCollisionCandidates(AliVEvent::kMB);
  80
  81     TString containerName3 = mgr->GetCommonFileName();
  82     containerName3 += ":PWGHF_hfeCalkMB";
  83
  84      AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  85      AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_EMCalMB", TList::Class(),AliAnalysisManager::kOutputContainer, containerName3.Data());
  86      mgr->ConnectInput(hfetaskMB, 0, cinput);
  87      mgr->ConnectOutput(hfetaskMB, 1, coutput1);
  88      //}
  89
  90
  91   return NULL;
  92 }