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