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