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