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