2 /**************************************************************************
3 * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
19 /////////////////////////////////////////////////////////////
21 // Class to handle systematic errors for charm hadrons
24 // AliHFSystEff syst; // DECAY = 1 for D0, 2, for D+, 3 for D*, 5 for Lc->pKp, 6 for Lc->pK0S
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*
29 // syst.DrawErrors(); // to see a plot of the error contributions
30 // syst.GetTotalSystErr(pt); // to get the total err at pt
32 // Author: A.Dainese, andrea.dainese@pd.infn.it
33 /////////////////////////////////////////////////////////////
36 #include <TGraphAsymmErrors.h>
44 #include "AliHFSystErr.h"
47 ClassImp(AliHFSystErr)
49 //--------------------------------------------------------------------------
50 AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
62 fCentralityClass("0100"),
63 fRapidityRange("0101"),
66 fIsRapidityScan(false)
69 // Default Constructor
73 //--------------------------------------------------------------------------
74 AliHFSystErr::~AliHFSystErr() {
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; }
91 //--------------------------------------------------------------------------
92 void AliHFSystErr::Init(Int_t decay){
94 // Variables/histos initialization
97 if ((fRunNumber>11) && fIsLowEnergy==false) {
98 AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
103 if (fCollisionType==0) {
104 if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
105 else InitD0toKpi2010pp();
107 else if (fCollisionType==1) {
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");
118 else if (fRunNumber == 11){
119 if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
120 else if (fCentralityClass=="010") InitD0toKpi2011PbPb010();
121 else if (fCentralityClass=="3050InPlane") InitD0toKpi2011PbPb3050InPlane();
122 else if (fCentralityClass=="3050OutOfPlane")InitD0toKpi2011PbPb3050OutOfPlane();
123 else if (fCentralityClass == "3050")InitD0toKpi2011PbPb3050();
124 else if (fCentralityClass=="010" && fIsCentScan) InitD0toKpi2011PbPb010CentScan();
125 else if (fCentralityClass=="1020") InitD0toKpi2011PbPb1020CentScan();
126 else if (fCentralityClass=="2030") InitD0toKpi2011PbPb2030CentScan();
127 else if (fCentralityClass=="3040") InitD0toKpi2011PbPb3040CentScan();
128 else if (fCentralityClass=="4050") InitD0toKpi2011PbPb4050CentScan();
129 else if (fCentralityClass=="5080") InitD0toKpi2010PbPb5080CentScan();
130 else AliFatal("Not yet implemented");
133 else if (fCollisionType==2) {
134 if (fCentralityClass=="0100") InitD0toKpi2013pPb0100();
136 if (fCentralityClass=="020V0A") InitD0toKpi2013pPb020V0A();
137 if (fCentralityClass=="2040V0A") InitD0toKpi2013pPb2040V0A();
138 if (fCentralityClass=="4060V0A") InitD0toKpi2013pPb4060V0A();
139 if (fCentralityClass=="60100V0A") InitD0toKpi2013pPb60100V0A();
141 if (fCentralityClass=="020ZNA") InitD0toKpi2013pPb020ZNA();
142 if (fCentralityClass=="2040ZNA") InitD0toKpi2013pPb2040ZNA();
143 if (fCentralityClass=="4060ZNA") InitD0toKpi2013pPb4060ZNA();
144 if (fCentralityClass=="60100ZNA") InitD0toKpi2013pPb60100ZNA();
146 if (fIsRapidityScan) {
147 if (fRapidityRange == "0804") InitD0toKpi2013pPb0100RapScan0804();
148 if (fRapidityRange == "0401") InitD0toKpi2013pPb0100RapScan0401();
149 if (fRapidityRange == "0101") InitD0toKpi2013pPb0100RapScan0101();
150 if (fRapidityRange == "0104") InitD0toKpi2013pPb0100RapScan0104();
151 if (fRapidityRange == "0408") InitD0toKpi2013pPb0100RapScan0408();
154 else AliFatal("Not yet implemented");
158 if (fCollisionType==0) {
159 if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
160 else InitDplustoKpipi2010pp();
162 else if (fCollisionType==1) {
163 if (fRunNumber == 10){
164 if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
165 else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
166 else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
167 else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
168 else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
169 else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
170 else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
171 else AliFatal("Not yet implemented");
173 if(fRunNumber == 11){
174 if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
175 else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010();
176 else if (fCentralityClass=="010" && fIsCentScan) InitDplustoKpipi2011PbPb010CentScan();
177 else if (fCentralityClass=="1020") InitDplustoKpipi2011PbPb1020CentScan();
178 else if (fCentralityClass=="2030") InitDplustoKpipi2011PbPb2030CentScan();
179 else if (fCentralityClass=="3040") InitDplustoKpipi2011PbPb3040CentScan();
180 else if (fCentralityClass=="4050") InitDplustoKpipi2011PbPb4050CentScan();
181 else if (fCentralityClass=="5080") InitDplustoKpipi2010PbPb5080CentScan();
182 else if (fCentralityClass=="3050") InitDplustoKpipi2011PbPb3050();
183 else AliFatal("Not yet implemented");
186 else if (fCollisionType==2) {
187 if (fCentralityClass=="0100") InitDplustoKpipi2013pPb0100();
188 if (fIsRapidityScan) {
189 if (fRapidityRange == "0804") InitDplustoKpipi2013pPb0100RapScan0804();
190 if (fRapidityRange == "0401") InitDplustoKpipi2013pPb0100RapScan0401();
191 if (fRapidityRange == "0101") InitDplustoKpipi2013pPb0100RapScan0101();
192 if (fRapidityRange == "0104") InitDplustoKpipi2013pPb0100RapScan0104();
193 if (fRapidityRange == "0408") InitDplustoKpipi2013pPb0100RapScan0408();
196 else AliFatal("Not yet implemented");
199 if (fCollisionType==0) {
200 if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
201 else InitDstartoD0pi2010pp();
203 else if (fCollisionType==1) {
204 if (fRunNumber == 10){
205 if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
206 else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
207 else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
208 else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
209 else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
210 else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
211 else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
212 else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
213 else AliFatal("Not yet implemented");
215 if (fRunNumber == 11){
216 if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
217 else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010();
218 else if (fCentralityClass=="010" && fIsCentScan) InitDstartoD0pi2011PbPb010CentScan();
219 else if (fCentralityClass=="1020") InitDstartoD0pi2011PbPb1020CentScan();
220 else if (fCentralityClass=="2030") InitDstartoD0pi2011PbPb2030CentScan();
221 else if (fCentralityClass=="3040") InitDstartoD0pi2011PbPb3040CentScan();
222 else if (fCentralityClass=="4050") InitDstartoD0pi2011PbPb4050CentScan();
223 else if (fCentralityClass=="5080") InitDstartoD0pi2010PbPb5080CentScan();
224 else if (fCentralityClass=="3050") InitDstartoD0pi2011PbPb3050();
225 else AliFatal("Not yet implemented");
228 else if (fCollisionType==2) {
229 if (fCentralityClass=="0100") InitDstartoD0pi2013pPb0100();
231 if (fCentralityClass=="020V0A") InitDstartoD0pi2013pPb020V0A();
232 if (fCentralityClass=="2040V0A") InitDstartoD0pi2013pPb2040V0A();
233 if (fCentralityClass=="4060V0A") InitDstartoD0pi2013pPb4060V0A();
234 if (fCentralityClass=="60100V0A") InitDstartoD0pi2013pPb60100V0A();
236 if (fCentralityClass=="020ZNA") InitDstartoD0pi2013pPb020ZNA();
237 if (fCentralityClass=="2040ZNA") InitDstartoD0pi2013pPb2040ZNA();
238 if (fCentralityClass=="4060ZNA") InitDstartoD0pi2013pPb4060ZNA();
239 if (fCentralityClass=="60100ZNA") InitDstartoD0pi2013pPb60100ZNA();
241 if (fIsRapidityScan) {
242 if (fRapidityRange == "0804") InitDstartoD0pi2013pPb0100RapScan0804();
243 if (fRapidityRange == "0401") InitDstartoD0pi2013pPb0100RapScan0401();
244 if (fRapidityRange == "0101") InitDstartoD0pi2013pPb0100RapScan0101();
245 if (fRapidityRange == "0104") InitDstartoD0pi2013pPb0100RapScan0104();
246 if (fRapidityRange == "0408") InitDstartoD0pi2013pPb0100RapScan0408();
249 else AliFatal("Not yet implemented");
252 if (fCollisionType==0) InitDstoKKpi2010pp();
253 else if (fCollisionType==1) {
254 if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
255 else AliFatal("Not yet implemented");
257 else if (fCollisionType==2) {
258 if (fCentralityClass=="0100") InitDstoKKpi2013pPb0100();
260 else AliFatal("Not yet implemented");
263 if (fCollisionType==0) InitLctopKpi2010pp();
264 else AliFatal("Not yet implemented");
267 if (fCollisionType==0) InitLctopK0S2010pp();
268 else AliFatal("Not yet implemented");
272 printf("Invalid decay type: %d\n",decay);
278 //--------------------------------------------------------------------------
279 void AliHFSystErr::InitD0toKpi2010pp() {
281 // D0->Kpi syst errors. Responsible: A. Rossi
286 fNorm = new TH1F("fNorm","fNorm",24,0,24);
287 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
290 fBR = new TH1F("fBR","fBR",24,0,24);
291 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
293 // Tracking efficiency
294 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
295 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
297 // Raw yield extraction
298 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
299 fRawYield->SetBinContent(1,1);
300 fRawYield->SetBinContent(2,0.22);
301 fRawYield->SetBinContent(3,0.1);
302 for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
303 for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
304 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
305 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
307 // Cuts efficiency (from cuts variation)
308 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
309 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
311 // PID efficiency (from PID/noPID)
312 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
313 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
314 fPIDEff->SetBinContent(2,0.05); // 5%
317 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
318 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
319 fMCPtShape->SetBinContent(1,0.03);
320 fMCPtShape->SetBinContent(2,0.03);
322 // particle-antiparticle
323 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",24,0,24);
324 // fPartAntipart->SetBinContent(1,1);
325 // for(Int_t i=2;i<=24;i++) fPartAntipart->SetBinContent(i,0.05);
329 //--------------------------------------------------------------------------
330 void AliHFSystErr::InitD0toKpi2010PbPb020() {
332 // D0->Kpi syst errors. Responsible: A. Rossi
333 // 2010 PbPb sample, 0-20 CC
337 fNorm = new TH1F("fNorm","fNorm",20,0,20);
338 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
341 fBR = new TH1F("fBR","fBR",20,0,20);
342 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
344 // Tracking efficiency
345 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
346 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
348 // Raw yield extraction
349 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
350 fRawYield->SetBinContent(1,0);
351 fRawYield->SetBinContent(2,0);
352 fRawYield->SetBinContent(3,0.08);
353 for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
354 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
356 // Cuts efficiency (from cuts variation)
357 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
358 fCutsEff->SetBinContent(1,0.);
359 fCutsEff->SetBinContent(2,0.);
360 fCutsEff->SetBinContent(3,0.13);
361 fCutsEff->SetBinContent(4,0.11);
362 for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
363 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
365 // PID efficiency (from PID/noPID)
366 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
367 for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
370 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
371 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
372 fMCPtShape->SetBinContent(3,0.04);
373 fMCPtShape->SetBinContent(4,0.02);
374 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
376 // // particle-antiparticle
377 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
378 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
379 // fPartAntipart->SetBinContent(3,0.10);
380 // fPartAntipart->SetBinContent(4,0.10);
381 // fPartAntipart->SetBinContent(7,0.10);
382 // fPartAntipart->SetBinContent(8,0.10);
386 //--------------------------------------------------------------------------
387 //--------------------------------------------------------------------------
388 void AliHFSystErr::InitD0toKpi2011PbPb07half() {
390 // D0->Kpi syst errors. Responsible: A. Rossi
391 // 2011 PbPb sample, 0-7.5 CC
395 fNorm = new TH1F("fNorm","fNorm",36,0,36);
396 for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
399 fBR = new TH1F("fBR","fBR",36,0,36);
400 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
402 // Tracking efficiency
403 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
404 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
405 for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
407 // Raw yield extraction
408 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
409 fRawYield->SetBinContent(1,0);
410 fRawYield->SetBinContent(2,0.1);
411 fRawYield->SetBinContent(3,0.05);
412 for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
413 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
414 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
415 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
417 // Cuts efficiency (from cuts variation)
418 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
419 fCutsEff->SetBinContent(1,0.);
420 fCutsEff->SetBinContent(2,0.15);
421 fCutsEff->SetBinContent(3,0.13);
422 fCutsEff->SetBinContent(4,0.11);
423 fCutsEff->SetBinContent(5,0.08);
424 for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
425 for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
427 // PID efficiency (from PID/noPID)
428 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
429 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
430 for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
433 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
434 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
435 fMCPtShape->SetBinContent(2,0.06);
436 fMCPtShape->SetBinContent(3,0.04);
437 fMCPtShape->SetBinContent(4,0.02);
438 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
439 for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
440 for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
442 // // particle-antiparticle
443 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",36,0,36);
444 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
445 // fPartAntipart->SetBinContent(3,0.10);
446 // fPartAntipart->SetBinContent(4,0.10);
447 // fPartAntipart->SetBinContent(7,0.10);
448 // fPartAntipart->SetBinContent(8,0.10);
453 //--------------------------------------------------------------------------
454 void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
456 // D0->Kpi syst errors. Responsible: D. Caffarri
457 // 2011 PbPb sample, 30-50 CC InPlane
459 InitD0toKpi2011PbPb07half();
460 // Raw yield extraction
461 // fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
462 fRawYield->SetBinContent(1,0);
463 fRawYield->SetBinContent(2,0);
464 fRawYield->SetBinContent(3,0.05);
465 fRawYield->SetBinContent(4,0.08);
466 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
467 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
468 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
469 for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
471 // Cuts efficiency (from cuts variation)
472 //fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
473 fCutsEff->SetBinContent(1,0.);
474 fCutsEff->SetBinContent(2,0.0);
475 fCutsEff->SetBinContent(3,0.10);
476 fCutsEff->SetBinContent(4,0.10);
477 fCutsEff->SetBinContent(5,0.10);
478 fCutsEff->SetBinContent(6,0.10);
479 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
480 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
481 for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
484 //--------------------------------------------------------------------------
485 void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
487 // D0->Kpi syst errors. Responsible: D. Caffarri
488 // 2011 PbPb sample, 30-50 CC OutOfPlane
490 InitD0toKpi2011PbPb07half();
491 // Raw yield extraction
492 //fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
493 fRawYield->SetBinContent(1,0);
494 fRawYield->SetBinContent(2,0.);
495 fRawYield->SetBinContent(3,0.05);
496 for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.07);
497 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
498 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
499 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
500 for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);// OUT OF MEASUREMENT RANGE
502 // Cuts efficiency (from cuts variation)
503 // fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
504 fCutsEff->SetBinContent(1,0.);
505 fCutsEff->SetBinContent(2,0.);
506 fCutsEff->SetBinContent(3,0.10);
507 fCutsEff->SetBinContent(4,0.10);
508 fCutsEff->SetBinContent(5,0.10);
509 fCutsEff->SetBinContent(6,0.10);
510 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
511 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
512 for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.0);// OUT OF MEASUREMENT RANGE
515 //--------------------------------------------------------------------------
516 void AliHFSystErr::InitD0toKpi2010PbPb4080() {
518 // D0->Kpi syst errors. Responsible: A. Rossi
519 // 2010 PbPb sample, 40-80 CC
523 fNorm = new TH1F("fNorm","fNorm",20,0,20);
524 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
527 fBR = new TH1F("fBR","fBR",20,0,20);
528 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
530 // Tracking efficiency
531 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
532 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
535 // Raw yield extraction
536 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
537 fRawYield->SetBinContent(1,0);
538 fRawYield->SetBinContent(2,0);
539 for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
540 //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
542 // Cuts efficiency (from cuts variation)
543 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
544 fCutsEff->SetBinContent(1,0.);
545 fCutsEff->SetBinContent(2,0.);
546 fCutsEff->SetBinContent(3,0.13);
547 fCutsEff->SetBinContent(4,0.11);
548 for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
549 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
551 // PID efficiency (from PID/noPID)
552 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
553 // for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
554 // for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
555 for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
558 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
559 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
560 // fMCPtShape->SetBinContent(3,0.04); Not set for peripherals (Raa Vs pt is flat)
561 // fMCPtShape->SetBinContent(4,0.02);
562 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
564 // // particle-antiparticle
565 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
566 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
571 //--------------------------------------------------------------------------
572 void AliHFSystErr::InitD0toKpi2011PbPb010() {
574 // D0->Kpi syst errors. Responsible: A. Rossi
575 // 2011 PbPb sample, 0-10 CC
578 AliInfo(" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
581 fNorm = new TH1F("fNorm","fNorm",24,0,24);
582 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
585 fBR = new TH1F("fBR","fBR",24,0,24);
586 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
588 // Tracking efficiency
589 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
590 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
593 // Raw yield extraction
594 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
595 fRawYield->SetBinContent(1,0);
596 fRawYield->SetBinContent(2,0.15);
597 fRawYield->SetBinContent(3,0.09);
598 fRawYield->SetBinContent(4,0.06);
599 fRawYield->SetBinContent(5,0.06);
600 fRawYield->SetBinContent(6,0.06);
601 for (Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.06);
602 for (Int_t i=9; i<=12; i++)fRawYield->SetBinContent(i,0.07);
603 for (Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.11);
604 for (Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.30);
606 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
607 fCutsEff->SetBinContent(1,0.);
608 fCutsEff->SetBinContent(2,0.15);
609 fCutsEff->SetBinContent(3,0.12);
610 fCutsEff->SetBinContent(4,0.05);
611 fCutsEff->SetBinContent(5,0.05);
612 fCutsEff->SetBinContent(6,0.05);
613 for (Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i,0.05);
614 for (Int_t i=9; i<=12; i++)fCutsEff->SetBinContent(i,0.05);
615 for (Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i,0.05);
616 for (Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i,0.05);
619 // PID efficiency (from PID/noPID)
620 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
621 // for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
622 // for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
623 for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
626 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
627 fMCPtShape->SetBinContent(1,0.);
628 fMCPtShape->SetBinContent(2,0.15);
629 fMCPtShape->SetBinContent(3,0.12);
630 fMCPtShape->SetBinContent(4,0.02);
631 for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
633 // // particle-antiparticle
634 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
635 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
640 //--------------------------------------------------------------------------
641 void AliHFSystErr::InitD0toKpi2011PbPb3050() {
643 // D0->Kpi syst errors. Responsible: D. Caffarri
644 // 2011 PbPb sample, 30-50 CC
648 fNorm = new TH1F("fNorm","fNorm",20,0,20);
649 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
652 fBR = new TH1F("fBR","fBR",20,0,20);
653 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
655 // Tracking efficiency
656 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
657 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
660 // Raw yield extraction
661 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
662 fRawYield->SetBinContent(1,0);
663 fRawYield->SetBinContent(2,0.102);
664 fRawYield->SetBinContent(3,0.071);
665 fRawYield->SetBinContent(4,0.052);
666 fRawYield->SetBinContent(5,0.063);
667 fRawYield->SetBinContent(6,0.071);
668 fRawYield->SetBinContent(7,0.071);
669 fRawYield->SetBinContent(8,0.071);
670 fRawYield->SetBinContent(9,0.071);
671 fRawYield->SetBinContent(10,0.071);
672 fRawYield->SetBinContent(11,0.071);
673 fRawYield->SetBinContent(12,0.071);
675 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.112);
676 for(Int_t i=17;i<=20;i++) fRawYield->SetBinContent(i,0);
678 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
679 fCutsEff->SetBinContent(1,0.);
680 fCutsEff->SetBinContent(2,0.10);
681 fCutsEff->SetBinContent(3,0.10);
682 fCutsEff->SetBinContent(4,0.10);
683 fCutsEff->SetBinContent(5,0.10);
684 fCutsEff->SetBinContent(6,0.10);
685 for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
686 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
689 // PID efficiency (from PID/noPID)
690 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
691 // for(Int_t i=3;i<=6;i++) fPIDEff->SetBinContent(i,0.10);
692 // for(Int_t i=7;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
693 for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
696 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",0,0,20);
697 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
698 fMCPtShape->SetBinContent(2,0.06);
699 fMCPtShape->SetBinContent(3,0.04);
700 fMCPtShape->SetBinContent(4,0.02);
701 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
702 for(Int_t i=17;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
703 // // particle-antiparticle
704 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
705 // for(Int_t i=3;i<=12;i++) fPartAntipart->SetBinContent(i,0.05);
711 //--------------------------------------------------------------------------
712 void AliHFSystErr::InitD0toKpi2010ppLowEn() {
714 // D0->Kpi syst errors. Low energy run
715 // 2011 2.76 TeV pp sample
717 AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
720 fNorm = new TH1F("fNorm","fNorm",20,0,20);
721 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 4% error on sigmaV0and
724 fBR = new TH1F("fBR","fBR",20,0,20);
725 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
727 // Tracking efficiency
728 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
729 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
731 // Raw yield extraction
732 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
733 fRawYield->SetBinContent(1,1);
734 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
736 // Cuts efficiency (from cuts variation)
737 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
738 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
739 fCutsEff->SetBinContent(2,0.20);
740 for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
743 // PID efficiency (from PID/noPID)
744 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
745 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15); // 10%
746 // fPIDEff->SetBinContent(2,0.20);
747 for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
750 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
751 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
752 fMCPtShape->SetBinContent(1,0.03);
753 fMCPtShape->SetBinContent(2,0.03);
755 // // particle-antiparticle
756 // fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
757 // fPartAntipart->SetBinContent(1,1);
758 // fPartAntipart->SetBinContent(2,1);
759 // for(Int_t i=3;i<=6;i++) fPartAntipart->SetBinContent(i,0.08);
760 // for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.);
764 //_________________________________________________________________________
765 void AliHFSystErr::InitD0toKpi2013pPb0100(){
767 // D0->Kpi syst errors. p-Pb data sample
769 AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
771 fNorm = new TH1F("fNorm","fNorm",24,0,24);
772 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
775 fBR = new TH1F("fBR","fBR",24,0,24);
776 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013); // 1.2% PDG2010
778 // Tracking efficiency
779 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
780 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06); // 3% per track
782 // Raw yield extraction
783 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
784 fRawYield->SetBinContent(1,0.0);
785 fRawYield->SetBinContent(2,0.08);
786 for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.04);
787 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.07);
788 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.11);
790 // Cuts efficiency (from cuts variation)
791 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
792 fCutsEff->SetBinContent(1,0.0);
793 fCutsEff->SetBinContent(2,0.08);
794 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
796 // PID efficiency (from PID/noPID)
797 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
798 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
801 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
802 fMCPtShape->SetBinContent(1,0.0);
803 fMCPtShape->SetBinContent(2,0.02);
804 for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
806 // particle-antiparticle
808 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
809 for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
817 //--------------------------------------------------------------------------
818 void AliHFSystErr::InitDplustoKpipi2010pp() {
820 // D+->Kpipi syst errors. Responsible: R. Bala
826 fNorm = new TH1F("fNorm","fNorm",24,0,24);
827 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
830 fBR = new TH1F("fBR","fBR",24,0,24);
831 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
833 // Tracking efficiency
834 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
835 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
838 // Raw yield extraction
839 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
840 fRawYield->SetBinContent(1,1);
841 fRawYield->SetBinContent(2,0.25);
842 fRawYield->SetBinContent(3,0.25);
843 fRawYield->SetBinContent(4,0.20);
844 fRawYield->SetBinContent(5,0.09);
845 fRawYield->SetBinContent(6,0.09);
846 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
847 for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
849 // Cuts efficiency (from cuts variation)
850 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
851 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
853 // PID efficiency (from PID/noPID)
854 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
855 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 5%
856 fPIDEff->SetBinContent(1,0.15); // 15%
857 fPIDEff->SetBinContent(2,0.15); // 15%
858 fPIDEff->SetBinContent(3,0.15); // 15%
859 fPIDEff->SetBinContent(4,0.15); // 15%
860 for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10); // 5%
863 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
864 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
865 fMCPtShape->SetBinContent(1,0.03);
866 fMCPtShape->SetBinContent(2,0.03);
869 // particle-antiparticle
871 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
872 fPartAntipart->SetBinContent(1,1);
873 fPartAntipart->SetBinContent(2,1);
874 fPartAntipart->SetBinContent(3,0.12);
875 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
880 //--------------------------------------------------------------------------
881 void AliHFSystErr::InitDstoKKpi2010pp() {
883 // D+s->KKpi syst errors. Responsible: G.M. Innocenti
889 fNorm = new TH1F("fNorm","fNorm",12,0,12);
890 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
893 fBR = new TH1F("fBR","fBR",12,0,12);
894 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
896 // Tracking efficiency
897 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
898 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
901 // Raw yield extraction
902 fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
903 fRawYield->SetBinContent(1,1);
904 fRawYield->SetBinContent(2,1);
905 fRawYield->SetBinContent(3,0.20);
906 fRawYield->SetBinContent(4,0.20);
907 fRawYield->SetBinContent(5,0.15);
908 fRawYield->SetBinContent(6,0.15);
909 fRawYield->SetBinContent(7,0.15);
910 fRawYield->SetBinContent(8,0.15);
911 fRawYield->SetBinContent(9,0.20);
912 fRawYield->SetBinContent(10,0.20);
913 fRawYield->SetBinContent(11,0.20);
914 fRawYield->SetBinContent(12,0.20);
916 // Cuts efficiency (from cuts variation)
917 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
918 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15); // 15%
920 // PID efficiency (from PID/noPID)
921 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
922 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07); // 7%
925 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
926 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
927 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
928 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
929 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
930 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
933 // particle-antiparticle
935 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
936 fPartAntipart->SetBinContent(1,1);
937 fPartAntipart->SetBinContent(2,1);
938 fPartAntipart->SetBinContent(3,0.12);
939 for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
945 //--------------------------------------------------------------------------
946 void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
948 // D+->Kpipi syst errors. Responsible: E. Bruna
949 // 2011 PbPb sample, 0-7.5% CC
953 fNorm = new TH1F("fNorm","fNorm",36,0,36);
954 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
957 fBR = new TH1F("fBR","fBR",36,0,36);
958 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
960 // Tracking efficiency
961 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
962 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
964 // Raw yield extraction
965 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
966 for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
967 fRawYield->SetBinContent(4,0.30);
968 fRawYield->SetBinContent(5,0.20);
969 fRawYield->SetBinContent(6,0.20);
970 for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
971 for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
972 for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
973 for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
974 for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
977 // Cuts efficiency (from cuts variation)
978 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
979 for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
980 for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);
983 // PID efficiency (from PID/noPID)
984 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
985 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
987 // MC dN/dpt (24/7/2012)
988 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
989 for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
990 for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
991 // for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
994 // particle-antiparticle
996 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
997 fPartAntipart->SetBinContent(1,1);
998 fPartAntipart->SetBinContent(2,1);
999 fPartAntipart->SetBinContent(3,0.12);
1000 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1005 //_________________________________________________________________
1006 void AliHFSystErr::InitDplustoKpipi2011PbPb010() {
1008 // D+->Kpipi syst errors. Responsible: E. Bruna, C. Bedda
1009 // 2011 PbPb sample, 0-10% CC
1013 fNorm = new TH1F("fNorm","fNorm",36,0,36);
1014 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1017 fBR = new TH1F("fBR","fBR",36,0,36);
1018 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1020 // Tracking efficiency
1021 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
1022 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1024 // Raw yield extraction
1025 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1026 for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i,0.); //5 to 10%
1027 fRawYield->SetBinContent(4,0.10);
1028 fRawYield->SetBinContent(5,0.08);
1029 fRawYield->SetBinContent(6,0.08);
1030 for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.08);
1031 for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
1032 for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.08);
1033 for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
1034 for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.08);
1037 // Cuts efficiency (from cuts variation)
1038 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1039 for(Int_t i=1; i<=36; i++) fCutsEff->SetBinContent(i,0.10);
1042 // PID efficiency (from PID/noPID)
1043 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1044 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1047 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1048 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.);
1049 fMCPtShape->SetBinContent(4,0.06);
1050 fMCPtShape->SetBinContent(5,0.02);
1051 for (Int_t i=6; i<=36; i++) fMCPtShape->SetBinContent(i,0.01);
1052 // for(Int_t iBin=13; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1055 // particle-antiparticle
1057 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1058 fPartAntipart->SetBinContent(1,1);
1059 fPartAntipart->SetBinContent(2,1);
1060 fPartAntipart->SetBinContent(3,0.12);
1061 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1067 //___________________________________________________________________
1068 void AliHFSystErr::InitDplustoKpipi2011PbPb3050() {
1070 // D+->Kpipi syst errors. Responsible: R. Bala
1071 // 2011 PbPb sample, 30-50% CC
1073 AliInfo(" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
1075 fNorm = new TH1F("fNorm","fNorm",16,0,16);
1076 for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1079 fBR = new TH1F("fBR","fBR",16,0,16);
1080 for(Int_t i=1;i<=16;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1082 // Tracking efficiency
1083 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
1084 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1086 // Raw yield extraction
1087 fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
1088 for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.09); //5 to 10%
1089 for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.12);
1092 // Cuts efficiency (from cuts variation)
1093 fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
1094 for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
1095 for(Int_t i=12; i<=16; i++) fCutsEff->SetBinContent(i,0.14);
1097 // PID efficiency (from PID/noPID)
1098 fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
1099 for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1101 // MC dN/dpt (24/7/2012)
1102 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
1103 fMCPtShape->SetBinContent(3,0.10);
1104 fMCPtShape->SetBinContent(4,0.04);
1105 for(Int_t iBin=5; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.02);
1108 // particle-antiparticle
1110 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1111 fPartAntipart->SetBinContent(1,1);
1112 fPartAntipart->SetBinContent(2,1);
1113 fPartAntipart->SetBinContent(3,0.12);
1114 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1121 //--------------------------------------------------------------------------
1122 void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
1124 // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
1125 // 2011 PbPb sample, 0-7.5% CC
1129 fNorm = new TH1F("fNorm","fNorm",12,0,12);
1130 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1133 fBR = new TH1F("fBR","fBR",12,0,12);
1134 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
1136 // Tracking efficiency
1137 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
1138 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1140 // Raw yield extraction
1141 fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
1142 for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30);
1143 for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
1147 // Cuts efficiency (from cuts variation)
1148 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
1149 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
1151 // PID efficiency (from PID/noPID)
1152 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
1153 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
1156 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
1157 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
1158 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
1159 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
1160 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
1161 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
1163 // particle-antiparticle
1165 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
1166 fPartAntipart->SetBinContent(1,1);
1167 fPartAntipart->SetBinContent(2,1);
1168 fPartAntipart->SetBinContent(3,0.12);
1169 for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1175 //_________________________________________________________________________
1176 void AliHFSystErr::InitDstoKKpi2013pPb0100(){
1179 AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
1182 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1183 for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); //
1186 fBR = new TH1F("fBR","fBR",12,0,12);
1187 for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
1189 // Tracking efficiency
1190 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1191 for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1193 // Raw yield extraction
1194 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1195 for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
1196 for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
1197 for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
1198 for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
1200 // Cuts efficiency (from cuts variation)
1201 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1202 for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //
1204 // PID efficiency (from PID/noPID)
1205 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1206 for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10);
1207 for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
1210 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1211 for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
1213 // particle-antiparticle
1215 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1216 fPartAntipart->SetBinContent(1,1);
1217 fPartAntipart->SetBinContent(2,1);
1218 fPartAntipart->SetBinContent(3,0.12);
1219 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1225 //--------------------------------------------------------------------------
1226 void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
1228 // D+->Kpipi syst errors. Responsible: ??
1229 // 2010 PbPb sample, 0-20 CC
1233 fNorm = new TH1F("fNorm","fNorm",20,0,20);
1234 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
1237 fBR = new TH1F("fBR","fBR",20,0,20);
1238 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1240 // Tracking efficiency
1241 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1242 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1244 // Raw yield extraction
1245 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1246 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10); //5 to 10%
1247 // fRawYield->SetBinContent(5,0.23);
1248 //fRawYield->SetBinContent(6,0.23);
1249 fRawYield->SetBinContent(7,0.20);
1250 fRawYield->SetBinContent(8,0.20);
1251 fRawYield->SetBinContent(9,0.15);
1252 fRawYield->SetBinContent(10,0.15);
1253 fRawYield->SetBinContent(11,0.15);
1254 fRawYield->SetBinContent(12,0.15);
1256 // Cuts efficiency (from cuts variation)
1257 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1258 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
1260 // PID efficiency (from PID/noPID)
1261 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1262 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1264 // MC dN/dpt (2/2/2012)
1265 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1266 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
1267 for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1268 for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1269 for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
1272 // particle-antiparticle
1274 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1275 fPartAntipart->SetBinContent(1,1);
1276 fPartAntipart->SetBinContent(2,1);
1277 fPartAntipart->SetBinContent(3,0.12);
1278 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1284 //--------------------------------------------------------------------------
1285 void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
1287 // D+->Kpipi syst errors. Responsible: ??
1288 // 2010 PbPb sample, 40-80 CC
1293 fNorm = new TH1F("fNorm","fNorm",20,0,20);
1294 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1297 fBR = new TH1F("fBR","fBR",20,0,20);
1298 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
1300 // Tracking efficiency
1301 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1302 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1305 // Raw yield extraction
1306 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1307 fRawYield->SetBinContent(1,1);
1308 fRawYield->SetBinContent(2,1);
1309 fRawYield->SetBinContent(3,1);
1310 fRawYield->SetBinContent(4,0.15);
1311 fRawYield->SetBinContent(5,0.05);
1312 fRawYield->SetBinContent(6,0.05);
1313 fRawYield->SetBinContent(7,0.15);
1314 fRawYield->SetBinContent(8,0.15);
1315 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
1316 for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
1318 // Cuts efficiency (from cuts variation)
1319 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1320 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1322 // PID efficiency (from PID/noPID)
1323 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1324 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1325 fPIDEff->SetBinContent(3,0.13); // 13%
1328 // MC dN/dpt (2/2/2012)
1329 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1330 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
1331 for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
1332 for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1333 for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
1336 // particle-antiparticle
1338 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1339 fPartAntipart->SetBinContent(1,1);
1340 fPartAntipart->SetBinContent(2,1);
1341 fPartAntipart->SetBinContent(3,0.12);
1342 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1347 //--------------------------------------------------------------------------
1348 void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
1351 // D+->Kpipi syst errors. Responsible: R. Bala
1352 // 2011 2.76 TeV pp sample
1354 AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1357 fNorm = new TH1F("fNorm","fNorm",20,0,20);
1358 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
1361 fBR = new TH1F("fBR","fBR",20,0,20);
1362 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1364 // Tracking efficiency
1365 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1366 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
1368 // Raw yield extraction
1369 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1370 fRawYield->SetBinContent(1,1);
1371 fRawYield->SetBinContent(2,1);
1372 for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10); //5 to 10%
1373 fRawYield->SetBinContent(7,0.15);
1374 fRawYield->SetBinContent(8,0.15);
1375 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055); //5 to 10%
1377 // Cuts efficiency (from cuts variation)
1378 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1379 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
1381 // PID efficiency (from PID/noPID)
1382 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1383 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 5%
1384 fPIDEff->SetBinContent(3,0.10); // 13%
1385 fPIDEff->SetBinContent(4,0.10); // 13%
1387 // MC dN/dpt (copied from D0 : will update later)
1388 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1389 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1390 fMCPtShape->SetBinContent(1,0.03);
1391 fMCPtShape->SetBinContent(2,0.03);
1396 //_________________________________________________________________________
1397 void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
1400 AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
1402 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1403 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
1405 fBR = new TH1F("fBR","fBR",20,0,24);
1406 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
1408 // Tracking efficiency
1409 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1410 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1412 // Raw yield extraction
1413 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1414 fRawYield->SetBinContent(1,0.);
1415 fRawYield->SetBinContent(2,0.1);
1416 fRawYield->SetBinContent(3,0.08);
1417 fRawYield->SetBinContent(4,0.05);
1418 fRawYield->SetBinContent(5,0.05);
1419 fRawYield->SetBinContent(6,0.05);
1420 fRawYield->SetBinContent(7,0.05);
1421 fRawYield->SetBinContent(8,0.05);
1422 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1423 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
1424 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
1426 // Cuts efficiency (from cuts variation)
1427 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1428 fCutsEff->SetBinContent(1,0.);
1429 fCutsEff->SetBinContent(2,0.1);
1430 fCutsEff->SetBinContent(3,0.1);
1431 fCutsEff->SetBinContent(4,0.06);
1432 fCutsEff->SetBinContent(5,0.06);
1433 fCutsEff->SetBinContent(6,0.06);
1434 fCutsEff->SetBinContent(7,0.05);
1435 fCutsEff->SetBinContent(8,0.05);
1436 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
1437 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
1438 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
1440 // PID efficiency (from PID/noPID)
1441 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1442 fPIDEff->SetBinContent(1,1);
1443 fPIDEff->SetBinContent(2,0.1);
1444 for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0); //
1447 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1448 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
1450 // particle-antiparticle
1452 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1453 fPartAntipart->SetBinContent(1,1);
1454 fPartAntipart->SetBinContent(2,1);
1455 fPartAntipart->SetBinContent(3,0.12);
1456 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
1463 //--------------------------------------------------------------------------
1464 void AliHFSystErr::InitDstartoD0pi2010pp() {
1466 // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
1471 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1472 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
1475 fBR = new TH1F("fBR","fBR",24,0,24);
1476 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1478 // Tracking efficiency
1479 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1480 fTrackingEff->SetBinContent(1,1.0);
1481 fTrackingEff->SetBinContent(2,0.13); // 10% (ITSsa) \oplus 8% (4% per ITSTPC track)
1482 fTrackingEff->SetBinContent(3,0.12);
1483 fTrackingEff->SetBinContent(3,0.12);
1484 for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
1487 // Raw yield extraction
1488 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1489 fRawYield->SetBinContent(1,1.0);
1490 fRawYield->SetBinContent(2,0.10);
1491 fRawYield->SetBinContent(3,0.04);
1492 fRawYield->SetBinContent(4,0.03);
1493 fRawYield->SetBinContent(5,0.03);
1494 fRawYield->SetBinContent(6,0.05);
1495 fRawYield->SetBinContent(7,0.05);
1496 fRawYield->SetBinContent(8,0.05);
1497 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04); //4%
1498 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09); //4%
1499 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
1501 // Cuts efficiency (from cuts variation)
1502 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1503 fCutsEff->SetBinContent(2,0.22);
1504 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1506 // PID efficiency (from PID/noPID)
1507 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1508 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1511 // MC dN/dpt (copied from D0 : will update later)
1512 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1513 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
1514 fMCPtShape->SetBinContent(1,0.03);
1515 fMCPtShape->SetBinContent(2,0.03);
1521 //--------------------------------------------------------------------------
1522 void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
1525 // D+->Kpipi syst errors. Responsible: A. Grelli
1526 // 2011 2.76 TeV pp sample
1528 AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1531 fNorm = new TH1F("fNorm","fNorm",20,0,20);
1532 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
1535 fBR = new TH1F("fBR","fBR",20,0,20);
1536 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1538 // Tracking efficiency
1539 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
1540 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
1542 // Raw yield extraction
1543 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1544 fRawYield->SetBinContent(1,1);
1545 fRawYield->SetBinContent(2,1);
1546 fRawYield->SetBinContent(3,0.14);
1547 fRawYield->SetBinContent(4,0.14);
1548 fRawYield->SetBinContent(5,0.12);
1549 fRawYield->SetBinContent(6,0.12);
1550 fRawYield->SetBinContent(7,0.06);
1551 fRawYield->SetBinContent(8,0.06);
1552 fRawYield->SetBinContent(9,0.08);
1553 fRawYield->SetBinContent(10,0.08);
1554 fRawYield->SetBinContent(11,0.08);
1555 fRawYield->SetBinContent(12,0.08);
1556 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
1558 // Cuts efficiency (from cuts variation)
1559 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
1560 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
1561 fCutsEff->SetBinContent(3,0.15);
1562 fCutsEff->SetBinContent(4,0.15);
1563 fCutsEff->SetBinContent(5,0.15);
1564 fCutsEff->SetBinContent(6,0.15);
1565 fCutsEff->SetBinContent(7,0.10);
1566 fCutsEff->SetBinContent(8,0.10);
1567 fCutsEff->SetBinContent(9,0.10);
1568 fCutsEff->SetBinContent(10,0.10);
1569 fCutsEff->SetBinContent(11,0.10);
1570 fCutsEff->SetBinContent(12,0.10);
1572 // PID efficiency (from PID/noPID)
1573 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
1574 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
1577 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
1578 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
1579 fMCPtShape->SetBinContent(1,0.03);
1580 fMCPtShape->SetBinContent(2,0.03);
1585 //------------------------------------------------------------------------
1586 void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
1588 // D*+->D0pi syst errors. Responsible: A. Grelli
1589 // 2010 PbPb sample, 0-20 CC
1592 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1595 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1596 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
1599 fBR = new TH1F("fBR","fBR",24,0,24);
1600 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1602 // Tracking efficiency
1603 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1604 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1607 // Raw yield extraction
1608 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1609 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
1610 fRawYield->SetBinContent(3,0.2);
1611 fRawYield->SetBinContent(4,0.2);
1612 fRawYield->SetBinContent(5,0.2);
1613 fRawYield->SetBinContent(6,0.2);
1615 // Cuts efficiency (from cuts variation)
1616 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1617 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1618 fCutsEff->SetBinContent(4,0.15);
1619 fCutsEff->SetBinContent(5,0.15);
1620 fCutsEff->SetBinContent(6,0.15);
1622 // PID efficiency (from PID/noPID)
1623 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1624 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1627 // MC dN/dpt (from study on D* pt shape)
1628 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1629 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
1630 fMCPtShape->SetBinContent(4,0.025);
1631 fMCPtShape->SetBinContent(5,0.025);
1632 fMCPtShape->SetBinContent(6,0.025);
1633 fMCPtShape->SetBinContent(7,0.04);
1634 fMCPtShape->SetBinContent(8,0.04);
1635 fMCPtShape->SetBinContent(9,0.03);
1636 fMCPtShape->SetBinContent(10,0.03);
1637 fMCPtShape->SetBinContent(11,0.03);
1638 fMCPtShape->SetBinContent(12,0.03);
1643 // ----------------------------- 2011 ----------------------------------
1644 void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
1646 // D*+->D0pi syst errors. Responsible: A. Grelli
1647 // 2011 PbPb sample, 0-7.5 CC
1650 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1654 fNorm = new TH1F("fNorm","fNorm",36,0,36);
1655 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1658 fBR = new TH1F("fBR","fBR",36,0,36);
1659 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1661 // Tracking efficiency
1662 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1663 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1666 // Raw yield extraction
1667 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1668 for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05); //4%
1669 fRawYield->SetBinContent(4,0.2);
1670 fRawYield->SetBinContent(5,0.10);
1671 fRawYield->SetBinContent(6,0.10);
1672 fRawYield->SetBinContent(7,0.08);
1673 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15); //4%
1675 // Cuts efficiency (from cuts variation)
1676 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1677 for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1679 // PID efficiency (from PID/noPID)
1680 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1681 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1682 fPIDEff->SetBinContent(4,0.09);
1684 // MC dN/dpt (from study on D* pt shape)
1685 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1686 for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
1687 fMCPtShape->SetBinContent(4,0.015);
1688 fMCPtShape->SetBinContent(5,0.015);
1689 fMCPtShape->SetBinContent(6,0.015);
1690 fMCPtShape->SetBinContent(7,0.02);
1691 fMCPtShape->SetBinContent(8,0.02);
1692 fMCPtShape->SetBinContent(9,0.03);
1693 fMCPtShape->SetBinContent(10,0.03);
1694 fMCPtShape->SetBinContent(11,0.03);
1695 fMCPtShape->SetBinContent(12,0.03);
1702 //-------------------------------------------------------------------------
1703 void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
1705 // D*+->D0pi syst errors. Responsible: A. Grelli
1706 // 2010 PbPb sample, 20-40 CC
1709 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1712 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1713 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10); // 10% error on sigmaV0and
1716 fBR = new TH1F("fBR","fBR",24,0,24);
1717 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1719 // Tracking efficiency
1720 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1721 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1724 // Raw yield extraction
1725 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1726 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15); //4%
1728 // Cuts efficiency (from cuts variation)
1729 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1730 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1732 // PID efficiency (from PID/noPID)
1733 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1734 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04); // 3%
1737 // MC dN/dpt (copied from D0 : will update later)
1738 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1739 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
1740 fMCPtShape->SetBinContent(1,0.03);
1741 fMCPtShape->SetBinContent(2,0.03);
1747 //--------------------------------------------------------------------------
1748 void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
1750 // D*+->D0pi syst errors. Responsible: A. Grelli
1751 // 2010 PbPb sample, 40-80 CC
1754 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1757 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1758 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1761 fBR = new TH1F("fBR","fBR",24,0,24);
1762 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1764 // Tracking efficiency
1765 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1766 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1769 // Raw yield extraction
1770 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1771 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2); //4%
1772 fRawYield->SetBinContent(1,1);
1773 fRawYield->SetBinContent(2,0.15);
1774 fRawYield->SetBinContent(3,0.15);
1775 fRawYield->SetBinContent(4,0.15);
1776 fRawYield->SetBinContent(5,0.15);
1777 fRawYield->SetBinContent(6,0.10);
1778 fRawYield->SetBinContent(7,0.10);
1779 fRawYield->SetBinContent(8,0.10);
1780 fRawYield->SetBinContent(9,0.11);
1781 fRawYield->SetBinContent(10,0.11);
1782 fRawYield->SetBinContent(11,0.11);
1783 fRawYield->SetBinContent(12,0.11);
1784 fRawYield->SetBinContent(13,0.08);
1785 fRawYield->SetBinContent(14,0.08);
1786 fRawYield->SetBinContent(15,0.08);
1787 fRawYield->SetBinContent(16,0.08);
1790 // Cuts efficiency (from cuts variation)
1791 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1792 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1794 // PID efficiency (from PID/noPID)
1795 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1796 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1799 // MC dN/dpt (copied from D0 : will update later)
1800 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1801 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1802 fMCPtShape->SetBinContent(2,0.05);
1803 fMCPtShape->SetBinContent(3,0.05);
1804 fMCPtShape->SetBinContent(4,0.05);
1805 fMCPtShape->SetBinContent(5,0.04);
1806 fMCPtShape->SetBinContent(6,0.02);
1807 fMCPtShape->SetBinContent(7,0.04);
1808 fMCPtShape->SetBinContent(8,0.04);
1814 // ----------------------------- 2011 ----------------------------------
1815 void AliHFSystErr::InitDstartoD0pi2011PbPb010() {
1817 // D*+->D0pi syst errors. Responsible: A. Grelli
1818 // 2011 PbPb sample, 0-10% CC
1821 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
1825 fNorm = new TH1F("fNorm","fNorm",36,0,36);
1826 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
1829 fBR = new TH1F("fBR","fBR",36,0,36);
1830 for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1832 // Tracking efficiency
1833 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
1834 for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1837 // Raw yield extraction
1838 fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
1839 for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i, 0.);
1840 fRawYield->SetBinContent(4,0.10);
1841 fRawYield->SetBinContent(5,0.08);
1842 fRawYield->SetBinContent(6,0.05);
1843 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
1844 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1845 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
1846 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
1847 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10); //4%
1849 // Cuts efficiency (from cuts variation)
1850 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
1851 for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i, 0.);
1852 fCutsEff->SetBinContent(4,0.15);
1853 for(Int_t i=5;i<=36;i++) fCutsEff->SetBinContent(i,0.10); // 10%
1855 // PID efficiency (from PID/noPID)
1856 fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
1857 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1859 // MC dN/dpt (from study on D* pt shape)
1860 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
1861 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i, 0.);
1862 fMCPtShape->SetBinContent(4,0.06);
1863 fMCPtShape->SetBinContent(5,0.03);
1864 for(Int_t i=5;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
1870 //----------------------------------------------------------------------
1871 void AliHFSystErr::InitDstartoD0pi2011PbPb3050() {
1873 // D*+->D0pi syst errors. Responsible: A. Grelli
1874 // 2011 PbPb sample, 30-50 CC
1877 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
1880 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1881 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
1884 fBR = new TH1F("fBR","fBR",24,0,24);
1885 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1887 // Tracking efficiency
1888 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
1889 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
1892 // Raw yield extraction
1893 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1894 for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1); //4%
1895 fRawYield->SetBinContent(1,0.12);
1896 fRawYield->SetBinContent(2,0.12);
1897 fRawYield->SetBinContent(3,0.12);
1898 fRawYield->SetBinContent(4,0.06);
1899 fRawYield->SetBinContent(5,0.06);
1900 fRawYield->SetBinContent(6,0.05);
1901 fRawYield->SetBinContent(7,0.08);
1902 fRawYield->SetBinContent(8,0.05);
1903 fRawYield->SetBinContent(9,0.05);
1904 fRawYield->SetBinContent(10,0.05);
1905 fRawYield->SetBinContent(11,0.05);
1906 fRawYield->SetBinContent(12,0.05);
1907 fRawYield->SetBinContent(13,0.05);
1908 fRawYield->SetBinContent(14,0.05);
1909 fRawYield->SetBinContent(15,0.05);
1910 fRawYield->SetBinContent(16,0.05);
1913 // Cuts efficiency (from cuts variation)
1914 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1915 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.05); // 10%
1916 fCutsEff->SetBinContent(1,0.16);
1917 fCutsEff->SetBinContent(2,0.16);
1918 fCutsEff->SetBinContent(3,0.16);
1919 fCutsEff->SetBinContent(4,0.10);
1920 fCutsEff->SetBinContent(5,0.10);
1921 fCutsEff->SetBinContent(6,0.10);
1922 // PID efficiency (from PID/noPID)
1923 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1924 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
1927 // MC dN/dpt (copied from D0 : will update later)
1928 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1929 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1930 fMCPtShape->SetBinContent(2,0.11);
1931 fMCPtShape->SetBinContent(3,0.11);
1932 fMCPtShape->SetBinContent(4,0.03);
1933 fMCPtShape->SetBinContent(5,0.02);
1934 fMCPtShape->SetBinContent(6,0.01);
1935 fMCPtShape->SetBinContent(7,0.01);
1936 fMCPtShape->SetBinContent(8,0.01);
1940 //_________________________________________________________________________
1941 void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
1942 // Dstar syst in pPb 2013 MB
1944 AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
1947 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1948 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); //
1951 fBR = new TH1F("fBR","fBR",24,0,24);
1952 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1954 // Tracking efficiency
1955 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1956 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
1958 // Raw yield extraction
1959 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1960 fRawYield->SetBinContent(1,1.0);
1961 fRawYield->SetBinContent(2,0.08);
1962 fRawYield->SetBinContent(3,0.02);
1963 fRawYield->SetBinContent(4,0.02);
1964 fRawYield->SetBinContent(5,0.05);
1965 fRawYield->SetBinContent(6,0.02);
1966 fRawYield->SetBinContent(7,0.04);
1967 fRawYield->SetBinContent(8,0.04);
1968 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
1969 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
1970 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
1972 // Cuts efficiency (from cuts variation)
1973 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
1974 fCutsEff->SetBinContent(1,0.);
1975 fCutsEff->SetBinContent(2,0.10);
1976 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
1978 // PID efficiency (from PID/noPID)
1979 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
1980 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
1983 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
1984 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
1985 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
1986 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
1988 // particle-antiparticle
1990 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
1991 fPartAntipart->SetBinContent(1,1);
1992 fPartAntipart->SetBinContent(2,1);
1993 fPartAntipart->SetBinContent(3,0.12);
1994 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2001 //--------------------------------------------------------------------------
2002 void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
2003 // define errors for RAA vs. centrality
2004 InitD0toKpi2010PbPb020();
2005 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2006 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
2007 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2009 //--------------------------------------------------------------------------
2010 void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
2011 // define errors for RAA vs. centrality
2012 InitD0toKpi2010PbPb020();
2013 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2014 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
2015 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2017 //--------------------------------------------------------------------------
2018 void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
2019 // define errors for RAA vs. centrality
2020 InitD0toKpi2010PbPb4080();
2021 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2022 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.14);
2023 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2025 //--------------------------------------------------------------------------
2026 void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
2027 // define errors for RAA vs. centrality
2028 InitD0toKpi2010PbPb4080();
2029 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
2030 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.11);
2031 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2033 //--------------------------------------------------------------------------
2034 void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
2035 // define errors for RAA vs. centrality
2036 InitD0toKpi2010PbPb4080();
2037 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2038 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
2039 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2041 //--------------------------------------------------------------------------
2042 void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
2043 // define errors for RAA vs. centrality
2044 InitDplustoKpipi2010PbPb020();
2045 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
2046 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
2049 //--------------------------------------------------------------------------
2050 void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
2051 // define errors for RAA vs. centrality
2052 InitDplustoKpipi2010PbPb020();
2053 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
2054 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2056 //--------------------------------------------------------------------------
2057 void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
2058 // define errors for RAA vs. centrality
2059 InitDplustoKpipi2010PbPb020();
2060 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2061 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
2063 //--------------------------------------------------------------------------
2064 void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
2065 // define errors for RAA vs. centrality
2066 InitDplustoKpipi2010PbPb4080();
2067 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2068 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2070 //--------------------------------------------------------------------------
2071 void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
2072 // define errors for RAA vs. centrality
2073 InitDplustoKpipi2010PbPb4080();
2074 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
2075 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
2078 //--------------------------------------------------------------------------
2079 void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
2080 // define errors for RAA vs. centrality
2081 InitDstartoD0pi2010PbPb020();
2082 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
2083 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
2085 //--------------------------------------------------------------------------
2086 void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
2087 // define errors for RAA vs. centrality
2088 InitDstartoD0pi2010PbPb020();
2089 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
2090 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
2092 //--------------------------------------------------------------------------
2093 void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
2094 // define errors for RAA vs. centrality
2095 InitDstartoD0pi2010PbPb2040();
2096 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2097 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
2099 //--------------------------------------------------------------------------
2100 void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
2101 // define errors for RAA vs. centrality
2102 InitDstartoD0pi2010PbPb4080();
2103 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2104 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
2106 //--------------------------------------------------------------------------
2107 void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
2108 // define errors for RAA vs. centrality
2109 InitDstartoD0pi2010PbPb4080();
2110 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
2111 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
2114 //--------------------------------------------------------------------------
2115 void AliHFSystErr::InitD0toKpi2011PbPb010CentScan(){
2116 // define errors for RAA vs. centrality 2011
2117 InitD0toKpi2011PbPb07half();
2120 //--------------------------------------------------------------------------
2121 void AliHFSystErr::InitD0toKpi2011PbPb1020CentScan(){
2122 // define errors for RAA vs. centrality 2011
2123 InitD0toKpi2011PbPb07half();
2126 //--------------------------------------------------------------------------
2127 void AliHFSystErr::InitD0toKpi2011PbPb2030CentScan(){
2128 // define errors for RAA vs. centrality 2011
2129 InitD0toKpi2011PbPb07half();
2132 //--------------------------------------------------------------------------
2133 void AliHFSystErr::InitD0toKpi2011PbPb3040CentScan(){
2134 // define errors for RAA vs. centrality 2011
2135 InitD0toKpi2011PbPb07half();
2138 //--------------------------------------------------------------------------
2139 void AliHFSystErr::InitD0toKpi2011PbPb4050CentScan(){
2140 // define errors for RAA vs. centrality 2011
2141 InitD0toKpi2011PbPb07half();
2144 //--------------------------------------------------------------------------
2145 void AliHFSystErr::InitD0toKpi2010PbPb5080CentScan(){
2146 // define errors for RAA vs. centrality 2011
2147 InitD0toKpi2011PbPb07half();
2151 //--------------------------------------------------------------------------
2152 void AliHFSystErr::InitDplustoKpipi2011PbPb010CentScan(){
2153 // define errors for RAA vs. centrality 2011
2154 InitDplustoKpipi2011PbPb07half();
2156 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.15);
2157 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2158 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.15);
2160 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2161 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2162 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2164 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
2165 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2166 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.10);
2169 //--------------------------------------------------------------------------
2170 void AliHFSystErr::InitDplustoKpipi2011PbPb1020CentScan(){
2171 // define errors for RAA vs. centrality 2011
2172 InitDplustoKpipi2011PbPb07half();
2174 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2175 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2176 for(Int_t i=9;i<=17;i++)fRawYield->SetBinContent(i,0.15);
2178 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2179 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2180 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2182 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
2183 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2184 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.15);
2187 //--------------------------------------------------------------------------
2188 void AliHFSystErr::InitDplustoKpipi2011PbPb2030CentScan(){
2189 // define errors for RAA vs. centrality 2011
2190 InitDplustoKpipi2011PbPb07half();
2192 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2193 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2194 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.08);
2196 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2197 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2198 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2200 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2201 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2202 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2205 //--------------------------------------------------------------------------
2206 void AliHFSystErr::InitDplustoKpipi2011PbPb3040CentScan(){
2207 // define errors for RAA vs. centrality 2011
2208 InitDplustoKpipi2011PbPb07half();
2210 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.20);
2211 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2212 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.20);
2214 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2215 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2216 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2218 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2219 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
2220 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2223 //--------------------------------------------------------------------------
2224 void AliHFSystErr::InitDplustoKpipi2011PbPb4050CentScan(){
2225 // define errors for RAA vs. centrality 2011
2226 InitDplustoKpipi2011PbPb07half();
2228 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
2229 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
2230 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.05);
2232 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2233 for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2234 for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2236 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
2237 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
2238 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
2241 //--------------------------------------------------------------------------
2242 void AliHFSystErr::InitDplustoKpipi2010PbPb5080CentScan(){
2243 // define errors for RAA vs. centrality 2011
2244 InitDplustoKpipi2010PbPb4080();
2246 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2247 for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2248 for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.10);
2250 for(Int_t i=1;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
2251 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20); // 10%
2252 for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10); // 10%
2253 for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20); // 10%
2256 //--------------------------------------------------------------------------
2257 void AliHFSystErr::InitDstartoD0pi2011PbPb010CentScan(){
2258 // define errors for RAA vs. centrality 2011
2259 InitDstartoD0pi2011PbPb07half();
2262 //--------------------------------------------------------------------------
2263 void AliHFSystErr::InitDstartoD0pi2011PbPb1020CentScan(){
2264 // define errors for RAA vs. centrality 2011
2265 InitDstartoD0pi2011PbPb07half();
2268 //--------------------------------------------------------------------------
2269 void AliHFSystErr::InitDstartoD0pi2011PbPb2030CentScan(){
2270 // define errors for RAA vs. centrality 2011
2271 InitDstartoD0pi2011PbPb07half();
2274 //--------------------------------------------------------------------------
2275 void AliHFSystErr::InitDstartoD0pi2011PbPb3040CentScan(){
2276 // define errors for RAA vs. centrality 2011
2277 InitDstartoD0pi2011PbPb07half();
2280 //--------------------------------------------------------------------------
2281 void AliHFSystErr::InitDstartoD0pi2011PbPb4050CentScan(){
2282 // define errors for RAA vs. centrality 2011
2283 InitDstartoD0pi2011PbPb07half();
2286 //--------------------------------------------------------------------------
2287 void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
2288 // define errors for RAA vs. centrality 2011
2289 InitDstartoD0pi2011PbPb07half();
2293 //_________________________________________________________________________
2294 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
2295 InitD0toKpi2013pPb0100();
2297 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2299 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2300 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.06);
2301 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2303 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2305 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2306 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2307 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2309 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2310 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2311 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2315 //_________________________________________________________________________
2316 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
2317 InitD0toKpi2013pPb0100();
2319 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2321 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2322 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.08);
2323 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2325 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2327 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2328 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2329 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2331 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2332 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2333 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2335 //_________________________________________________________________________
2336 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
2337 InitD0toKpi2013pPb0100();
2339 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2341 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2342 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.12);
2343 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2345 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2347 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2348 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2349 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2351 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2352 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2353 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2356 //_________________________________________________________________________
2357 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
2358 InitD0toKpi2013pPb0100();
2360 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2362 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2363 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2364 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2366 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2368 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2369 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2370 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2372 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2373 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2374 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2377 //_________________________________________________________________________
2378 void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
2379 InitD0toKpi2013pPb0100();
2381 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2383 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2384 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2385 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2387 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2389 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2390 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2391 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2393 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2394 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2395 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2399 //_________________________________________________________________________
2400 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
2401 InitDplustoKpipi2013pPb0100();
2403 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2405 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2406 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2407 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2409 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2411 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2412 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2413 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2415 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2416 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.10);
2417 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.10);
2420 //_________________________________________________________________________
2421 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
2422 InitDplustoKpipi2013pPb0100();
2423 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2425 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2426 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2427 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2429 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2431 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2432 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2433 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2435 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2436 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2437 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2440 //_________________________________________________________________________
2441 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
2442 InitDplustoKpipi2013pPb0100();
2443 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2445 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2446 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2447 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2449 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2451 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2452 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2453 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2455 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2456 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2457 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2460 //_________________________________________________________________________
2461 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
2462 InitDplustoKpipi2013pPb0100();
2463 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2465 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2466 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2467 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2469 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2471 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2472 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2473 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2475 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2476 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2477 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2480 //_________________________________________________________________________
2481 void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
2482 InitDplustoKpipi2013pPb0100();
2484 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2486 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
2487 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2488 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2490 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2492 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
2493 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
2494 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
2496 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2497 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2498 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2503 //_________________________________________________________________________
2504 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
2505 InitDstartoD0pi2013pPb0100();
2507 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2509 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2510 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2511 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2513 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2515 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2516 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2517 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2519 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2520 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2521 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2524 //_________________________________________________________________________
2525 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
2526 InitDstartoD0pi2013pPb0100();
2528 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2530 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2531 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
2532 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2534 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2536 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2537 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2538 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2540 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2541 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2542 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2545 //_________________________________________________________________________
2546 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
2547 InitDstartoD0pi2013pPb0100();
2549 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2551 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
2552 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2553 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2555 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2557 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2558 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2559 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2561 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2562 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2563 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2566 //_________________________________________________________________________
2567 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
2568 InitDstartoD0pi2013pPb0100();
2570 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2572 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
2573 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
2574 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2576 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2578 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2579 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2580 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2582 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2583 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2584 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2586 //_________________________________________________________________________
2587 void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
2588 InitDstartoD0pi2013pPb0100();
2590 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2592 for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
2593 for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.07);
2594 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2596 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2598 for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
2599 for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
2600 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
2602 for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
2603 for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
2604 for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
2608 //_________________________________________________________________________
2609 void AliHFSystErr::InitD0toKpi2013pPb020V0A(){
2610 InitD0toKpi2013pPb0100();
2612 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2615 // Raw yield extraction
2616 fRawYield->SetBinContent(1,0.0);
2617 fRawYield->SetBinContent(2,0.085);
2618 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2619 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2620 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2621 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2622 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2623 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2625 // Cuts efficiency (from cuts variation)
2626 fCutsEff->SetBinContent(1,0.0);
2627 fCutsEff->SetBinContent(2,0.10);
2628 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2630 // PID efficiency (from PID/noPID)
2631 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2636 //_________________________________________________________________________
2637 void AliHFSystErr::InitD0toKpi2013pPb2040V0A(){
2638 InitD0toKpi2013pPb0100();
2640 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2643 // Raw yield extraction
2644 fRawYield->SetBinContent(1,0.0);
2645 fRawYield->SetBinContent(2,0.085);
2646 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2647 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2648 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2649 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2650 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2651 //for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2654 // Cuts efficiency (from cuts variation)
2655 fCutsEff->SetBinContent(1,0.0);
2656 fCutsEff->SetBinContent(2,0.10);
2657 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2660 // PID efficiency (from PID/noPID)
2661 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2665 //_________________________________________________________________________
2666 void AliHFSystErr::InitD0toKpi2013pPb4060V0A(){
2667 InitD0toKpi2013pPb0100();
2669 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2672 // Raw yield extraction
2673 fRawYield->SetBinContent(1,0.0);
2674 fRawYield->SetBinContent(2,0.085);
2675 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2676 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2677 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2678 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2679 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2680 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2682 // Cuts efficiency (from cuts variation)
2683 fCutsEff->SetBinContent(1,0.0);
2684 fCutsEff->SetBinContent(2,0.10);
2685 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2688 // PID efficiency (from PID/noPID)
2689 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2693 //_________________________________________________________________________
2694 void AliHFSystErr::InitD0toKpi2013pPb60100V0A(){
2695 InitD0toKpi2013pPb0100();
2697 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2699 // Raw yield extraction
2700 fRawYield->SetBinContent(1,0.0);
2701 fRawYield->SetBinContent(2,0.085);
2702 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2703 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2704 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2705 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2706 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2709 // Cuts efficiency (from cuts variation)
2710 fCutsEff->SetBinContent(1,0.0);
2711 fCutsEff->SetBinContent(2,0.10);
2712 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2715 // PID efficiency (from PID/noPID)
2716 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
2720 //_________________________________________________________________________
2721 void AliHFSystErr::InitD0toKpi2013pPb020ZNA(){
2722 InitD0toKpi2013pPb0100();
2724 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
2726 // Raw yield extraction
2727 fRawYield->SetBinContent(1,0.0);
2728 fRawYield->SetBinContent(2,0.085);
2729 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2730 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2731 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2732 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2733 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2734 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2736 // Cuts efficiency (from cuts variation)
2737 fCutsEff->SetBinContent(1,0.0);
2738 fCutsEff->SetBinContent(2,0.10);
2739 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2742 // PID efficiency (from PID/noPID)
2743 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2747 //_________________________________________________________________________
2748 void AliHFSystErr::InitD0toKpi2013pPb2040ZNA(){
2749 InitD0toKpi2013pPb0100();
2751 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
2753 // Raw yield extraction
2754 fRawYield->SetBinContent(1,0.0);
2755 fRawYield->SetBinContent(2,0.085);
2756 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2757 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2758 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2759 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2760 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2761 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2764 // Cuts efficiency (from cuts variation)
2765 fCutsEff->SetBinContent(1,0.0);
2766 fCutsEff->SetBinContent(2,0.10);
2767 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2769 // PID efficiency (from PID/noPID)
2770 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2774 //_________________________________________________________________________
2775 void AliHFSystErr::InitD0toKpi2013pPb4060ZNA(){
2776 InitD0toKpi2013pPb0100();
2778 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
2780 // Raw yield extraction
2781 fRawYield->SetBinContent(1,0.0);
2782 fRawYield->SetBinContent(2,0.085);
2783 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2784 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2785 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2786 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2787 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2788 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2790 // Cuts efficiency (from cuts variation)
2791 fCutsEff->SetBinContent(1,0.0);
2792 fCutsEff->SetBinContent(2,0.10);
2793 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2795 // PID efficiency (from PID/noPID)
2796 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2800 //_________________________________________________________________________
2801 void AliHFSystErr::InitD0toKpi2013pPb60100ZNA(){
2802 InitD0toKpi2013pPb0100();
2804 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
2806 // Raw yield extraction
2807 fRawYield->SetBinContent(1,0.0);
2808 fRawYield->SetBinContent(2,0.085);
2809 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
2810 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
2811 for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
2812 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
2813 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
2814 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2816 // Cuts efficiency (from cuts variation)
2817 fCutsEff->SetBinContent(1,0.0);
2818 fCutsEff->SetBinContent(2,0.10);
2819 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2821 // PID efficiency (from PID/noPID)
2822 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
2826 //_________________________________________________________________________
2827 void AliHFSystErr::InitDstartoD0pi2013pPb020V0A(){
2828 InitDstartoD0pi2013pPb0100();
2830 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2832 // Raw yield extraction
2833 fRawYield->SetBinContent(1,0.0);
2834 fRawYield->SetBinContent(2,0.10);
2835 fRawYield->SetBinContent(3,0.06);
2836 fRawYield->SetBinContent(4,0.06);
2837 fRawYield->SetBinContent(5,0.05);
2838 fRawYield->SetBinContent(6,0.05);
2839 fRawYield->SetBinContent(7,0.05);
2840 fRawYield->SetBinContent(8,0.05);
2841 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2842 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2843 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2845 // Cuts efficiency (from cuts variation)
2846 fCutsEff->SetBinContent(1,0.0);
2847 fCutsEff->SetBinContent(2,0.10);
2848 fCutsEff->SetBinContent(3,0.06);
2849 fCutsEff->SetBinContent(4,0.06);
2850 fCutsEff->SetBinContent(5,0.06);
2851 fCutsEff->SetBinContent(6,0.06);
2852 fCutsEff->SetBinContent(7,0.06);
2853 fCutsEff->SetBinContent(8,0.06);
2854 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
2855 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2857 // PID efficiency (from PID/noPID)
2858 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2863 //_________________________________________________________________________
2864 void AliHFSystErr::InitDstartoD0pi2013pPb2040V0A(){
2865 InitDstartoD0pi2013pPb0100();
2867 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2869 // Raw yield extraction
2870 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2871 fRawYield->SetBinContent(1,0.0);
2872 fRawYield->SetBinContent(2,0.10);
2873 fRawYield->SetBinContent(3,0.06);
2874 fRawYield->SetBinContent(4,0.06);
2875 fRawYield->SetBinContent(5,0.05);
2876 fRawYield->SetBinContent(6,0.05);
2877 fRawYield->SetBinContent(7,0.05);
2878 fRawYield->SetBinContent(8,0.05);
2879 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2880 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2881 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
2883 // Cuts efficiency (from cuts variation)
2884 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2885 fCutsEff->SetBinContent(1,0.);
2886 fCutsEff->SetBinContent(2,0.10);
2887 fCutsEff->SetBinContent(3,0.06);
2888 fCutsEff->SetBinContent(4,0.06);
2889 fCutsEff->SetBinContent(5,0.06);
2890 fCutsEff->SetBinContent(6,0.06);
2891 fCutsEff->SetBinContent(7,0.06);
2892 fCutsEff->SetBinContent(8,0.06);
2893 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2894 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
2895 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2897 // PID efficiency (from PID/noPID)
2898 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2899 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2902 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2903 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2904 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2905 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2907 // particle-antiparticle
2909 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2910 fPartAntipart->SetBinContent(1,1);
2911 fPartAntipart->SetBinContent(2,1);
2912 fPartAntipart->SetBinContent(3,0.12);
2913 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2918 //_________________________________________________________________________
2919 void AliHFSystErr::InitDstartoD0pi2013pPb4060V0A(){
2920 InitDstartoD0pi2013pPb0100();
2922 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2924 // Raw yield extraction
2925 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2926 fRawYield->SetBinContent(1,0.0);
2927 fRawYield->SetBinContent(2,0.0);
2928 fRawYield->SetBinContent(3,0.06);
2929 fRawYield->SetBinContent(4,0.06);
2930 fRawYield->SetBinContent(5,0.05);
2931 fRawYield->SetBinContent(6,0.05);
2932 fRawYield->SetBinContent(7,0.05);
2933 fRawYield->SetBinContent(8,0.05);
2934 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2935 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2936 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2938 // Cuts efficiency (from cuts variation)
2939 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2940 fCutsEff->SetBinContent(1,0.);
2941 fCutsEff->SetBinContent(2,0.);
2942 fCutsEff->SetBinContent(3,0.06);
2943 fCutsEff->SetBinContent(4,0.06);
2944 fCutsEff->SetBinContent(5,0.06);
2945 fCutsEff->SetBinContent(6,0.06);
2946 fCutsEff->SetBinContent(7,0.06);
2947 fCutsEff->SetBinContent(8,0.06);
2948 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
2949 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
2950 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2952 // PID efficiency (from PID/noPID)
2953 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
2954 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2957 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
2958 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
2959 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
2960 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
2962 // particle-antiparticle
2964 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
2965 fPartAntipart->SetBinContent(1,1);
2966 fPartAntipart->SetBinContent(2,1);
2967 fPartAntipart->SetBinContent(3,0.12);
2968 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
2974 //_________________________________________________________________________
2975 void AliHFSystErr::InitDstartoD0pi2013pPb60100V0A(){
2976 InitDstartoD0pi2013pPb0100();
2978 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2980 // Raw yield extraction
2981 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
2982 fRawYield->SetBinContent(1,0.0);
2983 fRawYield->SetBinContent(2,0.0);
2984 fRawYield->SetBinContent(3,0.06);
2985 fRawYield->SetBinContent(4,0.06);
2986 fRawYield->SetBinContent(5,0.05);
2987 fRawYield->SetBinContent(6,0.05);
2988 fRawYield->SetBinContent(7,0.05);
2989 fRawYield->SetBinContent(8,0.05);
2990 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
2991 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2992 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2994 // Cuts efficiency (from cuts variation)
2995 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
2996 fCutsEff->SetBinContent(1,0.);
2997 fCutsEff->SetBinContent(2,0.);
2998 fCutsEff->SetBinContent(3,0.06);
2999 fCutsEff->SetBinContent(4,0.06);
3000 fCutsEff->SetBinContent(5,0.06);
3001 fCutsEff->SetBinContent(6,0.06);
3002 fCutsEff->SetBinContent(7,0.06);
3003 fCutsEff->SetBinContent(8,0.06);
3004 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3005 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3006 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3008 // PID efficiency (from PID/noPID)
3009 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3010 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3013 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3014 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3015 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3016 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3018 // particle-antiparticle
3020 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3021 fPartAntipart->SetBinContent(1,1);
3022 fPartAntipart->SetBinContent(2,1);
3023 fPartAntipart->SetBinContent(3,0.12);
3024 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3030 //_________________________________________________________________________
3031 void AliHFSystErr::InitDstartoD0pi2013pPb020ZNA(){
3032 InitDstartoD0pi2013pPb0100();
3034 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3036 // Raw yield extraction
3037 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3038 fRawYield->SetBinContent(1,0.);
3039 fRawYield->SetBinContent(2,0.10);
3040 fRawYield->SetBinContent(3,0.06);
3041 fRawYield->SetBinContent(4,0.06);
3042 fRawYield->SetBinContent(5,0.05);
3043 fRawYield->SetBinContent(6,0.05);
3044 fRawYield->SetBinContent(7,0.05);
3045 fRawYield->SetBinContent(8,0.05);
3046 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3047 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3048 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3050 // Cuts efficiency (from cuts variation)
3051 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3052 fCutsEff->SetBinContent(1,0.);
3053 fCutsEff->SetBinContent(2,0.10);
3054 fCutsEff->SetBinContent(3,0.06);
3055 fCutsEff->SetBinContent(4,0.06);
3056 fCutsEff->SetBinContent(5,0.06);
3057 fCutsEff->SetBinContent(6,0.06);
3058 fCutsEff->SetBinContent(7,0.06);
3059 fCutsEff->SetBinContent(8,0.06);
3060 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3061 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3062 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3064 // PID efficiency (from PID/noPID)
3065 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3066 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
3069 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3070 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3071 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3072 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3074 // particle-antiparticle
3076 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3077 fPartAntipart->SetBinContent(1,1);
3078 fPartAntipart->SetBinContent(2,1);
3079 fPartAntipart->SetBinContent(3,0.12);
3080 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3086 //_________________________________________________________________________
3087 void AliHFSystErr::InitDstartoD0pi2013pPb2040ZNA(){
3088 InitDstartoD0pi2013pPb0100();
3090 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3092 // Raw yield extraction
3093 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3094 fRawYield->SetBinContent(1,0.0);
3095 fRawYield->SetBinContent(2,0.0);
3096 fRawYield->SetBinContent(3,0.06);
3097 fRawYield->SetBinContent(4,0.06);
3098 fRawYield->SetBinContent(5,0.05);
3099 fRawYield->SetBinContent(6,0.05);
3100 fRawYield->SetBinContent(7,0.05);
3101 fRawYield->SetBinContent(8,0.05);
3102 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3103 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3104 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3106 // Cuts efficiency (from cuts variation)
3107 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3108 fCutsEff->SetBinContent(1,0.);
3109 fCutsEff->SetBinContent(2,0.);
3110 fCutsEff->SetBinContent(3,0.06);
3111 fCutsEff->SetBinContent(4,0.06);
3112 fCutsEff->SetBinContent(5,0.06);
3113 fCutsEff->SetBinContent(6,0.06);
3114 fCutsEff->SetBinContent(7,0.06);
3115 fCutsEff->SetBinContent(8,0.06);
3116 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3117 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3118 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3120 // PID efficiency (from PID/noPID)
3121 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3122 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3125 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3126 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3127 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3128 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3130 // particle-antiparticle
3132 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3133 fPartAntipart->SetBinContent(1,1);
3134 fPartAntipart->SetBinContent(2,1);
3135 fPartAntipart->SetBinContent(3,0.12);
3136 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3142 //_________________________________________________________________________
3143 void AliHFSystErr::InitDstartoD0pi2013pPb4060ZNA(){
3144 InitDstartoD0pi2013pPb0100();
3146 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3148 // Raw yield extraction
3149 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3150 fRawYield->SetBinContent(1,0.0);
3151 fRawYield->SetBinContent(2,0.0);
3152 fRawYield->SetBinContent(3,0.06);
3153 fRawYield->SetBinContent(4,0.06);
3154 fRawYield->SetBinContent(5,0.05);
3155 fRawYield->SetBinContent(6,0.05);
3156 fRawYield->SetBinContent(7,0.05);
3157 fRawYield->SetBinContent(8,0.05);
3158 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3159 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3160 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3162 // Cuts efficiency (from cuts variation)
3163 // fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3164 fCutsEff->SetBinContent(1,0.);
3165 fCutsEff->SetBinContent(2,0.);
3166 fCutsEff->SetBinContent(3,0.06);
3167 fCutsEff->SetBinContent(4,0.06);
3168 fCutsEff->SetBinContent(5,0.06);
3169 fCutsEff->SetBinContent(6,0.06);
3170 fCutsEff->SetBinContent(7,0.06);
3171 fCutsEff->SetBinContent(8,0.06);
3172 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3173 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3174 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3176 // PID efficiency (from PID/noPID)
3177 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3178 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3181 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3182 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3183 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3184 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3186 // particle-antiparticle
3188 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3189 fPartAntipart->SetBinContent(1,1);
3190 fPartAntipart->SetBinContent(2,1);
3191 fPartAntipart->SetBinContent(3,0.12);
3192 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3199 //_________________________________________________________________________
3200 void AliHFSystErr::InitDstartoD0pi2013pPb60100ZNA(){
3201 InitDstartoD0pi2013pPb0100();
3203 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3205 // Raw yield extraction
3206 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3207 fRawYield->SetBinContent(1,0.);
3208 fRawYield->SetBinContent(2,0.);
3209 fRawYield->SetBinContent(3,0.06);
3210 fRawYield->SetBinContent(4,0.06);
3211 fRawYield->SetBinContent(5,0.05);
3212 fRawYield->SetBinContent(6,0.05);
3213 fRawYield->SetBinContent(7,0.05);
3214 fRawYield->SetBinContent(8,0.05);
3215 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
3216 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
3217 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
3219 // Cuts efficiency (from cuts variation)
3220 //fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3221 fCutsEff->SetBinContent(1,0.);
3222 fCutsEff->SetBinContent(2,0.);
3223 fCutsEff->SetBinContent(3,0.06);
3224 fCutsEff->SetBinContent(4,0.06);
3225 fCutsEff->SetBinContent(5,0.06);
3226 fCutsEff->SetBinContent(6,0.06);
3227 fCutsEff->SetBinContent(7,0.06);
3228 fCutsEff->SetBinContent(8,0.06);
3229 for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
3230 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
3231 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3233 // PID efficiency (from PID/noPID)
3234 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3235 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
3238 /* fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3239 for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
3240 for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
3241 for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);*/
3243 // particle-antiparticle
3245 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
3246 fPartAntipart->SetBinContent(1,1);
3247 fPartAntipart->SetBinContent(2,1);
3248 fPartAntipart->SetBinContent(3,0.12);
3249 for(Int_t i=4;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
3257 //--------------------------------------------------------------------------
3258 void AliHFSystErr::InitLctopKpi2010pp() {
3260 // Lc->pKpi syst errors. Responsible: R. Romita
3265 fNorm = new TH1F("fNorm","fNorm",4,2,6);
3266 for(Int_t i=1;i<=4;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
3268 // Tracking efficiency
3269 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",4,2,6);
3270 for(Int_t i=1;i<=4;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
3272 // Raw yield extraction
3273 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3274 fRawYield->SetBinContent(1,0.20);
3275 fRawYield->SetBinContent(2,0.15);
3276 fRawYield->SetBinContent(3,0.10);
3277 fRawYield->SetBinContent(4,0.15);
3279 fCutsEff = new TH1F("fCutsEff","fCutsEff",4,2,6);
3280 fCutsEff->SetBinContent(1,0.50);
3281 fCutsEff->SetBinContent(2,0.07);
3282 fCutsEff->SetBinContent(3,0.30);
3283 fCutsEff->SetBinContent(4,0.40);
3285 // PID efficiency (from PID/noPID)
3286 fPIDEff = new TH1F("fPIDEff","fPIDEff",4,2,4);
3287 fPIDEff->SetBinContent(1,0.13); // 15%
3288 fPIDEff->SetBinContent(2,0.30); // 15%
3289 fPIDEff->SetBinContent(3,0.5); // 15%
3290 fPIDEff->SetBinContent(4,0.30); // 15%
3293 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
3294 for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
3295 for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
3296 for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
3297 for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
3298 for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
3301 //--------------------------------------------------------------------------
3302 void AliHFSystErr::InitLctopK0S2010pp() {
3304 // Lc->pK0S syst errors. Responsible: A. De Caro
3306 // Only tool preparation: no physical mean for these values!
3310 fNorm = new TH1F("fNorm","fNorm",12,0.,12.);
3311 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05); // 5% error on sigmaV0and // only a try
3313 // Tracking efficiency
3314 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0.,12.);
3315 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.1); // 10% (4% per track) // only a try
3317 // Raw yield extraction
3318 fRawYield = new TH1F("fRawYield","fRawYield",12,0.,12.);
3319 for(Int_t i=1;i<=2;i++) fRawYield->SetBinContent(i,1.); // only a try
3320 for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.1); // only a try
3322 fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0.,12.);
3323 for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.1); // only a try
3325 // PID efficiency (from PID/noPID)
3326 fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0.,12.);
3327 for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // only a try
3330 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0.,12.);
3331 for(Int_t i=1; i<=2;i++) fMCPtShape->SetBinContent(i,1.); // only a try
3332 for(Int_t i=3; i<=12;i++) fMCPtShape->SetBinContent(i,0.1); // only a try
3336 //--------------------------------------------------------------------------
3337 Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
3342 Int_t bin=fCutsEff->FindBin(pt);
3344 return fCutsEff->GetBinContent(bin);
3346 //--------------------------------------------------------------------------
3347 Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
3352 Int_t bin=fMCPtShape->FindBin(pt);
3354 return fMCPtShape->GetBinContent(bin);
3356 //--------------------------------------------------------------------------
3357 Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
3362 Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3364 return TMath::Sqrt(err);
3366 //--------------------------------------------------------------------------
3367 Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
3372 Int_t bin=fPIDEff->FindBin(pt);
3374 return fPIDEff->GetBinContent(bin);
3376 //--------------------------------------------------------------------------
3377 Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
3382 Int_t bin=fTrackingEff->FindBin(pt);
3384 return fTrackingEff->GetBinContent(bin);
3386 //--------------------------------------------------------------------------
3387 Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
3392 Int_t bin=fRawYield->FindBin(pt);
3394 return fRawYield->GetBinContent(bin);
3396 //--------------------------------------------------------------------------
3397 Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
3402 Int_t bin=fPartAntipart->FindBin(pt);
3404 return fPartAntipart->GetBinContent(bin);
3406 //--------------------------------------------------------------------------
3407 Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
3409 // Get total syst error (except norm. error)
3414 if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
3415 if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
3416 // if(fBR) err += GetBRErr()*GetBRErr();
3417 if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
3418 if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
3419 if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3420 if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
3422 err += feeddownErr*feeddownErr;
3424 return TMath::Sqrt(err);
3426 //---------------------------------------------------------------------------
3427 void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
3431 gStyle->SetOptStat(0);
3433 TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
3434 cSystErr->Range(0.20,-0.5,18.4,0.34);
3435 cSystErr->SetRightMargin(0.318);
3436 cSystErr->SetFillColor(0);
3438 TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
3439 hFrame->SetAxisRange(1.,35.9,"X");
3440 hFrame->SetAxisRange(-0.5,0.5,"Y");
3443 TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
3444 leg->SetTextSize(0.03601695);
3445 leg->SetFillStyle(0);
3446 leg->SetBorderSize(0);
3448 TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
3449 Int_t nbins = fNorm->GetNbinsX();
3450 TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
3451 for(Int_t i=1;i<=36;i++) {
3452 Double_t pt = hTotErr->GetBinCenter(i);
3453 Double_t ptwidth = hTotErr->GetBinWidth(i);
3456 Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
3457 Double_t toterryl=0., toterryh=0.;
3458 for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
3459 grErrFeeddown->GetPoint(j,x,y);
3460 errxh = grErrFeeddown->GetErrorXhigh(j);
3461 errxl = grErrFeeddown->GetErrorXlow(j);
3462 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
3463 erryh = grErrFeeddown->GetErrorYhigh(j);
3464 erryl = grErrFeeddown->GetErrorYlow(j);
3467 if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
3468 else toterryl = GetTotalSystErr(pt);
3469 if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
3470 else toterryh = GetTotalSystErr(pt);
3472 hTotErr->SetBinContent(i,toterryh);
3473 gTotErr->SetPoint(i,pt,0.);
3474 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh); // i, exl, exh, eyl, eyh
3477 hTotErr->SetBinContent(i,GetTotalSystErr(pt));
3478 gTotErr->SetPoint(i,pt,0.);
3479 gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt)); // i, exl, exh, eyl, eyh
3483 gTotErr->SetLineColor(kBlack);
3484 gTotErr->SetFillColor(kRed);
3485 gTotErr->SetFillStyle(3002);
3487 leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
3488 // hTotErr->SetLineColor(1);
3489 // hTotErr->SetLineWidth(3);
3490 // hTotErr->Draw("same");
3491 // leg->AddEntry(hTotErr,"Total (excl. norm.)","l");
3494 fNorm->SetFillColor(1);
3495 fNorm->SetFillStyle(3002);
3496 //fNorm->Draw("same");
3497 //TH1F *hNormRefl = ReflectHisto(fNorm);
3498 //hNormRefl->Draw("same");
3499 Double_t norm = fNorm->GetBinContent(1)*100;
3500 leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
3503 grErrFeeddown->SetFillColor(kTeal-8);
3504 grErrFeeddown->SetFillStyle(3001);
3505 grErrFeeddown->Draw("2");
3506 leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
3509 fTrackingEff->SetFillColor(4);
3510 fTrackingEff->SetFillStyle(3006);
3511 fTrackingEff->Draw("same");
3512 TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
3513 hTrackingEffRefl->Draw("same");
3514 leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
3517 fBR->SetFillColor(6);
3518 fBR->SetFillStyle(3005);
3519 //fBR->SetFillStyle(3020);
3521 TH1F *hBRRefl = ReflectHisto(fBR);
3522 hBRRefl->Draw("same");
3523 leg->AddEntry(fBR,"Branching ratio","f");
3526 Int_t ci; // for color index setting
3527 ci = TColor::GetColor("#00cc00");
3528 fRawYield->SetLineColor(ci);
3529 // fRawYield->SetLineColor(3);
3530 fRawYield->SetLineWidth(3);
3531 fRawYield->Draw("same");
3532 TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
3533 hRawYieldRefl->Draw("same");
3534 leg->AddEntry(fRawYield,"Yield extraction","l");
3537 fCutsEff->SetLineColor(4);
3538 fCutsEff->SetLineWidth(3);
3539 fCutsEff->Draw("same");
3540 TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
3541 hCutsEffRefl->Draw("same");
3542 leg->AddEntry(fCutsEff,"Cuts efficiency","l");
3545 fPIDEff->SetLineColor(7);
3546 fPIDEff->SetLineWidth(3);
3547 fPIDEff->Draw("same");
3548 TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
3549 hPIDEffRefl->Draw("same");
3550 leg->AddEntry(fPIDEff,"PID efficiency","l");
3553 Int_t ci = TColor::GetColor("#9933ff");
3554 fMCPtShape->SetLineColor(ci);
3555 // fMCPtShape->SetLineColor(8);
3556 fMCPtShape->SetLineWidth(3);
3557 fMCPtShape->Draw("same");
3558 TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
3559 hMCPtShapeRefl->Draw("same");
3560 leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
3563 Int_t ci = TColor::GetColor("#ff6600");
3564 fPartAntipart->SetLineColor(ci);
3565 // fPartAntipart->SetLineColor(9);
3566 fPartAntipart->SetLineWidth(3);
3567 fPartAntipart->Draw("same");
3568 TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
3569 hPartAntipartRefl->Draw("same");
3570 leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
3576 cSystErr->SaveAs("RelativeSystematics.eps");
3580 //-------------------------------------------------------------------------
3581 TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
3583 // Clones and reflects histogram
3585 TH1F *hout=(TH1F*)hin->Clone("hout");