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