]>
Commit | Line | Data |
---|---|---|
218da8e5 | 1 | AliAnalysisTaskSEDs *AddTaskDs(Bool_t readMC=kTRUE) |
25c94fc8 | 2 | { |
3 | // | |
4 | // Test macro for the AliAnalysisTaskSE for Ds candidates | |
5 | ||
6 | ||
7 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); | |
8 | if (!mgr) { | |
9 | ::Error("AddTaskDs", "No analysis manager to connect to."); | |
10 | return NULL; | |
11 | } | |
12 | ||
13 | ||
4c230f6d | 14 | |
15 | //ptbins | |
16 | const Int_t nptbins=4; | |
17 | Float_t ptbins[nptbins+1]={0.,2.,3.,5.,99999.}; | |
18 | ||
19 | //setting cut values | |
20 | Int_t nvars=14; | |
21 | Float_t** ancutsval; | |
22 | ancutsval=new Float_t*[nvars]; | |
23 | for(Int_t ic=0;ic<nvars;ic++){ancutsval[ic]=new Float_t[nptbins];} | |
24 | for(Int_t ipt=0;ipt<nptbins;ipt++){ | |
25 | ancutsval[0][ipt]=0.2; | |
26 | ancutsval[1][ipt]=0.4; | |
27 | ancutsval[2][ipt]=0.4; | |
28 | ancutsval[3][ipt]=0.; | |
29 | ancutsval[4][ipt]=0.; | |
30 | ancutsval[5][ipt]=0.005; | |
31 | ancutsval[6][ipt]=0.038; | |
32 | ancutsval[7][ipt]=0.; | |
33 | ancutsval[8][ipt]=0.; | |
34 | ancutsval[9][ipt]=0.95; | |
35 | ancutsval[10][ipt]=0.; | |
36 | ancutsval[11][ipt]=0.1; | |
37 | ancutsval[12][ipt]=0.004; | |
38 | ancutsval[13][ipt]=0.035; | |
39 | } | |
40 | Float_t** prcutsval; | |
41 | prcutsval=new Float_t*[nvars]; | |
42 | for(Int_t ic=0;ic<nvars;ic++){prcutsval[ic]=new Float_t[nptbins];} | |
43 | for(Int_t ipt=0;ipt<nptbins;ipt++){ | |
44 | prcutsval[0][ipt]=0.2; | |
45 | prcutsval[1][ipt]=0.4; | |
46 | prcutsval[2][ipt]=0.4; | |
47 | prcutsval[3][ipt]=0.; | |
48 | prcutsval[4][ipt]=0.; | |
49 | prcutsval[5][ipt]=0.005; | |
50 | prcutsval[6][ipt]=0.06; | |
51 | prcutsval[7][ipt]=0.; | |
52 | prcutsval[8][ipt]=0.; | |
53 | prcutsval[9][ipt]=0.85; | |
54 | prcutsval[10][ipt]=0.; | |
55 | prcutsval[11][ipt]=0.1; | |
56 | prcutsval[12][ipt]=0.1; | |
57 | prcutsval[13][ipt]=0.1; | |
58 | } | |
59 | ||
60 | //Analysis cuts | |
61 | AliRDHFCutsDstoKKpi *analysiscuts = new AliRDHFCutsDstoKKpi(); | |
62 | analysiscuts->SetPtBins(nptbins+1,ptbins); | |
63 | analysiscuts->SetCuts(nvars,nptbins,ancutsval); | |
64 | AliRDHFCutsDstoKKpi *prodcuts = new AliRDHFCutsDstoKKpi(); | |
65 | prodcuts->SetPtBins(nptbins+1,ptbins); | |
66 | prodcuts->SetCuts(nvars,nptbins,prcutsval); | |
67 | ||
68 | // Analysis task | |
69 | AliAnalysisTaskSEDs *dsTask = new AliAnalysisTaskSEDs("DsAnalysis",prodcuts,analysiscuts); | |
218da8e5 | 70 | dsTask->SetReadMC(readMC); |
25c94fc8 | 71 | dsTask->SetDebugLevel(0); |
72 | mgr->AddTask(dsTask); | |
73 | ||
74 | // | |
75 | // Create containers for input/output | |
76 | AliAnalysisDataContainer *cinputDs = mgr->CreateContainer("cinputDs",TChain::Class(), | |
77 | AliAnalysisManager::kInputContainer); | |
78 | ||
4c230f6d | 79 | TString outputfile = AliAnalysisManager::GetCommonFileName(); |
80 | outputfile += ":PWG3_D2H_InvMassDs"; | |
81 | ||
82 | AliAnalysisDataContainer *coutputDsCuts = mgr->CreateContainer("coutputDsCuts",TList::Class(), | |
83 | AliAnalysisManager::kOutputContainer, | |
84 | outputfile.Data()); | |
85 | ||
86 | ||
25c94fc8 | 87 | AliAnalysisDataContainer *coutputDs = mgr->CreateContainer("coutputDs",TList::Class(), |
88 | AliAnalysisManager::kOutputContainer, | |
89 | outputfile.Data()); | |
90 | ||
4c230f6d | 91 | |
25c94fc8 | 92 | mgr->ConnectInput(dsTask,0,mgr->GetCommonInputContainer()); |
93 | ||
94 | mgr->ConnectOutput(dsTask,1,coutputDs); | |
4c230f6d | 95 | |
96 | mgr->ConnectOutput(dsTask,2,coutputDsCuts); | |
25c94fc8 | 97 | return dsTask; |
98 | } |