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