]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGHF/vertexingHF/AliHFSystErr.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Add systematic uncertainty values for Raa vs pt in the 0-10% class for 2011 data...
[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / AliHFSystErr.cxx
CommitLineData
ada850b0 1
0233abe6 2/**************************************************************************
3 * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
4 * *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
7 * *
8 * Permission to use, copy, modify and distribute this software and its *
9 * documentation strictly for non-commercial purposes is hereby granted *
10 * without fee, provided that the above copyright notice appears in all *
11 * copies and that both the copyright notice and this permission notice *
12 * appear in the supporting documentation. The authors make no claims *
13 * about the suitability of this software for any purpose. It is *
14 * provided "as is" without express or implied warranty. *
15 **************************************************************************/
16
27de2dfb 17/* $Id$ */
18
0233abe6 19/////////////////////////////////////////////////////////////
20//
21// Class to handle systematic errors for charm hadrons
22//
23// Usage:
53c91ed7 24// AliHFSystEff syst; // DECAY = 1 for D0, 2, for D+, 3 for D*, 5 for Lc->pKp, 6 for Lc->pK0S
5541b811 25// syst.SetRunNumber(YEAR); // YEAR = two last numbers of the year (is 10 for 2010)
26// syst.SetCollisionType(TYPE); // TYPE = 0 is pp, 1 is PbPb
27// syst.SetCentrality(CENT); // CENT is centrality, 0100 for MB, 020 (4080) for 0-20 (40-80) CC...
28// syst.Init(DECAY); // DECAY = 1 for D0, 2, for D+, 3 for D*
0233abe6 29// syst.DrawErrors(); // to see a plot of the error contributions
30// syst.GetTotalSystErr(pt); // to get the total err at pt
31//
32// Author: A.Dainese, andrea.dainese@pd.infn.it
33/////////////////////////////////////////////////////////////
34
35#include <TStyle.h>
36#include <TGraphAsymmErrors.h>
37#include <TMath.h>
38#include <TCanvas.h>
39#include <TH2F.h>
40#include <TLegend.h>
5ad8ad40 41#include <TColor.h>
0233abe6 42
5541b811 43#include "AliLog.h"
0233abe6 44#include "AliHFSystErr.h"
45
46
47ClassImp(AliHFSystErr)
48
49//--------------------------------------------------------------------------
50AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
51TNamed(name,title),
52fNorm(0),
53fRawYield(0),
54fTrackingEff(0),
55fBR(0),
56fCutsEff(0),
57fPIDEff(0),
58fMCPtShape(0),
5541b811 59fPartAntipart(0),
60fRunNumber(10),
61fCollisionType(0),
a8c3295d 62fCentralityClass("0100"),
44dd56b2 63fRapidityRange("0101"),
cede55f5 64fIsLowEnergy(false),
44dd56b2 65fIsCentScan(false),
66fIsRapidityScan(false)
0233abe6 67{
68 //
69 // Default Constructor
70 //
71}
5541b811 72
0233abe6 73//--------------------------------------------------------------------------
5541b811 74AliHFSystErr::~AliHFSystErr() {
75 //
76 // Default Destructor
0233abe6 77 //
a8c3295d 78 /*
5541b811 79
80 if(fNorm) { delete fNorm; fNorm=0; }
81 if(fRawYield) { delete fRawYield; fRawYield=0; }
82 if(fTrackingEff) { delete fTrackingEff; fTrackingEff=0; }
83 if(fBR) { delete fBR; fBR=0; }
84 if(fCutsEff) { delete fCutsEff; fCutsEff=0; }
85 if(fPIDEff) { delete fPIDEff; fPIDEff=0; }
86 if(fMCPtShape) { delete fMCPtShape; fMCPtShape=0; }
87 if(fPartAntipart) { delete fPartAntipart; fPartAntipart=0; }
a8c3295d 88 */
5541b811 89}
90
91//--------------------------------------------------------------------------
92void AliHFSystErr::Init(Int_t decay){
93 //
94 // Variables/histos initialization
0233abe6 95 //
96
9341a8fa 97 if ((fRunNumber>11) && fIsLowEnergy==false) {
27360a0e 98 AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
5541b811 99 }
100
0233abe6 101 switch(decay) {
102 case 1: // D0->Kpi
a8c3295d 103 if (fCollisionType==0) {
104 if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
105 else InitD0toKpi2010pp();
9341a8fa 106 }
107 else if (fCollisionType==1) {
3db225bf 108 if (fRunNumber == 10){
109 if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
110 else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
111 else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
112 else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
113 else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
114 else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
115 else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
116 else AliFatal("Not yet implemented");
117 }
9341a8fa 118 else if (fRunNumber == 11){
3db225bf 119 if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
c8617efe 120 else if (fCentralityClass=="010") InitD0toKpi2011PbPb010();
3db225bf 121 else if (fCentralityClass=="3050InPlane") InitD0toKpi2011PbPb3050InPlane();
830baf70 122 else if (fCentralityClass=="3050OutOfPlane")InitD0toKpi2011PbPb3050OutOfPlane();
123 else if (fCentralityClass == "3050")InitD0toKpi2011PbPb3050();
c8617efe 124 else if (fCentralityClass=="010" && fIsCentScan) InitD0toKpi2011PbPb010CentScan();
3db225bf 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");
131 }
9341a8fa 132 }
133 else if (fCollisionType==2) {
86164e35 134 if (fCentralityClass=="0100") InitD0toKpi2013pPb0100();
fa90243b 135
136 if (fCentralityClass=="020V0A") InitD0toKpi2013pPb020V0A();
137 if (fCentralityClass=="2040V0A") InitD0toKpi2013pPb2040V0A();
138 if (fCentralityClass=="4060V0A") InitD0toKpi2013pPb4060V0A();
139 if (fCentralityClass=="60100V0A") InitD0toKpi2013pPb60100V0A();
140
141 if (fCentralityClass=="020ZNA") InitD0toKpi2013pPb020ZNA();
142 if (fCentralityClass=="2040ZNA") InitD0toKpi2013pPb2040ZNA();
143 if (fCentralityClass=="4060ZNA") InitD0toKpi2013pPb4060ZNA();
144 if (fCentralityClass=="60100ZNA") InitD0toKpi2013pPb60100ZNA();
145
44dd56b2 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();
152 }
9341a8fa 153 }
154 else AliFatal("Not yet implemented");
0233abe6 155 break;
3db225bf 156
0233abe6 157 case 2: // D+->Kpipi
a8c3295d 158 if (fCollisionType==0) {
159 if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
160 else InitDplustoKpipi2010pp();
9341a8fa 161 }
162 else if (fCollisionType==1) {
3db225bf 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");
172 }
173 if(fRunNumber == 11){
174 if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
c8617efe 175 else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010();
176 else if (fCentralityClass=="010" && fIsCentScan) InitDplustoKpipi2011PbPb010CentScan();
3db225bf 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();
830baf70 182 else if (fCentralityClass=="3050") InitDplustoKpipi2011PbPb3050();
3db225bf 183 else AliFatal("Not yet implemented");
184 }
9341a8fa 185 }
186 else if (fCollisionType==2) {
86164e35 187 if (fCentralityClass=="0100") InitDplustoKpipi2013pPb0100();
44dd56b2 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();
194 }
9341a8fa 195 }
196 else AliFatal("Not yet implemented");
0233abe6 197 break;
198 case 3: // D*->D0pi
fa2be30c 199 if (fCollisionType==0) {
200 if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
201 else InitDstartoD0pi2010pp();
9341a8fa 202 }
203 else if (fCollisionType==1) {
3db225bf 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");
214 }
215 if (fRunNumber == 11){
216 if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
c8617efe 217 else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010();
218 else if (fCentralityClass=="010" && fIsCentScan) InitDstartoD0pi2011PbPb010CentScan();
3db225bf 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();
830baf70 224 else if (fCentralityClass=="3050") InitDstartoD0pi2011PbPb3050();
3db225bf 225 else AliFatal("Not yet implemented");
226 }
227 }
228 else if (fCollisionType==2) {
86164e35 229 if (fCentralityClass=="0100") InitDstartoD0pi2013pPb0100();
fa90243b 230
231 if (fCentralityClass=="020V0A") InitDstartoD0pi2013pPb020V0A();
232 if (fCentralityClass=="2040V0A") InitDstartoD0pi2013pPb2040V0A();
233 if (fCentralityClass=="4060V0A") InitDstartoD0pi2013pPb4060V0A();
234 if (fCentralityClass=="60100V0A") InitDstartoD0pi2013pPb60100V0A();
235
236 if (fCentralityClass=="020ZNA") InitDstartoD0pi2013pPb020ZNA();
237 if (fCentralityClass=="2040ZNA") InitDstartoD0pi2013pPb2040ZNA();
238 if (fCentralityClass=="4060ZNA") InitDstartoD0pi2013pPb4060ZNA();
239 if (fCentralityClass=="60100ZNA") InitDstartoD0pi2013pPb60100ZNA();
240
44dd56b2 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();
247 }
5541b811 248 }
9341a8fa 249 else AliFatal("Not yet implemented");
0233abe6 250 break;
66247513 251 case 4: // D+s->KKpi
252 if (fCollisionType==0) InitDstoKKpi2010pp();
46988859 253 else if (fCollisionType==1) {
254 if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
255 else AliFatal("Not yet implemented");
9341a8fa 256 }
257 else if (fCollisionType==2) {
86164e35 258 if (fCentralityClass=="0100") InitDstoKKpi2013pPb0100();
46988859 259 }
27360a0e 260 else AliFatal("Not yet implemented");
66247513 261 break;
9ca30701 262 case 5: // Lc->pKpi
263 if (fCollisionType==0) InitLctopKpi2010pp();
264 else AliFatal("Not yet implemented");
265 break;
53c91ed7 266 case 6: // Lc->pK0S
267 if (fCollisionType==0) InitLctopK0S2010pp();
268 else AliFatal("Not yet implemented");
269 break;
66247513 270
0233abe6 271 default:
272 printf("Invalid decay type: %d\n",decay);
273 break;
274 }
275
0233abe6 276}
5541b811 277
0233abe6 278//--------------------------------------------------------------------------
5541b811 279void AliHFSystErr::InitD0toKpi2010pp() {
280 //
281 // D0->Kpi syst errors. Responsible: A. Rossi
282 // 2010 pp sample
0233abe6 283 //
284
5541b811 285 // Normalization
f1d72175 286 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 287 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
5541b811 288
289 // Branching ratio
f1d72175 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
5541b811 292
293 // Tracking efficiency
f1d72175 294 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
f8bc42d3 295 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08); // 8% (4% per track)
5541b811 296
297 // Raw yield extraction
f1d72175 298 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
5541b811 299 fRawYield->SetBinContent(1,1);
d2e5959f 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);
f1d72175 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);
0233abe6 306
5541b811 307 // Cuts efficiency (from cuts variation)
f1d72175 308 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
309 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
5541b811 310
311 // PID efficiency (from PID/noPID)
f1d72175 312 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
313 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03); // 3%
dc80fe6b 314 fPIDEff->SetBinContent(2,0.05); // 5%
5541b811 315
316 // MC dN/dpt
f1d72175 317 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 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);
5541b811 321
322 // particle-antiparticle
f1d72175 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);
5541b811 326
327 return;
0233abe6 328}
329//--------------------------------------------------------------------------
5541b811 330void AliHFSystErr::InitD0toKpi2010PbPb020() {
0233abe6 331 //
aebba8df 332 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 333 // 2010 PbPb sample, 0-20 CC
0233abe6 334 //
5541b811 335
336 // Normalization
337 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 338 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
5541b811 339
340 // Branching ratio
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
343
344 // Tracking efficiency
345 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 346 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);// Jacek, 5% per track
5541b811 347
348 // Raw yield extraction
349 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
350 fRawYield->SetBinContent(1,0);
351 fRawYield->SetBinContent(2,0);
6c388dee 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);
5541b811 355
356 // Cuts efficiency (from cuts variation)
357 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 358 fCutsEff->SetBinContent(1,0.);
359 fCutsEff->SetBinContent(2,0.);
cc4bb301 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);
6c388dee 363 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
5541b811 364
365 // PID efficiency (from PID/noPID)
366 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
6c388dee 367 for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
5541b811 368
369 // MC dN/dpt
370 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
cc4bb301 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);
b797b51b 374 for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
578bfbba 375
6c388dee 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);
0233abe6 383
5541b811 384 return;
385}
386//--------------------------------------------------------------------------
9a8b5080 387//--------------------------------------------------------------------------
388void AliHFSystErr::InitD0toKpi2011PbPb07half() {
389 //
390 // D0->Kpi syst errors. Responsible: A. Rossi
391 // 2011 PbPb sample, 0-7.5 CC
392 //
393
394 // Normalization
395 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 396 for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
9a8b5080 397
398 // Branching ratio
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
401
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
406
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);
fd8ae434 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);
9a8b5080 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
416
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
426
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
431
432 // MC dN/dpt
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
441
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);
449
450 return;
451}
93da8165 452
453//--------------------------------------------------------------------------
454void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
455 //
456 // D0->Kpi syst errors. Responsible: D. Caffarri
457 // 2011 PbPb sample, 30-50 CC InPlane
458 //
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
470
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
482}
483
484//--------------------------------------------------------------------------
485void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
486 //
487 // D0->Kpi syst errors. Responsible: D. Caffarri
488 // 2011 PbPb sample, 30-50 CC OutOfPlane
489 //
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
501
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
513}
514
9a8b5080 515//--------------------------------------------------------------------------
5541b811 516void AliHFSystErr::InitD0toKpi2010PbPb4080() {
93da8165 517 //
aebba8df 518 // D0->Kpi syst errors. Responsible: A. Rossi
5541b811 519 // 2010 PbPb sample, 40-80 CC
520 //
521
0233abe6 522 // Normalization
523 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 524 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
0233abe6 525
526 // Branching ratio
527 fBR = new TH1F("fBR","fBR",20,0,20);
844405bc 528 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
0233abe6 529
aebba8df 530 // Tracking efficiency
531 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 532 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
aebba8df 533
534
0233abe6 535 // Raw yield extraction
536 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
5541b811 537 fRawYield->SetBinContent(1,0);
538 fRawYield->SetBinContent(2,0);
6c388dee 539 for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
666cf569 540 //for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,0);
0233abe6 541
542 // Cuts efficiency (from cuts variation)
543 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
aebba8df 544 fCutsEff->SetBinContent(1,0.);
545 fCutsEff->SetBinContent(2,0.);
cc4bb301 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);
6c388dee 549 for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
0233abe6 550
551 // PID efficiency (from PID/noPID)
552 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
6c388dee 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);
0233abe6 556
557 // MC dN/dpt
558 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
cc4bb301 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);
563
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);
0233abe6 567
568 return;
569}
a8c3295d 570
c8617efe 571//--------------------------------------------------------------------------
572void AliHFSystErr::InitD0toKpi2011PbPb010() {
573 //
574 // D0->Kpi syst errors. Responsible: A. Rossi
575 // 2011 PbPb sample, 0-10 CC
576 //
577
578 AliInfo(" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
579
580 // Normalization
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
583
584 // Branching ratio
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
587
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
591
592
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);
605
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);
617
618
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);
624
625 // MC dN/dpt
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.);
632
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);
636
637 return;
638}
639
830baf70 640//--------------------------------------------------------------------------
641void AliHFSystErr::InitD0toKpi2011PbPb3050() {
642 //
643 // D0->Kpi syst errors. Responsible: D. Caffarri
644 // 2011 PbPb sample, 30-50 CC
645 //
646
647 // Normalization
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
650
651 // Branching ratio
652 fBR = new TH1F("fBR","fBR",20,0,20);
653 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012); // 1.2% PDG2010
654
655 // Tracking efficiency
656 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
657 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); // Jacek, 5% per track
658
659
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);
674
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);
677
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);
687
688
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);
694
695 // MC dN/dpt
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);
706
707 return;
708}
709
710
a8c3295d 711//--------------------------------------------------------------------------
712void AliHFSystErr::InitD0toKpi2010ppLowEn() {
713 //
714 // D0->Kpi syst errors. Low energy run
715 // 2011 2.76 TeV pp sample
716 //
717 AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
718
719 // Normalization
720 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 721 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 4% error on sigmaV0and
a8c3295d 722
723 // Branching ratio
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
726
727 // Tracking efficiency
728 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 729 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10); //10% (5% per track)
a8c3295d 730
731 // Raw yield extraction
732 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
733 fRawYield->SetBinContent(1,1);
f9a85c82 734 for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
a8c3295d 735
736 // Cuts efficiency (from cuts variation)
737 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
6c388dee 738 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10); // 10%
f9a85c82 739 fCutsEff->SetBinContent(2,0.20);
740 for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
741
a8c3295d 742
743 // PID efficiency (from PID/noPID)
744 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
f9a85c82 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%
a8c3295d 748
749 // MC dN/dpt
750 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 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);
a8c3295d 754
f9a85c82 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.);
a8c3295d 761
762 return;
763}
86164e35 764//_________________________________________________________________________
765void AliHFSystErr::InitD0toKpi2013pPb0100(){
bd332272 766 //
767 // D0->Kpi syst errors. p-Pb data sample
768 //
769 AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
86164e35 770
771 fNorm = new TH1F("fNorm","fNorm",24,0,24);
9341a8fa 772 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
86164e35 773
774 // Branching ratio
9341a8fa 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
86164e35 777
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
781
782 // Raw yield extraction
783 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
9341a8fa 784 fRawYield->SetBinContent(1,0.0);
bc95efa5 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);
86164e35 789
790 // Cuts efficiency (from cuts variation)
791 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
9341a8fa 792 fCutsEff->SetBinContent(1,0.0);
e5f8509c 793 fCutsEff->SetBinContent(2,0.08);
9341a8fa 794 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
86164e35 795
796 // PID efficiency (from PID/noPID)
797 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
bc95efa5 798 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
86164e35 799
800 // MC dN/dpt
801 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
9341a8fa 802 fMCPtShape->SetBinContent(1,0.0);
803 fMCPtShape->SetBinContent(2,0.02);
e5f8509c 804 for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
86164e35 805
806 // particle-antiparticle
807 /*
808 fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",20,0,20);
9341a8fa 809 for(Int_t i=1;i<=20;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
86164e35 810 */
811 return;
812
813
814}
a8c3295d 815
44dd56b2 816
0233abe6 817//--------------------------------------------------------------------------
5541b811 818void AliHFSystErr::InitDplustoKpipi2010pp() {
0233abe6 819 //
820 // D+->Kpipi syst errors. Responsible: R. Bala
5541b811 821 // 2010 pp sample
0233abe6 822 //
823
1f91ef56 824
825// Normalization
826 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 827 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
626bcf57 828
829 // Branching ratio
1f91ef56 830 fBR = new TH1F("fBR","fBR",24,0,24);
f7c8a2f1 831 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
626bcf57 832
833 // Tracking efficiency
1f91ef56 834 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
f8bc42d3 835 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
626bcf57 836
837
838 // Raw yield extraction
1f91ef56 839 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
d2e5959f 840 fRawYield->SetBinContent(1,1);
1f91ef56 841 fRawYield->SetBinContent(2,0.25);
842 fRawYield->SetBinContent(3,0.25);
d2e5959f 843 fRawYield->SetBinContent(4,0.20);
1f91ef56 844 fRawYield->SetBinContent(5,0.09);
845 fRawYield->SetBinContent(6,0.09);
d2e5959f 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);
1f91ef56 848
626bcf57 849 // Cuts efficiency (from cuts variation)
1f91ef56 850 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
851 for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10); // 10%
626bcf57 852
853 // PID efficiency (from PID/noPID)
1f91ef56 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%
861
862 // MC dN/dpt
863 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
864 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
578bfbba 865 fMCPtShape->SetBinContent(1,0.03);
866 fMCPtShape->SetBinContent(2,0.03);
626bcf57 867
868
869 // particle-antiparticle
870 /*
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%
876 */
0233abe6 877 return;
878}
1f91ef56 879
66247513 880//--------------------------------------------------------------------------
881void AliHFSystErr::InitDstoKKpi2010pp() {
882 //
883 // D+s->KKpi syst errors. Responsible: G.M. Innocenti
884 // 2010 pp sample
885 //
886
887
888// Normalization
889 fNorm = new TH1F("fNorm","fNorm",12,0,12);
494408cc 890 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035); // 3.5% error on sigmaV0and
66247513 891
892 // Branching ratio
893 fBR = new TH1F("fBR","fBR",12,0,12);
494408cc 894 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06); // 0.14/2.32 PDG2010
66247513 895
896 // Tracking efficiency
897 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
494408cc 898 for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12); // 12% (4% per track)
66247513 899
900
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);
915
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%
919
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%
923
494408cc 924 // MC dN/dpt
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);
66247513 931
932
933 // particle-antiparticle
934 /*
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%
940 */
941 return;
942}
943
944
436b185f 945//--------------------------------------------------------------------------
946void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
947 //
948 // D+->Kpipi syst errors. Responsible: E. Bruna
949 // 2011 PbPb sample, 0-7.5% CC
950 //
951
952 // Normalization
953 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 954 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
436b185f 955
956 // Branching ratio
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
959
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
963
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%
eb4886c7 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);
436b185f 972 for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
eb4886c7 973 for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
436b185f 974 for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
975
976
977 // Cuts efficiency (from cuts variation)
978 fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
a7cfa9ed 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);
eb4886c7 981
436b185f 982
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%
986
a7cfa9ed 987 // MC dN/dpt (24/7/2012)
436b185f 988 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
3d446ac8 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);
436b185f 992
993
994 // particle-antiparticle
995 /*
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%
1001 */
1002
1003 return;
1004}
c8617efe 1005//_________________________________________________________________
1006void AliHFSystErr::InitDplustoKpipi2011PbPb010() {
1007 //
1008 // D+->Kpipi syst errors. Responsible: E. Bruna, C. Bedda
1009 // 2011 PbPb sample, 0-10% CC
1010 //
1011
1012 // Normalization
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
1015
1016 // Branching ratio
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
1019
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
1023
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);
1035
1036
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);
1040
1041
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%
1045
1046 // MC dN/dpt
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);
1053
1054
1055 // particle-antiparticle
1056 /*
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%
1062 */
1063
1064 return;
1065}
1066
830baf70 1067//___________________________________________________________________
1068void AliHFSystErr::InitDplustoKpipi2011PbPb3050() {
1069 //
1070 // D+->Kpipi syst errors. Responsible: R. Bala
1071 // 2011 PbPb sample, 30-50% CC
1072 //
1073 AliInfo(" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
1074 // Normalization
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
1077
1078 // Branching ratio
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
1081
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
1085
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);
1090
1091
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);
1096
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%
1100
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);
1106
1107
1108 // particle-antiparticle
1109 /*
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%
1115 */
1116
1117 return;
1118}
1119
436b185f 1120
46988859 1121//--------------------------------------------------------------------------
1122void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
1123 //
1124 // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
1125 // 2011 PbPb sample, 0-7.5% CC
1126 //
1127
1128 // Normalization
1129 fNorm = new TH1F("fNorm","fNorm",12,0,12);
a713683e 1130 for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
46988859 1131
1132 // Branching ratio
1133 fBR = new TH1F("fBR","fBR",12,0,12);
d92cd593 1134 for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053); // 0.12/2.28 PDG2012
46988859 1135
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
1139
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);
1144
1145
1146
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%
1150
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%
1154
1155 // MC dN/dpt
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);
1162
1163 // particle-antiparticle
1164 /*
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%
1170 */
1171
1172 return;
1173}
1174
86164e35 1175//_________________________________________________________________________
1176void AliHFSystErr::InitDstoKKpi2013pPb0100(){
bd332272 1177 //
86164e35 1178 // Dstar syst
bd332272 1179 AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
86164e35 1180
1181 //Normalization
1182 fNorm = new TH1F("fNorm","fNorm",24,0,24);
1183 for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00); //
1184
1185 // Branching ratio
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
1188
1189 // Tracking efficiency
1190 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
59604f68 1191 for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09); // 3% per track
86164e35 1192
1193 // Raw yield extraction
1194 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
59604f68 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);
1199
86164e35 1200 // Cuts efficiency (from cuts variation)
1201 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
59604f68 1202 for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15); //
86164e35 1203
1204 // PID efficiency (from PID/noPID)
1205 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
59604f68 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);
86164e35 1208
1209 // MC dN/dpt
1210 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
59604f68 1211 for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
86164e35 1212
1213 // particle-antiparticle
1214 /*
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%
1220 */
1221 return;
86164e35 1222}
1223
1224
0233abe6 1225//--------------------------------------------------------------------------
5541b811 1226void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
1227 //
1228 // D+->Kpipi syst errors. Responsible: ??
1229 // 2010 PbPb sample, 0-20 CC
1230 //
aebba8df 1231
aebba8df 1232 // Normalization
1233 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 1234 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
aebba8df 1235
1236 // Branching ratio
1237 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1238 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
aebba8df 1239
1240 // Tracking efficiency
1241 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 1242 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
aebba8df 1243
1244 // Raw yield extraction
1245 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
6c388dee 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);
cb142b90 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);
aebba8df 1255
1256 // Cuts efficiency (from cuts variation)
1257 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
666cf569 1258 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
aebba8df 1259
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%
aebba8df 1263
f7c8a2f1 1264 // MC dN/dpt (2/2/2012)
aebba8df 1265 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 1266 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
f7c8a2f1 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);
aebba8df 1270
1271
1272 // particle-antiparticle
1273 /*
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%
1279 */
6c388dee 1280
5541b811 1281 return;
1282}
1283
1284//--------------------------------------------------------------------------
1285void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
1286 //
1287 // D+->Kpipi syst errors. Responsible: ??
1288 // 2010 PbPb sample, 40-80 CC
1289 //
1290
aebba8df 1291
1292 // Normalization
1293 fNorm = new TH1F("fNorm","fNorm",20,0,20);
d5e15690 1294 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
aebba8df 1295
1296 // Branching ratio
1297 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1298 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1%
aebba8df 1299
1300 // Tracking efficiency
1301 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
69bcc748 1302 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
aebba8df 1303
1304
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);
6970b225 1310 fRawYield->SetBinContent(4,0.15);
cb142b90 1311 fRawYield->SetBinContent(5,0.05);
1312 fRawYield->SetBinContent(6,0.05);
1313 fRawYield->SetBinContent(7,0.15);
1314 fRawYield->SetBinContent(8,0.15);
6970b225 1315 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
aebba8df 1316 for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1); //5 to 10%
1317
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%
1321
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%
1326
1327
f7c8a2f1 1328 // MC dN/dpt (2/2/2012)
aebba8df 1329 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
578bfbba 1330 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
f7c8a2f1 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);
aebba8df 1334
1335
1336 // particle-antiparticle
1337 /*
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%
1343 */
5541b811 1344 return;
1345}
1346
a8c3295d 1347//--------------------------------------------------------------------------
1348void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
1349
1350 //
1351 // D+->Kpipi syst errors. Responsible: R. Bala
1352 // 2011 2.76 TeV pp sample
1353 //
1354 AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
1355
1356 // Normalization
1357 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 1358 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
a8c3295d 1359
1360 // Branching ratio
1361 fBR = new TH1F("fBR","fBR",20,0,20);
f7c8a2f1 1362 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021); // 2.1% PDG2010
a8c3295d 1363
1364 // Tracking efficiency
1365 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 1366 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); // 3% (1% per track)
a8c3295d 1367
1368 // Raw yield extraction
1369 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
1370 fRawYield->SetBinContent(1,1);
1371 fRawYield->SetBinContent(2,1);
f9a85c82 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%
a8c3295d 1376
1377 // Cuts efficiency (from cuts variation)
1378 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
f9a85c82 1379 for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15); // 10%
a8c3295d 1380
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%
f9a85c82 1384 fPIDEff->SetBinContent(3,0.10); // 13%
1385 fPIDEff->SetBinContent(4,0.10); // 13%
a8c3295d 1386
1387 // MC dN/dpt (copied from D0 : will update later)
1388 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 1389 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
578bfbba 1390 fMCPtShape->SetBinContent(1,0.03);
1391 fMCPtShape->SetBinContent(2,0.03);
a8c3295d 1392
1393 return;
1394}
1395
86164e35 1396//_________________________________________________________________________
1397void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
bd332272 1398 //
1399 // Dplus syst
1400 AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
86164e35 1401
1402 fNorm = new TH1F("fNorm","fNorm",24,0,24);
3e92d3ee 1403 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02); //
86164e35 1404
2d325332 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
86164e35 1407
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
4c1e4a80 1411
86164e35 1412 // Raw yield extraction
1413 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
3e92d3ee 1414 fRawYield->SetBinContent(1,0.);
1415 fRawYield->SetBinContent(2,0.1);
1416 fRawYield->SetBinContent(3,0.08);
1417 fRawYield->SetBinContent(4,0.05);
4c1e4a80 1418 fRawYield->SetBinContent(5,0.05);
3e92d3ee 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);
4c1e4a80 1425
86164e35 1426 // Cuts efficiency (from cuts variation)
1427 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
3e92d3ee 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);
86164e35 1439
1440 // PID efficiency (from PID/noPID)
1441 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
4c1e4a80 1442 fPIDEff->SetBinContent(1,1);
3e92d3ee 1443 fPIDEff->SetBinContent(2,0.1);
1444 for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0); //
86164e35 1445
4c1e4a80 1446 // MC dN/dpt
86164e35 1447 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
3e92d3ee 1448 for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
86164e35 1449
1450 // particle-antiparticle
1451 /*
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%
1457 */
1458 return;
1459
1460
1461}
1462
5541b811 1463//--------------------------------------------------------------------------
1464void AliHFSystErr::InitDstartoD0pi2010pp() {
0233abe6 1465 //
2cf98d0d 1466 // D*+->D0pi syst errors. Responsible: A. Grelli, Y. Wang
5541b811 1467 // 2010 pp sample
0233abe6 1468 //
1469
2cf98d0d 1470 // Normalization
fa2be30c 1471 fNorm = new TH1F("fNorm","fNorm",24,0,24);
27360a0e 1472 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035); // 4% error on sigmaV0and
fa2be30c 1473
1474 // Branching ratio
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
1477
1478 // Tracking efficiency
1479 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
1480 fTrackingEff->SetBinContent(1,1.0);
f8bc42d3 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)
fa2be30c 1485
1486
1487 // Raw yield extraction
1488 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
1489 fRawYield->SetBinContent(1,1.0);
d2e5959f 1490 fRawYield->SetBinContent(2,0.10);
fa2be30c 1491 fRawYield->SetBinContent(3,0.04);
d2e5959f 1492 fRawYield->SetBinContent(4,0.03);
fa2be30c 1493 fRawYield->SetBinContent(5,0.03);
d2e5959f 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%
fa2be30c 1500
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%
1505
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%
1509
1510
1511 // MC dN/dpt (copied from D0 : will update later)
1512 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 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);
fa2be30c 1516
1517 return;
1518
1519
1520}
1521//--------------------------------------------------------------------------
1522void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
1523
1524 //
1525 // D+->Kpipi syst errors. Responsible: A. Grelli
1526 // 2011 2.76 TeV pp sample
1527 //
1528 AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
1529
1530// Normalization
2cf98d0d 1531 fNorm = new TH1F("fNorm","fNorm",20,0,20);
f9a85c82 1532 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031); // 10% error on sigmaV0and
2cf98d0d 1533
1534 // Branching ratio
1535 fBR = new TH1F("fBR","fBR",20,0,20);
8accac87 1536 for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
2cf98d0d 1537
1538 // Tracking efficiency
1539 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
27360a0e 1540 for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15); //10% (to be checked!!)
2cf98d0d 1541
1542 // Raw yield extraction
1543 fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
8accac87 1544 fRawYield->SetBinContent(1,1);
fa2be30c 1545 fRawYield->SetBinContent(2,1);
f9a85c82 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);
fa2be30c 1556 for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
2cf98d0d 1557
1558 // Cuts efficiency (from cuts variation)
1559 fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
f9a85c82 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);
2cf98d0d 1571
1572 // PID efficiency (from PID/noPID)
1573 fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
fa2be30c 1574 for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05); // 10%
2cf98d0d 1575
fa2be30c 1576 // MC dN/dpt
2cf98d0d 1577 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
27360a0e 1578 for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
578bfbba 1579 fMCPtShape->SetBinContent(1,0.03);
1580 fMCPtShape->SetBinContent(2,0.03);
2cf98d0d 1581
0233abe6 1582 return;
1583}
6c388dee 1584
fa2be30c 1585//------------------------------------------------------------------------
5541b811 1586void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
1587 //
fa2be30c 1588 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 1589 // 2010 PbPb sample, 0-20 CC
1590 //
1591
fa2be30c 1592 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
1593
1594 // Normalization
1595 fNorm = new TH1F("fNorm","fNorm",24,0,24);
d5e15690 1596 for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
fa2be30c 1597
1598 // Branching ratio
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
1601
1602 // Tracking efficiency
1603 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1604 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1605
1606
1607 // Raw yield extraction
1608 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
6c388dee 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);
fa2be30c 1614
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%
6fd587a4 1618 fCutsEff->SetBinContent(4,0.15);
1619 fCutsEff->SetBinContent(5,0.15);
1620 fCutsEff->SetBinContent(6,0.15);
fa2be30c 1621
1622 // PID efficiency (from PID/noPID)
1623 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
6c388dee 1624 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
fa2be30c 1625
1626
6fd587a4 1627 // MC dN/dpt (from study on D* pt shape)
fa2be30c 1628 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
6fd587a4 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);
1639
9a8b5080 1640 return;
1641
1642}
1643// ----------------------------- 2011 ----------------------------------
1644void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
1645 //
1646 // D*+->D0pi syst errors. Responsible: A. Grelli
1647 // 2011 PbPb sample, 0-7.5 CC
1648 //
1649
1650 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
1651
1652
1653 // Normalization
1654 fNorm = new TH1F("fNorm","fNorm",36,0,36);
a713683e 1655 for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048); // TAA and pp norm
9a8b5080 1656
1657 // Branching ratio
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
1660
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
1664
1665
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%
ada850b0 1669 fRawYield->SetBinContent(4,0.2);
9a8b5080 1670 fRawYield->SetBinContent(5,0.10);
ada850b0 1671 fRawYield->SetBinContent(6,0.10);
1672 fRawYield->SetBinContent(7,0.08);
9a8b5080 1673 for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15); //4%
1674
9a8b5080 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%
9a8b5080 1678
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%
ada850b0 1682 fPIDEff->SetBinContent(4,0.09);
9a8b5080 1683
1684 // MC dN/dpt (from study on D* pt shape)
1685 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
273a7638 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);
9a8b5080 1692 fMCPtShape->SetBinContent(9,0.03);
1693 fMCPtShape->SetBinContent(10,0.03);
1694 fMCPtShape->SetBinContent(11,0.03);
1695 fMCPtShape->SetBinContent(12,0.03);
1696
6fd587a4 1697
fa2be30c 1698
1699 return;
1700
1701}
1702//-------------------------------------------------------------------------
1703void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
1704 //
1705 // D*+->D0pi syst errors. Responsible: A. Grelli
1706 // 2010 PbPb sample, 20-40 CC
1707 //
1708
1709 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
1710
1711 // Normalization
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
1714
1715 // Branching ratio
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
1718
1719 // Tracking efficiency
1720 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1721 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1722
1723
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%
1727
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%
1731
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%
1735
1736
1737 // MC dN/dpt (copied from D0 : will update later)
1738 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
578bfbba 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);
fa2be30c 1742
5541b811 1743 return;
fa2be30c 1744
5541b811 1745}
6c388dee 1746
5541b811 1747//--------------------------------------------------------------------------
1748void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
1749 //
fa2be30c 1750 // D*+->D0pi syst errors. Responsible: A. Grelli
5541b811 1751 // 2010 PbPb sample, 40-80 CC
1752 //
1753
fa2be30c 1754 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
1755
1756 // Normalization
1757 fNorm = new TH1F("fNorm","fNorm",24,0,24);
d5e15690 1758 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); // TAA and pp norm
fa2be30c 1759
1760 // Branching ratio
1761 fBR = new TH1F("fBR","fBR",24,0,24);
1762 for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015); // 1.5% PDG2010
1763
1764 // Tracking efficiency
1765 fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
69bcc748 1766 for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
fa2be30c 1767
1768
1769 // Raw yield extraction
1770 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
6c388dee 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);
1788
1789
fa2be30c 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%
1793
1794 // PID efficiency (from PID/noPID)
1795 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
6c388dee 1796 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05); // 3%
fa2be30c 1797
1798
1799 // MC dN/dpt (copied from D0 : will update later)
1800 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
6fd587a4 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);
1809
5541b811 1810 return;
fa2be30c 1811
5541b811 1812}
6c388dee 1813
c8617efe 1814// ----------------------------- 2011 ----------------------------------
1815void AliHFSystErr::InitDstartoD0pi2011PbPb010() {
1816 //
1817 // D*+->D0pi syst errors. Responsible: A. Grelli
1818 // 2011 PbPb sample, 0-10% CC
1819 //
1820
1821 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
1822
1823
1824 // Normalization
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
1827
1828 // Branching ratio
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
1831
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
1835
1836
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%
1848
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%
1854
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%
1858
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);
1865
1866 return;
1867
1868}
1869
830baf70 1870//----------------------------------------------------------------------
1871void AliHFSystErr::InitDstartoD0pi2011PbPb3050() {
1872 //
1873 // D*+->D0pi syst errors. Responsible: A. Grelli
1874 // 2011 PbPb sample, 30-50 CC
1875 //
1876
1877 AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
1878
1879 // Normalization
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
1882
1883 // Branching ratio
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
1886
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
1890
1891
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);
1911
1912
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%
1925
1926
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);
1937
1938 return;
1939}
86164e35 1940//_________________________________________________________________________
1941void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
f21cf4a7 1942 // Dstar syst in pPb 2013 MB
1943 //
bd332272 1944 AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
86164e35 1945
86164e35 1946 //Normalization
1947 fNorm = new TH1F("fNorm","fNorm",24,0,24);
fa90243b 1948 for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07); //
86164e35 1949
f21cf4a7 1950 // Branching ratio
86164e35 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
1953
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
1957
1958 // Raw yield extraction
1959 fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
f21cf4a7 1960 fRawYield->SetBinContent(1,1.0);
30dc63e6 1961 fRawYield->SetBinContent(2,0.08);
1962 fRawYield->SetBinContent(3,0.02);
f21cf4a7 1963 fRawYield->SetBinContent(4,0.02);
30dc63e6 1964 fRawYield->SetBinContent(5,0.05);
1965 fRawYield->SetBinContent(6,0.02);
f21cf4a7 1966 fRawYield->SetBinContent(7,0.04);
1967 fRawYield->SetBinContent(8,0.04);
30dc63e6 1968 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
f21cf4a7 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);
1971
86164e35 1972 // Cuts efficiency (from cuts variation)
1973 fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
f21cf4a7 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);
86164e35 1977
1978 // PID efficiency (from PID/noPID)
1979 fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
3e92d3ee 1980 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
86164e35 1981
1982 // MC dN/dpt
1983 fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
c9569eb5 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);
86164e35 1987
1988 // particle-antiparticle
1989 /*
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%
1995 */
fa90243b 1996
1997 return;
86164e35 1998
1999}
2000
414bc6a5 2001//--------------------------------------------------------------------------
2002void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
2003 // define errors for RAA vs. centrality
2004 InitD0toKpi2010PbPb020();
41cfcd39 2005 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 2006 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
41cfcd39 2007 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2008}
2009//--------------------------------------------------------------------------
2010void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
2011 // define errors for RAA vs. centrality
2012 InitD0toKpi2010PbPb020();
41cfcd39 2013 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 2014 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
41cfcd39 2015 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2016}
2017//--------------------------------------------------------------------------
2018void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
2019 // define errors for RAA vs. centrality
41cfcd39 2020 InitD0toKpi2010PbPb4080();
2021 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
27360a0e 2022 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.14);
41cfcd39 2023 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2024}
2025//--------------------------------------------------------------------------
2026void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
41cfcd39 2027 // define errors for RAA vs. centrality
414bc6a5 2028 InitD0toKpi2010PbPb4080();
41cfcd39 2029 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
27360a0e 2030 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.11);
41cfcd39 2031 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2032}
2033//--------------------------------------------------------------------------
2034void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
41cfcd39 2035 // define errors for RAA vs. centrality
414bc6a5 2036 InitD0toKpi2010PbPb4080();
41cfcd39 2037 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
27360a0e 2038 for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
41cfcd39 2039 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2040}
2041//--------------------------------------------------------------------------
2042void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
2043 // define errors for RAA vs. centrality
2044 InitDplustoKpipi2010PbPb020();
7b03f1ec 2045 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
eefab142 2046 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
7b03f1ec 2047
414bc6a5 2048}
2049//--------------------------------------------------------------------------
2050void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
2051 // define errors for RAA vs. centrality
2052 InitDplustoKpipi2010PbPb020();
7b03f1ec 2053 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
eefab142 2054 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2055}
2056//--------------------------------------------------------------------------
2057void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
2058 // define errors for RAA vs. centrality
2059 InitDplustoKpipi2010PbPb020();
7b03f1ec 2060 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
eefab142 2061 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
414bc6a5 2062}
2063//--------------------------------------------------------------------------
2064void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
2065 // define errors for RAA vs. centrality
2066 InitDplustoKpipi2010PbPb4080();
7b03f1ec 2067 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
eefab142 2068 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2069}
2070//--------------------------------------------------------------------------
2071void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
2072 // define errors for RAA vs. centrality
2073 InitDplustoKpipi2010PbPb4080();
7b03f1ec 2074 for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
eefab142 2075 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
414bc6a5 2076}
2077
2078//--------------------------------------------------------------------------
2079void 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);
eefab142 2083 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
414bc6a5 2084}
2085//--------------------------------------------------------------------------
2086void 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);
eefab142 2090 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
414bc6a5 2091}
2092//--------------------------------------------------------------------------
2093void 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);
eefab142 2097 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
414bc6a5 2098}
2099//--------------------------------------------------------------------------
2100void 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);
eefab142 2104 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
414bc6a5 2105}
2106//--------------------------------------------------------------------------
2107void 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);
eefab142 2111 for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
414bc6a5 2112}
2113
3db225bf 2114//--------------------------------------------------------------------------
2115void AliHFSystErr::InitD0toKpi2011PbPb010CentScan(){
2116 // define errors for RAA vs. centrality 2011
2117 InitD0toKpi2011PbPb07half();
2118
2119}
2120//--------------------------------------------------------------------------
2121void AliHFSystErr::InitD0toKpi2011PbPb1020CentScan(){
2122 // define errors for RAA vs. centrality 2011
2123 InitD0toKpi2011PbPb07half();
2124
2125}
2126//--------------------------------------------------------------------------
2127void AliHFSystErr::InitD0toKpi2011PbPb2030CentScan(){
2128 // define errors for RAA vs. centrality 2011
2129 InitD0toKpi2011PbPb07half();
2130
2131}
2132//--------------------------------------------------------------------------
2133void AliHFSystErr::InitD0toKpi2011PbPb3040CentScan(){
2134 // define errors for RAA vs. centrality 2011
2135 InitD0toKpi2011PbPb07half();
2136
2137}
2138//--------------------------------------------------------------------------
2139void AliHFSystErr::InitD0toKpi2011PbPb4050CentScan(){
2140 // define errors for RAA vs. centrality 2011
2141 InitD0toKpi2011PbPb07half();
2142
2143}
2144//--------------------------------------------------------------------------
2145void AliHFSystErr::InitD0toKpi2010PbPb5080CentScan(){
2146 // define errors for RAA vs. centrality 2011
2147 InitD0toKpi2011PbPb07half();
2148
2149}
2eab2b1d 2150
3db225bf 2151//--------------------------------------------------------------------------
2152void AliHFSystErr::InitDplustoKpipi2011PbPb010CentScan(){
2153 // define errors for RAA vs. centrality 2011
2154 InitDplustoKpipi2011PbPb07half();
2eab2b1d 2155
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);
2159
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);
2163
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);
3db225bf 2167}
2eab2b1d 2168
3db225bf 2169//--------------------------------------------------------------------------
2170void AliHFSystErr::InitDplustoKpipi2011PbPb1020CentScan(){
2171 // define errors for RAA vs. centrality 2011
2172 InitDplustoKpipi2011PbPb07half();
2eab2b1d 2173
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);
2177
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);
3db225bf 2181
2eab2b1d 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);
3db225bf 2185}
2eab2b1d 2186
3db225bf 2187//--------------------------------------------------------------------------
2188void AliHFSystErr::InitDplustoKpipi2011PbPb2030CentScan(){
2189 // define errors for RAA vs. centrality 2011
2190 InitDplustoKpipi2011PbPb07half();
2191
2eab2b1d 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);
2195
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);
2199
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);
3db225bf 2203}
2eab2b1d 2204
3db225bf 2205//--------------------------------------------------------------------------
2206void AliHFSystErr::InitDplustoKpipi2011PbPb3040CentScan(){
2207 // define errors for RAA vs. centrality 2011
2208 InitDplustoKpipi2011PbPb07half();
2eab2b1d 2209
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);
2213
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);
2217
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);
3db225bf 2221}
2eab2b1d 2222
3db225bf 2223//--------------------------------------------------------------------------
2224void AliHFSystErr::InitDplustoKpipi2011PbPb4050CentScan(){
2225 // define errors for RAA vs. centrality 2011
2226 InitDplustoKpipi2011PbPb07half();
2eab2b1d 2227
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);
2231
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);
2235
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);
3db225bf 2239}
2eab2b1d 2240
3db225bf 2241//--------------------------------------------------------------------------
2242void AliHFSystErr::InitDplustoKpipi2010PbPb5080CentScan(){
2243 // define errors for RAA vs. centrality 2011
2eab2b1d 2244 InitDplustoKpipi2010PbPb4080();
2245
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);
2249
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%
3db225bf 2254}
2eab2b1d 2255
3db225bf 2256//--------------------------------------------------------------------------
2257void AliHFSystErr::InitDstartoD0pi2011PbPb010CentScan(){
2258 // define errors for RAA vs. centrality 2011
2259 InitDstartoD0pi2011PbPb07half();
2260
2261}
2262//--------------------------------------------------------------------------
2263void AliHFSystErr::InitDstartoD0pi2011PbPb1020CentScan(){
2264 // define errors for RAA vs. centrality 2011
2265 InitDstartoD0pi2011PbPb07half();
2266
2267}
2268//--------------------------------------------------------------------------
2269void AliHFSystErr::InitDstartoD0pi2011PbPb2030CentScan(){
2270 // define errors for RAA vs. centrality 2011
2271 InitDstartoD0pi2011PbPb07half();
2272
2273}
2274//--------------------------------------------------------------------------
2275void AliHFSystErr::InitDstartoD0pi2011PbPb3040CentScan(){
2276 // define errors for RAA vs. centrality 2011
2277 InitDstartoD0pi2011PbPb07half();
2278
2279}
2280//--------------------------------------------------------------------------
2281void AliHFSystErr::InitDstartoD0pi2011PbPb4050CentScan(){
2282 // define errors for RAA vs. centrality 2011
2283 InitDstartoD0pi2011PbPb07half();
2284
2285}
2286//--------------------------------------------------------------------------
2287void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
2288 // define errors for RAA vs. centrality 2011
2289 InitDstartoD0pi2011PbPb07half();
2290
2291}
2292
44dd56b2 2293//_________________________________________________________________________
2294void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
2295 InitD0toKpi2013pPb0100();
2296
a1f0af09 2297 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2298
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);
44dd56b2 2301 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2302
44dd56b2 2303 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2304
a1f0af09 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);
2308
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);
44dd56b2 2312
2313
2314};
2315//_________________________________________________________________________
2316void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
2317 InitD0toKpi2013pPb0100();
a1f0af09 2318
2319 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
44dd56b2 2320
a1f0af09 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);
44dd56b2 2323 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
2324
44dd56b2 2325 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2326
a1f0af09 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);
2330
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);
44dd56b2 2334};
2335//_________________________________________________________________________
2336void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
a1f0af09 2337 InitD0toKpi2013pPb0100();
2338
2339 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
44dd56b2 2340
a1f0af09 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);
44dd56b2 2343 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
2344
44dd56b2 2345 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2346
a1f0af09 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);
2350
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);
44dd56b2 2354
2355};
2356//_________________________________________________________________________
2357void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
2358 InitD0toKpi2013pPb0100();
2359
a1f0af09 2360 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
44dd56b2 2361
a1f0af09 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);
2365
2366 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2367
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);
2371
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);
44dd56b2 2375
2376};
2377//_________________________________________________________________________
2378void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
2379 InitD0toKpi2013pPb0100();
2380
a1f0af09 2381 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2382
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);
44dd56b2 2386
44dd56b2 2387 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2388
a1f0af09 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);
2392
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);
44dd56b2 2396
2397};
2398
2399//_________________________________________________________________________
2400void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
2401 InitDplustoKpipi2013pPb0100();
2402
a1f0af09 2403 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
44dd56b2 2404
a1f0af09 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);
2408
44dd56b2 2409 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2410
a1f0af09 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);
2414
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);
44dd56b2 2418
2419};
2420//_________________________________________________________________________
2421void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
2422 InitDplustoKpipi2013pPb0100();
a1f0af09 2423 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
44dd56b2 2424
a1f0af09 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);
44dd56b2 2427 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
2428
44dd56b2 2429 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2430
a1f0af09 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);
2434
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);
44dd56b2 2438
2439};
2440//_________________________________________________________________________
2441void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
2442 InitDplustoKpipi2013pPb0100();
a1f0af09 2443 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
44dd56b2 2444
a1f0af09 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);
44dd56b2 2447 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2448
44dd56b2 2449 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2450
a1f0af09 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);
2454
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);
44dd56b2 2458
2459};
2460//_________________________________________________________________________
2461void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
2462 InitDplustoKpipi2013pPb0100();
a1f0af09 2463 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
2464
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);
44dd56b2 2467 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
a1f0af09 2468
44dd56b2 2469 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2470
a1f0af09 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);
44dd56b2 2474
a1f0af09 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);
2478
44dd56b2 2479};
2480//_________________________________________________________________________
2481void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
2482 InitDplustoKpipi2013pPb0100();
2483
a1f0af09 2484 AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
2485
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);
44dd56b2 2489
44dd56b2 2490 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2491
a1f0af09 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);
44dd56b2 2495
a1f0af09 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);
2499
44dd56b2 2500};
2501
2502
2503//_________________________________________________________________________
2504void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
2505 InitDstartoD0pi2013pPb0100();
2506
a1f0af09 2507 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
2508
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);
44dd56b2 2511 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2512
44dd56b2 2513 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2514
a1f0af09 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);
2518
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);
44dd56b2 2522
2523};
2524//_________________________________________________________________________
2525void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
2526 InitDstartoD0pi2013pPb0100();
2527
a1f0af09 2528 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
2529
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);
44dd56b2 2532 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2533
44dd56b2 2534 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2535
a1f0af09 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);
2539
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);
44dd56b2 2543
2544};
2545//_________________________________________________________________________
2546void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
2547InitDstartoD0pi2013pPb0100();
2548
a1f0af09 2549 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
2550
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);
44dd56b2 2553 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2554
44dd56b2 2555 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2556
a1f0af09 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);
44dd56b2 2560
a1f0af09 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);
2564
44dd56b2 2565};
2566//_________________________________________________________________________
2567void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
2568 InitDstartoD0pi2013pPb0100();
a1f0af09 2569
2570 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
44dd56b2 2571
a1f0af09 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);
44dd56b2 2574 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
2575
44dd56b2 2576 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2577
a1f0af09 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);
2581
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);
44dd56b2 2585};
2586//_________________________________________________________________________
2587void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
2588 InitDstartoD0pi2013pPb0100();
a1f0af09 2589
2590 AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
44dd56b2 2591
a1f0af09 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);
44dd56b2 2594 for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
2595
44dd56b2 2596 for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
2597
a1f0af09 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);
44dd56b2 2601
a1f0af09 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);
2605
44dd56b2 2606};
3db225bf 2607
fa90243b 2608//_________________________________________________________________________
2609void AliHFSystErr::InitD0toKpi2013pPb020V0A(){
2610 InitD0toKpi2013pPb0100();
2611
2612 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2613
2614
2615 // Raw yield extraction
2616 fRawYield->SetBinContent(1,0.0);
751d48b6 2617 fRawYield->SetBinContent(2,0.085);
fa90243b 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);
751d48b6 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);
fa90243b 2623 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
2624
2625 // Cuts efficiency (from cuts variation)
2626 fCutsEff->SetBinContent(1,0.0);
2627 fCutsEff->SetBinContent(2,0.10);
751d48b6 2628 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
fa90243b 2629
2630 // PID efficiency (from PID/noPID)
2631 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2632
2633
2634}
2635
2636//_________________________________________________________________________
2637void AliHFSystErr::InitD0toKpi2013pPb2040V0A(){
2638 InitD0toKpi2013pPb0100();
2639
2640 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2641
2642
2643 // Raw yield extraction
2644 fRawYield->SetBinContent(1,0.0);
751d48b6 2645 fRawYield->SetBinContent(2,0.085);
2646 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
fa90243b 2647 for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 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);
2652
fa90243b 2653
2654 // Cuts efficiency (from cuts variation)
2655 fCutsEff->SetBinContent(1,0.0);
2656 fCutsEff->SetBinContent(2,0.10);
751d48b6 2657 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2658
fa90243b 2659
2660 // PID efficiency (from PID/noPID)
2661 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2662
2663}
2664
2665//_________________________________________________________________________
2666void AliHFSystErr::InitD0toKpi2013pPb4060V0A(){
2667 InitD0toKpi2013pPb0100();
2668
2669 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2670
2671
2672 // Raw yield extraction
2673 fRawYield->SetBinContent(1,0.0);
751d48b6 2674 fRawYield->SetBinContent(2,0.085);
fa90243b 2675 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 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);
fa90243b 2680 // for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
2681
2682 // Cuts efficiency (from cuts variation)
2683 fCutsEff->SetBinContent(1,0.0);
2684 fCutsEff->SetBinContent(2,0.10);
751d48b6 2685 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2686
fa90243b 2687
2688 // PID efficiency (from PID/noPID)
2689 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2690
2691}
2692
2693//_________________________________________________________________________
2694void AliHFSystErr::InitD0toKpi2013pPb60100V0A(){
2695 InitD0toKpi2013pPb0100();
2696
2697 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2698
2699 // Raw yield extraction
2700 fRawYield->SetBinContent(1,0.0);
751d48b6 2701 fRawYield->SetBinContent(2,0.085);
fa90243b 2702 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 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);
2707
fa90243b 2708
2709 // Cuts efficiency (from cuts variation)
2710 fCutsEff->SetBinContent(1,0.0);
751d48b6 2711 fCutsEff->SetBinContent(2,0.10);
2712 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2713
fa90243b 2714
2715 // PID efficiency (from PID/noPID)
751d48b6 2716 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 2717
2718}
2719
2720//_________________________________________________________________________
2721void AliHFSystErr::InitD0toKpi2013pPb020ZNA(){
2722 InitD0toKpi2013pPb0100();
2723
2724 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
2725
2726 // Raw yield extraction
2727 fRawYield->SetBinContent(1,0.0);
751d48b6 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);
fa90243b 2735
2736 // Cuts efficiency (from cuts variation)
2737 fCutsEff->SetBinContent(1,0.0);
751d48b6 2738 fCutsEff->SetBinContent(2,0.10);
2739 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2740
fa90243b 2741
2742 // PID efficiency (from PID/noPID)
2743 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2744
2745}
2746
2747//_________________________________________________________________________
2748void AliHFSystErr::InitD0toKpi2013pPb2040ZNA(){
2749 InitD0toKpi2013pPb0100();
2750
2751 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
2752
2753 // Raw yield extraction
2754 fRawYield->SetBinContent(1,0.0);
751d48b6 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);
2762
fa90243b 2763
2764 // Cuts efficiency (from cuts variation)
2765 fCutsEff->SetBinContent(1,0.0);
2766 fCutsEff->SetBinContent(2,0.10);
751d48b6 2767 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2768
fa90243b 2769 // PID efficiency (from PID/noPID)
2770 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2771
2772}
2773
2774//_________________________________________________________________________
2775void AliHFSystErr::InitD0toKpi2013pPb4060ZNA(){
2776 InitD0toKpi2013pPb0100();
2777
2778 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
2779
2780 // Raw yield extraction
2781 fRawYield->SetBinContent(1,0.0);
751d48b6 2782 fRawYield->SetBinContent(2,0.085);
fa90243b 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);
751d48b6 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);
fa90243b 2789
2790 // Cuts efficiency (from cuts variation)
2791 fCutsEff->SetBinContent(1,0.0);
751d48b6 2792 fCutsEff->SetBinContent(2,0.10);
2793 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2794
fa90243b 2795 // PID efficiency (from PID/noPID)
2796 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
2797
2798}
2799
2800//_________________________________________________________________________
2801void AliHFSystErr::InitD0toKpi2013pPb60100ZNA(){
2802 InitD0toKpi2013pPb0100();
2803
2804 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
2805
2806 // Raw yield extraction
2807 fRawYield->SetBinContent(1,0.0);
751d48b6 2808 fRawYield->SetBinContent(2,0.085);
fa90243b 2809 for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
751d48b6 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);
fa90243b 2815
2816 // Cuts efficiency (from cuts variation)
2817 fCutsEff->SetBinContent(1,0.0);
751d48b6 2818 fCutsEff->SetBinContent(2,0.10);
2819 for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
2820
fa90243b 2821 // PID efficiency (from PID/noPID)
751d48b6 2822 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 2823
2824}
2825
2826//_________________________________________________________________________
2827void AliHFSystErr::InitDstartoD0pi2013pPb020V0A(){
2828 InitDstartoD0pi2013pPb0100();
2829
2830 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
2831
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);
751d48b6 2841 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
fa90243b 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);
2844
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);
751d48b6 2854 for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2855 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
2856
2857 // PID efficiency (from PID/noPID)
751d48b6 2858 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
fa90243b 2859
2860
2861}
2862
2863//_________________________________________________________________________
2864void AliHFSystErr::InitDstartoD0pi2013pPb2040V0A(){
2865 InitDstartoD0pi2013pPb0100();
2866
2867 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
2868
2869 // Raw yield extraction
2870 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2871 fRawYield->SetBinContent(1,0.0);
2872 fRawYield->SetBinContent(2,0.10);
fa90243b 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);
751d48b6 2879 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
fa90243b 2880 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
751d48b6 2881 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
fa90243b 2882
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);
751d48b6 2894 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 2895 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2896
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);
2900
2901 // MC dN/dpt
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);*/
2906
2907 // particle-antiparticle
2908 /*
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%
2914 */
2915
2916}
2917
2918//_________________________________________________________________________
2919void AliHFSystErr::InitDstartoD0pi2013pPb4060V0A(){
2920 InitDstartoD0pi2013pPb0100();
2921
2922 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
2923
2924 // Raw yield extraction
2925 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2926 fRawYield->SetBinContent(1,0.0);
2927 fRawYield->SetBinContent(2,0.0);
2928 fRawYield->SetBinContent(3,0.06);
2929 fRawYield->SetBinContent(4,0.06);
fa90243b 2930 fRawYield->SetBinContent(5,0.05);
2931 fRawYield->SetBinContent(6,0.05);
2932 fRawYield->SetBinContent(7,0.05);
2933 fRawYield->SetBinContent(8,0.05);
751d48b6 2934 for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
fa90243b 2935 for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
751d48b6 2936 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
fa90243b 2937
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);
751d48b6 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);
fa90243b 2950 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
2951
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);
2955
2956 // MC dN/dpt
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);*/
2961
2962 // particle-antiparticle
2963 /*
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%
2969 */
2970
2971
2972}
2973
2974//_________________________________________________________________________
2975void AliHFSystErr::InitDstartoD0pi2013pPb60100V0A(){
2976 InitDstartoD0pi2013pPb0100();
2977
2978 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
2979
2980 // Raw yield extraction
2981 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 2982 fRawYield->SetBinContent(1,0.0);
2983 fRawYield->SetBinContent(2,0.0);
fa90243b 2984 fRawYield->SetBinContent(3,0.06);
2985 fRawYield->SetBinContent(4,0.06);
751d48b6 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);
fa90243b 2992 for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
2993
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.);
751d48b6 2998 fCutsEff->SetBinContent(3,0.06);
2999 fCutsEff->SetBinContent(4,0.06);
fa90243b 3000 fCutsEff->SetBinContent(5,0.06);
3001 fCutsEff->SetBinContent(6,0.06);
3002 fCutsEff->SetBinContent(7,0.06);
3003 fCutsEff->SetBinContent(8,0.06);
751d48b6 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);
fa90243b 3006 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3007
3008 // PID efficiency (from PID/noPID)
3009 //fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 3010 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 3011
3012 // MC dN/dpt
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);*/
3017
3018 // particle-antiparticle
3019 /*
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%
3025 */
3026
3027
3028}
3029
3030//_________________________________________________________________________
3031void AliHFSystErr::InitDstartoD0pi2013pPb020ZNA(){
3032 InitDstartoD0pi2013pPb0100();
3033
3034 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
3035
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);
751d48b6 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);
fa90243b 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.);
3049
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);
751d48b6 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);
fa90243b 3062 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3063
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.);
3067
3068 // MC dN/dpt
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);*/
3073
3074 // particle-antiparticle
3075 /*
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%
3081 */
3082
3083
3084}
3085
3086//_________________________________________________________________________
3087void AliHFSystErr::InitDstartoD0pi2013pPb2040ZNA(){
3088 InitDstartoD0pi2013pPb0100();
3089
3090 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
3091
3092 // Raw yield extraction
3093 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 3094 fRawYield->SetBinContent(1,0.0);
3095 fRawYield->SetBinContent(2,0.0);
3096 fRawYield->SetBinContent(3,0.06);
3097 fRawYield->SetBinContent(4,0.06);
fa90243b 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.);
3105
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);
751d48b6 3117 for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
fa90243b 3118 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3119
3120 // PID efficiency (from PID/noPID)
3121 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 3122 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 3123
3124 // MC dN/dpt
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);*/
3129
3130 // particle-antiparticle
3131 /*
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%
3137 */
3138
3139
3140}
3141
3142//_________________________________________________________________________
3143void AliHFSystErr::InitDstartoD0pi2013pPb4060ZNA(){
3144 InitDstartoD0pi2013pPb0100();
3145
3146 AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
3147
3148 // Raw yield extraction
3149 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 3150 fRawYield->SetBinContent(1,0.0);
3151 fRawYield->SetBinContent(2,0.0);
3152 fRawYield->SetBinContent(3,0.06);
3153 fRawYield->SetBinContent(4,0.06);
fa90243b 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.);
3161
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.);
751d48b6 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);
fa90243b 3174 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3175
3176 // PID efficiency (from PID/noPID)
3177 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 3178 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 3179
3180 // MC dN/dpt
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);*/
3185
3186 // particle-antiparticle
3187 /*
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%
3193 */
3194
3195
3196
3197}
3198
3199//_________________________________________________________________________
3200void AliHFSystErr::InitDstartoD0pi2013pPb60100ZNA(){
3201 InitDstartoD0pi2013pPb0100();
3202
3203 AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
3204
3205 // Raw yield extraction
3206 //fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
751d48b6 3207 fRawYield->SetBinContent(1,0.);
3208 fRawYield->SetBinContent(2,0.);
fa90243b 3209 fRawYield->SetBinContent(3,0.06);
3210 fRawYield->SetBinContent(4,0.06);
3211 fRawYield->SetBinContent(5,0.05);
3212 fRawYield->SetBinContent(6,0.05);
751d48b6 3213 fRawYield->SetBinContent(7,0.05);
3214 fRawYield->SetBinContent(8,0.05);
fa90243b 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.);
3218
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.);
751d48b6 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);
fa90243b 3231 for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
3232
3233 // PID efficiency (from PID/noPID)
3234 // fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
751d48b6 3235 for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fa90243b 3236
3237 // MC dN/dpt
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);*/
3242
3243 // particle-antiparticle
3244 /*
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%
3250 */
3251
3252}
3253
3254
3255
414bc6a5 3256
9ca30701 3257//--------------------------------------------------------------------------
3258void AliHFSystErr::InitLctopKpi2010pp() {
3259 //
3260 // Lc->pKpi syst errors. Responsible: R. Romita
3261 // 2010 pp sample
3262 //
3263
3264 // Normalization
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
3267
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)
3271
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);
3278
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);
3284
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%
3291
3292 // MC dN/dpt
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);
3299}
3300
53c91ed7 3301//--------------------------------------------------------------------------
3302void AliHFSystErr::InitLctopK0S2010pp() {
3303 //
3304 // Lc->pK0S syst errors. Responsible: A. De Caro
3305 // 2010 pp sample
3306 // Only tool preparation: no physical mean for these values!
3307 //
3308
3309 // Normalization
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
3312
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
3316
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
3321
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
3324
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
3328
3329 // MC dN/dpt
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
3333
3334}
3335
5541b811 3336//--------------------------------------------------------------------------
0233abe6 3337Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
3338 //
3339 // Get error
3340 //
3341
3342 Int_t bin=fCutsEff->FindBin(pt);
3343
3344 return fCutsEff->GetBinContent(bin);
3345}
3346//--------------------------------------------------------------------------
3347Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
3348 //
3349 // Get error
3350 //
3351
3352 Int_t bin=fMCPtShape->FindBin(pt);
3353
3354 return fMCPtShape->GetBinContent(bin);
3355}
3356//--------------------------------------------------------------------------
3357Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
3358 //
3359 // Get error
3360 //
3361
3362 Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
3363
3364 return TMath::Sqrt(err);
3365}
3366//--------------------------------------------------------------------------
3367Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
3368 //
3369 // Get error
3370 //
3371
3372 Int_t bin=fPIDEff->FindBin(pt);
3373
3374 return fPIDEff->GetBinContent(bin);
3375}
3376//--------------------------------------------------------------------------
3377Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
3378 //
3379 // Get error
3380 //
3381
3382 Int_t bin=fTrackingEff->FindBin(pt);
3383
3384 return fTrackingEff->GetBinContent(bin);
3385}
3386//--------------------------------------------------------------------------
3387Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
3388 //
3389 // Get error
3390 //
3391
3392 Int_t bin=fRawYield->FindBin(pt);
3393
3394 return fRawYield->GetBinContent(bin);
3395}
3396//--------------------------------------------------------------------------
3397Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
3398 //
3399 // Get error
3400 //
3401
3402 Int_t bin=fPartAntipart->FindBin(pt);
3403
3404 return fPartAntipart->GetBinContent(bin);
3405}
3406//--------------------------------------------------------------------------
3407Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
3408 //
3409 // Get total syst error (except norm. error)
3410 //
3411
3412 Double_t err=0.;
3413
3414 if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
3415 if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
f9a85c82 3416 // if(fBR) err += GetBRErr()*GetBRErr();
0233abe6 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);
3421
3422 err += feeddownErr*feeddownErr;
3423
3424 return TMath::Sqrt(err);
3425}
3426//---------------------------------------------------------------------------
5ad8ad40 3427void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
0233abe6 3428 //
3429 // Draw errors
3430 //
3431 gStyle->SetOptStat(0);
3432
a07118d0 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);
0233abe6 3436 cSystErr->SetFillColor(0);
3437
436b185f 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");
a07118d0 3440 hFrame->SetAxisRange(-0.5,0.5,"Y");
0233abe6 3441 hFrame->Draw();
3442
a07118d0 3443 TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
3444 leg->SetTextSize(0.03601695);
0233abe6 3445 leg->SetFillStyle(0);
3446 leg->SetBorderSize(0);
2ee3afe2 3447
436b185f 3448 TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
2ee3afe2 3449 Int_t nbins = fNorm->GetNbinsX();
3450 TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
436b185f 3451 for(Int_t i=1;i<=36;i++) {
2ee3afe2 3452 Double_t pt = hTotErr->GetBinCenter(i);
3453 Double_t ptwidth = hTotErr->GetBinWidth(i);
3454
3455 if(grErrFeeddown) {
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);
a07118d0 3462 if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
2ee3afe2 3463 erryh = grErrFeeddown->GetErrorYhigh(j);
3464 erryl = grErrFeeddown->GetErrorYlow(j);
3465 }
3466 }
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);
3471
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
3475 }
3476 else {
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
3480 }
3481
3482 }
f968ef30 3483 gTotErr->SetLineColor(kBlack);
3484 gTotErr->SetFillColor(kRed);
3485 gTotErr->SetFillStyle(3002);
2ee3afe2 3486 gTotErr->Draw("2");
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");
3492
0233abe6 3493
e373b9f0 3494 fNorm->SetFillColor(1);
3495 fNorm->SetFillStyle(3002);
3496 //fNorm->Draw("same");
3497 //TH1F *hNormRefl = ReflectHisto(fNorm);
3498 //hNormRefl->Draw("same");
27360a0e 3499 Double_t norm = fNorm->GetBinContent(1)*100;
3500 leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
e373b9f0 3501
5ad8ad40 3502 if(grErrFeeddown) {
3503 grErrFeeddown->SetFillColor(kTeal-8);
3504 grErrFeeddown->SetFillStyle(3001);
3505 grErrFeeddown->Draw("2");
3506 leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
3507 }
0233abe6 3508 if(fTrackingEff) {
5ad8ad40 3509 fTrackingEff->SetFillColor(4);
3510 fTrackingEff->SetFillStyle(3006);
0233abe6 3511 fTrackingEff->Draw("same");
f968ef30 3512 TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
3513 hTrackingEffRefl->Draw("same");
0233abe6 3514 leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
3515 }
3516 if(fBR) {
3517 fBR->SetFillColor(6);
3518 fBR->SetFillStyle(3005);
5ad8ad40 3519 //fBR->SetFillStyle(3020);
0233abe6 3520 fBR->Draw("same");
f968ef30 3521 TH1F *hBRRefl = ReflectHisto(fBR);
3522 hBRRefl->Draw("same");
0233abe6 3523 leg->AddEntry(fBR,"Branching ratio","f");
3524 }
3525 if(fRawYield) {
5ad8ad40 3526 Int_t ci; // for color index setting
3527 ci = TColor::GetColor("#00cc00");
3528 fRawYield->SetLineColor(ci);
3529 // fRawYield->SetLineColor(3);
0233abe6 3530 fRawYield->SetLineWidth(3);
3531 fRawYield->Draw("same");
f968ef30 3532 TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
3533 hRawYieldRefl->Draw("same");
844405bc 3534 leg->AddEntry(fRawYield,"Yield extraction","l");
0233abe6 3535 }
3536 if(fCutsEff) {
3537 fCutsEff->SetLineColor(4);
3538 fCutsEff->SetLineWidth(3);
3539 fCutsEff->Draw("same");
f968ef30 3540 TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
3541 hCutsEffRefl->Draw("same");
0233abe6 3542 leg->AddEntry(fCutsEff,"Cuts efficiency","l");
3543 }
3544 if(fPIDEff) {
3545 fPIDEff->SetLineColor(7);
3546 fPIDEff->SetLineWidth(3);
3547 fPIDEff->Draw("same");
f968ef30 3548 TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
3549 hPIDEffRefl->Draw("same");
0233abe6 3550 leg->AddEntry(fPIDEff,"PID efficiency","l");
3551 }
3552 if(fMCPtShape) {
5ad8ad40 3553 Int_t ci = TColor::GetColor("#9933ff");
3554 fMCPtShape->SetLineColor(ci);
3555 // fMCPtShape->SetLineColor(8);
0233abe6 3556 fMCPtShape->SetLineWidth(3);
3557 fMCPtShape->Draw("same");
f968ef30 3558 TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
3559 hMCPtShapeRefl->Draw("same");
0233abe6 3560 leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
3561 }
3562 if(fPartAntipart) {
5ad8ad40 3563 Int_t ci = TColor::GetColor("#ff6600");
3564 fPartAntipart->SetLineColor(ci);
3565 // fPartAntipart->SetLineColor(9);
0233abe6 3566 fPartAntipart->SetLineWidth(3);
3567 fPartAntipart->Draw("same");
f968ef30 3568 TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
3569 hPartAntipartRefl->Draw("same");
0233abe6 3570 leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
3571 }
0233abe6 3572
0233abe6 3573
3574 leg->Draw();
3575
5ad8ad40 3576 cSystErr->SaveAs("RelativeSystematics.eps");
3577
0233abe6 3578 return;
3579}
f968ef30 3580//-------------------------------------------------------------------------
3581TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
3582 //
3583 // Clones and reflects histogram
3584 //
3585 TH1F *hout=(TH1F*)hin->Clone("hout");
3586 hout->Scale(-1.);
3587
3588 return hout;
3589}
9a8b5080 3590
ada850b0 3591
3592