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