fRunNumber(10),
fCollisionType(0),
fCentralityClass("0100"),
-fIsLowEnergy(false)
+fIsLowEnergy(false),
+fIsCentScan(false)
{
//
// Default Constructor
if (fIsLowEnergy) InitD0toKpi2010ppLowEn();
else InitD0toKpi2010pp();
} else if (fCollisionType==1) {
- if (fCentralityClass=="020") InitD0toKpi2010PbPb020();
+ 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 AliError("Not yet implemented");
}
if (fIsLowEnergy) InitDplustoKpipi2010ppLowEn();
else InitDplustoKpipi2010pp();
} else if (fCollisionType==1) {
- if (fCentralityClass=="020") InitDplustoKpipi2010PbPb020();
+ 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 AliError("Not yet implemented");
}
if(fIsLowEnergy) InitDstartoD0pi2010ppLowEn();
else InitDstartoD0pi2010pp();
}else if (fCollisionType==1) {
- if (fCentralityClass=="020") InitDstartoD0pi2010PbPb020();
- if (fCentralityClass=="2040") InitDstartoD0pi2010PbPb2040();
- if (fCentralityClass=="4080") InitDstartoD0pi2010PbPb4080();
- if (fCentralityClass!="4080" && fCentralityClass!="2040" && fCentralityClass!="020") AliError("Not yet implemented");
+ 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 AliError("Not yet implemented");
}
break;
case 4: // D+s->KKpi
fIsLowEnergy = flag;
AliInfo(" Settings for the low energy run");
}
+ void SetIsPbPb2010EnergyScan(Bool_t flag) {
+ fIsCentScan = flag;
+ AliInfo(" Settings for the PbPb 2010 energy scan");
+ }
// Function to initialize the variables/histograms
void Init(Int_t decay);
void InitDplustoKpipi2010ppLowEn();
void InitDstartoD0pi2010ppLowEn();
- void InitDstartoD0pi2010PbPb2040();
+ void InitDstartoD0pi2010PbPb2040();
TH1F* ReflectHisto(TH1F *hin) const;
TString fCentralityClass; // Centrality class
// MB:0100, 0-10:010, 0-20:020 ...40-80:4080...
Bool_t fIsLowEnergy; // flag for the low energy (2.76TeV) run
+ Bool_t fIsCentScan; // flag fot the PbPb centrality scan
- ClassDef(AliHFSystErr,3); // class for systematic errors of charm hadrons
+ ClassDef(AliHFSystErr,4); // class for systematic errors of charm hadrons
};
#endif