PWGJE/EMCALJetTasks/macros/AddTaskJetMassResponseDet.C

   1 AliAnalysisTaskJetMassResponseDet* AddTaskJetMassResponseDet(const char * njetsPart,
   2                                                              const char * njetsDet,
   3                                                              const Double_t R,
   4                                                              const char *type,
   5                                                              Int_t       pSel,
   6                                                              TString     kEmcalTriggers = "",
   7                                                              TString     tag            = "") {
   8
   9   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  10   if (!mgr)
  11     {
  12       Error("AddTaskEmcalJetMass","No analysis manager found.");
  13       return 0;
  14     }
  15   Bool_t ismc=kFALSE;
  16   ismc = (mgr->GetMCtruthEventHandler())?kTRUE:kFALSE;
  17
  18   // Check the analysis type using the event handlers connected to the analysis manager.
  19   //==============================================================================
  20   if (!mgr->GetInputEventHandler())
  21     {
  22       ::Error("AddTaskEmcalJetMass", "This task requires an input event handler");
  23       return NULL;
  24     }
  25
  26   TString wagonName = Form("JetMassResponseDet_%s%s",njetsDet,tag.Data());
  27   TString strType(type);
  28
  29   //Configure jet mass detector response task
  30   AliAnalysisTaskJetMassResponseDet *task = new AliAnalysisTaskJetMassResponseDet(wagonName.Data());
  31
  32   task->SetNCentBins(1);
  33   //task->SetVzRange(-10.,10.);
  34
  35   task->SetJetContainerPart(0);
  36   task->SetJetContainerDet(1);
  37
  38   AliJetContainer *jetContPart = task->AddJetContainer(njetsPart,strType,R);
  39   if(jetContPart) {
  40     // jetContPart->SetPercAreaCut(0.6);
  41   }
  42
  43   AliJetContainer *jetContDet = task->AddJetContainer(njetsDet,strType,R);
  44   if(jetContDet) {
  45     jetContDet->SetPercAreaCut(0.6);
  46   }
  47
  48   task->SetCaloTriggerPatchInfoName(kEmcalTriggers.Data());
  49   task->SelectCollisionCandidates(pSel);
  50   task->SetUseAliAnaUtils(kFALSE);
  51
  52   mgr->AddTask(task);
  53
  54   //Connnect input
  55   mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer() );
  56
  57   //Connect output
  58   TString contName(wagonName);
  59   TString outputfile = Form("%s",AliAnalysisManager::GetCommonFileName());
  60   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(contName.Data(), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
  61   mgr->ConnectOutput(task,1,coutput1);
  62
  63   return task;
  64 }
  65