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