]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG3/vertexingHF/AliHFSystErr.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
memory usage optimised
[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliHFSystErr.cxx
CommitLineData
0233abe6 1/**************************************************************************
2 * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
27de2dfb 16/* $Id$ */
17
0233abe6 18/////////////////////////////////////////////////////////////
19//
20// Class to handle systematic errors for charm hadrons
21//
22// Usage:
5541b811 23// AliHFSystEff syst; // DECAY = 1 for D0, 2, for D+, 3 for D*
24// syst.SetRunNumber(YEAR); // YEAR = two last numbers of the year (is 10 for 2010)
25// syst.SetCollisionType(TYPE); // TYPE = 0 is pp, 1 is PbPb
26// syst.SetCentrality(CENT); // CENT is centrality, 0100 for MB, 020 (4080) for 0-20 (40-80) CC...
27// syst.Init(DECAY); // DECAY = 1 for D0, 2, for D+, 3 for D*
0233abe6 28// syst.DrawErrors(); // to see a plot of the error contributions
29// syst.GetTotalSystErr(pt); // to get the total err at pt
30//
31// Author: A.Dainese, andrea.dainese@pd.infn.it
32/////////////////////////////////////////////////////////////
33
34#include <TStyle.h>
35#include <TGraphAsymmErrors.h>
36#include <TMath.h>
37#include <TCanvas.h>
38#include <TH2F.h>
39#include <TLegend.h>
5ad8ad40 40#include <TColor.h>
0233abe6 41
5541b811 42#include "AliLog.h"
0233abe6 43#include "AliHFSystErr.h"
44
45
46ClassImp(AliHFSystErr)
47
48//--------------------------------------------------------------------------
49AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
50TNamed(name,title),
51fNorm(0),
52fRawYield(0),
53fTrackingEff(0),
54fBR(0),
55fCutsEff(0),
56fPIDEff(0),
57fMCPtShape(0),
5541b811 58fPartAntipart(0),
59fRunNumber(10),
60fCollisionType(0),
a8c3295d 61fCentralityClass("0100"),
62fIsLowEnergy(false)
0233abe6 63{
64 //
65 // Default Constructor
66 //
67}
5541b811 68
0233abe6 69//--------------------------------------------------------------------------
5541b811 70AliHFSystErr::~AliHFSystErr() {
71 //
72 // Default Destructor
0233abe6 73 //
a8c3295d 74 /*
5541b811 75
76 if(fNorm) { delete fNorm; fNorm=0; }
77 if(fRawYield) { delete fRawYield; fRawYield=0; }
78 if(fTrackingEff) { delete fTrackingEff; fTrackingEff=0; }
79 if(fBR) { delete fBR; fBR=0; }
80 if(fCutsEff) { delete fCutsEff; fCutsEff=0; }
81 if(fPIDEff) { delete fPIDEff; fPIDEff=0; }
82 if(fMCPtShape) { delete fMCPtShape; fMCPtShape=0; }
83 if(fPartAntipart) { delete fPartAntipart; fPartAntipart=0; }
a8c3295d 84 */
5541b811 85}
86
87//--------------------------------------------------------------------------
88void AliHFSystErr::Init(Int_t decay){
89 //
90 // Variables/histos initialization
0233abe6 91 //
92
a8c3295d 93 if (fRunNumber!=10 && fIsLowEnergy==false) {
94 AliError("Only settings for 2010 and the low energy runs are implemented so far");
5541b811 95 }
aebba8df 96 if (fCentralityClass!="020" && fCentralityClass!="4080" && fCentralityClass!="0100"){
5541b811 97 AliError("Only settings for MB2010 are implemented so far");
98 }
99
0233abe6 100 switch(decay) {
101 case 1: // D0->Kpi
a8c3295d 102 if (fCollisionType==0) {
103 if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
104 else InitD0toKpi2010pp();
105 } else if (fCollisionType==1) {
aebba8df 106 if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
107 else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
5541b811 108 else AliError("Not yet implemented");
109 }
a8c3295d 110 // else if (fCollisionType==2) InitD0toKpi2010ppLowEn();
0233abe6 111 break;
112 case 2: // D+->Kpipi
a8c3295d 113 if (fCollisionType==0) {
114 if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
115 else InitDplustoKpipi2010pp();
116 } else if (fCollisionType==1) {
aebba8df 117 if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
118 else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
5541b811 119 else AliError("Not yet implemented");
120 }
0233abe6 121 break;
122 case 3: // D*->D0pi
fa2be30c 123 if (fCollisionType==0) {
124 if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
125 else InitDstartoD0pi2010pp();
126 }else if (fCollisionType==1) {
127 if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
128 if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
129 if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
130 if (fCentralityClass!="4080" && fCentralityClass!="2040" && fCentralityClass!="020") AliError("Not yet implemented");
5541b811 131 }
0233abe6 132 break;
133 default:
134 printf("Invalid decay type: %d\n",decay);
135 break;
136 }
137
0233abe6 138}
5541b811 139
0233abe6 140//--------------------------------------------------------------------------
5541b811 141void AliHFSystErr::InitD0toKpi2010pp() {
142 //
143 // D0->Kpi syst errors. Responsible: A. Rossi
144 // 2010 pp sample
0233abe6 145 //
146
5541b811 147 // Normalization
f1d72175 148 fNorm = new TH1F("fNorm","fNorm",24,0,24);
149 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
5541b811 150
151 // Branching ratio
f1d72175 152 fBR = new TH1F("fBR","fBR",24,0,24);
153 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
5541b811 154
155 // Tracking efficiency
f1d72175 156 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
157 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.02); // 2% (1% per track)
5541b811 158
159 // Raw yield extraction
f1d72175 160 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
5541b811 161 fRawYield->SetBinContent(1,1);
f1d72175 162 fRawYield->SetBinContent(2,0.2);
163 for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
164 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
165 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
0233abe6 166
5541b811 167 // Cuts efficiency (from cuts variation)
f1d72175 168 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
169 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
5541b811 170
171 // PID efficiency (from PID/noPID)
f1d72175 172 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
173 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
174 fPIDEff->SetBinContent(2,0.10); // 10%
5541b811 175
176 // MC dN/dpt
f1d72175 177 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 178 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
179 fMCPtShape->SetBinContent(1,0.03);
180 fMCPtShape->SetBinContent(2,0.03);
5541b811 181
182 // particle-antiparticle
f1d72175 183 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
184 // fPartAntipart->SetBinContent(1,1);
185 // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
5541b811 186
187 return;
0233abe6 188}
189//--------------------------------------------------------------------------
5541b811 190void AliHFSystErr::InitD0toKpi2010PbPb020() {
0233abe6 191 //
aebba8df 192 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 193 // 2010 PbPb sample, 0-20 CC
0233abe6 194 //
5541b811 195
196 // Normalization
197 fNorm = new TH1F("fNorm","fNorm",20,0,20);
198 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
199
200 // Branching ratio
201 fBR = new TH1F("fBR","fBR",20,0,20);
202 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
203
204 // Tracking efficiency
205 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
aebba8df 206 for(Int_t i=2;i<=12;i++) fTrackingEff->SetBinContent(i,0.05+0.005*(Float_t)i);
666cf569 207
5541b811 208
209 // Raw yield extraction
210 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
211 fRawYield->SetBinContent(1,0);
212 fRawYield->SetBinContent(2,0);
aebba8df 213 fRawYield->SetBinContent(3,0.05);
214 fRawYield->SetBinContent(4,0.05);
215 fRawYield->SetBinContent(5,0.10);
216 fRawYield->SetBinContent(6,0.10);
217 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
666cf569 218 // for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
5541b811 219
220 // Cuts efficiency (from cuts variation)
221 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 222 fCutsEff->SetBinContent(1,0.);
223 fCutsEff->SetBinContent(2,0.);
666cf569 224 fCutsEff->SetBinContent(3,0.18);
aebba8df 225 fCutsEff->SetBinContent(4,0.18);
226 fCutsEff->SetBinContent(5,0.18);
227 fCutsEff->SetBinContent(6,0.15);
228 fCutsEff->SetBinContent(7,0.15);
229 fCutsEff->SetBinContent(8,0.15);
666cf569 230fCutsEff->SetBinContent(9,0.15);
231 fCutsEff->SetBinContent(10,0.15);
aebba8df 232 fCutsEff->SetBinContent(11,0.15);
233 fCutsEff->SetBinContent(12,0.15);
666cf569 234 // for(Int_t i=13;i<=20;i++) fCutsEff->SetBinContent(i,0.);
5541b811 235
236 // PID efficiency (from PID/noPID)
237 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
aebba8df 238 for(Int_t i=3;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
239 fPIDEff->SetBinContent(5,0.10);
240 fPIDEff->SetBinContent(6,0.10);
5541b811 241
242 // MC dN/dpt
243 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 244 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
245 fMCPtShape->SetBinContent(1,0.03);
246 fMCPtShape->SetBinContent(2,0.03);
247
5541b811 248
249 // particle-antiparticle
250 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
aebba8df 251 for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
666cf569 252 fPartAntipart->SetBinContent(3,0.10);
aebba8df 253 fPartAntipart->SetBinContent(4,0.10);
254 fPartAntipart->SetBinContent(7,0.10);
255 fPartAntipart->SetBinContent(8,0.10);
0233abe6 256
5541b811 257 return;
258}
259//--------------------------------------------------------------------------
260void AliHFSystErr::InitD0toKpi2010PbPb4080() {
261 //
aebba8df 262 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 263 // 2010 PbPb sample, 40-80 CC
264 //
265
0233abe6 266 // Normalization
267 fNorm = new TH1F("fNorm","fNorm",20,0,20);
844405bc 268 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
0233abe6 269
270 // Branching ratio
271 fBR = new TH1F("fBR","fBR",20,0,20);
844405bc 272 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
0233abe6 273
aebba8df 274 // Tracking efficiency
275 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
276 for(Int_t i=2;i<=12;i++) fTrackingEff->SetBinContent(i,0.5*(0.05+0.005*(Float_t)i));
277
278
0233abe6 279 // Raw yield extraction
280 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
5541b811 281 fRawYield->SetBinContent(1,0);
282 fRawYield->SetBinContent(2,0);
46cfcb37 283 for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.05);
666cf569 284 //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
0233abe6 285
286 // Cuts efficiency (from cuts variation)
287 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 288 fCutsEff->SetBinContent(1,0.);
289 fCutsEff->SetBinContent(2,0.);
666cf569 290 fCutsEff->SetBinContent(3,0.15);
46cfcb37 291 fCutsEff->SetBinContent(4,0.15);
292 fCutsEff->SetBinContent(5,0.15);
aebba8df 293 fCutsEff->SetBinContent(6,0.15);
294 fCutsEff->SetBinContent(7,0.15);
295 fCutsEff->SetBinContent(8,0.15);
296 fCutsEff->SetBinContent(9,0.15);
297 fCutsEff->SetBinContent(10,0.15);
298 fCutsEff->SetBinContent(11,0.15);
299 fCutsEff->SetBinContent(12,0.15);
666cf569 300 // for(Int_t i=13;i<=20;i++) fCutsEff->SetBinContent(i,0.);
0233abe6 301
302 // PID efficiency (from PID/noPID)
303 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
aebba8df 304 for(Int_t i=3;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
0233abe6 305
306 // MC dN/dpt
307 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 308 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
309 fMCPtShape->SetBinContent(1,0.03);
310 fMCPtShape->SetBinContent(2,0.03);
0233abe6 311
312 // particle-antiparticle
313 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
aebba8df 314 for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
0233abe6 315
316 return;
317}
a8c3295d 318
319//--------------------------------------------------------------------------
320void AliHFSystErr::InitD0toKpi2010ppLowEn() {
321 //
322 // D0->Kpi syst errors. Low energy run
323 // 2011 2.76 TeV pp sample
324 //
325 AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
326
327 // Normalization
328 fNorm = new TH1F("fNorm","fNorm",20,0,20);
329 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
330
331 // Branching ratio
332 fBR = new TH1F("fBR","fBR",20,0,20);
333 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
334
335 // Tracking efficiency
336 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
337 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
338
339 // Raw yield extraction
340 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
341 fRawYield->SetBinContent(1,1);
342 fRawYield->SetBinContent(2,1);
343 fRawYield->SetBinContent(3,0.2);
344 fRawYield->SetBinContent(4,0.2);
345 fRawYield->SetBinContent(5,0.1);
346 fRawYield->SetBinContent(6,0.1);
347 fRawYield->SetBinContent(7,0.2);
348 fRawYield->SetBinContent(8,0.2);
349 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
350
351 // Cuts efficiency (from cuts variation)
352 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
353 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.18); // 10%
354
355 // PID efficiency (from PID/noPID)
356 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
357 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.20); // 10%
358 fPIDEff->SetBinContent(3,0.25); // 10%
359 fPIDEff->SetBinContent(4,0.25); // 10%
360 fPIDEff->SetBinContent(7,0.25); // 10%
361 fPIDEff->SetBinContent(8,0.25); // 10%
362
363 // MC dN/dpt
364 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 365 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
366 fMCPtShape->SetBinContent(1,0.03);
367 fMCPtShape->SetBinContent(2,0.03);
a8c3295d 368
369 // particle-antiparticle
370 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
371 fPartAntipart->SetBinContent(1,1);
372 fPartAntipart->SetBinContent(2,1);
373 for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
374
375 return;
376}
377
0233abe6 378//--------------------------------------------------------------------------
5541b811 379void AliHFSystErr::InitDplustoKpipi2010pp() {
0233abe6 380 //
381 // D+->Kpipi syst errors. Responsible: R. Bala
5541b811 382 // 2010 pp sample
0233abe6 383 //
384
626bcf57 385 // Normalization
386 fNorm = new TH1F("fNorm","fNorm",20,0,20);
844405bc 387 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
626bcf57 388
389 // Branching ratio
390 fBR = new TH1F("fBR","fBR",20,0,20);
844405bc 391 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.04); // 4% PDG2010
626bcf57 392
393 // Tracking efficiency
394 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
844405bc 395 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.03); // 3% (1% per track)
626bcf57 396
397
398 // Raw yield extraction
399 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
400 fRawYield->SetBinContent(1,1);
401 fRawYield->SetBinContent(2,1);
5c73c2f5 402 fRawYield->SetBinContent(3,0.20);
5ad8ad40 403 for(Int_t i=4;i<=20;i++) fRawYield->SetBinContent(i,0.055); //5 to 10%
626bcf57 404
405 // Cuts efficiency (from cuts variation)
406 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
407 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
408
409 // PID efficiency (from PID/noPID)
410 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
411 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
412 fPIDEff->SetBinContent(3,0.13); // 13%
413
414
415 // MC dN/dpt (copied from D0 : will update later)
416 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 417 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
418 fMCPtShape->SetBinContent(1,0.03);
419 fMCPtShape->SetBinContent(2,0.03);
626bcf57 420
421
422 // particle-antiparticle
423 /*
424 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
425 fPartAntipart->SetBinContent(1,1);
426 fPartAntipart->SetBinContent(2,1);
427 fPartAntipart->SetBinContent(3,0.12);
428 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
429 */
0233abe6 430 return;
431}
432//--------------------------------------------------------------------------
5541b811 433void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
434 //
435 // D+->Kpipi syst errors. Responsible: ??
436 // 2010 PbPb sample, 0-20 CC
437 //
aebba8df 438
aebba8df 439 // Normalization
440 fNorm = new TH1F("fNorm","fNorm",20,0,20);
441 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
442
443 // Branching ratio
444 fBR = new TH1F("fBR","fBR",20,0,20);
445 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.04); // 4% PDG2010
446
447 // Tracking efficiency
448 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
666cf569 449 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,i,0.075+0.005*(Float_t)i);
aebba8df 450
451
452 // Raw yield extraction
453 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
454 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
666cf569 455// for(Int_t i=6;i<=7;i++) fRawYield->SetBinContent(i,0.32);
456// for(Int_t i=7;i<=9;i++) fRawYield->SetBinContent(i,0.32);
457// for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);//31./156.);
a8c3295d 458 fRawYield->SetBinContent(6,0.12);
459 fRawYield->SetBinContent(7,0.17);
460 fRawYield->SetBinContent(8,0.17);
461 fRawYield->SetBinContent(9,0.08);
462 fRawYield->SetBinContent(10,0.08);
463 fRawYield->SetBinContent(11,0.08);
464 fRawYield->SetBinContent(12,0.08);
aebba8df 465
466 // Cuts efficiency (from cuts variation)
467 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
666cf569 468 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
aebba8df 469
470 // PID efficiency (from PID/noPID)
471 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
472 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
473 fPIDEff->SetBinContent(3,0.13); // 13%
474
475
476 // MC dN/dpt (copied from D0 : will update later)
477 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 478 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
479 fMCPtShape->SetBinContent(1,0.03);
480 fMCPtShape->SetBinContent(2,0.03);
aebba8df 481
482
483 // particle-antiparticle
484 /*
485 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
486 fPartAntipart->SetBinContent(1,1);
487 fPartAntipart->SetBinContent(2,1);
488 fPartAntipart->SetBinContent(3,0.12);
489 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
490 */
5541b811 491 return;
492}
493
494//--------------------------------------------------------------------------
495void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
496 //
497 // D+->Kpipi syst errors. Responsible: ??
498 // 2010 PbPb sample, 40-80 CC
499 //
500
aebba8df 501
502 // Normalization
503 fNorm = new TH1F("fNorm","fNorm",20,0,20);
504 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
505
506 // Branching ratio
507 fBR = new TH1F("fBR","fBR",20,0,20);
508 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.04); // 4% PDG2010
509
510 // Tracking efficiency
511 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
512 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.03); // 3% (1% per track)
513
514
515 // Raw yield extraction
516 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
517 fRawYield->SetBinContent(1,1);
518 fRawYield->SetBinContent(2,1);
519 fRawYield->SetBinContent(3,1);
520 fRawYield->SetBinContent(4,1.);
521 fRawYield->SetBinContent(5,4./23.);
522 fRawYield->SetBinContent(6,4./23.);
523 fRawYield->SetBinContent(7,5./20.);
524 fRawYield->SetBinContent(8,5./20.);
525 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,1./16.);
526 for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
527
528 // Cuts efficiency (from cuts variation)
529 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
530 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
531
532 // PID efficiency (from PID/noPID)
533 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
534 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
535 fPIDEff->SetBinContent(3,0.13); // 13%
536
537
538 // MC dN/dpt (copied from D0 : will update later)
539 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 540 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
541 fMCPtShape->SetBinContent(1,0.03);
542 fMCPtShape->SetBinContent(2,0.03);
aebba8df 543
544
545 // particle-antiparticle
546 /*
547 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
548 fPartAntipart->SetBinContent(1,1);
549 fPartAntipart->SetBinContent(2,1);
550 fPartAntipart->SetBinContent(3,0.12);
551 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
552 */
5541b811 553 return;
554}
555
a8c3295d 556//--------------------------------------------------------------------------
557void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
558
559 //
560 // D+->Kpipi syst errors. Responsible: R. Bala
561 // 2011 2.76 TeV pp sample
562 //
563 AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
564
565 // Normalization
566 fNorm = new TH1F("fNorm","fNorm",20,0,20);
567 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
568
569 // Branching ratio
570 fBR = new TH1F("fBR","fBR",20,0,20);
571 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.04); // 4% PDG2010
572
573 // Tracking efficiency
574 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
575 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
576
577 // Raw yield extraction
578 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
579 fRawYield->SetBinContent(1,1);
580 fRawYield->SetBinContent(2,1);
581 fRawYield->SetBinContent(3,0.20);
582 fRawYield->SetBinContent(4,0.20);
583 for(Int_t i=5;i<=20;i++) fRawYield->SetBinContent(i,0.055); //5 to 10%
584
585 // Cuts efficiency (from cuts variation)
586 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
587 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
588
589 // PID efficiency (from PID/noPID)
590 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
591 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
592 fPIDEff->SetBinContent(3,0.15); // 13%
593 fPIDEff->SetBinContent(4,0.15); // 13%
594
595 // MC dN/dpt (copied from D0 : will update later)
596 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 597 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
598 fMCPtShape->SetBinContent(1,0.03);
599 fMCPtShape->SetBinContent(2,0.03);
a8c3295d 600
601 return;
602}
603
5541b811 604//--------------------------------------------------------------------------
605void AliHFSystErr::InitDstartoD0pi2010pp() {
0233abe6 606 //
2cf98d0d 607 // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
5541b811 608 // 2010 pp sample
0233abe6 609 //
610
2cf98d0d 611 // Normalization
fa2be30c 612 fNorm = new TH1F("fNorm","fNorm",24,0,24);
613 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
614
615 // Branching ratio
616 fBR = new TH1F("fBR","fBR",24,0,24);
617 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
618
619 // Tracking efficiency
620 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
621 fTrackingEff->SetBinContent(1,1.0);
622 fTrackingEff->SetBinContent(2,0.12);
623 fTrackingEff->SetBinContent(3,0.08);
624 fTrackingEff->SetBinContent(3,0.05);
625 for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.03); // 3% (1% per track)
626
627
628 // Raw yield extraction
629 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
630 fRawYield->SetBinContent(1,1.0);
631 fRawYield->SetBinContent(2,0.09);
632 fRawYield->SetBinContent(3,0.04);
633 fRawYield->SetBinContent(4,0.02);
634 fRawYield->SetBinContent(5,0.03);
08ae44d3 635 fRawYield->SetBinContent(6,0.06);
636 fRawYield->SetBinContent(7,0.04);
637 fRawYield->SetBinContent(8,0.11);
fa2be30c 638 for(Int_t i=5;i<=24;i++) fRawYield->SetBinContent(i,0.03); //4%
639 fRawYield->SetBinContent(13,0.09);
640 fRawYield->SetBinContent(14,0.09);
641 fRawYield->SetBinContent(15,0.09);
642 fRawYield->SetBinContent(16,0.09);
643 fRawYield->SetBinContent(17,0.24);
644 fRawYield->SetBinContent(18,0.24);
645 fRawYield->SetBinContent(19,0.24);
646 fRawYield->SetBinContent(20,0.24);
647 fRawYield->SetBinContent(21,0.24);
648 fRawYield->SetBinContent(22,0.24);
649 fRawYield->SetBinContent(23,0.24);
650 fRawYield->SetBinContent(24,0.24);
651
652 // Cuts efficiency (from cuts variation)
653 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
654 fCutsEff->SetBinContent(2,0.22);
655 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
656
657 // PID efficiency (from PID/noPID)
658 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
659 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
660
661
662 // MC dN/dpt (copied from D0 : will update later)
663 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 664 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
665 fMCPtShape->SetBinContent(1,0.03);
666 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 667
668 return;
669
670
671}
672//--------------------------------------------------------------------------
673void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
674
675 //
676 // D+->Kpipi syst errors. Responsible: A. Grelli
677 // 2011 2.76 TeV pp sample
678 //
679 AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
680
681// Normalization
2cf98d0d 682 fNorm = new TH1F("fNorm","fNorm",20,0,20);
683 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
684
685 // Branching ratio
686 fBR = new TH1F("fBR","fBR",20,0,20);
8accac87 687 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
2cf98d0d 688
689 // Tracking efficiency
690 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
08ae44d3 691 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
2cf98d0d 692
693 // Raw yield extraction
694 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
8accac87 695 fRawYield->SetBinContent(1,1);
fa2be30c 696 fRawYield->SetBinContent(2,1);
697 fRawYield->SetBinContent(3,0.2);
698 fRawYield->SetBinContent(4,0.2);
08ae44d3 699 fRawYield->SetBinContent(5,0.08);
700 fRawYield->SetBinContent(6,0.08);
701 fRawYield->SetBinContent(7,0.1);
702 fRawYield->SetBinContent(8,0.1);
fa2be30c 703 fRawYield->SetBinContent(9,0.2);
704 fRawYield->SetBinContent(10,0.2);
705 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
2cf98d0d 706
707 // Cuts efficiency (from cuts variation)
708 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
fa2be30c 709 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.18);
710 fCutsEff->SetBinContent(3,0.2);
711 fCutsEff->SetBinContent(4,0.2);
712 fCutsEff->SetBinContent(5,0.1);
713 fCutsEff->SetBinContent(6,0.1);
714 fCutsEff->SetBinContent(7,0.12);
715 fCutsEff->SetBinContent(8,0.12);
716 fCutsEff->SetBinContent(9,0.2);
717 fCutsEff->SetBinContent(10,0.2);
066c6b9b 718 fCutsEff->SetBinContent(11,0.2);
719 fCutsEff->SetBinContent(12,0.2);
2cf98d0d 720
721 // PID efficiency (from PID/noPID)
722 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
fa2be30c 723 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
2cf98d0d 724
fa2be30c 725 // MC dN/dpt
2cf98d0d 726 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 727 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
728 fMCPtShape->SetBinContent(1,0.03);
729 fMCPtShape->SetBinContent(2,0.03);
2cf98d0d 730
fa2be30c 731
0233abe6 732 return;
733}
fa2be30c 734//------------------------------------------------------------------------
5541b811 735void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
736 //
fa2be30c 737 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 738 // 2010 PbPb sample, 0-20 CC
739 //
740
fa2be30c 741 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
742
743 // Normalization
744 fNorm = new TH1F("fNorm","fNorm",24,0,24);
745 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
746
747 // Branching ratio
748 fBR = new TH1F("fBR","fBR",24,0,24);
749 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
750
751 // Tracking efficiency
752 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
753 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 3% (1% per track)
754
755
756 // Raw yield extraction
757 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
758 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
759
760 // Cuts efficiency (from cuts variation)
761 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
762 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
763
764 // PID efficiency (from PID/noPID)
765 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
766 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
767
768
769 // MC dN/dpt (copied from D0 : will update later)
770 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 771 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
772 fMCPtShape->SetBinContent(1,0.03);
773 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 774
775 return;
776
777}
778//-------------------------------------------------------------------------
779void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
780 //
781 // D*+->D0pi syst errors. Responsible: A. Grelli
782 // 2010 PbPb sample, 20-40 CC
783 //
784
785 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
786
787 // Normalization
788 fNorm = new TH1F("fNorm","fNorm",24,0,24);
789 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
790
791 // Branching ratio
792 fBR = new TH1F("fBR","fBR",24,0,24);
793 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
794
795 // Tracking efficiency
796 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
797 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 3% (1% per track)
798
799
800 // Raw yield extraction
801 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
802 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
803
804 // Cuts efficiency (from cuts variation)
805 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
806 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
807
808 // PID efficiency (from PID/noPID)
809 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
810 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
811
812
813 // MC dN/dpt (copied from D0 : will update later)
814 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 815 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
816 fMCPtShape->SetBinContent(1,0.03);
817 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 818
5541b811 819 return;
fa2be30c 820
5541b811 821}
822//--------------------------------------------------------------------------
823void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
824 //
fa2be30c 825 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 826 // 2010 PbPb sample, 40-80 CC
827 //
828
fa2be30c 829 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
830
831 // Normalization
832 fNorm = new TH1F("fNorm","fNorm",24,0,24);
833 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
834
835 // Branching ratio
836 fBR = new TH1F("fBR","fBR",24,0,24);
837 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
838
839 // Tracking efficiency
840 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
841 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 3% (1% per track)
842
843
844 // Raw yield extraction
845 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
846 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
847
848 // Cuts efficiency (from cuts variation)
849 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
850 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
851
852 // PID efficiency (from PID/noPID)
853 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
854 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
855
856
857 // MC dN/dpt (copied from D0 : will update later)
858 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 859 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
860 fMCPtShape->SetBinContent(1,0.03);
861 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 862
5541b811 863 return;
fa2be30c 864
5541b811 865}
866//--------------------------------------------------------------------------
0233abe6 867Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
868 //
869 // Get error
870 //
871
872 Int_t bin=fCutsEff->FindBin(pt);
873
874 return fCutsEff->GetBinContent(bin);
875}
876//--------------------------------------------------------------------------
877Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
878 //
879 // Get error
880 //
881
882 Int_t bin=fMCPtShape->FindBin(pt);
883
884 return fMCPtShape->GetBinContent(bin);
885}
886//--------------------------------------------------------------------------
887Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
888 //
889 // Get error
890 //
891
892 Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
893
894 return TMath::Sqrt(err);
895}
896//--------------------------------------------------------------------------
897Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
898 //
899 // Get error
900 //
901
902 Int_t bin=fPIDEff->FindBin(pt);
903
904 return fPIDEff->GetBinContent(bin);
905}
906//--------------------------------------------------------------------------
907Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
908 //
909 // Get error
910 //
911
912 Int_t bin=fTrackingEff->FindBin(pt);
913
914 return fTrackingEff->GetBinContent(bin);
915}
916//--------------------------------------------------------------------------
917Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
918 //
919 // Get error
920 //
921
922 Int_t bin=fRawYield->FindBin(pt);
923
924 return fRawYield->GetBinContent(bin);
925}
926//--------------------------------------------------------------------------
927Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
928 //
929 // Get error
930 //
931
932 Int_t bin=fPartAntipart->FindBin(pt);
933
934 return fPartAntipart->GetBinContent(bin);
935}
936//--------------------------------------------------------------------------
937Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
938 //
939 // Get total syst error (except norm. error)
940 //
941
942 Double_t err=0.;
943
944 if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
945 if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
946 if(fBR) err += GetBRErr()*GetBRErr();
947 if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
948 if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
949 if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
950 if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
951
952 err += feeddownErr*feeddownErr;
953
954 return TMath::Sqrt(err);
955}
956//---------------------------------------------------------------------------
5ad8ad40 957void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
0233abe6 958 //
959 // Draw errors
960 //
961 gStyle->SetOptStat(0);
962
a07118d0 963 TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
964 cSystErr->Range(0.20,-0.5,18.4,0.34);
965 cSystErr->SetRightMargin(0.318);
0233abe6 966 cSystErr->SetFillColor(0);
967
2ee3afe2 968 TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} [GeV/c]; Relative Error",20,0,20,100,-1,+1);
f1d72175 969 hFrame->SetAxisRange(1.,15.9,"X");
a07118d0 970 hFrame->SetAxisRange(-0.5,0.5,"Y");
0233abe6 971 hFrame->Draw();
972
a07118d0 973 TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
974 leg->SetTextSize(0.03601695);
0233abe6 975 leg->SetFillStyle(0);
976 leg->SetBorderSize(0);
2ee3afe2 977
f1d72175 978 TH1F *hTotErr=new TH1F("hTotErr","",24,0,24);
2ee3afe2 979 Int_t nbins = fNorm->GetNbinsX();
980 TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
f1d72175 981 for(Int_t i=1;i<=24;i++) {
2ee3afe2 982 Double_t pt = hTotErr->GetBinCenter(i);
983 Double_t ptwidth = hTotErr->GetBinWidth(i);
984
985 if(grErrFeeddown) {
986 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
987 Double_t toterryl=0., toterryh=0.;
988 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
989 grErrFeeddown->GetPoint(j,x,y);
990 errxh = grErrFeeddown->GetErrorXhigh(j);
991 errxl = grErrFeeddown->GetErrorXlow(j);
a07118d0 992 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
2ee3afe2 993 erryh = grErrFeeddown->GetErrorYhigh(j);
994 erryl = grErrFeeddown->GetErrorYlow(j);
995 }
996 }
997 if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
998 else toterryl = GetTotalSystErr(pt);
999 if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
1000 else toterryh = GetTotalSystErr(pt);
1001
1002 hTotErr->SetBinContent(i,toterryh);
1003 gTotErr->SetPoint(i,pt,0.);
1004 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
1005 }
1006 else {
1007 hTotErr->SetBinContent(i,GetTotalSystErr(pt));
1008 gTotErr->SetPoint(i,pt,0.);
1009 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
1010 }
1011
1012 }
f968ef30 1013 gTotErr->SetLineColor(kBlack);
1014 gTotErr->SetFillColor(kRed);
1015 gTotErr->SetFillStyle(3002);
2ee3afe2 1016 gTotErr->Draw("2");
1017 leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
1018// hTotErr->SetLineColor(1);
1019// hTotErr->SetLineWidth(3);
1020// hTotErr->Draw("same");
1021// leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
1022
0233abe6 1023
e373b9f0 1024 fNorm->SetFillColor(1);
1025 fNorm->SetFillStyle(3002);
1026 //fNorm->Draw("same");
1027 //TH1F *hNormRefl = ReflectHisto(fNorm);
1028 //hNormRefl->Draw("same");
1029 leg->AddEntry(fNorm,"Normalization (10%)","");
1030
5ad8ad40 1031 if(grErrFeeddown) {
1032 grErrFeeddown->SetFillColor(kTeal-8);
1033 grErrFeeddown->SetFillStyle(3001);
1034 grErrFeeddown->Draw("2");
1035 leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
1036 }
0233abe6 1037 if(fTrackingEff) {
5ad8ad40 1038 fTrackingEff->SetFillColor(4);
1039 fTrackingEff->SetFillStyle(3006);
0233abe6 1040 fTrackingEff->Draw("same");
f968ef30 1041 TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
1042 hTrackingEffRefl->Draw("same");
0233abe6 1043 leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
1044 }
1045 if(fBR) {
1046 fBR->SetFillColor(6);
1047 fBR->SetFillStyle(3005);
5ad8ad40 1048 //fBR->SetFillStyle(3020);
0233abe6 1049 fBR->Draw("same");
f968ef30 1050 TH1F *hBRRefl = ReflectHisto(fBR);
1051 hBRRefl->Draw("same");
0233abe6 1052 leg->AddEntry(fBR,"Branching ratio","f");
1053 }
1054 if(fRawYield) {
5ad8ad40 1055 Int_t ci; // for color index setting
1056 ci = TColor::GetColor("#00cc00");
1057 fRawYield->SetLineColor(ci);
1058 // fRawYield->SetLineColor(3);
0233abe6 1059 fRawYield->SetLineWidth(3);
1060 fRawYield->Draw("same");
f968ef30 1061 TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
1062 hRawYieldRefl->Draw("same");
844405bc 1063 leg->AddEntry(fRawYield,"Yield extraction","l");
0233abe6 1064 }
1065 if(fCutsEff) {
1066 fCutsEff->SetLineColor(4);
1067 fCutsEff->SetLineWidth(3);
1068 fCutsEff->Draw("same");
f968ef30 1069 TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
1070 hCutsEffRefl->Draw("same");
0233abe6 1071 leg->AddEntry(fCutsEff,"Cuts efficiency","l");
1072 }
1073 if(fPIDEff) {
1074 fPIDEff->SetLineColor(7);
1075 fPIDEff->SetLineWidth(3);
1076 fPIDEff->Draw("same");
f968ef30 1077 TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
1078 hPIDEffRefl->Draw("same");
0233abe6 1079 leg->AddEntry(fPIDEff,"PID efficiency","l");
1080 }
1081 if(fMCPtShape) {
5ad8ad40 1082 Int_t ci = TColor::GetColor("#9933ff");
1083 fMCPtShape->SetLineColor(ci);
1084 // fMCPtShape->SetLineColor(8);
0233abe6 1085 fMCPtShape->SetLineWidth(3);
1086 fMCPtShape->Draw("same");
f968ef30 1087 TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
1088 hMCPtShapeRefl->Draw("same");
0233abe6 1089 leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
1090 }
1091 if(fPartAntipart) {
5ad8ad40 1092 Int_t ci = TColor::GetColor("#ff6600");
1093 fPartAntipart->SetLineColor(ci);
1094 // fPartAntipart->SetLineColor(9);
0233abe6 1095 fPartAntipart->SetLineWidth(3);
1096 fPartAntipart->Draw("same");
f968ef30 1097 TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
1098 hPartAntipartRefl->Draw("same");
0233abe6 1099 leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
1100 }
0233abe6 1101
0233abe6 1102
1103 leg->Draw();
1104
5ad8ad40 1105 cSystErr->SaveAs("RelativeSystematics.eps");
1106
0233abe6 1107 return;
1108}
f968ef30 1109//-------------------------------------------------------------------------
1110TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
1111 //
1112 // Clones and reflects histogram
1113 //
1114 TH1F *hout=(TH1F*)hin->Clone("hout");
1115 hout->Scale(-1.);
1116
1117 return hout;
1118}