PWG3/vertexingHF/AddTaskDStarJets.C

   1 //DEFINITION OF A FEW CONSTANTS
   2 const Int_t    mintrackrefsTPC = 0 ;
   3 const Int_t    mintrackrefsITS = 3 ;
   4 const Int_t    PDG = 421;
   5 const Int_t    minclustersTPC = 0 ;
   6 const Int_t    minITSClusters = 4;
   7 // ANALYSIS TYPE D*+ or D*-
   8 const Bool_t computeD0 = kFALSE;
   9 const Bool_t topologicalCut = kFALSE;
  10 //----------------------------------------------------
  11
  12 AliAnalysisTaskSEDStarJets *AddTaskDStarJets(Bool_t theMCon=kTRUE)
  13 {
  14
  15   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  16   if (!mgr) {
  17     ::Error("AddTaskDStarJets", "No analysis manager to connect to.");
  18     return NULL;
  19   }
  20
  21   //CREATE THE TASK
  22   printf("CREATE TASK\n");
  23   // create the task
  24   AliAnalysisTaskSEDStarJets *task = new AliAnalysisTaskSEDStarJets("AliAnalysisTaskSEDStarJets");
  25   task->SetMinITSClusters(minITSClusters);
  26   task->SetAnalType(computeD0);
  27   task->SetMC(theMCon);
  28   task->SetCutType(topologicalCut);
  29   // Create and connect containers for input/output
  30
  31   TString outputfile = AliAnalysisManager::GetCommonFileName();
  32   outputfile += ":PWG3_D2H_DStarJet";
  33
  34   // ------ input data ------
  35   AliAnalysisDataContainer *cinput0  = mgr->GetCommonInputContainer();
  36
  37   // ----- output data -----
  38
  39   // output TH1I for event counting
  40   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist0", TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
  41
  42   mgr->AddTask(task);
  43
  44   mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
  45   mgr->ConnectOutput(task,1,coutput1);
  46
  47   return task ;
  48 }
  49