Adding methods for pPb rapidity measurements

diff --git a/PWGHF/vertexingHF/AliHFSystErr.cxx b/PWGHF/vertexingHF/AliHFSystErr.cxx
index 5b617376daad5524deb5f1e4757f32b0ed5b25aa..a6a956e2fc14184d38daa64e342b32dde3c4ed72 100644 (file)
--- a/PWGHF/vertexingHF/AliHFSystErr.cxx
+++ b/PWGHF/vertexingHF/AliHFSystErr.cxx
@@ -60,8 +60,10 @@ fPartAntipart(0),
 fRunNumber(10),
 fCollisionType(0),
 fCentralityClass("0100"),
+fRapidityRange("0101"),
 fIsLowEnergy(false),
-fIsCentScan(false)
+fIsCentScan(false),
+fIsRapidityScan(false)
 {
   //
   // Default Constructor
@@ -128,6 +130,13 @@ void AliHFSystErr::Init(Int_t decay){
     } 
     else if (fCollisionType==2) { 
       if (fCentralityClass=="0100") InitD0toKpi2013pPb0100();
+      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;
@@ -161,6 +170,13 @@ void AliHFSystErr::Init(Int_t decay){
     } 
     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;
@@ -194,6 +210,13 @@ void AliHFSystErr::Init(Int_t decay){
     }
     else if (fCollisionType==2) { 
       if (fCentralityClass=="0100") InitDstartoD0pi2013pPb0100();
+      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;
@@ -625,6 +648,7 @@ void AliHFSystErr::InitD0toKpi2013pPb0100(){
 
 }
 
+
 //--------------------------------------------------------------------------
 void AliHFSystErr::InitDplustoKpipi2010pp() {
   // 
@@ -1865,6 +1889,249 @@ void AliHFSystErr::InitDstartoD0pi2010PbPb5080CentScan(){
   
 }
 
+//_________________________________________________________________________
+void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0804(){
+  InitD0toKpi2013pPb0100();
+  
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.04);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.06);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
+  
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+ 
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.13);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.07);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
+  
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0401(){
+  InitD0toKpi2013pPb0100();
+  
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.04);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.08);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
+  
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+   for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.20);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0101(){
+InitD0toKpi2013pPb0100();
+
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.04);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.12);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
+
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+   for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.20);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0104(){
+  InitD0toKpi2013pPb0100();
+
+   for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.04);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.03);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.08);
+
+   fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+   for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+   for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.20);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.05);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.08);
+
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitD0toKpi2013pPb0100RapScan0408(){
+  InitD0toKpi2013pPb0100();
+
+   for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.04);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.03);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.04);
+
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+  
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.13);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.07);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
+  
+};
+
+//_________________________________________________________________________
+void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0804(){
+  InitDplustoKpipi2013pPb0100();
+
+   for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.07);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.10);
+
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+  
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.22);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.16);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0401(){
+  InitDplustoKpipi2013pPb0100();
+  
+   for(Int_t i=2;i<=4;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<=16;i++) fRawYield->SetBinContent(i,0.10);
+
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+   for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.24);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.06);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.03);
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0101(){
+  InitDplustoKpipi2013pPb0100();
+
+   for(Int_t i=2;i<=4;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<=16;i++) fRawYield->SetBinContent(i,0.05);
+
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+   for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.24);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.06);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.03);
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0104(){
+  InitDplustoKpipi2013pPb0100();
+
+   for(Int_t i=2;i<=4;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<=16;i++) fRawYield->SetBinContent(i,0.05);
+
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+   for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.24);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.06);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.03);
+  
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDplustoKpipi2013pPb0100RapScan0408(){
+  InitDplustoKpipi2013pPb0100();
+
+   for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.07);
+   for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
+
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+  
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.22);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.16);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
+
+};
+
+
+//_________________________________________________________________________
+void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0804(){
+  InitDstartoD0pi2013pPb0100();
+
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.06);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.03);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
+  
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.11);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.06);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0401(){
+  InitDstartoD0pi2013pPb0100();
+
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.04);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
+  
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.15);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.07);
+
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0101(){
+InitDstartoD0pi2013pPb0100();
+
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.05);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.03);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
+
+ fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.15);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.07);
+  
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0104(){
+  InitDstartoD0pi2013pPb0100();
+  
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.04);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.03);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.05);
+  
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+  
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.15);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.10);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.07);
+};
+//_________________________________________________________________________
+void AliHFSystErr::InitDstartoD0pi2013pPb0100RapScan0408(){
+  InitDstartoD0pi2013pPb0100();
+  
+  for(Int_t i=2;i<=4;i++) fRawYield->SetBinContent(i,0.06);
+  for(Int_t i=5;i<=8;i++) fRawYield->SetBinContent(i,0.07);
+  for(Int_t i=9;i<=16;i++) fRawYield->SetBinContent(i,0.07);
+  
+  fPIDEff = new TH1F("fPIDEff","fPIDEff",36,0,36);
+  for(Int_t i=1;i<=36;i++) fPIDEff->SetBinContent(i,0.0); // 0%
+  
+  for(Int_t i=2;i<=4;i++) fCutsEff->SetBinContent(i,0.11);
+  for(Int_t i=5;i<=8;i++) fCutsEff->SetBinContent(i,0.06);
+  for(Int_t i=9;i<=16;i++) fCutsEff->SetBinContent(i,0.06);
+  
+};
 
 
 //--------------------------------------------------------------------------

