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