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