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