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