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