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