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