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