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