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