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