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