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