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