]>
Commit | Line | Data |
---|---|---|
2da38f7e | 1 | //if like define a different number of signal for TPC PID\r |
2 | //by default the task is anyway computing 1, 2 and 3 sigmas\r | |
d962650e | 3 | const Bool_t theRareOn = kFALSE;\r |
4 | const Bool_t anaType = 1;//0 HD; 1 UU;\r | |
dfb051fd | 5 | //----------------------------------------------------\r |
6 | \r | |
7 | AliAnalysisTaskSEDStarSpectra *AddTaskDStarSpectra(Bool_t theMCon=kFALSE)\r | |
8 | {\r | |
9 | \r | |
10 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r | |
11 | if (!mgr) {\r | |
12 | ::Error("AddTaskDStarSpectra", "No analysis manager to connect to.");\r | |
13 | return NULL;\r | |
14 | } \r | |
15 | \r | |
2da38f7e | 16 | // cuts are stored in a TFile generated by makeTFile4CutsDStartoKpipi.C in ./macros/\r |
17 | // set there the cuts!!!!!\r | |
ae39ad32 | 18 | Bool_t stdcuts=kFALSE;\r |
2da38f7e | 19 | TFile* filecuts=new TFile("DStartoKpipiCuts.root");\r |
20 | if(!filecuts->IsOpen()){\r | |
21 | cout<<"Input file not found: exit"<<endl;\r | |
ae39ad32 | 22 | stdcuts=kTRUE;\r |
dfb051fd | 23 | }\r |
24 | \r | |
2da38f7e | 25 | AliRDHFCutsDStartoKpipi* RDHFDStartoKpipi=new AliRDHFCutsDStartoKpipi();\r |
ae39ad32 | 26 | if(stdcuts) RDHFDStartoKpipi->SetStandardCutsPP2010();\r |
27 | else RDHFDStartoKpipi = (AliRDHFCutsDStartoKpipi*)filecuts->Get("DStartoKpipiCuts");\r | |
2da38f7e | 28 | RDHFDStartoKpipi->SetName("DStartoKpipiCuts");\r |
dfb051fd | 29 | \r |
2da38f7e | 30 | // mm let's see if everything is ok\r |
31 | if(!RDHFDStartoKpipi){\r | |
32 | cout<<"Specific AliRDHFCuts not found"<<endl;\r | |
33 | return;\r | |
dfb051fd | 34 | }\r |
dfb051fd | 35 | \r |
36 | //CREATE THE TASK\r | |
37 | printf("CREATE TASK\n");\r | |
38 | // create the task\r | |
39 | AliAnalysisTaskSEDStarSpectra *task = new AliAnalysisTaskSEDStarSpectra("AliAnalysisTaskSEDStarSpectra",RDHFDStartoKpipi);\r | |
40 | task->SetAnalysisType(anaType);\r | |
dfb051fd | 41 | task->SetMC(theMCon);\r |
d962650e | 42 | task->SetRareSearch(theRareOn);\r |
dfb051fd | 43 | task->SetDebugLevel(0);\r |
44 | \r | |
45 | mgr->AddTask(task);\r | |
46 | \r | |
47 | // Create and connect containers for input/output\r | |
48 | \r | |
49 | TString outputfile = AliAnalysisManager::GetCommonFileName();\r | |
2da38f7e | 50 | outputfile += ":PWG3_D2H_DStarSpectra";\r |
dfb051fd | 51 | \r |
52 | // ------ input data ------\r | |
53 | \r | |
54 | //AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();\r | |
55 | AliAnalysisDataContainer *cinput0 = mgr->CreateContainer("indstar",TChain::Class(), \r | |
56 | AliAnalysisManager::kInputContainer);\r | |
57 | // ----- output data -----\r | |
58 | \r | |
59 | AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1",TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());\r | |
d962650e | 60 | AliAnalysisDataContainer *coutputDStar1 = mgr->CreateContainer("DStarAll",TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());\r |
61 | AliAnalysisDataContainer *coutputDStar2 = mgr->CreateContainer("DStarPID",TList::Class(),AliAnalysisManager::kOutputContainer, outputfile.Data());\r | |
62 | AliAnalysisDataContainer *coutputDStar3 = mgr->CreateContainer("cuts",AliRDHFCutsDStartoKpipi::Class(),AliAnalysisManager::kOutputContainer, outputfile.Data()); //cuts\r | |
63 | AliAnalysisDataContainer *coutputDstarNorm = mgr->CreateContainer("coutputDstarNorm",AliNormalizationCounter::Class(),AliAnalysisManager::kOutputContainer, outputfile.Data());\r | |
dfb051fd | 64 | \r |
65 | mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());\r | |
66 | \r | |
67 | mgr->ConnectOutput(task,1,coutput1);\r | |
68 | mgr->ConnectOutput(task,2,coutputDStar1);\r | |
69 | mgr->ConnectOutput(task,3,coutputDStar2);\r | |
70 | mgr->ConnectOutput(task,4,coutputDStar3);\r | |
d962650e | 71 | mgr->ConnectOutput(task,5,coutputDstarNorm);\r |
dfb051fd | 72 | \r |
73 | return task;\r | |
74 | }\r | |
75 | \r |