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