AliPhysics  e59a9ba (e59a9ba)
 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 //--------------------------------------------------------------------------
51 AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
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 //--------------------------------------------------------------------------
4687 Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
4688  //
4689  // Get error
4690  //
4691 
4692  Int_t bin=fCutsEff->FindBin(pt);
4693 
4694  return fCutsEff->GetBinContent(bin);
4695 }
4696 //--------------------------------------------------------------------------
4697 Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
4698  //
4699  // Get error
4700  //
4701 
4702  Int_t bin=fMCPtShape->FindBin(pt);
4703 
4704  return fMCPtShape->GetBinContent(bin);
4705 }
4706 //--------------------------------------------------------------------------
4707 Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
4708  //
4709  // Get error
4710  //
4711 
4712  Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
4713 
4714  return TMath::Sqrt(err);
4715 }
4716 //--------------------------------------------------------------------------
4717 Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
4718  //
4719  // Get error
4720  //
4721 
4722  Int_t bin=fPIDEff->FindBin(pt);
4723 
4724  return fPIDEff->GetBinContent(bin);
4725 }
4726 //--------------------------------------------------------------------------
4727 Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
4728  //
4729  // Get error
4730  //
4731 
4732  Int_t bin=fTrackingEff->FindBin(pt);
4733 
4734  return fTrackingEff->GetBinContent(bin);
4735 }
4736 //--------------------------------------------------------------------------
4737 Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
4738  //
4739  // Get error
4740  //
4741 
4742  Int_t bin=fRawYield->FindBin(pt);
4743 
4744  return fRawYield->GetBinContent(bin);
4745 }
4746 //--------------------------------------------------------------------------
4747 Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
4748  //
4749  // Get error
4750  //
4751 
4752  Int_t bin=fPartAntipart->FindBin(pt);
4753 
4754  return fPartAntipart->GetBinContent(bin);
4755 }
4756 //--------------------------------------------------------------------------
4757 Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
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 {
4828  hTotErr->SetBinContent(i,GetTotalSystErr(pt));
4829  gTotErr->SetPoint(i,pt,0.);
4830  gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
4831  }
4832 
4833  }
4834  gTotErr->SetLineColor(kBlack);
4835  gTotErr->SetFillColor(kRed);
4836  gTotErr->SetFillStyle(3002);
4837  gTotErr->Draw("2");
4838  leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
4839 // hTotErr->SetLineColor(1);
4840 // hTotErr->SetLineWidth(3);
4841 // hTotErr->Draw("same");
4842 // leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
4843 
4844 
4845  fNorm->SetFillColor(1);
4846  fNorm->SetFillStyle(3002);
4847  //fNorm->Draw("same");
4848  //TH1F *hNormRefl = ReflectHisto(fNorm);
4849  //hNormRefl->Draw("same");
4850  Double_t norm = fNorm->GetBinContent(1)*100;
4851  leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
4852 
4853  if(grErrFeeddown) {
4854  grErrFeeddown->SetFillColor(kTeal-8);
4855  grErrFeeddown->SetFillStyle(3001);
4856  grErrFeeddown->Draw("2");
4857  leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
4858  }
4859  if(fTrackingEff) {
4860  fTrackingEff->SetFillColor(4);
4861  fTrackingEff->SetFillStyle(3006);
4862  fTrackingEff->Draw("same");
4863  TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
4864  hTrackingEffRefl->Draw("same");
4865  leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
4866  }
4867  if(fBR) {
4868  fBR->SetFillColor(6);
4869  fBR->SetFillStyle(3005);
4870  //fBR->SetFillStyle(3020);
4871  fBR->Draw("same");
4872  TH1F *hBRRefl = ReflectHisto(fBR);
4873  hBRRefl->Draw("same");
4874  leg->AddEntry(fBR,"Branching ratio","f");
4875  }
4876  if(fRawYield) {
4877  Int_t ci; // for color index setting
4878  ci = TColor::GetColor("#00cc00");
4879  fRawYield->SetLineColor(ci);
4880  // fRawYield->SetLineColor(3);
4881  fRawYield->SetLineWidth(3);
4882  fRawYield->Draw("same");
4883  TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
4884  hRawYieldRefl->Draw("same");
4885  leg->AddEntry(fRawYield,"Yield extraction","l");
4886  }
4887  if(fCutsEff) {
4888  fCutsEff->SetLineColor(4);
4889  fCutsEff->SetLineWidth(3);
4890  fCutsEff->Draw("same");
4891  TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
4892  hCutsEffRefl->Draw("same");
4893  leg->AddEntry(fCutsEff,"Cuts efficiency","l");
4894  }
4895  if(fPIDEff) {
4896  fPIDEff->SetLineColor(7);
4897  fPIDEff->SetLineWidth(3);
4898  fPIDEff->Draw("same");
4899  TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
4900  hPIDEffRefl->Draw("same");
4901  leg->AddEntry(fPIDEff,"PID efficiency","l");
4902  }
4903  if(fMCPtShape) {
4904  Int_t ci = TColor::GetColor("#9933ff");
4905  fMCPtShape->SetLineColor(ci);
4906  // fMCPtShape->SetLineColor(8);
4907  fMCPtShape->SetLineWidth(3);
4908  fMCPtShape->Draw("same");
4909  TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
4910  hMCPtShapeRefl->Draw("same");
4911  leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
4912  }
4913  if(fPartAntipart) {
4914  Int_t ci = TColor::GetColor("#ff6600");
4915  fPartAntipart->SetLineColor(ci);
4916  // fPartAntipart->SetLineColor(9);
4917  fPartAntipart->SetLineWidth(3);
4918  fPartAntipart->Draw("same");
4919  TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
4920  hPartAntipartRefl->Draw("same");
4921  leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
4922  }
4923 
4924 
4925  leg->Draw();
4926 
4927  cSystErr->SaveAs("RelativeSystematics.eps");
4928 
4929  return;
4930 }
4931 //-------------------------------------------------------------------------
4932 TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
4933  //
4934  // Clones and reflects histogram
4935  //
4936  TH1F *hout=(TH1F*)hin->Clone("hout");
4937  hout->Scale(-1.);
4938 
4939  return hout;
4940 }
4941 
4942 
4943 
void InitDstartoD0pi2010PbPb2040CentScan()
void InitDstartoD0pi2011PbPb2030CentScan()
void InitDplustoKpipi2011PbPb3040CentScan()
void InitDstartoD0pi2013pPb2040ZNA()
void InitD0toKpi2010PbPb6080CentScan()
void InitDplustoKpipi2010ppLowEn()
void InitDplustoKpipi2013pPb020ZNA()
void InitD0toKpi2013pPb0100()
void InitD0toKpi2013pPb60100CL1()
void InitDplustoKpipi2013pPb60100CL1()
void InitDstartoD0pi2013pPb2040V0A()
const char * title
Definition: MakeQAPdf.C:26
TString fRapidityRange
MB:0100, 0-10:010, 0-20:020 ...40-80:4080...
Definition: AliHFSystErr.h:265
void InitDstartoD0pi2010PbPb6080CentScan()
void InitDstartoD0pi2013pPb0100RapScan0104()
void InitDplustoKpipi2011PbPb3050()
void InitDplustoKpipi2013pPb0100RapScan0408()
void InitD0toKpi2013pPb2040CL1()
void InitD0toKpi2013pPb4060CL1()
TH1F * fPIDEff
cuts efficiency
Definition: AliHFSystErr.h:257
void InitDplustoKpipi2011PbPb010()
void InitDstartoD0pi2013pPb60100CL1()
void InitDplustoKpipi2011PbPb010CentScan()
void InitD0toKpi2011PbPb1020CentScan()
Bool_t fIsLowEnergy
Rapidity range fot y measurements.
Definition: AliHFSystErr.h:267
void InitDplustoKpipi2013pPb0100RapScan0401()
void InitDstartoD0pi2011PbPb3040CentScan()
void InitD0toKpi2013pPb0100RapScan0408()
TH1F * fBR
tracking efficiency
Definition: AliHFSystErr.h:255
void InitD0toKpi2010PbPb020()
void InitDstartoD0pi2013pPb0100()
void InitDstartoD0pi2012pp()
TH1F * fMCPtShape
PID efficiency.
Definition: AliHFSystErr.h:258
void InitDplustoKpipi2013pPb4060V0A()
void InitLctopK0S2010pp()
void InitDstartoD0pi2010PbPb2040()
void InitDstartoD0pi2010PbPb4060CentScan()
void InitDplustoKpipi2010PbPb1020CentScan()
void InitDstoKKpi2011PbPb2050()
Double_t GetCutsEffErr(Double_t pt) const
void InitDstartoD0pi2011PbPb010()
void InitDstartoD0pi2010PbPb020()
void InitDstartoD0pi2013pPb60100ZNA()
void InitDplustoKpipi2013pPb0100RapScan0101()
void InitDplustoKpipi2010PbPb020()
void InitDstartoD0pi2013pPb0100RapScan0408()
TH1F * ReflectHisto(TH1F *hin) const
void InitD0toKpi2010pp()
void InitD0toKpi2011PbPb010CentScan()
void InitDplustoKpipi2013pPb60100ZNA()
void InitDplustoKpipi2013pPb0100RapScan0804()
void InitD0toKpi2011PbPb3050()
void InitDstartoD0pi2011PbPb1020CentScan()
Double_t GetPartAntipartErr(Double_t pt) const
void InitDstartoD0pi2011PbPb4050CentScan()
void InitDplustoKpipi2010PbPb4060CentScan()
void InitD0toKpi2013pPb2040V0A()
void InitDplustoKpipi2010PbPb4080()
void InitD0toKpi2013pPb4060ZNA()
void InitD0toKpi2010PbPb2040CentScan()
Double_t GetSeleEffErr(Double_t pt) const
void InitD0toKpi2013pPb0100RapScan0401()
void InitD0toKpi2013pPb020V0A()
void InitD0toKpi2010PbPb010CentScan()
void InitD0toKpi2013pPb0100RapScan0804()
void InitD0toKpi2013pPb2040ZNA()
void InitDstartoD0pi2010ppLowEn()
void InitLctopKpi2010pp()
void InitDstartoD0pi2013pPb020V0A()
void InitLctopK0S2013pPb()
void InitD0toKpi2010PbPb5080CentScan()
void InitD0toKpi2013pPb0100RapScan0104()
void InitDstartoD0pi2010PbPb1020CentScan()
void InitD0toKpi2013pPb020CL1()
void InitDstartoD0pi2010PbPb010CentScan()
Double_t GetTrackingEffErr(Double_t pt) const
void InitD0toKpi2013pPb60100V0A()
void InitDstartoD0pi2011PbPb07half()
void InitD0toKpi2011PbPb07half()
void InitDstartoD0pi2013pPb4060ZNA()
void InitDstartoD0pi2013pPb0100RapScan0401()
TH1F * fRawYield
normalization
Definition: AliHFSystErr.h:253
void InitDstartoD0pi2010PbPb5080CentScan()
void InitDplustoKpipi2013pPb2040CL1()
void InitDstartoD0pi2010PbPb4080()
void InitDplustoKpipi2013pPb020V0A()
void InitDstartoD0pi2013pPb4060V0A()
AliHFSystErr(const Char_t *name="HFSystErr", const Char_t *title="")
void InitDstartoD0pi2013pPb0100RapScan0101()
void InitDplustoKpipi2010PbPb6080CentScan()
TH1F * fPartAntipart
MC dNdpt.
Definition: AliHFSystErr.h:259
void InitDstartoD0pi2013pPb020ZNA()
Bool_t fIsRapidityScan
flag fot the PbPb centrality scan
Definition: AliHFSystErr.h:271
void InitD0toKpi2013pPb4060V0A()
void InitDstoKKpi2013pPb0100()
void InitDplustoKpipi2013pPb4060ZNA()
TH1F * fTrackingEff
raw yield
Definition: AliHFSystErr.h:254
void InitDplustoKpipi2013pPb0100()
TString fCentralityClass
Collision type: pp=0, PbPb=1.
Definition: AliHFSystErr.h:263
void InitDplustoKpipi2010pp()
Double_t GetTotalSystErr(Double_t pt, Double_t feeddownErr=0) const
void InitD0toKpi2011PbPb3040CentScan()
void InitDplustoKpipi2010PbPb010CentScan()
void InitDplustoKpipi2011PbPb4050CentScan()
void InitDstartoD0pi2013pPb020CL1()
void InitD0toKpi2013pPb60100ZNA()
void InitD0toKpi2013pPb020ZNA()
void InitD0toKpi2010ppLowEn()
Bool_t fIsLowPtAnalysis
flag for the low energy (2.76TeV) run
Definition: AliHFSystErr.h:268
void InitDstoKKpi2011PbPb07half()
void InitD0toKpi2010PbPb4060CentScan()
void InitDplustoKpipi2010ppPass4()
decay
Definition: HFPtSpectrum.C:41
void Init(Int_t decay)
Function to initialize the variables/histograms.
void InitD0toKpi2010PbPb1020CentScan()
void InitDstartoD0pi2011PbPb010CentScan()
Double_t GetRawYieldErr(Double_t pt) const
ClassImp(AliAnalysisTaskCRC) AliAnalysisTaskCRC
void InitDstoKKpi2010pp()
void InitDplustoKpipi2011PbPb1020CentScan()
void InitDplustoKpipi2013pPb020CL1()
void InitD0toKpi2010ppLowPtAn()
void InitD0toKpi2010PbPb4080()
void InitDstartoD0pi2013pPb2040CL1()
void InitD0toKpi2011PbPb2030CentScan()
void InitDplustoKpipi2011PbPb07half()
void InitDstartoD0pi2010pp()
void InitDplustoKpipi2013pPb2040ZNA()
void InitD0toKpi2011PbPb010()
void InitDplustoKpipi2013pPb2040V0A()
void InitD0toKpi2011PbPb4050CentScan()
void InitDstartoD0pi2013pPb4060CL1()
void InitDplustoKpipi2013pPb4060CL1()
void InitDplustoKpipi2010PbPb2040CentScan()
void InitD0toKpi2013pPb0100RapScan0101()
const Int_t nbins
Bool_t fIsPass4Analysis
flag for the low pt analysis (no topological cuts)
Definition: AliHFSystErr.h:269
void DrawErrors(TGraphAsymmErrors *grErrFeeddown=0) const
TH1F * fCutsEff
branching ratio
Definition: AliHFSystErr.h:256
void InitDstoKKpi2011PbPb010()
Double_t GetMCPtShapeErr(Double_t pt) const
Int_t fCollisionType
Run Number (year)
Definition: AliHFSystErr.h:262
void InitDplustoKpipi2013pPb0100RapScan0104()
void InitDstartoD0pi2013pPb60100V0A()
void InitD0toKpi2011PbPb3050OutOfPlane()
void InitDstartoD0pi2013pPb0100RapScan0804()
void InitD0toKpi2013pPb0100LowPtAn()
void InitDplustoKpipi2011PbPb2030CentScan()
Int_t fRunNumber
particle=antiparticle
Definition: AliHFSystErr.h:261
void InitDstartoD0pi2011PbPb3050()
virtual ~AliHFSystErr()
Bool_t fIsCentScan
flag for the pass4 analysis
Definition: AliHFSystErr.h:270
void InitDplustoKpipi2012pp()
void InitDplustoKpipi2013pPb60100V0A()
void InitD0toKpi2011PbPb3050InPlane()
void InitDplustoKpipi2010PbPb5080CentScan()
Double_t GetPIDEffEr