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