#include <TStyle.h>
#include <TGraphAsymmErrors.h>
#include <TMath.h>
#include <TCanvas.h>
#include <TH2F.h>
#include <TLegend.h>
#include <TColor.h>
#include "AliLog.h"
#include "AliHFSystErr.h"
ClassImp(AliHFSystErr)
AliHFSystErr::AliHFSystErr(const Char_t* name, const Char_t* title) :
TNamed(name,title),
fNorm(0),
fRawYield(0),
fTrackingEff(0),
fBR(0),
fCutsEff(0),
fPIDEff(0),
fMCPtShape(0),
fPartAntipart(0),
fRunNumber(10),
fCollisionType(0),
fCentralityClass("0100"),
fRapidityRange("0101"),
fIsLowEnergy(false),
fIsCentScan(false),
fIsRapidityScan(false)
{
}
AliHFSystErr::~AliHFSystErr() {
}
void AliHFSystErr::Init(Int_t decay){
if ((fRunNumber>11) && fIsLowEnergy==false) {
AliFatal("Only settings for 2010 and the low energy runs are implemented so far");
}
switch(decay) {
case 1:
if (fCollisionType==0) {
if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
else InitD0toKpi2010pp();
}
else if (fCollisionType==1) {
if (fRunNumber == 10){
if (fCentralityClass=="010") InitD0toKpi2010PbPb010CentScan();
else if (fCentralityClass=="1020") InitD0toKpi2010PbPb1020CentScan();
else if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
else if (fCentralityClass=="2040") InitD0toKpi2010PbPb2040CentScan();
else if (fCentralityClass=="4060") InitD0toKpi2010PbPb4060CentScan();
else if (fCentralityClass=="6080") InitD0toKpi2010PbPb6080CentScan();
else if (fCentralityClass=="4080") InitD0toKpi2010PbPb4080();
else AliFatal("Not yet implemented");
}
else if (fRunNumber == 11){
if (fCentralityClass=="07half") InitD0toKpi2011PbPb07half();
else if (fCentralityClass=="010") InitD0toKpi2011PbPb010();
else if (fCentralityClass=="3050InPlane") InitD0toKpi2011PbPb3050InPlane();
else if (fCentralityClass=="3050OutOfPlane")InitD0toKpi2011PbPb3050OutOfPlane();
else if (fCentralityClass == "3050")InitD0toKpi2011PbPb3050();
else if (fCentralityClass=="010" && fIsCentScan) InitD0toKpi2011PbPb010CentScan();
else if (fCentralityClass=="1020") InitD0toKpi2011PbPb1020CentScan();
else if (fCentralityClass=="2030") InitD0toKpi2011PbPb2030CentScan();
else if (fCentralityClass=="3040") InitD0toKpi2011PbPb3040CentScan();
else if (fCentralityClass=="4050") InitD0toKpi2011PbPb4050CentScan();
else if (fCentralityClass=="5080") InitD0toKpi2010PbPb5080CentScan();
else AliFatal("Not yet implemented");
}
}
else if (fCollisionType==2) {
if (fCentralityClass=="0100") InitD0toKpi2013pPb0100();
if (fCentralityClass=="020V0A") InitD0toKpi2013pPb020V0A();
if (fCentralityClass=="2040V0A") InitD0toKpi2013pPb2040V0A();
if (fCentralityClass=="4060V0A") InitD0toKpi2013pPb4060V0A();
if (fCentralityClass=="60100V0A") InitD0toKpi2013pPb60100V0A();
if (fCentralityClass=="020ZNA") InitD0toKpi2013pPb020ZNA();
if (fCentralityClass=="2040ZNA") InitD0toKpi2013pPb2040ZNA();
if (fCentralityClass=="4060ZNA") InitD0toKpi2013pPb4060ZNA();
if (fCentralityClass=="60100ZNA") InitD0toKpi2013pPb60100ZNA();
if (fIsRapidityScan) {
if (fRapidityRange == "0804") InitD0toKpi2013pPb0100RapScan0804();
if (fRapidityRange == "0401") InitD0toKpi2013pPb0100RapScan0401();
if (fRapidityRange == "0101") InitD0toKpi2013pPb0100RapScan0101();
if (fRapidityRange == "0104") InitD0toKpi2013pPb0100RapScan0104();
if (fRapidityRange == "0408") InitD0toKpi2013pPb0100RapScan0408();
}
}
else AliFatal("Not yet implemented");
break;
case 2:
if (fCollisionType==0) {
if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
else InitDplustoKpipi2010pp();
}
else if (fCollisionType==1) {
if (fRunNumber == 10){
if (fCentralityClass=="010") InitDplustoKpipi2010PbPb010CentScan();
else if (fCentralityClass=="1020") InitDplustoKpipi2010PbPb1020CentScan();
else if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
else if (fCentralityClass=="2040") InitDplustoKpipi2010PbPb2040CentScan();
else if (fCentralityClass=="4060") InitDplustoKpipi2010PbPb4060CentScan();
else if (fCentralityClass=="6080") InitDplustoKpipi2010PbPb6080CentScan();
else if (fCentralityClass=="4080") InitDplustoKpipi2010PbPb4080();
else AliFatal("Not yet implemented");
}
if(fRunNumber == 11){
if (fCentralityClass=="07half") InitDplustoKpipi2011PbPb07half();
else if (fCentralityClass=="010") InitDplustoKpipi2011PbPb010();
else if (fCentralityClass=="010" && fIsCentScan) InitDplustoKpipi2011PbPb010CentScan();
else if (fCentralityClass=="1020") InitDplustoKpipi2011PbPb1020CentScan();
else if (fCentralityClass=="2030") InitDplustoKpipi2011PbPb2030CentScan();
else if (fCentralityClass=="3040") InitDplustoKpipi2011PbPb3040CentScan();
else if (fCentralityClass=="4050") InitDplustoKpipi2011PbPb4050CentScan();
else if (fCentralityClass=="5080") InitDplustoKpipi2010PbPb5080CentScan();
else if (fCentralityClass=="3050") InitDplustoKpipi2011PbPb3050();
else AliFatal("Not yet implemented");
}
}
else if (fCollisionType==2) {
if (fCentralityClass=="0100") InitDplustoKpipi2013pPb0100();
if (fIsRapidityScan) {
if (fRapidityRange == "0804") InitDplustoKpipi2013pPb0100RapScan0804();
if (fRapidityRange == "0401") InitDplustoKpipi2013pPb0100RapScan0401();
if (fRapidityRange == "0101") InitDplustoKpipi2013pPb0100RapScan0101();
if (fRapidityRange == "0104") InitDplustoKpipi2013pPb0100RapScan0104();
if (fRapidityRange == "0408") InitDplustoKpipi2013pPb0100RapScan0408();
}
}
else AliFatal("Not yet implemented");
break;
case 3:
if (fCollisionType==0) {
if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
else InitDstartoD0pi2010pp();
}
else if (fCollisionType==1) {
if (fRunNumber == 10){
if (fCentralityClass=="010") InitDstartoD0pi2010PbPb010CentScan();
else if (fCentralityClass=="1020") InitDstartoD0pi2010PbPb1020CentScan();
else if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
else if (fCentralityClass=="2040" && fIsCentScan) InitDstartoD0pi2010PbPb2040CentScan();
else if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
else if (fCentralityClass=="4060") InitDstartoD0pi2010PbPb4060CentScan();
else if (fCentralityClass=="6080") InitDstartoD0pi2010PbPb6080CentScan();
else if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
else AliFatal("Not yet implemented");
}
if (fRunNumber == 11){
if (fCentralityClass=="07half") InitDstartoD0pi2011PbPb07half();
else if (fCentralityClass=="010") InitDstartoD0pi2011PbPb010();
else if (fCentralityClass=="010" && fIsCentScan) InitDstartoD0pi2011PbPb010CentScan();
else if (fCentralityClass=="1020") InitDstartoD0pi2011PbPb1020CentScan();
else if (fCentralityClass=="2030") InitDstartoD0pi2011PbPb2030CentScan();
else if (fCentralityClass=="3040") InitDstartoD0pi2011PbPb3040CentScan();
else if (fCentralityClass=="4050") InitDstartoD0pi2011PbPb4050CentScan();
else if (fCentralityClass=="5080") InitDstartoD0pi2010PbPb5080CentScan();
else if (fCentralityClass=="3050") InitDstartoD0pi2011PbPb3050();
else AliFatal("Not yet implemented");
}
}
else if (fCollisionType==2) {
if (fCentralityClass=="0100") InitDstartoD0pi2013pPb0100();
if (fCentralityClass=="020V0A") InitDstartoD0pi2013pPb020V0A();
if (fCentralityClass=="2040V0A") InitDstartoD0pi2013pPb2040V0A();
if (fCentralityClass=="4060V0A") InitDstartoD0pi2013pPb4060V0A();
if (fCentralityClass=="60100V0A") InitDstartoD0pi2013pPb60100V0A();
if (fCentralityClass=="020ZNA") InitDstartoD0pi2013pPb020ZNA();
if (fCentralityClass=="2040ZNA") InitDstartoD0pi2013pPb2040ZNA();
if (fCentralityClass=="4060ZNA") InitDstartoD0pi2013pPb4060ZNA();
if (fCentralityClass=="60100ZNA") InitDstartoD0pi2013pPb60100ZNA();
if (fIsRapidityScan) {
if (fRapidityRange == "0804") InitDstartoD0pi2013pPb0100RapScan0804();
if (fRapidityRange == "0401") InitDstartoD0pi2013pPb0100RapScan0401();
if (fRapidityRange == "0101") InitDstartoD0pi2013pPb0100RapScan0101();
if (fRapidityRange == "0104") InitDstartoD0pi2013pPb0100RapScan0104();
if (fRapidityRange == "0408") InitDstartoD0pi2013pPb0100RapScan0408();
}
}
else AliFatal("Not yet implemented");
break;
case 4:
if (fCollisionType==0) InitDstoKKpi2010pp();
else if (fCollisionType==1) {
if (fCentralityClass=="07half") InitDstoKKpi2011PbPb07half();
else AliFatal("Not yet implemented");
}
else if (fCollisionType==2) {
if (fCentralityClass=="0100") InitDstoKKpi2013pPb0100();
}
else AliFatal("Not yet implemented");
break;
case 5:
if (fCollisionType==0) InitLctopKpi2010pp();
else AliFatal("Not yet implemented");
break;
case 6:
if (fCollisionType==0) InitLctopK0S2010pp();
else AliFatal("Not yet implemented");
break;
default:
printf("Invalid decay type: %d\n",decay);
break;
}
}
void AliHFSystErr::InitD0toKpi2010pp() {
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.08);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,1);
fRawYield->SetBinContent(2,0.22);
fRawYield->SetBinContent(3,0.1);
for(Int_t i=4;i<=7;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=8;i<=12;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,1);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.03);
fPIDEff->SetBinContent(2,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
fMCPtShape->SetBinContent(1,0.03);
fMCPtShape->SetBinContent(2,0.03);
return;
}
void AliHFSystErr::InitD0toKpi2010PbPb020() {
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
fRawYield->SetBinContent(1,0);
fRawYield->SetBinContent(2,0);
fRawYield->SetBinContent(3,0.08);
for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.13);
fCutsEff->SetBinContent(4,0.11);
for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(3,0.04);
fMCPtShape->SetBinContent(4,0.02);
for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
return;
}
void AliHFSystErr::InitD0toKpi2011PbPb07half() {
fNorm = new TH1F("fNorm","fNorm",36,0,36);
for(Int_t i=1;i<36;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",36,0,36);
for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.012);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);
for(Int_t i=25;i<=36;i++) fTrackingEff->SetBinContent(i,0.);
fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
fRawYield->SetBinContent(1,0);
fRawYield->SetBinContent(2,0.1);
fRawYield->SetBinContent(3,0.05);
for(Int_t i=4;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.30);
for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.);
fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.15);
fCutsEff->SetBinContent(3,0.13);
fCutsEff->SetBinContent(4,0.11);
fCutsEff->SetBinContent(5,0.08);
for(Int_t i=6;i<=24;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);
fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
for(Int_t i=25;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(2,0.06);
fMCPtShape->SetBinContent(3,0.04);
fMCPtShape->SetBinContent(4,0.02);
for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=17;i<=24;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=25;i<=36;i++) fMCPtShape->SetBinContent(i,0.0);
return;
}
void AliHFSystErr::InitD0toKpi2011PbPb3050InPlane() {
InitD0toKpi2011PbPb07half();
fRawYield->SetBinContent(1,0);
fRawYield->SetBinContent(2,0);
fRawYield->SetBinContent(3,0.05);
fRawYield->SetBinContent(4,0.08);
for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.0);
fCutsEff->SetBinContent(3,0.10);
fCutsEff->SetBinContent(4,0.10);
fCutsEff->SetBinContent(5,0.10);
fCutsEff->SetBinContent(6,0.10);
for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=25;i<=36;i++) fCutsEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitD0toKpi2011PbPb3050OutOfPlane() {
InitD0toKpi2011PbPb07half();
fRawYield->SetBinContent(1,0);
fRawYield->SetBinContent(2,0.);
fRawYield->SetBinContent(3,0.05);
for(Int_t i=4;i<=6;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.15);
for(Int_t i=17;i<=36;i++) fRawYield->SetBinContent(i,0.);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.10);
fCutsEff->SetBinContent(4,0.10);
fCutsEff->SetBinContent(5,0.10);
fCutsEff->SetBinContent(6,0.10);
for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=17;i<=36;i++) fCutsEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitD0toKpi2010PbPb4080() {
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
fRawYield->SetBinContent(1,0);
fRawYield->SetBinContent(2,0);
for(Int_t i=3;i<=16;i++) fRawYield->SetBinContent(i,0.05);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.13);
fCutsEff->SetBinContent(4,0.11);
for(Int_t i=5;i<=16;i++) fCutsEff->SetBinContent(i,0.10);
for(Int_t i=17;i<=20;i++) fCutsEff->SetBinContent(i,0.);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=3;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
return;
}
void AliHFSystErr::InitD0toKpi2011PbPb010() {
AliInfo(" Settings for D0 --> Kpi Pb-Pb collisions at 2.76 TeV - 0-10 centrality - newww!!");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.012);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.10);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,0);
fRawYield->SetBinContent(2,0.15);
fRawYield->SetBinContent(3,0.09);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.06);
fRawYield->SetBinContent(6,0.06);
for (Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.06);
for (Int_t i=9; i<=12; i++)fRawYield->SetBinContent(i,0.07);
for (Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.11);
for (Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.30);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.15);
fCutsEff->SetBinContent(3,0.12);
fCutsEff->SetBinContent(4,0.05);
fCutsEff->SetBinContent(5,0.05);
fCutsEff->SetBinContent(6,0.05);
for (Int_t i=7; i<=8; i++) fCutsEff->SetBinContent(i,0.05);
for (Int_t i=9; i<=12; i++)fCutsEff->SetBinContent(i,0.05);
for (Int_t i=13; i<=16; i++) fCutsEff->SetBinContent(i,0.05);
for (Int_t i=17; i<=24; i++) fCutsEff->SetBinContent(i,0.05);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
fMCPtShape->SetBinContent(1,0.);
fMCPtShape->SetBinContent(2,0.15);
fMCPtShape->SetBinContent(3,0.12);
fMCPtShape->SetBinContent(4,0.02);
for(Int_t i=5;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
return;
}
void AliHFSystErr::InitD0toKpi2011PbPb3050() {
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
fRawYield->SetBinContent(1,0);
fRawYield->SetBinContent(2,0.102);
fRawYield->SetBinContent(3,0.071);
fRawYield->SetBinContent(4,0.052);
fRawYield->SetBinContent(5,0.063);
fRawYield->SetBinContent(6,0.071);
fRawYield->SetBinContent(7,0.071);
fRawYield->SetBinContent(8,0.071);
fRawYield->SetBinContent(9,0.071);
fRawYield->SetBinContent(10,0.071);
fRawYield->SetBinContent(11,0.071);
fRawYield->SetBinContent(12,0.071);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.112);
for(Int_t i=17;i<=20;i++) fRawYield->SetBinContent(i,0);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.10);
fCutsEff->SetBinContent(3,0.10);
fCutsEff->SetBinContent(4,0.10);
fCutsEff->SetBinContent(5,0.10);
fCutsEff->SetBinContent(6,0.10);
for(Int_t i=7;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.15);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=2;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",0,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(2,0.06);
fMCPtShape->SetBinContent(3,0.04);
fMCPtShape->SetBinContent(4,0.02);
for(Int_t i=13;i<=16;i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=17;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
return;
}
void AliHFSystErr::InitD0toKpi2010ppLowEn() {
AliInfo(" Settings for D0 --> K pi, p-p collisions at 2.76 TeV");
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.012);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.10);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
fRawYield->SetBinContent(1,1);
for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,0.15);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
fCutsEff->SetBinContent(2,0.20);
for(Int_t i=7;i<=20;i++) fCutsEff->SetBinContent(i,0.15);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.15);
for(Int_t i=7;i<=20;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(1,0.03);
fMCPtShape->SetBinContent(2,0.03);
return;
}
void AliHFSystErr::InitD0toKpi2013pPb0100(){
AliInfo(" Settings for D0 --> K pi, p-Pb collisions at 5.023 TeV");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.013);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.06);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.08);
for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.11);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.08);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
fMCPtShape->SetBinContent(1,0.0);
fMCPtShape->SetBinContent(2,0.02);
for(Int_t i=3;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
return;
}
void AliHFSystErr::InitDplustoKpipi2010pp() {
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.12);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,1);
fRawYield->SetBinContent(2,0.25);
fRawYield->SetBinContent(3,0.25);
fRawYield->SetBinContent(4,0.20);
fRawYield->SetBinContent(5,0.09);
fRawYield->SetBinContent(6,0.09);
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=12;i<=24;i++) fRawYield->SetBinContent(i,0.10);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
fPIDEff->SetBinContent(1,0.15);
fPIDEff->SetBinContent(2,0.15);
fPIDEff->SetBinContent(3,0.15);
fPIDEff->SetBinContent(4,0.15);
for(Int_t i=12;i<=16;i++) fPIDEff->SetBinContent(i,0.10);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
fMCPtShape->SetBinContent(1,0.03);
fMCPtShape->SetBinContent(2,0.03);
return;
}
void AliHFSystErr::InitDstoKKpi2010pp() {
fNorm = new TH1F("fNorm","fNorm",12,0,12);
for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.035);
fBR = new TH1F("fBR","fBR",12,0,12);
for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.06);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.12);
fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
fRawYield->SetBinContent(1,1);
fRawYield->SetBinContent(2,1);
fRawYield->SetBinContent(3,0.20);
fRawYield->SetBinContent(4,0.20);
fRawYield->SetBinContent(5,0.15);
fRawYield->SetBinContent(6,0.15);
fRawYield->SetBinContent(7,0.15);
fRawYield->SetBinContent(8,0.15);
fRawYield->SetBinContent(9,0.20);
fRawYield->SetBinContent(10,0.20);
fRawYield->SetBinContent(11,0.20);
fRawYield->SetBinContent(12,0.20);
fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.15);
fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.07);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
return;
}
void AliHFSystErr::InitDplustoKpipi2011PbPb07half() {
fNorm = new TH1F("fNorm","fNorm",36,0,36);
for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",36,0,36);
for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,.10);
fRawYield->SetBinContent(4,0.30);
fRawYield->SetBinContent(5,0.20);
fRawYield->SetBinContent(6,0.20);
for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.20);
fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
for(Int_t i=13; i<=36; i++) fCutsEff->SetBinContent(i,0.08);
fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
for(Int_t iBin=1; iBin<=5; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
for(Int_t iBin=6; iBin<=36; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
return;
}
void AliHFSystErr::InitDplustoKpipi2011PbPb010() {
fNorm = new TH1F("fNorm","fNorm",36,0,36);
for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",36,0,36);
for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);
for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i,0.);
fRawYield->SetBinContent(4,0.10);
fRawYield->SetBinContent(5,0.08);
fRawYield->SetBinContent(6,0.08);
for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=17; i<=24; i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=25; i<=36; i++) fRawYield->SetBinContent(i,0.08);
fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
for(Int_t i=1; i<=36; i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.);
fMCPtShape->SetBinContent(4,0.06);
fMCPtShape->SetBinContent(5,0.02);
for (Int_t i=6; i<=36; i++) fMCPtShape->SetBinContent(i,0.01);
return;
}
void AliHFSystErr::InitDplustoKpipi2011PbPb3050() {
AliInfo(" Settings for D+ --> K pi pi, Pb-Pb collisions for 30-50 CC");
fNorm = new TH1F("fNorm","fNorm",16,0,16);
for(Int_t i=1;i<=16;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",16,0,16);
for(Int_t i=1;i<=16;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",16,0,16);
for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",16,0,16);
for(Int_t i=1;i<=12;i++) fRawYield->SetBinContent(i,0.09);
for(Int_t i=13; i<=16; i++) fRawYield->SetBinContent(i,0.12);
fCutsEff = new TH1F("fCutsEff","fCutsEff",16,0,16);
for(Int_t i=1; i<=12; i++) fCutsEff->SetBinContent(i,0.10);
for(Int_t i=12; i<=16; i++) fCutsEff->SetBinContent(i,0.14);
fPIDEff = new TH1F("fPIDEff","fPIDEff",16,0,16);
for(Int_t i=1;i<=16;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",16,0,16);
fMCPtShape->SetBinContent(3,0.10);
fMCPtShape->SetBinContent(4,0.04);
for(Int_t iBin=5; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.02);
return;
}
void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
fNorm = new TH1F("fNorm","fNorm",12,0,12);
for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",12,0,12);
for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,0.053);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",12,0,12);
for(Int_t i=1;i<=6;i++) fRawYield->SetBinContent(i,.30);
for(Int_t i=7; i<=12; i++) fRawYield->SetBinContent(i,0.20);
fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20);
fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
return;
}
void AliHFSystErr::InitDstoKKpi2013pPb0100(){
AliInfo(" Settings for Ds --> K K pi, p-Pb collisions at 5.023 TeV");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=3;i<=12;i++) fNorm->SetBinContent(i,0.00);
fBR = new TH1F("fBR","fBR",12,0,12);
for(Int_t i=3;i<=12;i++) fBR->SetBinContent(i,0.053);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
for(Int_t i=3;i<=12;i++) fTrackingEff->SetBinContent(i,0.09);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
for(Int_t i=3; i<=4; i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=5; i<=6; i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=7; i<=8; i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9; i<=12; i++) fRawYield->SetBinContent(i,0.08);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
for(Int_t i=3;i<=12;i++) fCutsEff->SetBinContent(i,0.15);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=3;i<=4;i++) fPIDEff->SetBinContent(i,0.10);
for(Int_t i=5;i<=12;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=3;i<=12;i++) fMCPtShape->SetBinContent(i,0.04);
return;
}
void AliHFSystErr::InitDplustoKpipi2010PbPb020() {
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
for(Int_t i=1;i<=20;i++) fRawYield->SetBinContent(i,.10);
fRawYield->SetBinContent(7,0.20);
fRawYield->SetBinContent(8,0.20);
fRawYield->SetBinContent(9,0.15);
fRawYield->SetBinContent(10,0.15);
fRawYield->SetBinContent(11,0.15);
fRawYield->SetBinContent(12,0.15);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.);
for(Int_t iBin=7; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.05);
return;
}
void AliHFSystErr::InitDplustoKpipi2010PbPb4080() {
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
fRawYield->SetBinContent(1,1);
fRawYield->SetBinContent(2,1);
fRawYield->SetBinContent(3,1);
fRawYield->SetBinContent(4,0.15);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.15);
fRawYield->SetBinContent(8,0.15);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.15);
for(Int_t i=13;i<=20;i++) fRawYield->SetBinContent(i,1);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05);
fPIDEff->SetBinContent(3,0.13);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0);
for(Int_t iBin=4; iBin<=8; iBin++) fMCPtShape->SetBinContent(iBin,0.01);
for(Int_t iBin=9; iBin<=12; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
for(Int_t iBin=13; iBin<=16; iBin++) fMCPtShape->SetBinContent(iBin,0.03);
return;
}
void AliHFSystErr::InitDplustoKpipi2010ppLowEn() {
AliInfo(" Settings for D+ --> K pi pi p-p collisions at 2.76 TeV");
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
fRawYield->SetBinContent(1,1);
fRawYield->SetBinContent(2,1);
for(Int_t i=3;i<=6;i++) fRawYield->SetBinContent(i,0.10);
fRawYield->SetBinContent(7,0.15);
fRawYield->SetBinContent(8,0.15);
for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.055);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.15);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05);
fPIDEff->SetBinContent(3,0.10);
fPIDEff->SetBinContent(4,0.10);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(1,0.03);
fMCPtShape->SetBinContent(2,0.03);
return;
}
void AliHFSystErr::InitDplustoKpipi2013pPb0100(){
AliInfo(" Settings for D+ --> K pi pi, p-Pb collisions at 5.023 TeV");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.02);
fBR = new TH1F("fBR","fBR",20,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.021);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,0.);
fRawYield->SetBinContent(2,0.1);
fRawYield->SetBinContent(3,0.08);
fRawYield->SetBinContent(4,0.05);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.08);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.1);
fCutsEff->SetBinContent(3,0.1);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.05);
fCutsEff->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
fPIDEff->SetBinContent(1,1);
fPIDEff->SetBinContent(2,0.1);
for(Int_t i=3;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.0);
return;
}
void AliHFSystErr::InitDstartoD0pi2010pp() {
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.035);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
fTrackingEff->SetBinContent(1,1.0);
fTrackingEff->SetBinContent(2,0.13);
fTrackingEff->SetBinContent(3,0.12);
fTrackingEff->SetBinContent(3,0.12);
for(Int_t i=4;i<=24;i++) fTrackingEff->SetBinContent(i,0.12);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,1.0);
fRawYield->SetBinContent(2,0.10);
fRawYield->SetBinContent(3,0.04);
fRawYield->SetBinContent(4,0.03);
fRawYield->SetBinContent(5,0.03);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.2);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
fCutsEff->SetBinContent(2,0.22);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0);
fMCPtShape->SetBinContent(1,0.03);
fMCPtShape->SetBinContent(2,0.03);
return;
}
void AliHFSystErr::InitDstartoD0pi2010ppLowEn() {
AliInfo(" Settings for D*+ --> D0 pi p-p collisions at 2.76 TeV");
fNorm = new TH1F("fNorm","fNorm",20,0,20);
for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.031);
fBR = new TH1F("fBR","fBR",20,0,20);
for(Int_t i=1;i<=20;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",20,0,20);
for(Int_t i=1;i<=20;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",20,0,20);
fRawYield->SetBinContent(1,1);
fRawYield->SetBinContent(2,1);
fRawYield->SetBinContent(3,0.14);
fRawYield->SetBinContent(4,0.14);
fRawYield->SetBinContent(5,0.12);
fRawYield->SetBinContent(6,0.12);
fRawYield->SetBinContent(7,0.06);
fRawYield->SetBinContent(8,0.06);
fRawYield->SetBinContent(9,0.08);
fRawYield->SetBinContent(10,0.08);
fRawYield->SetBinContent(11,0.08);
fRawYield->SetBinContent(12,0.08);
for(Int_t i=9;i<=20;i++) fRawYield->SetBinContent(i,0.065);
fCutsEff = new TH1F("fCutsEff","fCutsEff",20,0,20);
for(Int_t i=1;i<=20;i++) fCutsEff->SetBinContent(i,0.10);
fCutsEff->SetBinContent(3,0.15);
fCutsEff->SetBinContent(4,0.15);
fCutsEff->SetBinContent(5,0.15);
fCutsEff->SetBinContent(6,0.15);
fCutsEff->SetBinContent(7,0.10);
fCutsEff->SetBinContent(8,0.10);
fCutsEff->SetBinContent(9,0.10);
fCutsEff->SetBinContent(10,0.10);
fCutsEff->SetBinContent(11,0.10);
fCutsEff->SetBinContent(12,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",20,0,20);
for(Int_t i=1;i<=20;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",20,0,20);
for(Int_t i=1;i<=20;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(1,0.03);
fMCPtShape->SetBinContent(2,0.03);
return;
}
void AliHFSystErr::InitDstartoD0pi2010PbPb020() {
AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-20 centrality - DUMMY");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=20;i++) fNorm->SetBinContent(i,0.05);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);
fRawYield->SetBinContent(3,0.2);
fRawYield->SetBinContent(4,0.2);
fRawYield->SetBinContent(5,0.2);
fRawYield->SetBinContent(6,0.2);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10);
fCutsEff->SetBinContent(4,0.15);
fCutsEff->SetBinContent(5,0.15);
fCutsEff->SetBinContent(6,0.15);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.045);
fMCPtShape->SetBinContent(4,0.025);
fMCPtShape->SetBinContent(5,0.025);
fMCPtShape->SetBinContent(6,0.025);
fMCPtShape->SetBinContent(7,0.04);
fMCPtShape->SetBinContent(8,0.04);
fMCPtShape->SetBinContent(9,0.03);
fMCPtShape->SetBinContent(10,0.03);
fMCPtShape->SetBinContent(11,0.03);
fMCPtShape->SetBinContent(12,0.03);
return;
}
void AliHFSystErr::InitDstartoD0pi2011PbPb07half() {
AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-7.5 centrality - DUMMY");
fNorm = new TH1F("fNorm","fNorm",36,0,36);
for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",36,0,36);
for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
for(Int_t i=1;i<=36;i++) fRawYield->SetBinContent(i,0.05);
fRawYield->SetBinContent(4,0.2);
fRawYield->SetBinContent(5,0.10);
fRawYield->SetBinContent(6,0.10);
fRawYield->SetBinContent(7,0.08);
for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.15);
fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
for(Int_t i=1;i<=36;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05);
fPIDEff->SetBinContent(4,0.09);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
for(Int_t i=1;i<=36;i++) fMCPtShape->SetBinContent(i,0.035);
fMCPtShape->SetBinContent(4,0.015);
fMCPtShape->SetBinContent(5,0.015);
fMCPtShape->SetBinContent(6,0.015);
fMCPtShape->SetBinContent(7,0.02);
fMCPtShape->SetBinContent(8,0.02);
fMCPtShape->SetBinContent(9,0.03);
fMCPtShape->SetBinContent(10,0.03);
fMCPtShape->SetBinContent(11,0.03);
fMCPtShape->SetBinContent(12,0.03);
return;
}
void AliHFSystErr::InitDstartoD0pi2010PbPb2040() {
AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 20-40 centrality - DUMMY");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.10);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.15);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.04);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.);
fMCPtShape->SetBinContent(1,0.03);
fMCPtShape->SetBinContent(2,0.03);
return;
}
void AliHFSystErr::InitDstartoD0pi2010PbPb4080() {
AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 40-80 centrality - DUMMY");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.2);
fRawYield->SetBinContent(1,1);
fRawYield->SetBinContent(2,0.15);
fRawYield->SetBinContent(3,0.15);
fRawYield->SetBinContent(4,0.15);
fRawYield->SetBinContent(5,0.15);
fRawYield->SetBinContent(6,0.10);
fRawYield->SetBinContent(7,0.10);
fRawYield->SetBinContent(8,0.10);
fRawYield->SetBinContent(9,0.11);
fRawYield->SetBinContent(10,0.11);
fRawYield->SetBinContent(11,0.11);
fRawYield->SetBinContent(12,0.11);
fRawYield->SetBinContent(13,0.08);
fRawYield->SetBinContent(14,0.08);
fRawYield->SetBinContent(15,0.08);
fRawYield->SetBinContent(16,0.08);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(2,0.05);
fMCPtShape->SetBinContent(3,0.05);
fMCPtShape->SetBinContent(4,0.05);
fMCPtShape->SetBinContent(5,0.04);
fMCPtShape->SetBinContent(6,0.02);
fMCPtShape->SetBinContent(7,0.04);
fMCPtShape->SetBinContent(8,0.04);
return;
}
void AliHFSystErr::InitDstartoD0pi2011PbPb010() {
AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 0-10% centrality ");
fNorm = new TH1F("fNorm","fNorm",36,0,36);
for(Int_t i=1;i<=36;i++) fNorm->SetBinContent(i,0.048);
fBR = new TH1F("fBR","fBR",36,0,36);
for(Int_t i=1;i<=36;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",36,0,36);;
for(Int_t i=1;i<=36;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",36,0,36);
for(Int_t i=1;i<=3;i++) fRawYield->SetBinContent(i, 0.);
fRawYield->SetBinContent(4,0.10);
fRawYield->SetBinContent(5,0.08);
fRawYield->SetBinContent(6,0.05);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=25;i<=36;i++) fRawYield->SetBinContent(i,0.10);
fCutsEff = new TH1F("fCutsEff","fCutsEff",36,0,36);
for(Int_t i=1;i<=3;i++) fCutsEff->SetBinContent(i, 0.);
fCutsEff->SetBinContent(4,0.15);
for(Int_t i=5;i<=36;i++) fCutsEff->SetBinContent(i,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",36,0,36);
for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i, 0.);
fMCPtShape->SetBinContent(4,0.06);
fMCPtShape->SetBinContent(5,0.03);
for(Int_t i=5;i<=36;i++) fMCPtShape->SetBinContent(i,0.01);
return;
}
void AliHFSystErr::InitDstartoD0pi2011PbPb3050() {
AliInfo(" Settings for D*+ --> D0pi Pb-Pb collisions at 2.76 TeV - 30-50 centrality - DUMMY");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);;
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.15);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
for(Int_t i=1;i<=24;i++) fRawYield->SetBinContent(i,0.1);
fRawYield->SetBinContent(1,0.12);
fRawYield->SetBinContent(2,0.12);
fRawYield->SetBinContent(3,0.12);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.06);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.08);
fRawYield->SetBinContent(8,0.05);
fRawYield->SetBinContent(9,0.05);
fRawYield->SetBinContent(10,0.05);
fRawYield->SetBinContent(11,0.05);
fRawYield->SetBinContent(12,0.05);
fRawYield->SetBinContent(13,0.05);
fRawYield->SetBinContent(14,0.05);
fRawYield->SetBinContent(15,0.05);
fRawYield->SetBinContent(16,0.05);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
for(Int_t i=1;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
fCutsEff->SetBinContent(1,0.16);
fCutsEff->SetBinContent(2,0.16);
fCutsEff->SetBinContent(3,0.16);
fCutsEff->SetBinContent(4,0.10);
fCutsEff->SetBinContent(5,0.10);
fCutsEff->SetBinContent(6,0.10);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.05);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
fMCPtShape->SetBinContent(2,0.11);
fMCPtShape->SetBinContent(3,0.11);
fMCPtShape->SetBinContent(4,0.03);
fMCPtShape->SetBinContent(5,0.02);
fMCPtShape->SetBinContent(6,0.01);
fMCPtShape->SetBinContent(7,0.01);
fMCPtShape->SetBinContent(8,0.01);
return;
}
void AliHFSystErr::InitDstartoD0pi2013pPb0100(){
AliInfo(" Settings for D* --> D0 pi, p-Pb collisions at 5.023 TeV");
fNorm = new TH1F("fNorm","fNorm",24,0,24);
for(Int_t i=1;i<=24;i++) fNorm->SetBinContent(i,0.07);
fBR = new TH1F("fBR","fBR",24,0,24);
for(Int_t i=1;i<=24;i++) fBR->SetBinContent(i,0.015);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",24,0,24);
for(Int_t i=1;i<=24;i++) fTrackingEff->SetBinContent(i,0.09);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,1.0);
fRawYield->SetBinContent(2,0.08);
fRawYield->SetBinContent(3,0.02);
fRawYield->SetBinContent(4,0.02);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.02);
fRawYield->SetBinContent(7,0.04);
fRawYield->SetBinContent(8,0.04);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.06);
fCutsEff = new TH1F("fCutsEff","fCutsEff",24,0,24);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
fPIDEff = new TH1F("fPIDEff","fPIDEff",24,0,24);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",24,0,24);
for(Int_t i=1;i<=3;i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=4;i<=6;i++) fMCPtShape->SetBinContent(i,0.);
for(Int_t i=7;i<=24;i++) fMCPtShape->SetBinContent(i,0.01);
return;
}
void AliHFSystErr::InitD0toKpi2010PbPb010CentScan(){
InitD0toKpi2010PbPb020();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitD0toKpi2010PbPb1020CentScan(){
InitD0toKpi2010PbPb020();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.17);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitD0toKpi2010PbPb2040CentScan(){
InitD0toKpi2010PbPb4080();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.14);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitD0toKpi2010PbPb4060CentScan(){
InitD0toKpi2010PbPb4080();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.06);
for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.11);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitD0toKpi2010PbPb6080CentScan(){
InitD0toKpi2010PbPb4080();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=3;i<=5;i++) fMCPtShape->SetBinContent(i,0.08);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitDplustoKpipi2010PbPb010CentScan(){
InitDplustoKpipi2010PbPb020();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.18);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.09);
}
void AliHFSystErr::InitDplustoKpipi2010PbPb1020CentScan(){
InitDplustoKpipi2010PbPb020();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.23);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitDplustoKpipi2010PbPb2040CentScan(){
InitDplustoKpipi2010PbPb020();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.095);
}
void AliHFSystErr::InitDplustoKpipi2010PbPb4060CentScan(){
InitDplustoKpipi2010PbPb4080();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitDplustoKpipi2010PbPb6080CentScan(){
InitDplustoKpipi2010PbPb4080();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.15);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.07);
}
void AliHFSystErr::InitDstartoD0pi2010PbPb010CentScan(){
InitDstartoD0pi2010PbPb020();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.16);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
}
void AliHFSystErr::InitDstartoD0pi2010PbPb1020CentScan(){
InitDstartoD0pi2010PbPb020();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.15);
}
void AliHFSystErr::InitDstartoD0pi2010PbPb2040CentScan(){
InitDstartoD0pi2010PbPb2040();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.08);
}
void AliHFSystErr::InitDstartoD0pi2010PbPb4060CentScan(){
InitDstartoD0pi2010PbPb4080();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
}
void AliHFSystErr::InitDstartoD0pi2010PbPb6080CentScan(){
InitDstartoD0pi2010PbPb4080();
for(Int_t i=7;i<=12;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=7;i<=12;i++) fMCPtShape->SetBinContent(i,0.045);
}
void AliHFSystErr::InitD0toKpi2011PbPb010CentScan(){
InitD0toKpi2011PbPb07half();
}
void AliHFSystErr::InitD0toKpi2011PbPb1020CentScan(){
InitD0toKpi2011PbPb07half();
}
void AliHFSystErr::InitD0toKpi2011PbPb2030CentScan(){
InitD0toKpi2011PbPb07half();
}
void AliHFSystErr::InitD0toKpi2011PbPb3040CentScan(){
InitD0toKpi2011PbPb07half();
}
void AliHFSystErr::InitD0toKpi2011PbPb4050CentScan(){
InitD0toKpi2011PbPb07half();
}
void AliHFSystErr::InitD0toKpi2010PbPb5080CentScan(){
InitD0toKpi2011PbPb07half();
}
void AliHFSystErr::InitDplustoKpipi2011PbPb010CentScan(){
InitDplustoKpipi2011PbPb07half();
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.15);
for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.15);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.10);
}
void AliHFSystErr::InitDplustoKpipi2011PbPb1020CentScan(){
InitDplustoKpipi2011PbPb07half();
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=9;i<=17;i++)fRawYield->SetBinContent(i,0.15);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.15);
}
void AliHFSystErr::InitDplustoKpipi2011PbPb2030CentScan(){
InitDplustoKpipi2011PbPb07half();
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
}
void AliHFSystErr::InitDplustoKpipi2011PbPb3040CentScan(){
InitDplustoKpipi2011PbPb07half();
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.20);
for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.20);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.20);
for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
}
void AliHFSystErr::InitDplustoKpipi2011PbPb4050CentScan(){
InitDplustoKpipi2011PbPb07half();
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=5;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.15);
for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
}
void AliHFSystErr::InitDplustoKpipi2010PbPb5080CentScan(){
InitDplustoKpipi2010PbPb4080();
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=17;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=1;i<=17;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.20);
for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10);
for(Int_t i=9;i<=17;i++) fCutsEff->SetBinContent(i,0.20);
}
void AliHFSystErr::InitDstartoD0pi2011PbPb010CentScan(){
InitDstartoD0pi2011PbPb07half();
}
void AliHFSystErr::InitDstartoD0pi2011PbPb1020CentScan(){
InitDstartoD0pi2011PbPb07half();
}
void AliHFSystErr::InitDstartoD0pi2011PbPb2030CentScan(){
InitDstartoD0pi2011PbPb07half();
}
void AliHFSystErr::InitDstartoD0pi2011PbPb3040CentScan(){
InitDstartoD0pi2011PbPb07half();
}
void AliHFSystErr::InitDstartoD0pi2011PbPb4050CentScan(){
InitDstartoD0pi2011PbPb07half();
}
void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
InitDstartoD0pi2011PbPb07half();
}
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.06);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.12);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
InitDplustoKpipi2013pPb0100();
AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.10);
};
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
InitDplustoKpipi2013pPb0100();
AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
InitDplustoKpipi2013pPb0100();
AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
InitDplustoKpipi2013pPb0100();
AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
InitDplustoKpipi2013pPb0100();
AliInfo("Settings for D+ --> K pi pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.07);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.8 < y < -0.4");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.4 < y < -0.1");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan -0.1 < y < 0.1");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.1 < y < 0.4");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.03);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D*+ --> D0 pi, pPb collisions at 5.023 TeV, rapidity scan 0.4 < y < 0.8");
for(Int_t i=2;i<=5;i++) fRawYield->SetBinContent(i,0.06);
for(Int_t i=6;i<=8;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0);
for(Int_t i=2;i<=5;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=6;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=2;i<=5;i++) fMCPtShape->SetBinContent(i,0.10);
for(Int_t i=6;i<=8;i++) fMCPtShape->SetBinContent(i,0.05);
for(Int_t i=9;i<=16;i++) fMCPtShape->SetBinContent(i,0.05);
};
void AliHFSystErr::InitD0toKpi2013pPb020V0A(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitD0toKpi2013pPb2040V0A(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitD0toKpi2013pPb4060V0A(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitD0toKpi2013pPb60100V0A(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitD0toKpi2013pPb020ZNA(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 0-20% mult ZNA");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitD0toKpi2013pPb2040ZNA(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 20-40% mult ZNA");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitD0toKpi2013pPb4060ZNA(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult ZNA");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitD0toKpi2013pPb60100ZNA(){
InitD0toKpi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 60-100% mult ZNA");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.085);
for(Int_t i=3;i<=4;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=5;i<=6;i++) fRawYield->SetBinContent(i,0.04);
for(Int_t i=7;i<=8;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.058);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.09);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
for(Int_t i=3;i<=24;i++) fCutsEff->SetBinContent(i,0.05);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitDstartoD0pi2013pPb020V0A(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.10);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
fCutsEff->SetBinContent(1,0.0);
fCutsEff->SetBinContent(2,0.10);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.0);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitDstartoD0pi2013pPb2040V0A(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.10);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.00);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.10);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitDstartoD0pi2013pPb4060V0A(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.0);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.0);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.00);
}
void AliHFSystErr::InitDstartoD0pi2013pPb60100V0A(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.0);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitDstartoD0pi2013pPb020ZNA(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 0-20% mult V0A");
fRawYield->SetBinContent(1,0.);
fRawYield->SetBinContent(2,0.10);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.10);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.);
}
void AliHFSystErr::InitDstartoD0pi2013pPb2040ZNA(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 20-40% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.0);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitDstartoD0pi2013pPb4060ZNA(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D0 --> K pi, pPb collisions at 5.023 TeV, 40-60% mult V0A");
fRawYield->SetBinContent(1,0.0);
fRawYield->SetBinContent(2,0.0);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitDstartoD0pi2013pPb60100ZNA(){
InitDstartoD0pi2013pPb0100();
AliInfo("Settings for D* --> D0 pi, pPb collisions at 5.023 TeV, 60-100% mult V0A");
fRawYield->SetBinContent(1,0.);
fRawYield->SetBinContent(2,0.);
fRawYield->SetBinContent(3,0.06);
fRawYield->SetBinContent(4,0.06);
fRawYield->SetBinContent(5,0.05);
fRawYield->SetBinContent(6,0.05);
fRawYield->SetBinContent(7,0.05);
fRawYield->SetBinContent(8,0.05);
for(Int_t i=9;i<=12;i++) fRawYield->SetBinContent(i,0.08);
for(Int_t i=13;i<=16;i++) fRawYield->SetBinContent(i,0.10);
for(Int_t i=17;i<=24;i++) fRawYield->SetBinContent(i,0.);
fCutsEff->SetBinContent(1,0.);
fCutsEff->SetBinContent(2,0.);
fCutsEff->SetBinContent(3,0.06);
fCutsEff->SetBinContent(4,0.06);
fCutsEff->SetBinContent(5,0.06);
fCutsEff->SetBinContent(6,0.06);
fCutsEff->SetBinContent(7,0.06);
fCutsEff->SetBinContent(8,0.06);
for(Int_t i=9;i<=12;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=13;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
for(Int_t i=17;i<=24;i++) fCutsEff->SetBinContent(i,0.);
for(Int_t i=1;i<=24;i++) fPIDEff->SetBinContent(i,0.0);
}
void AliHFSystErr::InitLctopKpi2010pp() {
fNorm = new TH1F("fNorm","fNorm",4,2,6);
for(Int_t i=1;i<=4;i++) fNorm->SetBinContent(i,0.035);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",4,2,6);
for(Int_t i=1;i<=4;i++) fTrackingEff->SetBinContent(i,0.12);
fRawYield = new TH1F("fRawYield","fRawYield",24,0,24);
fRawYield->SetBinContent(1,0.20);
fRawYield->SetBinContent(2,0.15);
fRawYield->SetBinContent(3,0.10);
fRawYield->SetBinContent(4,0.15);
fCutsEff = new TH1F("fCutsEff","fCutsEff",4,2,6);
fCutsEff->SetBinContent(1,0.50);
fCutsEff->SetBinContent(2,0.07);
fCutsEff->SetBinContent(3,0.30);
fCutsEff->SetBinContent(4,0.40);
fPIDEff = new TH1F("fPIDEff","fPIDEff",4,2,4);
fPIDEff->SetBinContent(1,0.13);
fPIDEff->SetBinContent(2,0.30);
fPIDEff->SetBinContent(3,0.5);
fPIDEff->SetBinContent(4,0.30);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0,12);
for(Int_t i=1; i<=2; i++) fMCPtShape->SetBinContent(i,1.);
for(Int_t i=3; i<=4; i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=5; i<=6; i++) fMCPtShape->SetBinContent(i,0.03);
for(Int_t i=7; i<=8; i++) fMCPtShape->SetBinContent(i,0.02);
for(Int_t i=9; i<=12; i++) fMCPtShape->SetBinContent(i,0.02);
}
void AliHFSystErr::InitLctopK0S2010pp() {
fNorm = new TH1F("fNorm","fNorm",12,0.,12.);
for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05);
fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0.,12.);
for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.1);
fRawYield = new TH1F("fRawYield","fRawYield",12,0.,12.);
for(Int_t i=1;i<=2;i++) fRawYield->SetBinContent(i,1.);
for(Int_t i=3;i<=12;i++) fRawYield->SetBinContent(i,0.1);
fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0.,12.);
for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.1);
fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0.,12.);
for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1);
fMCPtShape = new TH1F("fMCPtShape","fMCPtShape",12,0.,12.);
for(Int_t i=1; i<=2;i++) fMCPtShape->SetBinContent(i,1.);
for(Int_t i=3; i<=12;i++) fMCPtShape->SetBinContent(i,0.1);
}
Double_t AliHFSystErr::GetCutsEffErr(Double_t pt) const {
Int_t bin=fCutsEff->FindBin(pt);
return fCutsEff->GetBinContent(bin);
}
Double_t AliHFSystErr::GetMCPtShapeErr(Double_t pt) const {
Int_t bin=fMCPtShape->FindBin(pt);
return fMCPtShape->GetBinContent(bin);
}
Double_t AliHFSystErr::GetSeleEffErr(Double_t pt) const {
Double_t err=GetCutsEffErr(pt)*GetCutsEffErr(pt)+GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
return TMath::Sqrt(err);
}
Double_t AliHFSystErr::GetPIDEffErr(Double_t pt) const {
Int_t bin=fPIDEff->FindBin(pt);
return fPIDEff->GetBinContent(bin);
}
Double_t AliHFSystErr::GetTrackingEffErr(Double_t pt) const {
Int_t bin=fTrackingEff->FindBin(pt);
return fTrackingEff->GetBinContent(bin);
}
Double_t AliHFSystErr::GetRawYieldErr(Double_t pt) const {
Int_t bin=fRawYield->FindBin(pt);
return fRawYield->GetBinContent(bin);
}
Double_t AliHFSystErr::GetPartAntipartErr(Double_t pt) const {
Int_t bin=fPartAntipart->FindBin(pt);
return fPartAntipart->GetBinContent(bin);
}
Double_t AliHFSystErr::GetTotalSystErr(Double_t pt,Double_t feeddownErr) const {
Double_t err=0.;
if(fRawYield) err += GetRawYieldErr(pt)*GetRawYieldErr(pt);
if(fTrackingEff) err += GetTrackingEffErr(pt)*GetTrackingEffErr(pt);
if(fCutsEff) err += GetCutsEffErr(pt)*GetCutsEffErr(pt);
if(fPIDEff) err += GetPIDEffErr(pt)*GetPIDEffErr(pt);
if(fMCPtShape) err += GetMCPtShapeErr(pt)*GetMCPtShapeErr(pt);
if(fPartAntipart) err += GetPartAntipartErr(pt)*GetPartAntipartErr(pt);
err += feeddownErr*feeddownErr;
return TMath::Sqrt(err);
}
void AliHFSystErr::DrawErrors(TGraphAsymmErrors *grErrFeeddown) const {
gStyle->SetOptStat(0);
TCanvas *cSystErr = new TCanvas("cSystErr","Systematic Errors",300,80,640,500);
cSystErr->Range(0.20,-0.5,18.4,0.34);
cSystErr->SetRightMargin(0.318);
cSystErr->SetFillColor(0);
TH2F *hFrame = new TH2F("hFrame","Systematic errors; p_{t} (GeV/c); Relative Error",40,0,40,100,-1,+1);
hFrame->SetAxisRange(1.,35.9,"X");
hFrame->SetAxisRange(-0.5,0.5,"Y");
hFrame->Draw();
TLegend *leg = new TLegend(0.69,0.44,0.98,0.86,NULL,"brNDC");
leg->SetTextSize(0.03601695);
leg->SetFillStyle(0);
leg->SetBorderSize(0);
TH1F *hTotErr=new TH1F("hTotErr","",36,0,36);
Int_t nbins = fNorm->GetNbinsX();
TGraphAsymmErrors *gTotErr = new TGraphAsymmErrors(nbins);
for(Int_t i=1;i<=36;i++) {
Double_t pt = hTotErr->GetBinCenter(i);
Double_t ptwidth = hTotErr->GetBinWidth(i);
if(grErrFeeddown) {
Double_t x=0., y=0., errxl=0., errxh=0., erryl=0., erryh=0.;
Double_t toterryl=0., toterryh=0.;
for(Int_t j=0; j<grErrFeeddown->GetN(); j++) {
grErrFeeddown->GetPoint(j,x,y);
errxh = grErrFeeddown->GetErrorXhigh(j);
errxl = grErrFeeddown->GetErrorXlow(j);
if ( ( (x-errxl) <= pt) && ( (x+errxl) >= pt) ) {
erryh = grErrFeeddown->GetErrorYhigh(j);
erryl = grErrFeeddown->GetErrorYlow(j);
}
}
if (erryl>=1e-3) toterryl = GetTotalSystErr(pt,erryl);
else toterryl = GetTotalSystErr(pt);
if (erryh>=1e-3) toterryh = GetTotalSystErr(pt,erryh);
else toterryh = GetTotalSystErr(pt);
hTotErr->SetBinContent(i,toterryh);
gTotErr->SetPoint(i,pt,0.);
gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,toterryl,toterryh);
}
else {
hTotErr->SetBinContent(i,GetTotalSystErr(pt));
gTotErr->SetPoint(i,pt,0.);
gTotErr->SetPointError(i,ptwidth/2.,ptwidth/2.,GetTotalSystErr(pt),GetTotalSystErr(pt));
}
}
gTotErr->SetLineColor(kBlack);
gTotErr->SetFillColor(kRed);
gTotErr->SetFillStyle(3002);
gTotErr->Draw("2");
leg->AddEntry(gTotErr,"Total (excl. norm.)","f");
fNorm->SetFillColor(1);
fNorm->SetFillStyle(3002);
Double_t norm = fNorm->GetBinContent(1)*100;
leg->AddEntry(fNorm,Form("Normalization (%.1f%s)",norm,"%"),"");
if(grErrFeeddown) {
grErrFeeddown->SetFillColor(kTeal-8);
grErrFeeddown->SetFillStyle(3001);
grErrFeeddown->Draw("2");
leg->AddEntry(grErrFeeddown,"Feed-down from B","f");
}
if(fTrackingEff) {
fTrackingEff->SetFillColor(4);
fTrackingEff->SetFillStyle(3006);
fTrackingEff->Draw("same");
TH1F *hTrackingEffRefl = ReflectHisto(fTrackingEff);
hTrackingEffRefl->Draw("same");
leg->AddEntry(fTrackingEff,"Tracking efficiency","f");
}
if(fBR) {
fBR->SetFillColor(6);
fBR->SetFillStyle(3005);
fBR->Draw("same");
TH1F *hBRRefl = ReflectHisto(fBR);
hBRRefl->Draw("same");
leg->AddEntry(fBR,"Branching ratio","f");
}
if(fRawYield) {
Int_t ci;
ci = TColor::GetColor("#00cc00");
fRawYield->SetLineColor(ci);
fRawYield->SetLineWidth(3);
fRawYield->Draw("same");
TH1F *hRawYieldRefl = ReflectHisto(fRawYield);
hRawYieldRefl->Draw("same");
leg->AddEntry(fRawYield,"Yield extraction","l");
}
if(fCutsEff) {
fCutsEff->SetLineColor(4);
fCutsEff->SetLineWidth(3);
fCutsEff->Draw("same");
TH1F *hCutsEffRefl = ReflectHisto(fCutsEff);
hCutsEffRefl->Draw("same");
leg->AddEntry(fCutsEff,"Cuts efficiency","l");
}
if(fPIDEff) {
fPIDEff->SetLineColor(7);
fPIDEff->SetLineWidth(3);
fPIDEff->Draw("same");
TH1F *hPIDEffRefl = ReflectHisto(fPIDEff);
hPIDEffRefl->Draw("same");
leg->AddEntry(fPIDEff,"PID efficiency","l");
}
if(fMCPtShape) {
Int_t ci = TColor::GetColor("#9933ff");
fMCPtShape->SetLineColor(ci);
fMCPtShape->SetLineWidth(3);
fMCPtShape->Draw("same");
TH1F *hMCPtShapeRefl = ReflectHisto(fMCPtShape);
hMCPtShapeRefl->Draw("same");
leg->AddEntry(fMCPtShape,"MC p_{t} shape","l");
}
if(fPartAntipart) {
Int_t ci = TColor::GetColor("#ff6600");
fPartAntipart->SetLineColor(ci);
fPartAntipart->SetLineWidth(3);
fPartAntipart->Draw("same");
TH1F *hPartAntipartRefl = ReflectHisto(fPartAntipart);
hPartAntipartRefl->Draw("same");
leg->AddEntry(fPartAntipart,"D = #bar{D}","l");
}
leg->Draw();
cSystErr->SaveAs("RelativeSystematics.eps");
return;
}
TH1F* AliHFSystErr::ReflectHisto(TH1F *hin) const {
TH1F *hout=(TH1F*)hin->Clone("hout");
hout->Scale(-1.);
return hout;
}