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