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