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