AliPhysics  vAN-20150822 (d56cf94)
 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 
1274  // Normalization
1275  fNorm = new TH1F("fNorm","fNorm",12,0,12);
1276  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1277 
1278  // Branching ratio
1279  fBR = new TH1F("fBR","fBR",12,0,12);
1280  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
1281 
1282  // Tracking efficiency
1283  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
1284  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1285 
1286  // Raw yield extraction
1287  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
1288  for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30);
1289  for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
1290 
1291 
1292 
1293  // Cuts efficiency (from cuts variation)
1294  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
1295  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
1296 
1297  // PID efficiency (from PID/noPID)
1298  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
1299  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
1300 
1301  // MC dN/dpt
1302  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
1303  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
1304  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
1305  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
1306  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
1307  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
1308 
1309  // particle-antiparticle
1310  /*
1311  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
1312  fPartAntipart->SetBinContent(1,1);
1313  fPartAntipart->SetBinContent(2,1);
1314  fPartAntipart->SetBinContent(3,0.12);
1315  for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1316  */
1317 
1318  return;
1319 }
1320 
1321 //--------------------------------------------------------------------------
1323  //
1324  // D+s->Kpipi syst errors. Responsible: A.Barbano
1325  // 2011 PbPb sample, 0-10% CC
1326  //
1327 
1328  // Normalization
1329  fNorm = new TH1F("fNorm","fNorm",12,0,12);
1330  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.047); // TAA and pp norm
1331 
1332  // Branching ratio
1333  fBR = new TH1F("fBR","fBR",12,0,12);
1334  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.045); // PDG2015
1335 
1336  // Tracking efficiency
1337  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
1338  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // 5% per track
1339 
1340  // Raw yield extraction
1341  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
1342  for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,0.15); // 15%
1343  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10); // 10%
1344 
1345 
1346  // Cuts efficiency (from cuts variation + phi mass selection cuts)
1347  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
1348  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
1349 
1350  // PID efficiency (from PID/noPID)
1351  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
1352  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.10); // 10%
1353 
1354  // MC dN/dpt
1355  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
1356  for(Int_t i=1; i<=6; i++) fMCPtShape->SetBinContent(i,0.02); //2%
1357  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.01); //1%
1358  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.01); //1%
1359 
1360  return;
1361 }
1362 
1363 //--------------------------------------------------------------------------
1365  //
1366  // D+s->Kpipi syst errors. Responsible: A.Barbano
1367  // 2011 PbPb sample, 20-50% CC
1368  //
1369 
1370  // Normalization
1371  fNorm = new TH1F("fNorm","fNorm",12,0,12);
1372  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
1373 
1374  // Branching ratio
1375  fBR = new TH1F("fBR","fBR",12,0,12);
1376  for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.045); // PDG2015
1377 
1378  // Tracking efficiency
1379  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
1380  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // 5% per track
1381 
1382  // Raw yield extraction
1383  fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
1384  for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,.15); //15%
1385 
1386 
1387  // Cuts efficiency (from cuts variation + phi selection cuts)
1388  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
1389  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
1390 
1391  // PID efficiency (from PID/noPID)
1392  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
1393  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.10); // 10%
1394 
1395  // MC dN/dpt
1396  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
1397  for(Int_t i=1; i<=12; i++) fMCPtShape->SetBinContent(i,0.01); //1%
1398 
1399  return;
1400 }
1401 
1402 //_________________________________________________________________________
1404  //
1405  // Dstar syst
1406  AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
1407 
1408  //Normalization
1409  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1410  for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); //
1411 
1412  // Branching ratio
1413  fBR = new TH1F("fBR","fBR",12,0,12);
1414  for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
1415 
1416  // Tracking efficiency
1417  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1418  for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1419 
1420  // Raw yield extraction
1421  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1422  for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
1423  for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
1424  for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
1425  for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
1426 
1427  // Cuts efficiency (from cuts variation)
1428  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1429  for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //
1430 
1431  // PID efficiency (from PID/noPID)
1432  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1433  for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10);
1434  for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
1435 
1436  // MC dN/dpt
1437  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1438  for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
1439 
1440  // particle-antiparticle
1441  /*
1442  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1443  fPartAntipart->SetBinContent(1,1);
1444  fPartAntipart->SetBinContent(2,1);
1445  fPartAntipart->SetBinContent(3,0.12);
1446  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1447  */
1448  return;
1449 }
1450 
1451 
1452 //--------------------------------------------------------------------------
1454  //
1455  // D+->Kpipi syst errors. Responsible: ??
1456  // 2010 PbPb sample, 0-20 CC
1457  //
1458 
1459  // Normalization
1460  fNorm = new TH1F("fNorm","fNorm",20,0,20);
1461  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
1462 
1463  // Branching ratio
1464  fBR = new TH1F("fBR","fBR",20,0,20);
1465  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1466 
1467  // Tracking efficiency
1468  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1469  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1470 
1471  // Raw yield extraction
1472  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1473  for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
1474  // fRawYield->SetBinContent(5,0.23);
1475  //fRawYield->SetBinContent(6,0.23);
1476  fRawYield->SetBinContent(7,0.20);
1477  fRawYield->SetBinContent(8,0.20);
1478  fRawYield->SetBinContent(9,0.15);
1479  fRawYield->SetBinContent(10,0.15);
1480  fRawYield->SetBinContent(11,0.15);
1481  fRawYield->SetBinContent(12,0.15);
1482 
1483  // Cuts efficiency (from cuts variation)
1484  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1485  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
1486 
1487  // PID efficiency (from PID/noPID)
1488  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1489  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1490 
1491  // MC dN/dpt (2/2/2012)
1492  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1493  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
1494  for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1495  for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1496  for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1497 
1498 
1499  // particle-antiparticle
1500  /*
1501  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1502  fPartAntipart->SetBinContent(1,1);
1503  fPartAntipart->SetBinContent(2,1);
1504  fPartAntipart->SetBinContent(3,0.12);
1505  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1506  */
1507 
1508  return;
1509 }
1510 
1511 //--------------------------------------------------------------------------
1513  //
1514  // D+->Kpipi syst errors. Responsible: ??
1515  // 2010 PbPb sample, 40-80 CC
1516  //
1517 
1518 
1519  // Normalization
1520  fNorm = new TH1F("fNorm","fNorm",20,0,20);
1521  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1522 
1523  // Branching ratio
1524  fBR = new TH1F("fBR","fBR",20,0,20);
1525  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
1526 
1527  // Tracking efficiency
1528  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1529  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1530 
1531 
1532  // Raw yield extraction
1533  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1534  fRawYield->SetBinContent(1,1);
1535  fRawYield->SetBinContent(2,1);
1536  fRawYield->SetBinContent(3,1);
1537  fRawYield->SetBinContent(4,0.15);
1538  fRawYield->SetBinContent(5,0.05);
1539  fRawYield->SetBinContent(6,0.05);
1540  fRawYield->SetBinContent(7,0.15);
1541  fRawYield->SetBinContent(8,0.15);
1542  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
1543  for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
1544 
1545  // Cuts efficiency (from cuts variation)
1546  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1547  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1548 
1549  // PID efficiency (from PID/noPID)
1550  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1551  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1552  fPIDEff->SetBinContent(3,0.13); // 13%
1553 
1554 
1555  // MC dN/dpt (2/2/2012)
1556  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1557  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
1558  for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1559  for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1560  for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1561 
1562 
1563  // particle-antiparticle
1564  /*
1565  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1566  fPartAntipart->SetBinContent(1,1);
1567  fPartAntipart->SetBinContent(2,1);
1568  fPartAntipart->SetBinContent(3,0.12);
1569  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1570  */
1571  return;
1572 }
1573 
1574 //--------------------------------------------------------------------------
1576 
1577  //
1578  // D+->Kpipi syst errors. Responsible: R. Bala
1579  // 2011 2.76 TeV pp sample
1580  //
1581  AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1582 
1583  // Normalization
1584  fNorm = new TH1F("fNorm","fNorm",20,0,20);
1585  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); // error on sigmaV0and
1586 
1587  // Branching ratio
1588  fBR = new TH1F("fBR","fBR",20,0,20);
1589  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1590 
1591  // Tracking efficiency
1592  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1593  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
1594 
1595  // Raw yield extraction
1596  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1597  fRawYield->SetBinContent(1,1);
1598  fRawYield->SetBinContent(2,1);
1599  for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10); //5 to 10%
1600  fRawYield->SetBinContent(7,0.15);
1601  fRawYield->SetBinContent(8,0.15);
1602  for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055); //5 to 10%
1603 
1604  // Cuts efficiency (from cuts variation)
1605  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1606  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
1607 
1608  // PID efficiency (from PID/noPID)
1609  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1610  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1611  fPIDEff->SetBinContent(3,0.10); // 13%
1612  fPIDEff->SetBinContent(4,0.10); // 13%
1613 
1614  // MC dN/dpt (copied from D0 : will update later)
1615  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1616  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1617  fMCPtShape->SetBinContent(1,0.03);
1618  fMCPtShape->SetBinContent(2,0.03);
1619 
1620  return;
1621 }
1622 
1623 //_________________________________________________________________________
1625  //
1626  // Dplus syst
1627  AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
1628 
1629  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1630  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
1631 
1632  fBR = new TH1F("fBR","fBR",20,0,24);
1633  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1634 
1635  // Tracking efficiency
1636  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1637  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1638 
1639  // Raw yield extraction
1640  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1641  fRawYield->SetBinContent(1,0.);
1642  fRawYield->SetBinContent(2,0.1);
1643  fRawYield->SetBinContent(3,0.08);
1644  fRawYield->SetBinContent(4,0.05);
1645  fRawYield->SetBinContent(5,0.05);
1646  fRawYield->SetBinContent(6,0.05);
1647  fRawYield->SetBinContent(7,0.05);
1648  fRawYield->SetBinContent(8,0.05);
1649  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1650  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
1651  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
1652 
1653  // Cuts efficiency (from cuts variation)
1654  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1655  fCutsEff->SetBinContent(1,0.);
1656  fCutsEff->SetBinContent(2,0.1);
1657  fCutsEff->SetBinContent(3,0.1);
1658  fCutsEff->SetBinContent(4,0.06);
1659  fCutsEff->SetBinContent(5,0.06);
1660  fCutsEff->SetBinContent(6,0.06);
1661  fCutsEff->SetBinContent(7,0.05);
1662  fCutsEff->SetBinContent(8,0.05);
1663  for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
1664  for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1665  for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
1666 
1667  // PID efficiency (from PID/noPID)
1668  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1669  fPIDEff->SetBinContent(1,1);
1670  fPIDEff->SetBinContent(2,0.1);
1671  for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0); //
1672 
1673  // MC dN/dpt
1674  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1675  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
1676 
1677  // particle-antiparticle
1678  /*
1679  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1680  fPartAntipart->SetBinContent(1,1);
1681  fPartAntipart->SetBinContent(2,1);
1682  fPartAntipart->SetBinContent(3,0.12);
1683  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1684  */
1685  return;
1686 
1687 
1688 }
1689 
1690 //--------------------------------------------------------------------------
1692  //
1693  // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
1694  // 2010 pp sample
1695  //
1696 
1697  // Normalization
1698  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1699  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
1700 
1701  // Branching ratio
1702  fBR = new TH1F("fBR","fBR",24,0,24);
1703  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1704 
1705  // Tracking efficiency
1706  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1707  fTrackingEff->SetBinContent(1,1.0);
1708  fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
1709  fTrackingEff->SetBinContent(3,0.12);
1710  fTrackingEff->SetBinContent(3,0.12);
1711  for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
1712 
1713 
1714  // Raw yield extraction
1715  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1716  fRawYield->SetBinContent(1,1.0);
1717  fRawYield->SetBinContent(2,0.10);
1718  fRawYield->SetBinContent(3,0.04);
1719  fRawYield->SetBinContent(4,0.03);
1720  fRawYield->SetBinContent(5,0.03);
1721  fRawYield->SetBinContent(6,0.05);
1722  fRawYield->SetBinContent(7,0.05);
1723  fRawYield->SetBinContent(8,0.05);
1724  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04); //4%
1725  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09); //4%
1726  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
1727 
1728  // Cuts efficiency (from cuts variation)
1729  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1730  fCutsEff->SetBinContent(2,0.22);
1731  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1732 
1733  // PID efficiency (from PID/noPID)
1734  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1735  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1736 
1737 
1738  // MC dN/dpt (copied from D0 : will update later)
1739  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1740  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1741  fMCPtShape->SetBinContent(1,0.03);
1742  fMCPtShape->SetBinContent(2,0.03);
1743 
1744  return;
1745 
1746 
1747 }
1748 //--------------------------------------------------------------------------
1750 
1751  //
1752  // D+->Kpipi syst errors. Responsible: A. Grelli
1753  // 2011 2.76 TeV pp sample
1754  //
1755  AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1756 
1757 // Normalization
1758  fNorm = new TH1F("fNorm","fNorm",20,0,20);
1759  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.019); // error on sigmaV0and
1760 
1761  // Branching ratio
1762  fBR = new TH1F("fBR","fBR",20,0,20);
1763  for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1764 
1765  // Tracking efficiency
1766  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1767  for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
1768 
1769  // Raw yield extraction
1770  fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1771  fRawYield->SetBinContent(1,1);
1772  fRawYield->SetBinContent(2,1);
1773  fRawYield->SetBinContent(3,0.14);
1774  fRawYield->SetBinContent(4,0.14);
1775  fRawYield->SetBinContent(5,0.12);
1776  fRawYield->SetBinContent(6,0.12);
1777  fRawYield->SetBinContent(7,0.06);
1778  fRawYield->SetBinContent(8,0.06);
1779  fRawYield->SetBinContent(9,0.08);
1780  fRawYield->SetBinContent(10,0.08);
1781  fRawYield->SetBinContent(11,0.08);
1782  fRawYield->SetBinContent(12,0.08);
1783  for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
1784 
1785  // Cuts efficiency (from cuts variation)
1786  fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1787  for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
1788  fCutsEff->SetBinContent(3,0.15);
1789  fCutsEff->SetBinContent(4,0.15);
1790  fCutsEff->SetBinContent(5,0.15);
1791  fCutsEff->SetBinContent(6,0.15);
1792  fCutsEff->SetBinContent(7,0.10);
1793  fCutsEff->SetBinContent(8,0.10);
1794  fCutsEff->SetBinContent(9,0.10);
1795  fCutsEff->SetBinContent(10,0.10);
1796  fCutsEff->SetBinContent(11,0.10);
1797  fCutsEff->SetBinContent(12,0.10);
1798 
1799  // PID efficiency (from PID/noPID)
1800  fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1801  for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
1802 
1803  // MC dN/dpt
1804  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1805  for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1806  fMCPtShape->SetBinContent(1,0.03);
1807  fMCPtShape->SetBinContent(2,0.03);
1808 
1809  return;
1810 }
1811 
1812 //------------------------------------------------------------------------
1814  //
1815  // D*+->D0pi syst errors. Responsible: A. Grelli
1816  // 2010 PbPb sample, 0-20 CC
1817  //
1818 
1819  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1820 
1821  // Normalization
1822  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1823  for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
1824 
1825  // Branching ratio
1826  fBR = new TH1F("fBR","fBR",24,0,24);
1827  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1828 
1829  // Tracking efficiency
1830  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1831  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1832 
1833 
1834  // Raw yield extraction
1835  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1836  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
1837  fRawYield->SetBinContent(3,0.2);
1838  fRawYield->SetBinContent(4,0.2);
1839  fRawYield->SetBinContent(5,0.2);
1840  fRawYield->SetBinContent(6,0.2);
1841 
1842  // Cuts efficiency (from cuts variation)
1843  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1844  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1845  fCutsEff->SetBinContent(4,0.15);
1846  fCutsEff->SetBinContent(5,0.15);
1847  fCutsEff->SetBinContent(6,0.15);
1848 
1849  // PID efficiency (from PID/noPID)
1850  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1851  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1852 
1853 
1854  // MC dN/dpt (from study on D* pt shape)
1855  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1856  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
1857  fMCPtShape->SetBinContent(4,0.025);
1858  fMCPtShape->SetBinContent(5,0.025);
1859  fMCPtShape->SetBinContent(6,0.025);
1860  fMCPtShape->SetBinContent(7,0.04);
1861  fMCPtShape->SetBinContent(8,0.04);
1862  fMCPtShape->SetBinContent(9,0.03);
1863  fMCPtShape->SetBinContent(10,0.03);
1864  fMCPtShape->SetBinContent(11,0.03);
1865  fMCPtShape->SetBinContent(12,0.03);
1866 
1867  return;
1868 
1869 }
1870 // ----------------------------- 2011 ----------------------------------
1872  //
1873  // D*+->D0pi syst errors. Responsible: A. Grelli
1874  // 2011 PbPb sample, 0-7.5 CC
1875  //
1876 
1877  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1878 
1879 
1880  // Normalization
1881  fNorm = new TH1F("fNorm","fNorm",36,0,36);
1882  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1883 
1884  // Branching ratio
1885  fBR = new TH1F("fBR","fBR",36,0,36);
1886  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1887 
1888  // Tracking efficiency
1889  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1890  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1891 
1892 
1893  // Raw yield extraction
1894  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1895  for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05); //4%
1896  fRawYield->SetBinContent(4,0.2);
1897  fRawYield->SetBinContent(5,0.10);
1898  fRawYield->SetBinContent(6,0.10);
1899  fRawYield->SetBinContent(7,0.08);
1900  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15); //4%
1901 
1902  // Cuts efficiency (from cuts variation)
1903  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1904  for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1905 
1906  // PID efficiency (from PID/noPID)
1907  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1908  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1909  fPIDEff->SetBinContent(4,0.09);
1910 
1911  // MC dN/dpt (from study on D* pt shape)
1912  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1913  for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
1914  fMCPtShape->SetBinContent(4,0.015);
1915  fMCPtShape->SetBinContent(5,0.015);
1916  fMCPtShape->SetBinContent(6,0.015);
1917  fMCPtShape->SetBinContent(7,0.02);
1918  fMCPtShape->SetBinContent(8,0.02);
1919  fMCPtShape->SetBinContent(9,0.03);
1920  fMCPtShape->SetBinContent(10,0.03);
1921  fMCPtShape->SetBinContent(11,0.03);
1922  fMCPtShape->SetBinContent(12,0.03);
1923 
1924 
1925 
1926  return;
1927 
1928 }
1929 //-------------------------------------------------------------------------
1931  //
1932  // D*+->D0pi syst errors. Responsible: A. Grelli
1933  // 2010 PbPb sample, 20-40 CC
1934  //
1935 
1936  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1937 
1938  // Normalization
1939  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1940  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
1941 
1942  // Branching ratio
1943  fBR = new TH1F("fBR","fBR",24,0,24);
1944  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1945 
1946  // Tracking efficiency
1947  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1948  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1949 
1950 
1951  // Raw yield extraction
1952  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1953  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
1954 
1955  // Cuts efficiency (from cuts variation)
1956  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1957  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1958 
1959  // PID efficiency (from PID/noPID)
1960  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1961  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1962 
1963 
1964  // MC dN/dpt (copied from D0 : will update later)
1965  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1966  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
1967  fMCPtShape->SetBinContent(1,0.03);
1968  fMCPtShape->SetBinContent(2,0.03);
1969 
1970  return;
1971 
1972 }
1973 
1974 //--------------------------------------------------------------------------
1976  //
1977  // D*+->D0pi syst errors. Responsible: A. Grelli
1978  // 2010 PbPb sample, 40-80 CC
1979  //
1980 
1981  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1982 
1983  // Normalization
1984  fNorm = new TH1F("fNorm","fNorm",24,0,24);
1985  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1986 
1987  // Branching ratio
1988  fBR = new TH1F("fBR","fBR",24,0,24);
1989  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1990 
1991  // Tracking efficiency
1992  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1993  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1994 
1995 
1996  // Raw yield extraction
1997  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1998  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
1999  fRawYield->SetBinContent(1,1);
2000  fRawYield->SetBinContent(2,0.15);
2001  fRawYield->SetBinContent(3,0.15);
2002  fRawYield->SetBinContent(4,0.15);
2003  fRawYield->SetBinContent(5,0.15);
2004  fRawYield->SetBinContent(6,0.10);
2005  fRawYield->SetBinContent(7,0.10);
2006  fRawYield->SetBinContent(8,0.10);
2007  fRawYield->SetBinContent(9,0.11);
2008  fRawYield->SetBinContent(10,0.11);
2009  fRawYield->SetBinContent(11,0.11);
2010  fRawYield->SetBinContent(12,0.11);
2011  fRawYield->SetBinContent(13,0.08);
2012  fRawYield->SetBinContent(14,0.08);
2013  fRawYield->SetBinContent(15,0.08);
2014  fRawYield->SetBinContent(16,0.08);
2015 
2016 
2017  // Cuts efficiency (from cuts variation)
2018  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2019  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
2020 
2021  // PID efficiency (from PID/noPID)
2022  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2023  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
2024 
2025 
2026  // MC dN/dpt (copied from D0 : will update later)
2027  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2028  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
2029  fMCPtShape->SetBinContent(2,0.05);
2030  fMCPtShape->SetBinContent(3,0.05);
2031  fMCPtShape->SetBinContent(4,0.05);
2032  fMCPtShape->SetBinContent(5,0.04);
2033  fMCPtShape->SetBinContent(6,0.02);
2034  fMCPtShape->SetBinContent(7,0.04);
2035  fMCPtShape->SetBinContent(8,0.04);
2036 
2037  return;
2038 
2039 }
2040 
2041 // ----------------------------- 2011 ----------------------------------
2043  //
2044  // D*+->D0pi syst errors. Responsible: A. Grelli
2045  // 2011 PbPb sample, 0-10% CC
2046  //
2047 
2048  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
2049 
2050 
2051  // Normalization
2052  fNorm = new TH1F("fNorm","fNorm",36,0,36);
2053  for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
2054 
2055  // Branching ratio
2056  fBR = new TH1F("fBR","fBR",36,0,36);
2057  for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
2058 
2059  // Tracking efficiency
2060  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
2061  for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
2062 
2063 
2064  // Raw yield extraction
2065  fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
2066  for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i, 0.);
2067  fRawYield->SetBinContent(4,0.10);
2068  fRawYield->SetBinContent(5,0.08);
2069  fRawYield->SetBinContent(6,0.05);
2070  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2071  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2072  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2073  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
2074  for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10); //4%
2075 
2076  // Cuts efficiency (from cuts variation)
2077  fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
2078  for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i, 0.);
2079  fCutsEff->SetBinContent(4,0.15);
2080  for(Int_t i=5;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
2081 
2082  // PID efficiency (from PID/noPID)
2083  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
2084  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
2085 
2086  // MC dN/dpt (from study on D* pt shape)
2087  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
2088  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i, 0.);
2089  fMCPtShape->SetBinContent(4,0.06);
2090  fMCPtShape->SetBinContent(5,0.03);
2091  for(Int_t i=5;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
2092 
2093  return;
2094 
2095 }
2096 
2097 //----------------------------------------------------------------------
2099  //
2100  // D*+->D0pi syst errors. Responsible: A. Grelli
2101  // 2011 PbPb sample, 30-50 CC
2102  //
2103 
2104  AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
2105 
2106  // Normalization
2107  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2108  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
2109 
2110  // Branching ratio
2111  fBR = new TH1F("fBR","fBR",24,0,24);
2112  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
2113 
2114  // Tracking efficiency
2115  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
2116  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
2117 
2118 
2119  // Raw yield extraction
2120  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2121  for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
2122  fRawYield->SetBinContent(1,0.12);
2123  fRawYield->SetBinContent(2,0.12);
2124  fRawYield->SetBinContent(3,0.12);
2125  fRawYield->SetBinContent(4,0.06);
2126  fRawYield->SetBinContent(5,0.06);
2127  fRawYield->SetBinContent(6,0.05);
2128  fRawYield->SetBinContent(7,0.08);
2129  fRawYield->SetBinContent(8,0.05);
2130  fRawYield->SetBinContent(9,0.05);
2131  fRawYield->SetBinContent(10,0.05);
2132  fRawYield->SetBinContent(11,0.05);
2133  fRawYield->SetBinContent(12,0.05);
2134  fRawYield->SetBinContent(13,0.05);
2135  fRawYield->SetBinContent(14,0.05);
2136  fRawYield->SetBinContent(15,0.05);
2137  fRawYield->SetBinContent(16,0.05);
2138 
2139 
2140  // Cuts efficiency (from cuts variation)
2141  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2142  for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
2143  fCutsEff->SetBinContent(1,0.16);
2144  fCutsEff->SetBinContent(2,0.16);
2145  fCutsEff->SetBinContent(3,0.16);
2146  fCutsEff->SetBinContent(4,0.10);
2147  fCutsEff->SetBinContent(5,0.10);
2148  fCutsEff->SetBinContent(6,0.10);
2149  // PID efficiency (from PID/noPID)
2150  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2151  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
2152 
2153 
2154  // MC dN/dpt (copied from D0 : will update later)
2155  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2156  for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
2157  fMCPtShape->SetBinContent(2,0.11);
2158  fMCPtShape->SetBinContent(3,0.11);
2159  fMCPtShape->SetBinContent(4,0.03);
2160  fMCPtShape->SetBinContent(5,0.02);
2161  fMCPtShape->SetBinContent(6,0.01);
2162  fMCPtShape->SetBinContent(7,0.01);
2163  fMCPtShape->SetBinContent(8,0.01);
2164 
2165  return;
2166 }
2167 //_________________________________________________________________________
2169  // Dstar syst in pPb 2013 MB
2170  //
2171  AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
2172 
2173  //Normalization
2174  fNorm = new TH1F("fNorm","fNorm",24,0,24);
2175  for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); //
2176 
2177  // Branching ratio
2178  fBR = new TH1F("fBR","fBR",24,0,24);
2179  for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
2180 
2181  // Tracking efficiency
2182  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
2183  for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
2184 
2185  // Raw yield extraction
2186  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2187  fRawYield->SetBinContent(1,1.0);
2188  fRawYield->SetBinContent(2,0.08);
2189  fRawYield->SetBinContent(3,0.02);
2190  fRawYield->SetBinContent(4,0.02);
2191  fRawYield->SetBinContent(5,0.05);
2192  fRawYield->SetBinContent(6,0.02);
2193  fRawYield->SetBinContent(7,0.04);
2194  fRawYield->SetBinContent(8,0.04);
2195  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2196  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
2197  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
2198 
2199  // Cuts efficiency (from cuts variation)
2200  fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2201  fCutsEff->SetBinContent(1,0.);
2202  fCutsEff->SetBinContent(2,0.10);
2203  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2204 
2205  // PID efficiency (from PID/noPID)
2206  fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2207  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2208 
2209  // MC dN/dpt
2210  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2211  for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2212  for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2213  for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
2214 
2215  // particle-antiparticle
2216  /*
2217  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2218  fPartAntipart->SetBinContent(1,1);
2219  fPartAntipart->SetBinContent(2,1);
2220  fPartAntipart->SetBinContent(3,0.12);
2221  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2222  */
2223 
2224  return;
2225 
2226 }
2227 
2228 //--------------------------------------------------------------------------
2230  // define errors for RAA vs. centrality
2232  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2233  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
2234  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2235 }
2236 //--------------------------------------------------------------------------
2238  // define errors for RAA vs. centrality
2240  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2241  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
2242  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2243 }
2244 //--------------------------------------------------------------------------
2246  // define errors for RAA vs. centrality
2248  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2249  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.14);
2250  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2251 }
2252 //--------------------------------------------------------------------------
2254  // define errors for RAA vs. centrality
2256  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2257  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.11);
2258  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2259 }
2260 //--------------------------------------------------------------------------
2262  // define errors for RAA vs. centrality
2264  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2265  for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
2266  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2267 }
2268 //--------------------------------------------------------------------------
2270  // define errors for RAA vs. centrality
2272  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
2273  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
2274 
2275 }
2276 //--------------------------------------------------------------------------
2278  // define errors for RAA vs. centrality
2280  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
2281  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2282 }
2283 //--------------------------------------------------------------------------
2285  // define errors for RAA vs. centrality
2287  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2288  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
2289 }
2290 //--------------------------------------------------------------------------
2292  // define errors for RAA vs. centrality
2294  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2295  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2296 }
2297 //--------------------------------------------------------------------------
2299  // define errors for RAA vs. centrality
2301  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
2302  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
2303 }
2304 
2305 //--------------------------------------------------------------------------
2307  // define errors for RAA vs. centrality
2309  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
2310  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
2311 }
2312 //--------------------------------------------------------------------------
2314  // define errors for RAA vs. centrality
2316  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2317  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
2318 }
2319 //--------------------------------------------------------------------------
2321  // define errors for RAA vs. centrality
2323  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2324  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2325 }
2326 //--------------------------------------------------------------------------
2328  // define errors for RAA vs. centrality
2330  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2331  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
2332 }
2333 //--------------------------------------------------------------------------
2335  // define errors for RAA vs. centrality
2337  for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2338  for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
2339 }
2340 
2341 //--------------------------------------------------------------------------
2343  // define errors for RAA vs. centrality 2011
2345 
2346 }
2347 //--------------------------------------------------------------------------
2349  // define errors for RAA vs. centrality 2011
2351 
2352 }
2353 //--------------------------------------------------------------------------
2355  // define errors for RAA vs. centrality 2011
2357 
2358 }
2359 //--------------------------------------------------------------------------
2361  // define errors for RAA vs. centrality 2011
2363 
2364 }
2365 //--------------------------------------------------------------------------
2367  // define errors for RAA vs. centrality 2011
2369 
2370 }
2371 //--------------------------------------------------------------------------
2373  // define errors for RAA vs. centrality 2011
2375 
2376 }
2377 
2378 //--------------------------------------------------------------------------
2380  // define errors for RAA vs. centrality 2011
2382 
2383  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.15);
2384  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2385  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.15);
2386 
2387  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2388  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2389  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2390 
2391  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
2392  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2393  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.10);
2394 }
2395 
2396 //--------------------------------------------------------------------------
2398  // define errors for RAA vs. centrality 2011
2400 
2401  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2402  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2403  for(Int_t i=9;i<=17;i++)fRawYield->SetBinContent(i,0.15);
2404 
2405  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2406  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2407  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2408 
2409  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
2410  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2411  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.15);
2412 }
2413 
2414 //--------------------------------------------------------------------------
2416  // define errors for RAA vs. centrality 2011
2418 
2419  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2420  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2421  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.08);
2422 
2423  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2424  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2425  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2426 
2427  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2428  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2429  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2430 }
2431 
2432 //--------------------------------------------------------------------------
2434  // define errors for RAA vs. centrality 2011
2436 
2437  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.20);
2438  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2439  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.20);
2440 
2441  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2442  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2443  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2444 
2445  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2446  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2447  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2448 }
2449 
2450 //--------------------------------------------------------------------------
2452  // define errors for RAA vs. centrality 2011
2454 
2455  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2456  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2457  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.05);
2458 
2459  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2460  for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2461  for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2462 
2463  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2464  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2465  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2466 }
2467 
2468 //--------------------------------------------------------------------------
2470  // define errors for RAA vs. centrality 2011
2472 
2473  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2474  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2475  for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.10);
2476 
2477  for(Int_t i=1;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2478  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20); // 10%
2479  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10); // 10%
2480  for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20); // 10%
2481 }
2482 
2483 //--------------------------------------------------------------------------
2485  // define errors for RAA vs. centrality 2011
2487 
2488 }
2489 //--------------------------------------------------------------------------
2491  // define errors for RAA vs. centrality 2011
2493 
2494 }
2495 //--------------------------------------------------------------------------
2497  // define errors for RAA vs. centrality 2011
2499 
2500 }
2501 //--------------------------------------------------------------------------
2503  // define errors for RAA vs. centrality 2011
2505 
2506 }
2507 //--------------------------------------------------------------------------
2509  // define errors for RAA vs. centrality 2011
2511 
2512 }
2513 //--------------------------------------------------------------------------
2515  // define errors for RAA vs. centrality 2011
2517 
2518 }
2519 
2520 //_________________________________________________________________________
2523 
2524  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2525 
2526  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2527  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.06);
2528  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2529 
2530  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2531 
2532  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2533  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2534  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2535 
2536  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2537  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2538  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2539 
2540 
2541 };
2542 //_________________________________________________________________________
2545 
2546  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2547 
2548  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2549  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.08);
2550  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2551 
2552  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2553 
2554  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2555  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2556  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2557 
2558  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2559  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2560  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2561 };
2562 //_________________________________________________________________________
2565 
2566  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2567 
2568  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2569  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.12);
2570  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2571 
2572  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2573 
2574  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2575  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2576  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2577 
2578  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2579  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2580  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2581 
2582 };
2583 //_________________________________________________________________________
2586 
2587  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2588 
2589  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2590  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2591  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2592 
2593  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2594 
2595  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2596  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2597  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2598 
2599  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2600  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2601  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2602 
2603 };
2604 //_________________________________________________________________________
2607 
2608  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2609 
2610  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2611  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2612  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2613 
2614  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2615 
2616  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2617  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2618  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2619 
2620  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2621  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2622  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2623 
2624 };
2625 
2626 //_________________________________________________________________________
2629 
2630  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2631 
2632  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2633  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2634  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2635 
2636  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2637 
2638  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2639  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2640  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2641 
2642  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2643  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.10);
2644  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.10);
2645 
2646 };
2647 //_________________________________________________________________________
2650  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2651 
2652  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2653  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2654  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2655 
2656  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2657 
2658  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2659  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2660  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2661 
2662  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2663  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2664  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2665 
2666 };
2667 //_________________________________________________________________________
2670  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2671 
2672  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2673  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2674  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2675 
2676  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2677 
2678  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2679  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2680  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2681 
2682  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2683  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2684  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2685 
2686 };
2687 //_________________________________________________________________________
2690  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2691 
2692  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2693  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2694  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2695 
2696  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2697 
2698  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2699  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2700  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2701 
2702  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2703  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2704  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2705 
2706 };
2707 //_________________________________________________________________________
2710 
2711  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2712 
2713  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2714  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2715  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2716 
2717  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2718 
2719  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2720  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2721  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2722 
2723  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2724  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2725  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2726 
2727 };
2728 
2729 
2730 //_________________________________________________________________________
2733 
2734  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2735 
2736  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2737  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2738  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2739 
2740  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2741 
2742  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2743  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2744  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2745 
2746  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2747  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2748  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2749 
2750 };
2751 //_________________________________________________________________________
2754 
2755  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2756 
2757  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2758  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2759  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2760 
2761  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2762 
2763  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2764  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2765  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2766 
2767  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2768  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2769  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2770 
2771 };
2772 //_________________________________________________________________________
2775 
2776  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2777 
2778  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2779  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2780  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2781 
2782  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2783 
2784  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2785  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2786  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2787 
2788  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2789  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2790  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2791 
2792 };
2793 //_________________________________________________________________________
2796 
2797  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2798 
2799  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2800  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2801  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2802 
2803  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2804 
2805  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2806  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2807  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2808 
2809  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2810  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2811  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2812 };
2813 //_________________________________________________________________________
2816 
2817  AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2818 
2819  for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2820  for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.07);
2821  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2822 
2823  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2824 
2825  for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2826  for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2827  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2828 
2829  for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2830  for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2831  for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2832 
2833 };
2834 
2838 
2839  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2840 
2841 
2842  // Raw yield extraction
2843  fRawYield->SetBinContent(1,0.0);
2844  fRawYield->SetBinContent(2,0.085);
2845  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
2846  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
2847  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2848  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2849  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2850  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2851 
2852  // Cuts efficiency (from cuts variation)
2853  fCutsEff->SetBinContent(1,0.0);
2854  fCutsEff->SetBinContent(2,0.10);
2855  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2856 
2857  // PID efficiency (from PID/noPID)
2858  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2859 
2860  fMCPtShape->SetBinContent(1,0.0);
2861  fMCPtShape->SetBinContent(2,0.028);
2862  fMCPtShape->SetBinContent(3,0.02);
2863  fMCPtShape->SetBinContent(4,0.02);
2864  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
2865 
2866 
2867 
2868 }
2869 
2870 //_________________________________________________________________________
2873 
2874  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2875 
2876 
2877  // Raw yield extraction
2878  // Raw yield extraction
2879  fRawYield->SetBinContent(1,0.0);
2880  fRawYield->SetBinContent(2,0.085);
2881  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
2882  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
2883  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2884  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2885  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2886  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2887 
2888  // Cuts efficiency (from cuts variation)
2889  fCutsEff->SetBinContent(1,0.0);
2890  fCutsEff->SetBinContent(2,0.10);
2891  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2892 
2893 
2894  // PID efficiency (from PID/noPID)
2895  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2896 
2897  fMCPtShape->SetBinContent(1,0.0);
2898  fMCPtShape->SetBinContent(2,0.028);
2899  fMCPtShape->SetBinContent(3,0.02);
2900  fMCPtShape->SetBinContent(4,0.02);
2901  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
2902 
2903 }
2904 
2905 //_________________________________________________________________________
2908 
2909  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2910 
2911 
2912  // Raw yield extraction
2913  fRawYield->SetBinContent(1,0.0);
2914  fRawYield->SetBinContent(2,0.085);
2915  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
2916  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
2917  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2918  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2919  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2920  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2921 
2922  // Cuts efficiency (from cuts variation)
2923  fCutsEff->SetBinContent(1,0.0);
2924  fCutsEff->SetBinContent(2,0.10);
2925  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2926 
2927 
2928  // PID efficiency (from PID/noPID)
2929  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2930 
2931  fMCPtShape->SetBinContent(1,0.0);
2932  fMCPtShape->SetBinContent(2,0.028);
2933  fMCPtShape->SetBinContent(3,0.02);
2934  fMCPtShape->SetBinContent(4,0.02);
2935  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
2936 
2937 
2938 }
2939 
2940 //_________________________________________________________________________
2943 
2944  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2945 
2946  // Raw yield extraction
2947  fRawYield->SetBinContent(1,0.0);
2948  fRawYield->SetBinContent(2,0.085);
2949  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
2950  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
2951  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2952  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2953  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2954  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2955 
2956 
2957  // Cuts efficiency (from cuts variation)
2958  fCutsEff->SetBinContent(1,0.0);
2959  fCutsEff->SetBinContent(2,0.10);
2960  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2961 
2962 
2963  // PID efficiency (from PID/noPID)
2964  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
2965 
2966  fMCPtShape->SetBinContent(1,0.0);
2967  fMCPtShape->SetBinContent(2,0.072);
2968  fMCPtShape->SetBinContent(3,0.06);
2969  fMCPtShape->SetBinContent(4,0.06);
2970  fMCPtShape->SetBinContent(5,0.03);
2971  fMCPtShape->SetBinContent(6,0.03);
2972  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.02);
2973 
2974 
2975 
2976 }
2977 
2978 //_________________________________________________________________________
2981 
2982  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
2983 
2984  // Raw yield extraction
2985  fRawYield->SetBinContent(1,0.0);
2986  fRawYield->SetBinContent(2,0.085);
2987  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
2988  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
2989  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2990  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2991  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2992  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2993 
2994  // Cuts efficiency (from cuts variation)
2995  fCutsEff->SetBinContent(1,0.0);
2996  fCutsEff->SetBinContent(2,0.10);
2997  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2998 
2999 
3000  // PID efficiency (from PID/noPID)
3001  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3002 
3003 
3004  fMCPtShape->SetBinContent(1,0.0);
3005  fMCPtShape->SetBinContent(2,0.053);
3006  fMCPtShape->SetBinContent(3,0.05);
3007  fMCPtShape->SetBinContent(4,0.05);
3008  fMCPtShape->SetBinContent(5,0.03);
3009  fMCPtShape->SetBinContent(6,0.03);
3010  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3011 
3012 
3013 }
3014 
3015 //_________________________________________________________________________
3018 
3019  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
3020 
3021  // Raw yield extraction
3022  fRawYield->SetBinContent(1,0.0);
3023  fRawYield->SetBinContent(2,0.085);
3024  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
3025  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
3026  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
3027  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
3028  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
3029  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
3030 
3031 
3032  // Cuts efficiency (from cuts variation)
3033  fCutsEff->SetBinContent(1,0.0);
3034  fCutsEff->SetBinContent(2,0.10);
3035  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3036 
3037  // PID efficiency (from PID/noPID)
3038  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3039 
3040 
3041  fMCPtShape->SetBinContent(1,0.0);
3042  fMCPtShape->SetBinContent(2,0.053);
3043  fMCPtShape->SetBinContent(3,0.05);
3044  fMCPtShape->SetBinContent(4,0.05);
3045  fMCPtShape->SetBinContent(5,0.03);
3046  fMCPtShape->SetBinContent(6,0.03);
3047  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3048 
3049 }
3050 
3051 //_________________________________________________________________________
3054 
3055  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
3056 
3057  // Raw yield extraction
3058  fRawYield->SetBinContent(1,0.0);
3059  fRawYield->SetBinContent(2,0.085);
3060  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
3061  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
3062  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
3063  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
3064  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
3065  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
3066 
3067  // Cuts efficiency (from cuts variation)
3068  fCutsEff->SetBinContent(1,0.0);
3069  fCutsEff->SetBinContent(2,0.10);
3070  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3071 
3072  // PID efficiency (from PID/noPID)
3073  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3074 
3075 
3076  fMCPtShape->SetBinContent(1,0.0);
3077  fMCPtShape->SetBinContent(2,0.053);
3078  fMCPtShape->SetBinContent(3,0.05);
3079  fMCPtShape->SetBinContent(4,0.05);
3080  fMCPtShape->SetBinContent(5,0.03);
3081  fMCPtShape->SetBinContent(6,0.03);
3082  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3083 
3084 
3085 
3086 }
3087 
3088 //_________________________________________________________________________
3091 
3092  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
3093 
3094  // Raw yield extraction
3095  fRawYield->SetBinContent(1,0.0);
3096  fRawYield->SetBinContent(2,0.085);
3097  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
3098  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
3099  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
3100  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
3101  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
3102  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
3103 
3104  // Cuts efficiency (from cuts variation)
3105  fCutsEff->SetBinContent(1,0.0);
3106  fCutsEff->SetBinContent(2,0.10);
3107  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3108 
3109  // PID efficiency (from PID/noPID)
3110  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3111 
3112  fMCPtShape->SetBinContent(1,0.0);
3113  fMCPtShape->SetBinContent(2,0.010);
3114  fMCPtShape->SetBinContent(3,0.010);
3115  fMCPtShape->SetBinContent(4,0.010);
3116  fMCPtShape->SetBinContent(5,0.05);
3117  fMCPtShape->SetBinContent(6,0.05);
3118  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3119 
3120 
3121 }
3122 
3123 //_________________________________________________________________________
3126 
3127  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
3128 
3129 
3130  // Raw yield extraction
3131  fRawYield->SetBinContent(1,0.0);
3132  fRawYield->SetBinContent(2,0.085);
3133  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
3134  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
3135  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
3136  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
3137  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
3138  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
3139 
3140  // Cuts efficiency (from cuts variation)
3141  fCutsEff->SetBinContent(1,0.0);
3142  fCutsEff->SetBinContent(2,0.10);
3143  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3144 
3145  // PID efficiency (from PID/noPID)
3146  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3147 
3148  fMCPtShape->SetBinContent(1,0.0);
3149  fMCPtShape->SetBinContent(2,0.028);
3150  fMCPtShape->SetBinContent(3,0.02);
3151  fMCPtShape->SetBinContent(4,0.02);
3152  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3153 
3154 
3155 }
3156 
3157 //_________________________________________________________________________
3160 
3161  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3162 
3163  // Raw yield extraction
3164  fRawYield->SetBinContent(1,0.0);
3165  fRawYield->SetBinContent(2,0.085);
3166  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
3167  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
3168  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
3169  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
3170  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
3171  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
3172 
3173 
3174  // Cuts efficiency (from cuts variation)
3175  fCutsEff->SetBinContent(1,0.0);
3176  fCutsEff->SetBinContent(2,0.10);
3177  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3178 
3179 
3180  // PID efficiency (from PID/noPID)
3181  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3182 
3183  fMCPtShape->SetBinContent(1,0.0);
3184  fMCPtShape->SetBinContent(2,0.028);
3185  fMCPtShape->SetBinContent(3,0.02);
3186  fMCPtShape->SetBinContent(4,0.02);
3187  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3188 }
3189 
3190 //_________________________________________________________________________
3193 
3194  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3195 
3196 
3197  // Raw yield extraction
3198  fRawYield->SetBinContent(1,0.0);
3199  fRawYield->SetBinContent(2,0.085);
3200  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
3201  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
3202  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
3203  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
3204  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
3205  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
3206 
3207  // Cuts efficiency (from cuts variation)
3208  fCutsEff->SetBinContent(1,0.0);
3209  fCutsEff->SetBinContent(2,0.10);
3210  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3211 
3212 
3213  // PID efficiency (from PID/noPID)
3214  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3215 
3216  fMCPtShape->SetBinContent(1,0.0);
3217  fMCPtShape->SetBinContent(2,0.028);
3218  fMCPtShape->SetBinContent(3,0.02);
3219  fMCPtShape->SetBinContent(4,0.02);
3220  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);}
3221 
3222 //_________________________________________________________________________
3225 
3226  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3227 
3228  // Raw yield extraction
3229  fRawYield->SetBinContent(1,0.0);
3230  fRawYield->SetBinContent(2,0.085);
3231  for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.042);
3232  for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.042);
3233  for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
3234  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
3235  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
3236  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
3237 
3238 
3239  // Cuts efficiency (from cuts variation)
3240  fCutsEff->SetBinContent(1,0.0);
3241  fCutsEff->SetBinContent(2,0.10);
3242  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3243 
3244 
3245  // PID efficiency (from PID/noPID)
3246  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3247 
3248  fMCPtShape->SetBinContent(1,0.0);
3249  fMCPtShape->SetBinContent(2,0.072);
3250  fMCPtShape->SetBinContent(3,0.06);
3251  fMCPtShape->SetBinContent(4,0.06);
3252  fMCPtShape->SetBinContent(5,0.03);
3253  fMCPtShape->SetBinContent(6,0.03);
3254  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.02);
3255 
3256 }
3257 
3258 //_________________________________________________________________________
3261 
3262  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3263 
3264  // Raw yield extraction
3265  fRawYield->SetBinContent(1,0.0);
3266  fRawYield->SetBinContent(2,0.10);
3267  fRawYield->SetBinContent(3,0.05);
3268  fRawYield->SetBinContent(4,0.05);
3269  fRawYield->SetBinContent(5,0.04);
3270  fRawYield->SetBinContent(6,0.04);
3271  fRawYield->SetBinContent(7,0.04);
3272  fRawYield->SetBinContent(8,0.04);
3273  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3274  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3275  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
3276 
3277  // Cuts efficiency (from cuts variation)
3278  fCutsEff->SetBinContent(1,0.0);
3279  fCutsEff->SetBinContent(2,0.10);
3280  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3281 
3282  // PID efficiency (from PID/noPID)
3283  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3284 
3285 
3286  fMCPtShape->SetBinContent(1,0.0);
3287  fMCPtShape->SetBinContent(2,0.028);
3288  fMCPtShape->SetBinContent(3,0.02);
3289  fMCPtShape->SetBinContent(4,0.02);
3290  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3291 
3292 
3293 
3294 }
3295 
3296 //_________________________________________________________________________
3299 
3300  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3301 
3302  // Raw yield extraction
3303  fRawYield->SetBinContent(1,0.0);
3304  fRawYield->SetBinContent(2,0.10);
3305  fRawYield->SetBinContent(3,0.05);
3306  fRawYield->SetBinContent(4,0.05);
3307  fRawYield->SetBinContent(5,0.04);
3308  fRawYield->SetBinContent(6,0.04);
3309  fRawYield->SetBinContent(7,0.04);
3310  fRawYield->SetBinContent(8,0.04);
3311  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3312  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3313  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3314 
3315 
3316  // Cuts efficiency (from cuts variation)
3317  fCutsEff->SetBinContent(1,0.0);
3318  fCutsEff->SetBinContent(2,0.10);
3319  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3320 
3321  // PID efficiency (from PID/noPID)
3322  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3323  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3324 
3325  fMCPtShape->SetBinContent(1,0.0);
3326  fMCPtShape->SetBinContent(2,0.028);
3327  fMCPtShape->SetBinContent(3,0.02);
3328  fMCPtShape->SetBinContent(4,0.02);
3329  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3330 
3331  // particle-antiparticle
3332  /*
3333  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3334  fPartAntipart->SetBinContent(1,1);
3335  fPartAntipart->SetBinContent(2,1);
3336  fPartAntipart->SetBinContent(3,0.12);
3337  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3338  */
3339 
3340 }
3341 
3342 //_________________________________________________________________________
3345 
3346  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3347 
3348  // Raw yield extraction
3349  fRawYield->SetBinContent(1,0.0);
3350  fRawYield->SetBinContent(2,0.10);
3351  fRawYield->SetBinContent(3,0.05);
3352  fRawYield->SetBinContent(4,0.05);
3353  fRawYield->SetBinContent(5,0.04);
3354  fRawYield->SetBinContent(6,0.04);
3355  fRawYield->SetBinContent(7,0.04);
3356  fRawYield->SetBinContent(8,0.04);
3357  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3358  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3359  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3360 
3361  // Cuts efficiency (from cuts variation)
3362  fCutsEff->SetBinContent(1,0.0);
3363  fCutsEff->SetBinContent(2,0.10);
3364  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3365 
3366 
3367  // PID efficiency (from PID/noPID)
3368  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3369  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3370 
3371  fMCPtShape->SetBinContent(1,0.0);
3372  fMCPtShape->SetBinContent(2,0.028);
3373  fMCPtShape->SetBinContent(3,0.02);
3374  fMCPtShape->SetBinContent(4,0.02);
3375  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3376 
3377  // particle-antiparticle
3378  /*
3379  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3380  fPartAntipart->SetBinContent(1,1);
3381  fPartAntipart->SetBinContent(2,1);
3382  fPartAntipart->SetBinContent(3,0.12);
3383  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3384  */
3385 
3386 
3387 }
3388 
3389 //_________________________________________________________________________
3392 
3393  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3394 
3395  // Raw yield extraction
3396  fRawYield->SetBinContent(1,0.0);
3397  fRawYield->SetBinContent(2,0.0);
3398  fRawYield->SetBinContent(3,0.05);
3399  fRawYield->SetBinContent(4,0.05);
3400  fRawYield->SetBinContent(5,0.04);
3401  fRawYield->SetBinContent(6,0.04);
3402  fRawYield->SetBinContent(7,0.04);
3403  fRawYield->SetBinContent(8,0.04);
3404  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3405  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3406  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3407 
3408 
3409  // Cuts efficiency (from cuts variation)
3410  fCutsEff->SetBinContent(1,0.0);
3411  fCutsEff->SetBinContent(2,0.10);
3412  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3413 
3414  // PID efficiency (from PID/noPID)
3415  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3416  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3417 
3418 
3419  fMCPtShape->SetBinContent(1,0.0);
3420  fMCPtShape->SetBinContent(2,0.072);
3421  fMCPtShape->SetBinContent(3,0.06);
3422  fMCPtShape->SetBinContent(4,0.06);
3423  fMCPtShape->SetBinContent(5,0.03);
3424  fMCPtShape->SetBinContent(6,0.03);
3425  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.02);
3426  // particle-antiparticle
3427  /*
3428  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3429  fPartAntipart->SetBinContent(1,1);
3430  fPartAntipart->SetBinContent(2,1);
3431  fPartAntipart->SetBinContent(3,0.12);
3432  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3433  */
3434 
3435 
3436 }
3437 
3438 //_________________________________________________________________________
3441 
3442  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3443 
3444  // Raw yield extraction
3445  fRawYield->SetBinContent(1,0.0);
3446  fRawYield->SetBinContent(2,0.10);
3447  fRawYield->SetBinContent(3,0.05);
3448  fRawYield->SetBinContent(4,0.05);
3449  fRawYield->SetBinContent(5,0.04);
3450  fRawYield->SetBinContent(6,0.04);
3451  fRawYield->SetBinContent(7,0.04);
3452  fRawYield->SetBinContent(8,0.04);
3453  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3454  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3455  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3456 
3457  // Cuts efficiency (from cuts variation)
3458  fCutsEff->SetBinContent(1,0.0);
3459  fCutsEff->SetBinContent(2,0.10);
3460  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3461 
3462  // PID efficiency (from PID/noPID)
3463  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3464  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
3465 
3466 
3467  fMCPtShape->SetBinContent(1,0.0);
3468  fMCPtShape->SetBinContent(2,0.053);
3469  fMCPtShape->SetBinContent(3,0.05);
3470  fMCPtShape->SetBinContent(4,0.05);
3471  fMCPtShape->SetBinContent(5,0.03);
3472  fMCPtShape->SetBinContent(6,0.03);
3473  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3474 
3475 
3476  // particle-antiparticle
3477  /*
3478  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3479  fPartAntipart->SetBinContent(1,1);
3480  fPartAntipart->SetBinContent(2,1);
3481  fPartAntipart->SetBinContent(3,0.12);
3482  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3483  */
3484 
3485 
3486 }
3487 
3488 //_________________________________________________________________________
3491 
3492  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3493 
3494  // Raw yield extraction
3495  fRawYield->SetBinContent(1,0.0);
3496  fRawYield->SetBinContent(2,0.10);
3497  fRawYield->SetBinContent(3,0.05);
3498  fRawYield->SetBinContent(4,0.05);
3499  fRawYield->SetBinContent(5,0.04);
3500  fRawYield->SetBinContent(6,0.04);
3501  fRawYield->SetBinContent(7,0.04);
3502  fRawYield->SetBinContent(8,0.04);
3503  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3504  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3505  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
3506 
3507  // Cuts efficiency (from cuts variation)
3508  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3509  fCutsEff->SetBinContent(1,0.0);
3510  fCutsEff->SetBinContent(2,0.10);
3511  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3512 
3513  // PID efficiency (from PID/noPID)
3514  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3515  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3516 
3517 
3518  fMCPtShape->SetBinContent(1,0.0);
3519  fMCPtShape->SetBinContent(2,0.053);
3520  fMCPtShape->SetBinContent(3,0.05);
3521  fMCPtShape->SetBinContent(4,0.05);
3522  fMCPtShape->SetBinContent(5,0.03);
3523  fMCPtShape->SetBinContent(6,0.03);
3524  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3525 
3526  // particle-antiparticle
3527  /*
3528  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3529  fPartAntipart->SetBinContent(1,1);
3530  fPartAntipart->SetBinContent(2,1);
3531  fPartAntipart->SetBinContent(3,0.12);
3532  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3533  */
3534 
3535 
3536 }
3537 
3538 //_________________________________________________________________________
3541 
3542  AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3543 
3544  // Raw yield extraction
3545  fRawYield->SetBinContent(1,0.0);
3546  fRawYield->SetBinContent(2,0.10);
3547  fRawYield->SetBinContent(3,0.05);
3548  fRawYield->SetBinContent(4,0.05);
3549  fRawYield->SetBinContent(5,0.04);
3550  fRawYield->SetBinContent(6,0.04);
3551  fRawYield->SetBinContent(7,0.04);
3552  fRawYield->SetBinContent(8,0.04);
3553  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3554  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3555  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3556 
3557 
3558  // Cuts efficiency (from cuts variation)
3559  // fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3560  fCutsEff->SetBinContent(1,0.0);
3561  fCutsEff->SetBinContent(2,0.10);
3562  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3563 
3564  // PID efficiency (from PID/noPID)
3565  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3566  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3567 
3568 
3569  fMCPtShape->SetBinContent(1,0.0);
3570  fMCPtShape->SetBinContent(2,0.053);
3571  fMCPtShape->SetBinContent(3,0.05);
3572  fMCPtShape->SetBinContent(4,0.05);
3573  fMCPtShape->SetBinContent(5,0.03);
3574  fMCPtShape->SetBinContent(6,0.03);
3575  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3576 
3577 
3578 
3579  // particle-antiparticle
3580  /*
3581  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3582  fPartAntipart->SetBinContent(1,1);
3583  fPartAntipart->SetBinContent(2,1);
3584  fPartAntipart->SetBinContent(3,0.12);
3585  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3586  */
3587 
3588 
3589 
3590 }
3591 
3592 //_________________________________________________________________________
3595 
3596  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3597 
3598  // Raw yield extraction
3599  fRawYield->SetBinContent(1,0.0);
3600  fRawYield->SetBinContent(2,0.00);
3601  fRawYield->SetBinContent(3,0.05);
3602  fRawYield->SetBinContent(4,0.05);
3603  fRawYield->SetBinContent(5,0.04);
3604  fRawYield->SetBinContent(6,0.04);
3605  fRawYield->SetBinContent(7,0.04);
3606  fRawYield->SetBinContent(8,0.04);
3607  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3608  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3609  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3610 
3611  // Cuts efficiency (from cuts variation)
3612  fCutsEff->SetBinContent(1,0.0);
3613  fCutsEff->SetBinContent(2,0.10);
3614  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3615 
3616  // PID efficiency (from PID/noPID)
3617  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3618  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3619 
3620  fMCPtShape->SetBinContent(1,0.0);
3621  fMCPtShape->SetBinContent(2,0.010);
3622  fMCPtShape->SetBinContent(3,0.010);
3623  fMCPtShape->SetBinContent(4,0.010);
3624  fMCPtShape->SetBinContent(5,0.05);
3625  fMCPtShape->SetBinContent(6,0.05);
3626  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
3627 
3628  // particle-antiparticle
3629  /*
3630  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3631  fPartAntipart->SetBinContent(1,1);
3632  fPartAntipart->SetBinContent(2,1);
3633  fPartAntipart->SetBinContent(3,0.12);
3634  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3635  */
3636 
3637 }
3638 
3639 
3640 //_________________________________________________________________________
3643 
3644  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
3645 
3646  // Raw yield extraction
3647  fRawYield->SetBinContent(1,0.0);
3648  fRawYield->SetBinContent(2,0.10);
3649  fRawYield->SetBinContent(3,0.05);
3650  fRawYield->SetBinContent(4,0.05);
3651  fRawYield->SetBinContent(5,0.04);
3652  fRawYield->SetBinContent(6,0.04);
3653  fRawYield->SetBinContent(7,0.04);
3654  fRawYield->SetBinContent(8,0.04);
3655  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3656  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3657  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
3658 
3659  // Cuts efficiency (from cuts variation)
3660  fCutsEff->SetBinContent(1,0.0);
3661  fCutsEff->SetBinContent(2,0.10);
3662  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3663 
3664  // PID efficiency (from PID/noPID)
3665  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3666 
3667  fMCPtShape->SetBinContent(1,0.0);
3668  fMCPtShape->SetBinContent(2,0.028);
3669  fMCPtShape->SetBinContent(3,0.02);
3670  fMCPtShape->SetBinContent(4,0.02);
3671  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3672 
3673 
3674 
3675 }
3676 
3677 //_________________________________________________________________________
3680 
3681  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
3682 
3683  // Raw yield extraction
3684  fRawYield->SetBinContent(1,0.0);
3685  fRawYield->SetBinContent(2,0.10);
3686  fRawYield->SetBinContent(3,0.05);
3687  fRawYield->SetBinContent(4,0.05);
3688  fRawYield->SetBinContent(5,0.04);
3689  fRawYield->SetBinContent(6,0.04);
3690  fRawYield->SetBinContent(7,0.04);
3691  fRawYield->SetBinContent(8,0.04);
3692  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3693  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3694  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3695 
3696  // Cuts efficiency (from cuts variation)
3697  fCutsEff->SetBinContent(1,0.0);
3698  fCutsEff->SetBinContent(2,0.10);
3699  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3700 
3701  // PID efficiency (from PID/noPID)
3702  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3703  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3704 
3705 
3706  fMCPtShape->SetBinContent(1,0.0);
3707  fMCPtShape->SetBinContent(2,0.028);
3708  fMCPtShape->SetBinContent(3,0.02);
3709  fMCPtShape->SetBinContent(4,0.02);
3710  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3711 
3712  // particle-antiparticle
3713  /*
3714  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3715  fPartAntipart->SetBinContent(1,1);
3716  fPartAntipart->SetBinContent(2,1);
3717  fPartAntipart->SetBinContent(3,0.12);
3718  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3719  */
3720 
3721 }
3722 
3723 //_________________________________________________________________________
3726 
3727  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
3728 
3729  // Raw yield extraction
3730  fRawYield->SetBinContent(1,0.0);
3731  fRawYield->SetBinContent(2,0.10);
3732  fRawYield->SetBinContent(3,0.05);
3733  fRawYield->SetBinContent(4,0.05);
3734  fRawYield->SetBinContent(5,0.04);
3735  fRawYield->SetBinContent(6,0.04);
3736  fRawYield->SetBinContent(7,0.04);
3737  fRawYield->SetBinContent(8,0.04);
3738  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3739  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3740  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
3741 
3742  // Cuts efficiency (from cuts variation)
3743  fCutsEff->SetBinContent(1,0.0);
3744  fCutsEff->SetBinContent(2,0.10);
3745  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3746 
3747  // PID efficiency (from PID/noPID)
3748  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3749  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3750 
3751  // MC dN/dpt
3752 
3753  fMCPtShape->SetBinContent(1,0.0);
3754  fMCPtShape->SetBinContent(2,0.028);
3755  fMCPtShape->SetBinContent(3,0.02);
3756  fMCPtShape->SetBinContent(4,0.02);
3757  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3758 
3759  // particle-antiparticle
3760  /*
3761  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3762  fPartAntipart->SetBinContent(1,1);
3763  fPartAntipart->SetBinContent(2,1);
3764  fPartAntipart->SetBinContent(3,0.12);
3765  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3766  */
3767 
3768 
3769 }
3770 
3771 //_________________________________________________________________________
3774 
3775  AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
3776 
3777  // Raw yield extraction
3778  fRawYield->SetBinContent(1,0.0);
3779  fRawYield->SetBinContent(2,0.00);
3780  fRawYield->SetBinContent(3,0.05);
3781  fRawYield->SetBinContent(4,0.05);
3782  fRawYield->SetBinContent(5,0.04);
3783  fRawYield->SetBinContent(6,0.04);
3784  fRawYield->SetBinContent(7,0.04);
3785  fRawYield->SetBinContent(8,0.04);
3786  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3787  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.00);
3788  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
3789 
3790  // Cuts efficiency (from cuts variation)
3791  fCutsEff->SetBinContent(1,0.0);
3792  fCutsEff->SetBinContent(2,0.10);
3793  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3794 
3795  // PID efficiency (from PID/noPID)
3796  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3797  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3798 
3799  fMCPtShape->SetBinContent(1,0.0);
3800  fMCPtShape->SetBinContent(2,0.072);
3801  fMCPtShape->SetBinContent(3,0.06);
3802  fMCPtShape->SetBinContent(4,0.06);
3803  fMCPtShape->SetBinContent(5,0.03);
3804  fMCPtShape->SetBinContent(6,0.03);
3805  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.02);
3806  // particle-antiparticle
3807  /*
3808  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3809  fPartAntipart->SetBinContent(1,1);
3810  fPartAntipart->SetBinContent(2,1);
3811  fPartAntipart->SetBinContent(3,0.12);
3812  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3813  */
3814 
3815 
3816 }
3817 
3818 //_________________________________________________________________________
3821 
3822  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3823 
3824  // Raw yield extraction
3825  fRawYield->SetBinContent(1,0.0);
3826  fRawYield->SetBinContent(2,0.08);
3827  fRawYield->SetBinContent(3,0.03);
3828  fRawYield->SetBinContent(4,0.03);
3829  fRawYield->SetBinContent(5,0.03);
3830  fRawYield->SetBinContent(6,0.03);
3831  fRawYield->SetBinContent(7,0.06);
3832  fRawYield->SetBinContent(8,0.06);
3833  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3834  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3835  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
3836 
3837  // Cuts efficiency (from cuts variation)
3838  fCutsEff->SetBinContent(1,0.0);
3839  fCutsEff->SetBinContent(2,0.10);
3840  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3841 
3842  //MC ptshape + multiplicty weights
3843  fMCPtShape->SetBinContent(1,0.0);
3844  fMCPtShape->SetBinContent(2,0.028);
3845  fMCPtShape->SetBinContent(3,0.02);
3846  fMCPtShape->SetBinContent(4,0.02);
3847  for(Int_t i=4;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3848 
3849  // PID efficiency (from PID/noPID)
3850  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3851 
3852 
3853 }
3854 
3855 //_________________________________________________________________________
3858 
3859  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3860 
3861  // Raw yield extraction
3862  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3863  fRawYield->SetBinContent(1,0.0);
3864  fRawYield->SetBinContent(2,0.08);
3865  fRawYield->SetBinContent(3,0.03);
3866  fRawYield->SetBinContent(4,0.03);
3867  fRawYield->SetBinContent(5,0.03);
3868  fRawYield->SetBinContent(6,0.03);
3869  fRawYield->SetBinContent(7,0.06);
3870  fRawYield->SetBinContent(8,0.06);
3871  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3872  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3873  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
3874 
3875  // Cuts efficiency (from cuts variation)
3876  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3877  fCutsEff->SetBinContent(1,0.0);
3878  fCutsEff->SetBinContent(2,0.10);
3879  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
3880 
3881  // PID efficiency (from PID/noPID)
3882  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3883  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3884 
3885  //MC ptshape + multiplicty weights
3886  fMCPtShape->SetBinContent(1,0.0);
3887  fMCPtShape->SetBinContent(2,0.028);
3888  fMCPtShape->SetBinContent(3,0.02);
3889  fMCPtShape->SetBinContent(4,0.02);
3890  for(Int_t i=4;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3891 
3892  // particle-antiparticle
3893  /*
3894  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3895  fPartAntipart->SetBinContent(1,1);
3896  fPartAntipart->SetBinContent(2,1);
3897  fPartAntipart->SetBinContent(3,0.12);
3898  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3899  */
3900 
3901 }
3902 
3903 //_________________________________________________________________________
3906 
3907  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3908 
3909  // Raw yield extraction
3910  // Raw yield extraction
3911  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3912  fRawYield->SetBinContent(1,0.0);
3913  fRawYield->SetBinContent(2,0.08);
3914  fRawYield->SetBinContent(3,0.03);
3915  fRawYield->SetBinContent(4,0.03);
3916  fRawYield->SetBinContent(5,0.03);
3917  fRawYield->SetBinContent(6,0.03);
3918  fRawYield->SetBinContent(7,0.06);
3919  fRawYield->SetBinContent(8,0.06);
3920  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3921  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3922  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
3923 
3924 
3925 
3926  // Cuts efficiency (from cuts variation)
3927  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3928  fCutsEff->SetBinContent(1,0.0);
3929  fCutsEff->SetBinContent(2,0.10);
3930  for(Int_t i=3;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
3931 
3932  // PID efficiency (from PID/noPID)
3933  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3934  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3935 
3936  //MC ptshape + multiplicty weights
3937  fMCPtShape->SetBinContent(1,0.0);
3938  fMCPtShape->SetBinContent(2,0.028);
3939  fMCPtShape->SetBinContent(3,0.02);
3940  fMCPtShape->SetBinContent(4,0.02);
3941  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
3942 
3943  // particle-antiparticle
3944  /*
3945  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3946  fPartAntipart->SetBinContent(1,1);
3947  fPartAntipart->SetBinContent(2,1);
3948  fPartAntipart->SetBinContent(3,0.12);
3949  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3950  */
3951 
3952 
3953 }
3954 
3955 //_________________________________________________________________________
3958 
3959  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3960 
3961  // Raw yield extraction
3962  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3963  fRawYield->SetBinContent(1,0.0);
3964  fRawYield->SetBinContent(2,0.08);
3965  fRawYield->SetBinContent(3,0.03);
3966  fRawYield->SetBinContent(4,0.03);
3967  fRawYield->SetBinContent(5,0.03);
3968  fRawYield->SetBinContent(6,0.03);
3969  fRawYield->SetBinContent(7,0.06);
3970  fRawYield->SetBinContent(8,0.06);
3971  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
3972  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
3973  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
3974 
3975 
3976  // Cuts efficiency (from cuts variation)
3977  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3978  fCutsEff->SetBinContent(1,0.0);
3979  fCutsEff->SetBinContent(2,0.10);
3980  for(Int_t i=3;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
3981 
3982  // PID efficiency (from PID/noPID)
3983  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3984  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3985 
3986  //MC ptshape + multiplicty weights
3987  fMCPtShape->SetBinContent(1,0.0);
3988  fMCPtShape->SetBinContent(2,0.072);
3989  fMCPtShape->SetBinContent(3,0.06);
3990  fMCPtShape->SetBinContent(4,0.06);
3991  fMCPtShape->SetBinContent(5,0.03);
3992  fMCPtShape->SetBinContent(6,0.03);
3993  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.02);
3994 
3995  // particle-antiparticle
3996  /*
3997  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3998  fPartAntipart->SetBinContent(1,1);
3999  fPartAntipart->SetBinContent(2,1);
4000  fPartAntipart->SetBinContent(3,0.12);
4001  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4002  */
4003 
4004 
4005 }
4006 
4007 //_________________________________________________________________________
4010 
4011  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
4012 
4013  // Raw yield extraction
4014  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4015  fRawYield->SetBinContent(1,0.);
4016  fRawYield->SetBinContent(2,0.08);
4017  fRawYield->SetBinContent(3,0.03);
4018  fRawYield->SetBinContent(4,0.03);
4019  fRawYield->SetBinContent(5,0.03);
4020  fRawYield->SetBinContent(6,0.03);
4021  fRawYield->SetBinContent(7,0.06);
4022  fRawYield->SetBinContent(8,0.06);
4023  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4024  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4025  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
4026 
4027  // Cuts efficiency (from cuts variation)
4028  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4029  fCutsEff->SetBinContent(1,0.0);
4030  fCutsEff->SetBinContent(2,0.10);
4031  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
4032 
4033  // PID efficiency (from PID/noPID)
4034  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4035  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
4036 
4037  fMCPtShape->SetBinContent(1,0.0);
4038  fMCPtShape->SetBinContent(2,0.053);
4039  fMCPtShape->SetBinContent(3,0.05);
4040  fMCPtShape->SetBinContent(4,0.05);
4041  fMCPtShape->SetBinContent(5,0.03);
4042  fMCPtShape->SetBinContent(6,0.03);
4043  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
4044  // particle-antiparticle
4045  /*
4046  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4047  fPartAntipart->SetBinContent(1,1);
4048  fPartAntipart->SetBinContent(2,1);
4049  fPartAntipart->SetBinContent(3,0.12);
4050  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4051  */
4052 
4053 
4054 }
4055 
4056 //_________________________________________________________________________
4059 
4060  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
4061 
4062  // Raw yield extraction
4063  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4064  fRawYield->SetBinContent(1,0.);
4065  fRawYield->SetBinContent(2,0.08);
4066  fRawYield->SetBinContent(3,0.03);
4067  fRawYield->SetBinContent(4,0.03);
4068  fRawYield->SetBinContent(5,0.03);
4069  fRawYield->SetBinContent(6,0.03);
4070  fRawYield->SetBinContent(7,0.06);
4071  fRawYield->SetBinContent(8,0.06);
4072  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4073  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4074  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
4075 
4076  // Cuts efficiency (from cuts variation)
4077  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4078  fCutsEff->SetBinContent(1,0.0);
4079  fCutsEff->SetBinContent(2,0.10);
4080  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
4081 
4082 
4083  // PID efficiency (from PID/noPID)
4084  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4085  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
4086 
4087  fMCPtShape->SetBinContent(1,0.0);
4088  fMCPtShape->SetBinContent(2,0.053);
4089  fMCPtShape->SetBinContent(3,0.05);
4090  fMCPtShape->SetBinContent(4,0.05);
4091  fMCPtShape->SetBinContent(5,0.03);
4092  fMCPtShape->SetBinContent(6,0.03);
4093  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
4094 
4095 
4096  // particle-antiparticle
4097  /*
4098  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4099  fPartAntipart->SetBinContent(1,1);
4100  fPartAntipart->SetBinContent(2,1);
4101  fPartAntipart->SetBinContent(3,0.12);
4102  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4103  */
4104 
4105 
4106 }
4107 
4108 //_________________________________________________________________________
4111 
4112  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
4113 
4114  // Raw yield extraction
4115  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4116  fRawYield->SetBinContent(1,0.);
4117  fRawYield->SetBinContent(2,0.08);
4118  fRawYield->SetBinContent(3,0.03);
4119  fRawYield->SetBinContent(4,0.03);
4120  fRawYield->SetBinContent(5,0.03);
4121  fRawYield->SetBinContent(6,0.03);
4122  fRawYield->SetBinContent(7,0.06);
4123  fRawYield->SetBinContent(8,0.06);
4124  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4125  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4126  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
4127 
4128 
4129  // Cuts efficiency (from cuts variation)
4130  // fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4131  fCutsEff->SetBinContent(1,0.0);
4132  fCutsEff->SetBinContent(2,0.10);
4133  for(Int_t i=3;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
4134 
4135 
4136  // PID efficiency (from PID/noPID)
4137  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4138  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
4139 
4140  fMCPtShape->SetBinContent(1,0.0);
4141  fMCPtShape->SetBinContent(2,0.053);
4142  fMCPtShape->SetBinContent(3,0.05);
4143  fMCPtShape->SetBinContent(4,0.05);
4144  fMCPtShape->SetBinContent(5,0.03);
4145  fMCPtShape->SetBinContent(6,0.03);
4146  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
4147 
4148  // particle-antiparticle
4149  /*
4150  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4151  fPartAntipart->SetBinContent(1,1);
4152  fPartAntipart->SetBinContent(2,1);
4153  fPartAntipart->SetBinContent(3,0.12);
4154  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4155  */
4156 
4157 
4158 
4159 }
4160 
4161 //_________________________________________________________________________
4164 
4165  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
4166 
4167  // Raw yield extraction
4168  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4169  fRawYield->SetBinContent(1,0.);
4170  fRawYield->SetBinContent(2,0.08);
4171  fRawYield->SetBinContent(3,0.03);
4172  fRawYield->SetBinContent(4,0.03);
4173  fRawYield->SetBinContent(5,0.03);
4174  fRawYield->SetBinContent(6,0.03);
4175  fRawYield->SetBinContent(7,0.06);
4176  fRawYield->SetBinContent(8,0.06);
4177  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4178  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4179  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
4180 
4181  // Cuts efficiency (from cuts variation)
4182  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4183  fCutsEff->SetBinContent(1,0.0);
4184  fCutsEff->SetBinContent(2,0.10);
4185  for(Int_t i=3;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
4186 
4187  // PID efficiency (from PID/noPID)
4188  // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4189  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
4190 
4191  fMCPtShape->SetBinContent(1,0.0);
4192  fMCPtShape->SetBinContent(2,0.010);
4193  fMCPtShape->SetBinContent(3,0.010);
4194  fMCPtShape->SetBinContent(4,0.010);
4195  fMCPtShape->SetBinContent(5,0.05);
4196  fMCPtShape->SetBinContent(6,0.05);
4197  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
4198 
4199  // particle-antiparticle
4200  /*
4201  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4202  fPartAntipart->SetBinContent(1,1);
4203  fPartAntipart->SetBinContent(2,1);
4204  fPartAntipart->SetBinContent(3,0.12);
4205  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4206  */
4207 
4208 }
4209 
4210 
4211 //_________________________________________________________________________
4214 
4215  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 0-20% mult CL1");
4216 
4217  // Raw yield extraction
4218  fRawYield->SetBinContent(1,0.0);
4219  fRawYield->SetBinContent(2,0.08);
4220  fRawYield->SetBinContent(3,0.03);
4221  fRawYield->SetBinContent(4,0.03);
4222  fRawYield->SetBinContent(5,0.03);
4223  fRawYield->SetBinContent(6,0.03);
4224  fRawYield->SetBinContent(7,0.06);
4225  fRawYield->SetBinContent(8,0.06);
4226  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4227  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4228  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
4229 
4230  // Cuts efficiency (from cuts variation)
4231  fCutsEff->SetBinContent(1,0.0);
4232  fCutsEff->SetBinContent(2,0.10);
4233  for(Int_t i=3;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
4234 
4235  fMCPtShape->SetBinContent(1,0.0);
4236  fMCPtShape->SetBinContent(2,0.028);
4237  fMCPtShape->SetBinContent(3,0.02);
4238  fMCPtShape->SetBinContent(4,0.02);
4239  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
4240 
4241  // PID efficiency (from PID/noPID)
4242  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
4243 
4244 
4245 }
4246 
4247 //_________________________________________________________________________
4250 
4251  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 20-40% mult CL1");
4252 
4253  // Raw yield extraction
4254  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4255  fRawYield->SetBinContent(1,0.0);
4256  fRawYield->SetBinContent(2,0.08);
4257  fRawYield->SetBinContent(3,0.03);
4258  fRawYield->SetBinContent(4,0.03);
4259  fRawYield->SetBinContent(5,0.03);
4260  fRawYield->SetBinContent(6,0.03);
4261  fRawYield->SetBinContent(7,0.06);
4262  fRawYield->SetBinContent(8,0.06);
4263  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4264  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4265  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
4266 
4267 
4268  // Cuts efficiency (from cuts variation)
4269  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4270  fCutsEff->SetBinContent(1,0.0);
4271  fCutsEff->SetBinContent(2,0.10);
4272  for(Int_t i=3;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
4273 
4274  // PID efficiency (from PID/noPID)
4275  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4276  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
4277 
4278  fMCPtShape->SetBinContent(1,0.0);
4279  fMCPtShape->SetBinContent(2,0.028);
4280  fMCPtShape->SetBinContent(3,0.02);
4281  fMCPtShape->SetBinContent(4,0.02);
4282  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
4283 
4284  // particle-antiparticle
4285  /*
4286  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4287  fPartAntipart->SetBinContent(1,1);
4288  fPartAntipart->SetBinContent(2,1);
4289  fPartAntipart->SetBinContent(3,0.12);
4290  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4291  */
4292 
4293 }
4294 
4295 //_________________________________________________________________________
4298 
4299  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 40-60% mult CL1");
4300 
4301  // Raw yield extraction
4302  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4303  fRawYield->SetBinContent(1,0.0);
4304  fRawYield->SetBinContent(2,0.08);
4305  fRawYield->SetBinContent(3,0.03);
4306  fRawYield->SetBinContent(4,0.03);
4307  fRawYield->SetBinContent(5,0.03);
4308  fRawYield->SetBinContent(6,0.03);
4309  fRawYield->SetBinContent(7,0.06);
4310  fRawYield->SetBinContent(8,0.06);
4311  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4312  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4313  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
4314 
4315  // Cuts efficiency (from cuts variation)
4316  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4317  fCutsEff->SetBinContent(1,0.0);
4318  fCutsEff->SetBinContent(2,0.10);
4319  for(Int_t i=3;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
4320 
4321  // PID efficiency (from PID/noPID)
4322  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4323  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
4324 
4325  fMCPtShape->SetBinContent(1,0.0);
4326  fMCPtShape->SetBinContent(2,0.028);
4327  fMCPtShape->SetBinContent(3,0.02);
4328  fMCPtShape->SetBinContent(4,0.02);
4329  for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
4330 
4331 
4332  // particle-antiparticle
4333  /*
4334  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4335  fPartAntipart->SetBinContent(1,1);
4336  fPartAntipart->SetBinContent(2,1);
4337  fPartAntipart->SetBinContent(3,0.12);
4338  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4339  */
4340 
4341 
4342 }
4343 
4344 //_________________________________________________________________________
4347 
4348  AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, 60-100% mult CL1");
4349 
4350  // Raw yield extraction
4351  //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4352  fRawYield->SetBinContent(1,0.0);
4353  fRawYield->SetBinContent(2,0.00);
4354  fRawYield->SetBinContent(3,0.03);
4355  fRawYield->SetBinContent(4,0.03);
4356  fRawYield->SetBinContent(5,0.03);
4357  fRawYield->SetBinContent(6,0.03);
4358  fRawYield->SetBinContent(7,0.06);
4359  fRawYield->SetBinContent(8,0.06);
4360  for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
4361  for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
4362  for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
4363 
4364 
4365  // Cuts efficiency (from cuts variation)
4366  //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4367  fCutsEff->SetBinContent(1,0.0);
4368  fCutsEff->SetBinContent(2,0.0);
4369  for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
4370 
4371  // PID efficiency (from PID/noPID)
4372  //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4373  for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
4374 
4375  //MC ptshape + multiplicty weights
4376  fMCPtShape->SetBinContent(1,0.0);
4377  fMCPtShape->SetBinContent(2,0.072);
4378  fMCPtShape->SetBinContent(3,0.06);
4379  fMCPtShape->SetBinContent(4,0.06);
4380  fMCPtShape->SetBinContent(5,0.03);
4381  fMCPtShape->SetBinContent(6,0.03);
4382  for(Int_t i=7;i<=16;i++) fMCPtShape->SetBinContent(i,0.02);
4383 
4384  // particle-antiparticle
4385  /*
4386  fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
4387  fPartAntipart->SetBinContent(1,1);
4388  fPartAntipart->SetBinContent(2,1);
4389  fPartAntipart->SetBinContent(3,0.12);
4390  for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
4391  */
4392 
4393 
4394 }
4395 
4396 //--------------------------------------------------------------------------
4398  //
4399  // Lc->pKpi syst errors. Responsible: R. Romita
4400  // 2010 pp sample
4401  //
4402 
4403  // Normalization
4404  fNorm = new TH1F("fNorm","fNorm",4,2,6);
4405  for(Int_t i=1;i<=4;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
4406 
4407  // Tracking efficiency
4408  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",4,2,6);
4409  for(Int_t i=1;i<=4;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
4410 
4411  // Raw yield extraction
4412  fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4413  fRawYield->SetBinContent(1,0.20);
4414  fRawYield->SetBinContent(2,0.15);
4415  fRawYield->SetBinContent(3,0.10);
4416  fRawYield->SetBinContent(4,0.15);
4417 
4418  fCutsEff = new TH1F("fCutsEff","fCutsEff",4,2,6);
4419  fCutsEff->SetBinContent(1,0.50);
4420  fCutsEff->SetBinContent(2,0.07);
4421  fCutsEff->SetBinContent(3,0.30);
4422  fCutsEff->SetBinContent(4,0.40);
4423 
4424  // PID efficiency (from PID/noPID)
4425  fPIDEff = new TH1F("fPIDEff","fPIDEff",4,2,4);
4426  fPIDEff->SetBinContent(1,0.13); // 15%
4427  fPIDEff->SetBinContent(2,0.30); // 15%
4428  fPIDEff->SetBinContent(3,0.5); // 15%
4429  fPIDEff->SetBinContent(4,0.30); // 15%
4430 
4431  // MC dN/dpt
4432  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
4433  for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
4434  for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
4435  for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
4436  for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
4437  for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
4438 }
4439 
4440 //--------------------------------------------------------------------------
4442  //
4443  // Lc->pK0S syst errors. Responsible: A. De Caro
4444  // 2010 pp sample
4445  // Only tool preparation: no physical mean for these values!
4446  //
4447 
4448  // Normalization
4449  fNorm = new TH1F("fNorm","fNorm",12,0.,12.);
4450  for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05); // 5% error on sigmaV0and // only a try
4451 
4452  // Tracking efficiency
4453  fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0.,12.);
4454  for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.1); // 10% (4% per track) // only a try
4455 
4456  // Raw yield extraction
4457  fRawYield = new TH1F("fRawYield","fRawYield",12,0.,12.);
4458  for(Int_t i=1;i<=2;i++) fRawYield->SetBinContent(i,1.); // only a try
4459  for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.1); // only a try
4460 
4461  fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0.,12.);
4462  for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.1); // only a try
4463 
4464  // PID efficiency (from PID/noPID)
4465  fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0.,12.);
4466  for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // only a try
4467 
4468  // MC dN/dpt
4469  fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0.,12.);
4470  for(Int_t i=1; i<=2;i++) fMCPtShape->SetBinContent(i,1.); // only a try
4471  for(Int_t i=3; i<=12;i++) fMCPtShape->SetBinContent(i,0.1); // only a try
4472 
4473 }
4474 
4475 //--------------------------------------------------------------------------
4476 Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
4477  //
4478  // Get error
4479  //
4480 
4481  Int_t bin=fCutsEff->FindBin(pt);
4482 
4483  return fCutsEff->GetBinContent(bin);
4484 }
4485 //--------------------------------------------------------------------------
4486 Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
4487  //
4488  // Get error
4489  //
4490 
4491  Int_t bin=fMCPtShape->FindBin(pt);
4492 
4493  return fMCPtShape->GetBinContent(bin);
4494 }
4495 //--------------------------------------------------------------------------
4496 Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
4497  //
4498  // Get error
4499  //
4500 
4501  Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
4502 
4503  return TMath::Sqrt(err);
4504 }
4505 //--------------------------------------------------------------------------
4506 Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
4507  //
4508  // Get error
4509  //
4510 
4511  Int_t bin=fPIDEff->FindBin(pt);
4512 
4513  return fPIDEff->GetBinContent(bin);
4514 }
4515 //--------------------------------------------------------------------------
4516 Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
4517  //
4518  // Get error
4519  //
4520 
4521  Int_t bin=fTrackingEff->FindBin(pt);
4522 
4523  return fTrackingEff->GetBinContent(bin);
4524 }
4525 //--------------------------------------------------------------------------
4526 Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
4527  //
4528  // Get error
4529  //
4530 
4531  Int_t bin=fRawYield->FindBin(pt);
4532 
4533  return fRawYield->GetBinContent(bin);
4534 }
4535 //--------------------------------------------------------------------------
4536 Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
4537  //
4538  // Get error
4539  //
4540 
4541  Int_t bin=fPartAntipart->FindBin(pt);
4542 
4543  return fPartAntipart->GetBinContent(bin);
4544 }
4545 //--------------------------------------------------------------------------
4546 Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
4547  //
4548  // Get total syst error (except norm. error)
4549  //
4550 
4551  Double_t err=0.;
4552 
4553  if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
4555  // if(fBR) err += GetBRErr()*GetBRErr();
4556  if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
4557  if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
4558  if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
4560 
4561  err += feeddownErr*feeddownErr;
4562 
4563  return TMath::Sqrt(err);
4564 }
4565 //---------------------------------------------------------------------------
4566 void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
4567  //
4568  // Draw errors
4569  //
4570  gStyle->SetOptStat(0);
4571 
4572  TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
4573  cSystErr->Range(0.20,-0.5,18.4,0.34);
4574  cSystErr->SetRightMargin(0.318);
4575  cSystErr->SetFillColor(0);
4576 
4577  TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
4578  hFrame->SetAxisRange(1.,35.9,"X");
4579  if(fIsLowPtAnalysis) hFrame->SetAxisRange(0.,24.,"X");
4580  hFrame->SetAxisRange(-0.5,0.5,"Y");
4581  hFrame->Draw();
4582 
4583  TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
4584  leg->SetTextSize(0.03601695);
4585  leg->SetFillStyle(0);
4586  leg->SetBorderSize(0);
4587 
4588  TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
4589  Int_t nbins = fNorm->GetNbinsX();
4590  TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
4591  for(Int_t i=1;i<=36;i++) {
4592  Double_t pt = hTotErr->GetBinCenter(i);
4593  Double_t ptwidth = hTotErr->GetBinWidth(i);
4594 
4595  if(grErrFeeddown) {
4596  Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
4597  Double_t toterryl=0., toterryh=0.;
4598  for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
4599  grErrFeeddown->GetPoint(j,x,y);
4600  errxh = grErrFeeddown->GetErrorXhigh(j);
4601  errxl = grErrFeeddown->GetErrorXlow(j);
4602  if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
4603  erryh = grErrFeeddown->GetErrorYhigh(j);
4604  erryl = grErrFeeddown->GetErrorYlow(j);
4605  }
4606  }
4607  if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
4608  else toterryl = GetTotalSystErr(pt);
4609  if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
4610  else toterryh = GetTotalSystErr(pt);
4611 
4612  hTotErr->SetBinContent(i,toterryh);
4613  gTotErr->SetPoint(i,pt,0.);
4614  gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
4615  }
4616  else {
4617  hTotErr->SetBinContent(i,GetTotalSystErr(pt));
4618  gTotErr->SetPoint(i,pt,0.);
4619  gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
4620  }
4621 
4622  }
4623  gTotErr->SetLineColor(kBlack);
4624  gTotErr->SetFillColor(kRed);
4625  gTotErr->SetFillStyle(3002);
4626  gTotErr->Draw("2");
4627  leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
4628 // hTotErr->SetLineColor(1);
4629 // hTotErr->SetLineWidth(3);
4630 // hTotErr->Draw("same");
4631 // leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
4632 
4633 
4634  fNorm->SetFillColor(1);
4635  fNorm->SetFillStyle(3002);
4636  //fNorm->Draw("same");
4637  //TH1F *hNormRefl = ReflectHisto(fNorm);
4638  //hNormRefl->Draw("same");
4639  Double_t norm = fNorm->GetBinContent(1)*100;
4640  leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
4641 
4642  if(grErrFeeddown) {
4643  grErrFeeddown->SetFillColor(kTeal-8);
4644  grErrFeeddown->SetFillStyle(3001);
4645  grErrFeeddown->Draw("2");
4646  leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
4647  }
4648  if(fTrackingEff) {
4649  fTrackingEff->SetFillColor(4);
4650  fTrackingEff->SetFillStyle(3006);
4651  fTrackingEff->Draw("same");
4652  TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
4653  hTrackingEffRefl->Draw("same");
4654  leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
4655  }
4656  if(fBR) {
4657  fBR->SetFillColor(6);
4658  fBR->SetFillStyle(3005);
4659  //fBR->SetFillStyle(3020);
4660  fBR->Draw("same");
4661  TH1F *hBRRefl = ReflectHisto(fBR);
4662  hBRRefl->Draw("same");
4663  leg->AddEntry(fBR,"Branching ratio","f");
4664  }
4665  if(fRawYield) {
4666  Int_t ci; // for color index setting
4667  ci = TColor::GetColor("#00cc00");
4668  fRawYield->SetLineColor(ci);
4669  // fRawYield->SetLineColor(3);
4670  fRawYield->SetLineWidth(3);
4671  fRawYield->Draw("same");
4672  TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
4673  hRawYieldRefl->Draw("same");
4674  leg->AddEntry(fRawYield,"Yield extraction","l");
4675  }
4676  if(fCutsEff) {
4677  fCutsEff->SetLineColor(4);
4678  fCutsEff->SetLineWidth(3);
4679  fCutsEff->Draw("same");
4680  TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
4681  hCutsEffRefl->Draw("same");
4682  leg->AddEntry(fCutsEff,"Cuts efficiency","l");
4683  }
4684  if(fPIDEff) {
4685  fPIDEff->SetLineColor(7);
4686  fPIDEff->SetLineWidth(3);
4687  fPIDEff->Draw("same");
4688  TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
4689  hPIDEffRefl->Draw("same");
4690  leg->AddEntry(fPIDEff,"PID efficiency","l");
4691  }
4692  if(fMCPtShape) {
4693  Int_t ci = TColor::GetColor("#9933ff");
4694  fMCPtShape->SetLineColor(ci);
4695  // fMCPtShape->SetLineColor(8);
4696  fMCPtShape->SetLineWidth(3);
4697  fMCPtShape->Draw("same");
4698  TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
4699  hMCPtShapeRefl->Draw("same");
4700  leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
4701  }
4702  if(fPartAntipart) {
4703  Int_t ci = TColor::GetColor("#ff6600");
4704  fPartAntipart->SetLineColor(ci);
4705  // fPartAntipart->SetLineColor(9);
4706  fPartAntipart->SetLineWidth(3);
4707  fPartAntipart->Draw("same");
4708  TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
4709  hPartAntipartRefl->Draw("same");
4710  leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
4711  }
4712 
4713 
4714  leg->Draw();
4715 
4716  cSystErr->SaveAs("RelativeSystematics.eps");
4717 
4718  return;
4719 }
4720 //-------------------------------------------------------------------------
4721 TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
4722  //
4723  // Clones and reflects histogram
4724  //
4725  TH1F *hout=(TH1F*)hin->Clone("hout");
4726  hout->Scale(-1.);
4727 
4728  return hout;
4729 }
4730 
4731 
4732 
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