diff --git a/PWGHF/vertexingHF/AliHFSystErr.h b/PWGHF/vertexingHF/AliHFSystErr.h
index 1ac0847542ecf2bd69d920a8cbd0793049d5044a..f6b0e7c02923f033a6aa16098f90a4f7c7893223 100644 (file)
--- a/PWGHF/vertexingHF/AliHFSystErr.h
+++ b/PWGHF/vertexingHF/AliHFSystErr.h
@@ -70,6 +70,12 @@ class AliHFSystErr : public TNamed
     AliInfo(" Settings for the PbPb 2010 energy scan");
   }
 
+  void SetIspPb2011RapidityScan(Bool_t flag){
+    fIsRapidityScan = flag; 
+    AliInfo("Settings for the pPb vs y measurement");
+  }
+
+
   // Function to initialize the variables/histograms
   void Init(Int_t decay);
 
@@ -115,6 +121,24 @@ class AliHFSystErr : public TNamed
   void InitDstartoD0pi2011PbPb4050CentScan();
   void InitDstartoD0pi2010PbPb5080CentScan();
 
+  void InitD0toKpi2013pPb0100RapScan0804();
+  void InitD0toKpi2013pPb0100RapScan0401();
+  void InitD0toKpi2013pPb0100RapScan0101();
+  void InitD0toKpi2013pPb0100RapScan0104();
+  void InitD0toKpi2013pPb0100RapScan0408();
+
+  void InitDplustoKpipi2013pPb0100RapScan0804();
+  void InitDplustoKpipi2013pPb0100RapScan0401();
+  void InitDplustoKpipi2013pPb0100RapScan0101();
+  void InitDplustoKpipi2013pPb0100RapScan0104();
+  void InitDplustoKpipi2013pPb0100RapScan0408();
+
+  void InitDstartoD0pi2013pPb0100RapScan0804();
+  void InitDstartoD0pi2013pPb0100RapScan0401();
+  void InitDstartoD0pi2013pPb0100RapScan0101();
+  void InitDstartoD0pi2013pPb0100RapScan0104();
+  void InitDstartoD0pi2013pPb0100RapScan0408();
+
  private:
 
   AliHFSystErr(const AliHFSystErr& source);
@@ -165,10 +189,13 @@ class AliHFSystErr : public TNamed
   Int_t fCollisionType;    // Collision type: pp=0, PbPb=1
   TString fCentralityClass;  // Centrality class
                            // MB:0100, 0-10:010, 0-20:020 ...40-80:4080...
+  TString fRapidityRange;  // Rapidity range fot y measurements 
+
   Bool_t fIsLowEnergy;     // flag for the low energy (2.76TeV) run
   Bool_t fIsCentScan;      // flag fot the PbPb centrality scan
-
-  ClassDef(AliHFSystErr,4);  // class for systematic errors of charm hadrons
+  Bool_t fIsRapidityScan;  // flag for the pPb vs y measurement  
+ 
+  ClassDef(AliHFSystErr,5);  // class for systematic errors of charm hadrons
 };
 
 #endif