]> 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
Merge remote-tracking branch 'origin/master' into flatdev
[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:
53c91ed7 24// AliHFSystEff syst; // DECAY = 1 for D0, 2, for D+, 3 for D*, 5 for Lc->pKp, 6 for Lc->pK0S
5541b811 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();
830baf70 121 else if (fCentralityClass=="3050OutOfPlane")InitD0toKpi2011PbPb3050OutOfPlane();
122 else if (fCentralityClass == "3050")InitD0toKpi2011PbPb3050();
3db225bf 123 else if (fCentralityClass=="010") InitD0toKpi2011PbPb010CentScan();
124 else if (fCentralityClass=="1020") InitD0toKpi2011PbPb1020CentScan();
125 else if (fCentralityClass=="2030") InitD0toKpi2011PbPb2030CentScan();
126 else if (fCentralityClass=="3040") InitD0toKpi2011PbPb3040CentScan();
127 else if (fCentralityClass=="4050") InitD0toKpi2011PbPb4050CentScan();
128 else if (fCentralityClass=="5080") InitD0toKpi2010PbPb5080CentScan();
129 else AliFatal("Not yet implemented");
130 }
9341a8fa 131 }
132 else if (fCollisionType==2) {
86164e35 133 if (fCentralityClass=="0100") InitD0toKpi2013pPb0100();
fa90243b 134
135 if (fCentralityClass=="020V0A") InitD0toKpi2013pPb020V0A();
136 if (fCentralityClass=="2040V0A") InitD0toKpi2013pPb2040V0A();
137 if (fCentralityClass=="4060V0A") InitD0toKpi2013pPb4060V0A();
138 if (fCentralityClass=="60100V0A") InitD0toKpi2013pPb60100V0A();
139
140 if (fCentralityClass=="020ZNA") InitD0toKpi2013pPb020ZNA();
141 if (fCentralityClass=="2040ZNA") InitD0toKpi2013pPb2040ZNA();
142 if (fCentralityClass=="4060ZNA") InitD0toKpi2013pPb4060ZNA();
143 if (fCentralityClass=="60100ZNA") InitD0toKpi2013pPb60100ZNA();
144
44dd56b2 145 if (fIsRapidityScan) {
146 if (fRapidityRange == "0804") InitD0toKpi2013pPb0100RapScan0804();
147 if (fRapidityRange == "0401") InitD0toKpi2013pPb0100RapScan0401();
148 if (fRapidityRange == "0101") InitD0toKpi2013pPb0100RapScan0101();
149 if (fRapidityRange == "0104") InitD0toKpi2013pPb0100RapScan0104();
150 if (fRapidityRange == "0408") InitD0toKpi2013pPb0100RapScan0408();
151 }
9341a8fa 152 }
153 else AliFatal("Not yet implemented");
0233abe6 154 break;
3db225bf 155
0233abe6 156 case 2: // D+->Kpipi
a8c3295d 157 if (fCollisionType==0) {
158 if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
159 else InitDplustoKpipi2010pp();
9341a8fa 160 }
161 else if (fCollisionType==1) {
3db225bf 162 if (fRunNumber == 10){
163 if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
164 else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
165 else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
166 else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
167 else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
168 else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
169 else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
170 else AliFatal("Not yet implemented");
171 }
172 if(fRunNumber == 11){
173 if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
174 else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010CentScan();
175 else if (fCentralityClass=="1020") InitDplustoKpipi2011PbPb1020CentScan();
176 else if (fCentralityClass=="2030") InitDplustoKpipi2011PbPb2030CentScan();
177 else if (fCentralityClass=="3040") InitDplustoKpipi2011PbPb3040CentScan();
178 else if (fCentralityClass=="4050") InitDplustoKpipi2011PbPb4050CentScan();
179 else if (fCentralityClass=="5080") InitDplustoKpipi2010PbPb5080CentScan();
830baf70 180 else if (fCentralityClass=="3050") InitDplustoKpipi2011PbPb3050();
3db225bf 181 else AliFatal("Not yet implemented");
182 }
9341a8fa 183 }
184 else if (fCollisionType==2) {
86164e35 185 if (fCentralityClass=="0100") InitDplustoKpipi2013pPb0100();
44dd56b2 186 if (fIsRapidityScan) {
187 if (fRapidityRange == "0804") InitDplustoKpipi2013pPb0100RapScan0804();
188 if (fRapidityRange == "0401") InitDplustoKpipi2013pPb0100RapScan0401();
189 if (fRapidityRange == "0101") InitDplustoKpipi2013pPb0100RapScan0101();
190 if (fRapidityRange == "0104") InitDplustoKpipi2013pPb0100RapScan0104();
191 if (fRapidityRange == "0408") InitDplustoKpipi2013pPb0100RapScan0408();
192 }
9341a8fa 193 }
194 else AliFatal("Not yet implemented");
0233abe6 195 break;
196 case 3: // D*->D0pi
fa2be30c 197 if (fCollisionType==0) {
198 if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
199 else InitDstartoD0pi2010pp();
9341a8fa 200 }
201 else if (fCollisionType==1) {
3db225bf 202 if (fRunNumber == 10){
203 if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
204 else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
205 else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
206 else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
207 else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
208 else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
209 else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
210 else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
211 else AliFatal("Not yet implemented");
212 }
213 if (fRunNumber == 11){
214 if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
215 else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010CentScan();
216 else if (fCentralityClass=="1020") InitDstartoD0pi2011PbPb1020CentScan();
217 else if (fCentralityClass=="2030") InitDstartoD0pi2011PbPb2030CentScan();
218 else if (fCentralityClass=="3040") InitDstartoD0pi2011PbPb3040CentScan();
219 else if (fCentralityClass=="4050") InitDstartoD0pi2011PbPb4050CentScan();
220 else if (fCentralityClass=="5080") InitDstartoD0pi2010PbPb5080CentScan();
830baf70 221 else if (fCentralityClass=="3050") InitDstartoD0pi2011PbPb3050();
3db225bf 222 else AliFatal("Not yet implemented");
223 }
224 }
225 else if (fCollisionType==2) {
86164e35 226 if (fCentralityClass=="0100") InitDstartoD0pi2013pPb0100();
fa90243b 227
228 if (fCentralityClass=="020V0A") InitDstartoD0pi2013pPb020V0A();
229 if (fCentralityClass=="2040V0A") InitDstartoD0pi2013pPb2040V0A();
230 if (fCentralityClass=="4060V0A") InitDstartoD0pi2013pPb4060V0A();
231 if (fCentralityClass=="60100V0A") InitDstartoD0pi2013pPb60100V0A();
232
233 if (fCentralityClass=="020ZNA") InitDstartoD0pi2013pPb020ZNA();
234 if (fCentralityClass=="2040ZNA") InitDstartoD0pi2013pPb2040ZNA();
235 if (fCentralityClass=="4060ZNA") InitDstartoD0pi2013pPb4060ZNA();
236 if (fCentralityClass=="60100ZNA") InitDstartoD0pi2013pPb60100ZNA();
237
44dd56b2 238 if (fIsRapidityScan) {
239 if (fRapidityRange == "0804") InitDstartoD0pi2013pPb0100RapScan0804();
240 if (fRapidityRange == "0401") InitDstartoD0pi2013pPb0100RapScan0401();
241 if (fRapidityRange == "0101") InitDstartoD0pi2013pPb0100RapScan0101();
242 if (fRapidityRange == "0104") InitDstartoD0pi2013pPb0100RapScan0104();
243 if (fRapidityRange == "0408") InitDstartoD0pi2013pPb0100RapScan0408();
244 }
5541b811 245 }
9341a8fa 246 else AliFatal("Not yet implemented");
0233abe6 247 break;
66247513 248 case 4: // D+s->KKpi
249 if (fCollisionType==0) InitDstoKKpi2010pp();
46988859 250 else if (fCollisionType==1) {
251 if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
252 else AliFatal("Not yet implemented");
9341a8fa 253 }
254 else if (fCollisionType==2) {
86164e35 255 if (fCentralityClass=="0100") InitDstoKKpi2013pPb0100();
46988859 256 }
27360a0e 257 else AliFatal("Not yet implemented");
66247513 258 break;
9ca30701 259 case 5: // Lc->pKpi
260 if (fCollisionType==0) InitLctopKpi2010pp();
261 else AliFatal("Not yet implemented");
262 break;
53c91ed7 263 case 6: // Lc->pK0S
264 if (fCollisionType==0) InitLctopK0S2010pp();
265 else AliFatal("Not yet implemented");
266 break;
66247513 267
0233abe6 268 default:
269 printf("Invalid decay type: %d\n",decay);
270 break;
271 }
272
0233abe6 273}
5541b811 274
0233abe6 275//--------------------------------------------------------------------------
5541b811 276void AliHFSystErr::InitD0toKpi2010pp() {
277 //
278 // D0->Kpi syst errors. Responsible: A. Rossi
279 // 2010 pp sample
0233abe6 280 //
281
5541b811 282 // Normalization
f1d72175 283 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 284 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
5541b811 285
286 // Branching ratio
f1d72175 287 fBR = new TH1F("fBR","fBR",24,0,24);
288 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
5541b811 289
290 // Tracking efficiency
f1d72175 291 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
f8bc42d3 292 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
5541b811 293
294 // Raw yield extraction
f1d72175 295 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
5541b811 296 fRawYield->SetBinContent(1,1);
d2e5959f 297 fRawYield->SetBinContent(2,0.22);
298 fRawYield->SetBinContent(3,0.1);
299 for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
300 for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
f1d72175 301 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
302 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
0233abe6 303
5541b811 304 // Cuts efficiency (from cuts variation)
f1d72175 305 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
306 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
5541b811 307
308 // PID efficiency (from PID/noPID)
f1d72175 309 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
310 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
dc80fe6b 311 fPIDEff->SetBinContent(2,0.05); // 5%
5541b811 312
313 // MC dN/dpt
f1d72175 314 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 315 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
316 fMCPtShape->SetBinContent(1,0.03);
317 fMCPtShape->SetBinContent(2,0.03);
5541b811 318
319 // particle-antiparticle
f1d72175 320 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
321 // fPartAntipart->SetBinContent(1,1);
322 // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
5541b811 323
324 return;
0233abe6 325}
326//--------------------------------------------------------------------------
5541b811 327void AliHFSystErr::InitD0toKpi2010PbPb020() {
0233abe6 328 //
aebba8df 329 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 330 // 2010 PbPb sample, 0-20 CC
0233abe6 331 //
5541b811 332
333 // Normalization
334 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 335 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
5541b811 336
337 // Branching ratio
338 fBR = new TH1F("fBR","fBR",20,0,20);
339 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
340
341 // Tracking efficiency
342 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 343 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
5541b811 344
345 // Raw yield extraction
346 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
347 fRawYield->SetBinContent(1,0);
348 fRawYield->SetBinContent(2,0);
6c388dee 349 fRawYield->SetBinContent(3,0.08);
350 for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
351 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
5541b811 352
353 // Cuts efficiency (from cuts variation)
354 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 355 fCutsEff->SetBinContent(1,0.);
356 fCutsEff->SetBinContent(2,0.);
cc4bb301 357 fCutsEff->SetBinContent(3,0.13);
358 fCutsEff->SetBinContent(4,0.11);
359 for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
6c388dee 360 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
5541b811 361
362 // PID efficiency (from PID/noPID)
363 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
6c388dee 364 for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
5541b811 365
366 // MC dN/dpt
367 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
cc4bb301 368 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
369 fMCPtShape->SetBinContent(3,0.04);
370 fMCPtShape->SetBinContent(4,0.02);
b797b51b 371 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
578bfbba 372
6c388dee 373// // particle-antiparticle
374// fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
375// for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
376// fPartAntipart->SetBinContent(3,0.10);
377// fPartAntipart->SetBinContent(4,0.10);
378// fPartAntipart->SetBinContent(7,0.10);
379// fPartAntipart->SetBinContent(8,0.10);
0233abe6 380
5541b811 381 return;
382}
383//--------------------------------------------------------------------------
9a8b5080 384//--------------------------------------------------------------------------
385void AliHFSystErr::InitD0toKpi2011PbPb07half() {
386 //
387 // D0->Kpi syst errors. Responsible: A. Rossi
388 // 2011 PbPb sample, 0-7.5 CC
389 //
390
391 // Normalization
392 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 393 for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
9a8b5080 394
395 // Branching ratio
396 fBR = new TH1F("fBR","fBR",36,0,36);
397 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
398
399 // Tracking efficiency
400 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
401 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
402 for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
403
404 // Raw yield extraction
405 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
406 fRawYield->SetBinContent(1,0);
407 fRawYield->SetBinContent(2,0.1);
408 fRawYield->SetBinContent(3,0.05);
fd8ae434 409 for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
410 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
9a8b5080 411 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
412 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
413
414 // Cuts efficiency (from cuts variation)
415 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
416 fCutsEff->SetBinContent(1,0.);
417 fCutsEff->SetBinContent(2,0.15);
418 fCutsEff->SetBinContent(3,0.13);
419 fCutsEff->SetBinContent(4,0.11);
420 fCutsEff->SetBinContent(5,0.08);
421 for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
422 for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
423
424 // PID efficiency (from PID/noPID)
425 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
426 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
427 for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
428
429 // MC dN/dpt
430 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
431 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
432 fMCPtShape->SetBinContent(2,0.06);
433 fMCPtShape->SetBinContent(3,0.04);
434 fMCPtShape->SetBinContent(4,0.02);
435 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
436 for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
437 for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
438
439// // particle-antiparticle
440// fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
441// for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
442// fPartAntipart->SetBinContent(3,0.10);
443// fPartAntipart->SetBinContent(4,0.10);
444// fPartAntipart->SetBinContent(7,0.10);
445// fPartAntipart->SetBinContent(8,0.10);
446
447 return;
448}
93da8165 449
450//--------------------------------------------------------------------------
451void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
452 //
453 // D0->Kpi syst errors. Responsible: D. Caffarri
454 // 2011 PbPb sample, 30-50 CC InPlane
455 //
456 InitD0toKpi2011PbPb07half();
457 // Raw yield extraction
458 // fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
459 fRawYield->SetBinContent(1,0);
460 fRawYield->SetBinContent(2,0);
461 fRawYield->SetBinContent(3,0.05);
462 fRawYield->SetBinContent(4,0.08);
463 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
464 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
465 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
466 for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
467
468 // Cuts efficiency (from cuts variation)
469 //fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
470 fCutsEff->SetBinContent(1,0.);
471 fCutsEff->SetBinContent(2,0.0);
472 fCutsEff->SetBinContent(3,0.10);
473 fCutsEff->SetBinContent(4,0.10);
474 fCutsEff->SetBinContent(5,0.10);
475 fCutsEff->SetBinContent(6,0.10);
476 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
477 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
478 for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
479}
480
481//--------------------------------------------------------------------------
482void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
483 //
484 // D0->Kpi syst errors. Responsible: D. Caffarri
485 // 2011 PbPb sample, 30-50 CC OutOfPlane
486 //
487 InitD0toKpi2011PbPb07half();
488 // Raw yield extraction
489 //fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
490 fRawYield->SetBinContent(1,0);
491 fRawYield->SetBinContent(2,0.);
492 fRawYield->SetBinContent(3,0.05);
493 for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.07);
494 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
495 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
496 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
497 for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
498
499 // Cuts efficiency (from cuts variation)
500 // fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
501 fCutsEff->SetBinContent(1,0.);
502 fCutsEff->SetBinContent(2,0.);
503 fCutsEff->SetBinContent(3,0.10);
504 fCutsEff->SetBinContent(4,0.10);
505 fCutsEff->SetBinContent(5,0.10);
506 fCutsEff->SetBinContent(6,0.10);
507 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
508 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
509 for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
510}
511
9a8b5080 512//--------------------------------------------------------------------------
5541b811 513void AliHFSystErr::InitD0toKpi2010PbPb4080() {
93da8165 514 //
aebba8df 515 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 516 // 2010 PbPb sample, 40-80 CC
517 //
518
0233abe6 519 // Normalization
520 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 521 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
0233abe6 522
523 // Branching ratio
524 fBR = new TH1F("fBR","fBR",20,0,20);
844405bc 525 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
0233abe6 526
aebba8df 527 // Tracking efficiency
528 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 529 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
aebba8df 530
531
0233abe6 532 // Raw yield extraction
533 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
5541b811 534 fRawYield->SetBinContent(1,0);
535 fRawYield->SetBinContent(2,0);
6c388dee 536 for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
666cf569 537 //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
0233abe6 538
539 // Cuts efficiency (from cuts variation)
540 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 541 fCutsEff->SetBinContent(1,0.);
542 fCutsEff->SetBinContent(2,0.);
cc4bb301 543 fCutsEff->SetBinContent(3,0.13);
544 fCutsEff->SetBinContent(4,0.11);
545 for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
6c388dee 546 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
0233abe6 547
548 // PID efficiency (from PID/noPID)
549 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
6c388dee 550// for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
551// for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
552 for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
0233abe6 553
554 // MC dN/dpt
555 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
cc4bb301 556 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
557 // fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
558 // fMCPtShape->SetBinContent(4,0.02);
559 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
560
561 // // particle-antiparticle
562 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
563 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
0233abe6 564
565 return;
566}
a8c3295d 567
830baf70 568//--------------------------------------------------------------------------
569void AliHFSystErr::InitD0toKpi2011PbPb3050() {
570 //
571 // D0->Kpi syst errors. Responsible: D. Caffarri
572 // 2011 PbPb sample, 30-50 CC
573 //
574
575 // Normalization
576 fNorm = new TH1F("fNorm","fNorm",20,0,20);
577 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
578
579 // Branching ratio
580 fBR = new TH1F("fBR","fBR",20,0,20);
581 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
582
583 // Tracking efficiency
584 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
585 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
586
587
588 // Raw yield extraction
589 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
590 fRawYield->SetBinContent(1,0);
591 fRawYield->SetBinContent(2,0.102);
592 fRawYield->SetBinContent(3,0.071);
593 fRawYield->SetBinContent(4,0.052);
594 fRawYield->SetBinContent(5,0.063);
595 fRawYield->SetBinContent(6,0.071);
596 fRawYield->SetBinContent(7,0.071);
597 fRawYield->SetBinContent(8,0.071);
598 fRawYield->SetBinContent(9,0.071);
599 fRawYield->SetBinContent(10,0.071);
600 fRawYield->SetBinContent(11,0.071);
601 fRawYield->SetBinContent(12,0.071);
602
603 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.112);
604 for(Int_t i=17;i<=20;i++) fRawYield->SetBinContent(i,0);
605
606 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
607 fCutsEff->SetBinContent(1,0.);
608 fCutsEff->SetBinContent(2,0.10);
609 fCutsEff->SetBinContent(3,0.10);
610 fCutsEff->SetBinContent(4,0.10);
611 fCutsEff->SetBinContent(5,0.10);
612 fCutsEff->SetBinContent(6,0.10);
613 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
614 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
615
616
617 // PID efficiency (from PID/noPID)
618 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
619 // for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
620 // for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
621 for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
622
623 // MC dN/dpt
624 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",0,0,20);
625 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
626 fMCPtShape->SetBinContent(2,0.06);
627 fMCPtShape->SetBinContent(3,0.04);
628 fMCPtShape->SetBinContent(4,0.02);
629 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
630 for(Int_t i=17;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
631 // // particle-antiparticle
632 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
633 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
634
635 return;
636}
637
638
a8c3295d 639//--------------------------------------------------------------------------
640void AliHFSystErr::InitD0toKpi2010ppLowEn() {
641 //
642 // D0->Kpi syst errors. Low energy run
643 // 2011 2.76 TeV pp sample
644 //
645 AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
646
647 // Normalization
648 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 649 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 4% error on sigmaV0and
a8c3295d 650
651 // Branching ratio
652 fBR = new TH1F("fBR","fBR",20,0,20);
653 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
654
655 // Tracking efficiency
656 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 657 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
a8c3295d 658
659 // Raw yield extraction
660 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
661 fRawYield->SetBinContent(1,1);
f9a85c82 662 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
a8c3295d 663
664 // Cuts efficiency (from cuts variation)
665 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
6c388dee 666 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
f9a85c82 667 fCutsEff->SetBinContent(2,0.20);
668 for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
669
a8c3295d 670
671 // PID efficiency (from PID/noPID)
672 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
f9a85c82 673 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
674 // fPIDEff->SetBinContent(2,0.20);
675 for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
a8c3295d 676
677 // MC dN/dpt
678 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 679 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
680 fMCPtShape->SetBinContent(1,0.03);
681 fMCPtShape->SetBinContent(2,0.03);
a8c3295d 682
f9a85c82 683// // particle-antiparticle
684// fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
685// fPartAntipart->SetBinContent(1,1);
686// fPartAntipart->SetBinContent(2,1);
687// for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
688// for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
a8c3295d 689
690 return;
691}
86164e35 692//_________________________________________________________________________
693void AliHFSystErr::InitD0toKpi2013pPb0100(){
bd332272 694 //
695 // D0->Kpi syst errors. p-Pb data sample
696 //
697 AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
86164e35 698
699 fNorm = new TH1F("fNorm","fNorm",24,0,24);
9341a8fa 700 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
86164e35 701
702 // Branching ratio
9341a8fa 703 fBR = new TH1F("fBR","fBR",24,0,24);
704 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.2% PDG2010
86164e35 705
706 // Tracking efficiency
707 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
708 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track
709
710 // Raw yield extraction
711 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
9341a8fa 712 fRawYield->SetBinContent(1,0.0);
bc95efa5 713 fRawYield->SetBinContent(2,0.08);
714 for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.04);
715 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.07);
716 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.11);
86164e35 717
718 // Cuts efficiency (from cuts variation)
719 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
9341a8fa 720 fCutsEff->SetBinContent(1,0.0);
e5f8509c 721 fCutsEff->SetBinContent(2,0.08);
9341a8fa 722 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
86164e35 723
724 // PID efficiency (from PID/noPID)
725 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
bc95efa5 726 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
86164e35 727
728 // MC dN/dpt
729 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
9341a8fa 730 fMCPtShape->SetBinContent(1,0.0);
731 fMCPtShape->SetBinContent(2,0.02);
e5f8509c 732 for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
86164e35 733
734 // particle-antiparticle
735 /*
736 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
9341a8fa 737 for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
86164e35 738 */
739 return;
740
741
742}
a8c3295d 743
44dd56b2 744
0233abe6 745//--------------------------------------------------------------------------
5541b811 746void AliHFSystErr::InitDplustoKpipi2010pp() {
0233abe6 747 //
748 // D+->Kpipi syst errors. Responsible: R. Bala
5541b811 749 // 2010 pp sample
0233abe6 750 //
751
1f91ef56 752
753// Normalization
754 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 755 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
626bcf57 756
757 // Branching ratio
1f91ef56 758 fBR = new TH1F("fBR","fBR",24,0,24);
f7c8a2f1 759 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
626bcf57 760
761 // Tracking efficiency
1f91ef56 762 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
f8bc42d3 763 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
626bcf57 764
765
766 // Raw yield extraction
1f91ef56 767 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
d2e5959f 768 fRawYield->SetBinContent(1,1);
1f91ef56 769 fRawYield->SetBinContent(2,0.25);
770 fRawYield->SetBinContent(3,0.25);
d2e5959f 771 fRawYield->SetBinContent(4,0.20);
1f91ef56 772 fRawYield->SetBinContent(5,0.09);
773 fRawYield->SetBinContent(6,0.09);
d2e5959f 774 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
775 for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
1f91ef56 776
626bcf57 777 // Cuts efficiency (from cuts variation)
1f91ef56 778 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
779 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
626bcf57 780
781 // PID efficiency (from PID/noPID)
1f91ef56 782 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
783 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
784 fPIDEff->SetBinContent(1,0.15); // 15%
785 fPIDEff->SetBinContent(2,0.15); // 15%
786 fPIDEff->SetBinContent(3,0.15); // 15%
787 fPIDEff->SetBinContent(4,0.15); // 15%
788 for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%
789
790 // MC dN/dpt
791 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
792 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
578bfbba 793 fMCPtShape->SetBinContent(1,0.03);
794 fMCPtShape->SetBinContent(2,0.03);
626bcf57 795
796
797 // particle-antiparticle
798 /*
799 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
800 fPartAntipart->SetBinContent(1,1);
801 fPartAntipart->SetBinContent(2,1);
802 fPartAntipart->SetBinContent(3,0.12);
803 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
804 */
0233abe6 805 return;
806}
1f91ef56 807
66247513 808//--------------------------------------------------------------------------
809void AliHFSystErr::InitDstoKKpi2010pp() {
810 //
811 // D+s->KKpi syst errors. Responsible: G.M. Innocenti
812 // 2010 pp sample
813 //
814
815
816// Normalization
817 fNorm = new TH1F("fNorm","fNorm",12,0,12);
494408cc 818 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
66247513 819
820 // Branching ratio
821 fBR = new TH1F("fBR","fBR",12,0,12);
494408cc 822 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
66247513 823
824 // Tracking efficiency
825 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
494408cc 826 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
66247513 827
828
829 // Raw yield extraction
830 fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
831 fRawYield->SetBinContent(1,1);
832 fRawYield->SetBinContent(2,1);
833 fRawYield->SetBinContent(3,0.20);
834 fRawYield->SetBinContent(4,0.20);
835 fRawYield->SetBinContent(5,0.15);
836 fRawYield->SetBinContent(6,0.15);
837 fRawYield->SetBinContent(7,0.15);
838 fRawYield->SetBinContent(8,0.15);
839 fRawYield->SetBinContent(9,0.20);
840 fRawYield->SetBinContent(10,0.20);
841 fRawYield->SetBinContent(11,0.20);
842 fRawYield->SetBinContent(12,0.20);
843
844 // Cuts efficiency (from cuts variation)
845 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
846 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%
847
848 // PID efficiency (from PID/noPID)
849 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
850 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
851
494408cc 852 // MC dN/dpt
853 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
854 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
855 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
856 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
857 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
858 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
66247513 859
860
861 // particle-antiparticle
862 /*
863 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
864 fPartAntipart->SetBinContent(1,1);
865 fPartAntipart->SetBinContent(2,1);
866 fPartAntipart->SetBinContent(3,0.12);
867 for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
868 */
869 return;
870}
871
872
436b185f 873//--------------------------------------------------------------------------
874void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
875 //
876 // D+->Kpipi syst errors. Responsible: E. Bruna
877 // 2011 PbPb sample, 0-7.5% CC
878 //
879
880 // Normalization
881 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 882 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
436b185f 883
884 // Branching ratio
885 fBR = new TH1F("fBR","fBR",36,0,36);
886 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
887
888 // Tracking efficiency
889 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
890 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
891
892 // Raw yield extraction
893 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
894 for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
eb4886c7 895 fRawYield->SetBinContent(4,0.30);
896 fRawYield->SetBinContent(5,0.20);
897 fRawYield->SetBinContent(6,0.20);
898 for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
899 for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
436b185f 900 for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
eb4886c7 901 for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
436b185f 902 for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
903
904
905 // Cuts efficiency (from cuts variation)
906 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
a7cfa9ed 907 for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
908 for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);
eb4886c7 909
436b185f 910
911 // PID efficiency (from PID/noPID)
912 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
913 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
914
a7cfa9ed 915 // MC dN/dpt (24/7/2012)
436b185f 916 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
3d446ac8 917 for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
918 for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
919 // for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
436b185f 920
921
922 // particle-antiparticle
923 /*
924 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
925 fPartAntipart->SetBinContent(1,1);
926 fPartAntipart->SetBinContent(2,1);
927 fPartAntipart->SetBinContent(3,0.12);
928 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
929 */
930
931 return;
932}
830baf70 933//___________________________________________________________________
934void AliHFSystErr::InitDplustoKpipi2011PbPb3050() {
935 //
936 // D+->Kpipi syst errors. Responsible: R. Bala
937 // 2011 PbPb sample, 30-50% CC
938 //
939 AliInfo(" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
940 // Normalization
941 fNorm = new TH1F("fNorm","fNorm",16,0,16);
942 for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
943
944 // Branching ratio
945 fBR = new TH1F("fBR","fBR",16,0,16);
946 for(Int_t i=1;i<=16;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
947
948 // Tracking efficiency
949 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
950 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
951
952 // Raw yield extraction
953 fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
954 for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.09); //5 to 10%
955 for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.12);
956
957
958 // Cuts efficiency (from cuts variation)
959 fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
960 for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
961 for(Int_t i=12; i<=16; i++) fCutsEff->SetBinContent(i,0.14);
962
963 // PID efficiency (from PID/noPID)
964 fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
965 for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.05); // 5%
966
967 // MC dN/dpt (24/7/2012)
968 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
969 fMCPtShape->SetBinContent(3,0.10);
970 fMCPtShape->SetBinContent(4,0.04);
971 for(Int_t iBin=5; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.02);
972
973
974 // particle-antiparticle
975 /*
976 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
977 fPartAntipart->SetBinContent(1,1);
978 fPartAntipart->SetBinContent(2,1);
979 fPartAntipart->SetBinContent(3,0.12);
980 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
981 */
982
983 return;
984}
985
436b185f 986
46988859 987//--------------------------------------------------------------------------
988void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
989 //
990 // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
991 // 2011 PbPb sample, 0-7.5% CC
992 //
993
994 // Normalization
995 fNorm = new TH1F("fNorm","fNorm",12,0,12);
a713683e 996 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
46988859 997
998 // Branching ratio
999 fBR = new TH1F("fBR","fBR",12,0,12);
d92cd593 1000 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
46988859 1001
1002 // Tracking efficiency
1003 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
1004 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1005
1006 // Raw yield extraction
1007 fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
1008 for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30);
1009 for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
1010
1011
1012
1013 // Cuts efficiency (from cuts variation)
1014 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
1015 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
1016
1017 // PID efficiency (from PID/noPID)
1018 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
1019 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
1020
1021 // MC dN/dpt
1022 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
1023 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
1024 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
1025 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
1026 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
1027 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
1028
1029 // particle-antiparticle
1030 /*
1031 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
1032 fPartAntipart->SetBinContent(1,1);
1033 fPartAntipart->SetBinContent(2,1);
1034 fPartAntipart->SetBinContent(3,0.12);
1035 for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1036 */
1037
1038 return;
1039}
1040
86164e35 1041//_________________________________________________________________________
1042void AliHFSystErr::InitDstoKKpi2013pPb0100(){
bd332272 1043 //
86164e35 1044 // Dstar syst
bd332272 1045 AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
86164e35 1046
1047 //Normalization
1048 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1049 for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); //
1050
1051 // Branching ratio
1052 fBR = new TH1F("fBR","fBR",12,0,12);
1053 for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
1054
1055 // Tracking efficiency
1056 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
59604f68 1057 for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
86164e35 1058
1059 // Raw yield extraction
1060 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
59604f68 1061 for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
1062 for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
1063 for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
1064 for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
1065
86164e35 1066 // Cuts efficiency (from cuts variation)
1067 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
59604f68 1068 for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //
86164e35 1069
1070 // PID efficiency (from PID/noPID)
1071 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
59604f68 1072 for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10);
1073 for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
86164e35 1074
1075 // MC dN/dpt
1076 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
59604f68 1077 for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
86164e35 1078
1079 // particle-antiparticle
1080 /*
1081 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1082 fPartAntipart->SetBinContent(1,1);
1083 fPartAntipart->SetBinContent(2,1);
1084 fPartAntipart->SetBinContent(3,0.12);
1085 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1086 */
1087 return;
86164e35 1088}
1089
1090
0233abe6 1091//--------------------------------------------------------------------------
5541b811 1092void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
1093 //
1094 // D+->Kpipi syst errors. Responsible: ??
1095 // 2010 PbPb sample, 0-20 CC
1096 //
aebba8df 1097
aebba8df 1098 // Normalization
1099 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 1100 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
aebba8df 1101
1102 // Branching ratio
1103 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1104 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
aebba8df 1105
1106 // Tracking efficiency
1107 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 1108 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
aebba8df 1109
1110 // Raw yield extraction
1111 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
6c388dee 1112 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
1113 // fRawYield->SetBinContent(5,0.23);
1114 //fRawYield->SetBinContent(6,0.23);
cb142b90 1115 fRawYield->SetBinContent(7,0.20);
1116 fRawYield->SetBinContent(8,0.20);
1117 fRawYield->SetBinContent(9,0.15);
1118 fRawYield->SetBinContent(10,0.15);
1119 fRawYield->SetBinContent(11,0.15);
1120 fRawYield->SetBinContent(12,0.15);
aebba8df 1121
1122 // Cuts efficiency (from cuts variation)
1123 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
666cf569 1124 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
aebba8df 1125
1126 // PID efficiency (from PID/noPID)
1127 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1128 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
aebba8df 1129
f7c8a2f1 1130 // MC dN/dpt (2/2/2012)
aebba8df 1131 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 1132 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
f7c8a2f1 1133 for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1134 for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1135 for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
aebba8df 1136
1137
1138 // particle-antiparticle
1139 /*
1140 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1141 fPartAntipart->SetBinContent(1,1);
1142 fPartAntipart->SetBinContent(2,1);
1143 fPartAntipart->SetBinContent(3,0.12);
1144 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1145 */
6c388dee 1146
5541b811 1147 return;
1148}
1149
1150//--------------------------------------------------------------------------
1151void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
1152 //
1153 // D+->Kpipi syst errors. Responsible: ??
1154 // 2010 PbPb sample, 40-80 CC
1155 //
1156
aebba8df 1157
1158 // Normalization
1159 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 1160 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
aebba8df 1161
1162 // Branching ratio
1163 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1164 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
aebba8df 1165
1166 // Tracking efficiency
1167 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 1168 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
aebba8df 1169
1170
1171 // Raw yield extraction
1172 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1173 fRawYield->SetBinContent(1,1);
1174 fRawYield->SetBinContent(2,1);
1175 fRawYield->SetBinContent(3,1);
6970b225 1176 fRawYield->SetBinContent(4,0.15);
cb142b90 1177 fRawYield->SetBinContent(5,0.05);
1178 fRawYield->SetBinContent(6,0.05);
1179 fRawYield->SetBinContent(7,0.15);
1180 fRawYield->SetBinContent(8,0.15);
6970b225 1181 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
aebba8df 1182 for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
1183
1184 // Cuts efficiency (from cuts variation)
1185 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1186 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1187
1188 // PID efficiency (from PID/noPID)
1189 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1190 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1191 fPIDEff->SetBinContent(3,0.13); // 13%
1192
1193
f7c8a2f1 1194 // MC dN/dpt (2/2/2012)
aebba8df 1195 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 1196 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
f7c8a2f1 1197 for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1198 for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1199 for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
aebba8df 1200
1201
1202 // particle-antiparticle
1203 /*
1204 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1205 fPartAntipart->SetBinContent(1,1);
1206 fPartAntipart->SetBinContent(2,1);
1207 fPartAntipart->SetBinContent(3,0.12);
1208 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1209 */
5541b811 1210 return;
1211}
1212
a8c3295d 1213//--------------------------------------------------------------------------
1214void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
1215
1216 //
1217 // D+->Kpipi syst errors. Responsible: R. Bala
1218 // 2011 2.76 TeV pp sample
1219 //
1220 AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1221
1222 // Normalization
1223 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 1224 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
a8c3295d 1225
1226 // Branching ratio
1227 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1228 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
a8c3295d 1229
1230 // Tracking efficiency
1231 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 1232 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
a8c3295d 1233
1234 // Raw yield extraction
1235 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1236 fRawYield->SetBinContent(1,1);
1237 fRawYield->SetBinContent(2,1);
f9a85c82 1238 for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10); //5 to 10%
1239 fRawYield->SetBinContent(7,0.15);
1240 fRawYield->SetBinContent(8,0.15);
1241 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055); //5 to 10%
a8c3295d 1242
1243 // Cuts efficiency (from cuts variation)
1244 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
f9a85c82 1245 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
a8c3295d 1246
1247 // PID efficiency (from PID/noPID)
1248 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1249 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
f9a85c82 1250 fPIDEff->SetBinContent(3,0.10); // 13%
1251 fPIDEff->SetBinContent(4,0.10); // 13%
a8c3295d 1252
1253 // MC dN/dpt (copied from D0 : will update later)
1254 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 1255 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
578bfbba 1256 fMCPtShape->SetBinContent(1,0.03);
1257 fMCPtShape->SetBinContent(2,0.03);
a8c3295d 1258
1259 return;
1260}
1261
86164e35 1262//_________________________________________________________________________
1263void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
bd332272 1264 //
1265 // Dplus syst
1266 AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
86164e35 1267
1268 fNorm = new TH1F("fNorm","fNorm",24,0,24);
3e92d3ee 1269 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
86164e35 1270
2d325332 1271 fBR = new TH1F("fBR","fBR",20,0,24);
1272 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
86164e35 1273
1274 // Tracking efficiency
1275 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1276 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
4c1e4a80 1277
86164e35 1278 // Raw yield extraction
1279 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3e92d3ee 1280 fRawYield->SetBinContent(1,0.);
1281 fRawYield->SetBinContent(2,0.1);
1282 fRawYield->SetBinContent(3,0.08);
1283 fRawYield->SetBinContent(4,0.05);
4c1e4a80 1284 fRawYield->SetBinContent(5,0.05);
3e92d3ee 1285 fRawYield->SetBinContent(6,0.05);
1286 fRawYield->SetBinContent(7,0.05);
1287 fRawYield->SetBinContent(8,0.05);
1288 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1289 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
1290 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
4c1e4a80 1291
86164e35 1292 // Cuts efficiency (from cuts variation)
1293 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3e92d3ee 1294 fCutsEff->SetBinContent(1,0.);
1295 fCutsEff->SetBinContent(2,0.1);
1296 fCutsEff->SetBinContent(3,0.1);
1297 fCutsEff->SetBinContent(4,0.06);
1298 fCutsEff->SetBinContent(5,0.06);
1299 fCutsEff->SetBinContent(6,0.06);
1300 fCutsEff->SetBinContent(7,0.05);
1301 fCutsEff->SetBinContent(8,0.05);
1302 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
1303 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1304 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
86164e35 1305
1306 // PID efficiency (from PID/noPID)
1307 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4c1e4a80 1308 fPIDEff->SetBinContent(1,1);
3e92d3ee 1309 fPIDEff->SetBinContent(2,0.1);
1310 for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0); //
86164e35 1311
4c1e4a80 1312 // MC dN/dpt
86164e35 1313 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3e92d3ee 1314 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
86164e35 1315
1316 // particle-antiparticle
1317 /*
1318 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1319 fPartAntipart->SetBinContent(1,1);
1320 fPartAntipart->SetBinContent(2,1);
1321 fPartAntipart->SetBinContent(3,0.12);
1322 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1323 */
1324 return;
1325
1326
1327}
1328
5541b811 1329//--------------------------------------------------------------------------
1330void AliHFSystErr::InitDstartoD0pi2010pp() {
0233abe6 1331 //
2cf98d0d 1332 // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
5541b811 1333 // 2010 pp sample
0233abe6 1334 //
1335
2cf98d0d 1336 // Normalization
fa2be30c 1337 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 1338 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
fa2be30c 1339
1340 // Branching ratio
1341 fBR = new TH1F("fBR","fBR",24,0,24);
1342 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1343
1344 // Tracking efficiency
1345 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1346 fTrackingEff->SetBinContent(1,1.0);
f8bc42d3 1347 fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
1348 fTrackingEff->SetBinContent(3,0.12);
1349 fTrackingEff->SetBinContent(3,0.12);
1350 for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
fa2be30c 1351
1352
1353 // Raw yield extraction
1354 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1355 fRawYield->SetBinContent(1,1.0);
d2e5959f 1356 fRawYield->SetBinContent(2,0.10);
fa2be30c 1357 fRawYield->SetBinContent(3,0.04);
d2e5959f 1358 fRawYield->SetBinContent(4,0.03);
fa2be30c 1359 fRawYield->SetBinContent(5,0.03);
d2e5959f 1360 fRawYield->SetBinContent(6,0.05);
1361 fRawYield->SetBinContent(7,0.05);
1362 fRawYield->SetBinContent(8,0.05);
1363 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04); //4%
1364 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09); //4%
1365 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
fa2be30c 1366
1367 // Cuts efficiency (from cuts variation)
1368 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1369 fCutsEff->SetBinContent(2,0.22);
1370 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1371
1372 // PID efficiency (from PID/noPID)
1373 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1374 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1375
1376
1377 // MC dN/dpt (copied from D0 : will update later)
1378 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 1379 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1380 fMCPtShape->SetBinContent(1,0.03);
1381 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 1382
1383 return;
1384
1385
1386}
1387//--------------------------------------------------------------------------
1388void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
1389
1390 //
1391 // D+->Kpipi syst errors. Responsible: A. Grelli
1392 // 2011 2.76 TeV pp sample
1393 //
1394 AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1395
1396// Normalization
2cf98d0d 1397 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 1398 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
2cf98d0d 1399
1400 // Branching ratio
1401 fBR = new TH1F("fBR","fBR",20,0,20);
8accac87 1402 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
2cf98d0d 1403
1404 // Tracking efficiency
1405 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 1406 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
2cf98d0d 1407
1408 // Raw yield extraction
1409 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
8accac87 1410 fRawYield->SetBinContent(1,1);
fa2be30c 1411 fRawYield->SetBinContent(2,1);
f9a85c82 1412 fRawYield->SetBinContent(3,0.14);
1413 fRawYield->SetBinContent(4,0.14);
1414 fRawYield->SetBinContent(5,0.12);
1415 fRawYield->SetBinContent(6,0.12);
1416 fRawYield->SetBinContent(7,0.06);
1417 fRawYield->SetBinContent(8,0.06);
1418 fRawYield->SetBinContent(9,0.08);
1419 fRawYield->SetBinContent(10,0.08);
1420 fRawYield->SetBinContent(11,0.08);
1421 fRawYield->SetBinContent(12,0.08);
fa2be30c 1422 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
2cf98d0d 1423
1424 // Cuts efficiency (from cuts variation)
1425 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
f9a85c82 1426 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
1427 fCutsEff->SetBinContent(3,0.15);
1428 fCutsEff->SetBinContent(4,0.15);
1429 fCutsEff->SetBinContent(5,0.15);
1430 fCutsEff->SetBinContent(6,0.15);
1431 fCutsEff->SetBinContent(7,0.10);
1432 fCutsEff->SetBinContent(8,0.10);
1433 fCutsEff->SetBinContent(9,0.10);
1434 fCutsEff->SetBinContent(10,0.10);
1435 fCutsEff->SetBinContent(11,0.10);
1436 fCutsEff->SetBinContent(12,0.10);
2cf98d0d 1437
1438 // PID efficiency (from PID/noPID)
1439 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
fa2be30c 1440 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
2cf98d0d 1441
fa2be30c 1442 // MC dN/dpt
2cf98d0d 1443 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 1444 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
578bfbba 1445 fMCPtShape->SetBinContent(1,0.03);
1446 fMCPtShape->SetBinContent(2,0.03);
2cf98d0d 1447
0233abe6 1448 return;
1449}
6c388dee 1450
fa2be30c 1451//------------------------------------------------------------------------
5541b811 1452void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
1453 //
fa2be30c 1454 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 1455 // 2010 PbPb sample, 0-20 CC
1456 //
1457
fa2be30c 1458 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1459
1460 // Normalization
1461 fNorm = new TH1F("fNorm","fNorm",24,0,24);
d5e15690 1462 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
fa2be30c 1463
1464 // Branching ratio
1465 fBR = new TH1F("fBR","fBR",24,0,24);
1466 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1467
1468 // Tracking efficiency
1469 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1470 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1471
1472
1473 // Raw yield extraction
1474 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
6c388dee 1475 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
1476 fRawYield->SetBinContent(3,0.2);
1477 fRawYield->SetBinContent(4,0.2);
1478 fRawYield->SetBinContent(5,0.2);
1479 fRawYield->SetBinContent(6,0.2);
fa2be30c 1480
1481 // Cuts efficiency (from cuts variation)
1482 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1483 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
6fd587a4 1484 fCutsEff->SetBinContent(4,0.15);
1485 fCutsEff->SetBinContent(5,0.15);
1486 fCutsEff->SetBinContent(6,0.15);
fa2be30c 1487
1488 // PID efficiency (from PID/noPID)
1489 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
6c388dee 1490 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
fa2be30c 1491
1492
6fd587a4 1493 // MC dN/dpt (from study on D* pt shape)
fa2be30c 1494 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
6fd587a4 1495 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
1496 fMCPtShape->SetBinContent(4,0.025);
1497 fMCPtShape->SetBinContent(5,0.025);
1498 fMCPtShape->SetBinContent(6,0.025);
1499 fMCPtShape->SetBinContent(7,0.04);
1500 fMCPtShape->SetBinContent(8,0.04);
1501 fMCPtShape->SetBinContent(9,0.03);
1502 fMCPtShape->SetBinContent(10,0.03);
1503 fMCPtShape->SetBinContent(11,0.03);
1504 fMCPtShape->SetBinContent(12,0.03);
1505
9a8b5080 1506 return;
1507
1508}
1509// ----------------------------- 2011 ----------------------------------
1510void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
1511 //
1512 // D*+->D0pi syst errors. Responsible: A. Grelli
1513 // 2011 PbPb sample, 0-7.5 CC
1514 //
1515
1516 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1517
1518
1519 // Normalization
1520 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 1521 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
9a8b5080 1522
1523 // Branching ratio
1524 fBR = new TH1F("fBR","fBR",36,0,36);
1525 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1526
1527 // Tracking efficiency
1528 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1529 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1530
1531
1532 // Raw yield extraction
1533 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1534 for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05); //4%
ada850b0 1535 fRawYield->SetBinContent(4,0.2);
9a8b5080 1536 fRawYield->SetBinContent(5,0.10);
ada850b0 1537 fRawYield->SetBinContent(6,0.10);
1538 fRawYield->SetBinContent(7,0.08);
9a8b5080 1539 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15); //4%
1540
9a8b5080 1541 // Cuts efficiency (from cuts variation)
1542 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1543 for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
9a8b5080 1544
1545 // PID efficiency (from PID/noPID)
1546 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1547 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
ada850b0 1548 fPIDEff->SetBinContent(4,0.09);
9a8b5080 1549
1550 // MC dN/dpt (from study on D* pt shape)
1551 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
273a7638 1552 for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
1553 fMCPtShape->SetBinContent(4,0.015);
1554 fMCPtShape->SetBinContent(5,0.015);
1555 fMCPtShape->SetBinContent(6,0.015);
1556 fMCPtShape->SetBinContent(7,0.02);
1557 fMCPtShape->SetBinContent(8,0.02);
9a8b5080 1558 fMCPtShape->SetBinContent(9,0.03);
1559 fMCPtShape->SetBinContent(10,0.03);
1560 fMCPtShape->SetBinContent(11,0.03);
1561 fMCPtShape->SetBinContent(12,0.03);
1562
6fd587a4 1563
fa2be30c 1564
1565 return;
1566
1567}
1568//-------------------------------------------------------------------------
1569void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
1570 //
1571 // D*+->D0pi syst errors. Responsible: A. Grelli
1572 // 2010 PbPb sample, 20-40 CC
1573 //
1574
1575 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1576
1577 // Normalization
1578 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1579 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
1580
1581 // Branching ratio
1582 fBR = new TH1F("fBR","fBR",24,0,24);
1583 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1584
1585 // Tracking efficiency
1586 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1587 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1588
1589
1590 // Raw yield extraction
1591 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1592 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
1593
1594 // Cuts efficiency (from cuts variation)
1595 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1596 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1597
1598 // PID efficiency (from PID/noPID)
1599 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1600 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1601
1602
1603 // MC dN/dpt (copied from D0 : will update later)
1604 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 1605 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
1606 fMCPtShape->SetBinContent(1,0.03);
1607 fMCPtShape->SetBinContent(2,0.03);
fa2be30c 1608
5541b811 1609 return;
fa2be30c 1610
5541b811 1611}
6c388dee 1612
5541b811 1613//--------------------------------------------------------------------------
1614void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
1615 //
fa2be30c 1616 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 1617 // 2010 PbPb sample, 40-80 CC
1618 //
1619
fa2be30c 1620 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1621
1622 // Normalization
1623 fNorm = new TH1F("fNorm","fNorm",24,0,24);
d5e15690 1624 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
fa2be30c 1625
1626 // Branching ratio
1627 fBR = new TH1F("fBR","fBR",24,0,24);
1628 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1629
1630 // Tracking efficiency
1631 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1632 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1633
1634
1635 // Raw yield extraction
1636 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
6c388dee 1637 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
1638 fRawYield->SetBinContent(1,1);
1639 fRawYield->SetBinContent(2,0.15);
1640 fRawYield->SetBinContent(3,0.15);
1641 fRawYield->SetBinContent(4,0.15);
1642 fRawYield->SetBinContent(5,0.15);
1643 fRawYield->SetBinContent(6,0.10);
1644 fRawYield->SetBinContent(7,0.10);
1645 fRawYield->SetBinContent(8,0.10);
1646 fRawYield->SetBinContent(9,0.11);
1647 fRawYield->SetBinContent(10,0.11);
1648 fRawYield->SetBinContent(11,0.11);
1649 fRawYield->SetBinContent(12,0.11);
1650 fRawYield->SetBinContent(13,0.08);
1651 fRawYield->SetBinContent(14,0.08);
1652 fRawYield->SetBinContent(15,0.08);
1653 fRawYield->SetBinContent(16,0.08);
1654
1655
fa2be30c 1656 // Cuts efficiency (from cuts variation)
1657 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1658 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1659
1660 // PID efficiency (from PID/noPID)
1661 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
6c388dee 1662 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
fa2be30c 1663
1664
1665 // MC dN/dpt (copied from D0 : will update later)
1666 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
6fd587a4 1667 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1668 fMCPtShape->SetBinContent(2,0.05);
1669 fMCPtShape->SetBinContent(3,0.05);
1670 fMCPtShape->SetBinContent(4,0.05);
1671 fMCPtShape->SetBinContent(5,0.04);
1672 fMCPtShape->SetBinContent(6,0.02);
1673 fMCPtShape->SetBinContent(7,0.04);
1674 fMCPtShape->SetBinContent(8,0.04);
1675
5541b811 1676 return;
fa2be30c 1677
5541b811 1678}
6c388dee 1679
830baf70 1680//----------------------------------------------------------------------
1681void AliHFSystErr::InitDstartoD0pi2011PbPb3050() {
1682 //
1683 // D*+->D0pi syst errors. Responsible: A. Grelli
1684 // 2011 PbPb sample, 30-50 CC
1685 //
1686
1687 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
1688
1689 // Normalization
1690 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1691 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1692
1693 // Branching ratio
1694 fBR = new TH1F("fBR","fBR",24,0,24);
1695 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1696
1697 // Tracking efficiency
1698 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1699 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1700
1701
1702 // Raw yield extraction
1703 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1704 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
1705 fRawYield->SetBinContent(1,0.12);
1706 fRawYield->SetBinContent(2,0.12);
1707 fRawYield->SetBinContent(3,0.12);
1708 fRawYield->SetBinContent(4,0.06);
1709 fRawYield->SetBinContent(5,0.06);
1710 fRawYield->SetBinContent(6,0.05);
1711 fRawYield->SetBinContent(7,0.08);
1712 fRawYield->SetBinContent(8,0.05);
1713 fRawYield->SetBinContent(9,0.05);
1714 fRawYield->SetBinContent(10,0.05);
1715 fRawYield->SetBinContent(11,0.05);
1716 fRawYield->SetBinContent(12,0.05);
1717 fRawYield->SetBinContent(13,0.05);
1718 fRawYield->SetBinContent(14,0.05);
1719 fRawYield->SetBinContent(15,0.05);
1720 fRawYield->SetBinContent(16,0.05);
1721
1722
1723 // Cuts efficiency (from cuts variation)
1724 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1725 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
1726 fCutsEff->SetBinContent(1,0.16);
1727 fCutsEff->SetBinContent(2,0.16);
1728 fCutsEff->SetBinContent(3,0.16);
1729 fCutsEff->SetBinContent(4,0.10);
1730 fCutsEff->SetBinContent(5,0.10);
1731 fCutsEff->SetBinContent(6,0.10);
1732 // PID efficiency (from PID/noPID)
1733 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1734 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1735
1736
1737 // MC dN/dpt (copied from D0 : will update later)
1738 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1739 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1740 fMCPtShape->SetBinContent(2,0.11);
1741 fMCPtShape->SetBinContent(3,0.11);
1742 fMCPtShape->SetBinContent(4,0.03);
1743 fMCPtShape->SetBinContent(5,0.02);
1744 fMCPtShape->SetBinContent(6,0.01);
1745 fMCPtShape->SetBinContent(7,0.01);
1746 fMCPtShape->SetBinContent(8,0.01);
1747
1748 return;
1749}
86164e35 1750//_________________________________________________________________________
1751void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
f21cf4a7 1752 // Dstar syst in pPb 2013 MB
1753 //
bd332272 1754 AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
86164e35 1755
86164e35 1756 //Normalization
1757 fNorm = new TH1F("fNorm","fNorm",24,0,24);
fa90243b 1758 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); //
86164e35 1759
f21cf4a7 1760 // Branching ratio
86164e35 1761 fBR = new TH1F("fBR","fBR",24,0,24);
1762 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1763
1764 // Tracking efficiency
1765 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1766 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1767
1768 // Raw yield extraction
1769 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
f21cf4a7 1770 fRawYield->SetBinContent(1,1.0);
30dc63e6 1771 fRawYield->SetBinContent(2,0.08);
1772 fRawYield->SetBinContent(3,0.02);
f21cf4a7 1773 fRawYield->SetBinContent(4,0.02);
30dc63e6 1774 fRawYield->SetBinContent(5,0.05);
1775 fRawYield->SetBinContent(6,0.02);
f21cf4a7 1776 fRawYield->SetBinContent(7,0.04);
1777 fRawYield->SetBinContent(8,0.04);
30dc63e6 1778 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
f21cf4a7 1779 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
1780 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
1781
86164e35 1782 // Cuts efficiency (from cuts variation)
1783 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
f21cf4a7 1784 fCutsEff->SetBinContent(1,0.);
1785 fCutsEff->SetBinContent(2,0.10);
1786 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
86164e35 1787
1788 // PID efficiency (from PID/noPID)
1789 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3e92d3ee 1790 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
86164e35 1791
1792 // MC dN/dpt
1793 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
c9569eb5 1794 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
1795 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
1796 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
86164e35 1797
1798 // particle-antiparticle
1799 /*
1800 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1801 fPartAntipart->SetBinContent(1,1);
1802 fPartAntipart->SetBinContent(2,1);
1803 fPartAntipart->SetBinContent(3,0.12);
1804 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1805 */
fa90243b 1806
1807 return;
86164e35 1808
1809}
1810
414bc6a5 1811//--------------------------------------------------------------------------
1812void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
1813 // define errors for RAA vs. centrality
1814 InitD0toKpi2010PbPb020();
41cfcd39 1815 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 1816 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
41cfcd39 1817 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1818}
1819//--------------------------------------------------------------------------
1820void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
1821 // define errors for RAA vs. centrality
1822 InitD0toKpi2010PbPb020();
41cfcd39 1823 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 1824 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
41cfcd39 1825 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1826}
1827//--------------------------------------------------------------------------
1828void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
1829 // define errors for RAA vs. centrality
41cfcd39 1830 InitD0toKpi2010PbPb4080();
1831 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 1832 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.14);
41cfcd39 1833 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1834}
1835//--------------------------------------------------------------------------
1836void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
41cfcd39 1837 // define errors for RAA vs. centrality
414bc6a5 1838 InitD0toKpi2010PbPb4080();
41cfcd39 1839 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
27360a0e 1840 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.11);
41cfcd39 1841 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1842}
1843//--------------------------------------------------------------------------
1844void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
41cfcd39 1845 // define errors for RAA vs. centrality
414bc6a5 1846 InitD0toKpi2010PbPb4080();
41cfcd39 1847 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
27360a0e 1848 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
41cfcd39 1849 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1850}
1851//--------------------------------------------------------------------------
1852void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
1853 // define errors for RAA vs. centrality
1854 InitDplustoKpipi2010PbPb020();
7b03f1ec 1855 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
eefab142 1856 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
7b03f1ec 1857
414bc6a5 1858}
1859//--------------------------------------------------------------------------
1860void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
1861 // define errors for RAA vs. centrality
1862 InitDplustoKpipi2010PbPb020();
7b03f1ec 1863 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
eefab142 1864 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1865}
1866//--------------------------------------------------------------------------
1867void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
1868 // define errors for RAA vs. centrality
1869 InitDplustoKpipi2010PbPb020();
7b03f1ec 1870 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
eefab142 1871 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
414bc6a5 1872}
1873//--------------------------------------------------------------------------
1874void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
1875 // define errors for RAA vs. centrality
1876 InitDplustoKpipi2010PbPb4080();
7b03f1ec 1877 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
eefab142 1878 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1879}
1880//--------------------------------------------------------------------------
1881void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
1882 // define errors for RAA vs. centrality
1883 InitDplustoKpipi2010PbPb4080();
7b03f1ec 1884 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
eefab142 1885 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
414bc6a5 1886}
1887
1888//--------------------------------------------------------------------------
1889void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
1890 // define errors for RAA vs. centrality
1891 InitDstartoD0pi2010PbPb020();
1892 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
eefab142 1893 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
414bc6a5 1894}
1895//--------------------------------------------------------------------------
1896void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
1897 // define errors for RAA vs. centrality
1898 InitDstartoD0pi2010PbPb020();
1899 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
eefab142 1900 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
414bc6a5 1901}
1902//--------------------------------------------------------------------------
1903void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
1904 // define errors for RAA vs. centrality
1905 InitDstartoD0pi2010PbPb2040();
1906 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
eefab142 1907 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 1908}
1909//--------------------------------------------------------------------------
1910void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
1911 // define errors for RAA vs. centrality
1912 InitDstartoD0pi2010PbPb4080();
1913 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
eefab142 1914 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
414bc6a5 1915}
1916//--------------------------------------------------------------------------
1917void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
1918 // define errors for RAA vs. centrality
1919 InitDstartoD0pi2010PbPb4080();
1920 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
eefab142 1921 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
414bc6a5 1922}
1923
3db225bf 1924//--------------------------------------------------------------------------
1925void AliHFSystErr::InitD0toKpi2011PbPb010CentScan(){
1926 // define errors for RAA vs. centrality 2011
1927 InitD0toKpi2011PbPb07half();
1928
1929}
1930//--------------------------------------------------------------------------
1931void AliHFSystErr::InitD0toKpi2011PbPb1020CentScan(){
1932 // define errors for RAA vs. centrality 2011
1933 InitD0toKpi2011PbPb07half();
1934
1935}
1936//--------------------------------------------------------------------------
1937void AliHFSystErr::InitD0toKpi2011PbPb2030CentScan(){
1938 // define errors for RAA vs. centrality 2011
1939 InitD0toKpi2011PbPb07half();
1940
1941}
1942//--------------------------------------------------------------------------
1943void AliHFSystErr::InitD0toKpi2011PbPb3040CentScan(){
1944 // define errors for RAA vs. centrality 2011
1945 InitD0toKpi2011PbPb07half();
1946
1947}
1948//--------------------------------------------------------------------------
1949void AliHFSystErr::InitD0toKpi2011PbPb4050CentScan(){
1950 // define errors for RAA vs. centrality 2011
1951 InitD0toKpi2011PbPb07half();
1952
1953}
1954//--------------------------------------------------------------------------
1955void AliHFSystErr::InitD0toKpi2010PbPb5080CentScan(){
1956 // define errors for RAA vs. centrality 2011
1957 InitD0toKpi2011PbPb07half();
1958
1959}
2eab2b1d 1960
3db225bf 1961//--------------------------------------------------------------------------
1962void AliHFSystErr::InitDplustoKpipi2011PbPb010CentScan(){
1963 // define errors for RAA vs. centrality 2011
1964 InitDplustoKpipi2011PbPb07half();
2eab2b1d 1965
1966 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.15);
1967 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
1968 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.15);
1969
1970 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1971 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1972 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
1973
1974 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
1975 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
1976 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.10);
3db225bf 1977}
2eab2b1d 1978
3db225bf 1979//--------------------------------------------------------------------------
1980void AliHFSystErr::InitDplustoKpipi2011PbPb1020CentScan(){
1981 // define errors for RAA vs. centrality 2011
1982 InitDplustoKpipi2011PbPb07half();
2eab2b1d 1983
1984 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
1985 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
1986 for(Int_t i=9;i<=17;i++)fRawYield->SetBinContent(i,0.15);
1987
1988 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
1989 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
1990 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
3db225bf 1991
2eab2b1d 1992 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
1993 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
1994 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.15);
3db225bf 1995}
2eab2b1d 1996
3db225bf 1997//--------------------------------------------------------------------------
1998void AliHFSystErr::InitDplustoKpipi2011PbPb2030CentScan(){
1999 // define errors for RAA vs. centrality 2011
2000 InitDplustoKpipi2011PbPb07half();
2001
2eab2b1d 2002 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2003 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2004 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.08);
2005
2006 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2007 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2008 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2009
2010 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2011 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2012 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
3db225bf 2013}
2eab2b1d 2014
3db225bf 2015//--------------------------------------------------------------------------
2016void AliHFSystErr::InitDplustoKpipi2011PbPb3040CentScan(){
2017 // define errors for RAA vs. centrality 2011
2018 InitDplustoKpipi2011PbPb07half();
2eab2b1d 2019
2020 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.20);
2021 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2022 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.20);
2023
2024 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2025 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2026 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2027
2028 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2029 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2030 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
3db225bf 2031}
2eab2b1d 2032
3db225bf 2033//--------------------------------------------------------------------------
2034void AliHFSystErr::InitDplustoKpipi2011PbPb4050CentScan(){
2035 // define errors for RAA vs. centrality 2011
2036 InitDplustoKpipi2011PbPb07half();
2eab2b1d 2037
2038 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2039 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2040 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.05);
2041
2042 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2043 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2044 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2045
2046 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2047 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2048 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
3db225bf 2049}
2eab2b1d 2050
3db225bf 2051//--------------------------------------------------------------------------
2052void AliHFSystErr::InitDplustoKpipi2010PbPb5080CentScan(){
2053 // define errors for RAA vs. centrality 2011
2eab2b1d 2054 InitDplustoKpipi2010PbPb4080();
2055
2056 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2057 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2058 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.10);
2059
2060 for(Int_t i=1;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2061 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20); // 10%
2062 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10); // 10%
2063 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20); // 10%
3db225bf 2064}
2eab2b1d 2065
3db225bf 2066//--------------------------------------------------------------------------
2067void AliHFSystErr::InitDstartoD0pi2011PbPb010CentScan(){
2068 // define errors for RAA vs. centrality 2011
2069 InitDstartoD0pi2011PbPb07half();
2070
2071}
2072//--------------------------------------------------------------------------
2073void AliHFSystErr::InitDstartoD0pi2011PbPb1020CentScan(){
2074 // define errors for RAA vs. centrality 2011
2075 InitDstartoD0pi2011PbPb07half();
2076
2077}
2078//--------------------------------------------------------------------------
2079void AliHFSystErr::InitDstartoD0pi2011PbPb2030CentScan(){
2080 // define errors for RAA vs. centrality 2011
2081 InitDstartoD0pi2011PbPb07half();
2082
2083}
2084//--------------------------------------------------------------------------
2085void AliHFSystErr::InitDstartoD0pi2011PbPb3040CentScan(){
2086 // define errors for RAA vs. centrality 2011
2087 InitDstartoD0pi2011PbPb07half();
2088
2089}
2090//--------------------------------------------------------------------------
2091void AliHFSystErr::InitDstartoD0pi2011PbPb4050CentScan(){
2092 // define errors for RAA vs. centrality 2011
2093 InitDstartoD0pi2011PbPb07half();
2094
2095}
2096//--------------------------------------------------------------------------
2097void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
2098 // define errors for RAA vs. centrality 2011
2099 InitDstartoD0pi2011PbPb07half();
2100
2101}
2102
44dd56b2 2103//_________________________________________________________________________
2104void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
2105 InitD0toKpi2013pPb0100();
2106
a1f0af09 2107 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2108
2109 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2110 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.06);
44dd56b2 2111 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2112
44dd56b2 2113 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2114
a1f0af09 2115 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2116 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2117 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2118
2119 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2120 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2121 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2122
2123
2124};
2125//_________________________________________________________________________
2126void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
2127 InitD0toKpi2013pPb0100();
a1f0af09 2128
2129 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
44dd56b2 2130
a1f0af09 2131 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2132 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.08);
44dd56b2 2133 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2134
44dd56b2 2135 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2136
a1f0af09 2137 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2138 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2139 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2140
2141 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2142 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2143 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2144};
2145//_________________________________________________________________________
2146void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
a1f0af09 2147 InitD0toKpi2013pPb0100();
2148
2149 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
44dd56b2 2150
a1f0af09 2151 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2152 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.12);
44dd56b2 2153 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2154
44dd56b2 2155 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2156
a1f0af09 2157 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2158 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2159 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2160
2161 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2162 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2163 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2164
2165};
2166//_________________________________________________________________________
2167void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
2168 InitD0toKpi2013pPb0100();
2169
a1f0af09 2170 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
44dd56b2 2171
a1f0af09 2172 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2173 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2174 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2175
2176 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2177
2178 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2179 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2180 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2181
2182 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2183 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2184 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2185
2186};
2187//_________________________________________________________________________
2188void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
2189 InitD0toKpi2013pPb0100();
2190
a1f0af09 2191 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2192
2193 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2194 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2195 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
44dd56b2 2196
44dd56b2 2197 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2198
a1f0af09 2199 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2200 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2201 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2202
2203 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2204 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2205 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2206
2207};
2208
2209//_________________________________________________________________________
2210void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
2211 InitDplustoKpipi2013pPb0100();
2212
a1f0af09 2213 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
44dd56b2 2214
a1f0af09 2215 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2216 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2217 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2218
44dd56b2 2219 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2220
a1f0af09 2221 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2222 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2223 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2224
2225 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2226 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.10);
2227 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.10);
44dd56b2 2228
2229};
2230//_________________________________________________________________________
2231void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
2232 InitDplustoKpipi2013pPb0100();
a1f0af09 2233 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
44dd56b2 2234
a1f0af09 2235 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2236 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2237 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2238
44dd56b2 2239 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2240
a1f0af09 2241 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2242 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2243 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2244
2245 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2246 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2247 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2248
2249};
2250//_________________________________________________________________________
2251void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
2252 InitDplustoKpipi2013pPb0100();
a1f0af09 2253 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
44dd56b2 2254
a1f0af09 2255 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2256 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2257 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2258
44dd56b2 2259 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2260
a1f0af09 2261 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2262 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2263 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2264
2265 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2266 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2267 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2268
2269};
2270//_________________________________________________________________________
2271void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
2272 InitDplustoKpipi2013pPb0100();
a1f0af09 2273 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2274
2275 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2276 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2277 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
a1f0af09 2278
44dd56b2 2279 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2280
a1f0af09 2281 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2282 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2283 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
44dd56b2 2284
a1f0af09 2285 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2286 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2287 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2288
44dd56b2 2289};
2290//_________________________________________________________________________
2291void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
2292 InitDplustoKpipi2013pPb0100();
2293
a1f0af09 2294 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2295
2296 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2297 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2298 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
44dd56b2 2299
44dd56b2 2300 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2301
a1f0af09 2302 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2303 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2304 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
44dd56b2 2305
a1f0af09 2306 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2307 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2308 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2309
44dd56b2 2310};
2311
2312
2313//_________________________________________________________________________
2314void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
2315 InitDstartoD0pi2013pPb0100();
2316
a1f0af09 2317 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2318
2319 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2320 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
44dd56b2 2321 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2322
44dd56b2 2323 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2324
a1f0af09 2325 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2326 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2327 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2328
2329 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2330 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2331 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2332
2333};
2334//_________________________________________________________________________
2335void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
2336 InitDstartoD0pi2013pPb0100();
2337
a1f0af09 2338 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2339
2340 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2341 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
44dd56b2 2342 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2343
44dd56b2 2344 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2345
a1f0af09 2346 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2347 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2348 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2349
2350 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2351 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2352 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2353
2354};
2355//_________________________________________________________________________
2356void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
2357InitDstartoD0pi2013pPb0100();
2358
a1f0af09 2359 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2360
2361 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2362 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
44dd56b2 2363 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2364
44dd56b2 2365 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2366
a1f0af09 2367 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2368 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2369 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
44dd56b2 2370
a1f0af09 2371 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2372 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2373 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2374
44dd56b2 2375};
2376//_________________________________________________________________________
2377void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
2378 InitDstartoD0pi2013pPb0100();
a1f0af09 2379
2380 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
44dd56b2 2381
a1f0af09 2382 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2383 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
44dd56b2 2384 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2385
44dd56b2 2386 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2387
a1f0af09 2388 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2389 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2390 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2391
2392 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2393 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2394 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
44dd56b2 2395};
2396//_________________________________________________________________________
2397void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
2398 InitDstartoD0pi2013pPb0100();
a1f0af09 2399
2400 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
44dd56b2 2401
a1f0af09 2402 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2403 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.07);
44dd56b2 2404 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2405
44dd56b2 2406 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2407
a1f0af09 2408 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2409 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2410 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
44dd56b2 2411
a1f0af09 2412 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2413 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2414 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2415
44dd56b2 2416};
3db225bf 2417
fa90243b 2418//_________________________________________________________________________
2419void AliHFSystErr::InitD0toKpi2013pPb020V0A(){
2420 InitD0toKpi2013pPb0100();
2421
2422 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2423
2424
2425 // Raw yield extraction
2426 fRawYield->SetBinContent(1,0.0);
751d48b6 2427 fRawYield->SetBinContent(2,0.085);
fa90243b 2428 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2429 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 2430 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2431 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2432 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
fa90243b 2433 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2434
2435 // Cuts efficiency (from cuts variation)
2436 fCutsEff->SetBinContent(1,0.0);
2437 fCutsEff->SetBinContent(2,0.10);
751d48b6 2438 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
fa90243b 2439
2440 // PID efficiency (from PID/noPID)
2441 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2442
2443
2444}
2445
2446//_________________________________________________________________________
2447void AliHFSystErr::InitD0toKpi2013pPb2040V0A(){
2448 InitD0toKpi2013pPb0100();
2449
2450 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2451
2452
2453 // Raw yield extraction
2454 fRawYield->SetBinContent(1,0.0);
751d48b6 2455 fRawYield->SetBinContent(2,0.085);
2456 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
fa90243b 2457 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 2458 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2459 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2460 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2461 //for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2462
fa90243b 2463
2464 // Cuts efficiency (from cuts variation)
2465 fCutsEff->SetBinContent(1,0.0);
2466 fCutsEff->SetBinContent(2,0.10);
751d48b6 2467 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2468
fa90243b 2469
2470 // PID efficiency (from PID/noPID)
2471 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2472
2473}
2474
2475//_________________________________________________________________________
2476void AliHFSystErr::InitD0toKpi2013pPb4060V0A(){
2477 InitD0toKpi2013pPb0100();
2478
2479 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2480
2481
2482 // Raw yield extraction
2483 fRawYield->SetBinContent(1,0.0);
751d48b6 2484 fRawYield->SetBinContent(2,0.085);
fa90243b 2485 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 2486 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2487 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2488 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2489 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
fa90243b 2490 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2491
2492 // Cuts efficiency (from cuts variation)
2493 fCutsEff->SetBinContent(1,0.0);
2494 fCutsEff->SetBinContent(2,0.10);
751d48b6 2495 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2496
fa90243b 2497
2498 // PID efficiency (from PID/noPID)
2499 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2500
2501}
2502
2503//_________________________________________________________________________
2504void AliHFSystErr::InitD0toKpi2013pPb60100V0A(){
2505 InitD0toKpi2013pPb0100();
2506
2507 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2508
2509 // Raw yield extraction
2510 fRawYield->SetBinContent(1,0.0);
751d48b6 2511 fRawYield->SetBinContent(2,0.085);
fa90243b 2512 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 2513 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2514 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2515 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2516 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2517
fa90243b 2518
2519 // Cuts efficiency (from cuts variation)
2520 fCutsEff->SetBinContent(1,0.0);
751d48b6 2521 fCutsEff->SetBinContent(2,0.10);
2522 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2523
fa90243b 2524
2525 // PID efficiency (from PID/noPID)
751d48b6 2526 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 2527
2528}
2529
2530//_________________________________________________________________________
2531void AliHFSystErr::InitD0toKpi2013pPb020ZNA(){
2532 InitD0toKpi2013pPb0100();
2533
2534 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
2535
2536 // Raw yield extraction
2537 fRawYield->SetBinContent(1,0.0);
751d48b6 2538 fRawYield->SetBinContent(2,0.085);
2539 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2540 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2541 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2542 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2543 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2544 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
fa90243b 2545
2546 // Cuts efficiency (from cuts variation)
2547 fCutsEff->SetBinContent(1,0.0);
751d48b6 2548 fCutsEff->SetBinContent(2,0.10);
2549 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2550
fa90243b 2551
2552 // PID efficiency (from PID/noPID)
2553 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2554
2555}
2556
2557//_________________________________________________________________________
2558void AliHFSystErr::InitD0toKpi2013pPb2040ZNA(){
2559 InitD0toKpi2013pPb0100();
2560
2561 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
2562
2563 // Raw yield extraction
2564 fRawYield->SetBinContent(1,0.0);
751d48b6 2565 fRawYield->SetBinContent(2,0.085);
2566 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2567 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2568 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2569 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2570 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2571 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2572
fa90243b 2573
2574 // Cuts efficiency (from cuts variation)
2575 fCutsEff->SetBinContent(1,0.0);
2576 fCutsEff->SetBinContent(2,0.10);
751d48b6 2577 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2578
fa90243b 2579 // PID efficiency (from PID/noPID)
2580 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2581
2582}
2583
2584//_________________________________________________________________________
2585void AliHFSystErr::InitD0toKpi2013pPb4060ZNA(){
2586 InitD0toKpi2013pPb0100();
2587
2588 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
2589
2590 // Raw yield extraction
2591 fRawYield->SetBinContent(1,0.0);
751d48b6 2592 fRawYield->SetBinContent(2,0.085);
fa90243b 2593 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2594 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 2595 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2596 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2597 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2598 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
fa90243b 2599
2600 // Cuts efficiency (from cuts variation)
2601 fCutsEff->SetBinContent(1,0.0);
751d48b6 2602 fCutsEff->SetBinContent(2,0.10);
2603 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2604
fa90243b 2605 // PID efficiency (from PID/noPID)
2606 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2607
2608}
2609
2610//_________________________________________________________________________
2611void AliHFSystErr::InitD0toKpi2013pPb60100ZNA(){
2612 InitD0toKpi2013pPb0100();
2613
2614 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
2615
2616 // Raw yield extraction
2617 fRawYield->SetBinContent(1,0.0);
751d48b6 2618 fRawYield->SetBinContent(2,0.085);
fa90243b 2619 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 2620 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2621 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2622 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2623 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2624 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
fa90243b 2625
2626 // Cuts efficiency (from cuts variation)
2627 fCutsEff->SetBinContent(1,0.0);
751d48b6 2628 fCutsEff->SetBinContent(2,0.10);
2629 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2630
fa90243b 2631 // PID efficiency (from PID/noPID)
751d48b6 2632 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 2633
2634}
2635
2636//_________________________________________________________________________
2637void AliHFSystErr::InitDstartoD0pi2013pPb020V0A(){
2638 InitDstartoD0pi2013pPb0100();
2639
2640 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2641
2642 // Raw yield extraction
2643 fRawYield->SetBinContent(1,0.0);
2644 fRawYield->SetBinContent(2,0.10);
2645 fRawYield->SetBinContent(3,0.06);
2646 fRawYield->SetBinContent(4,0.06);
2647 fRawYield->SetBinContent(5,0.05);
2648 fRawYield->SetBinContent(6,0.05);
2649 fRawYield->SetBinContent(7,0.05);
2650 fRawYield->SetBinContent(8,0.05);
751d48b6 2651 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
fa90243b 2652 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2653 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2654
2655 // Cuts efficiency (from cuts variation)
2656 fCutsEff->SetBinContent(1,0.0);
2657 fCutsEff->SetBinContent(2,0.10);
2658 fCutsEff->SetBinContent(3,0.06);
2659 fCutsEff->SetBinContent(4,0.06);
2660 fCutsEff->SetBinContent(5,0.06);
2661 fCutsEff->SetBinContent(6,0.06);
2662 fCutsEff->SetBinContent(7,0.06);
2663 fCutsEff->SetBinContent(8,0.06);
751d48b6 2664 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2665 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2666
2667 // PID efficiency (from PID/noPID)
751d48b6 2668 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
fa90243b 2669
2670
2671}
2672
2673//_________________________________________________________________________
2674void AliHFSystErr::InitDstartoD0pi2013pPb2040V0A(){
2675 InitDstartoD0pi2013pPb0100();
2676
2677 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2678
2679 // Raw yield extraction
2680 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2681 fRawYield->SetBinContent(1,0.0);
2682 fRawYield->SetBinContent(2,0.10);
fa90243b 2683 fRawYield->SetBinContent(3,0.06);
2684 fRawYield->SetBinContent(4,0.06);
2685 fRawYield->SetBinContent(5,0.05);
2686 fRawYield->SetBinContent(6,0.05);
2687 fRawYield->SetBinContent(7,0.05);
2688 fRawYield->SetBinContent(8,0.05);
751d48b6 2689 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
fa90243b 2690 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
751d48b6 2691 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
fa90243b 2692
2693 // Cuts efficiency (from cuts variation)
2694 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2695 fCutsEff->SetBinContent(1,0.);
2696 fCutsEff->SetBinContent(2,0.10);
2697 fCutsEff->SetBinContent(3,0.06);
2698 fCutsEff->SetBinContent(4,0.06);
2699 fCutsEff->SetBinContent(5,0.06);
2700 fCutsEff->SetBinContent(6,0.06);
2701 fCutsEff->SetBinContent(7,0.06);
2702 fCutsEff->SetBinContent(8,0.06);
2703 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
751d48b6 2704 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2705 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2706
2707 // PID efficiency (from PID/noPID)
2708 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2709 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2710
2711 // MC dN/dpt
2712 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2713 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2714 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2715 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2716
2717 // particle-antiparticle
2718 /*
2719 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2720 fPartAntipart->SetBinContent(1,1);
2721 fPartAntipart->SetBinContent(2,1);
2722 fPartAntipart->SetBinContent(3,0.12);
2723 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2724 */
2725
2726}
2727
2728//_________________________________________________________________________
2729void AliHFSystErr::InitDstartoD0pi2013pPb4060V0A(){
2730 InitDstartoD0pi2013pPb0100();
2731
2732 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2733
2734 // Raw yield extraction
2735 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2736 fRawYield->SetBinContent(1,0.0);
2737 fRawYield->SetBinContent(2,0.0);
2738 fRawYield->SetBinContent(3,0.06);
2739 fRawYield->SetBinContent(4,0.06);
fa90243b 2740 fRawYield->SetBinContent(5,0.05);
2741 fRawYield->SetBinContent(6,0.05);
2742 fRawYield->SetBinContent(7,0.05);
2743 fRawYield->SetBinContent(8,0.05);
751d48b6 2744 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
fa90243b 2745 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
751d48b6 2746 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
fa90243b 2747
2748 // Cuts efficiency (from cuts variation)
2749 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2750 fCutsEff->SetBinContent(1,0.);
2751 fCutsEff->SetBinContent(2,0.);
2752 fCutsEff->SetBinContent(3,0.06);
2753 fCutsEff->SetBinContent(4,0.06);
2754 fCutsEff->SetBinContent(5,0.06);
2755 fCutsEff->SetBinContent(6,0.06);
751d48b6 2756 fCutsEff->SetBinContent(7,0.06);
2757 fCutsEff->SetBinContent(8,0.06);
2758 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2759 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2760 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2761
2762 // PID efficiency (from PID/noPID)
2763 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2764 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2765
2766 // MC dN/dpt
2767 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2768 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2769 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2770 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2771
2772 // particle-antiparticle
2773 /*
2774 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2775 fPartAntipart->SetBinContent(1,1);
2776 fPartAntipart->SetBinContent(2,1);
2777 fPartAntipart->SetBinContent(3,0.12);
2778 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2779 */
2780
2781
2782}
2783
2784//_________________________________________________________________________
2785void AliHFSystErr::InitDstartoD0pi2013pPb60100V0A(){
2786 InitDstartoD0pi2013pPb0100();
2787
2788 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2789
2790 // Raw yield extraction
2791 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2792 fRawYield->SetBinContent(1,0.0);
2793 fRawYield->SetBinContent(2,0.0);
fa90243b 2794 fRawYield->SetBinContent(3,0.06);
2795 fRawYield->SetBinContent(4,0.06);
751d48b6 2796 fRawYield->SetBinContent(5,0.05);
2797 fRawYield->SetBinContent(6,0.05);
2798 fRawYield->SetBinContent(7,0.05);
2799 fRawYield->SetBinContent(8,0.05);
2800 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2801 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
fa90243b 2802 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2803
2804 // Cuts efficiency (from cuts variation)
2805 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2806 fCutsEff->SetBinContent(1,0.);
2807 fCutsEff->SetBinContent(2,0.);
751d48b6 2808 fCutsEff->SetBinContent(3,0.06);
2809 fCutsEff->SetBinContent(4,0.06);
fa90243b 2810 fCutsEff->SetBinContent(5,0.06);
2811 fCutsEff->SetBinContent(6,0.06);
2812 fCutsEff->SetBinContent(7,0.06);
2813 fCutsEff->SetBinContent(8,0.06);
751d48b6 2814 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2815 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2816 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2817
2818 // PID efficiency (from PID/noPID)
2819 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 2820 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 2821
2822 // MC dN/dpt
2823 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2824 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2825 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2826 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2827
2828 // particle-antiparticle
2829 /*
2830 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2831 fPartAntipart->SetBinContent(1,1);
2832 fPartAntipart->SetBinContent(2,1);
2833 fPartAntipart->SetBinContent(3,0.12);
2834 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2835 */
2836
2837
2838}
2839
2840//_________________________________________________________________________
2841void AliHFSystErr::InitDstartoD0pi2013pPb020ZNA(){
2842 InitDstartoD0pi2013pPb0100();
2843
2844 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2845
2846 // Raw yield extraction
2847 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2848 fRawYield->SetBinContent(1,0.);
2849 fRawYield->SetBinContent(2,0.10);
2850 fRawYield->SetBinContent(3,0.06);
2851 fRawYield->SetBinContent(4,0.06);
751d48b6 2852 fRawYield->SetBinContent(5,0.05);
2853 fRawYield->SetBinContent(6,0.05);
2854 fRawYield->SetBinContent(7,0.05);
2855 fRawYield->SetBinContent(8,0.05);
2856 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
fa90243b 2857 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2858 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2859
2860 // Cuts efficiency (from cuts variation)
2861 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2862 fCutsEff->SetBinContent(1,0.);
2863 fCutsEff->SetBinContent(2,0.10);
2864 fCutsEff->SetBinContent(3,0.06);
2865 fCutsEff->SetBinContent(4,0.06);
2866 fCutsEff->SetBinContent(5,0.06);
2867 fCutsEff->SetBinContent(6,0.06);
751d48b6 2868 fCutsEff->SetBinContent(7,0.06);
2869 fCutsEff->SetBinContent(8,0.06);
2870 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2871 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2872 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2873
2874 // PID efficiency (from PID/noPID)
2875 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2876 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
2877
2878 // MC dN/dpt
2879 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2880 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2881 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2882 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2883
2884 // particle-antiparticle
2885 /*
2886 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2887 fPartAntipart->SetBinContent(1,1);
2888 fPartAntipart->SetBinContent(2,1);
2889 fPartAntipart->SetBinContent(3,0.12);
2890 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2891 */
2892
2893
2894}
2895
2896//_________________________________________________________________________
2897void AliHFSystErr::InitDstartoD0pi2013pPb2040ZNA(){
2898 InitDstartoD0pi2013pPb0100();
2899
2900 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2901
2902 // Raw yield extraction
2903 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2904 fRawYield->SetBinContent(1,0.0);
2905 fRawYield->SetBinContent(2,0.0);
2906 fRawYield->SetBinContent(3,0.06);
2907 fRawYield->SetBinContent(4,0.06);
fa90243b 2908 fRawYield->SetBinContent(5,0.05);
2909 fRawYield->SetBinContent(6,0.05);
2910 fRawYield->SetBinContent(7,0.05);
2911 fRawYield->SetBinContent(8,0.05);
2912 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2913 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2914 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2915
2916 // Cuts efficiency (from cuts variation)
2917 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2918 fCutsEff->SetBinContent(1,0.);
2919 fCutsEff->SetBinContent(2,0.);
2920 fCutsEff->SetBinContent(3,0.06);
2921 fCutsEff->SetBinContent(4,0.06);
2922 fCutsEff->SetBinContent(5,0.06);
2923 fCutsEff->SetBinContent(6,0.06);
2924 fCutsEff->SetBinContent(7,0.06);
2925 fCutsEff->SetBinContent(8,0.06);
2926 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
751d48b6 2927 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2928 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2929
2930 // PID efficiency (from PID/noPID)
2931 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 2932 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 2933
2934 // MC dN/dpt
2935 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2936 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2937 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2938 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2939
2940 // particle-antiparticle
2941 /*
2942 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2943 fPartAntipart->SetBinContent(1,1);
2944 fPartAntipart->SetBinContent(2,1);
2945 fPartAntipart->SetBinContent(3,0.12);
2946 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2947 */
2948
2949
2950}
2951
2952//_________________________________________________________________________
2953void AliHFSystErr::InitDstartoD0pi2013pPb4060ZNA(){
2954 InitDstartoD0pi2013pPb0100();
2955
2956 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2957
2958 // Raw yield extraction
2959 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2960 fRawYield->SetBinContent(1,0.0);
2961 fRawYield->SetBinContent(2,0.0);
2962 fRawYield->SetBinContent(3,0.06);
2963 fRawYield->SetBinContent(4,0.06);
fa90243b 2964 fRawYield->SetBinContent(5,0.05);
2965 fRawYield->SetBinContent(6,0.05);
2966 fRawYield->SetBinContent(7,0.05);
2967 fRawYield->SetBinContent(8,0.05);
2968 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2969 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2970 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2971
2972 // Cuts efficiency (from cuts variation)
2973 // fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2974 fCutsEff->SetBinContent(1,0.);
2975 fCutsEff->SetBinContent(2,0.);
751d48b6 2976 fCutsEff->SetBinContent(3,0.06);
2977 fCutsEff->SetBinContent(4,0.06);
2978 fCutsEff->SetBinContent(5,0.06);
2979 fCutsEff->SetBinContent(6,0.06);
2980 fCutsEff->SetBinContent(7,0.06);
2981 fCutsEff->SetBinContent(8,0.06);
2982 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2983 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2984 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2985
2986 // PID efficiency (from PID/noPID)
2987 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 2988 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 2989
2990 // MC dN/dpt
2991 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2992 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2993 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2994 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2995
2996 // particle-antiparticle
2997 /*
2998 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2999 fPartAntipart->SetBinContent(1,1);
3000 fPartAntipart->SetBinContent(2,1);
3001 fPartAntipart->SetBinContent(3,0.12);
3002 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3003 */
3004
3005
3006
3007}
3008
3009//_________________________________________________________________________
3010void AliHFSystErr::InitDstartoD0pi2013pPb60100ZNA(){
3011 InitDstartoD0pi2013pPb0100();
3012
3013 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3014
3015 // Raw yield extraction
3016 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 3017 fRawYield->SetBinContent(1,0.);
3018 fRawYield->SetBinContent(2,0.);
fa90243b 3019 fRawYield->SetBinContent(3,0.06);
3020 fRawYield->SetBinContent(4,0.06);
3021 fRawYield->SetBinContent(5,0.05);
3022 fRawYield->SetBinContent(6,0.05);
751d48b6 3023 fRawYield->SetBinContent(7,0.05);
3024 fRawYield->SetBinContent(8,0.05);
fa90243b 3025 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3026 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3027 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3028
3029 // Cuts efficiency (from cuts variation)
3030 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3031 fCutsEff->SetBinContent(1,0.);
3032 fCutsEff->SetBinContent(2,0.);
751d48b6 3033 fCutsEff->SetBinContent(3,0.06);
3034 fCutsEff->SetBinContent(4,0.06);
3035 fCutsEff->SetBinContent(5,0.06);
3036 fCutsEff->SetBinContent(6,0.06);
3037 fCutsEff->SetBinContent(7,0.06);
3038 fCutsEff->SetBinContent(8,0.06);
3039 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3040 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 3041 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3042
3043 // PID efficiency (from PID/noPID)
3044 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 3045 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 3046
3047 // MC dN/dpt
3048 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3049 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3050 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3051 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3052
3053 // particle-antiparticle
3054 /*
3055 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3056 fPartAntipart->SetBinContent(1,1);
3057 fPartAntipart->SetBinContent(2,1);
3058 fPartAntipart->SetBinContent(3,0.12);
3059 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3060 */
3061
3062}
3063
3064
3065
414bc6a5 3066
9ca30701 3067//--------------------------------------------------------------------------
3068void AliHFSystErr::InitLctopKpi2010pp() {
3069 //
3070 // Lc->pKpi syst errors. Responsible: R. Romita
3071 // 2010 pp sample
3072 //
3073
3074 // Normalization
3075 fNorm = new TH1F("fNorm","fNorm",4,2,6);
3076 for(Int_t i=1;i<=4;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
3077
3078 // Tracking efficiency
3079 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",4,2,6);
3080 for(Int_t i=1;i<=4;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
3081
3082 // Raw yield extraction
3083 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3084 fRawYield->SetBinContent(1,0.20);
3085 fRawYield->SetBinContent(2,0.15);
3086 fRawYield->SetBinContent(3,0.10);
3087 fRawYield->SetBinContent(4,0.15);
3088
3089 fCutsEff = new TH1F("fCutsEff","fCutsEff",4,2,6);
3090 fCutsEff->SetBinContent(1,0.50);
3091 fCutsEff->SetBinContent(2,0.07);
3092 fCutsEff->SetBinContent(3,0.30);
3093 fCutsEff->SetBinContent(4,0.40);
3094
3095 // PID efficiency (from PID/noPID)
3096 fPIDEff = new TH1F("fPIDEff","fPIDEff",4,2,4);
3097 fPIDEff->SetBinContent(1,0.13); // 15%
3098 fPIDEff->SetBinContent(2,0.30); // 15%
3099 fPIDEff->SetBinContent(3,0.5); // 15%
3100 fPIDEff->SetBinContent(4,0.30); // 15%
3101
3102 // MC dN/dpt
3103 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
3104 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
3105 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
3106 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
3107 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
3108 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
3109}
3110
53c91ed7 3111//--------------------------------------------------------------------------
3112void AliHFSystErr::InitLctopK0S2010pp() {
3113 //
3114 // Lc->pK0S syst errors. Responsible: A. De Caro
3115 // 2010 pp sample
3116 // Only tool preparation: no physical mean for these values!
3117 //
3118
3119 // Normalization
3120 fNorm = new TH1F("fNorm","fNorm",12,0.,12.);
3121 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05); // 5% error on sigmaV0and // only a try
3122
3123 // Tracking efficiency
3124 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0.,12.);
3125 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.1); // 10% (4% per track) // only a try
3126
3127 // Raw yield extraction
3128 fRawYield = new TH1F("fRawYield","fRawYield",12,0.,12.);
3129 for(Int_t i=1;i<=2;i++) fRawYield->SetBinContent(i,1.); // only a try
3130 for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.1); // only a try
3131
3132 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0.,12.);
3133 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.1); // only a try
3134
3135 // PID efficiency (from PID/noPID)
3136 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0.,12.);
3137 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // only a try
3138
3139 // MC dN/dpt
3140 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0.,12.);
3141 for(Int_t i=1; i<=2;i++) fMCPtShape->SetBinContent(i,1.); // only a try
3142 for(Int_t i=3; i<=12;i++) fMCPtShape->SetBinContent(i,0.1); // only a try
3143
3144}
3145
5541b811 3146//--------------------------------------------------------------------------
0233abe6 3147Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
3148 //
3149 // Get error
3150 //
3151
3152 Int_t bin=fCutsEff->FindBin(pt);
3153
3154 return fCutsEff->GetBinContent(bin);
3155}
3156//--------------------------------------------------------------------------
3157Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
3158 //
3159 // Get error
3160 //
3161
3162 Int_t bin=fMCPtShape->FindBin(pt);
3163
3164 return fMCPtShape->GetBinContent(bin);
3165}
3166//--------------------------------------------------------------------------
3167Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
3168 //
3169 // Get error
3170 //
3171
3172 Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3173
3174 return TMath::Sqrt(err);
3175}
3176//--------------------------------------------------------------------------
3177Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
3178 //
3179 // Get error
3180 //
3181
3182 Int_t bin=fPIDEff->FindBin(pt);
3183
3184 return fPIDEff->GetBinContent(bin);
3185}
3186//--------------------------------------------------------------------------
3187Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
3188 //
3189 // Get error
3190 //
3191
3192 Int_t bin=fTrackingEff->FindBin(pt);
3193
3194 return fTrackingEff->GetBinContent(bin);
3195}
3196//--------------------------------------------------------------------------
3197Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
3198 //
3199 // Get error
3200 //
3201
3202 Int_t bin=fRawYield->FindBin(pt);
3203
3204 return fRawYield->GetBinContent(bin);
3205}
3206//--------------------------------------------------------------------------
3207Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
3208 //
3209 // Get error
3210 //
3211
3212 Int_t bin=fPartAntipart->FindBin(pt);
3213
3214 return fPartAntipart->GetBinContent(bin);
3215}
3216//--------------------------------------------------------------------------
3217Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
3218 //
3219 // Get total syst error (except norm. error)
3220 //
3221
3222 Double_t err=0.;
3223
3224 if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
3225 if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
f9a85c82 3226 // if(fBR) err += GetBRErr()*GetBRErr();
0233abe6 3227 if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
3228 if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
3229 if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3230 if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
3231
3232 err += feeddownErr*feeddownErr;
3233
3234 return TMath::Sqrt(err);
3235}
3236//---------------------------------------------------------------------------
5ad8ad40 3237void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
0233abe6 3238 //
3239 // Draw errors
3240 //
3241 gStyle->SetOptStat(0);
3242
a07118d0 3243 TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
3244 cSystErr->Range(0.20,-0.5,18.4,0.34);
3245 cSystErr->SetRightMargin(0.318);
0233abe6 3246 cSystErr->SetFillColor(0);
3247
436b185f 3248 TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
3249 hFrame->SetAxisRange(1.,35.9,"X");
a07118d0 3250 hFrame->SetAxisRange(-0.5,0.5,"Y");
0233abe6 3251 hFrame->Draw();
3252
a07118d0 3253 TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
3254 leg->SetTextSize(0.03601695);
0233abe6 3255 leg->SetFillStyle(0);
3256 leg->SetBorderSize(0);
2ee3afe2 3257
436b185f 3258 TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
2ee3afe2 3259 Int_t nbins = fNorm->GetNbinsX();
3260 TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
436b185f 3261 for(Int_t i=1;i<=36;i++) {
2ee3afe2 3262 Double_t pt = hTotErr->GetBinCenter(i);
3263 Double_t ptwidth = hTotErr->GetBinWidth(i);
3264
3265 if(grErrFeeddown) {
3266 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
3267 Double_t toterryl=0., toterryh=0.;
3268 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
3269 grErrFeeddown->GetPoint(j,x,y);
3270 errxh = grErrFeeddown->GetErrorXhigh(j);
3271 errxl = grErrFeeddown->GetErrorXlow(j);
a07118d0 3272 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
2ee3afe2 3273 erryh = grErrFeeddown->GetErrorYhigh(j);
3274 erryl = grErrFeeddown->GetErrorYlow(j);
3275 }
3276 }
3277 if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
3278 else toterryl = GetTotalSystErr(pt);
3279 if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
3280 else toterryh = GetTotalSystErr(pt);
3281
3282 hTotErr->SetBinContent(i,toterryh);
3283 gTotErr->SetPoint(i,pt,0.);
3284 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
3285 }
3286 else {
3287 hTotErr->SetBinContent(i,GetTotalSystErr(pt));
3288 gTotErr->SetPoint(i,pt,0.);
3289 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
3290 }
3291
3292 }
f968ef30 3293 gTotErr->SetLineColor(kBlack);
3294 gTotErr->SetFillColor(kRed);
3295 gTotErr->SetFillStyle(3002);
2ee3afe2 3296 gTotErr->Draw("2");
3297 leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
3298// hTotErr->SetLineColor(1);
3299// hTotErr->SetLineWidth(3);
3300// hTotErr->Draw("same");
3301// leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
3302
0233abe6 3303
e373b9f0 3304 fNorm->SetFillColor(1);
3305 fNorm->SetFillStyle(3002);
3306 //fNorm->Draw("same");
3307 //TH1F *hNormRefl = ReflectHisto(fNorm);
3308 //hNormRefl->Draw("same");
27360a0e 3309 Double_t norm = fNorm->GetBinContent(1)*100;
3310 leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
e373b9f0 3311
5ad8ad40 3312 if(grErrFeeddown) {
3313 grErrFeeddown->SetFillColor(kTeal-8);
3314 grErrFeeddown->SetFillStyle(3001);
3315 grErrFeeddown->Draw("2");
3316 leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
3317 }
0233abe6 3318 if(fTrackingEff) {
5ad8ad40 3319 fTrackingEff->SetFillColor(4);
3320 fTrackingEff->SetFillStyle(3006);
0233abe6 3321 fTrackingEff->Draw("same");
f968ef30 3322 TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
3323 hTrackingEffRefl->Draw("same");
0233abe6 3324 leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
3325 }
3326 if(fBR) {
3327 fBR->SetFillColor(6);
3328 fBR->SetFillStyle(3005);
5ad8ad40 3329 //fBR->SetFillStyle(3020);
0233abe6 3330 fBR->Draw("same");
f968ef30 3331 TH1F *hBRRefl = ReflectHisto(fBR);
3332 hBRRefl->Draw("same");
0233abe6 3333 leg->AddEntry(fBR,"Branching ratio","f");
3334 }
3335 if(fRawYield) {
5ad8ad40 3336 Int_t ci; // for color index setting
3337 ci = TColor::GetColor("#00cc00");
3338 fRawYield->SetLineColor(ci);
3339 // fRawYield->SetLineColor(3);
0233abe6 3340 fRawYield->SetLineWidth(3);
3341 fRawYield->Draw("same");
f968ef30 3342 TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
3343 hRawYieldRefl->Draw("same");
844405bc 3344 leg->AddEntry(fRawYield,"Yield extraction","l");
0233abe6 3345 }
3346 if(fCutsEff) {
3347 fCutsEff->SetLineColor(4);
3348 fCutsEff->SetLineWidth(3);
3349 fCutsEff->Draw("same");
f968ef30 3350 TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
3351 hCutsEffRefl->Draw("same");
0233abe6 3352 leg->AddEntry(fCutsEff,"Cuts efficiency","l");
3353 }
3354 if(fPIDEff) {
3355 fPIDEff->SetLineColor(7);
3356 fPIDEff->SetLineWidth(3);
3357 fPIDEff->Draw("same");
f968ef30 3358 TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
3359 hPIDEffRefl->Draw("same");
0233abe6 3360 leg->AddEntry(fPIDEff,"PID efficiency","l");
3361 }
3362 if(fMCPtShape) {
5ad8ad40 3363 Int_t ci = TColor::GetColor("#9933ff");
3364 fMCPtShape->SetLineColor(ci);
3365 // fMCPtShape->SetLineColor(8);
0233abe6 3366 fMCPtShape->SetLineWidth(3);
3367 fMCPtShape->Draw("same");
f968ef30 3368 TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
3369 hMCPtShapeRefl->Draw("same");
0233abe6 3370 leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
3371 }
3372 if(fPartAntipart) {
5ad8ad40 3373 Int_t ci = TColor::GetColor("#ff6600");
3374 fPartAntipart->SetLineColor(ci);
3375 // fPartAntipart->SetLineColor(9);
0233abe6 3376 fPartAntipart->SetLineWidth(3);
3377 fPartAntipart->Draw("same");
f968ef30 3378 TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
3379 hPartAntipartRefl->Draw("same");
0233abe6 3380 leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
3381 }
0233abe6 3382
0233abe6 3383
3384 leg->Draw();
3385
5ad8ad40 3386 cSystErr->SaveAs("RelativeSystematics.eps");
3387
0233abe6 3388 return;
3389}
f968ef30 3390//-------------------------------------------------------------------------
3391TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
3392 //
3393 // Clones and reflects histogram
3394 //
3395 TH1F *hout=(TH1F*)hin->Clone("hout");
3396 hout->Scale(-1.);
3397
3398 return hout;
3399}
9a8b5080 3400
ada850b0 3401
3402