AliPhysics  3b4a69f (3b4a69f)
AliHFSystErr.cxx
Go to the documentation of this file.
1 
2 /**************************************************************************
3  * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
4  * *
5  * Author: The ALICE Off-line Project. *
6  * Contributors are mentioned in the code where appropriate. *
7  * *
8  * Permission to use, copy, modify and distribute this software and its *
9  * documentation strictly for non-commercial purposes is hereby granted *
10  * without fee, provided that the above copyright notice appears in all *
11  * copies and that both the copyright notice and this permission notice *
12  * appear in the supporting documentation. The authors make no claims *
13  * about the suitability of this software for any purpose. It is *
14  * provided "as is" without express or implied warranty. *
15  **************************************************************************/
16 
17 /* $Id$ */
18 
20 //
21 // Class to handle systematic errors for charm hadrons
22 //
23 // Usage:
24 // AliHFSystEff syst; // DECAY = 1 for D0, 2, for D+, 3 for D*, 5 for Lc->pKp, 6 for Lc->pK0S
25 // syst.SetRunNumber(YEAR); // YEAR = two last numbers of the year (is 10 for 2010)
26 // syst.SetCollisionType(TYPE); // TYPE = 0 is pp, 1 is PbPb
27 // syst.SetCentrality(CENT); // CENT is centrality, 0100 for MB, 020 (4080) for 0-20 (40-80) CC...
28 // syst.Init(DECAY); // DECAY = 1 for D0, 2, for D+, 3 for D*
29 // syst.DrawErrors(); // to see a plot of the error contributions
30 // syst.GetTotalSystErr(pt); // to get the total err at pt
31 //
32 // Author: A.Dainese, andrea.dainese@pd.infn.it
34 
35 #include <TStyle.h>
36 #include <TGraphAsymmErrors.h>
37 #include <TMath.h>
38 #include <TCanvas.h>
39 #include <TH2F.h>
40 #include <TLegend.h>
41 #include <TColor.h>
42 
43 #include "AliLog.h"
44 #include "AliHFSystErr.h"
45 
47 ClassImp(AliHFSystErr);
49 
50 //--------------------------------------------------------------------------
51 AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
52  TNamed(name,title),
53  fNorm(0),
54  fRawYield(0),
55  fTrackingEff(0),
56  fBR(0),
57  fCutsEff(0),
58  fPIDEff(0),
59  fMCPtShape(0),
60  fPartAntipart(0),
61  fRunNumber(10),
62  fCollisionType(0),
63  fCentralityClass("0100"),
64  fRapidityRange("0101"),
65  fIsLowEnergy(false),
66  fIsLowPtAnalysis(false),
67  fIsPass4Analysis(false),
68  fIs5TeVAnalysis(false),
69  fIsBDTAnalysis(false),
70  fIsCentScan(false),
71  fStandardBins(false),
72  fIsRapidityScan(false)
73 {
74  //
76  //
77 }
78 
79 //--------------------------------------------------------------------------
80 AliHFSystErr::~AliHFSystErr() {
81  //
83  //
84  /*
85 
86  if(fNorm) { delete fNorm; fNorm=0; }
87  if(fRawYield) { delete fRawYield; fRawYield=0; }
88  if(fTrackingEff) { delete fTrackingEff; fTrackingEff=0; }
89  if(fBR) { delete fBR; fBR=0; }
90  if(fCutsEff) { delete fCutsEff; fCutsEff=0; }
91  if(fPIDEff) { delete fPIDEff; fPIDEff=0; }
92  if(fMCPtShape) { delete fMCPtShape; fMCPtShape=0; }
93  if(fPartAntipart) { delete fPartAntipart; fPartAntipart=0; }
94  */
95 }
96 
97 //--------------------------------------------------------------------------
98 void AliHFSystErr::Init(Int_t decay){
99  //
101  //
102 
103  // if ((fRunNumber>11) && fIsLowEnergy==false) {
104  // AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
105  // }
106 
107  switch(decay) {
108  case 1: // D0->Kpi
109  if (fCollisionType==0) {
110  if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
111  else if(fIs5TeVAnalysis){
112  if(fIsLowPtAnalysis){
113  if(fStandardBins) InitD0toKpi2017pp5TeVLowPtAn();
114  else InitD0toKpi2017pp5TeVLowPtAn_finebins();
115  }else{
116  if(fRunNumber==17){
117  if(fStandardBins)InitD0toKpi2017pp5TeV();
118  else InitD0toKpi2017pp5TeV_finebins();
119  }
120  else InitD0toKpi2015pp5TeV();
121  }
122  }
123  else if(fRunNumber == 10){
124  if(fIsLowPtAnalysis) InitD0toKpi2010ppLowPtAn();
125  else if(fIsPass4Analysis) InitD0toKpi2010ppPass4();
126  else InitD0toKpi2010pp();
127  }
128  else if(fRunNumber == 16 || fRunNumber==2016) InitD0toKpi2016pp13TeV();
129  else AliFatal("Not yet implemented");
130  }
131  else if (fCollisionType==1) {
132  if (fRunNumber == 10){
133  if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
134  else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
135  else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
136  else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
137  else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
138  else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
139  else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
140  else AliFatal("Not yet implemented");
141  }
142  else if (fRunNumber == 11){
143  if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
144  else if (fCentralityClass=="010") InitD0toKpi2011PbPb010();
145  else if (fCentralityClass=="3050InPlane") InitD0toKpi2011PbPb3050InPlane();
146  else if (fCentralityClass=="3050OutOfPlane")InitD0toKpi2011PbPb3050OutOfPlane();
147  else if (fCentralityClass == "3050")InitD0toKpi2011PbPb3050();
148  else if (fCentralityClass=="010" && fIsCentScan) InitD0toKpi2011PbPb010CentScan();
149  else if (fCentralityClass=="1020") InitD0toKpi2011PbPb1020CentScan();
150  else if (fCentralityClass=="2030") InitD0toKpi2011PbPb2030CentScan();
151  else if (fCentralityClass=="3040") InitD0toKpi2011PbPb3040CentScan();
152  else if (fCentralityClass=="4050") InitD0toKpi2011PbPb4050CentScan();
153  else if (fCentralityClass=="5080") InitD0toKpi2010PbPb5080CentScan();
154  else AliFatal("Not yet implemented");
155  }
156  else if (fRunNumber == 15){
157  if (fCentralityClass=="010") InitD0toKpi2015PbPb010();
158  else if (fCentralityClass=="3050") InitD0toKpi2015PbPb3050();
159  else if (fCentralityClass=="6080") InitD0toKpi2015PbPb6080();
160  else AliFatal("Not yet implemented");
161  }
162  }
163  else if (fCollisionType==2) {
164  if (fCentralityClass=="0100"){
165  if(fIsLowPtAnalysis){
166  if(fRunNumber==16 || fRunNumber==2016) InitD0toKpi2016pPb0100LowPtAn();
167  else InitD0toKpi2013pPb0100LowPtAn();
168  }else{
169  if(fRunNumber==16 || fRunNumber==2016) InitD0toKpi2016pPb0100();
170  else InitD0toKpi2013pPb0100();
171  }
172  }
173  if(fRunNumber==16 || fRunNumber==2016){
174  if (fCentralityClass=="010ZNA") InitD0toKpi2016pPb010ZNA();
175  else if(fCentralityClass=="60100ZNA") InitD0toKpi2016pPb60100ZNA();
176  }else{
177  if (fCentralityClass=="020V0A") InitD0toKpi2013pPb020V0A();
178  if (fCentralityClass=="2040V0A") InitD0toKpi2013pPb2040V0A();
179  if (fCentralityClass=="4060V0A") InitD0toKpi2013pPb4060V0A();
180  if (fCentralityClass=="60100V0A") InitD0toKpi2013pPb60100V0A();
181 
182  if (fCentralityClass=="020ZNA") InitD0toKpi2013pPb020ZNA();
183  if (fCentralityClass=="2040ZNA") InitD0toKpi2013pPb2040ZNA();
184  if (fCentralityClass=="4060ZNA") InitD0toKpi2013pPb4060ZNA();
185  if (fCentralityClass=="60100ZNA")InitD0toKpi2013pPb60100ZNA();
186 
187  if (fCentralityClass=="020CL1") InitD0toKpi2013pPb020CL1();
188  if (fCentralityClass=="2040CL1") InitD0toKpi2013pPb2040CL1();
189  if (fCentralityClass=="4060CL1") InitD0toKpi2013pPb4060CL1();
190  if (fCentralityClass=="60100CL1") InitD0toKpi2013pPb60100CL1();
191 
192  if (fIsRapidityScan) {
193  if (fRapidityRange == "0804") InitD0toKpi2013pPb0100RapScan0804();
194  if (fRapidityRange == "0401") InitD0toKpi2013pPb0100RapScan0401();
195  if (fRapidityRange == "0101") InitD0toKpi2013pPb0100RapScan0101();
196  if (fRapidityRange == "0104") InitD0toKpi2013pPb0100RapScan0104();
197  if (fRapidityRange == "0408") InitD0toKpi2013pPb0100RapScan0408();
198  }
199  }
200  }
201  else AliFatal("Not yet implemented");
202  break;
203 
204  case 2: // D+->Kpipi
205  if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
206  if (fCollisionType==0) {
207  if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
208  else if(fIs5TeVAnalysis){
209  if(fRunNumber==17){
210  if(fStandardBins) InitDplustoKpipi2017pp5TeV();
211  else InitDplustoKpipi2017pp5TeV_finebins();
212  }
213  else InitDplustoKpipi2015pp5TeV();
214  }
215  else if(fRunNumber == 10){
216  if(fIsPass4Analysis) InitDplustoKpipi2010ppPass4();
217  else InitDplustoKpipi2010pp();
218  } else if(fRunNumber == 12){
219  InitDplustoKpipi2012pp();
220  } else if(fRunNumber == 16 || fRunNumber == 2016){
221  InitDplustoKpipi2016pp13TeV();
222  } else AliFatal("Not yet implemented");
223  }
224  else if (fCollisionType==1) {
225  if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
226  if (fRunNumber == 10){
227  if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
228  else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
229  else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
230  else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
231  else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
232  else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
233  else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
234  else AliFatal("Not yet implemented");
235  }
236  if(fRunNumber == 11){
237  if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
238  else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010();
239  else if (fCentralityClass=="010" && fIsCentScan) InitDplustoKpipi2011PbPb010CentScan();
240  else if (fCentralityClass=="1020") InitDplustoKpipi2011PbPb1020CentScan();
241  else if (fCentralityClass=="2030") InitDplustoKpipi2011PbPb2030CentScan();
242  else if (fCentralityClass=="3040") InitDplustoKpipi2011PbPb3040CentScan();
243  else if (fCentralityClass=="4050") InitDplustoKpipi2011PbPb4050CentScan();
244  else if (fCentralityClass=="5080") InitDplustoKpipi2010PbPb5080CentScan();
245  else if (fCentralityClass=="3050") InitDplustoKpipi2011PbPb3050();
246  else AliFatal("Not yet implemented");
247  }
248  if(fRunNumber == 15){
249  if (fCentralityClass=="010") InitDplustoKpipi2015PbPb010();
250  else if (fCentralityClass=="3050") InitDplustoKpipi2015PbPb3050();
251  else if (fCentralityClass=="6080") InitDplustoKpipi2015PbPb6080();
252  else AliFatal("Not yet implemented");
253  }
254  }
255  else if (fCollisionType==2) {
256  if(fRunNumber==16 || fRunNumber==2016) {
257  if (fCentralityClass=="0100") InitDplustoKpipi2016pPb0100();
258 
259  if (fCentralityClass=="010ZNA") InitDplustoKpipi2016pPb010ZNA();
260  if (fCentralityClass=="60100ZNA") InitDplustoKpipi2016pPb60100ZNA();
261 
262  if (fCentralityClass=="140trkl") InitDplustoKpipi2016pPb140trkl();
263  if (fCentralityClass=="4070trkl") InitDplustoKpipi2016pPb4070trkl();
264  if (fCentralityClass=="70200trkl") InitDplustoKpipi2016pPb70200trkl();
265  }
266  else {
267  if (fCentralityClass=="0100"){InitDplustoKpipi2013pPb0100();}
268 
269  if (fCentralityClass=="020V0A") InitDplustoKpipi2013pPb020V0A();
270  if (fCentralityClass=="2040V0A") InitDplustoKpipi2013pPb2040V0A();
271  if (fCentralityClass=="4060V0A") InitDplustoKpipi2013pPb4060V0A();
272  if (fCentralityClass=="60100V0A") InitDplustoKpipi2013pPb60100V0A();
273 
274  if (fCentralityClass=="020ZNA") InitDplustoKpipi2013pPb020ZNA();
275  if (fCentralityClass=="2040ZNA") InitDplustoKpipi2013pPb2040ZNA();
276  if (fCentralityClass=="4060ZNA") InitDplustoKpipi2013pPb4060ZNA();
277  if (fCentralityClass=="60100ZNA") InitDplustoKpipi2013pPb60100ZNA();
278 
279 
280  if (fCentralityClass=="020CL1") InitDplustoKpipi2013pPb020CL1();
281  if (fCentralityClass=="2040CL1") InitDplustoKpipi2013pPb2040CL1();
282  if (fCentralityClass=="4060CL1") InitDplustoKpipi2013pPb4060CL1();
283  if (fCentralityClass=="60100CL1") InitDplustoKpipi2013pPb60100CL1();
284 
285  if (fIsRapidityScan) {
286  if (fRapidityRange == "0804") InitDplustoKpipi2013pPb0100RapScan0804();
287  if (fRapidityRange == "0401") InitDplustoKpipi2013pPb0100RapScan0401();
288  if (fRapidityRange == "0101") InitDplustoKpipi2013pPb0100RapScan0101();
289  if (fRapidityRange == "0104") InitDplustoKpipi2013pPb0100RapScan0104();
290  if (fRapidityRange == "0408") InitDplustoKpipi2013pPb0100RapScan0408();
291  }
292  }
293  }
294  else AliFatal("Not yet implemented");
295  break;
296  case 3: // D*->D0pi
297  if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
298  if (fCollisionType==0) {
299  if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
300  else if(fRunNumber == 10 || fRunNumber==2010){
301  if(fIsPass4Analysis) InitDstartoD0pi2010ppPass4();
302  else InitDstartoD0pi2010pp();
303  } else if(fRunNumber == 12 || fRunNumber==2012){
304  InitDstartoD0pi2012pp();
305  } else if(fRunNumber == 16 || fRunNumber==2016){
306  InitDstartoKpipi2016pp13TeV();
307  } else if(fRunNumber == 17 || fRunNumber == 2017){
308  if(fIs5TeVAnalysis){
309  if(fStandardBins) InitDstartoD0pi2017pp5TeV();
310  else InitDstartoD0pi2017pp5TeV_finebins();
311  }
312  } else AliFatal("Not yet implemented");
313  }
314  else if (fCollisionType==1) {
315  if (fRunNumber == 10 || fRunNumber==2010){
316  if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
317  else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
318  else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
319  else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
320  else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
321  else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
322  else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
323  else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
324  else AliFatal("Not yet implemented");
325  }
326  if (fRunNumber == 11 || fRunNumber==2011){
327  if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
328  else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010();
329  else if (fCentralityClass=="010" && fIsCentScan) InitDstartoD0pi2011PbPb010CentScan();
330  else if (fCentralityClass=="1020") InitDstartoD0pi2011PbPb1020CentScan();
331  else if (fCentralityClass=="2030") InitDstartoD0pi2011PbPb2030CentScan();
332  else if (fCentralityClass=="3040") InitDstartoD0pi2011PbPb3040CentScan();
333  else if (fCentralityClass=="4050") InitDstartoD0pi2011PbPb4050CentScan();
334  else if (fCentralityClass=="5080") InitDstartoD0pi2010PbPb5080CentScan();
335  else if (fCentralityClass=="3050") InitDstartoD0pi2011PbPb3050();
336  else AliFatal("Not yet implemented");
337  }
338  else if (fRunNumber == 15 || fRunNumber==2015){
339  if (fCentralityClass=="010") InitDstartoD0pi2015PbPb010();
340  else if (fCentralityClass=="3050") InitDstartoD0pi2015PbPb3050();
341  else if (fCentralityClass=="6080") InitDstartoD0pi2015PbPb6080();
342  else AliFatal("Not yet implemented");
343  }
344  }
345  else if (fCollisionType==2) {
346  if (fRunNumber == 16 || fRunNumber==2016){
347  if (fCentralityClass=="0100")InitDstartoD0pi2016pPb0100();
348  else if (fCentralityClass=="010ZNA")InitDstartoD0pi2016pPb010ZNA();
349  else if (fCentralityClass=="60100ZNA")InitDstartoD0pi2016pPb60100ZNA();
350  else AliFatal("Not yet implemented");
351 
352 
353  }
354  else if (fRunNumber == 13 || fRunNumber==2013){
355 
356  if (fCentralityClass=="020V0A") InitDstartoD0pi2013pPb020V0A();
357  if (fCentralityClass=="2040V0A") InitDstartoD0pi2013pPb2040V0A();
358  if (fCentralityClass=="4060V0A") InitDstartoD0pi2013pPb4060V0A();
359  if (fCentralityClass=="60100V0A") InitDstartoD0pi2013pPb60100V0A();
360 
361  if (fCentralityClass=="020ZNA") InitDstartoD0pi2013pPb020ZNA();
362  if (fCentralityClass=="2040ZNA") InitDstartoD0pi2013pPb2040ZNA();
363  if (fCentralityClass=="4060ZNA") InitDstartoD0pi2013pPb4060ZNA();
364  if (fCentralityClass=="60100ZNA") InitDstartoD0pi2013pPb60100ZNA();
365 
366  if (fCentralityClass=="020CL1") InitDstartoD0pi2013pPb020CL1();
367  if (fCentralityClass=="2040CL1") InitDstartoD0pi2013pPb2040CL1();
368  if (fCentralityClass=="4060CL1") InitDstartoD0pi2013pPb4060CL1();
369  if (fCentralityClass=="60100CL1") InitDstartoD0pi2013pPb60100CL1();
370 
371  if (fIsRapidityScan) {
372  if (fRapidityRange == "0804") InitDstartoD0pi2013pPb0100RapScan0804();
373  if (fRapidityRange == "0401") InitDstartoD0pi2013pPb0100RapScan0401();
374  if (fRapidityRange == "0101") InitDstartoD0pi2013pPb0100RapScan0101();
375  if (fRapidityRange == "0104") InitDstartoD0pi2013pPb0100RapScan0104();
376  if (fRapidityRange == "0408") InitDstartoD0pi2013pPb0100RapScan0408();
377  }
378  }
379  }
380  else AliFatal("Not yet implemented");
381  break;
382  case 4: // D+s->KKpi
383  if(fIsLowPtAnalysis) AliFatal("Not yet implemented");
384  if (fCollisionType==0) {
385  if(fIsPass4Analysis) InitDstoKKpi2010ppPass4();
386  else if (fRunNumber==16 || fRunNumber==2016) InitDstoKKpi2016pp13TeV();
387  else if (fRunNumber==17 || fRunNumber==2017){
388  if(fIsBDTAnalysis)
389  InitDstoKKpi2017pp5TeVBDT();
390  else
391  InitDstoKKpi2017pp5TeV();
392  }
393  else InitDstoKKpi2010pp();
394  }
395  else if (fCollisionType==1) {
396  if (fRunNumber == 15){
397  if (fCentralityClass=="010") InitDstoKKpi2015PbPb010();
398  else if (fCentralityClass=="3050") InitDstoKKpi2015PbPb3050();
399  else if (fCentralityClass=="6080") InitDstoKKpi2015PbPb6080();
400  else AliFatal("Not yet implemented");
401  }
402  else{
403  if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
404  else if (fCentralityClass=="010") InitDstoKKpi2011PbPb010();
405  else if (fCentralityClass=="2050") InitDstoKKpi2011PbPb2050();
406  else AliFatal("Not yet implemented");
407  }
408  }
409  else if (fCollisionType==2) {
410  if(fRunNumber==13 || fRunNumber==2013) {
411  if (fCentralityClass=="0100") InitDstoKKpi2013pPb0100();
412  }
413  if(fRunNumber==16 || fRunNumber==2016) {
414  if (fCentralityClass=="0100") InitDstoKKpi2016pPb0100();
415  else if (fCentralityClass=="140trkl") InitDstoKKpi2016pPb140trkl();
416  else if (fCentralityClass=="4070trkl") InitDstoKKpi2016pPb4070trkl();
417  else if (fCentralityClass=="70200trkl") InitDstoKKpi2016pPb70200trkl();
418  }
419  }
420  else AliFatal("Not yet implemented");
421  break;
422 
423  case 5: // Lc->pKpi
424  if (fCollisionType==0) {
425  if (fIsBDTAnalysis) InitLctopKpi2010ppBDT();
426  else InitLctopKpi2010pp();
427  }
428  else if (fCollisionType==2) {
429  if(fRunNumber==13 || fRunNumber==2013) {
430  if (fIsBDTAnalysis) InitLctopKpi2013pPbBDT();
431  else InitLctopKpi2013pPb();
432  }
433  else if(fRunNumber==16 || fRunNumber==2016) {
434  InitLctopKpi2016pPb();
435  }
436  }
437  else AliFatal("Not yet implemented");
438  break;
439  case 6: // Lc->pK0S
440  if (fCollisionType==0) InitLctopK0S2010pp();
441  else if (fCollisionType==2) {
442  if(fRunNumber==13 || fRunNumber==2013) {
443  if (fIsBDTAnalysis) InitLctopK0S2013pPbBDT();
444  else InitLctopK0S2013pPb();
445  }
446  else if(fRunNumber==16 || fRunNumber==2016) {
447  if (fIsBDTAnalysis) InitLctopK0S2016pPbBDT();
448  else InitLctopK0S2016pPb();
449  }
450  }
451  else AliFatal("Not yet implemented");
452  break;
453 
454  default:
455  printf("Invalid decay type: %d\n",decay);
456  break;
457  }
458 
459 }
460 
461 //--------------------------------------------------------------------------
462 void AliHFSystErr::InitD0toKpi2010pp() {
463  //
464  // D0->Kpi syst errors. Responsible: A. Rossi
465  // 2010 pp sample
466  //
467 
468  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010pp");
469 
470  // Normalization
471  fNorm = new TH1F("fNorm","fNorm",24,0,24);
472  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
473 
474  // Branching ratio
475  fBR = new TH1F("fBR","fBR",24,0,24);
476  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
477 
478  // Tracking efficiency
479  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
480  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
481 
482  // Raw yield extraction
483  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
484  fRawYield->SetBinContent(1,1);
485  fRawYield->SetBinContent(2,0.22);
486  fRawYield->SetBinContent(3,0.1);
487  for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
488  for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
489  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
490  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
491 
492  // Cuts efficiency (from cuts variation)
493  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
494  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
495 
496  // PID efficiency (from PID/noPID)
497  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
498  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
499  fPIDEff->SetBinContent(2,0.05); // 5%
500 
501  // MC dN/dpt
502  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
503  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
504  fMCPtShape->SetBinContent(1,0.03);
505  fMCPtShape->SetBinContent(2,0.03);
506 
507  // particle-antiparticle
508  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
509  // fPartAntipart->SetBinContent(1,1);
510  // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
511 
512  return;
513 }
514 //-------------------------------------------------------------------------
515 void AliHFSystErr::InitD0toKpi2010ppPass4(){
516  //
517  // D0->Kpi syst errors. Responsible: C. Terrevoli
518  // 2010 pp sample - pass4
519  //
520  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010ppPass4");
521 
522  // Normalization
523  fNorm = new TH1F("fNorm","fNorm",36,0,36);
524  fNorm->SetBinContent(1,0.0); // bin 0.-1.
525  for(Int_t i=2;i<=36;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
526 
527  // Branching ratio
528  fBR = new TH1F("fBR","fBR",36,0,36);
529  fBR->SetBinContent(1,0.0); // 1.2% PDG2010
530  for(Int_t i=2;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
531 
532  // Tracking efficiency
533  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
534  fTrackingEff->SetBinContent(1,0.0); //
535  fTrackingEff->SetBinContent(2,0.03); //
536  fTrackingEff->SetBinContent(3,0.035); //
537  for(Int_t i=4;i<=5;i++) fTrackingEff->SetBinContent(i,0.04); //
538  for(Int_t i=6;i<=9;i++) fTrackingEff->SetBinContent(i,0.045); //
539  for(Int_t i=10;i<=14;i++) fTrackingEff->SetBinContent(i,0.04); //
540  for(Int_t i=15;i<=36;i++) fTrackingEff->SetBinContent(i,0.035); //
541 
542  // Raw yield extraction
543  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
544  fRawYield->SetBinContent(1,0);
545  fRawYield->SetBinContent(2,0.04);
546  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.03);
547  for(Int_t i=13;i<=36;i++) fRawYield->SetBinContent(i,0.05);
548 
549  // Cuts efficiency (from cuts variation)
550  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
551  fCutsEff->SetBinContent(1,0.0);
552  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.05);
553 
554  // PID efficiency (from PID/noPID)
555  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
556  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
557 
558  // MC dN/dpt
559  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
560  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
561  fMCPtShape->SetBinContent(2,0.03);
562 
563  // particle-antiparticle
564  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
565  // fPartAntipart->SetBinContent(1,1);
566  // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
567 
568  return;
569 }
570 //-------------------------------------------------------------------------
571 void AliHFSystErr::InitD0toKpi2015pp5TeV(){
572  //
573  // D0->Kpi syst errors. Responsible: C. Terrevoli
574  // 2015 pp sample - 5TeV
575  //
576  SetNameTitle("AliHFSystErr","SystErrD0toKpi2015pp5TeV");
577 
578  // Normalization
579  fNorm = new TH1F("fNorm","fNorm",36,0,36);
580  fNorm->SetBinContent(1,0.0); // bin 0.-1.
581  for(Int_t i=2;i<=36;i++) fNorm->SetBinContent(i,0.023);
582 
583  // Branching ratio
584  fBR = new TH1F("fBR","fBR",36,0,36);
585  fBR->SetBinContent(1,0.0); // 1.2% PDG2010
586  for(Int_t i=2;i<=36;i++) fBR->SetBinContent(i,0.01); //
587 
588  // Tracking efficiency
589  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
590  fTrackingEff->SetBinContent(1,0.0); // 0-1
591  fTrackingEff->SetBinContent(2,0.045); // 1-2
592  fTrackingEff->SetBinContent(3,0.04); // 2-3
593  fTrackingEff->SetBinContent(4,0.035); // 3-4
594  for(Int_t i=5;i<=10;i++) fTrackingEff->SetBinContent(i,0.03); //
595  for(Int_t i=11;i<=12;i++) fTrackingEff->SetBinContent(i,0.035); //
596  for(Int_t i=13;i<=16;i++) fTrackingEff->SetBinContent(i,0.045); //
597  for(Int_t i=17;i<=36;i++) fTrackingEff->SetBinContent(i,0.05); //
598 
599  // Raw yield extraction
600  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
601  fRawYield->SetBinContent(1,0.03);
602  for(Int_t i=2;i<=7;i++) fRawYield->SetBinContent(i,0.03);
603  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.04);
604  for(Int_t i=12;i<=16;i++) fRawYield->SetBinContent(i,0.05);
605  for(Int_t i=16;i<=36;i++) fRawYield->SetBinContent(i,0.06);
606 
607  // Cuts efficiency (from cuts variation)
608  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
609  fCutsEff->SetBinContent(1,0.0);
610  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.05);
611 
612  // PID efficiency (from PID/noPID)
613  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
614  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
615 
616  // MC dN/dpt
617  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
618  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
619  fMCPtShape->SetBinContent(2,0.03);
620 
621  // particle-antiparticle
622  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
623  // fPartAntipart->SetBinContent(1,1);
624  // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
625 
626  return;
627 }
628 
629 //--------------------------------------------------------------------------
630 void AliHFSystErr::InitD0toKpi2010PbPb020() {
631  //
632  // D0->Kpi syst errors. Responsible: A. Rossi
633  // 2010 PbPb sample, 0-20 CC
634  //
635 
636  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb020");
637 
638  // Normalization
639  fNorm = new TH1F("fNorm","fNorm",20,0,20);
640  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
641 
642  // Branching ratio
643  fBR = new TH1F("fBR","fBR",20,0,20);
644  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
645 
646  // Tracking efficiency
647  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
648  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
649 
650  // Raw yield extraction
651  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
652  fRawYield->SetBinContent(1,0);
653  fRawYield->SetBinContent(2,0);
654  fRawYield->SetBinContent(3,0.08);
655  for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
656  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
657 
658  // Cuts efficiency (from cuts variation)
659  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
660  fCutsEff->SetBinContent(1,0.);
661  fCutsEff->SetBinContent(2,0.);
662  fCutsEff->SetBinContent(3,0.13);
663  fCutsEff->SetBinContent(4,0.11);
664  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
665  for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
666 
667  // PID efficiency (from PID/noPID)
668  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
669  for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
670 
671  // MC dN/dpt
672  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
673  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
674  fMCPtShape->SetBinContent(3,0.04);
675  fMCPtShape->SetBinContent(4,0.02);
676  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
677 
678  // // particle-antiparticle
679  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
680  // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
681  // fPartAntipart->SetBinContent(3,0.10);
682  // fPartAntipart->SetBinContent(4,0.10);
683  // fPartAntipart->SetBinContent(7,0.10);
684  // fPartAntipart->SetBinContent(8,0.10);
685 
686  return;
687 }
688 //--------------------------------------------------------------------------
689 //--------------------------------------------------------------------------
690 void AliHFSystErr::InitD0toKpi2011PbPb07half() {
691  //
692  // D0->Kpi syst errors. Responsible: A. Rossi
693  // 2011 PbPb sample, 0-7.5 CC
694  //
695 
696  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb07half");
697 
698  // Normalization
699  fNorm = new TH1F("fNorm","fNorm",36,0,36);
700  for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
701 
702  // Branching ratio
703  fBR = new TH1F("fBR","fBR",36,0,36);
704  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
705 
706  // Tracking efficiency
707  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
708  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
709  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
710 
711  // Raw yield extraction
712  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
713  fRawYield->SetBinContent(1,0);
714  fRawYield->SetBinContent(2,0.1);
715  fRawYield->SetBinContent(3,0.05);
716  for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
717  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
718  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
719  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
720 
721  // Cuts efficiency (from cuts variation)
722  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
723  fCutsEff->SetBinContent(1,0.);
724  fCutsEff->SetBinContent(2,0.15);
725  fCutsEff->SetBinContent(3,0.13);
726  fCutsEff->SetBinContent(4,0.11);
727  fCutsEff->SetBinContent(5,0.08);
728  for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
729  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
730 
731  // PID efficiency (from PID/noPID)
732  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
733  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
734  for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
735 
736  // MC dN/dpt
737  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
738  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
739  fMCPtShape->SetBinContent(2,0.06);
740  fMCPtShape->SetBinContent(3,0.04);
741  fMCPtShape->SetBinContent(4,0.02);
742  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
743  for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
744  for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
745 
746  // // particle-antiparticle
747  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
748  // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
749  // fPartAntipart->SetBinContent(3,0.10);
750  // fPartAntipart->SetBinContent(4,0.10);
751  // fPartAntipart->SetBinContent(7,0.10);
752  // fPartAntipart->SetBinContent(8,0.10);
753 
754  return;
755 }
756 
757 //--------------------------------------------------------------------------
758 void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
759  //
760  // D0->Kpi syst errors. Responsible: D. Caffarri
761  // 2011 PbPb sample, 30-50 CC InPlane
762  //
763  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb3050InPlane");
764 
765  InitD0toKpi2011PbPb07half();
766  // Raw yield extraction
767  // fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
768  fRawYield->SetBinContent(1,0);
769  fRawYield->SetBinContent(2,0);
770  fRawYield->SetBinContent(3,0.05);
771  fRawYield->SetBinContent(4,0.08);
772  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
773  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
774  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
775  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
776 
777  // Cuts efficiency (from cuts variation)
778  //fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
779  fCutsEff->SetBinContent(1,0.);
780  fCutsEff->SetBinContent(2,0.0);
781  fCutsEff->SetBinContent(3,0.10);
782  fCutsEff->SetBinContent(4,0.10);
783  fCutsEff->SetBinContent(5,0.10);
784  fCutsEff->SetBinContent(6,0.10);
785  for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
786  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
787  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
788 }
789 
790 //--------------------------------------------------------------------------
791 void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
792  //
793  // D0->Kpi syst errors. Responsible: D. Caffarri
794  // 2011 PbPb sample, 30-50 CC OutOfPlane
795  //
796  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb3050OutOfPlane");
797  InitD0toKpi2011PbPb07half();
798  // Raw yield extraction
799  //fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
800  fRawYield->SetBinContent(1,0);
801  fRawYield->SetBinContent(2,0.);
802  fRawYield->SetBinContent(3,0.05);
803  for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.07);
804  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
805  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
806  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
807  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
808 
809  // Cuts efficiency (from cuts variation)
810  // fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
811  fCutsEff->SetBinContent(1,0.);
812  fCutsEff->SetBinContent(2,0.);
813  fCutsEff->SetBinContent(3,0.10);
814  fCutsEff->SetBinContent(4,0.10);
815  fCutsEff->SetBinContent(5,0.10);
816  fCutsEff->SetBinContent(6,0.10);
817  for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
818  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
819  for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
820 }
821 
822 //--------------------------------------------------------------------------
823 void AliHFSystErr::InitD0toKpi2010PbPb4080() {
824  //
825  // D0->Kpi syst errors. Responsible: A. Rossi
826  // 2010 PbPb sample, 40-80 CC
827  //
828  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010PbPb4080");
829 
830  // Normalization
831  fNorm = new TH1F("fNorm","fNorm",20,0,20);
832  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
833 
834  // Branching ratio
835  fBR = new TH1F("fBR","fBR",20,0,20);
836  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
837 
838  // Tracking efficiency
839  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
840  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
841 
842 
843  // Raw yield extraction
844  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
845  fRawYield->SetBinContent(1,0);
846  fRawYield->SetBinContent(2,0);
847  for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
848  //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
849 
850  // Cuts efficiency (from cuts variation)
851  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
852  fCutsEff->SetBinContent(1,0.);
853  fCutsEff->SetBinContent(2,0.);
854  fCutsEff->SetBinContent(3,0.13);
855  fCutsEff->SetBinContent(4,0.11);
856  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
857  for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
858 
859  // PID efficiency (from PID/noPID)
860  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
861  // for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
862  // for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
863  for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
864 
865  // MC dN/dpt
866  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
867  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
868  // fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
869  // fMCPtShape->SetBinContent(4,0.02);
870  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
871 
872  // // particle-antiparticle
873  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
874  // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
875 
876  return;
877 }
878 
879 //--------------------------------------------------------------------------
880 void AliHFSystErr::InitD0toKpi2011PbPb010() {
881  //
882  // D0->Kpi syst errors. Responsible: A. Rossi
883  // 2011 PbPb sample, 0-10 CC
884  //
885 
886  AliInfo(" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
887  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb010");
888 
889  // Normalization
890  fNorm = new TH1F("fNorm","fNorm",24,0,24);
891  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
892 
893  // Branching ratio
894  fBR = new TH1F("fBR","fBR",24,0,24);
895  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
896 
897  // Tracking efficiency
898  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
899  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
900 
901 
902  // Raw yield extraction
903  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
904  fRawYield->SetBinContent(1,0);
905  fRawYield->SetBinContent(2,0.15);
906  fRawYield->SetBinContent(3,0.09);
907  fRawYield->SetBinContent(4,0.06);
908  fRawYield->SetBinContent(5,0.06);
909  fRawYield->SetBinContent(6,0.06);
910  for (Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.06);
911  for (Int_t i=9; i<=12; i++)fRawYield->SetBinContent(i,0.07);
912  for (Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.11);
913  for (Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.30);
914 
915  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
916  fCutsEff->SetBinContent(1,0.);
917  fCutsEff->SetBinContent(2,0.15);
918  fCutsEff->SetBinContent(3,0.12);
919  fCutsEff->SetBinContent(4,0.05);
920  fCutsEff->SetBinContent(5,0.05);
921  fCutsEff->SetBinContent(6,0.05);
922  for (Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i,0.05);
923  for (Int_t i=9; i<=12; i++)fCutsEff->SetBinContent(i,0.05);
924  for (Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i,0.05);
925  for (Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i,0.05);
926 
927 
928  // PID efficiency (from PID/noPID)
929  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
930  // for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
931  // for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
932  for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
933 
934  // MC dN/dpt
935  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
936  fMCPtShape->SetBinContent(1,0.);
937  fMCPtShape->SetBinContent(2,0.15);
938  fMCPtShape->SetBinContent(3,0.12);
939  fMCPtShape->SetBinContent(4,0.02);
940  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
941 
942  // // particle-antiparticle
943  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
944  // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
945 
946  return;
947 }
948 
949 //--------------------------------------------------------------------------
950 void AliHFSystErr::InitD0toKpi2011PbPb3050() {
951  //
952  // D0->Kpi syst errors. Responsible: D. Caffarri
953  // 2011 PbPb sample, 30-50 CC
954  //
955 
956  SetNameTitle("AliHFSystErr","SystErrD0toKpi2011PbPb3050");
957  // Normalization
958  fNorm = new TH1F("fNorm","fNorm",20,0,20);
959  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
960 
961  // Branching ratio
962  fBR = new TH1F("fBR","fBR",20,0,20);
963  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
964 
965  // Tracking efficiency
966  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
967  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
968 
969 
970  // Raw yield extraction
971  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
972  fRawYield->SetBinContent(1,0);
973  fRawYield->SetBinContent(2,0.102);
974  fRawYield->SetBinContent(3,0.071);
975  fRawYield->SetBinContent(4,0.052);
976  fRawYield->SetBinContent(5,0.063);
977  fRawYield->SetBinContent(6,0.071);
978  fRawYield->SetBinContent(7,0.071);
979  fRawYield->SetBinContent(8,0.071);
980  fRawYield->SetBinContent(9,0.071);
981  fRawYield->SetBinContent(10,0.071);
982  fRawYield->SetBinContent(11,0.071);
983  fRawYield->SetBinContent(12,0.071);
984 
985  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.112);
986  for(Int_t i=17;i<=20;i++) fRawYield->SetBinContent(i,0);
987 
988  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
989  fCutsEff->SetBinContent(1,0.);
990  fCutsEff->SetBinContent(2,0.10);
991  fCutsEff->SetBinContent(3,0.10);
992  fCutsEff->SetBinContent(4,0.10);
993  fCutsEff->SetBinContent(5,0.10);
994  fCutsEff->SetBinContent(6,0.10);
995  for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
996  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
997 
998 
999  // PID efficiency (from PID/noPID)
1000  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1001  // for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
1002  // for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
1003  for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
1004 
1005  // MC dN/dpt
1006  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",0,0,20);
1007  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1008  fMCPtShape->SetBinContent(2,0.06);
1009  fMCPtShape->SetBinContent(3,0.04);
1010  fMCPtShape->SetBinContent(4,0.02);
1011  for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
1012  for(Int_t i=17;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
1013  // // particle-antiparticle
1014  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1015  // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
1016 
1017  return;
1018 }
1019 
1020 
1021 //--------------------------------------------------------------------------
1022 void AliHFSystErr::InitD0toKpi2010ppLowEn() {
1023  //
1024  // D0->Kpi syst errors. Low energy run
1025  // 2011 2.76 TeV pp sample
1026  //
1027  AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
1028  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010ppLowEn");
1029 
1030  // Normalization
1031  fNorm = new TH1F("fNorm","fNorm",20,0,20);
1032  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); // error on sigmaV0and
1033 
1034  // Branching ratio
1035  fBR = new TH1F("fBR","fBR",20,0,20);
1036  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
1037 
1038  // Tracking efficiency
1039  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1040  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
1041 
1042  // Raw yield extraction
1043  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1044  fRawYield->SetBinContent(1,1);
1045  for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
1046 
1047  // Cuts efficiency (from cuts variation)
1048  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1049  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1050  fCutsEff->SetBinContent(2,0.20);
1051  for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
1052 
1053 
1054  // PID efficiency (from PID/noPID)
1055  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1056  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
1057  // fPIDEff->SetBinContent(2,0.20);
1058  for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
1059 
1060  // MC dN/dpt
1061  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1062  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1063  fMCPtShape->SetBinContent(1,0.03);
1064  fMCPtShape->SetBinContent(2,0.03);
1065 
1066  // // particle-antiparticle
1067  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1068  // fPartAntipart->SetBinContent(1,1);
1069  // fPartAntipart->SetBinContent(2,1);
1070  // for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
1071  // for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
1072 
1073  return;
1074 }
1075 //--------------------------------------------------------------------------
1076 void AliHFSystErr::InitD0toKpi2010ppLowPtAn() {
1077  //
1078  // D0->Kpi syst errors. Responsible:
1079  // 2010 pp sample, analysis without topological cuts
1080  //
1081 
1082  AliInfo(" Settings for D0 --> K pi, pp collisions at 7 TeV, analysis without topological cuts");
1083  SetNameTitle("AliHFSystErr","SystErrD0toKpi2010ppLowPtAn");
1084 
1085  // Normalization
1086  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1087  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
1088 
1089  // Branching ratio
1090  fBR = new TH1F("fBR","fBR",24,0,24);
1091  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88)
1092 
1093  // Tracking efficiency
1094  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1095  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
1096 
1097  // Raw yield extraction
1098  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1099  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
1100  fRawYield->SetBinContent(1,0.14);
1101  fRawYield->SetBinContent(2,0.14);
1102  fRawYield->SetBinContent(3,0.10);
1103  fRawYield->SetBinContent(4,0.10);
1104  fRawYield->SetBinContent(5,0.10);
1105  fRawYield->SetBinContent(6,0.14);
1106  fRawYield->SetBinContent(7,0.18);
1107  fRawYield->SetBinContent(8,0.18);
1108  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.2);
1109 
1110  // Cuts efficiency (from cuts variation)
1111  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1112  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.);
1113 
1114  // PID efficiency (from PID/noPID)
1115  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1116  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
1117  fPIDEff->SetBinContent(1,0.05); // 5%
1118  fPIDEff->SetBinContent(2,0.05); // 5%
1119 
1120  // MC dN/dpt
1121  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1122  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1123 
1124  // particle-antiparticle
1125  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
1126  // for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);
1127 
1128  return;
1129 }
1130 //-------------------------------------------------------------------------
1131 void AliHFSystErr::InitD0toKpi2017pp5TeV(){
1132  // TO BE DONE
1133  // D0->Kpi syst errors. Responsible: C. Terrevoli
1134  // 2017 pp sample - 5TeV
1135  //
1136  printf("D0 in pp@5TeV 2017 \n");
1137  SetNameTitle("AliHFSystErr","SystErrD0toKpi2017pp5TeV");
1138 
1139  // Normalization
1140  fNorm = new TH1F("fNorm","fNorm",50,0,50);
1141  // fNorm->SetBinContent(1,0.0); // bin 0.-1.
1142  for(Int_t i=1;i<=50;i++) fNorm->SetBinContent(i,0.021);
1143 
1144  // Branching ratio
1145  fBR = new TH1F("fBR","fBR",50,0,50);
1146  //fBR->SetBinContent(1,0.0); // bin 0.-1.
1147  for(Int_t i=1;i<=50;i++) fBR->SetBinContent(i,0.0129); //
1148 
1149  // Tracking efficiency
1150  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",50,0,50);
1151  fTrackingEff->SetBinContent(1,0.03); // 0-1
1152  fTrackingEff->SetBinContent(2,0.035); // 1-2
1153  fTrackingEff->SetBinContent(3,0.04); // 2-3
1154  for(Int_t i=4;i<=6;i++) fTrackingEff->SetBinContent(i,0.045);
1155  for(Int_t i=7;i<=50;i++) fTrackingEff->SetBinContent(i,0.05);
1156 
1157  // Raw yield extraction
1158  fRawYield = new TH1F("fRawYield","fRawYield",50,0,50);
1159  fRawYield->SetBinContent(1,0.10);
1160  for(Int_t i=2;i<4;i++) fRawYield->SetBinContent(i,0.03);
1161  for(Int_t i=4;i<=13;i++) fRawYield->SetBinContent(i,0.02);
1162  for(Int_t i=14;i<=50;i++) fRawYield->SetBinContent(i,0.05);
1163 
1164  // Cuts efficiency (from cuts variation)
1165  fCutsEff = new TH1F("fCutsEff","fCutsEff",50,0,50);
1166  // fCutsEff->SetBinContent(1,0.0);
1167  for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i,0.05);
1168  for(Int_t i=4;i<=8;i++) fCutsEff->SetBinContent(i,0.03);
1169  for(Int_t i=9;i<=50;i++) fCutsEff->SetBinContent(i,0.03);
1170 
1171  // PID efficiency (from PID/noPID)
1172  fPIDEff = new TH1F("fPIDEff","fPIDEff",50,0,50);
1173  for(Int_t i=1;i<=50;i++) fPIDEff->SetBinContent(i,0.0);
1174 
1175  // MC dN/dpt
1176  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",50,0,50);
1177  fMCPtShape->SetBinContent(1,0.02);
1178  fMCPtShape->SetBinContent(2,0.02);
1179  for(Int_t i=3;i<=50;i++) fMCPtShape->SetBinContent(i,0);
1180 
1181  // particle-antiparticle
1182  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
1183  // fPartAntipart->SetBinContent(1,1);
1184  // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
1185 
1186  return;
1187 }
1188 //-------------------------------------------------------------------------
1189 void AliHFSystErr::InitD0toKpi2017pp5TeV_finebins(){
1190  Float_t xbins[24]={0.,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5,7,7.5,8,9,10,12,16,24,36,50};
1191  // D0->Kpi syst errors. Responsible: C. Terrevoli
1192  // 2017 pp sample - 5TeV with finer pt bins
1193  //
1194  SetNameTitle("AliHFSystErr","SystErrD0toKpi2017pp5TeVFineBins");
1195 
1196  // Normalization
1197  fNorm = new TH1F("fNorm","fNorm",23,xbins);
1198  // fNorm->SetBinContent(1,0.0); // bin 0.-1.
1199  for(Int_t i=1;i<=23;i++) fNorm->SetBinContent(i,0.021);
1200 
1201  // Branching ratio
1202  fBR = new TH1F("fBR","fBR",23,xbins);
1203  //fBR->SetBinContent(1,0.0); // bin 0.-1.
1204  for(Int_t i=1;i<=23;i++) fBR->SetBinContent(i,0.0129); //
1205 
1206  // Tracking efficiency
1207  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",23,xbins);
1208  fTrackingEff->SetBinContent(1,0.03); // 0-0.5
1209  fTrackingEff->SetBinContent(2,0.03); // 0.5-1.
1210  fTrackingEff->SetBinContent(3,0.035); // 1-1.5
1211  fTrackingEff->SetBinContent(4,0.035); // 1.5-2
1212  fTrackingEff->SetBinContent(5,0.04); // 2-2.5
1213  fTrackingEff->SetBinContent(6,0.04); // 2.5-3.
1214  for(Int_t i=7;i<=12;i++) fTrackingEff->SetBinContent(i,0.045);//3-3.5 3.5-4 4-4.5 4.5-5 5-5.5 5.5-6
1215  for(Int_t i=13;i<=23;i++) fTrackingEff->SetBinContent(i,0.05);// >6
1216 
1217  // Raw yield extraction
1218 
1219  fRawYield = new TH1F("fRawYield","fRawYield",23,xbins);//
1220  fRawYield->SetBinContent(1,0.1);//pt bin not used for this analysis 0-0.5, bins
1221  fRawYield->SetBinContent(2,0.1);//pt bin not used for this analysis 0.5-1 bins
1222  fRawYield->SetBinContent(3,0.07);// 1-1.5
1223  for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.03);//(1.5-2 2-2.5 2.5-3)
1224  for(Int_t i=7;i<=19;i++) fRawYield->SetBinContent(i,0.02);//3-3.5 --> 12
1225  for(Int_t i=20;i<=23;i++) fRawYield->SetBinContent(i,0.05);//>12
1226 
1227  // Cuts efficiency (from cuts variation)
1228  fCutsEff = new TH1F("fCutsEff","fCutsEff",23,xbins);
1229  for(Int_t i=1;i<=6;i++) fCutsEff->SetBinContent(i,0.05);
1230  for(Int_t i=7;i<=23;i++) fCutsEff->SetBinContent(i,0.03);
1231 
1232  // PID efficiency (from PID/noPID)
1233  fPIDEff = new TH1F("fPIDEff","fPIDEff",23,xbins);
1234  for(Int_t i=1;i<=23;i++) fPIDEff->SetBinContent(i,0.0);
1235 
1236  // MC dN/dpt
1237  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",23,xbins);
1238  fMCPtShape->SetBinContent(1,0.01);
1239  fMCPtShape->SetBinContent(2,0.01);
1240  fMCPtShape->SetBinContent(3,0.01);
1241  for(Int_t i=4;i<=23;i++) fMCPtShape->SetBinContent(i,0);
1242 
1243  // particle-antiparticle
1244  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
1245  // fPartAntipart->SetBinContent(1,1);
1246  // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
1247 
1248  return;
1249 }
1250 
1251 //--------------------------------------------------------------------------
1252 void AliHFSystErr::InitD0toKpi2017pp5TeVLowPtAn() {
1253  //
1254  // D0->Kpi syst errors. Responsible:
1255  // 2017 pp sample at 5 TeV, analysis without topological cuts
1256  //
1257 
1258  AliInfo(" Settings for D0 --> K pi, pp collisions at 5 TeV, analysis without topological cuts");
1259  SetNameTitle("AliHFSystErr","SystErrD0toKpi2017pp5TeVLowPtAn");
1260 
1261  // Normalization
1262  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1263  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.021); // 2.1% error on sigmaV0and
1264 
1265  // Branching ratio
1266  fBR = new TH1F("fBR","fBR",24,0,24);
1267  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88)
1268 
1269  // Tracking efficiency
1270  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1271  fTrackingEff->SetBinContent(1,0.03); // 0-1
1272  fTrackingEff->SetBinContent(2,0.035); // 1-2
1273  fTrackingEff->SetBinContent(3,0.04); // 2-3
1274  for(Int_t i=4;i<=6;i++) fTrackingEff->SetBinContent(i,0.045);
1275  for(Int_t i=7;i<=24;i++) fTrackingEff->SetBinContent(i,0.05);
1276 
1277  // Raw yield extraction
1278  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1279  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
1280  fRawYield->SetBinContent(1,0.05);
1281  for(Int_t i=2;i<=12;i++) fRawYield->SetBinContent(i,0.04);
1282  // fRawYield->SetBinContent(4,0.05); //5<pt<6
1283 
1284  // Cuts efficiency (from cuts variation)
1285  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1286  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.);
1287 
1288  // PID efficiency (from PID/noPID)
1289  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1290  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
1291 
1292  // MC dN/dpt
1293  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1294  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1295 
1296  // particle-antiparticle
1297  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
1298  // for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);
1299 
1300  return;
1301 }
1302 //--------------------------------------------------------------------------
1303 void AliHFSystErr::InitD0toKpi2017pp5TeVLowPtAn_finebins() {
1304  //
1305  // D0->Kpi syst errors. Responsible:
1306  // 2017 pp sample at 5 TeV, analysis without topological cuts
1307  //
1308 
1309  AliInfo(" Settings for D0 --> K pi, pp collisions at 5 TeV, analysis without topological cuts - 0.5 GeV bins");
1310  SetNameTitle("AliHFSystErr","SystErrD0toKpi2017pp5TeVLowPtAnFineBins");
1311 
1312  // Normalization
1313  fNorm = new TH1F("fNorm","fNorm",32,0,16);
1314  for(Int_t i=1;i<=32;i++) fNorm->SetBinContent(i,0.021); // 2.1% error on sigmaV0and
1315 
1316  // Branching ratio
1317  fBR = new TH1F("fBR","fBR",32,0,16);
1318  for(Int_t i=1;i<=32;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88)
1319 
1320  // Tracking efficiency
1321  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",32,0,16);
1322  fTrackingEff->SetBinContent(1,0.03); // 0-0.5
1323  fTrackingEff->SetBinContent(2,0.03); // 0.5-1.
1324  fTrackingEff->SetBinContent(3,0.035); // 1-1.5
1325  fTrackingEff->SetBinContent(4,0.035); // 1.5-2
1326  fTrackingEff->SetBinContent(5,0.04); // 2-2.5
1327  fTrackingEff->SetBinContent(6,0.04); // 2.5-3.
1328  for(Int_t i=7;i<=12;i++) fTrackingEff->SetBinContent(i,0.045);//3-6
1329  for(Int_t i=13;i<=32;i++) fTrackingEff->SetBinContent(i,0.05);// >6
1330 
1331  // Raw yield extraction
1332  fRawYield = new TH1F("fRawYield","fRawYield",32,0,16);
1333  for(Int_t i=1;i<=32;i++) fRawYield->SetBinContent(i,0.2);
1334  fRawYield->SetBinContent(1,0.09); // 0-0.5
1335  fRawYield->SetBinContent(2,0.06); // 0.5-1
1336  for(Int_t i=3;i<=24;i++) fRawYield->SetBinContent(i,0.04);
1337  fRawYield->SetBinContent(4,0.05); //5<pt<6
1338 
1339  // Cuts efficiency (from cuts variation)
1340  fCutsEff = new TH1F("fCutsEff","fCutsEff",32,0,16);
1341  for(Int_t i=1;i<=32;i++) fCutsEff->SetBinContent(i,0.);
1342 
1343  // PID efficiency (from PID/noPID)
1344  fPIDEff = new TH1F("fPIDEff","fPIDEff",32,0,16);
1345  for(Int_t i=1;i<=32;i++) fPIDEff->SetBinContent(i,0.);
1346 
1347  // MC dN/dpt
1348  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",32,0,16);
1349  for(Int_t i=1;i<=32;i++) fMCPtShape->SetBinContent(i,0);
1350 
1351  // particle-antiparticle
1352  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
1353  // for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);
1354 
1355  return;
1356 }
1357 //------------------------------------------------------------------------
1358 void AliHFSystErr::InitD0toKpi2016pp13TeV(){
1359  //
1360  // D0->Kpi syst errors
1361  // 2016 pp sample (Susanna)
1362  //
1363 
1364  AliInfo(" Settings for D0 --> K pi, pp collisions for 13 TeV");
1365  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pp13TeV");
1366 
1367  // Normalization
1368  fNorm = new TH1F("fNorm","fNorm",36,0,36);
1369  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); //
1370 
1371  // Branching ratio
1372  fBR = new TH1F("fBR","fBR",36,0,36);
1373  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88)
1374 
1375  // Tracking efficiency
1376  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1377  // for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.055); // 9% (3% per track)
1378  fTrackingEff->SetBinContent(2,0.03);
1379  fTrackingEff->SetBinContent(3,0.03);
1380  for(Int_t i=4;i<=36;i++) fTrackingEff->SetBinContent(i,0.04);
1381 
1382  // Raw yield extraction
1383  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1384  fRawYield->SetBinContent(1,1);
1385  fRawYield->SetBinContent(2,0.05); //1-2
1386  fRawYield->SetBinContent(3,0.02); //2-3
1387  fRawYield->SetBinContent(4,0.02); //3-4
1388  fRawYield->SetBinContent(5,0.02); //4-5
1389  fRawYield->SetBinContent(6,0.03); //5-6
1390  fRawYield->SetBinContent(7,0.04); //6-7
1391  fRawYield->SetBinContent(8,0.04); //7-8
1392  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04); //8-10 and 10-12
1393  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06); //12-16
1394  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.15); //16-24
1395 
1396  // Cuts efficiency (from cuts variation)
1397  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1398  fCutsEff->SetBinContent(1,1.);
1399  fCutsEff->SetBinContent(2,0.05);
1400  for(Int_t i=3;i<=36;i++) fCutsEff->SetBinContent(i,0.03);
1401 
1402  // PID efficiency (from PID/noPID)
1403  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1404  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00); // 0%
1405 
1406  // MC dN/dpt
1407  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1408  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
1409  fMCPtShape->SetBinContent(2,0.02);
1410 
1411  return;
1412 }
1413 //_________________________________________________________________________
1414 void AliHFSystErr::InitD0toKpi2013pPb0100(){
1415  //
1416  // D0->Kpi syst errors. p-Pb data sample
1417  //
1418  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
1419  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100");
1420 
1421  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1422  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
1423 
1424  // Branching ratio
1425  fBR = new TH1F("fBR","fBR",24,0,24);
1426  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.2% PDG2010
1427 
1428  // Tracking efficiency
1429  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1430  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track
1431 
1432  // Raw yield extraction
1433  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1434  fRawYield->SetBinContent(1,0.0);
1435  fRawYield->SetBinContent(2,0.08);
1436  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.04);
1437  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.07);
1438  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.11);
1439 
1440  // Cuts efficiency (from cuts variation)
1441  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1442  fCutsEff->SetBinContent(1,0.0);
1443  fCutsEff->SetBinContent(2,0.08);
1444  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
1445 
1446  // PID efficiency (from PID/noPID)
1447  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1448  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
1449 
1450  // MC dN/dpt
1451  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1452  fMCPtShape->SetBinContent(1,0.0);
1453  fMCPtShape->SetBinContent(2,0.02);
1454  for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
1455 
1456  // particle-antiparticle
1457  /*
1458  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1459  for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1460  */
1461  return;
1462 
1463 
1464 }
1465 
1466 //_________________________________________________________________________
1467 void AliHFSystErr::InitD0toKpi2013pPb0100LowPtAn(){
1468  //
1469  // D0->Kpi syst errors. p-Pb data sample
1470  // analysis without topological cuts
1471 
1472  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV, analysis without topological cuts");
1473  SetNameTitle("AliHFSystErr","SystErrD0toKpi2013pPb0100LowPtAn");
1474 
1475  // Normalization
1476  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1477  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
1478 
1479  // Branching ratio
1480  fBR = new TH1F("fBR","fBR",24,0,24);
1481  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88)
1482 
1483  // Tracking efficiency
1484  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1485  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track
1486 
1487  // Raw yield extraction
1488  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1489  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
1490  fRawYield->SetBinContent(1,0.15);
1491  fRawYield->SetBinContent(2,0.15);
1492  fRawYield->SetBinContent(3,0.10);
1493  fRawYield->SetBinContent(4,0.10);
1494  fRawYield->SetBinContent(5,0.10);
1495  fRawYield->SetBinContent(6,0.10);
1496  fRawYield->SetBinContent(7,0.10);
1497  fRawYield->SetBinContent(8,0.10);
1498  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.11);
1499 
1500  // Cuts efficiency (from cuts variation)
1501  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1502  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.);
1503 
1504  // PID efficiency (from PID/noPID)
1505  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1506  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
1507 
1508  // MC dN/dpt
1509  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1510  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1511 
1512  // particle-antiparticle
1513  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
1514  // for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);
1515 
1516  return;
1517 
1518 
1519 }
1520 //_________________________________________________________________________
1521 void AliHFSystErr::InitD0toKpi2016pPb0100(){
1522  //
1523  // D0->Kpi syst errors. p-Pb data sample
1524  //
1525  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV - 2016 run");
1526  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100");
1527 
1528  // Normalization
1529  fNorm = new TH1F("fNorm","fNorm",36,0,36);
1530  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); //
1531 
1532  // Branching ratio
1533  fBR = new TH1F("fBR","fBR",36,0,36);
1534  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // (0.05/3.88)
1535 
1536 
1537  // Tracking efficiency
1538  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1539  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); //
1540  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); //
1541  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); //
1542  for(Int_t i=24;i<=36;i++) fTrackingEff->SetBinContent(i,0.04); //
1543 
1544  // Raw yield extraction
1545  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1546  fRawYield->SetBinContent(1,0.0);
1547  fRawYield->SetBinContent(2,0.05);
1548  for(Int_t i=3;i<=17;i++) fRawYield->SetBinContent(i,0.02);
1549  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.03);
1550  for(Int_t i=24;i<=36;i++) fRawYield->SetBinContent(i,0.04);
1551 
1552  // Cuts efficiency (from cuts variation)
1553  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1554  fCutsEff->SetBinContent(1,0.0);
1555  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.03);
1556 
1557  // PID efficiency (from PID/noPID)
1558  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1559  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00);
1560 
1561  // MC dN/dpt
1562  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1563  fMCPtShape->SetBinContent(1,0.0);
1564  fMCPtShape->SetBinContent(2,0.00);
1565  for(Int_t i=3;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);//to be updated
1566 
1567  // particle-antiparticle
1568  /*
1569  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1570  for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1571  */
1572  return;
1573 
1574 
1575 }
1576 //_______________________________________________________________________
1577 void AliHFSystErr::InitD0toKpi2016pPb010ZNA(){
1578  //
1579  // D0->Kpi syst errors. p-Pb data sample
1580  //
1581  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV, 010 ZNA - 2016 run");
1582  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100");
1583 
1584  // Normalization
1585  fNorm = new TH1F("fNorm","fNorm",36,0,36);
1586  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); //
1587 
1588  // Branching ratio
1589  fBR = new TH1F("fBR","fBR",36,0,36);
1590  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // (0.05/3.88)
1591 
1592 
1593  // Tracking efficiency
1594  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1595  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); //
1596  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); //
1597  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); //
1598  for(Int_t i=24;i<=36;i++) fTrackingEff->SetBinContent(i,0.04); //
1599 
1600  // Raw yield extraction //aggiorna
1601  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1602  fRawYield->SetBinContent(1,0.0);
1603  fRawYield->SetBinContent(2,0.08);
1604  for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.04);
1605  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.03);
1606  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.035);
1607  for(Int_t i=16;i<=24;i++) fRawYield->SetBinContent(i,0.06);
1608  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10);
1609 
1610  // Cuts efficiency (from cuts variation)
1611  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1612  fCutsEff->SetBinContent(1,0.0);
1613  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.03);
1614 
1615  // PID efficiency (from PID/noPID)
1616  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1617  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00);
1618 
1619  // MC dN/dpt
1620  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1621  fMCPtShape->SetBinContent(1,0.0);
1622  fMCPtShape->SetBinContent(2,0.00);
1623  for(Int_t i=3;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);//to be updated
1624 
1625  // particle-antiparticle
1626  /*
1627  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1628  for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1629  */
1630  return;
1631 }
1632 //________________________________________________________________________
1633 void AliHFSystErr::InitD0toKpi2016pPb60100ZNA(){
1634  //
1635  // D0->Kpi syst errors. p-Pb data sample
1636  //
1637  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV 60100% ZNA - 2016 run");
1638  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100");
1639 
1640  // Normalization
1641  fNorm = new TH1F("fNorm","fNorm",36,0,36);
1642  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); //
1643 
1644  // Branching ratio
1645  fBR = new TH1F("fBR","fBR",36,0,36);
1646  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // (0.05/3.88)
1647 
1648 
1649  // Tracking efficiency
1650  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1651  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); //
1652  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); //
1653  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); //
1654  for(Int_t i=24;i<=36;i++) fTrackingEff->SetBinContent(i,0.04); //
1655 
1656  // Raw yield extraction
1657  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1658  fRawYield->SetBinContent(1,0.0);
1659  fRawYield->SetBinContent(2,0.06);
1660  for(Int_t i=3;i<=10;i++) fRawYield->SetBinContent(i,0.02);
1661  for(Int_t i=11;i<=16;i++) fRawYield->SetBinContent(i,0.03);
1662  for(Int_t i=16;i<=36;i++) fRawYield->SetBinContent(i,0.06);
1663 
1664  // Cuts efficiency (from cuts variation)
1665  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1666  fCutsEff->SetBinContent(1,0.0);
1667  for(Int_t i=2;i<=36;i++) fCutsEff->SetBinContent(i,0.03);
1668 
1669 
1670  // PID efficiency (from PID/noPID)
1671  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1672  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00);
1673 
1674  // MC dN/dpt
1675  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1676  fMCPtShape->SetBinContent(1,0.0);
1677  fMCPtShape->SetBinContent(2,0.00);
1678  for(Int_t i=3;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);//to be updated
1679 
1680  // particle-antiparticle
1681  /*
1682  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1683  for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1684  */
1685  return;
1686 }
1687 //_________________________________________________________________________
1688 void AliHFSystErr::InitD0toKpi2016pPb0100LowPtAn(){
1689  //
1690  // D0->Kpi syst errors. p-Pb 2016 data sample
1691  // analysis without topological cuts
1692 
1693  AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV 2016, analysis without topological cuts");
1694  SetNameTitle("AliHFSystErr","SystErrD0toKpi2016pPb0100LowPtAn");
1695 
1696  // Normalization
1697  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1698  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.037); // 4% error on sigmaV0and
1699 
1700  // Branching ratio
1701  fBR = new TH1F("fBR","fBR",24,0,24);
1702  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.0129); // (0.05/3.88)
1703 
1704  // Tracking efficiency
1705  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1706  for(Int_t i=0;i<=12;i++) fTrackingEff->SetBinContent(i,0.025); //
1707  for(Int_t i=12;i<=16;i++) fTrackingEff->SetBinContent(i,0.03); //
1708  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.035); //
1709 
1710  // Raw yield extraction
1711  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1712  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
1713  fRawYield->SetBinContent(1,0.06);
1714  fRawYield->SetBinContent(2,0.06);
1715  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1716 
1717  // Cuts efficiency (from cuts variation)
1718  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1719  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.);
1720 
1721  // PID efficiency (from PID/noPID)
1722  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1723  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
1724 
1725  // MC dN/dpt
1726  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1727  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1728 
1729  // particle-antiparticle
1730  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
1731  // for(Int_t i=1;i<=24;i++) fPartAntipart->SetBinContent(i,0.);
1732 
1733  return;
1734 
1735 
1736 }
1737 
1738 //--------------------------------------------------------------------------
1739 void AliHFSystErr::InitDplustoKpipi2010pp() {
1740  //
1741  // D+->Kpipi syst errors. Responsible: R. Bala
1742  // 2010 pp sample
1743  //
1744 
1745  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010pp");
1746 
1747  // Normalization
1748  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1749  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
1750 
1751  // Branching ratio
1752  fBR = new TH1F("fBR","fBR",24,0,24);
1753  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1754 
1755  // Tracking efficiency
1756  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1757  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
1758 
1759 
1760  // Raw yield extraction
1761  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1762  fRawYield->SetBinContent(1,1);
1763  fRawYield->SetBinContent(2,0.25);
1764  fRawYield->SetBinContent(3,0.25);
1765  fRawYield->SetBinContent(4,0.20);
1766  fRawYield->SetBinContent(5,0.09);
1767  fRawYield->SetBinContent(6,0.09);
1768  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1769  for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
1770 
1771  // Cuts efficiency (from cuts variation)
1772  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1773  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1774 
1775  // PID efficiency (from PID/noPID)
1776  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1777  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1778  fPIDEff->SetBinContent(1,0.15); // 15%
1779  fPIDEff->SetBinContent(2,0.15); // 15%
1780  fPIDEff->SetBinContent(3,0.15); // 15%
1781  fPIDEff->SetBinContent(4,0.15); // 15%
1782  for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%
1783 
1784  // MC dN/dpt
1785  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1786  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1787  fMCPtShape->SetBinContent(1,0.03);
1788  fMCPtShape->SetBinContent(2,0.03);
1789 
1790 
1791  // particle-antiparticle
1792  /*
1793  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1794  fPartAntipart->SetBinContent(1,1);
1795  fPartAntipart->SetBinContent(2,1);
1796  fPartAntipart->SetBinContent(3,0.12);
1797  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1798  */
1799  return;
1800 }
1801 
1802 //--------------------------------------------------------------------------
1803 void AliHFSystErr::InitDplustoKpipi2010ppPass4(){
1804  //
1805  // D+->Kpipi syst errors
1806  // 2010 pp sample - pass4
1807  //
1808 
1809 
1810  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010ppPass4");
1811  // Normalization
1812  fNorm = new TH1F("fNorm","fNorm",36,0,36);
1813  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
1814 
1815  // Branching ratio
1816  fBR = new TH1F("fBR","fBR",36,0,36);
1817  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1818 
1819  // Tracking efficiency
1820  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1821  fTrackingEff->SetBinContent(1,0.04);
1822  fTrackingEff->SetBinContent(2,0.045);
1823  fTrackingEff->SetBinContent(3,0.045);
1824  fTrackingEff->SetBinContent(4,0.05);
1825  fTrackingEff->SetBinContent(5,0.05);
1826  fTrackingEff->SetBinContent(6,0.055);
1827  for(Int_t i=7; i<=16; i++) fTrackingEff->SetBinContent(i,0.06);
1828  for(Int_t i=17; i<=36; i++) fTrackingEff->SetBinContent(i,0.055);
1829 
1830 
1831  // Raw yield extraction
1832  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1833  fRawYield->SetBinContent(1,1);
1834  fRawYield->SetBinContent(2,0.06);
1835  fRawYield->SetBinContent(3,0.06);
1836  for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.03);
1837  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
1838  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
1839  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.20);
1840 
1841  // Cuts efficiency (from cuts variation)
1842  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1843  for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i,0.1);
1844  for(Int_t i=4;i<=36;i++) fCutsEff->SetBinContent(i,0.05);
1845 
1846  // PID efficiency (from PID/noPID)
1847  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1848  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.); // 0%
1849 
1850  // MC dN/dpt
1851  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1852  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
1853  fMCPtShape->SetBinContent(2,0.03);
1854  for(Int_t i=8;i<=24;i++) fMCPtShape->SetBinContent(i,0.02);
1855 
1856 
1857  // particle-antiparticle
1858  /*
1859  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1860  fPartAntipart->SetBinContent(1,1);
1861  fPartAntipart->SetBinContent(2,1);
1862  fPartAntipart->SetBinContent(3,0.12);
1863  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1864  */
1865  return;
1866 }
1867 
1868 //--------------------------------------------------------------------------
1869 void AliHFSystErr::InitDplustoKpipi2012pp(){
1870  //
1871  // D+->Kpipi syst errors
1872  // 2012 pp sample
1873  //
1874 
1875  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2012pp");
1876  // Normalization
1877  fNorm = new TH1F("fNorm","fNorm",36,0,36);
1878  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); // 4% error on sigmaV0and
1879 
1880  // Branching ratio
1881  fBR = new TH1F("fBR","fBR",36,0,36);
1882  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); // PDG2015
1883 
1884  // Tracking efficiency
1885  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1886  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.09); // 9% (3% per track)
1887 
1888  // Raw yield extraction
1889  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1890  fRawYield->SetBinContent(1,1);
1891  fRawYield->SetBinContent(2,0.07);
1892  fRawYield->SetBinContent(3,0.05);
1893  for(Int_t i=4;i<=36;i++) fRawYield->SetBinContent(i,0.04);
1894 
1895  // Cuts efficiency (from cuts variation)
1896  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1897  fCutsEff->SetBinContent(1,1.);
1898  fCutsEff->SetBinContent(2,0.07);
1899  for(Int_t i=3;i<=36;i++) fCutsEff->SetBinContent(i,0.05); // 10%
1900 
1901  // PID efficiency (from PID/noPID)
1902  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1903  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00); // 0%
1904 
1905  // MC dN/dpt
1906  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1907  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
1908  for(Int_t i=2;i<=3;i++) fMCPtShape->SetBinContent(i,0.02);
1909  for(Int_t i=4;i<=5;i++) fMCPtShape->SetBinContent(i,0.01);
1910 
1911  return;
1912 }
1913 
1914 
1915 
1916 
1917 //--------------------------------------------------------------------------
1918 void AliHFSystErr::InitDplustoKpipi2015pp5TeV(){
1919  //
1920  // D+->Kpipi syst errors
1921  // 2015 pp sample
1922  //Responsible R. Bala
1923  AliInfo(" Settings for D+ --> K pi pi, pp collisions for 5.02 TeV");
1924  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2015pp");
1925  // Normalization
1926  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1927  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.023); // 4% error on sigmaV0and
1928 
1929  // Branching ratio
1930  fBR = new TH1F("fBR","fBR",24,0,24);
1931  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.025); // PDG2015
1932 
1933  // Tracking efficiency
1934  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1935  fTrackingEff->SetBinContent(2,0.06);
1936  fTrackingEff->SetBinContent(3,0.06);
1937  fTrackingEff->SetBinContent(4,0.06);
1938  fTrackingEff->SetBinContent(5,0.06);
1939  fTrackingEff->SetBinContent(6,0.05);
1940  fTrackingEff->SetBinContent(7,0.05);
1941  fTrackingEff->SetBinContent(8,0.05);
1942  fTrackingEff->SetBinContent(9,0.04);
1943  fTrackingEff->SetBinContent(10,0.04);
1944  for(Int_t i=11;i<=16;i++) fTrackingEff->SetBinContent(i,0.05);
1945  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.06);
1946 
1947  // Raw yield extraction
1948  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1949  fRawYield->SetBinContent(2,0.08);
1950  fRawYield->SetBinContent(3,0.04);
1951  fRawYield->SetBinContent(4,0.02);
1952  fRawYield->SetBinContent(5,0.02);
1953  fRawYield->SetBinContent(6,0.02);
1954  fRawYield->SetBinContent(7,0.04);
1955  fRawYield->SetBinContent(8,0.04);
1956  fRawYield->SetBinContent(9,0.03);
1957  fRawYield->SetBinContent(10,0.03);
1958  fRawYield->SetBinContent(11,0.06);
1959  fRawYield->SetBinContent(12,0.06);
1960 
1961 
1962  for(Int_t i=13;i<=24;i++) fRawYield->SetBinContent(i,0.04);
1963 
1964  // Cuts efficiency (from cuts variation)
1965  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1966  fCutsEff->SetBinContent(2,0.07);
1967  fCutsEff->SetBinContent(3,0.02);
1968  fCutsEff->SetBinContent(4,0.02);
1969  fCutsEff->SetBinContent(5,0.02);
1970  fCutsEff->SetBinContent(6,0.02);
1971  fCutsEff->SetBinContent(7,0.03);
1972  fCutsEff->SetBinContent(8,0.03);
1973  fCutsEff->SetBinContent(9,0.04);
1974  fCutsEff->SetBinContent(10,0.04);
1975 
1976  for(Int_t i=11;i<=16;i++) fCutsEff->SetBinContent(i,0.03);
1977  for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.04);
1978 
1979  // PID efficiency (from PID/noPID)
1980  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1981  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00); // 0%
1982 
1983  // MC dN/dpt
1984  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1985  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1986  return;
1987 }
1988 
1989 
1990 
1991 //--------------------------------------------------------------------------
1992 void AliHFSystErr::InitDplustoKpipi2017pp5TeV(){
1993  //
1994  // D+->Kpipi syst errors
1995  // 2017 pp sample
1996  //Responsible R. Bala
1997  AliInfo(" Settings for D+ --> K pi pi, pp collisions for 5.02 TeV 2017");
1998  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2017pp");
1999  // Normalization
2000  fNorm = new TH1F("fNorm","fNorm",36,0,36);
2001  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.021); // 2.1% error on sigmaV0and
2002 
2003  // Branching ratio
2004  fBR = new TH1F("fBR","fBR",36,0,36);
2005  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); // PDG2015
2006 
2007  // Tracking efficiency
2008  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
2009  fTrackingEff->SetBinContent(2,0.035);
2010  fTrackingEff->SetBinContent(3,0.045);
2011  fTrackingEff->SetBinContent(4,0.05);
2012  fTrackingEff->SetBinContent(5,0.06);
2013  fTrackingEff->SetBinContent(6,0.06);
2014  fTrackingEff->SetBinContent(7,0.065);
2015  fTrackingEff->SetBinContent(8,0.065);
2016  fTrackingEff->SetBinContent(9,0.065);
2017  fTrackingEff->SetBinContent(10,0.065);
2018  for(Int_t i=11;i<=36;i++) fTrackingEff->SetBinContent(i,0.07);
2019 
2020  // Raw yield extraction
2021  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
2022  fRawYield->SetBinContent(2,0.04);
2023  fRawYield->SetBinContent(3,0.02);
2024  fRawYield->SetBinContent(4,0.02);
2025  fRawYield->SetBinContent(5,0.02);
2026  fRawYield->SetBinContent(6,0.02);
2027  fRawYield->SetBinContent(7,0.02);
2028  fRawYield->SetBinContent(8,0.02);
2029  fRawYield->SetBinContent(9,0.02);
2030  fRawYield->SetBinContent(10,0.02);
2031  fRawYield->SetBinContent(11,0.03);
2032  fRawYield->SetBinContent(12,0.03);
2033 
2034 
2035  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.03);
2036  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.04);
2037  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.07);
2038 
2039  // Cuts efficiency (from cuts variation)
2040  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
2041  fCutsEff->SetBinContent(2,0.05);
2042  fCutsEff->SetBinContent(3,0.04);
2043  for(Int_t i=4;i<=36;i++) fCutsEff->SetBinContent(i,0.03);
2044 
2045  // PID efficiency (from PID/noPID)
2046  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
2047  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00); // 0%
2048 
2049  // MC dN/dpt
2050  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
2051  fMCPtShape->SetBinContent(2,0.04);
2052  fMCPtShape->SetBinContent(3,0.01);
2053  for(Int_t i=4;i<=36;i++) fMCPtShape->SetBinContent(i,0);
2054 
2055  return;
2056 }
2057 //--------------------------------------------------------------------------
2058 void AliHFSystErr::InitDplustoKpipi2017pp5TeV_finebins(){
2059  // DplustoKpipi syst errors. Responsible: R. Bala
2060  // 2017 pp sample - 5TeV with finer pt bins
2061  //
2062 
2063  Float_t xbins[21]={1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5,7,7.5,8,9,10,12,16,24,36};
2064  AliInfo(" Settings for D+ --> K pi pi, pp collisions for 5.02 TeV 2017 in finer bins");
2065  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2017pp5TeVFineBins");
2066 
2067  // Normalization
2068  fNorm = new TH1F("fNorm","fNorm",20,xbins);
2069  // fNorm->SetBinContent(1,0.0); // bin 0.-1.
2070  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.021);
2071 
2072  // Branching ratio
2073  fBR = new TH1F("fBR","fBR",20,xbins);
2074  //fBR->SetBinContent(1,0.0); // bin 0.-1.
2075  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.025); //
2076 
2077  // Tracking efficiency
2078  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,xbins);
2079  fTrackingEff->SetBinContent(1,0.035); // 1-1.5
2080  fTrackingEff->SetBinContent(2,0.035); // 1.5-2.
2081  fTrackingEff->SetBinContent(3,0.045); // 2-2.5
2082  fTrackingEff->SetBinContent(4,0.045); // 2.5-3
2083  fTrackingEff->SetBinContent(5,0.05); // 3-3.5
2084  fTrackingEff->SetBinContent(6,0.05); // 3.5-4.
2085  fTrackingEff->SetBinContent(7,0.06); // 4-4.5
2086  fTrackingEff->SetBinContent(8,0.06); // 4.5-5.
2087  fTrackingEff->SetBinContent(9,0.06); // 5-5.5.
2088  fTrackingEff->SetBinContent(10,0.06); // 5.5-6.
2089  fTrackingEff->SetBinContent(11,0.065); // 6-6.5.
2090  fTrackingEff->SetBinContent(12,0.065); // 6.5-7.
2091  fTrackingEff->SetBinContent(13,0.065); // 7-7.5.
2092  fTrackingEff->SetBinContent(14,0.065); // 7.5-8.
2093  fTrackingEff->SetBinContent(15,0.065); // 8-9
2094  fTrackingEff->SetBinContent(16,0.065); // 9-10
2095  for(Int_t i=17;i<=20;i++) fTrackingEff->SetBinContent(i,0.07);
2096 
2097  // Raw yield extraction
2098 
2099  fRawYield = new TH1F("fRawYield","fRawYield",20,xbins);//
2100  fRawYield->SetBinContent(1,0.05);//1-1.5
2101  fRawYield->SetBinContent(2,0.03);//1.5-2
2102  fRawYield->SetBinContent(3,0.03);// 2-2.5
2103  fRawYield->SetBinContent(4,0.03);// 2.5-3
2104  for(Int_t i=5;i<=15;i++) fRawYield->SetBinContent(i,0.02);//3-9
2105  for(Int_t i=16;i<=18;i++) fRawYield->SetBinContent(i,0.03);//9-16
2106  fRawYield->SetBinContent(19,0.04);// 16-24
2107  fRawYield->SetBinContent(20,0.07);// 24-36
2108 
2109  // Cuts efficiency (from cuts variation)
2110  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,xbins);
2111  fCutsEff->SetBinContent(1,0.05);
2112  fCutsEff->SetBinContent(2,0.05);
2113  fCutsEff->SetBinContent(3,0.04);
2114  fCutsEff->SetBinContent(4,0.04);
2115  for(Int_t i=5;i<=20;i++) fCutsEff->SetBinContent(i,0.03);
2116 
2117  // PID efficiency (from PID/noPID)
2118  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,xbins);
2119  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.0);
2120 
2121  // MC dN/dpt
2122  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,xbins);
2123  fMCPtShape->SetBinContent(1,0.04);
2124  fMCPtShape->SetBinContent(2,0.04);
2125  fMCPtShape->SetBinContent(3,0.01);
2126  fMCPtShape->SetBinContent(4,0.01);
2127  for(Int_t i=5;i<=20;i++) fMCPtShape->SetBinContent(i,0);
2128 
2129  // particle-antiparticle
2130  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
2131  // fPartAntipart->SetBinContent(1,1);
2132  // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
2133 
2134  return;
2135 }
2136 //--------------------------------------------------------------------------
2137 void AliHFSystErr::InitDplustoKpipi2016pp13TeV(){
2138  //
2139  // D+->Kpipi syst errors
2140  // 2016 pp sample
2141  //
2142 
2143  AliInfo(" Settings for D+ --> K pi pi, pp collisions for 13 TeV");
2144  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pp13TeV");
2145  // Normalization
2146  fNorm = new TH1F("fNorm","fNorm",36,0,36);
2147  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); // 4% error on sigmaV0and
2148 
2149  // Branching ratio
2150  fBR = new TH1F("fBR","fBR",36,0,36);
2151  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); // PDG2015
2152 
2153  // Tracking efficiency
2154  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
2155  // for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.09); // 9% (3% per track)
2156  fTrackingEff->SetBinContent(2,0.055);
2157  fTrackingEff->SetBinContent(3,0.055);
2158  for(Int_t i=4;i<=6;i++) fTrackingEff->SetBinContent(i,0.06); //
2159  for(Int_t i=7;i<=24;i++) fTrackingEff->SetBinContent(i,0.07); //
2160  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.075); //
2161 
2162  // Raw yield extraction
2163  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
2164  fRawYield->SetBinContent(1,1);
2165  fRawYield->SetBinContent(2,0.09);
2166  fRawYield->SetBinContent(3,0.03);
2167  for(Int_t i=4;i<=10;i++) fRawYield->SetBinContent(i,0.02);
2168  for(Int_t i=11;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2169  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.07);
2170  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.09);
2171 
2172  // Cuts efficiency (from cuts variation)
2173  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
2174  fCutsEff->SetBinContent(1,1.);
2175  fCutsEff->SetBinContent(2,0.07);
2176  fCutsEff->SetBinContent(3,0.06);
2177  for(Int_t i=4;i<=24;i++) fCutsEff->SetBinContent(i,0.04); // 10%
2178  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.04); // 10%
2179  // PID efficiency (from PID/noPID)
2180  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
2181  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00); // 0%
2182 
2183  // MC dN/dpt
2184  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
2185  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
2186  for(Int_t i=2;i<=3;i++) fMCPtShape->SetBinContent(i,0.04);
2187 
2188  return;
2189 }
2190 
2191 
2192 //--------------------------------------------------------------------------
2193 void AliHFSystErr::InitDstoKKpi2010pp() {
2194  //
2195  // D+s->KKpi syst errors. Responsible: G.M. Innocenti
2196  // 2010 pp sample
2197  //
2198 
2199 
2200  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2010pp");
2201  // Normalization
2202  fNorm = new TH1F("fNorm","fNorm",12,0,12);
2203  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
2204 
2205  // Branching ratio
2206  fBR = new TH1F("fBR","fBR",12,0,12);
2207  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
2208 
2209  // Tracking efficiency
2210  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
2211  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
2212 
2213 
2214  // Raw yield extraction
2215  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
2216  fRawYield->SetBinContent(1,1);
2217  fRawYield->SetBinContent(2,1);
2218  fRawYield->SetBinContent(3,0.20);
2219  fRawYield->SetBinContent(4,0.20);
2220  fRawYield->SetBinContent(5,0.15);
2221  fRawYield->SetBinContent(6,0.15);
2222  fRawYield->SetBinContent(7,0.15);
2223  fRawYield->SetBinContent(8,0.15);
2224  fRawYield->SetBinContent(9,0.20);
2225  fRawYield->SetBinContent(10,0.20);
2226  fRawYield->SetBinContent(11,0.20);
2227  fRawYield->SetBinContent(12,0.20);
2228 
2229  // Cuts efficiency (from cuts variation)
2230  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
2231  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%
2232 
2233  // PID efficiency (from PID/noPID)
2234  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
2235  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
2236 
2237  // MC dN/dpt
2238  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
2239  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
2240  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
2241  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
2242  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
2243  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
2244 
2245 
2246  // particle-antiparticle
2247  /*
2248  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
2249  fPartAntipart->SetBinContent(1,1);
2250  fPartAntipart->SetBinContent(2,1);
2251  fPartAntipart->SetBinContent(3,0.12);
2252  for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2253  */
2254  return;
2255 }
2256 
2257 //--------------------------------------------------------------------------
2258 void AliHFSystErr::InitDstoKKpi2010ppPass4() {
2259  //
2260  // D+s->KKpi syst errors. Responsible: A. Barbano
2261  // 2010 pp sample -- pass4
2262  //
2263 
2264 
2265  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2010ppPass4");
2266  // Normalization
2267  fNorm = new TH1F("fNorm","fNorm",12,0,12);
2268  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
2269 
2270  // Branching ratio
2271  fBR = new TH1F("fBR","fBR",12,0,12);
2272  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.035); // 0.08/2.27 PDG2015
2273 
2274  // Tracking efficiency
2275  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
2276  for(Int_t i=1; i<=2; i++) fTrackingEff->SetBinContent(i,0.03);
2277  for(Int_t i=3; i<=4; i++) fTrackingEff->SetBinContent(i,0.05);
2278  for(Int_t i=5; i<6; i++) fTrackingEff->SetBinContent(i,0.055);
2279  for(Int_t i=6; i<=12; i++) fTrackingEff->SetBinContent(i,0.06);
2280 
2281  // Raw yield extraction
2282  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
2283  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.05); //5%
2284 
2285  // Cuts efficiency (from cuts variation)
2286  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
2287  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.07); // 7%
2288 
2289  // PID efficiency (from PID/noPID)
2290  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
2291  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
2292 
2293  // MC dN/dpt
2294  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
2295  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
2296  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
2297  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
2298  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
2299  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
2300 
2301  return;
2302 }
2303 
2304 //_________________________________________________________________________
2305 void AliHFSystErr::InitDstoKKpi2017pp5TeV(){
2306  // Ds systematics in MB pp 5 TeV (2017)
2307  // Responsible: F. Grosa
2308 
2309  AliInfo(" Settings for Ds --> K K pi, MB pp collisions at 5 TeV - 2017 runs (LHC17pq)");
2310  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2017pp5TeV");
2311 
2312  // Normalization
2313  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2314  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.021); // from V0
2315 
2316  // Branching ratio
2317  fBR = new TH1F("fBR","fBR",24,0,24);
2318  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
2319 
2320  // Tracking efficiency
2321  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2322  for(Int_t i=1; i<=2; i++) fTrackingEff->SetBinContent(i, 1.000); // [0-2]
2323  for(Int_t i=3; i<=3; i++) fTrackingEff->SetBinContent(i, 0.045); // [2-3]
2324  for(Int_t i=4; i<=4; i++) fTrackingEff->SetBinContent(i, 0.055); // [3-4]
2325  for(Int_t i=5; i<=5; i++) fTrackingEff->SetBinContent(i, 0.06); // [4-5]
2326  for(Int_t i=6; i<=6; i++) fTrackingEff->SetBinContent(i, 0.065); // [5-6]
2327  for(Int_t i=7; i<=24; i++) fTrackingEff->SetBinContent(i, 0.07); // [12-24]
2328 
2329  // Raw yield extraction
2330  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2331  for(Int_t i=1; i<=2; i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
2332  for(Int_t i=3; i<=3; i++) fRawYield->SetBinContent(i, 0.07); // [2-3]
2333  for(Int_t i=4; i<=12; i++) fRawYield->SetBinContent(i, 0.03); // [3-12]
2334  for(Int_t i=13; i<=24; i++) fRawYield->SetBinContent(i, 0.05); // [12-16]
2335 
2336  // Cuts efficiency
2337  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2338  for(Int_t i=1; i<=2; i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
2339  for(Int_t i=3; i<=3; i++) fCutsEff->SetBinContent(i, 0.08); // [2-3]
2340  for(Int_t i=4; i<=24; i++) fCutsEff->SetBinContent(i, 0.05); // [3-24]
2341 
2342  // PID efficiency
2343  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2344  for(Int_t i=1; i<=2; i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
2345  for(Int_t i=3; i<=4; i++) fPIDEff->SetBinContent(i, 0.025); // [2-4]
2346  for(Int_t i=5; i<=6; i++) fPIDEff->SetBinContent(i, 0.030); // [4-6]
2347  for(Int_t i=7; i<=24; i++) fPIDEff->SetBinContent(i, 0.00); // [6-24]
2348 
2349  // MC dN/dpt
2350  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2351  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.); // [0-2]
2352  for(Int_t i=3; i<=3; i++) fMCPtShape->SetBinContent(i,0.01); // [2-3]
2353  for(Int_t i=4; i<=24; i++) fMCPtShape->SetBinContent(i,0.00); // [3-24]
2354 
2355  return;
2356 }
2357 
2358 //_________________________________________________________________________
2359 void AliHFSystErr::InitDstoKKpi2017pp5TeVBDT(){
2360  // Ds systematics in MB pp 5 TeV (2017) with BDT selections
2361  // Responsible: F. Catalano
2362 
2363  AliInfo(" Settings for Ds --> K K pi, MB pp collisions at 5 TeV - 2017 runs (LHC17pq) with BDT selections");
2364  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2017pp5TeVBDT");
2365 
2366  // Normalization
2367  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2368  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.05); // from V0
2369 
2370  // Branching ratio
2371  fBR = new TH1F("fBR","fBR",24,0,24);
2372  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
2373 
2374  // Tracking efficiency
2375  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2376  for(Int_t i=1; i<=2; i++) fTrackingEff->SetBinContent(i, 0.045); // [0-2]
2377  for(Int_t i=3; i<=3; i++) fTrackingEff->SetBinContent(i, 0.045); // [2-3]
2378  for(Int_t i=4; i<=4; i++) fTrackingEff->SetBinContent(i, 0.055); // [3-4]
2379  for(Int_t i=5; i<=5; i++) fTrackingEff->SetBinContent(i, 0.06); // [4-5]
2380  for(Int_t i=6; i<=6; i++) fTrackingEff->SetBinContent(i, 0.065); // [5-6]
2381  for(Int_t i=7; i<=16; i++) fTrackingEff->SetBinContent(i, 0.07); // [12-16]
2382  for(Int_t i=17; i<=24; i++) fTrackingEff->SetBinContent(i, 0.07); // [16-24]
2383 
2384  // Raw yield extraction
2385  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2386  for(Int_t i=1; i<=2; i++) fRawYield->SetBinContent(i, 0.1); // [0-2]
2387  for(Int_t i=3; i<=16; i++) fRawYield->SetBinContent(i, 0.02); // [2-16]
2388  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.05); // [16-24]
2389 
2390  // Cuts efficiency
2391  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2392  for(Int_t i=1; i<=2; i++) fCutsEff->SetBinContent(i, 0.15); // [0-2]
2393  for(Int_t i=3; i<=3; i++) fCutsEff->SetBinContent(i, 0.10); // [2-3]
2394  for(Int_t i=4; i<=6; i++) fCutsEff->SetBinContent(i, 0.08); // [3-6]
2395  for(Int_t i=7; i<=24; i++) fCutsEff->SetBinContent(i, 0.05); // [6-24]
2396 
2397  // PID efficiency
2398  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2399  for(Int_t i=1; i<=2; i++) fPIDEff->SetBinContent(i, 0.00); // [0-2]
2400  for(Int_t i=3; i<=6; i++) fPIDEff->SetBinContent(i, 0.00); // [2-6]
2401  for(Int_t i=7; i<=16; i++) fPIDEff->SetBinContent(i, 0.00); // [6-16]
2402  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.00); // [16-24]
2403 
2404  // MC dN/dpt
2405  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2406  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,0.01);
2407  for(Int_t i=3; i<=3; i++) fMCPtShape->SetBinContent(i,0.01);
2408  for(Int_t i=4; i<=16; i++) fMCPtShape->SetBinContent(i,0.00);
2409  for(Int_t i=17; i<=24; i++) fMCPtShape->SetBinContent(i,0.00);
2410 
2411  return;
2412 }
2413 
2414 //--------------------------------------------------------------------------
2415 void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
2416  //
2417  // D+->Kpipi syst errors. Responsible: E. Bruna
2418  // 2011 PbPb sample, 0-7.5% CC
2419  //
2420 
2421  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb07half");
2422  // Normalization
2423  fNorm = new TH1F("fNorm","fNorm",36,0,36);
2424  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
2425 
2426  // Branching ratio
2427  fBR = new TH1F("fBR","fBR",36,0,36);
2428  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
2429 
2430  // Tracking efficiency
2431  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
2432  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
2433 
2434  // Raw yield extraction
2435  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
2436  for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
2437  fRawYield->SetBinContent(4,0.30);
2438  fRawYield->SetBinContent(5,0.20);
2439  fRawYield->SetBinContent(6,0.20);
2440  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
2441  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
2442  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
2443  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
2444  for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
2445 
2446 
2447  // Cuts efficiency (from cuts variation)
2448  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
2449  for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
2450  for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);
2451 
2452 
2453  // PID efficiency (from PID/noPID)
2454  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
2455  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
2456 
2457  // MC dN/dpt (24/7/2012)
2458  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
2459  for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
2460  for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
2461  // for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
2462 
2463 
2464  // particle-antiparticle
2465  /*
2466  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2467  fPartAntipart->SetBinContent(1,1);
2468  fPartAntipart->SetBinContent(2,1);
2469  fPartAntipart->SetBinContent(3,0.12);
2470  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2471  */
2472 
2473  return;
2474 }
2475 //_________________________________________________________________
2476 void AliHFSystErr::InitDplustoKpipi2011PbPb010() {
2477  //
2478  // D+->Kpipi syst errors. Responsible: E. Bruna, C. Bedda
2479  // 2011 PbPb sample, 0-10% CC
2480  //
2481 
2482  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb010");
2483  // Normalization
2484  fNorm = new TH1F("fNorm","fNorm",36,0,36);
2485  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
2486 
2487  // Branching ratio
2488  fBR = new TH1F("fBR","fBR",36,0,36);
2489  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
2490 
2491  // Tracking efficiency
2492  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
2493  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
2494 
2495  // Raw yield extraction
2496  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
2497  for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i,0.); //5 to 10%
2498  fRawYield->SetBinContent(4,0.10);
2499  fRawYield->SetBinContent(5,0.08);
2500  fRawYield->SetBinContent(6,0.08);
2501  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.08);
2502  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
2503  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.08);
2504  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
2505  for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.08);
2506 
2507 
2508  // Cuts efficiency (from cuts variation)
2509  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
2510  for(Int_t i=1; i<=36; i++) fCutsEff->SetBinContent(i,0.10);
2511 
2512 
2513  // PID efficiency (from PID/noPID)
2514  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
2515  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
2516 
2517  // MC dN/dpt
2518  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
2519  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.);
2520  fMCPtShape->SetBinContent(4,0.06);
2521  fMCPtShape->SetBinContent(5,0.02);
2522  for (Int_t i=6; i<=36; i++) fMCPtShape->SetBinContent(i,0.01);
2523  // for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
2524 
2525 
2526  // particle-antiparticle
2527  /*
2528  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2529  fPartAntipart->SetBinContent(1,1);
2530  fPartAntipart->SetBinContent(2,1);
2531  fPartAntipart->SetBinContent(3,0.12);
2532  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2533  */
2534 
2535  return;
2536 }
2537 
2538 //___________________________________________________________________
2539 void AliHFSystErr::InitDplustoKpipi2011PbPb3050() {
2540  //
2541  // D+->Kpipi syst errors. Responsible: R. Bala
2542  // 2011 PbPb sample, 30-50% CC
2543  //
2544  AliInfo(" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
2545  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2011PbPb3050");
2546 
2547  // Normalization
2548  fNorm = new TH1F("fNorm","fNorm",16,0,16);
2549  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
2550 
2551  // Branching ratio
2552  fBR = new TH1F("fBR","fBR",16,0,16);
2553  for(Int_t i=1;i<=16;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
2554 
2555  // Tracking efficiency
2556  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
2557  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
2558 
2559  // Raw yield extraction
2560  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
2561  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.09); //5 to 10%
2562  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.12);
2563 
2564 
2565  // Cuts efficiency (from cuts variation)
2566  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
2567  for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
2568  for(Int_t i=12; i<=16; i++) fCutsEff->SetBinContent(i,0.14);
2569 
2570  // PID efficiency (from PID/noPID)
2571  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
2572  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.05); // 5%
2573 
2574  // MC dN/dpt (24/7/2012)
2575  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
2576  fMCPtShape->SetBinContent(3,0.10);
2577  fMCPtShape->SetBinContent(4,0.04);
2578  for(Int_t iBin=5; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.02);
2579 
2580 
2581  // particle-antiparticle
2582  /*
2583  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2584  fPartAntipart->SetBinContent(1,1);
2585  fPartAntipart->SetBinContent(2,1);
2586  fPartAntipart->SetBinContent(3,0.12);
2587  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2588  */
2589 
2590  return;
2591 }
2592 
2593 
2594 //--------------------------------------------------------------------------
2595 void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
2596  //
2597  // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
2598  // 2011 PbPb sample, 0-7.5% CC
2599  //
2600  AliInfo(" Settings for Ds --> K K pi, Pb-Pb collisions for 0-7.5 CC");
2601  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2011PbPb07half");
2602 
2603  // Normalization
2604  fNorm = new TH1F("fNorm","fNorm",12,0,12);
2605  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
2606 
2607  // Branching ratio
2608  fBR = new TH1F("fBR","fBR",12,0,12);
2609  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
2610 
2611  // Tracking efficiency
2612  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
2613  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
2614 
2615  // Raw yield extraction
2616  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
2617  for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30);
2618  for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
2619 
2620 
2621 
2622  // Cuts efficiency (from cuts variation)
2623  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
2624  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
2625 
2626  // PID efficiency (from PID/noPID)
2627  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
2628  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
2629 
2630  // MC dN/dpt
2631  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
2632  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
2633  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
2634  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
2635  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
2636  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
2637 
2638  // particle-antiparticle
2639  /*
2640  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
2641  fPartAntipart->SetBinContent(1,1);
2642  fPartAntipart->SetBinContent(2,1);
2643  fPartAntipart->SetBinContent(3,0.12);
2644  for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2645  */
2646 
2647  return;
2648 }
2649 
2650 //--------------------------------------------------------------------------
2651 void AliHFSystErr::InitDstoKKpi2011PbPb010() {
2652  //
2653  // D+s->Kpipi syst errors. Responsible: A.Barbano
2654  // 2011 PbPb sample, 0-10% CC
2655  //
2656 
2657  AliInfo(" Settings for Ds --> K K pi, Pb-Pb collisions for 0-10 CC");
2658  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2011PbPb010");
2659 
2660  // Normalization
2661  fNorm = new TH1F("fNorm","fNorm",12,0,12);
2662  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
2663 
2664  // Branching ratio
2665  fBR = new TH1F("fBR","fBR",12,0,12);
2666  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.045); // PDG2015
2667 
2668  // Tracking efficiency
2669  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
2670  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // 5% per track
2671 
2672  // Raw yield extraction
2673  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
2674  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.08); // 8%
2675 
2676 
2677  // Cuts efficiency (from cuts variation + phi mass selection cuts)
2678  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
2679  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
2680 
2681  // PID efficiency (from PID/noPID)
2682  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
2683  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 10%
2684 
2685  // MC dN/dpt
2686  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
2687  for(Int_t i=1; i<=6; i++) fMCPtShape->SetBinContent(i,0.02); //2%
2688  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.01); //1%
2689  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.01); //1%
2690 
2691  return;
2692 }
2693 
2694 //--------------------------------------------------------------------------
2695 void AliHFSystErr::InitDstoKKpi2011PbPb2050() {
2696  //
2697  // D+s->Kpipi syst errors. Responsible: A.Barbano
2698  // 2011 PbPb sample, 20-50% CC
2699  //
2700 
2701  AliInfo(" Settings for Ds --> K K pi, Pb-Pb collisions for 20-50 CC");
2702  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2011PbPb2050");
2703 
2704  // Normalization
2705  fNorm = new TH1F("fNorm","fNorm",12,0,12);
2706  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.051); // TAA and pp norm
2707 
2708  // Branching ratio
2709  fBR = new TH1F("fBR","fBR",12,0,12);
2710  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.045); // PDG2015
2711 
2712  // Tracking efficiency
2713  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
2714  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // 5% per track
2715 
2716  // Raw yield extraction
2717  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
2718  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.08); //8%
2719 
2720 
2721  // Cuts efficiency (from cuts variation + phi selection cuts)
2722  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
2723  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
2724 
2725  // PID efficiency (from PID/noPID)
2726  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
2727  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
2728 
2729  // MC dN/dpt
2730  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
2731  for(Int_t i=1; i<=12; i++) fMCPtShape->SetBinContent(i,0.01); //1%
2732 
2733  return;
2734 }
2735 
2736 //_________________________________________________________________________
2737 void AliHFSystErr::InitDstoKKpi2013pPb0100(){
2738  //
2739  // Dstar syst
2740  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
2741  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2013pPb0100");
2742 
2743  //Normalization
2744  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2745  for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); //
2746 
2747  // Branching ratio
2748  fBR = new TH1F("fBR","fBR",12,0,12);
2749  for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
2750 
2751  // Tracking efficiency
2752  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2753  for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
2754 
2755  // Raw yield extraction
2756  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2757  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
2758  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
2759  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
2760  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
2761 
2762  // Cuts efficiency (from cuts variation)
2763  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2764  for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //
2765 
2766  // PID efficiency (from PID/noPID)
2767  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2768  for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10);
2769  for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
2770 
2771  // MC dN/dpt
2772  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2773  for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
2774 
2775  // particle-antiparticle
2776  /*
2777  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2778  fPartAntipart->SetBinContent(1,1);
2779  fPartAntipart->SetBinContent(2,1);
2780  fPartAntipart->SetBinContent(3,0.12);
2781  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2782  */
2783  return;
2784 }
2785 
2786 //_________________________________________________________________________
2787 void AliHFSystErr::InitDstoKKpi2016pPb0100(){
2788  // Ds systematics in pPb 2016 MB
2789  // Responsible: J. Hamon
2790 
2791  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run");
2792  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb0100");
2793 
2794  // Normalization
2795  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2796  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037); // Same as D0 low pt?
2797 
2798  // Branching ratio
2799  fBR = new TH1F("fBR","fBR",24,0,24);
2800  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
2801 
2802  // Tracking efficiency
2803  // From C. Terrevoli's studies: 3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
2804  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2805  for(Int_t i=1; i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
2806  for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);
2807 
2808  // Raw yield extraction
2809  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2810  for(Int_t i=1; i<=2; i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
2811  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i, 0.03); // [2-4]
2812  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i, 0.01); // [4-6]
2813  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i, 0.01); // [6-8]
2814  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i, 0.01); // [8-12]
2815  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.03); // [12-16]
2816  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]
2817 
2818  // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
2819  // ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
2820  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2821  for(Int_t i=1; i<=2; i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
2822  for(Int_t i=3; i<=4; i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
2823  for(Int_t i=5; i<=6; i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
2824  for(Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
2825  for(Int_t i=9; i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
2826  for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
2827  for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]
2828 
2829  // PID efficiency
2830  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2831  for(Int_t i=1; i<=2; i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
2832  for(Int_t i=3; i<=4; i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
2833  for(Int_t i=5; i<=6; i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
2834  for(Int_t i=7; i<=8; i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
2835  for(Int_t i=9; i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
2836  for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.00); // [12-16] - ConservativePID
2837  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.00); // [16-24] - ConservativePID
2838 
2839  // MC dN/dpt
2840  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2841  for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.00); // No assigned systematics
2842 
2843  return;
2844 }
2845 //_________________________________________________________________________
2846 void AliHFSystErr::InitDstoKKpi2016pPb140trkl(){
2847  // Ds systematics in pPb 2016, 1-40 tracklets
2848  // Responsible: A. Barbano
2849 
2850  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run, [1,40) trkls");
2851  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb140trkls");
2852 
2853  // Normalization
2854  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2855  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037);
2856 
2857  // Branching ratio
2858  fBR = new TH1F("fBR","fBR",24,0,24);
2859  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
2860 
2861  // Tracking efficiency
2862  // From C. Terrevoli's studies: 3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
2863  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2864  for(Int_t i=1; i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
2865  for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);
2866 
2867  // Raw yield extraction
2868  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2869  for(Int_t i=1; i<=2; i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
2870  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i, 0.04); // [2-4]
2871  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i, 0.04); // [4-6]
2872  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i, 0.03); // [6-8]
2873  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i, 0.03); // [8-12]
2874  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.09); // [12-16]
2875  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]
2876 
2877  // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
2878  // ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
2879  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2880  for(Int_t i=1; i<=2; i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
2881  for(Int_t i=3; i<=4; i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
2882  for(Int_t i=5; i<=6; i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
2883  for(Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
2884  for(Int_t i=9; i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
2885  for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
2886  for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]
2887 
2888  // PID efficiency
2889  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2890  for(Int_t i=1; i<=2; i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
2891  for(Int_t i=3; i<=4; i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
2892  for(Int_t i=5; i<=6; i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
2893  for(Int_t i=7; i<=8; i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
2894  for(Int_t i=9; i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
2895  for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.02); // [12-16] - StrongPID
2896  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.02); // [16-24] - StrongPID
2897 
2898  // MC dN/dpt
2899  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2900  for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.01); // this stays for weights with Cand or with D
2901 
2902  return;
2903 }
2904 
2905 //_________________________________________________________________________
2906 void AliHFSystErr::InitDstoKKpi2016pPb4070trkl(){
2907  // Ds systematics in pPb 2016, 40-70 tracklets
2908  // Responsible: A. Barbano
2909 
2910  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run, [40,70) trkls");
2911  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb4070trkls");
2912 
2913  // Normalization
2914  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2915  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037);
2916 
2917  // Branching ratio
2918  fBR = new TH1F("fBR","fBR",24,0,24);
2919  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
2920 
2921  // Tracking efficiency
2922  // From C. Terrevoli's studies: 3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
2923  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2924  for(Int_t i=1; i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
2925  for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);
2926 
2927  // Raw yield extraction
2928  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2929  for(Int_t i=1; i<=2; i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
2930  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i, 0.03); // [2-4]
2931  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i, 0.03); // [4-6]
2932  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i, 0.04); // [6-8]
2933  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i, 0.05); // [8-12]
2934  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.12); // [12-16]
2935  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]
2936 
2937  // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
2938  // ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
2939  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2940  for(Int_t i=1; i<=2; i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
2941  for(Int_t i=3; i<=4; i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
2942  for(Int_t i=5; i<=6; i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
2943  for(Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
2944  for(Int_t i=9; i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
2945  for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
2946  for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]
2947 
2948  // PID efficiency
2949  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2950  for(Int_t i=1; i<=2; i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
2951  for(Int_t i=3; i<=4; i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
2952  for(Int_t i=5; i<=6; i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
2953  for(Int_t i=7; i<=8; i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
2954  for(Int_t i=9; i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
2955  for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.02); // [12-16] - StrongPID
2956  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.02); // [16-24] - StrongPID
2957 
2958  // MC dN/dpt
2959  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2960  for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.01); // this stays for weights with Cand or with D
2961 
2962  return;
2963 }
2964 
2965 //_________________________________________________________________________
2966 void AliHFSystErr::InitDstoKKpi2016pPb70200trkl(){
2967  // Ds systematics in pPb 2016, 70-200 tracklets
2968  // Responsible: A. Barbano
2969 
2970  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV - 2016 run, [70,200) trkls");
2971  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pPb70200trkls");
2972 
2973  // Normalization
2974  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2975  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.037);
2976 
2977  // Branching ratio
2978  fBR = new TH1F("fBR","fBR",24,0,24);
2979  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
2980 
2981  // Tracking efficiency
2982  // From C. Terrevoli's studies: 3.7 in 2-12GeV/c, 4 in 12-24 GeV/c
2983  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2984  for(Int_t i=1; i<=12; i++) fTrackingEff->SetBinContent(i, 0.037);
2985  for(Int_t i=13; i<=24; i++) fTrackingEff->SetBinContent(i, 0.040);
2986 
2987  // Raw yield extraction
2988  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2989  for(Int_t i=1; i<=2; i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
2990  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i, 0.04); // [2-4]
2991  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i, 0.04); // [4-6]
2992  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i, 0.04); // [6-8]
2993  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i, 0.04); // [8-12]
2994  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.10); // [12-16]
2995  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 0.06); // [16-24]
2996 
2997  // Cuts efficiency (envelope of cuts variation syst. and ImproverTask syst.)
2998  // ImproverTask syst. = ImproverTask_d0Mean / ImproverTask_Old
2999  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3000  for(Int_t i=1; i<=2; i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
3001  for(Int_t i=3; i<=4; i++) fCutsEff->SetBinContent(i, 0.14); // [2-4]
3002  for(Int_t i=5; i<=6; i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
3003  for(Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
3004  for(Int_t i=9; i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
3005  for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.07); // [12-16]
3006  for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 0.07); // [16-24]
3007 
3008  // PID efficiency
3009  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3010  for(Int_t i=1; i<=2; i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
3011  for(Int_t i=3; i<=4; i++) fPIDEff->SetBinContent(i, 0.02); // [2-4] - StrongPID
3012  for(Int_t i=5; i<=6; i++) fPIDEff->SetBinContent(i, 0.02); // [4-6] - StrongPID
3013  for(Int_t i=7; i<=8; i++) fPIDEff->SetBinContent(i, 0.02); // [6-8] - StrongPID
3014  for(Int_t i=9; i<=12; i++) fPIDEff->SetBinContent(i, 0.02); // [8-12] - StrongPID
3015  for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.02); // [12-16] - StrongPID
3016  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 0.02); // [16-24] - StrongPID
3017 
3018  // MC dN/dpt
3019  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3020  for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.01); // this stays for weights with Cand or with D
3021 
3022  return;
3023 }
3024 
3025 //_________________________________________________________________________
3026 void AliHFSystErr::InitDstoKKpi2016pp13TeV(){
3027  // Ds systematics in MB pp 13 TeV (2016)
3028  // Responsible: J. Hamon
3029 
3030  AliInfo(" Settings for Ds --> K K pi, MB pp collisions at 13 TeV - 2016 runs (LHC16k,l)");
3031  SetNameTitle("AliHFSystErr","SystErrDstoKKpi2016pp13TeV");
3032 
3033  // Normalization
3034  fNorm = new TH1F("fNorm","fNorm",24,0,24);
3035  for(Int_t i=1; i<=24; i++) fNorm->SetBinContent(i, 0.05); // from V0
3036 
3037  // Branching ratio
3038  fBR = new TH1F("fBR","fBR",24,0,24);
3039  for(Int_t i=1; i<=24; i++) fBR->SetBinContent(i, 0.035); // PDG 2016: BR = 2.27 +/- 0.08
3040 
3041  // Tracking efficiency
3042  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
3043  for(Int_t i=1; i<=2; i++) fTrackingEff->SetBinContent(i, 1.000); // [0-2]
3044  for(Int_t i=3; i<=4; i++) fTrackingEff->SetBinContent(i, 0.055); // [2-4]
3045  for(Int_t i=5; i<=6; i++) fTrackingEff->SetBinContent(i, 0.060); // [4-6]
3046  for(Int_t i=7; i<=12; i++) fTrackingEff->SetBinContent(i, 0.070); // [6-8] + [8-12]
3047  for(Int_t i=12; i<=16; i++) fTrackingEff->SetBinContent(i, 0.075); // [12-16]
3048  for(Int_t i=17; i<=24; i++) fTrackingEff->SetBinContent(i, 1.000); // [16-24]
3049 
3050  // Raw yield extraction
3051  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3052  for(Int_t i=1; i<=2; i++) fRawYield->SetBinContent(i, 1.00); // [0-2]
3053  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i, 0.04); // [2-4]
3054  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i, 0.04); // [4-6]
3055  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i, 0.03); // [6-8]
3056  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i, 0.04); // [8-12]
3057  for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i, 0.08); // [12-16]
3058  for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i, 1.00); // [16-24]
3059 
3060  // Cuts efficiency
3061  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3062  for(Int_t i=1; i<=2; i++) fCutsEff->SetBinContent(i, 1.00); // [0-2]
3063  for(Int_t i=3; i<=4; i++) fCutsEff->SetBinContent(i, 0.18); // [2-4]
3064  for(Int_t i=5; i<=6; i++) fCutsEff->SetBinContent(i, 0.09); // [4-6]
3065  for(Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i, 0.08); // [6-8]
3066  for(Int_t i=9; i<=12; i++) fCutsEff->SetBinContent(i, 0.08); // [8-12]
3067  for(Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i, 0.08); // [12-16]
3068  for(Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i, 1.00); // [16-24]
3069 
3070  // PID efficiency
3071  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3072  for(Int_t i=1; i<=2; i++) fPIDEff->SetBinContent(i, 1.00); // [0-2]
3073  for(Int_t i=3; i<=4; i++) fPIDEff->SetBinContent(i, 0.00); // [2-4]
3074  for(Int_t i=5; i<=6; i++) fPIDEff->SetBinContent(i, 0.00); // [4-6]
3075  for(Int_t i=7; i<=8; i++) fPIDEff->SetBinContent(i, 0.00); // [6-8]
3076  for(Int_t i=9; i<=12; i++) fPIDEff->SetBinContent(i, 0.00); // [8-12]
3077  for(Int_t i=13; i<=16; i++) fPIDEff->SetBinContent(i, 0.00); // [12-16]
3078  for(Int_t i=17; i<=24; i++) fPIDEff->SetBinContent(i, 1.00); // [16-24]
3079 
3080  // MC dN/dpt
3081  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3082  for(Int_t i=1; i<=24; i++) fMCPtShape->SetBinContent(i, 0.00); // No assigned systematics
3083  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i, 0.01); // [2-4] from toy MC pt shape
3084 
3085  return;
3086 }
3087 //--------------------------------------------------------------------------
3088 void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
3089  //
3090  // D+->Kpipi syst errors. Responsible: ??
3091  // 2010 PbPb sample, 0-20 CC
3092  //
3093 
3094  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb020");
3095  // Normalization
3096  fNorm = new TH1F("fNorm","fNorm",20,0,20);
3097  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
3098 
3099  // Branching ratio
3100  fBR = new TH1F("fBR","fBR",20,0,20);
3101  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
3102 
3103  // Tracking efficiency
3104  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
3105  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
3106 
3107  // Raw yield extraction
3108  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
3109  for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
3110  // fRawYield->SetBinContent(5,0.23);
3111  //fRawYield->SetBinContent(6,0.23);
3112  fRawYield->SetBinContent(7,0.20);
3113  fRawYield->SetBinContent(8,0.20);
3114  fRawYield->SetBinContent(9,0.15);
3115  fRawYield->SetBinContent(10,0.15);
3116  fRawYield->SetBinContent(11,0.15);
3117  fRawYield->SetBinContent(12,0.15);
3118 
3119  // Cuts efficiency (from cuts variation)
3120  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
3121  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
3122 
3123  // PID efficiency (from PID/noPID)
3124  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
3125  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
3126 
3127  // MC dN/dpt (2/2/2012)
3128  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
3129  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
3130  for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
3131  for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
3132  for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
3133 
3134 
3135  // particle-antiparticle
3136  /*
3137  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3138  fPartAntipart->SetBinContent(1,1);
3139  fPartAntipart->SetBinContent(2,1);
3140  fPartAntipart->SetBinContent(3,0.12);
3141  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3142  */
3143 
3144  return;
3145 }
3146 
3147 //--------------------------------------------------------------------------
3148 void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
3149  //
3150  // D+->Kpipi syst errors. Responsible: ??
3151  // 2010 PbPb sample, 40-80 CC
3152  //
3153 
3154 
3155  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010PbPb4080");
3156  // Normalization
3157  fNorm = new TH1F("fNorm","fNorm",20,0,20);
3158  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
3159 
3160  // Branching ratio
3161  fBR = new TH1F("fBR","fBR",20,0,20);
3162  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
3163 
3164  // Tracking efficiency
3165  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
3166  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
3167 
3168 
3169  // Raw yield extraction
3170  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
3171  fRawYield->SetBinContent(1,1);
3172  fRawYield->SetBinContent(2,1);
3173  fRawYield->SetBinContent(3,1);
3174  fRawYield->SetBinContent(4,0.15);
3175  fRawYield->SetBinContent(5,0.05);
3176  fRawYield->SetBinContent(6,0.05);
3177  fRawYield->SetBinContent(7,0.15);
3178  fRawYield->SetBinContent(8,0.15);
3179  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
3180  for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
3181 
3182  // Cuts efficiency (from cuts variation)
3183  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
3184  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
3185 
3186  // PID efficiency (from PID/noPID)
3187  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
3188  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
3189  fPIDEff->SetBinContent(3,0.13); // 13%
3190 
3191 
3192  // MC dN/dpt (2/2/2012)
3193  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
3194  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
3195  for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
3196  for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
3197  for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
3198 
3199 
3200  // particle-antiparticle
3201  /*
3202  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3203  fPartAntipart->SetBinContent(1,1);
3204  fPartAntipart->SetBinContent(2,1);
3205  fPartAntipart->SetBinContent(3,0.12);
3206  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3207  */
3208  return;
3209 }
3210 
3211 //--------------------------------------------------------------------------
3212 void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
3213 
3214  //
3215  // D+->Kpipi syst errors. Responsible: R. Bala
3216  // 2011 2.76 TeV pp sample
3217  //
3218  AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
3219  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2010ppLowEn");
3220 
3221  // Normalization
3222  fNorm = new TH1F("fNorm","fNorm",20,0,20);
3223  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); // error on sigmaV0and
3224 
3225  // Branching ratio
3226  fBR = new TH1F("fBR","fBR",20,0,20);
3227  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
3228 
3229  // Tracking efficiency
3230  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
3231  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
3232 
3233  // Raw yield extraction
3234  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
3235  fRawYield->SetBinContent(1,1);
3236  fRawYield->SetBinContent(2,1);
3237  for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10); //5 to 10%
3238  fRawYield->SetBinContent(7,0.15);
3239  fRawYield->SetBinContent(8,0.15);
3240  for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055); //5 to 10%
3241 
3242  // Cuts efficiency (from cuts variation)
3243  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
3244  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
3245 
3246  // PID efficiency (from PID/noPID)
3247  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
3248  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
3249  fPIDEff->SetBinContent(3,0.10); // 13%
3250  fPIDEff->SetBinContent(4,0.10); // 13%
3251 
3252  // MC dN/dpt (copied from D0 : will update later)
3253  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
3254  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
3255  fMCPtShape->SetBinContent(1,0.03);
3256  fMCPtShape->SetBinContent(2,0.03);
3257 
3258  return;
3259 }
3260 
3261 //_________________________________________________________________________
3262 void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
3263  //
3264  // Dplus syst
3265  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
3266  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2013pPb0100");
3267 
3268  fNorm = new TH1F("fNorm","fNorm",24,0,24);
3269  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
3270 
3271  fBR = new TH1F("fBR","fBR",20,0,24);
3272  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
3273 
3274  // Tracking efficiency
3275  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
3276  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
3277 
3278  // Raw yield extraction
3279  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3280  fRawYield->SetBinContent(1,0.);
3281  fRawYield->SetBinContent(2,0.1);
3282  fRawYield->SetBinContent(3,0.08);
3283  fRawYield->SetBinContent(4,0.05);
3284  fRawYield->SetBinContent(5,0.05);
3285  fRawYield->SetBinContent(6,0.05);
3286  fRawYield->SetBinContent(7,0.05);
3287  fRawYield->SetBinContent(8,0.05);
3288  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
3289  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3290  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
3291 
3292  // Cuts efficiency (from cuts variation)
3293  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3294  fCutsEff->SetBinContent(1,0.);
3295  fCutsEff->SetBinContent(2,0.1);
3296  fCutsEff->SetBinContent(3,0.1);
3297  fCutsEff->SetBinContent(4,0.06);
3298  fCutsEff->SetBinContent(5,0.06);
3299  fCutsEff->SetBinContent(6,0.06);
3300  fCutsEff->SetBinContent(7,0.05);
3301  fCutsEff->SetBinContent(8,0.05);
3302  for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
3303  for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
3304  for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3305 
3306  // PID efficiency (from PID/noPID)
3307  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3308  fPIDEff->SetBinContent(1,1);
3309  fPIDEff->SetBinContent(2,0.1);
3310  for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0); //
3311 
3312  // MC dN/dpt
3313  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3314  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
3315 
3316  // particle-antiparticle
3317  /*
3318  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3319  fPartAntipart->SetBinContent(1,1);
3320  fPartAntipart->SetBinContent(2,1);
3321  fPartAntipart->SetBinContent(3,0.12);
3322  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3323  */
3324  return;
3325 
3326 
3327 }
3328 
3329 //_________________________________________________________________________
3330 void AliHFSystErr::InitDplustoKpipi2016pPb0100(){
3331  //
3332  // Dplus syst
3333  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV - 2016 run");
3334  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb0100");
3335 
3336  fNorm = new TH1F("fNorm","fNorm",36,0,36);
3337  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.02); //
3338 
3339  fBR = new TH1F("fBR","fBR",36,0,36);
3340  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
3341 
3342  // Tracking efficiency
3343  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
3344  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);
3345  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);
3346  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.05);
3347 
3348 
3349  // Raw yield extraction
3350  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
3351  fRawYield->SetBinContent(2,0.03);
3352  fRawYield->SetBinContent(3,0.01);
3353  fRawYield->SetBinContent(4,0.01);
3354  fRawYield->SetBinContent(5,0.01);
3355  for(Int_t i=6;i<=12;i++) fRawYield->SetBinContent(i,0.02);
3356  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.03);
3357  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.04);
3358 
3359 
3360  // Cuts efficiency (from cuts variation)
3361  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
3362  fCutsEff->SetBinContent(2,0.12);
3363  fCutsEff->SetBinContent(3,0.07);
3364  for(Int_t i=4;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
3365  for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.05);
3366 
3367 
3368  // PID efficiency (from PID/noPID)
3369  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
3370 
3371  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); //
3372 
3373  // MC dN/dpt
3374  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
3375  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);
3376 
3377  // particle-antiparticle
3378  /*
3379  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3380  fPartAntipart->SetBinContent(1,1);
3381  fPartAntipart->SetBinContent(2,1);
3382  fPartAntipart->SetBinContent(3,0.12);
3383  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3384  */
3385  return;
3386 
3387 }
3388 
3389 //--------------------------------------------------------------------------
3390 void AliHFSystErr::InitDplustoKpipi2016pPb010ZNA() {
3391  //
3392  // D+->Kpipi syst errors. Responsible: F. Grosa
3393  // 2016 pPb sample, 0-10% ZNA
3394  //
3395  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 0-10% ZNA - 2016 run");
3396  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb010ZNA");
3397 
3398  fNorm = new TH1F("fNorm","fNorm",24,0,24);
3399  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
3400 
3401  fBR = new TH1F("fBR","fBR",24,0,24);
3402  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.025); // 0.24/9.46 PDG2016
3403 
3404  // Tracking efficiency (from MB)
3405  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
3406  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);
3407  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);
3408  //for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.05);
3409 
3410  // Raw yield extraction
3411  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3412  fRawYield->SetBinContent(2,0.05);
3413  for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.02);
3414  for(Int_t i=7;i<=16;i++) fRawYield->SetBinContent(i,0.03);
3415  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
3416 
3417  // Cuts efficiency (from cuts variation)
3418  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3419  fCutsEff->SetBinContent(2,0.07);
3420  fCutsEff->SetBinContent(3,0.04);
3421  for(Int_t i=4;i<=24;i++) fCutsEff->SetBinContent(i,0.03);
3422 
3423  // PID efficiency (from PID/noPID)
3424  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3425  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.); //
3426 
3427  // MC dN/dpt
3428  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3429  for(Int_t i=1;i<=7;i++) fMCPtShape->SetBinContent(i,0.01);
3430 
3431  return;
3432 }
3433 
3434 //--------------------------------------------------------------------------
3435 void AliHFSystErr::InitDplustoKpipi2016pPb60100ZNA() {
3436  //
3437  // D+->Kpipi syst errors. Responsible: F. Grosa
3438  // 2016 pPb sample, 60-100% ZNA
3439  //
3440  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 60-100% ZNA - 2016 run");
3441  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb60100ZNA");
3442 
3443  fNorm = new TH1F("fNorm","fNorm",24,0,24);
3444  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
3445 
3446  fBR = new TH1F("fBR","fBR",24,0,24);
3447  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.025); // 0.24/9.46 PDG2016
3448 
3449  // Tracking efficiency (from MB)
3450  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
3451  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);
3452  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);
3453  //for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.05);
3454 
3455  // Raw yield extraction
3456  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3457  fRawYield->SetBinContent(2,0.08);
3458  for(Int_t i=3;i<=7;i++) fRawYield->SetBinContent(i,0.02);
3459  for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.03);
3460  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.04);
3461  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
3462 
3463  // Cuts efficiency (from cuts variation)
3464  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3465  fCutsEff->SetBinContent(2,0.07);
3466  fCutsEff->SetBinContent(3,0.04);
3467  for(Int_t i=4;i<=24;i++) fCutsEff->SetBinContent(i,0.03);
3468 
3469  // PID efficiency (from PID/noPID)
3470  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3471  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.); //
3472 
3473  // MC dN/dpt
3474  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3475  fMCPtShape->SetBinContent(2,0.03);
3476  fMCPtShape->SetBinContent(3,0.02);
3477  for(Int_t i=4;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3478 
3479  return;
3480 }
3481 
3482 //--------------------------------------------------------------------------
3483 void AliHFSystErr::InitDplustoKpipi2016pPb140trkl() {
3484 
3485  //
3486  // D+->Kpipi syst errors. Responsible: F. Grosa
3487  // 2016 pPb sample, 1-40 tracklets
3488  //
3489  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 1-40 tracklets- 2016 run");
3490  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb140trkl");
3491 
3492  fNorm = new TH1F("fNorm","fNorm",16,0,16);
3493  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.02); //
3494 
3495  fBR = new TH1F("fBR","fBR",16,0,16);
3496  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); // 0.24/9.46 PDG2016
3497 
3498  // Tracking efficiency
3499  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
3500  for(Int_t i=1;i<=16;i++) fTrackingEff->SetBinContent(i,0.04);
3501 
3502  // Raw yield extraction
3503  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
3504  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.02);
3505  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.03);
3506 
3507  // Cuts efficiency (from cuts variation)
3508  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
3509  for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.055);
3510  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
3511 
3512  // PID efficiency (from PID/noPID)
3513  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
3514  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.); //
3515 
3516  // MC dN/dpt
3517  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
3518  for(Int_t i=1;i<=16;i++) fMCPtShape->SetBinContent(i,0.);
3519 
3520  return;
3521 }
3522 
3523 //--------------------------------------------------------------------------
3524 void AliHFSystErr::InitDplustoKpipi2016pPb4070trkl() {
3525 
3526  //
3527  // D+->Kpipi syst errors. Responsible: F. Grosa
3528  // 2016 pPb sample, 40-70 tracklets
3529  //
3530  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 1-40 tracklets- 2016 run");
3531  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb140trkl");
3532 
3533  fNorm = new TH1F("fNorm","fNorm",16,0,16);
3534  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.02); //
3535 
3536  fBR = new TH1F("fBR","fBR",16,0,16);
3537  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); // 0.24/9.46 PDG2016
3538 
3539  // Tracking efficiency
3540  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
3541  for(Int_t i=1;i<=16;i++) fTrackingEff->SetBinContent(i,0.04);
3542 
3543  // Raw yield extraction
3544  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
3545  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.02);
3546  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.03);
3547 
3548  // Cuts efficiency (from cuts variation)
3549  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
3550  for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.055);
3551  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
3552 
3553  // PID efficiency (from PID/noPID)
3554  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
3555  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.); //
3556 
3557  // MC dN/dpt
3558  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
3559  for(Int_t i=1;i<=16;i++) fMCPtShape->SetBinContent(i,0.);
3560 
3561  return;
3562 }
3563 
3564 //--------------------------------------------------------------------------
3565 void AliHFSystErr::InitDplustoKpipi2016pPb70200trkl() {
3566 
3567  //
3568  // D+->Kpipi syst errors. Responsible: F. Grosa
3569  // 2016 pPb sample, 70-200 tracklets
3570  //
3571  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV 1-40 tracklets- 2016 run");
3572  SetNameTitle("AliHFSystErr","SystErrDplustoKpipi2016pPb140trkl");
3573 
3574  fNorm = new TH1F("fNorm","fNorm",16,0,16);
3575  for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.02); //
3576 
3577  fBR = new TH1F("fBR","fBR",16,0,16);
3578  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.025); // 0.24/9.46 PDG2016
3579 
3580  // Tracking efficiency
3581  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
3582  for(Int_t i=1;i<=16;i++) fTrackingEff->SetBinContent(i,0.04);
3583 
3584  // Raw yield extraction
3585  fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
3586  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.03);
3587  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.02);
3588  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.03);
3589  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.04);
3590 
3591  // Cuts efficiency (from cuts variation)
3592  fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
3593  for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.055);
3594  for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.04);
3595 
3596  // PID efficiency (from PID/noPID)
3597  fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
3598  for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.); //
3599 
3600  // MC dN/dpt
3601  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
3602  for(Int_t i=1;i<=16;i++) fMCPtShape->SetBinContent(i,0.);
3603 
3604  return;
3605 }
3606 
3607 //--------------------------------------------------------------------------
3608 void AliHFSystErr::InitDstartoD0pi2010pp() {
3609  //
3610  // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
3611  // 2010 pp sample
3612  //
3613 
3614  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010pp");
3615  // Normalization
3616  fNorm = new TH1F("fNorm","fNorm",24,0,24);
3617  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
3618 
3619  // Branching ratio
3620  fBR = new TH1F("fBR","fBR",24,0,24);
3621  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
3622 
3623  // Tracking efficiency
3624  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
3625  fTrackingEff->SetBinContent(1,1.0);
3626  fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
3627  fTrackingEff->SetBinContent(3,0.12);
3628  fTrackingEff->SetBinContent(3,0.12);
3629  for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
3630 
3631 
3632  // Raw yield extraction
3633  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3634  fRawYield->SetBinContent(1,1.0);
3635  fRawYield->SetBinContent(2,0.10);
3636  fRawYield->SetBinContent(3,0.04);
3637  fRawYield->SetBinContent(4,0.03);
3638  fRawYield->SetBinContent(5,0.03);
3639  fRawYield->SetBinContent(6,0.05);
3640  fRawYield->SetBinContent(7,0.05);
3641  fRawYield->SetBinContent(8,0.05);
3642  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04); //4%
3643  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09); //4%
3644  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
3645 
3646  // Cuts efficiency (from cuts variation)
3647  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3648  fCutsEff->SetBinContent(2,0.22);
3649  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
3650 
3651  // PID efficiency (from PID/noPID)
3652  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3653  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
3654 
3655 
3656  // MC dN/dpt (copied from D0 : will update later)
3657  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3658  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
3659  fMCPtShape->SetBinContent(1,0.03);
3660  fMCPtShape->SetBinContent(2,0.03);
3661 
3662  return;
3663 
3664 
3665 }
3667 void AliHFSystErr::InitDstartoD0pi2010ppPass4(){
3668  //
3669  // D*+->Kpipi syst errors. Responsible: A. Grelli
3670  // 2010 pp sample - pass4
3671  //
3672  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010ppPass4");
3673 
3674  // Normalization
3675  fNorm = new TH1F("fNorm","fNorm",24,0,24);
3676  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
3677 
3678  // Branching ratio
3679  fBR = new TH1F("fBR","fBR",24,0,24);
3680  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
3681 
3682  // Tracking efficiency
3683  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
3684  fTrackingEff->SetBinContent(1,1.0);
3685  fTrackingEff->SetBinContent(2,0.06); //
3686  fTrackingEff->SetBinContent(3,0.06);
3687  fTrackingEff->SetBinContent(3,0.06);
3688  for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); //
3689 
3690 
3691  // Raw yield extraction
3692  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3693  fRawYield->SetBinContent(1,1.0);
3694  fRawYield->SetBinContent(2,0.07);
3695  fRawYield->SetBinContent(3,0.02);
3696  fRawYield->SetBinContent(4,0.02);
3697  fRawYield->SetBinContent(5,0.01);
3698  fRawYield->SetBinContent(6,0.01);
3699  fRawYield->SetBinContent(7,0.01);
3700  fRawYield->SetBinContent(8,0.01);
3701  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.02); //4%
3702  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.02); //4%
3703  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.02); //4%
3704 
3705  // Cuts efficiency (from cuts variation)
3706  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3707  fCutsEff->SetBinContent(2,0.07);
3708  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
3709  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.0); //4%
3710  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0); //4%
3711 
3712 
3713  // PID efficiency (from PID/noPID)
3714  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3715  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0); // 3%
3716 
3717 
3718  // MC dN/dpt (copied from D0 : will update later)
3719  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3720  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
3721  fMCPtShape->SetBinContent(1,0.03);
3722  fMCPtShape->SetBinContent(2,0.02);
3723 
3724 
3725  // particle-antiparticle
3726  // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
3727  // fPartAntipart->SetBinContent(1,1);
3728  // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
3729 
3730  return;
3731 }
3732 //----------------------------------------------------
3733 void AliHFSystErr::InitDstartoD0pi2017pp5TeV(){
3734  //
3735  // D*+->Kpipi syst errors. Responsible: L. Vermunt
3736  // 2017 pp sample - 5TeV
3737  //
3738  AliInfo(" Settings for D*+ --> D0 pi --> K pi pi, pp collisions for 5.02 TeV 2017");
3739  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2017pp5TeV");
3740 
3741  // Normalization
3742  fNorm = new TH1F("fNorm","fNorm",36,0,36);
3743  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.021); // 2.1% error on sigmaV0and (will be updated after VdM scan)
3744  for(Int_t i=1;i<=36;i++) fNorm->SetBinError(i,0.);
3745 
3746  // Branching ratio
3747  fBR = new TH1F("fBR","fBR",36,0,36);
3748  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // 1.3%
3749  for(Int_t i=1;i<=36;i++) fBR->SetBinError(i,0.);
3750 
3751  // Tracking efficiency (TPC contribution: 0.5% per D0 daughter + 1.0% for soft pion)
3752  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
3753  fTrackingEff->SetBinContent(2,0.04);
3754  fTrackingEff->SetBinContent(3,0.04);
3755  fTrackingEff->SetBinContent(4,0.045);
3756  for(Int_t i=5;i<=36;i++) fTrackingEff->SetBinContent(i,0.05);
3757  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinError(i,0.);
3758 
3759  // Raw yield extraction
3760  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
3761  fRawYield->SetBinContent(2,0.07);
3762  fRawYield->SetBinContent(3,0.03);
3763  fRawYield->SetBinContent(4,0.02);
3764  for(Int_t i= 5;i<=16;i++) fRawYield->SetBinContent(i,0.01);
3765  for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.02);
3766  for(Int_t i=1;i<=36;i++) fRawYield->SetBinError(i,0.);
3767 
3768  // Cuts efficiency (from cuts variation)
3769  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
3770  fCutsEff->SetBinContent(2,0.05);
3771  fCutsEff->SetBinContent(3,0.04);
3772  fCutsEff->SetBinContent(4,0.02);
3773  for(Int_t i=5;i<=36;i++) fCutsEff->SetBinContent(i,0.01);
3774  for(Int_t i=1;i<=36;i++) fCutsEff->SetBinError(i,0.);
3775 
3776  // PID efficiency (from PID/noPID)
3777  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
3778  for(Int_t i=2;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
3779  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinError(i,0.);
3780 
3781  // MC pT shape
3782  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
3783  fMCPtShape->SetBinContent(2,0.05);
3784  fMCPtShape->SetBinContent(3,0.01);
3785  for(Int_t i=4;i<=36;i++) fMCPtShape->SetBinContent(i,0);
3786  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinError(i,0.);
3787 
3788  return;
3789 }
3790 //----------------------------------------------------
3791 void AliHFSystErr::InitDstartoD0pi2017pp5TeV_finebins(){
3792  //
3793  // D*+->Kpipi syst errors. Responsible: L. Vermunt
3794  // 2017 pp sample - 5TeV
3795  //
3796  // First bin is [1-2], 0.5 GeV bins afterwards
3797  // Note: *->SetBinContent(3,* should be deleted if first bin is changed: [1,2] -> [1.5,2]
3798 
3799  AliInfo(" Settings for D*+ --> D0 pi --> K pi pi, pp collisions for 5.02 TeV 2017 (fine binning)");
3800  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2017pp5TeVFineBins");
3801 
3802  // Normalization
3803  fNorm = new TH1F("fNorm","fNorm",72,0,36); //Finer binning uses 500 MeV, so 72 bins needed
3804  for(Int_t i=1;i<=72;i++) fNorm->SetBinContent(i,0.021); // 2.1% error on sigmaV0and
3805  for(Int_t i=1;i<=72;i++) fNorm->SetBinError(i,0.);
3806 
3807  // Branching ratio
3808  fBR = new TH1F("fBR","fBR",72,0,36);
3809  for(Int_t i=1;i<=72;i++) fBR->SetBinContent(i,0.013); // 1.3%
3810  for(Int_t i=1;i<=72;i++) fBR->SetBinError(i,0.);
3811 
3812  // Tracking efficiency (TPC contribution: 0.5% per D0 daughter + 1.0% for soft pion)
3813  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",72,0,36);
3814  fTrackingEff->SetBinContent(3,0.04); //<- Should be deleted when starting with bin 1.5-2.0
3815  fTrackingEff->SetBinContent(4,0.04);
3816  fTrackingEff->SetBinContent(5,0.04);
3817  fTrackingEff->SetBinContent(6,0.04);
3818  fTrackingEff->SetBinContent(7,0.045);
3819  fTrackingEff->SetBinContent(8,0.045);
3820  for(Int_t i=9;i<=72;i++) fTrackingEff->SetBinContent(i,0.05);
3821  for(Int_t i=1;i<=72;i++) fTrackingEff->SetBinError(i,0.);
3822 
3823  // Raw yield extraction
3824  fRawYield = new TH1F("fRawYield","fRawYield",72,0,36);
3825  fRawYield->SetBinContent(3,0.07); //<- Should be deleted when starting with bin 1.5-2.0
3826  fRawYield->SetBinContent(4,0.07);
3827  fRawYield->SetBinContent(5,0.03);
3828  fRawYield->SetBinContent(6,0.03);
3829  fRawYield->SetBinContent(7,0.03);
3830  for(Int_t i= 8;i<=32;i++) fRawYield->SetBinContent(i,0.01);
3831  for(Int_t i=33;i<=72;i++) fRawYield->SetBinContent(i,0.02);
3832  for(Int_t i=1;i<=72;i++) fRawYield->SetBinError(i,0.);
3833 
3834  // Cuts efficiency (from cuts variation)
3835  fCutsEff = new TH1F("fCutsEff","fCutsEff",72,0,36);
3836  fCutsEff->SetBinContent(3,0.05); //<- Should be deleted when starting with bin 1.5-2.0
3837  fCutsEff->SetBinContent(4,0.05);
3838  fCutsEff->SetBinContent(5,0.05);
3839  fCutsEff->SetBinContent(6,0.03);
3840  fCutsEff->SetBinContent(7,0.02);
3841  fCutsEff->SetBinContent(8,0.02);
3842  for(Int_t i=9;i<=72;i++) fCutsEff->SetBinContent(i,0.01);
3843  for(Int_t i=1;i<=72;i++) fCutsEff->SetBinError(i,0.);
3844 
3845  // PID efficiency (from PID/noPID)
3846  fPIDEff = new TH1F("fPIDEff","fPIDEff",72,0,36);
3847  fPIDEff->SetBinContent(3,0.0); //<- Should be deleted when starting with bin 1.5-2.0
3848  for(Int_t i=4;i<=72;i++) fPIDEff->SetBinContent(i,0.0);
3849  for(Int_t i=1;i<=72;i++) fPIDEff->SetBinError(i,0.);
3850 
3851  // MC pT shape
3852  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",72,0,36);
3853  fMCPtShape->SetBinContent(3,0.05); //<- Should be deleted when starting with bin 1.5-2.0
3854  fMCPtShape->SetBinContent(4,0.05); //<- Should be "fMCPtShape->SetBinContent(4,0.01);" when starting with bin 1.5-2.0
3855  fMCPtShape->SetBinContent(5,0.01);
3856  for(Int_t i=6;i<=72;i++) fMCPtShape->SetBinContent(i,0);
3857  for(Int_t i=1;i<=72;i++) fMCPtShape->SetBinError(i,0.);
3858 
3859  return;
3860 
3861 }
3862 //--------------------------------------------------------------------------
3863 void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
3864 
3865  //
3866  // D+->Kpipi syst errors. Responsible: A. Grelli
3867  // 2011 2.76 TeV pp sample
3868  //
3869  AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
3870  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010ppLowEn");
3871 
3872  // Normalization
3873  fNorm = new TH1F("fNorm","fNorm",20,0,20);
3874  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); // error on sigmaV0and
3875 
3876  // Branching ratio
3877  fBR = new TH1F("fBR","fBR",20,0,20);
3878  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
3879 
3880  // Tracking efficiency
3881  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
3882  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
3883 
3884  // Raw yield extraction
3885  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
3886  fRawYield->SetBinContent(1,1);
3887  fRawYield->SetBinContent(2,1);
3888  fRawYield->SetBinContent(3,0.14);
3889  fRawYield->SetBinContent(4,0.14);
3890  fRawYield->SetBinContent(5,0.12);
3891  fRawYield->SetBinContent(6,0.12);
3892  fRawYield->SetBinContent(7,0.06);
3893  fRawYield->SetBinContent(8,0.06);
3894  fRawYield->SetBinContent(9,0.08);
3895  fRawYield->SetBinContent(10,0.08);
3896  fRawYield->SetBinContent(11,0.08);
3897  fRawYield->SetBinContent(12,0.08);
3898  for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
3899 
3900  // Cuts efficiency (from cuts variation)
3901  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
3902  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
3903  fCutsEff->SetBinContent(3,0.15);
3904  fCutsEff->SetBinContent(4,0.15);
3905  fCutsEff->SetBinContent(5,0.15);
3906  fCutsEff->SetBinContent(6,0.15);
3907  fCutsEff->SetBinContent(7,0.10);
3908  fCutsEff->SetBinContent(8,0.10);
3909  fCutsEff->SetBinContent(9,0.10);
3910  fCutsEff->SetBinContent(10,0.10);
3911  fCutsEff->SetBinContent(11,0.10);
3912  fCutsEff->SetBinContent(12,0.10);
3913 
3914  // PID efficiency (from PID/noPID)
3915  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
3916  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
3917 
3918  // MC dN/dpt
3919  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
3920  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
3921  fMCPtShape->SetBinContent(1,0.03);
3922  fMCPtShape->SetBinContent(2,0.03);
3923 
3924  return;
3925 }
3926 
3927 
3928 //--------------------------------------------------------------------------
3929 void AliHFSystErr::InitDstartoD0pi2012pp() {
3930  //
3931  // D*+->D0pi syst errors.
3932  // 2012 pp sample
3933  //
3934 
3935  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2012pp");
3936 
3937  // Normalization
3938  fNorm = new TH1F("fNorm","fNorm",24,0,24);
3939  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.05); // 5% error on sigmaV0and
3940 
3941  // Branching ratio
3942  fBR = new TH1F("fBR","fBR",24,0,24);
3943  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.3% PDG2015
3944 
3945  // Tracking efficiency
3946  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
3947  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 9% (3% per track)
3948 
3949  // Raw yield extraction
3950  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3951  fRawYield->SetBinContent(1,1.0);
3952  fRawYield->SetBinContent(2,0.08);
3953  fRawYield->SetBinContent(3,0.06);
3954  for(Int_t i=4;i<=8;i++) fRawYield->SetBinContent(i,0.03);
3955  for(Int_t i=9;i<=24;i++) fRawYield->SetBinContent(i,0.02);
3956 
3957  // Cuts efficiency (from cuts variation)
3958  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3959  fCutsEff->SetBinContent(1,1.0);
3960  for(Int_t i=2;i<=3;i++) fCutsEff->SetBinContent(i,0.10);
3961  for(Int_t i=4;i<=7;i++) fCutsEff->SetBinContent(i,0.05);
3962  for(Int_t i=8;i<=24;i++) fCutsEff->SetBinContent(i,0.01);
3963 
3964  // PID efficiency (from PID/noPID)
3965  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3966  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3967 
3968  // MC dN/dpt (copied from D0 : will update later)
3969  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3970  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
3971  for(Int_t i=2;i<=3;i++) fMCPtShape->SetBinContent(i,0.02);
3972  for(Int_t i=4;i<=5;i++) fMCPtShape->SetBinContent(i,0.01);
3973 
3974  return;
3975 }
3976 
3977 
3978 //------------------------------------------------------------------------
3979 void AliHFSystErr::InitDstartoKpipi2016pp13TeV(){
3980  //
3981  // D+->Kpipi syst errors
3982  // 2016 pp sample (Annelies Venn)
3983  //
3984 
3985  AliInfo(" Settings for D*+ --> D0 pi, pp collisions for 13 TeV");
3986  SetNameTitle("AliHFSystErr","SystErrDstartoKpipi2016pp13TeV");
3987  // Normalization
3988  fNorm = new TH1F("fNorm","fNorm",36,0,36);
3989  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.05); //
3990 
3991  // Branching ratio
3992  fBR = new TH1F("fBR","fBR",36,0,36);
3993  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // PDG2015
3994 
3995  // Tracking efficiency
3996  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
3997  // for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.055); // 9% (3% per track)
3998  fTrackingEff->SetBinContent(2,0.055);
3999  fTrackingEff->SetBinContent(3,0.055);
4000  for(Int_t i=4;i<=5;i++) fTrackingEff->SetBinContent(i,0.06);
4001  for(Int_t i=6;i<=16;i++) fTrackingEff->SetBinContent(i,0.065); //
4002  for(Int_t i=17;i<=24;i++) fTrackingEff->SetBinContent(i,0.07); //
4003  for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.075); //
4004 
4005  // Raw yield extraction
4006  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
4007  fRawYield->SetBinContent(1,1);
4008  fRawYield->SetBinContent(2,0.09);
4009  fRawYield->SetBinContent(3,0.07);
4010  for(Int_t i=4;i<=5;i++) fRawYield->SetBinContent(i,0.05);
4011  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.03);
4012  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.02);
4013  for(Int_t i=13;i<=24;i++) fRawYield->SetBinContent(i,0.03);
4014  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.21);
4015 
4016  // Cuts efficiency (from cuts variation)
4017  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
4018  fCutsEff->SetBinContent(1,1.);
4019  fCutsEff->SetBinContent(2,0.11);
4020  fCutsEff->SetBinContent(3,0.05);
4021  for(Int_t i=4;i<=24;i++) fCutsEff->SetBinContent(i,0.03); // 10%
4022  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.03); // 10%
4023  // PID efficiency (from PID/noPID)
4024  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
4025  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.00); // 0%
4026 
4027  // MC dN/dpt
4028  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
4029  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0);
4030  for(Int_t i=2;i<=3;i++) fMCPtShape->SetBinContent(i,0.04);
4031 
4032  return;
4033 }
4034 
4035 //------------------------------------------------------------------------
4036 void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
4037  //
4038  // D*+->D0pi syst errors. Responsible: A. Grelli
4039  // 2010 PbPb sample, 0-20 CC
4040  //
4041 
4042  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
4043  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb020");
4044 
4045  // Normalization
4046  fNorm = new TH1F("fNorm","fNorm",24,0,24);
4047  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
4048 
4049  // Branching ratio
4050  fBR = new TH1F("fBR","fBR",24,0,24);
4051  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
4052 
4053  // Tracking efficiency
4054  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
4055  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
4056 
4057 
4058  // Raw yield extraction
4059  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4060  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
4061  fRawYield->SetBinContent(3,0.2);
4062  fRawYield->SetBinContent(4,0.2);
4063  fRawYield->SetBinContent(5,0.2);
4064  fRawYield->SetBinContent(6,0.2);
4065 
4066  // Cuts efficiency (from cuts variation)
4067  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4068  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
4069  fCutsEff->SetBinContent(4,0.15);
4070  fCutsEff->SetBinContent(5,0.15);
4071  fCutsEff->SetBinContent(6,0.15);
4072 
4073  // PID efficiency (from PID/noPID)
4074  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4075  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
4076 
4077 
4078  // MC dN/dpt (from study on D* pt shape)
4079  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
4080  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
4081  fMCPtShape->SetBinContent(4,0.025);
4082  fMCPtShape->SetBinContent(5,0.025);
4083  fMCPtShape->SetBinContent(6,0.025);
4084  fMCPtShape->SetBinContent(7,0.04);
4085  fMCPtShape->SetBinContent(8,0.04);
4086  fMCPtShape->SetBinContent(9,0.03);
4087  fMCPtShape->SetBinContent(10,0.03);
4088  fMCPtShape->SetBinContent(11,0.03);
4089  fMCPtShape->SetBinContent(12,0.03);
4090 
4091  return;
4092 
4093 }
4094 // ----------------------------- 2011 ----------------------------------
4095 void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
4096  //
4097  // D*+->D0pi syst errors. Responsible: A. Grelli
4098  // 2011 PbPb sample, 0-7.5 CC
4099  //
4100 
4101  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
4102  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb07half");
4103 
4104 
4105  // Normalization
4106  fNorm = new TH1F("fNorm","fNorm",36,0,36);
4107  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
4108 
4109  // Branching ratio
4110  fBR = new TH1F("fBR","fBR",36,0,36);
4111  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
4112 
4113  // Tracking efficiency
4114  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
4115  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
4116 
4117 
4118  // Raw yield extraction
4119  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
4120  for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05); //4%
4121  fRawYield->SetBinContent(4,0.2);
4122  fRawYield->SetBinContent(5,0.10);
4123  fRawYield->SetBinContent(6,0.10);
4124  fRawYield->SetBinContent(7,0.08);
4125  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15); //4%
4126 
4127  // Cuts efficiency (from cuts variation)
4128  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
4129  for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
4130 
4131  // PID efficiency (from PID/noPID)
4132  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
4133  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
4134  fPIDEff->SetBinContent(4,0.09);
4135 
4136  // MC dN/dpt (from study on D* pt shape)
4137  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
4138  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
4139  fMCPtShape->SetBinContent(4,0.015);
4140  fMCPtShape->SetBinContent(5,0.015);
4141  fMCPtShape->SetBinContent(6,0.015);
4142  fMCPtShape->SetBinContent(7,0.02);
4143  fMCPtShape->SetBinContent(8,0.02);
4144  fMCPtShape->SetBinContent(9,0.03);
4145  fMCPtShape->SetBinContent(10,0.03);
4146  fMCPtShape->SetBinContent(11,0.03);
4147  fMCPtShape->SetBinContent(12,0.03);
4148 
4149 
4150 
4151  return;
4152 
4153 }
4154 //-------------------------------------------------------------------------
4155 void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
4156  //
4157  // D*+->D0pi syst errors. Responsible: A. Grelli
4158  // 2010 PbPb sample, 20-40 CC
4159  //
4160 
4161  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
4162  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb2040");
4163 
4164  // Normalization
4165  fNorm = new TH1F("fNorm","fNorm",24,0,24);
4166  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
4167 
4168  // Branching ratio
4169  fBR = new TH1F("fBR","fBR",24,0,24);
4170  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
4171 
4172  // Tracking efficiency
4173  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
4174  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
4175 
4176 
4177  // Raw yield extraction
4178  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4179  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
4180 
4181  // Cuts efficiency (from cuts variation)
4182  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4183  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
4184 
4185  // PID efficiency (from PID/noPID)
4186  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4187  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
4188 
4189 
4190  // MC dN/dpt (copied from D0 : will update later)
4191  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
4192  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
4193  fMCPtShape->SetBinContent(1,0.03);
4194  fMCPtShape->SetBinContent(2,0.03);
4195 
4196  return;
4197 
4198 }
4199 
4200 //--------------------------------------------------------------------------
4201 void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
4202  //
4203  // D*+->D0pi syst errors. Responsible: A. Grelli
4204  // 2010 PbPb sample, 40-80 CC
4205  //
4206 
4207  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
4208  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2010PbPb4080");
4209 
4210  // Normalization
4211  fNorm = new TH1F("fNorm","fNorm",24,0,24);
4212  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
4213 
4214  // Branching ratio
4215  fBR = new TH1F("fBR","fBR",24,0,24);
4216  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
4217 
4218  // Tracking efficiency
4219  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
4220  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
4221 
4222 
4223  // Raw yield extraction
4224  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4225  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
4226  fRawYield->SetBinContent(1,1);
4227  fRawYield->SetBinContent(2,0.15);
4228  fRawYield->SetBinContent(3,0.15);
4229  fRawYield->SetBinContent(4,0.15);
4230  fRawYield->SetBinContent(5,0.15);
4231  fRawYield->SetBinContent(6,0.10);
4232  fRawYield->SetBinContent(7,0.10);
4233  fRawYield->SetBinContent(8,0.10);
4234  fRawYield->SetBinContent(9,0.11);
4235  fRawYield->SetBinContent(10,0.11);
4236  fRawYield->SetBinContent(11,0.11);
4237  fRawYield->SetBinContent(12,0.11);
4238  fRawYield->SetBinContent(13,0.08);
4239  fRawYield->SetBinContent(14,0.08);
4240  fRawYield->SetBinContent(15,0.08);
4241  fRawYield->SetBinContent(16,0.08);
4242 
4243 
4244  // Cuts efficiency (from cuts variation)
4245  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4246  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
4247 
4248  // PID efficiency (from PID/noPID)
4249  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4250  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
4251 
4252 
4253  // MC dN/dpt (copied from D0 : will update later)
4254  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
4255  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
4256  fMCPtShape->SetBinContent(2,0.05);
4257  fMCPtShape->SetBinContent(3,0.05);
4258  fMCPtShape->SetBinContent(4,0.05);
4259  fMCPtShape->SetBinContent(5,0.04);
4260  fMCPtShape->SetBinContent(6,0.02);
4261  fMCPtShape->SetBinContent(7,0.04);
4262  fMCPtShape->SetBinContent(8,0.04);
4263 
4264  return;
4265 
4266 }
4267 
4268 // ----------------------------- 2011 ----------------------------------
4269 void AliHFSystErr::InitDstartoD0pi2011PbPb010() {
4270  //
4271  // D*+->D0pi syst errors. Responsible: A. Grelli
4272  // 2011 PbPb sample, 0-10% CC
4273  //
4274 
4275  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
4276  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb010");
4277 
4278 
4279  // Normalization
4280  fNorm = new TH1F("fNorm","fNorm",36,0,36);
4281  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
4282 
4283  // Branching ratio
4284  fBR = new TH1F("fBR","fBR",36,0,36);
4285  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
4286 
4287  // Tracking efficiency
4288  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
4289  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
4290 
4291 
4292  // Raw yield extraction
4293  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
4294  for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i, 0.);
4295  fRawYield->SetBinContent(4,0.10);
4296  fRawYield->SetBinContent(5,0.08);
4297  fRawYield->SetBinContent(6,0.05);
4298  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
4299  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
4300  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
4301  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
4302  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10); //4%
4303 
4304  // Cuts efficiency (from cuts variation)
4305  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
4306  for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i, 0.);
4307  fCutsEff->SetBinContent(4,0.15);
4308  for(Int_t i=5;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
4309 
4310  // PID efficiency (from PID/noPID)
4311  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
4312  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
4313 
4314  // MC dN/dpt (from study on D* pt shape)
4315  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
4316  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i, 0.);
4317  fMCPtShape->SetBinContent(4,0.06);
4318  fMCPtShape->SetBinContent(5,0.03);
4319  for(Int_t i=5;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
4320 
4321  return;
4322 
4323 }
4324 
4325 //----------------------------------------------------------------------
4326 void AliHFSystErr::InitDstartoD0pi2011PbPb3050() {
4327  //
4328  // D*+->D0pi syst errors. Responsible: A. Grelli
4329  // 2011 PbPb sample, 30-50 CC
4330  //
4331 
4332  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
4333  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2011PbPb3050");
4334 
4335  // Normalization
4336  fNorm = new TH1F("fNorm","fNorm",24,0,24);
4337  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
4338 
4339  // Branching ratio
4340  fBR = new TH1F("fBR","fBR",24,0,24);
4341  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
4342 
4343  // Tracking efficiency
4344  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
4345  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
4346 
4347 
4348  // Raw yield extraction
4349  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4350  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
4351  fRawYield->SetBinContent(1,0.12);
4352  fRawYield->SetBinContent(2,0.12);
4353  fRawYield->SetBinContent(3,0.12);
4354  fRawYield->SetBinContent(4,0.06);
4355  fRawYield->SetBinContent(5,0.06);
4356  fRawYield->SetBinContent(6,0.05);
4357  fRawYield->SetBinContent(7,0.08);
4358  fRawYield->SetBinContent(8,0.05);
4359  fRawYield->SetBinContent(9,0.05);
4360  fRawYield->SetBinContent(10,0.05);
4361  fRawYield->SetBinContent(11,0.05);
4362  fRawYield->SetBinContent(12,0.05);
4363  fRawYield->SetBinContent(13,0.05);
4364  fRawYield->SetBinContent(14,0.05);
4365  fRawYield->SetBinContent(15,0.05);
4366  fRawYield->SetBinContent(16,0.05);
4367 
4368 
4369  // Cuts efficiency (from cuts variation)
4370  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4371  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
4372  fCutsEff->SetBinContent(1,0.16);
4373  fCutsEff->SetBinContent(2,0.16);
4374  fCutsEff->SetBinContent(3,0.16);
4375  fCutsEff->SetBinContent(4,0.10);
4376  fCutsEff->SetBinContent(5,0.10);
4377  fCutsEff->SetBinContent(6,0.10);
4378  // PID efficiency (from PID/noPID)
4379  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4380  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
4381 
4382 
4383  // MC dN/dpt (copied from D0 : will update later)
4384  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
4385  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
4386  fMCPtShape->SetBinContent(2,0.11);
4387  fMCPtShape->SetBinContent(3,0.11);
4388  fMCPtShape->SetBinContent(4,0.03);
4389  fMCPtShape->SetBinContent(5,0.02);
4390  fMCPtShape->SetBinContent(6,0.01);
4391  fMCPtShape->SetBinContent(7,0.01);
4392  fMCPtShape->SetBinContent(8,0.01);
4393 
4394  return;
4395 }
4396 //_________________________________________________________________________
4397 void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
4398  // Dstar syst in pPb 2013 MB
4399  //
4400  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
4401  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2013pPb0100");
4402 
4403  //Normalization
4404  fNorm = new TH1F("fNorm","fNorm",24,0,24);
4405  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); //
4406 
4407  // Branching ratio
4408  fBR = new TH1F("fBR","fBR",24,0,24);
4409  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
4410 
4411  // Tracking efficiency
4412  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
4413  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
4414 
4415  // Raw yield extraction
4416  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4417  fRawYield->SetBinContent(1,1.0);
4418  fRawYield->SetBinContent(2,0.08);
4419  fRawYield->SetBinContent(3,0.02);
4420  fRawYield->SetBinContent(4,0.02);
4421  fRawYield->SetBinContent(5,0.05);
4422  fRawYield->SetBinContent(6,0.02);
4423  fRawYield->SetBinContent(7,0.04);
4424  fRawYield->SetBinContent(8,0.04);
4425  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
4426  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
4427  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
4428 
4429  // Cuts efficiency (from cuts variation)
4430  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4431  fCutsEff->SetBinContent(1,0.);
4432  fCutsEff->SetBinContent(2,0.10);
4433  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
4434 
4435  // PID efficiency (from PID/noPID)
4436  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4437  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
4438 
4439  // MC dN/dpt
4440  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
4441  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
4442  for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
4443  for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
4444 
4445  // particle-antiparticle
4446  /*
4447  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4448  fPartAntipart->SetBinContent(1,1);
4449  fPartAntipart->SetBinContent(2,1);
4450  fPartAntipart->SetBinContent(3,0.12);
4451  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4452  */
4453 
4454  return;
4455 
4456 }
4457 
4458 //--------------------------------------------------------------------------
4459 void AliHFSystErr::InitDstartoD0pi2016pPb0100(){
4460  // Dstar syst in pPb 2016 MB
4461  // Responsible: A.M. Veen
4462 
4463  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV 2016");
4464  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2016pPb0100");
4465 
4466  //Normalization
4467  fNorm = new TH1F("fNorm","fNorm",36,0,36);
4468  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.037); // Same as D0 low pt?
4469 
4470  // Branching ratio
4471  fBR = new TH1F("fBR","fBR",36,0,36);
4472  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // 1.3% PDG2016 (Dstar->d0+d0->kPipi uncertainty)
4473 
4474  // Tracking efficiency
4475  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36); //Dstar meson: 3.7 in 1-10, 4 in 10-24, 4.5 in 24-36
4476  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);//new version
4477  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);//
4478  for(Int_t i=25;i<=36; i++) fTrackingEff->SetBinContent(i,0.045);//
4479 
4480 
4481  // Raw yield extraction
4482  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
4483  fRawYield->SetBinContent(1,1.0);
4484  fRawYield->SetBinContent(2,0.06);
4485  fRawYield->SetBinContent(3,0.04);
4486  fRawYield->SetBinContent(4,0.03);
4487  fRawYield->SetBinContent(5,0.03);
4488  fRawYield->SetBinContent(6,0.02);
4489  fRawYield->SetBinContent(7,0.02);
4490  fRawYield->SetBinContent(8,0.02);
4491  for(Int_t i=9;i<=24;i++) fRawYield->SetBinContent(i,0.02);
4492  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.04);
4493 
4494  // Cuts efficiency (from cuts variation)
4495  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
4496  // For now uncertainty added as if we use average with full improver as uncertainty
4497  fCutsEff->SetBinContent(2,0.04);// Original values are higher then effect improver in most bins
4498  fCutsEff->SetBinContent(3,0.02);
4499  for(Int_t i=4;i<=16;i++) fCutsEff->SetBinContent(i,0.02);// Very conservative first values rounded up
4500  for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.03);// a little higher due to variations in d0 only
4501  for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.02);
4502 
4503  // PID efficiency (from PID/noPID)
4504  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
4505  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.02); // PID for 2 sigma in TPC
4506 
4507  // MC dN/dpt
4508  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
4509  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0); //No systematic assigned final
4510 
4511 
4512 
4513  return;
4514 
4515 }
4516 
4517 
4518 
4519 //--------------------------------------------------------------------------
4520 void AliHFSystErr::InitDstartoD0pi2016pPb010ZNA(){
4521  // Dstar syst in pPb 2016 010ZNA
4522  // Responsible: C. Bedda
4523 
4524  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV 2016");
4525  SetNameTitle("AliHFSystErr","SystErrDstartoD0pi2016pPb010ZNA");
4526 
4527  //Normalization
4528  fNorm = new TH1F("fNorm","fNorm",36,0,36);
4529  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.037);
4530 
4531  // Branching ratio
4532  fBR = new TH1F("fBR","fBR",36,0,36);
4533  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.013); // 1.3% PDG2016 (Dstar->d0+d0->kPipi uncertainty)
4534 
4535  // Tracking efficiency
4536  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36); //Dstar meson: 3.7 in 1-10, 4 in 10-24, 4.5 in 24-36
4537  for(Int_t i=1;i<=10;i++) fTrackingEff->SetBinContent(i,0.037);//same as MB
4538  for(Int_t i=11;i<=24;i++) fTrackingEff->SetBinContent(i,0.04);//
4539  for(Int_t i=25;i<=36; i++) fTrackingEff->SetBinContent(i,0.045);//
4540 
4541 
4542  // Raw yield extraction
4543  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
4544  fRawYield->SetBinContent(2,0.14);
4545  fRawYield->SetBinContent(3,0.07);
4546 
4547  for(Int_t i=4;i<=15;i++) fRawYield->SetBinContent(i,0.03);
4548  for(Int_t i=16;i<=36;i++) fRawYield->SetBinContent(i,0.05);
4549 
4550  // Cuts efficiency (from cuts variation)
4551  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);//same as MB
4552 
4553  fCutsEff->SetBinContent(2,0.04);
4554 
4555  for(Int_t i=3;i<=36;i++) fCutsEff->SetBinContent(i,0.02);
4556 
4557  // PID efficiency (from PID/noPID)
4558  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
4559  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.02); // PID for 2 sigma in TPC
4560 
4561  // MC dN/dpt
4562  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
4563  fMCPtShape->SetBinContent(2,0.01);
4564  for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0); //Systematics on multiplicity weights (MC pt shape negligible)
4565  for(Int_t i=24;i<=36;i++) fMCPtShape->SetBinContent(i,0.01); //Systematics on multiplicity weights (MC pt shape negligible)