]> 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
Adding different ways to extract the systematic uncertainties from the envelopes...
[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);
2427 fRawYield->SetBinContent(2,0.07);
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);
2430 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.08);
2431 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04);
2432 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
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);
2438 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.10);
2439 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.10);
2440 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2441 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.08);
2442 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2443 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.10);
2444
2445 // PID efficiency (from PID/noPID)
2446 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2447
2448
2449}
2450
2451//_________________________________________________________________________
2452void AliHFSystErr::InitD0toKpi2013pPb2040V0A(){
2453 InitD0toKpi2013pPb0100();
2454
2455 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2456
2457
2458 // Raw yield extraction
2459 fRawYield->SetBinContent(1,0.0);
2460 fRawYield->SetBinContent(2,0.09);
2461 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.06);
2462 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2463 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.06);
2464 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2465 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2466 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2467
2468 // Cuts efficiency (from cuts variation)
2469 fCutsEff->SetBinContent(1,0.0);
2470 fCutsEff->SetBinContent(2,0.10);
2471 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.10);
2472 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.03);
2473 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2474 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.08);
2475 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
2476// for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2477
2478 // PID efficiency (from PID/noPID)
2479 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2480
2481}
2482
2483//_________________________________________________________________________
2484void AliHFSystErr::InitD0toKpi2013pPb4060V0A(){
2485 InitD0toKpi2013pPb0100();
2486
2487 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2488
2489
2490 // Raw yield extraction
2491 fRawYield->SetBinContent(1,0.0);
2492 fRawYield->SetBinContent(2,0.07);
2493 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2494 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.06);
2495 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.08);
2496 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2497 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.20);
2498 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2499
2500 // Cuts efficiency (from cuts variation)
2501 fCutsEff->SetBinContent(1,0.0);
2502 fCutsEff->SetBinContent(2,0.10);
2503 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.05);
2504 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.08);
2505 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.10);
2506 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.08);
2507 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.18);
2508 // for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2509
2510 // PID efficiency (from PID/noPID)
2511 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2512
2513}
2514
2515//_________________________________________________________________________
2516void AliHFSystErr::InitD0toKpi2013pPb60100V0A(){
2517 InitD0toKpi2013pPb0100();
2518
2519 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2520
2521 // Raw yield extraction
2522 fRawYield->SetBinContent(1,0.0);
2523 fRawYield->SetBinContent(2,0.10);
2524 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2525 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.08);
2526 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.08);
2527 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2528 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2529 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2530
2531 // Cuts efficiency (from cuts variation)
2532 fCutsEff->SetBinContent(1,0.0);
2533 fCutsEff->SetBinContent(2,0.15);
2534 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.05);
2535 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.10);
2536 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.10);
2537 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
2538 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2539 // for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2540
2541 // PID efficiency (from PID/noPID)
2542 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.08);
2543
2544}
2545
2546//_________________________________________________________________________
2547void AliHFSystErr::InitD0toKpi2013pPb020ZNA(){
2548 InitD0toKpi2013pPb0100();
2549
2550 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
2551
2552 // Raw yield extraction
2553 fRawYield->SetBinContent(1,0.0);
2554 fRawYield->SetBinContent(2,0.08);
2555 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.05);
2556 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.05);
2557 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2558 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.07);
2559 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.11);
2560 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.16);
2561
2562 // Cuts efficiency (from cuts variation)
2563 fCutsEff->SetBinContent(1,0.0);
2564 fCutsEff->SetBinContent(2,0.18);
2565 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.05);
2566 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.05);
2567 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2568 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
2569 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2570 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.18);
2571
2572 // PID efficiency (from PID/noPID)
2573 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2574
2575}
2576
2577//_________________________________________________________________________
2578void AliHFSystErr::InitD0toKpi2013pPb2040ZNA(){
2579 InitD0toKpi2013pPb0100();
2580
2581 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
2582
2583 // Raw yield extraction
2584 fRawYield->SetBinContent(1,0.0);
2585 fRawYield->SetBinContent(2,0.15);
2586 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.05);
2587 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.05);
2588 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2589 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2590 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.16);
2591 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.16);
2592
2593 // Cuts efficiency (from cuts variation)
2594 fCutsEff->SetBinContent(1,0.0);
2595 fCutsEff->SetBinContent(2,0.10);
2596 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.05);
2597 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.05);
2598 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2599 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
2600 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2601 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2602
2603 // PID efficiency (from PID/noPID)
2604 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2605
2606}
2607
2608//_________________________________________________________________________
2609void AliHFSystErr::InitD0toKpi2013pPb4060ZNA(){
2610 InitD0toKpi2013pPb0100();
2611
2612 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
2613
2614 // Raw yield extraction
2615 fRawYield->SetBinContent(1,0.0);
2616 fRawYield->SetBinContent(2,0.09);
2617 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2618 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2619 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.07);
2620 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.07);
2621 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.2);
2622 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2);
2623
2624 // Cuts efficiency (from cuts variation)
2625 fCutsEff->SetBinContent(1,0.0);
2626 fCutsEff->SetBinContent(2,0.20);
2627 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.1);
2628 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.1);
2629 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.1);
2630 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.1);
2631 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.2);
2632 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.1);
2633
2634 // PID efficiency (from PID/noPID)
2635 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2636
2637}
2638
2639//_________________________________________________________________________
2640void AliHFSystErr::InitD0toKpi2013pPb60100ZNA(){
2641 InitD0toKpi2013pPb0100();
2642
2643 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
2644
2645 // Raw yield extraction
2646 fRawYield->SetBinContent(1,0.0);
2647 fRawYield->SetBinContent(2,0.05);
2648 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2649 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.05);
2650 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2651 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2652 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.11);
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.18);
2658 for(Int_t i=3;i<=4;i++) fCutsEff->SetBinContent(i,0.1);
2659 for(Int_t i=5;i<=6;i++) fCutsEff->SetBinContent(i,0.1);
2660 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.1);
2661 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.1);
2662 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.1);
2663 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2664
2665 // PID efficiency (from PID/noPID)
2666 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.08);
2667
2668}
2669
2670//_________________________________________________________________________
2671void AliHFSystErr::InitDstartoD0pi2013pPb020V0A(){
2672 InitDstartoD0pi2013pPb0100();
2673
2674 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2675
2676 // Raw yield extraction
2677 fRawYield->SetBinContent(1,0.0);
2678 fRawYield->SetBinContent(2,0.10);
2679 fRawYield->SetBinContent(3,0.06);
2680 fRawYield->SetBinContent(4,0.06);
2681 fRawYield->SetBinContent(5,0.05);
2682 fRawYield->SetBinContent(6,0.05);
2683 fRawYield->SetBinContent(7,0.05);
2684 fRawYield->SetBinContent(8,0.05);
2685 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2686 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2687 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2688
2689 // Cuts efficiency (from cuts variation)
2690 fCutsEff->SetBinContent(1,0.0);
2691 fCutsEff->SetBinContent(2,0.10);
2692 fCutsEff->SetBinContent(3,0.06);
2693 fCutsEff->SetBinContent(4,0.06);
2694 fCutsEff->SetBinContent(5,0.06);
2695 fCutsEff->SetBinContent(6,0.06);
2696 fCutsEff->SetBinContent(7,0.06);
2697 fCutsEff->SetBinContent(8,0.06);
2698 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2699 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2700
2701 // PID efficiency (from PID/noPID)
2702 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
2703
2704
2705}
2706
2707//_________________________________________________________________________
2708void AliHFSystErr::InitDstartoD0pi2013pPb2040V0A(){
2709 InitDstartoD0pi2013pPb0100();
2710
2711 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2712
2713 // Raw yield extraction
2714 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2715 fRawYield->SetBinContent(1,1.0);
2716 fRawYield->SetBinContent(2,0.12);
2717 fRawYield->SetBinContent(3,0.06);
2718 fRawYield->SetBinContent(4,0.06);
2719 fRawYield->SetBinContent(5,0.05);
2720 fRawYield->SetBinContent(6,0.05);
2721 fRawYield->SetBinContent(7,0.05);
2722 fRawYield->SetBinContent(8,0.05);
2723 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2724 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2725 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2726
2727 // Cuts efficiency (from cuts variation)
2728 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2729 fCutsEff->SetBinContent(1,0.);
2730 fCutsEff->SetBinContent(2,0.10);
2731 fCutsEff->SetBinContent(3,0.06);
2732 fCutsEff->SetBinContent(4,0.06);
2733 fCutsEff->SetBinContent(5,0.06);
2734 fCutsEff->SetBinContent(6,0.06);
2735 fCutsEff->SetBinContent(7,0.06);
2736 fCutsEff->SetBinContent(8,0.06);
2737 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2738 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2739 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2740
2741 // PID efficiency (from PID/noPID)
2742 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2743 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2744
2745 // MC dN/dpt
2746 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2747 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2748 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2749 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2750
2751 // particle-antiparticle
2752 /*
2753 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2754 fPartAntipart->SetBinContent(1,1);
2755 fPartAntipart->SetBinContent(2,1);
2756 fPartAntipart->SetBinContent(3,0.12);
2757 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2758 */
2759
2760}
2761
2762//_________________________________________________________________________
2763void AliHFSystErr::InitDstartoD0pi2013pPb4060V0A(){
2764 InitDstartoD0pi2013pPb0100();
2765
2766 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2767
2768 // Raw yield extraction
2769 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2770 fRawYield->SetBinContent(1,1.0);
2771 fRawYield->SetBinContent(2,1.0);
2772 fRawYield->SetBinContent(3,0.08);
2773 fRawYield->SetBinContent(4,0.08);
2774 fRawYield->SetBinContent(5,0.05);
2775 fRawYield->SetBinContent(6,0.05);
2776 fRawYield->SetBinContent(7,0.05);
2777 fRawYield->SetBinContent(8,0.05);
2778 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2779 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2780 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2781
2782 // Cuts efficiency (from cuts variation)
2783 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2784 fCutsEff->SetBinContent(1,0.);
2785 fCutsEff->SetBinContent(2,0.);
2786 fCutsEff->SetBinContent(3,0.06);
2787 fCutsEff->SetBinContent(4,0.06);
2788 fCutsEff->SetBinContent(5,0.06);
2789 fCutsEff->SetBinContent(6,0.06);
2790 fCutsEff->SetBinContent(7,0.10);
2791 fCutsEff->SetBinContent(8,0.10);
2792 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.10);
2793 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2794 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2795
2796 // PID efficiency (from PID/noPID)
2797 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2798 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2799
2800 // MC dN/dpt
2801 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2802 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2803 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2804 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2805
2806 // particle-antiparticle
2807 /*
2808 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2809 fPartAntipart->SetBinContent(1,1);
2810 fPartAntipart->SetBinContent(2,1);
2811 fPartAntipart->SetBinContent(3,0.12);
2812 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2813 */
2814
2815
2816}
2817
2818//_________________________________________________________________________
2819void AliHFSystErr::InitDstartoD0pi2013pPb60100V0A(){
2820 InitDstartoD0pi2013pPb0100();
2821
2822 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2823
2824 // Raw yield extraction
2825 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2826 fRawYield->SetBinContent(1,1.0);
2827 fRawYield->SetBinContent(2,1.0);
2828 fRawYield->SetBinContent(3,0.06);
2829 fRawYield->SetBinContent(4,0.06);
2830 fRawYield->SetBinContent(5,0.06);
2831 fRawYield->SetBinContent(6,0.06);
2832 fRawYield->SetBinContent(7,0.06);
2833 fRawYield->SetBinContent(8,0.06);
2834 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2835 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.12);
2836 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2837
2838 // Cuts efficiency (from cuts variation)
2839 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2840 fCutsEff->SetBinContent(1,0.);
2841 fCutsEff->SetBinContent(2,0.);
2842 fCutsEff->SetBinContent(3,0.08);
2843 fCutsEff->SetBinContent(4,0.08);
2844 fCutsEff->SetBinContent(5,0.06);
2845 fCutsEff->SetBinContent(6,0.06);
2846 fCutsEff->SetBinContent(7,0.06);
2847 fCutsEff->SetBinContent(8,0.06);
2848 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.10);
2849 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2850 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2851
2852 // PID efficiency (from PID/noPID)
2853 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2854 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.08);
2855
2856 // MC dN/dpt
2857 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2858 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2859 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2860 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2861
2862 // particle-antiparticle
2863 /*
2864 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2865 fPartAntipart->SetBinContent(1,1);
2866 fPartAntipart->SetBinContent(2,1);
2867 fPartAntipart->SetBinContent(3,0.12);
2868 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2869 */
2870
2871
2872}
2873
2874//_________________________________________________________________________
2875void AliHFSystErr::InitDstartoD0pi2013pPb020ZNA(){
2876 InitDstartoD0pi2013pPb0100();
2877
2878 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2879
2880 // Raw yield extraction
2881 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2882 fRawYield->SetBinContent(1,0.);
2883 fRawYield->SetBinContent(2,0.10);
2884 fRawYield->SetBinContent(3,0.06);
2885 fRawYield->SetBinContent(4,0.06);
2886 fRawYield->SetBinContent(5,0.04);
2887 fRawYield->SetBinContent(6,0.04);
2888 fRawYield->SetBinContent(7,0.04);
2889 fRawYield->SetBinContent(8,0.04);
2890 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2891 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2892 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2893
2894 // Cuts efficiency (from cuts variation)
2895 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2896 fCutsEff->SetBinContent(1,0.);
2897 fCutsEff->SetBinContent(2,0.10);
2898 fCutsEff->SetBinContent(3,0.06);
2899 fCutsEff->SetBinContent(4,0.06);
2900 fCutsEff->SetBinContent(5,0.06);
2901 fCutsEff->SetBinContent(6,0.06);
2902 fCutsEff->SetBinContent(7,0.08);
2903 fCutsEff->SetBinContent(8,0.08);
2904 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.08);
2905 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2906 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2907
2908 // PID efficiency (from PID/noPID)
2909 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2910 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
2911
2912 // MC dN/dpt
2913 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2914 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2915 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2916 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2917
2918 // particle-antiparticle
2919 /*
2920 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2921 fPartAntipart->SetBinContent(1,1);
2922 fPartAntipart->SetBinContent(2,1);
2923 fPartAntipart->SetBinContent(3,0.12);
2924 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2925 */
2926
2927
2928}
2929
2930//_________________________________________________________________________
2931void AliHFSystErr::InitDstartoD0pi2013pPb2040ZNA(){
2932 InitDstartoD0pi2013pPb0100();
2933
2934 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2935
2936 // Raw yield extraction
2937 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2938 fRawYield->SetBinContent(1,1.0);
2939 fRawYield->SetBinContent(2,1.0);
2940 fRawYield->SetBinContent(3,0.08);
2941 fRawYield->SetBinContent(4,0.08);
2942 fRawYield->SetBinContent(5,0.05);
2943 fRawYield->SetBinContent(6,0.05);
2944 fRawYield->SetBinContent(7,0.05);
2945 fRawYield->SetBinContent(8,0.05);
2946 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2947 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2948 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2949
2950 // Cuts efficiency (from cuts variation)
2951 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2952 fCutsEff->SetBinContent(1,0.);
2953 fCutsEff->SetBinContent(2,0.);
2954 fCutsEff->SetBinContent(3,0.06);
2955 fCutsEff->SetBinContent(4,0.06);
2956 fCutsEff->SetBinContent(5,0.06);
2957 fCutsEff->SetBinContent(6,0.06);
2958 fCutsEff->SetBinContent(7,0.06);
2959 fCutsEff->SetBinContent(8,0.06);
2960 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2961 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
2962 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2963
2964 // PID efficiency (from PID/noPID)
2965 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2966 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2967
2968 // MC dN/dpt
2969 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2970 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2971 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2972 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2973
2974 // particle-antiparticle
2975 /*
2976 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2977 fPartAntipart->SetBinContent(1,1);
2978 fPartAntipart->SetBinContent(2,1);
2979 fPartAntipart->SetBinContent(3,0.12);
2980 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2981 */
2982
2983
2984}
2985
2986//_________________________________________________________________________
2987void AliHFSystErr::InitDstartoD0pi2013pPb4060ZNA(){
2988 InitDstartoD0pi2013pPb0100();
2989
2990 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2991
2992 // Raw yield extraction
2993 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2994 fRawYield->SetBinContent(1,1.0);
2995 fRawYield->SetBinContent(2,1.0);
2996 fRawYield->SetBinContent(3,0.09);
2997 fRawYield->SetBinContent(4,0.09);
2998 fRawYield->SetBinContent(5,0.05);
2999 fRawYield->SetBinContent(6,0.05);
3000 fRawYield->SetBinContent(7,0.05);
3001 fRawYield->SetBinContent(8,0.05);
3002 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3003 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3004 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3005
3006 // Cuts efficiency (from cuts variation)
3007 // fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3008 fCutsEff->SetBinContent(1,0.);
3009 fCutsEff->SetBinContent(2,0.);
3010 fCutsEff->SetBinContent(3,0.10);
3011 fCutsEff->SetBinContent(4,0.10);
3012 fCutsEff->SetBinContent(5,0.08);
3013 fCutsEff->SetBinContent(6,0.08);
3014 fCutsEff->SetBinContent(7,0.08);
3015 fCutsEff->SetBinContent(8,0.08);
3016 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.08);
3017 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.12);
3018 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3019
3020 // PID efficiency (from PID/noPID)
3021 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3022 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
3023
3024 // MC dN/dpt
3025 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3026 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3027 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3028 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3029
3030 // particle-antiparticle
3031 /*
3032 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3033 fPartAntipart->SetBinContent(1,1);
3034 fPartAntipart->SetBinContent(2,1);
3035 fPartAntipart->SetBinContent(3,0.12);
3036 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3037 */
3038
3039
3040
3041}
3042
3043//_________________________________________________________________________
3044void AliHFSystErr::InitDstartoD0pi2013pPb60100ZNA(){
3045 InitDstartoD0pi2013pPb0100();
3046
3047 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3048
3049 // Raw yield extraction
3050 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3051 fRawYield->SetBinContent(1,1.0);
3052 fRawYield->SetBinContent(2,1.0);
3053 fRawYield->SetBinContent(3,0.06);
3054 fRawYield->SetBinContent(4,0.06);
3055 fRawYield->SetBinContent(5,0.05);
3056 fRawYield->SetBinContent(6,0.05);
3057 fRawYield->SetBinContent(7,0.08);
3058 fRawYield->SetBinContent(8,0.08);
3059 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3060 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3061 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3062
3063 // Cuts efficiency (from cuts variation)
3064 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3065 fCutsEff->SetBinContent(1,0.);
3066 fCutsEff->SetBinContent(2,0.);
3067 fCutsEff->SetBinContent(3,0.10);
3068 fCutsEff->SetBinContent(4,0.10);
3069 fCutsEff->SetBinContent(5,0.10);
3070 fCutsEff->SetBinContent(6,0.10);
3071 fCutsEff->SetBinContent(7,0.08);
3072 fCutsEff->SetBinContent(8,0.08);
3073 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.10);
3074 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
3075 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3076
3077 // PID efficiency (from PID/noPID)
3078 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3079 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.08);
3080
3081 // MC dN/dpt
3082 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3083 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3084 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3085 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3086
3087 // particle-antiparticle
3088 /*
3089 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3090 fPartAntipart->SetBinContent(1,1);
3091 fPartAntipart->SetBinContent(2,1);
3092 fPartAntipart->SetBinContent(3,0.12);
3093 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3094 */
3095
3096}
3097
3098
3099
414bc6a5 3100
9ca30701 3101//--------------------------------------------------------------------------
3102void AliHFSystErr::InitLctopKpi2010pp() {
3103 //
3104 // Lc->pKpi syst errors. Responsible: R. Romita
3105 // 2010 pp sample
3106 //
3107
3108 // Normalization
3109 fNorm = new TH1F("fNorm","fNorm",4,2,6);
3110 for(Int_t i=1;i<=4;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
3111
3112 // Tracking efficiency
3113 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",4,2,6);
3114 for(Int_t i=1;i<=4;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
3115
3116 // Raw yield extraction
3117 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3118 fRawYield->SetBinContent(1,0.20);
3119 fRawYield->SetBinContent(2,0.15);
3120 fRawYield->SetBinContent(3,0.10);
3121 fRawYield->SetBinContent(4,0.15);
3122
3123 fCutsEff = new TH1F("fCutsEff","fCutsEff",4,2,6);
3124 fCutsEff->SetBinContent(1,0.50);
3125 fCutsEff->SetBinContent(2,0.07);
3126 fCutsEff->SetBinContent(3,0.30);
3127 fCutsEff->SetBinContent(4,0.40);
3128
3129 // PID efficiency (from PID/noPID)
3130 fPIDEff = new TH1F("fPIDEff","fPIDEff",4,2,4);
3131 fPIDEff->SetBinContent(1,0.13); // 15%
3132 fPIDEff->SetBinContent(2,0.30); // 15%
3133 fPIDEff->SetBinContent(3,0.5); // 15%
3134 fPIDEff->SetBinContent(4,0.30); // 15%
3135
3136 // MC dN/dpt
3137 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
3138 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
3139 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
3140 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
3141 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
3142 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
3143}
3144
53c91ed7 3145//--------------------------------------------------------------------------
3146void AliHFSystErr::InitLctopK0S2010pp() {
3147 //
3148 // Lc->pK0S syst errors. Responsible: A. De Caro
3149 // 2010 pp sample
3150 // Only tool preparation: no physical mean for these values!
3151 //
3152
3153 // Normalization
3154 fNorm = new TH1F("fNorm","fNorm",12,0.,12.);
3155 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05); // 5% error on sigmaV0and // only a try
3156
3157 // Tracking efficiency
3158 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0.,12.);
3159 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.1); // 10% (4% per track) // only a try
3160
3161 // Raw yield extraction
3162 fRawYield = new TH1F("fRawYield","fRawYield",12,0.,12.);
3163 for(Int_t i=1;i<=2;i++) fRawYield->SetBinContent(i,1.); // only a try
3164 for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.1); // only a try
3165
3166 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0.,12.);
3167 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.1); // only a try
3168
3169 // PID efficiency (from PID/noPID)
3170 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0.,12.);
3171 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // only a try
3172
3173 // MC dN/dpt
3174 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0.,12.);
3175 for(Int_t i=1; i<=2;i++) fMCPtShape->SetBinContent(i,1.); // only a try
3176 for(Int_t i=3; i<=12;i++) fMCPtShape->SetBinContent(i,0.1); // only a try
3177
3178}
3179
5541b811 3180//--------------------------------------------------------------------------
0233abe6 3181Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
3182 //
3183 // Get error
3184 //
3185
3186 Int_t bin=fCutsEff->FindBin(pt);
3187
3188 return fCutsEff->GetBinContent(bin);
3189}
3190//--------------------------------------------------------------------------
3191Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
3192 //
3193 // Get error
3194 //
3195
3196 Int_t bin=fMCPtShape->FindBin(pt);
3197
3198 return fMCPtShape->GetBinContent(bin);
3199}
3200//--------------------------------------------------------------------------
3201Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
3202 //
3203 // Get error
3204 //
3205
3206 Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3207
3208 return TMath::Sqrt(err);
3209}
3210//--------------------------------------------------------------------------
3211Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
3212 //
3213 // Get error
3214 //
3215
3216 Int_t bin=fPIDEff->FindBin(pt);
3217
3218 return fPIDEff->GetBinContent(bin);
3219}
3220//--------------------------------------------------------------------------
3221Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
3222 //
3223 // Get error
3224 //
3225
3226 Int_t bin=fTrackingEff->FindBin(pt);
3227
3228 return fTrackingEff->GetBinContent(bin);
3229}
3230//--------------------------------------------------------------------------
3231Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
3232 //
3233 // Get error
3234 //
3235
3236 Int_t bin=fRawYield->FindBin(pt);
3237
3238 return fRawYield->GetBinContent(bin);
3239}
3240//--------------------------------------------------------------------------
3241Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
3242 //
3243 // Get error
3244 //
3245
3246 Int_t bin=fPartAntipart->FindBin(pt);
3247
3248 return fPartAntipart->GetBinContent(bin);
3249}
3250//--------------------------------------------------------------------------
3251Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
3252 //
3253 // Get total syst error (except norm. error)
3254 //
3255
3256 Double_t err=0.;
3257
3258 if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
3259 if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
f9a85c82 3260 // if(fBR) err += GetBRErr()*GetBRErr();
0233abe6 3261 if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
3262 if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
3263 if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3264 if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
3265
3266 err += feeddownErr*feeddownErr;
3267
3268 return TMath::Sqrt(err);
3269}
3270//---------------------------------------------------------------------------
5ad8ad40 3271void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
0233abe6 3272 //
3273 // Draw errors
3274 //
3275 gStyle->SetOptStat(0);
3276
a07118d0 3277 TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
3278 cSystErr->Range(0.20,-0.5,18.4,0.34);
3279 cSystErr->SetRightMargin(0.318);
0233abe6 3280 cSystErr->SetFillColor(0);
3281
436b185f 3282 TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
3283 hFrame->SetAxisRange(1.,35.9,"X");
a07118d0 3284 hFrame->SetAxisRange(-0.5,0.5,"Y");
0233abe6 3285 hFrame->Draw();
3286
a07118d0 3287 TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
3288 leg->SetTextSize(0.03601695);
0233abe6 3289 leg->SetFillStyle(0);
3290 leg->SetBorderSize(0);
2ee3afe2 3291
436b185f 3292 TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
2ee3afe2 3293 Int_t nbins = fNorm->GetNbinsX();
3294 TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
436b185f 3295 for(Int_t i=1;i<=36;i++) {
2ee3afe2 3296 Double_t pt = hTotErr->GetBinCenter(i);
3297 Double_t ptwidth = hTotErr->GetBinWidth(i);
3298
3299 if(grErrFeeddown) {
3300 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
3301 Double_t toterryl=0., toterryh=0.;
3302 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
3303 grErrFeeddown->GetPoint(j,x,y);
3304 errxh = grErrFeeddown->GetErrorXhigh(j);
3305 errxl = grErrFeeddown->GetErrorXlow(j);
a07118d0 3306 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
2ee3afe2 3307 erryh = grErrFeeddown->GetErrorYhigh(j);
3308 erryl = grErrFeeddown->GetErrorYlow(j);
3309 }
3310 }
3311 if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
3312 else toterryl = GetTotalSystErr(pt);
3313 if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
3314 else toterryh = GetTotalSystErr(pt);
3315
3316 hTotErr->SetBinContent(i,toterryh);
3317 gTotErr->SetPoint(i,pt,0.);
3318 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
3319 }
3320 else {
3321 hTotErr->SetBinContent(i,GetTotalSystErr(pt));
3322 gTotErr->SetPoint(i,pt,0.);
3323 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
3324 }
3325
3326 }
f968ef30 3327 gTotErr->SetLineColor(kBlack);
3328 gTotErr->SetFillColor(kRed);
3329 gTotErr->SetFillStyle(3002);
2ee3afe2 3330 gTotErr->Draw("2");
3331 leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
3332// hTotErr->SetLineColor(1);
3333// hTotErr->SetLineWidth(3);
3334// hTotErr->Draw("same");
3335// leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
3336
0233abe6 3337
e373b9f0 3338 fNorm->SetFillColor(1);
3339 fNorm->SetFillStyle(3002);
3340 //fNorm->Draw("same");
3341 //TH1F *hNormRefl = ReflectHisto(fNorm);
3342 //hNormRefl->Draw("same");
27360a0e 3343 Double_t norm = fNorm->GetBinContent(1)*100;
3344 leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
e373b9f0 3345
5ad8ad40 3346 if(grErrFeeddown) {
3347 grErrFeeddown->SetFillColor(kTeal-8);
3348 grErrFeeddown->SetFillStyle(3001);
3349 grErrFeeddown->Draw("2");
3350 leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
3351 }
0233abe6 3352 if(fTrackingEff) {
5ad8ad40 3353 fTrackingEff->SetFillColor(4);
3354 fTrackingEff->SetFillStyle(3006);
0233abe6 3355 fTrackingEff->Draw("same");
f968ef30 3356 TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
3357 hTrackingEffRefl->Draw("same");
0233abe6 3358 leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
3359 }
3360 if(fBR) {
3361 fBR->SetFillColor(6);
3362 fBR->SetFillStyle(3005);
5ad8ad40 3363 //fBR->SetFillStyle(3020);
0233abe6 3364 fBR->Draw("same");
f968ef30 3365 TH1F *hBRRefl = ReflectHisto(fBR);
3366 hBRRefl->Draw("same");
0233abe6 3367 leg->AddEntry(fBR,"Branching ratio","f");
3368 }
3369 if(fRawYield) {
5ad8ad40 3370 Int_t ci; // for color index setting
3371 ci = TColor::GetColor("#00cc00");
3372 fRawYield->SetLineColor(ci);
3373 // fRawYield->SetLineColor(3);
0233abe6 3374 fRawYield->SetLineWidth(3);
3375 fRawYield->Draw("same");
f968ef30 3376 TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
3377 hRawYieldRefl->Draw("same");
844405bc 3378 leg->AddEntry(fRawYield,"Yield extraction","l");
0233abe6 3379 }
3380 if(fCutsEff) {
3381 fCutsEff->SetLineColor(4);
3382 fCutsEff->SetLineWidth(3);
3383 fCutsEff->Draw("same");
f968ef30 3384 TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
3385 hCutsEffRefl->Draw("same");
0233abe6 3386 leg->AddEntry(fCutsEff,"Cuts efficiency","l");
3387 }
3388 if(fPIDEff) {
3389 fPIDEff->SetLineColor(7);
3390 fPIDEff->SetLineWidth(3);
3391 fPIDEff->Draw("same");
f968ef30 3392 TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
3393 hPIDEffRefl->Draw("same");
0233abe6 3394 leg->AddEntry(fPIDEff,"PID efficiency","l");
3395 }
3396 if(fMCPtShape) {
5ad8ad40 3397 Int_t ci = TColor::GetColor("#9933ff");
3398 fMCPtShape->SetLineColor(ci);
3399 // fMCPtShape->SetLineColor(8);
0233abe6 3400 fMCPtShape->SetLineWidth(3);
3401 fMCPtShape->Draw("same");
f968ef30 3402 TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
3403 hMCPtShapeRefl->Draw("same");
0233abe6 3404 leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
3405 }
3406 if(fPartAntipart) {
5ad8ad40 3407 Int_t ci = TColor::GetColor("#ff6600");
3408 fPartAntipart->SetLineColor(ci);
3409 // fPartAntipart->SetLineColor(9);
0233abe6 3410 fPartAntipart->SetLineWidth(3);
3411 fPartAntipart->Draw("same");
f968ef30 3412 TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
3413 hPartAntipartRefl->Draw("same");
0233abe6 3414 leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
3415 }
0233abe6 3416
0233abe6 3417
3418 leg->Draw();
3419
5ad8ad40 3420 cSystErr->SaveAs("RelativeSystematics.eps");
3421
0233abe6 3422 return;
3423}
f968ef30 3424//-------------------------------------------------------------------------
3425TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
3426 //
3427 // Clones and reflects histogram
3428 //
3429 TH1F *hout=(TH1F*)hin->Clone("hout");
3430 hout->Scale(-1.);
3431
3432 return hout;
3433}
9a8b5080 3434
ada850b0 3435
3436