+//--------------------------------------------------------------------------
+void AliHFSystErr::InitDstoKKpi2011PbPb07half() {
+ //
+ // D+s->Kpipi syst errors. Responsible: G.M. Innocenti
+ // 2011 PbPb sample, 0-7.5% CC
+ //
+
+ // Normalization
+ fNorm = new TH1F("fNorm","fNorm",12,0,12);
+ for(Int_t i=1;i<=12;i++) fNorm->SetBinContent(i,0.05); // TAA and pp norm
+
+ // Branching ratio
+ fBR = new TH1F("fBR","fBR",12,0,12);
+ for(Int_t i=1;i<=12;i++) fBR->SetBinContent(i,06); // 0.14/2.32 PDG2010
+
+ // Tracking efficiency
+ fTrackingEff = new TH1F("fTrackingEff","fTrackingEff",12,0,12);
+ for(Int_t i=1;i<=12;i++) fTrackingEff->SetBinContent(i,0.15); // Jacek, 5% per track
+
+ // Raw yield extraction
+ 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);
+
+
+
+ // Cuts efficiency (from cuts variation)
+ fCutsEff = new TH1F("fCutsEff","fCutsEff",12,0,12);
+ for(Int_t i=1;i<=12;i++) fCutsEff->SetBinContent(i,0.20); // 20%
+
+ // PID efficiency (from PID/noPID)
+ fPIDEff = new TH1F("fPIDEff","fPIDEff",12,0,12);
+ for(Int_t i=1;i<=12;i++) fPIDEff->SetBinContent(i,0.1); // 10%
+
+ // MC dN/dpt
+ 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);
+
+ // particle-antiparticle
+ /*
+ fPartAntipart = new TH1F("fPartAntipart","fPartAntipart",12,0,12);
+ fPartAntipart->SetBinContent(1,1);
+ fPartAntipart->SetBinContent(2,1);
+ fPartAntipart->SetBinContent(3,0.12);
+ for(Int_t i=4;i<=12;i++) fPartAntipart->SetBinContent(i,0.05); //5 to 12%
+ */
+
+ return;
+}
+