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