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