]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGHF/vertexingHF/AliHFSystErr.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Update D0 pPb systematics (A.Festanti)
[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / AliHFSystErr.cxx
CommitLineData
ada850b0 1
0233abe6 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
27de2dfb 17/* $Id$ */
18
0233abe6 19/////////////////////////////////////////////////////////////
20//
21// Class to handle systematic errors for charm hadrons
22//
23// Usage:
5541b811 24// AliHFSystEff syst; // DECAY = 1 for D0, 2, for D+, 3 for D*
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*
0233abe6 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
33/////////////////////////////////////////////////////////////
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>
5ad8ad40 41#include <TColor.h>
0233abe6 42
5541b811 43#include "AliLog.h"
0233abe6 44#include "AliHFSystErr.h"
45
46
47ClassImp(AliHFSystErr)
48
49//--------------------------------------------------------------------------
50AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
51TNamed(name,title),
52fNorm(0),
53fRawYield(0),
54fTrackingEff(0),
55fBR(0),
56fCutsEff(0),
57fPIDEff(0),
58fMCPtShape(0),
5541b811 59fPartAntipart(0),
60fRunNumber(10),
61fCollisionType(0),
a8c3295d 62fCentralityClass("0100"),
44dd56b2 63fRapidityRange("0101"),
cede55f5 64fIsLowEnergy(false),
44dd56b2 65fIsCentScan(false),
66fIsRapidityScan(false)
0233abe6 67{
68 //
69 // Default Constructor
70 //
71}
5541b811 72
0233abe6 73//--------------------------------------------------------------------------
5541b811 74AliHFSystErr::~AliHFSystErr() {
75 //
76 // Default Destructor
0233abe6 77 //
a8c3295d 78 /*
5541b811 79
80 if(fNorm) { delete fNorm; fNorm=0; }
81 if(fRawYield) { delete fRawYield; fRawYield=0; }
82 if(fTrackingEff) { delete fTrackingEff; fTrackingEff=0; }
83 if(fBR) { delete fBR; fBR=0; }
84 if(fCutsEff) { delete fCutsEff; fCutsEff=0; }
85 if(fPIDEff) { delete fPIDEff; fPIDEff=0; }
86 if(fMCPtShape) { delete fMCPtShape; fMCPtShape=0; }
87 if(fPartAntipart) { delete fPartAntipart; fPartAntipart=0; }
a8c3295d 88 */
5541b811 89}
90
91//--------------------------------------------------------------------------
92void AliHFSystErr::Init(Int_t decay){
93 //
94 // Variables/histos initialization
0233abe6 95 //
96
9341a8fa 97 if ((fRunNumber>11) && fIsLowEnergy==false) {
27360a0e 98 AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
5541b811 99 }
100
0233abe6 101 switch(decay) {
102 case 1: // D0->Kpi
a8c3295d 103 if (fCollisionType==0) {
104 if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
105 else InitD0toKpi2010pp();
9341a8fa 106 }
107 else if (fCollisionType==1) {
3db225bf 108 if (fRunNumber == 10){
109 if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
110 else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
111 else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
112 else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
113 else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
114 else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
115 else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
116 else AliFatal("Not yet implemented");
117 }
9341a8fa 118 else if (fRunNumber == 11){
3db225bf 119 if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
120 else if (fCentralityClass=="3050InPlane") InitD0toKpi2011PbPb3050InPlane();
121 else if (fCentralityClass=="3050OutOfPlane") InitD0toKpi2011PbPb3050OutOfPlane();
122 else if (fCentralityClass=="010") InitD0toKpi2011PbPb010CentScan();
123 else if (fCentralityClass=="1020") InitD0toKpi2011PbPb1020CentScan();
124 else if (fCentralityClass=="2030") InitD0toKpi2011PbPb2030CentScan();
125 else if (fCentralityClass=="3040") InitD0toKpi2011PbPb3040CentScan();
126 else if (fCentralityClass=="4050") InitD0toKpi2011PbPb4050CentScan();
127 else if (fCentralityClass=="5080") InitD0toKpi2010PbPb5080CentScan();
128 else AliFatal("Not yet implemented");
129 }
9341a8fa 130 }
131 else if (fCollisionType==2) {
86164e35 132 if (fCentralityClass=="0100") InitD0toKpi2013pPb0100();
44dd56b2 133 if (fIsRapidityScan) {
134 if (fRapidityRange == "0804") InitD0toKpi2013pPb0100RapScan0804();
135 if (fRapidityRange == "0401") InitD0toKpi2013pPb0100RapScan0401();
136 if (fRapidityRange == "0101") InitD0toKpi2013pPb0100RapScan0101();
137 if (fRapidityRange == "0104") InitD0toKpi2013pPb0100RapScan0104();
138 if (fRapidityRange == "0408") InitD0toKpi2013pPb0100RapScan0408();
139 }
9341a8fa 140 }
141 else AliFatal("Not yet implemented");
0233abe6 142 break;
3db225bf 143
0233abe6 144 case 2: // D+->Kpipi
a8c3295d 145 if (fCollisionType==0) {
146 if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
147 else InitDplustoKpipi2010pp();
9341a8fa 148 }
149 else if (fCollisionType==1) {
3db225bf 150 if (fRunNumber == 10){
151 if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
152 else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
153 else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
154 else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
155 else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
156 else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
157 else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
158 else AliFatal("Not yet implemented");
159 }
160 if(fRunNumber == 11){
161 if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
162 else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010CentScan();
163 else if (fCentralityClass=="1020") InitDplustoKpipi2011PbPb1020CentScan();
164 else if (fCentralityClass=="2030") InitDplustoKpipi2011PbPb2030CentScan();
165 else if (fCentralityClass=="3040") InitDplustoKpipi2011PbPb3040CentScan();
166 else if (fCentralityClass=="4050") InitDplustoKpipi2011PbPb4050CentScan();
167 else if (fCentralityClass=="5080") InitDplustoKpipi2010PbPb5080CentScan();
168 else AliFatal("Not yet implemented");
169 }
9341a8fa 170 }
171 else if (fCollisionType==2) {
86164e35 172 if (fCentralityClass=="0100") InitDplustoKpipi2013pPb0100();
44dd56b2 173 if (fIsRapidityScan) {
174 if (fRapidityRange == "0804") InitDplustoKpipi2013pPb0100RapScan0804();
175 if (fRapidityRange == "0401") InitDplustoKpipi2013pPb0100RapScan0401();
176 if (fRapidityRange == "0101") InitDplustoKpipi2013pPb0100RapScan0101();
177 if (fRapidityRange == "0104") InitDplustoKpipi2013pPb0100RapScan0104();
178 if (fRapidityRange == "0408") InitDplustoKpipi2013pPb0100RapScan0408();
179 }
9341a8fa 180 }
181 else AliFatal("Not yet implemented");
0233abe6 182 break;
183 case 3: // D*->D0pi
fa2be30c 184 if (fCollisionType==0) {
185 if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
186 else InitDstartoD0pi2010pp();
9341a8fa 187 }
188 else if (fCollisionType==1) {
3db225bf 189 if (fRunNumber == 10){
190 if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
191 else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
192 else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
193 else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
194 else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
195 else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
196 else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
197 else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
198 else AliFatal("Not yet implemented");
199 }
200 if (fRunNumber == 11){
201 if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
202 else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010CentScan();
203 else if (fCentralityClass=="1020") InitDstartoD0pi2011PbPb1020CentScan();
204 else if (fCentralityClass=="2030") InitDstartoD0pi2011PbPb2030CentScan();
205 else if (fCentralityClass=="3040") InitDstartoD0pi2011PbPb3040CentScan();
206 else if (fCentralityClass=="4050") InitDstartoD0pi2011PbPb4050CentScan();
207 else if (fCentralityClass=="5080") InitDstartoD0pi2010PbPb5080CentScan();
208 else AliFatal("Not yet implemented");
209 }
210 }
211 else if (fCollisionType==2) {
86164e35 212 if (fCentralityClass=="0100") InitDstartoD0pi2013pPb0100();
44dd56b2 213 if (fIsRapidityScan) {
214 if (fRapidityRange == "0804") InitDstartoD0pi2013pPb0100RapScan0804();
215 if (fRapidityRange == "0401") InitDstartoD0pi2013pPb0100RapScan0401();
216 if (fRapidityRange == "0101") InitDstartoD0pi2013pPb0100RapScan0101();
217 if (fRapidityRange == "0104") InitDstartoD0pi2013pPb0100RapScan0104();
218 if (fRapidityRange == "0408") InitDstartoD0pi2013pPb0100RapScan0408();
219 }
5541b811 220 }
9341a8fa 221 else AliFatal("Not yet implemented");
0233abe6 222 break;
66247513 223 case 4: // D+s->KKpi
224 if (fCollisionType==0) InitDstoKKpi2010pp();
46988859 225 else if (fCollisionType==1) {
226 if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
227 else AliFatal("Not yet implemented");
9341a8fa 228 }
229 else if (fCollisionType==2) {
86164e35 230 if (fCentralityClass=="0100") InitDstoKKpi2013pPb0100();
46988859 231 }
27360a0e 232 else AliFatal("Not yet implemented");
66247513 233 break;
9ca30701 234 case 5: // Lc->pKpi
235 if (fCollisionType==0) InitLctopKpi2010pp();
236 else AliFatal("Not yet implemented");
237 break;
66247513 238
0233abe6 239 default:
240 printf("Invalid decay type: %d\n",decay);
241 break;
242 }
243
0233abe6 244}
5541b811 245
0233abe6 246//--------------------------------------------------------------------------
5541b811 247void AliHFSystErr::InitD0toKpi2010pp() {
248 //
249 // D0->Kpi syst errors. Responsible: A. Rossi
250 // 2010 pp sample
0233abe6 251 //
252
5541b811 253 // Normalization
f1d72175 254 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 255 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
5541b811 256
257 // Branching ratio
f1d72175 258 fBR = new TH1F("fBR","fBR",24,0,24);
259 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
5541b811 260
261 // Tracking efficiency
f1d72175 262 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
f8bc42d3 263 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
5541b811 264
265 // Raw yield extraction
f1d72175 266 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
5541b811 267 fRawYield->SetBinContent(1,1);
d2e5959f 268 fRawYield->SetBinContent(2,0.22);
269 fRawYield->SetBinContent(3,0.1);
270 for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
271 for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
f1d72175 272 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
273 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
0233abe6 274
5541b811 275 // Cuts efficiency (from cuts variation)
f1d72175 276 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
277 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
5541b811 278
279 // PID efficiency (from PID/noPID)
f1d72175 280 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
281 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
dc80fe6b 282 fPIDEff->SetBinContent(2,0.05); // 5%
5541b811 283
284 // MC dN/dpt
f1d72175 285 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 286 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
287 fMCPtShape->SetBinContent(1,0.03);
288 fMCPtShape->SetBinContent(2,0.03);
5541b811 289
290 // particle-antiparticle
f1d72175 291 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
292 // fPartAntipart->SetBinContent(1,1);
293 // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
5541b811 294
295 return;
0233abe6 296}
297//--------------------------------------------------------------------------
5541b811 298void AliHFSystErr::InitD0toKpi2010PbPb020() {
0233abe6 299 //
aebba8df 300 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 301 // 2010 PbPb sample, 0-20 CC
0233abe6 302 //
5541b811 303
304 // Normalization
305 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 306 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
5541b811 307
308 // Branching ratio
309 fBR = new TH1F("fBR","fBR",20,0,20);
310 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
311
312 // Tracking efficiency
313 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 314 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
5541b811 315
316 // Raw yield extraction
317 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
318 fRawYield->SetBinContent(1,0);
319 fRawYield->SetBinContent(2,0);
6c388dee 320 fRawYield->SetBinContent(3,0.08);
321 for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
322 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
5541b811 323
324 // Cuts efficiency (from cuts variation)
325 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 326 fCutsEff->SetBinContent(1,0.);
327 fCutsEff->SetBinContent(2,0.);
cc4bb301 328 fCutsEff->SetBinContent(3,0.13);
329 fCutsEff->SetBinContent(4,0.11);
330 for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
6c388dee 331 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
5541b811 332
333 // PID efficiency (from PID/noPID)
334 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
6c388dee 335 for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
5541b811 336
337 // MC dN/dpt
338 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
cc4bb301 339 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
340 fMCPtShape->SetBinContent(3,0.04);
341 fMCPtShape->SetBinContent(4,0.02);
b797b51b 342 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
578bfbba 343
6c388dee 344// // particle-antiparticle
345// fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
346// for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
347// fPartAntipart->SetBinContent(3,0.10);
348// fPartAntipart->SetBinContent(4,0.10);
349// fPartAntipart->SetBinContent(7,0.10);
350// fPartAntipart->SetBinContent(8,0.10);
0233abe6 351
5541b811 352 return;
353}
354//--------------------------------------------------------------------------
9a8b5080 355//--------------------------------------------------------------------------
356void AliHFSystErr::InitD0toKpi2011PbPb07half() {
357 //
358 // D0->Kpi syst errors. Responsible: A. Rossi
359 // 2011 PbPb sample, 0-7.5 CC
360 //
361
362 // Normalization
363 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 364 for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
9a8b5080 365
366 // Branching ratio
367 fBR = new TH1F("fBR","fBR",36,0,36);
368 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
369
370 // Tracking efficiency
371 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
372 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
373 for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
374
375 // Raw yield extraction
376 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
377 fRawYield->SetBinContent(1,0);
378 fRawYield->SetBinContent(2,0.1);
379 fRawYield->SetBinContent(3,0.05);
fd8ae434 380 for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
381 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
9a8b5080 382 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
383 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
384
385 // Cuts efficiency (from cuts variation)
386 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
387 fCutsEff->SetBinContent(1,0.);
388 fCutsEff->SetBinContent(2,0.15);
389 fCutsEff->SetBinContent(3,0.13);
390 fCutsEff->SetBinContent(4,0.11);
391 fCutsEff->SetBinContent(5,0.08);
392 for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
393 for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
394
395 // PID efficiency (from PID/noPID)
396 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
397 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
398 for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
399
400 // MC dN/dpt
401 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
402 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
403 fMCPtShape->SetBinContent(2,0.06);
404 fMCPtShape->SetBinContent(3,0.04);
405 fMCPtShape->SetBinContent(4,0.02);
406 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
407 for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
408 for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
409
410// // particle-antiparticle
411// fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
412// for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
413// fPartAntipart->SetBinContent(3,0.10);
414// fPartAntipart->SetBinContent(4,0.10);
415// fPartAntipart->SetBinContent(7,0.10);
416// fPartAntipart->SetBinContent(8,0.10);
417
418 return;
419}
93da8165 420
421//--------------------------------------------------------------------------
422void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
423 //
424 // D0->Kpi syst errors. Responsible: D. Caffarri
425 // 2011 PbPb sample, 30-50 CC InPlane
426 //
427 InitD0toKpi2011PbPb07half();
428 // Raw yield extraction
429 // fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
430 fRawYield->SetBinContent(1,0);
431 fRawYield->SetBinContent(2,0);
432 fRawYield->SetBinContent(3,0.05);
433 fRawYield->SetBinContent(4,0.08);
434 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
435 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
436 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
437 for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
438
439 // Cuts efficiency (from cuts variation)
440 //fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
441 fCutsEff->SetBinContent(1,0.);
442 fCutsEff->SetBinContent(2,0.0);
443 fCutsEff->SetBinContent(3,0.10);
444 fCutsEff->SetBinContent(4,0.10);
445 fCutsEff->SetBinContent(5,0.10);
446 fCutsEff->SetBinContent(6,0.10);
447 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
448 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
449 for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
450}
451
452//--------------------------------------------------------------------------
453void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
454 //
455 // D0->Kpi syst errors. Responsible: D. Caffarri
456 // 2011 PbPb sample, 30-50 CC OutOfPlane
457 //
458 InitD0toKpi2011PbPb07half();
459 // Raw yield extraction
460 //fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
461 fRawYield->SetBinContent(1,0);
462 fRawYield->SetBinContent(2,0.);
463 fRawYield->SetBinContent(3,0.05);
464 for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.07);
465 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
466 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
467 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
468 for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
469
470 // Cuts efficiency (from cuts variation)
471 // fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
472 fCutsEff->SetBinContent(1,0.);
473 fCutsEff->SetBinContent(2,0.);
474 fCutsEff->SetBinContent(3,0.10);
475 fCutsEff->SetBinContent(4,0.10);
476 fCutsEff->SetBinContent(5,0.10);
477 fCutsEff->SetBinContent(6,0.10);
478 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
479 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
480 for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
481}
482
9a8b5080 483//--------------------------------------------------------------------------
5541b811 484void AliHFSystErr::InitD0toKpi2010PbPb4080() {
93da8165 485 //
aebba8df 486 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 487 // 2010 PbPb sample, 40-80 CC
488 //
489
0233abe6 490 // Normalization
491 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 492 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
0233abe6 493
494 // Branching ratio
495 fBR = new TH1F("fBR","fBR",20,0,20);
844405bc 496 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
0233abe6 497
aebba8df 498 // Tracking efficiency
499 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 500 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
aebba8df 501
502
0233abe6 503 // Raw yield extraction
504 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
5541b811 505 fRawYield->SetBinContent(1,0);
506 fRawYield->SetBinContent(2,0);
6c388dee 507 for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
666cf569 508 //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
0233abe6 509
510 // Cuts efficiency (from cuts variation)
511 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 512 fCutsEff->SetBinContent(1,0.);
513 fCutsEff->SetBinContent(2,0.);
cc4bb301 514 fCutsEff->SetBinContent(3,0.13);
515 fCutsEff->SetBinContent(4,0.11);
516 for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
6c388dee 517 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
0233abe6 518
519 // PID efficiency (from PID/noPID)
520 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
6c388dee 521// for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
522// for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
523 for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
0233abe6 524
525 // MC dN/dpt
526 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
cc4bb301 527 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
528 // fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
529 // fMCPtShape->SetBinContent(4,0.02);
530 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
531
532 // // particle-antiparticle
533 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
534 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
0233abe6 535
536 return;
537}
a8c3295d 538
539//--------------------------------------------------------------------------
540void AliHFSystErr::InitD0toKpi2010ppLowEn() {
541 //
542 // D0->Kpi syst errors. Low energy run
543 // 2011 2.76 TeV pp sample
544 //
545 AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
546
547 // Normalization
548 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 549 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 4% error on sigmaV0and
a8c3295d 550
551 // Branching ratio
552 fBR = new TH1F("fBR","fBR",20,0,20);
553 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
554
555 // Tracking efficiency
556 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 557 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
a8c3295d 558
559 // Raw yield extraction
560 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
561 fRawYield->SetBinContent(1,1);
f9a85c82 562 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
a8c3295d 563
564 // Cuts efficiency (from cuts variation)
565 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
6c388dee 566 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
f9a85c82 567 fCutsEff->SetBinContent(2,0.20);
568 for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
569
a8c3295d 570
571 // PID efficiency (from PID/noPID)
572 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
f9a85c82 573 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
574 // fPIDEff->SetBinContent(2,0.20);
575 for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
a8c3295d 576
577 // MC dN/dpt
578 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 579 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
580 fMCPtShape->SetBinContent(1,0.03);
581 fMCPtShape->SetBinContent(2,0.03);
a8c3295d 582
f9a85c82 583// // particle-antiparticle
584// fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
585// fPartAntipart->SetBinContent(1,1);
586// fPartAntipart->SetBinContent(2,1);
587// for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
588// for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
a8c3295d 589
590 return;
591}
86164e35 592//_________________________________________________________________________
593void AliHFSystErr::InitD0toKpi2013pPb0100(){
bd332272 594 //
595 // D0->Kpi syst errors. p-Pb data sample
596 //
597 AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
86164e35 598
599 fNorm = new TH1F("fNorm","fNorm",24,0,24);
9341a8fa 600 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
86164e35 601
602 // Branching ratio
9341a8fa 603 fBR = new TH1F("fBR","fBR",24,0,24);
604 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.2% PDG2010
86164e35 605
606 // Tracking efficiency
607 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
608 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track
609
610 // Raw yield extraction
611 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
9341a8fa 612 fRawYield->SetBinContent(1,0.0);
e5f8509c
ZCV		 613 fRawYield->SetBinContent(2,0.07);
614 for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.02);
615 fRawYield->SetBinContent(7,0.03);
616 fRawYield->SetBinContent(8,0.05);
617 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.02);
618 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
619 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.1);
86164e35 620
621 // Cuts efficiency (from cuts variation)
622 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
9341a8fa 623 fCutsEff->SetBinContent(1,0.0);
e5f8509c 624 fCutsEff->SetBinContent(2,0.08);
9341a8fa 625 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
86164e35 626
627 // PID efficiency (from PID/noPID)
628 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
9341a8fa 629 fPIDEff->SetBinContent(1,0.0);
e5f8509c
ZCV		 630 fPIDEff->SetBinContent(2,0.1);
631 for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
86164e35 632
633 // MC dN/dpt
634 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
9341a8fa 635 fMCPtShape->SetBinContent(1,0.0);
636 fMCPtShape->SetBinContent(2,0.02);
e5f8509c 637 for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
86164e35 638
639 // particle-antiparticle
640 /*
641 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
9341a8fa 642 for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
86164e35 643 */
644 return;
645
646
647}
a8c3295d 648
44dd56b2 649
0233abe6 650//--------------------------------------------------------------------------
5541b811 651void AliHFSystErr::InitDplustoKpipi2010pp() {
0233abe6 652 //
653 // D+->Kpipi syst errors. Responsible: R. Bala
5541b811 654 // 2010 pp sample
0233abe6 655 //
656
1f91ef56 657
658// Normalization
659 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 660 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
626bcf57 661
662 // Branching ratio
1f91ef56 663 fBR = new TH1F("fBR","fBR",24,0,24);
f7c8a2f1 664 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
626bcf57 665
666 // Tracking efficiency
1f91ef56 667 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
f8bc42d3 668 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
626bcf57 669
670
671 // Raw yield extraction
1f91ef56 672 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
d2e5959f 673 fRawYield->SetBinContent(1,1);
1f91ef56 674 fRawYield->SetBinContent(2,0.25);
675 fRawYield->SetBinContent(3,0.25);
d2e5959f 676 fRawYield->SetBinContent(4,0.20);
1f91ef56 677 fRawYield->SetBinContent(5,0.09);
678 fRawYield->SetBinContent(6,0.09);
d2e5959f 679 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
680 for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
1f91ef56 681
626bcf57 682 // Cuts efficiency (from cuts variation)
1f91ef56 683 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
684 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
626bcf57 685
686 // PID efficiency (from PID/noPID)
1f91ef56 687 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
688 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
689 fPIDEff->SetBinContent(1,0.15); // 15%
690 fPIDEff->SetBinContent(2,0.15); // 15%
691 fPIDEff->SetBinContent(3,0.15); // 15%
692 fPIDEff->SetBinContent(4,0.15); // 15%
693 for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%
694
695 // MC dN/dpt
696 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
697 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
578bfbba 698 fMCPtShape->SetBinContent(1,0.03);
699 fMCPtShape->SetBinContent(2,0.03);
626bcf57 700
701
702 // particle-antiparticle
703 /*
704 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
705 fPartAntipart->SetBinContent(1,1);
706 fPartAntipart->SetBinContent(2,1);
707 fPartAntipart->SetBinContent(3,0.12);
708 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
709 */
0233abe6 710 return;
711}
1f91ef56 712
66247513 713//--------------------------------------------------------------------------
714void AliHFSystErr::InitDstoKKpi2010pp() {
715 //
716 // D+s->KKpi syst errors. Responsible: G.M. Innocenti
717 // 2010 pp sample
718 //
719
720
721// Normalization
722 fNorm = new TH1F("fNorm","fNorm",12,0,12);
494408cc 723 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
66247513 724
725 // Branching ratio
726 fBR = new TH1F("fBR","fBR",12,0,12);
494408cc 727 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
66247513 728
729 // Tracking efficiency
730 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
494408cc 731 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
66247513 732
733
734 // Raw yield extraction
735 fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
736 fRawYield->SetBinContent(1,1);
737 fRawYield->SetBinContent(2,1);
738 fRawYield->SetBinContent(3,0.20);
739 fRawYield->SetBinContent(4,0.20);
740 fRawYield->SetBinContent(5,0.15);
741 fRawYield->SetBinContent(6,0.15);
742 fRawYield->SetBinContent(7,0.15);
743 fRawYield->SetBinContent(8,0.15);
744 fRawYield->SetBinContent(9,0.20);
745 fRawYield->SetBinContent(10,0.20);
746 fRawYield->SetBinContent(11,0.20);
747 fRawYield->SetBinContent(12,0.20);
748
749 // Cuts efficiency (from cuts variation)
750 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
751 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%
752
753 // PID efficiency (from PID/noPID)
754 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
755 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
756
494408cc 757 // MC dN/dpt
758 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
759 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
760 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
761 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
762 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
763 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
66247513 764
765
766 // particle-antiparticle
767 /*
768 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
769 fPartAntipart->SetBinContent(1,1);
770 fPartAntipart->SetBinContent(2,1);
771 fPartAntipart->SetBinContent(3,0.12);
772 for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
773 */
774 return;
775}
776
777
436b185f 778//--------------------------------------------------------------------------
779void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
780 //
781 // D+->Kpipi syst errors. Responsible: E. Bruna
782 // 2011 PbPb sample, 0-7.5% CC
783 //
784
785 // Normalization
786 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 787 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
436b185f 788
789 // Branching ratio
790 fBR = new TH1F("fBR","fBR",36,0,36);
791 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
792
793 // Tracking efficiency
794 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
795 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
796
797 // Raw yield extraction
798 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
799 for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
eb4886c7 800 fRawYield->SetBinContent(4,0.30);
801 fRawYield->SetBinContent(5,0.20);
802 fRawYield->SetBinContent(6,0.20);
803 for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
804 for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
436b185f 805 for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
eb4886c7 806 for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
436b185f 807 for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
808
809
810 // Cuts efficiency (from cuts variation)
811 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
a7cfa9ed 812 for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
813 for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);
eb4886c7 814
436b185f 815
816 // PID efficiency (from PID/noPID)
817 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
818 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
819
a7cfa9ed 820 // MC dN/dpt (24/7/2012)
436b185f 821 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
3d446ac8 822 for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
823 for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
824 // for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
436b185f 825
826
827 // particle-antiparticle
828 /*
829 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
830 fPartAntipart->SetBinContent(1,1);
831 fPartAntipart->SetBinContent(2,1);
832 fPartAntipart->SetBinContent(3,0.12);
833 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
834 */
835
836 return;
837}
838
46988859 839//--------------------------------------------------------------------------
840void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
841 //
842 // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
843 // 2011 PbPb sample, 0-7.5% CC
844 //
845
846 // Normalization
847 fNorm = new TH1F("fNorm","fNorm",12,0,12);
a713683e 848 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
46988859 849
850 // Branching ratio
851 fBR = new TH1F("fBR","fBR",12,0,12);
d92cd593 852 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
46988859 853
854 // Tracking efficiency
855 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
856 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
857
858 // Raw yield extraction
859 fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
860 for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30);
861 for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
862
863
864
865 // Cuts efficiency (from cuts variation)
866 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
867 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
868
869 // PID efficiency (from PID/noPID)
870 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
871 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
872
873 // MC dN/dpt
874 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
875 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
876 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
877 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
878 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
879 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
880
881 // particle-antiparticle
882 /*
883 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
884 fPartAntipart->SetBinContent(1,1);
885 fPartAntipart->SetBinContent(2,1);
886 fPartAntipart->SetBinContent(3,0.12);
887 for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
888 */
889
890 return;
891}
892
86164e35 893//_________________________________________________________________________
894void AliHFSystErr::InitDstoKKpi2013pPb0100(){
bd332272 895 //
86164e35 896 // Dstar syst
bd332272 897 AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
86164e35 898
899 //Normalization
900 fNorm = new TH1F("fNorm","fNorm",24,0,24);
901 for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); //
902
903 // Branching ratio
904 fBR = new TH1F("fBR","fBR",12,0,12);
905 for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
906
907 // Tracking efficiency
908 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
59604f68 909 for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
86164e35 910
911 // Raw yield extraction
912 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
59604f68 913 for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
914 for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
915 for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
916 for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
917
86164e35 918 // Cuts efficiency (from cuts variation)
919 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
59604f68 920 for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //
86164e35 921
922 // PID efficiency (from PID/noPID)
923 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
59604f68 924 for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10);
925 for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
86164e35 926
927 // MC dN/dpt
928 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
59604f68 929 for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
86164e35 930
931 // particle-antiparticle
932 /*
933 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
934 fPartAntipart->SetBinContent(1,1);
935 fPartAntipart->SetBinContent(2,1);
936 fPartAntipart->SetBinContent(3,0.12);
937 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
938 */
939 return;
86164e35 940}
941
942
0233abe6 943//--------------------------------------------------------------------------
5541b811 944void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
945 //
946 // D+->Kpipi syst errors. Responsible: ??
947 // 2010 PbPb sample, 0-20 CC
948 //
aebba8df 949
aebba8df 950 // Normalization
951 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 952 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
aebba8df 953
954 // Branching ratio
955 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 956 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
aebba8df 957
958 // Tracking efficiency
959 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 960 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
aebba8df 961
962 // Raw yield extraction
963 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
6c388dee 964 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
965 // fRawYield->SetBinContent(5,0.23);
966 //fRawYield->SetBinContent(6,0.23);
cb142b90 967 fRawYield->SetBinContent(7,0.20);
968 fRawYield->SetBinContent(8,0.20);
969 fRawYield->SetBinContent(9,0.15);
970 fRawYield->SetBinContent(10,0.15);
971 fRawYield->SetBinContent(11,0.15);
972 fRawYield->SetBinContent(12,0.15);
aebba8df 973
974 // Cuts efficiency (from cuts variation)
975 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
666cf569 976 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
aebba8df 977
978 // PID efficiency (from PID/noPID)
979 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
980 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
aebba8df 981
f7c8a2f1 982 // MC dN/dpt (2/2/2012)
aebba8df 983 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 984 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
f7c8a2f1 985 for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
986 for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
987 for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
aebba8df 988
989
990 // particle-antiparticle
991 /*
992 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
993 fPartAntipart->SetBinContent(1,1);
994 fPartAntipart->SetBinContent(2,1);
995 fPartAntipart->SetBinContent(3,0.12);
996 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
997 */
6c388dee 998
5541b811 999 return;
1000}
1001
1002//--------------------------------------------------------------------------
1003void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
1004 //
1005 // D+->Kpipi syst errors. Responsible: ??
1006 // 2010 PbPb sample, 40-80 CC
1007 //
1008
aebba8df 1009
1010 // Normalization
1011 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 1012 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
aebba8df 1013
1014 // Branching ratio
1015 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1016 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
aebba8df 1017
1018 // Tracking efficiency
1019 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 1020 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
aebba8df 1021
1022
1023 // Raw yield extraction
1024 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1025 fRawYield->SetBinContent(1,1);
1026 fRawYield->SetBinContent(2,1);
1027 fRawYield->SetBinContent(3,1);
6970b225 1028 fRawYield->SetBinContent(4,0.15);
cb142b90 1029 fRawYield->SetBinContent(5,0.05);
1030 fRawYield->SetBinContent(6,0.05);
1031 fRawYield->SetBinContent(7,0.15);
1032 fRawYield->SetBinContent(8,0.15);
6970b225 1033 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
aebba8df 1034 for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
1035
1036 // Cuts efficiency (from cuts variation)
1037 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1038 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1039
1040 // PID efficiency (from PID/noPID)
1041 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1042 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1043 fPIDEff->SetBinContent(3,0.13); // 13%
1044
1045
f7c8a2f1 1046 // MC dN/dpt (2/2/2012)
aebba8df 1047 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 1048 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
f7c8a2f1 1049 for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1050 for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1051 for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
aebba8df 1052
1053
1054 // particle-antiparticle
1055 /*
1056 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1057 fPartAntipart->SetBinContent(1,1);
1058 fPartAntipart->SetBinContent(2,1);
1059 fPartAntipart->SetBinContent(3,0.12);
1060 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1061 */
5541b811 1062 return;
1063}
1064
a8c3295d 1065//--------------------------------------------------------------------------
1066void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
1067
1068 //
1069 // D+->Kpipi syst errors. Responsible: R. Bala
1070 // 2011 2.76 TeV pp sample
1071 //
1072 AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1073
1074 // Normalization
1075 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 1076 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
a8c3295d 1077
1078 // Branching ratio
1079 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1080 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
a8c3295d 1081
1082 // Tracking efficiency
1083 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 1084 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
a8c3295d 1085
1086 // Raw yield extraction
1087 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1088 fRawYield->SetBinContent(1,1);
1089 fRawYield->SetBinContent(2,1);
f9a85c82 1090 for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10); //5 to 10%
1091 fRawYield->SetBinContent(7,0.15);
1092 fRawYield->SetBinContent(8,0.15);
1093 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055); //5 to 10%
a8c3295d 1094
1095 // Cuts efficiency (from cuts variation)
1096 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
f9a85c82 1097 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
a8c3295d 1098
1099 // PID efficiency (from PID/noPID)
1100 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1101 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
f9a85c82 1102 fPIDEff->SetBinContent(3,0.10); // 13%
1103 fPIDEff->SetBinContent(4,0.10); // 13%
a8c3295d 1104
1105 // MC dN/dpt (copied from D0 : will update later)
1106 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 1107 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
578bfbba 1108 fMCPtShape->SetBinContent(1,0.03);
1109 fMCPtShape->SetBinContent(2,0.03);
a8c3295d 1110
1111 return;
1112}
1113
86164e35 1114//_________________________________________________________________________
1115void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
bd332272 1116 //
1117 // Dplus syst
1118 AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
86164e35 1119
1120 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1121 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.00); //
1122
2d325332 1123 fBR = new TH1F("fBR","fBR",20,0,24);
1124 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
86164e35 1125
1126 // Tracking efficiency
1127 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1128 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
4c1e4a80 1129
86164e35 1130 // Raw yield extraction
1131 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
4c1e4a80 1132 fRawYield->SetBinContent(1,1);
1133 fRawYield->SetBinContent(2,1);
1134 fRawYield->SetBinContent(3,0.06);
1135 fRawYield->SetBinContent(4,0.12);
1136 fRawYield->SetBinContent(5,0.05);
1137 fRawYield->SetBinContent(6,0.04);
1138 fRawYield->SetBinContent(7,0.04);
1139 fRawYield->SetBinContent(8,0.04);
1140 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.07);
1141 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.17);
1142 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.14);
1143
86164e35 1144 // Cuts efficiency (from cuts variation)
1145 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
4c1e4a80 1146 fCutsEff->SetBinContent(1,1);
1147 fCutsEff->SetBinContent(2,1);
1148 fCutsEff->SetBinContent(3,0.07);
1149 fCutsEff->SetBinContent(4,0.04);
1150 fCutsEff->SetBinContent(5,0.04);
1151 fCutsEff->SetBinContent(6,0.08);
1152 fCutsEff->SetBinContent(7,0.08);
1153 fCutsEff->SetBinContent(8,0.08);
1154 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.08);
1155 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
1156 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.1);
86164e35 1157
1158 // PID efficiency (from PID/noPID)
1159 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4c1e4a80 1160 fPIDEff->SetBinContent(1,1);
1161 fPIDEff->SetBinContent(2,1);
1162 for(Int_t i=3;i<=12;i++) fPIDEff->SetBinContent(i,0.0); //
294292d5 1163 for(Int_t i=13;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
86164e35 1164
4c1e4a80 1165 // MC dN/dpt
86164e35 1166 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
4c1e4a80 1167 for(Int_t i=1;i<=12;i++) fMCPtShape->SetBinContent(i,0.02);
1168 for(Int_t i=13;i<=24;i++) fMCPtShape->SetBinContent(i,0.04);
86164e35 1169
1170 // particle-antiparticle
1171 /*
1172 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1173 fPartAntipart->SetBinContent(1,1);
1174 fPartAntipart->SetBinContent(2,1);
1175 fPartAntipart->SetBinContent(3,0.12);
1176 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1177 */
1178 return;
1179
1180
1181}
1182
5541b811 1183//--------------------------------------------------------------------------
1184void AliHFSystErr::InitDstartoD0pi2010pp() {
0233abe6 1185 //
2cf98d0d 1186 // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
5541b811 1187 // 2010 pp sample
0233abe6 1188 //
1189
2cf98d0d 1190 // Normalization
fa2be30c 1191 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 1192 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
fa2be30c 1193
1194 // Branching ratio
1195 fBR = new TH1F("fBR","fBR",24,0,24);
1196 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1197
1198 // Tracking efficiency
1199 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1200 fTrackingEff->SetBinContent(1,1.0);
f8bc42d3 1201 fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
1202 fTrackingEff->SetBinContent(3,0.12);
1203 fTrackingEff->SetBinContent(3,0.12);
1204 for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
fa2be30c 1205
1206
1207 // Raw yield extraction
1208 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1209 fRawYield->SetBinContent(1,1.0);
d2e5959f 1210 fRawYield->SetBinContent(2,0.10);
fa2be30c 1211 fRawYield->SetBinContent(3,0.04);
d2e5959f 1212 fRawYield->SetBinContent(4,0.03);
fa2be30c 1213 fRawYield->SetBinContent(5,0.03);
d2e5959f 1214 fRawYield->SetBinContent(6,0.05);
1215 fRawYield->SetBinContent(7,0.05);
1216 fRawYield->SetBinContent(8,0.05);
1217 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04); //4%
1218 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09); //4%
1219 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
fa2be30c 1220
1221 // Cuts efficiency (from cuts variation)
1222 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1223 fCutsEff->SetBinContent(2,0.22);
1224 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1225
1226 // PID efficiency (from PID/noPID)
1227 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1228 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1229
1230
1231 // MC dN/dpt (copied from D0 : will update later)
1232 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 1233 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1234 fMCPtShape->SetBinContent(1,0.03);
1235 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 1236
1237 return;
1238
1239
1240}
1241//--------------------------------------------------------------------------
1242void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
1243
1244 //
1245 // D+->Kpipi syst errors. Responsible: A. Grelli
1246 // 2011 2.76 TeV pp sample
1247 //
1248 AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1249
1250// Normalization
2cf98d0d 1251 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 1252 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
2cf98d0d 1253
1254 // Branching ratio
1255 fBR = new TH1F("fBR","fBR",20,0,20);
8accac87 1256 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
2cf98d0d 1257
1258 // Tracking efficiency
1259 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 1260 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
2cf98d0d 1261
1262 // Raw yield extraction
1263 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
8accac87 1264 fRawYield->SetBinContent(1,1);
fa2be30c 1265 fRawYield->SetBinContent(2,1);
f9a85c82 1266 fRawYield->SetBinContent(3,0.14);
1267 fRawYield->SetBinContent(4,0.14);
1268 fRawYield->SetBinContent(5,0.12);
1269 fRawYield->SetBinContent(6,0.12);
1270 fRawYield->SetBinContent(7,0.06);
1271 fRawYield->SetBinContent(8,0.06);
1272 fRawYield->SetBinContent(9,0.08);
1273 fRawYield->SetBinContent(10,0.08);
1274 fRawYield->SetBinContent(11,0.08);
1275 fRawYield->SetBinContent(12,0.08);
fa2be30c 1276 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
2cf98d0d 1277
1278 // Cuts efficiency (from cuts variation)
1279 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
f9a85c82 1280 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
1281 fCutsEff->SetBinContent(3,0.15);
1282 fCutsEff->SetBinContent(4,0.15);
1283 fCutsEff->SetBinContent(5,0.15);
1284 fCutsEff->SetBinContent(6,0.15);
1285 fCutsEff->SetBinContent(7,0.10);
1286 fCutsEff->SetBinContent(8,0.10);
1287 fCutsEff->SetBinContent(9,0.10);
1288 fCutsEff->SetBinContent(10,0.10);
1289 fCutsEff->SetBinContent(11,0.10);
1290 fCutsEff->SetBinContent(12,0.10);
2cf98d0d 1291
1292 // PID efficiency (from PID/noPID)
1293 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
fa2be30c 1294 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
2cf98d0d 1295
fa2be30c 1296 // MC dN/dpt
2cf98d0d 1297 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 1298 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
578bfbba 1299 fMCPtShape->SetBinContent(1,0.03);
1300 fMCPtShape->SetBinContent(2,0.03);
2cf98d0d 1301
0233abe6 1302 return;
1303}
6c388dee 1304
fa2be30c 1305//------------------------------------------------------------------------
5541b811 1306void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
1307 //
fa2be30c 1308 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 1309 // 2010 PbPb sample, 0-20 CC
1310 //
1311
fa2be30c 1312 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1313
1314 // Normalization
1315 fNorm = new TH1F("fNorm","fNorm",24,0,24);
d5e15690 1316 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
fa2be30c 1317
1318 // Branching ratio
1319 fBR = new TH1F("fBR","fBR",24,0,24);
1320 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1321
1322 // Tracking efficiency
1323 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1324 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1325
1326
1327 // Raw yield extraction
1328 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
6c388dee 1329 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
1330 fRawYield->SetBinContent(3,0.2);
1331 fRawYield->SetBinContent(4,0.2);
1332 fRawYield->SetBinContent(5,0.2);
1333 fRawYield->SetBinContent(6,0.2);
fa2be30c 1334
1335 // Cuts efficiency (from cuts variation)
1336 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1337 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
6fd587a4 1338 fCutsEff->SetBinContent(4,0.15);
1339 fCutsEff->SetBinContent(5,0.15);
1340 fCutsEff->SetBinContent(6,0.15);
fa2be30c 1341
1342 // PID efficiency (from PID/noPID)
1343 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
6c388dee 1344 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
fa2be30c 1345
1346
6fd587a4 1347 // MC dN/dpt (from study on D* pt shape)
fa2be30c 1348 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
6fd587a4 1349 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
1350 fMCPtShape->SetBinContent(4,0.025);
1351 fMCPtShape->SetBinContent(5,0.025);
1352 fMCPtShape->SetBinContent(6,0.025);
1353 fMCPtShape->SetBinContent(7,0.04);
1354 fMCPtShape->SetBinContent(8,0.04);
1355 fMCPtShape->SetBinContent(9,0.03);
1356 fMCPtShape->SetBinContent(10,0.03);
1357 fMCPtShape->SetBinContent(11,0.03);
1358 fMCPtShape->SetBinContent(12,0.03);
1359
9a8b5080 1360 return;
1361
1362}
1363// ----------------------------- 2011 ----------------------------------
1364void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
1365 //
1366 // D*+->D0pi syst errors. Responsible: A. Grelli
1367 // 2011 PbPb sample, 0-7.5 CC
1368 //
1369
1370 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1371
1372
1373 // Normalization
1374 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 1375 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
9a8b5080 1376
1377 // Branching ratio
1378 fBR = new TH1F("fBR","fBR",36,0,36);
1379 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1380
1381 // Tracking efficiency
1382 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1383 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1384
1385
1386 // Raw yield extraction
1387 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1388 for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05); //4%
ada850b0 1389 fRawYield->SetBinContent(4,0.2);
9a8b5080 1390 fRawYield->SetBinContent(5,0.10);
ada850b0 1391 fRawYield->SetBinContent(6,0.10);
1392 fRawYield->SetBinContent(7,0.08);
9a8b5080 1393 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15); //4%
1394
9a8b5080 1395 // Cuts efficiency (from cuts variation)
1396 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1397 for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
9a8b5080 1398
1399 // PID efficiency (from PID/noPID)
1400 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1401 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
ada850b0 1402 fPIDEff->SetBinContent(4,0.09);
9a8b5080 1403
1404 // MC dN/dpt (from study on D* pt shape)
1405 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
273a7638 1406 for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
1407 fMCPtShape->SetBinContent(4,0.015);
1408 fMCPtShape->SetBinContent(5,0.015);
1409 fMCPtShape->SetBinContent(6,0.015);
1410 fMCPtShape->SetBinContent(7,0.02);
1411 fMCPtShape->SetBinContent(8,0.02);
9a8b5080 1412 fMCPtShape->SetBinContent(9,0.03);
1413 fMCPtShape->SetBinContent(10,0.03);
1414 fMCPtShape->SetBinContent(11,0.03);
1415 fMCPtShape->SetBinContent(12,0.03);
1416
6fd587a4 1417
fa2be30c 1418
1419 return;
1420
1421}
1422//-------------------------------------------------------------------------
1423void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
1424 //
1425 // D*+->D0pi syst errors. Responsible: A. Grelli
1426 // 2010 PbPb sample, 20-40 CC
1427 //
1428
1429 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1430
1431 // Normalization
1432 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1433 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
1434
1435 // Branching ratio
1436 fBR = new TH1F("fBR","fBR",24,0,24);
1437 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1438
1439 // Tracking efficiency
1440 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1441 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1442
1443
1444 // Raw yield extraction
1445 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1446 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
1447
1448 // Cuts efficiency (from cuts variation)
1449 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1450 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1451
1452 // PID efficiency (from PID/noPID)
1453 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1454 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1455
1456
1457 // MC dN/dpt (copied from D0 : will update later)
1458 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 1459 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
1460 fMCPtShape->SetBinContent(1,0.03);
1461 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 1462
5541b811 1463 return;
fa2be30c 1464
5541b811 1465}
6c388dee 1466
5541b811 1467//--------------------------------------------------------------------------
1468void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
1469 //
fa2be30c 1470 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 1471 // 2010 PbPb sample, 40-80 CC
1472 //
1473
fa2be30c 1474 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1475
1476 // Normalization
1477 fNorm = new TH1F("fNorm","fNorm",24,0,24);
d5e15690 1478 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
fa2be30c 1479
1480 // Branching ratio
1481 fBR = new TH1F("fBR","fBR",24,0,24);
1482 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1483
1484 // Tracking efficiency
1485 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1486 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1487
1488
1489 // Raw yield extraction
1490 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
6c388dee 1491 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
1492 fRawYield->SetBinContent(1,1);
1493 fRawYield->SetBinContent(2,0.15);
1494 fRawYield->SetBinContent(3,0.15);
1495 fRawYield->SetBinContent(4,0.15);
1496 fRawYield->SetBinContent(5,0.15);
1497 fRawYield->SetBinContent(6,0.10);
1498 fRawYield->SetBinContent(7,0.10);
1499 fRawYield->SetBinContent(8,0.10);
1500 fRawYield->SetBinContent(9,0.11);
1501 fRawYield->SetBinContent(10,0.11);
1502 fRawYield->SetBinContent(11,0.11);
1503 fRawYield->SetBinContent(12,0.11);
1504 fRawYield->SetBinContent(13,0.08);
1505 fRawYield->SetBinContent(14,0.08);
1506 fRawYield->SetBinContent(15,0.08);
1507 fRawYield->SetBinContent(16,0.08);
1508
1509
fa2be30c 1510 // Cuts efficiency (from cuts variation)
1511 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1512 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1513
1514 // PID efficiency (from PID/noPID)
1515 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
6c388dee 1516 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
fa2be30c 1517
1518
1519 // MC dN/dpt (copied from D0 : will update later)
1520 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
6fd587a4 1521 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1522 fMCPtShape->SetBinContent(2,0.05);
1523 fMCPtShape->SetBinContent(3,0.05);
1524 fMCPtShape->SetBinContent(4,0.05);
1525 fMCPtShape->SetBinContent(5,0.04);
1526 fMCPtShape->SetBinContent(6,0.02);
1527 fMCPtShape->SetBinContent(7,0.04);
1528 fMCPtShape->SetBinContent(8,0.04);
1529
5541b811 1530 return;
fa2be30c 1531
5541b811 1532}
6c388dee 1533
86164e35 1534//_________________________________________________________________________
1535void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
f21cf4a7 1536 // Dstar syst in pPb 2013 MB
1537 //
bd332272 1538 AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
86164e35 1539
86164e35 1540 //Normalization
1541 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1542 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.00); //
1543
f21cf4a7 1544 // Branching ratio
86164e35 1545 fBR = new TH1F("fBR","fBR",24,0,24);
1546 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1547
1548 // Tracking efficiency
1549 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1550 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1551
1552 // Raw yield extraction
1553 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
f21cf4a7 1554 fRawYield->SetBinContent(1,1.0);
30dc63e6 1555 fRawYield->SetBinContent(2,0.08);
1556 fRawYield->SetBinContent(3,0.02);
f21cf4a7 1557 fRawYield->SetBinContent(4,0.02);
30dc63e6 1558 fRawYield->SetBinContent(5,0.05);
1559 fRawYield->SetBinContent(6,0.02);
f21cf4a7 1560 fRawYield->SetBinContent(7,0.04);
1561 fRawYield->SetBinContent(8,0.04);
30dc63e6 1562 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
f21cf4a7 1563 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
1564 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
1565
86164e35 1566 // Cuts efficiency (from cuts variation)
1567 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
f21cf4a7 1568 fCutsEff->SetBinContent(1,0.);
1569 fCutsEff->SetBinContent(2,0.10);
1570 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
86164e35 1571
1572 // PID efficiency (from PID/noPID)
1573 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
f21cf4a7 1574 fPIDEff->SetBinContent(1,0.0);
86164e35 1575 fPIDEff->SetBinContent(2,0.05);
f21cf4a7 1576 for(Int_t i=3;i<=12;i++) fPIDEff->SetBinContent(i,0.00);
1577 for(Int_t i=12;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
86164e35 1578
1579 // MC dN/dpt
1580 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
c9569eb5 1581 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
1582 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
1583 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
86164e35 1584
1585 // particle-antiparticle
1586 /*
1587 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1588 fPartAntipart->SetBinContent(1,1);
1589 fPartAntipart->SetBinContent(2,1);
1590 fPartAntipart->SetBinContent(3,0.12);
1591 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1592 */
1593 return;
1594
1595
1596}
1597
414bc6a5 1598//--------------------------------------------------------------------------
1599void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
1600 // define errors for RAA vs. centrality
1601 InitD0toKpi2010PbPb020();
41cfcd39 1602 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 1603 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
41cfcd39 1604 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1605}
1606//--------------------------------------------------------------------------
1607void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
1608 // define errors for RAA vs. centrality
1609 InitD0toKpi2010PbPb020();
41cfcd39 1610 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 1611 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
41cfcd39 1612 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1613}
1614//--------------------------------------------------------------------------
1615void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
1616 // define errors for RAA vs. centrality
41cfcd39 1617 InitD0toKpi2010PbPb4080();
1618 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 1619 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.14);
41cfcd39 1620 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1621}
1622//--------------------------------------------------------------------------
1623void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
41cfcd39 1624 // define errors for RAA vs. centrality
414bc6a5 1625 InitD0toKpi2010PbPb4080();
41cfcd39 1626 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
27360a0e 1627 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.11);
41cfcd39 1628 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1629}
1630//--------------------------------------------------------------------------
1631void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
41cfcd39 1632 // define errors for RAA vs. centrality
414bc6a5 1633 InitD0toKpi2010PbPb4080();
41cfcd39 1634 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
27360a0e 1635 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
41cfcd39 1636 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1637}
1638//--------------------------------------------------------------------------
1639void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
1640 // define errors for RAA vs. centrality
1641 InitDplustoKpipi2010PbPb020();
7b03f1ec 1642 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
eefab142 1643 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
7b03f1ec 1644
414bc6a5 1645}
1646//--------------------------------------------------------------------------
1647void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
1648 // define errors for RAA vs. centrality
1649 InitDplustoKpipi2010PbPb020();
7b03f1ec 1650 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
eefab142 1651 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1652}
1653//--------------------------------------------------------------------------
1654void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
1655 // define errors for RAA vs. centrality
1656 InitDplustoKpipi2010PbPb020();
7b03f1ec 1657 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
eefab142 1658 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
414bc6a5 1659}
1660//--------------------------------------------------------------------------
1661void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
1662 // define errors for RAA vs. centrality
1663 InitDplustoKpipi2010PbPb4080();
7b03f1ec 1664 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
eefab142 1665 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1666}
1667//--------------------------------------------------------------------------
1668void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
1669 // define errors for RAA vs. centrality
1670 InitDplustoKpipi2010PbPb4080();
7b03f1ec 1671 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
eefab142 1672 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
414bc6a5 1673}
1674
1675//--------------------------------------------------------------------------
1676void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
1677 // define errors for RAA vs. centrality
1678 InitDstartoD0pi2010PbPb020();
1679 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
eefab142 1680 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
414bc6a5 1681}
1682//--------------------------------------------------------------------------
1683void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
1684 // define errors for RAA vs. centrality
1685 InitDstartoD0pi2010PbPb020();
1686 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
eefab142 1687 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
414bc6a5 1688}
1689//--------------------------------------------------------------------------
1690void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
1691 // define errors for RAA vs. centrality
1692 InitDstartoD0pi2010PbPb2040();
1693 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
eefab142 1694 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1695}
1696//--------------------------------------------------------------------------
1697void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
1698 // define errors for RAA vs. centrality
1699 InitDstartoD0pi2010PbPb4080();
1700 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
eefab142 1701 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
414bc6a5 1702}
1703//--------------------------------------------------------------------------
1704void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
1705 // define errors for RAA vs. centrality
1706 InitDstartoD0pi2010PbPb4080();
1707 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
eefab142 1708 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
414bc6a5 1709}
1710
3db225bf 1711//--------------------------------------------------------------------------
1712void AliHFSystErr::InitD0toKpi2011PbPb010CentScan(){
1713 // define errors for RAA vs. centrality 2011
1714 InitD0toKpi2011PbPb07half();
1715
1716}
1717//--------------------------------------------------------------------------
1718void AliHFSystErr::InitD0toKpi2011PbPb1020CentScan(){
1719 // define errors for RAA vs. centrality 2011
1720 InitD0toKpi2011PbPb07half();
1721
1722}
1723//--------------------------------------------------------------------------
1724void AliHFSystErr::InitD0toKpi2011PbPb2030CentScan(){
1725 // define errors for RAA vs. centrality 2011
1726 InitD0toKpi2011PbPb07half();
1727
1728}
1729//--------------------------------------------------------------------------
1730void AliHFSystErr::InitD0toKpi2011PbPb3040CentScan(){
1731 // define errors for RAA vs. centrality 2011
1732 InitD0toKpi2011PbPb07half();
1733
1734}
1735//--------------------------------------------------------------------------
1736void AliHFSystErr::InitD0toKpi2011PbPb4050CentScan(){
1737 // define errors for RAA vs. centrality 2011
1738 InitD0toKpi2011PbPb07half();
1739
1740}
1741//--------------------------------------------------------------------------
1742void AliHFSystErr::InitD0toKpi2010PbPb5080CentScan(){
1743 // define errors for RAA vs. centrality 2011
1744 InitD0toKpi2011PbPb07half();
1745
1746}
2eab2b1d 1747
3db225bf 1748//--------------------------------------------------------------------------
1749void AliHFSystErr::InitDplustoKpipi2011PbPb010CentScan(){
1750 // define errors for RAA vs. centrality 2011
1751 InitDplustoKpipi2011PbPb07half();
2eab2b1d 1752
1753 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.15);
1754 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
1755 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.15);
1756
1757 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1758 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1759 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
1760
1761 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
1762 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
1763 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.10);
3db225bf 1764}
2eab2b1d 1765
3db225bf 1766//--------------------------------------------------------------------------
1767void AliHFSystErr::InitDplustoKpipi2011PbPb1020CentScan(){
1768 // define errors for RAA vs. centrality 2011
1769 InitDplustoKpipi2011PbPb07half();
2eab2b1d 1770
1771 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
1772 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
1773 for(Int_t i=9;i<=17;i++)fRawYield->SetBinContent(i,0.15);
1774
1775 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1776 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1777 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
3db225bf 1778
2eab2b1d 1779 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
1780 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
1781 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.15);
3db225bf 1782}
2eab2b1d 1783
3db225bf 1784//--------------------------------------------------------------------------
1785void AliHFSystErr::InitDplustoKpipi2011PbPb2030CentScan(){
1786 // define errors for RAA vs. centrality 2011
1787 InitDplustoKpipi2011PbPb07half();
1788
2eab2b1d 1789 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
1790 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
1791 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.08);
1792
1793 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1794 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1795 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
1796
1797 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
1798 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
1799 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
3db225bf 1800}
2eab2b1d 1801
3db225bf 1802//--------------------------------------------------------------------------
1803void AliHFSystErr::InitDplustoKpipi2011PbPb3040CentScan(){
1804 // define errors for RAA vs. centrality 2011
1805 InitDplustoKpipi2011PbPb07half();
2eab2b1d 1806
1807 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.20);
1808 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
1809 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.20);
1810
1811 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1812 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1813 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
1814
1815 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
1816 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
1817 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
3db225bf 1818}
2eab2b1d 1819
3db225bf 1820//--------------------------------------------------------------------------
1821void AliHFSystErr::InitDplustoKpipi2011PbPb4050CentScan(){
1822 // define errors for RAA vs. centrality 2011
1823 InitDplustoKpipi2011PbPb07half();
2eab2b1d 1824
1825 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
1826 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
1827 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.05);
1828
1829 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1830 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1831 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
1832
1833 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
1834 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
1835 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
3db225bf 1836}
2eab2b1d 1837
3db225bf 1838//--------------------------------------------------------------------------
1839void AliHFSystErr::InitDplustoKpipi2010PbPb5080CentScan(){
1840 // define errors for RAA vs. centrality 2011
2eab2b1d 1841 InitDplustoKpipi2010PbPb4080();
1842
1843 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
1844 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
1845 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.10);
1846
1847 for(Int_t i=1;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
1848 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20); // 10%
1849 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1850 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20); // 10%
3db225bf 1851}
2eab2b1d 1852
3db225bf 1853//--------------------------------------------------------------------------
1854void AliHFSystErr::InitDstartoD0pi2011PbPb010CentScan(){
1855 // define errors for RAA vs. centrality 2011
1856 InitDstartoD0pi2011PbPb07half();
1857
1858}
1859//--------------------------------------------------------------------------
1860void AliHFSystErr::InitDstartoD0pi2011PbPb1020CentScan(){
1861 // define errors for RAA vs. centrality 2011
1862 InitDstartoD0pi2011PbPb07half();
1863
1864}
1865//--------------------------------------------------------------------------
1866void AliHFSystErr::InitDstartoD0pi2011PbPb2030CentScan(){
1867 // define errors for RAA vs. centrality 2011
1868 InitDstartoD0pi2011PbPb07half();
1869
1870}
1871//--------------------------------------------------------------------------
1872void AliHFSystErr::InitDstartoD0pi2011PbPb3040CentScan(){
1873 // define errors for RAA vs. centrality 2011
1874 InitDstartoD0pi2011PbPb07half();
1875
1876}
1877//--------------------------------------------------------------------------
1878void AliHFSystErr::InitDstartoD0pi2011PbPb4050CentScan(){
1879 // define errors for RAA vs. centrality 2011
1880 InitDstartoD0pi2011PbPb07half();
1881
1882}
1883//--------------------------------------------------------------------------
1884void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
1885 // define errors for RAA vs. centrality 2011
1886 InitDstartoD0pi2011PbPb07half();
1887
1888}
1889
44dd56b2 1890//_________________________________________________________________________
1891void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
1892 InitD0toKpi2013pPb0100();
1893
a1f0af09 1894 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
1895
1896 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
1897 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.06);
44dd56b2 1898 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
1899
44dd56b2 1900 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
1901
a1f0af09 1902 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
1903 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
1904 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1905
1906 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1907 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1908 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 1909
1910
1911};
1912//_________________________________________________________________________
1913void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
1914 InitD0toKpi2013pPb0100();
a1f0af09 1915
1916 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
44dd56b2 1917
a1f0af09 1918 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
1919 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.08);
44dd56b2 1920 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
1921
44dd56b2 1922 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
1923
a1f0af09 1924 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
1925 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
1926 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1927
1928 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1929 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1930 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 1931};
1932//_________________________________________________________________________
1933void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
a1f0af09 1934 InitD0toKpi2013pPb0100();
1935
1936 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
44dd56b2 1937
a1f0af09 1938 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
1939 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.12);
44dd56b2 1940 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
1941
44dd56b2 1942 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
1943
a1f0af09 1944 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
1945 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
1946 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1947
1948 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1949 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1950 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 1951
1952};
1953//_________________________________________________________________________
1954void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
1955 InitD0toKpi2013pPb0100();
1956
a1f0af09 1957 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
44dd56b2 1958
a1f0af09 1959 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
1960 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
1961 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
1962
1963 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
1964
1965 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
1966 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
1967 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1968
1969 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1970 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1971 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 1972
1973};
1974//_________________________________________________________________________
1975void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
1976 InitD0toKpi2013pPb0100();
1977
a1f0af09 1978 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
1979
1980 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
1981 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
1982 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
44dd56b2 1983
44dd56b2 1984 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
1985
a1f0af09 1986 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
1987 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
1988 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1989
1990 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1991 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1992 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 1993
1994};
1995
1996//_________________________________________________________________________
1997void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
1998 InitDplustoKpipi2013pPb0100();
1999
a1f0af09 2000 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
44dd56b2 2001
a1f0af09 2002 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2003 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2004 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2005
44dd56b2 2006 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2007
a1f0af09 2008 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2009 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2010 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2011
2012 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2013 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.10);
2014 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.10);
44dd56b2 2015
2016};
2017//_________________________________________________________________________
2018void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
2019 InitDplustoKpipi2013pPb0100();
a1f0af09 2020 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
44dd56b2 2021
a1f0af09 2022 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2023 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2024 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2025
44dd56b2 2026 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2027
a1f0af09 2028 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2029 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2030 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2031
2032 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2033 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2034 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2035
2036};
2037//_________________________________________________________________________
2038void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
2039 InitDplustoKpipi2013pPb0100();
a1f0af09 2040 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
44dd56b2 2041
a1f0af09 2042 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2043 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2044 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2045
44dd56b2 2046 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2047
a1f0af09 2048 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2049 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2050 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2051
2052 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2053 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2054 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2055
2056};
2057//_________________________________________________________________________
2058void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
2059 InitDplustoKpipi2013pPb0100();
a1f0af09 2060 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2061
2062 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2063 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2064 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
a1f0af09 2065
44dd56b2 2066 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2067
a1f0af09 2068 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2069 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2070 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
44dd56b2 2071
a1f0af09 2072 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2073 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2074 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2075
44dd56b2 2076};
2077//_________________________________________________________________________
2078void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
2079 InitDplustoKpipi2013pPb0100();
2080
a1f0af09 2081 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2082
2083 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2084 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2085 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
44dd56b2 2086
44dd56b2 2087 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2088
a1f0af09 2089 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2090 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2091 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
44dd56b2 2092
a1f0af09 2093 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2094 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2095 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2096
44dd56b2 2097};
2098
2099
2100//_________________________________________________________________________
2101void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
2102 InitDstartoD0pi2013pPb0100();
2103
a1f0af09 2104 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2105
2106 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2107 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
44dd56b2 2108 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2109
44dd56b2 2110 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2111
a1f0af09 2112 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2113 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2114 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2115
2116 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2117 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2118 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2119
2120};
2121//_________________________________________________________________________
2122void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
2123 InitDstartoD0pi2013pPb0100();
2124
a1f0af09 2125 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2126
2127 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2128 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2129 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2130
44dd56b2 2131 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2132
a1f0af09 2133 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2134 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2135 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2136
2137 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2138 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2139 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2140
2141};
2142//_________________________________________________________________________
2143void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
2144InitDstartoD0pi2013pPb0100();
2145
a1f0af09 2146 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2147
2148 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2149 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
44dd56b2 2150 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2151
44dd56b2 2152 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2153
a1f0af09 2154 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2155 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2156 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
44dd56b2 2157
a1f0af09 2158 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2159 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2160 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2161
44dd56b2 2162};
2163//_________________________________________________________________________
2164void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
2165 InitDstartoD0pi2013pPb0100();
a1f0af09 2166
2167 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
44dd56b2 2168
a1f0af09 2169 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2170 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
44dd56b2 2171 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2172
44dd56b2 2173 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2174
a1f0af09 2175 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2176 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2177 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2178
2179 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2180 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2181 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2182};
2183//_________________________________________________________________________
2184void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
2185 InitDstartoD0pi2013pPb0100();
a1f0af09 2186
2187 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
44dd56b2 2188
a1f0af09 2189 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2190 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.07);
44dd56b2 2191 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2192
44dd56b2 2193 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2194
a1f0af09 2195 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2196 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2197 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
44dd56b2 2198
a1f0af09 2199 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2200 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2201 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2202
44dd56b2 2203};
3db225bf 2204
414bc6a5 2205
9ca30701 2206//--------------------------------------------------------------------------
2207void AliHFSystErr::InitLctopKpi2010pp() {
2208 //
2209 // Lc->pKpi syst errors. Responsible: R. Romita
2210 // 2010 pp sample
2211 //
2212
2213 // Normalization
2214 fNorm = new TH1F("fNorm","fNorm",4,2,6);
2215 for(Int_t i=1;i<=4;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
2216
2217 // Tracking efficiency
2218 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",4,2,6);
2219 for(Int_t i=1;i<=4;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
2220
2221 // Raw yield extraction
2222 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2223 fRawYield->SetBinContent(1,0.20);
2224 fRawYield->SetBinContent(2,0.15);
2225 fRawYield->SetBinContent(3,0.10);
2226 fRawYield->SetBinContent(4,0.15);
2227
2228 fCutsEff = new TH1F("fCutsEff","fCutsEff",4,2,6);
2229 fCutsEff->SetBinContent(1,0.50);
2230 fCutsEff->SetBinContent(2,0.07);
2231 fCutsEff->SetBinContent(3,0.30);
2232 fCutsEff->SetBinContent(4,0.40);
2233
2234 // PID efficiency (from PID/noPID)
2235 fPIDEff = new TH1F("fPIDEff","fPIDEff",4,2,4);
2236 fPIDEff->SetBinContent(1,0.13); // 15%
2237 fPIDEff->SetBinContent(2,0.30); // 15%
2238 fPIDEff->SetBinContent(3,0.5); // 15%
2239 fPIDEff->SetBinContent(4,0.30); // 15%
2240
2241 // MC dN/dpt
2242 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
2243 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
2244 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
2245 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
2246 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
2247 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
2248}
2249
5541b811 2250//--------------------------------------------------------------------------
0233abe6 2251Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
2252 //
2253 // Get error
2254 //
2255
2256 Int_t bin=fCutsEff->FindBin(pt);
2257
2258 return fCutsEff->GetBinContent(bin);
2259}
2260//--------------------------------------------------------------------------
2261Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
2262 //
2263 // Get error
2264 //
2265
2266 Int_t bin=fMCPtShape->FindBin(pt);
2267
2268 return fMCPtShape->GetBinContent(bin);
2269}
2270//--------------------------------------------------------------------------
2271Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
2272 //
2273 // Get error
2274 //
2275
2276 Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
2277
2278 return TMath::Sqrt(err);
2279}
2280//--------------------------------------------------------------------------
2281Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
2282 //
2283 // Get error
2284 //
2285
2286 Int_t bin=fPIDEff->FindBin(pt);
2287
2288 return fPIDEff->GetBinContent(bin);
2289}
2290//--------------------------------------------------------------------------
2291Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
2292 //
2293 // Get error
2294 //
2295
2296 Int_t bin=fTrackingEff->FindBin(pt);
2297
2298 return fTrackingEff->GetBinContent(bin);
2299}
2300//--------------------------------------------------------------------------
2301Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
2302 //
2303 // Get error
2304 //
2305
2306 Int_t bin=fRawYield->FindBin(pt);
2307
2308 return fRawYield->GetBinContent(bin);
2309}
2310//--------------------------------------------------------------------------
2311Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
2312 //
2313 // Get error
2314 //
2315
2316 Int_t bin=fPartAntipart->FindBin(pt);
2317
2318 return fPartAntipart->GetBinContent(bin);
2319}
2320//--------------------------------------------------------------------------
2321Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
2322 //
2323 // Get total syst error (except norm. error)
2324 //
2325
2326 Double_t err=0.;
2327
2328 if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
2329 if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
f9a85c82 2330 // if(fBR) err += GetBRErr()*GetBRErr();
0233abe6 2331 if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
2332 if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
2333 if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
2334 if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
2335
2336 err += feeddownErr*feeddownErr;
2337
2338 return TMath::Sqrt(err);
2339}
2340//---------------------------------------------------------------------------
5ad8ad40 2341void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
0233abe6 2342 //
2343 // Draw errors
2344 //
2345 gStyle->SetOptStat(0);
2346
a07118d0 2347 TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
2348 cSystErr->Range(0.20,-0.5,18.4,0.34);
2349 cSystErr->SetRightMargin(0.318);
0233abe6 2350 cSystErr->SetFillColor(0);
2351
436b185f 2352 TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
2353 hFrame->SetAxisRange(1.,35.9,"X");
a07118d0 2354 hFrame->SetAxisRange(-0.5,0.5,"Y");
0233abe6 2355 hFrame->Draw();
2356
a07118d0 2357 TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
2358 leg->SetTextSize(0.03601695);
0233abe6 2359 leg->SetFillStyle(0);
2360 leg->SetBorderSize(0);
2ee3afe2 2361
436b185f 2362 TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
2ee3afe2 2363 Int_t nbins = fNorm->GetNbinsX();
2364 TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
436b185f 2365 for(Int_t i=1;i<=36;i++) {
2ee3afe2 2366 Double_t pt = hTotErr->GetBinCenter(i);
2367 Double_t ptwidth = hTotErr->GetBinWidth(i);
2368
2369 if(grErrFeeddown) {
2370 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
2371 Double_t toterryl=0., toterryh=0.;
2372 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
2373 grErrFeeddown->GetPoint(j,x,y);
2374 errxh = grErrFeeddown->GetErrorXhigh(j);
2375 errxl = grErrFeeddown->GetErrorXlow(j);
a07118d0 2376 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
2ee3afe2 2377 erryh = grErrFeeddown->GetErrorYhigh(j);
2378 erryl = grErrFeeddown->GetErrorYlow(j);
2379 }
2380 }
2381 if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
2382 else toterryl = GetTotalSystErr(pt);
2383 if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
2384 else toterryh = GetTotalSystErr(pt);
2385
2386 hTotErr->SetBinContent(i,toterryh);
2387 gTotErr->SetPoint(i,pt,0.);
2388 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
2389 }
2390 else {
2391 hTotErr->SetBinContent(i,GetTotalSystErr(pt));
2392 gTotErr->SetPoint(i,pt,0.);
2393 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
2394 }
2395
2396 }
f968ef30 2397 gTotErr->SetLineColor(kBlack);
2398 gTotErr->SetFillColor(kRed);
2399 gTotErr->SetFillStyle(3002);
2ee3afe2 2400 gTotErr->Draw("2");
2401 leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
2402// hTotErr->SetLineColor(1);
2403// hTotErr->SetLineWidth(3);
2404// hTotErr->Draw("same");
2405// leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
2406
0233abe6 2407
e373b9f0 2408 fNorm->SetFillColor(1);
2409 fNorm->SetFillStyle(3002);
2410 //fNorm->Draw("same");
2411 //TH1F *hNormRefl = ReflectHisto(fNorm);
2412 //hNormRefl->Draw("same");
27360a0e 2413 Double_t norm = fNorm->GetBinContent(1)*100;
2414 leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
e373b9f0 2415
5ad8ad40 2416 if(grErrFeeddown) {
2417 grErrFeeddown->SetFillColor(kTeal-8);
2418 grErrFeeddown->SetFillStyle(3001);
2419 grErrFeeddown->Draw("2");
2420 leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
2421 }
0233abe6 2422 if(fTrackingEff) {
5ad8ad40 2423 fTrackingEff->SetFillColor(4);
2424 fTrackingEff->SetFillStyle(3006);
0233abe6 2425 fTrackingEff->Draw("same");
f968ef30 2426 TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
2427 hTrackingEffRefl->Draw("same");
0233abe6 2428 leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
2429 }
2430 if(fBR) {
2431 fBR->SetFillColor(6);
2432 fBR->SetFillStyle(3005);
5ad8ad40 2433 //fBR->SetFillStyle(3020);
0233abe6 2434 fBR->Draw("same");
f968ef30 2435 TH1F *hBRRefl = ReflectHisto(fBR);
2436 hBRRefl->Draw("same");
0233abe6 2437 leg->AddEntry(fBR,"Branching ratio","f");
2438 }
2439 if(fRawYield) {
5ad8ad40 2440 Int_t ci; // for color index setting
2441 ci = TColor::GetColor("#00cc00");
2442 fRawYield->SetLineColor(ci);
2443 // fRawYield->SetLineColor(3);
0233abe6 2444 fRawYield->SetLineWidth(3);
2445 fRawYield->Draw("same");
f968ef30 2446 TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
2447 hRawYieldRefl->Draw("same");
844405bc 2448 leg->AddEntry(fRawYield,"Yield extraction","l");
0233abe6 2449 }
2450 if(fCutsEff) {
2451 fCutsEff->SetLineColor(4);
2452 fCutsEff->SetLineWidth(3);
2453 fCutsEff->Draw("same");
f968ef30 2454 TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
2455 hCutsEffRefl->Draw("same");
0233abe6 2456 leg->AddEntry(fCutsEff,"Cuts efficiency","l");
2457 }
2458 if(fPIDEff) {
2459 fPIDEff->SetLineColor(7);
2460 fPIDEff->SetLineWidth(3);
2461 fPIDEff->Draw("same");
f968ef30 2462 TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
2463 hPIDEffRefl->Draw("same");
0233abe6 2464 leg->AddEntry(fPIDEff,"PID efficiency","l");
2465 }
2466 if(fMCPtShape) {
5ad8ad40 2467 Int_t ci = TColor::GetColor("#9933ff");
2468 fMCPtShape->SetLineColor(ci);
2469 // fMCPtShape->SetLineColor(8);
0233abe6 2470 fMCPtShape->SetLineWidth(3);
2471 fMCPtShape->Draw("same");
f968ef30 2472 TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
2473 hMCPtShapeRefl->Draw("same");
0233abe6 2474 leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
2475 }
2476 if(fPartAntipart) {
5ad8ad40 2477 Int_t ci = TColor::GetColor("#ff6600");
2478 fPartAntipart->SetLineColor(ci);
2479 // fPartAntipart->SetLineColor(9);
0233abe6 2480 fPartAntipart->SetLineWidth(3);
2481 fPartAntipart->Draw("same");
f968ef30 2482 TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
2483 hPartAntipartRefl->Draw("same");
0233abe6 2484 leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
2485 }
0233abe6 2486
0233abe6 2487
2488 leg->Draw();
2489
5ad8ad40 2490 cSystErr->SaveAs("RelativeSystematics.eps");
2491
0233abe6 2492 return;
2493}
f968ef30 2494//-------------------------------------------------------------------------
2495TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
2496 //
2497 // Clones and reflects histogram
2498 //
2499 TH1F *hout=(TH1F*)hin->Clone("hout");
2500 hout->Scale(-1.);
2501
2502 return hout;
2503}
9a8b5080 2504
ada850b0 2505
2506