New task for D* pt-dep analysis (A. Grelli)
[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AddTaskDStarSpectra.C
1 //DEFINITION OF A FEW CONSTANTS
2 const Int_t    minITSClusters = 5;
3 const Int_t    minITSClustersSoft = 4;
4 const Int_t    numberOfSigmasPID = 3;
5 // ANALYSIS TYPE DATA/MC
6 const Bool_t usePIDforKaons = kFALSE;
7 //----------------------------------------------------
8
9 AliAnalysisTaskSEDStarSpectra *AddTaskDStarSpectra(Bool_t theMCon=kTRUE)
10 {
11
12   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
13   if (!mgr) {
14     ::Error("AddTaskDStarJets", "No analysis manager to connect to.");
15     return NULL;
16   }  
17   
18   //CREATE THE TASK
19   printf("CREATE TASK\n");
20   // create the task
21   AliAnalysisTaskSEDStarSpectra *task = new AliAnalysisTaskSEDStarSpectra("AliAnalysisTaskSEDStarSpectra");
22   task->SetMinITSClusters(minITSClusters);
23   task->SetMinITSClustersSoft(minITSClustersSoft);
24   task->SetMinITSClustersSoft(numberOfSigmasPID);
25   task->SetMC(theMCon);
26   task->SetMC(usePIDforKaons);
27
28   // Create and connect containers for input/output
29   
30   TString outputfile = AliAnalysisManager::GetCommonFileName();
31   outputfile += ":PWG3_D2H_DStarSpectra";
32
33   // ------ input data ------
34   AliAnalysisDataContainer *cinput0  = mgr->GetCommonInputContainer();
35   
36   // ----- output data -----
37   
38   // output TH1I for event counting
39   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist0", TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
40   
41   mgr->AddTask(task);
42   
43   mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
44   mgr->ConnectOutput(task,1,coutput1);
45   
46   return task ;
47 }
48