//const Int_t binBeginTPC[maxBins] = {-1, -1, -1, -1, -1, -1, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, -1, -1, -1, -1, -1};
//const Int_t binBeginTPCITS[maxBins] = {-1, -1, -1, -1, -1, -1, -1, 14, 10, 8, 7, 6, 5, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, -1, -1, -1, -1, -1, -1, -1};
+ const Int_t binBeginSPD[maxBins1] = {19, 18, 17, 15, 14, 12, 10, 9, 8, 7, 6, 6, 5, 5, 5, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4};
+ const Int_t binBeginTPC[maxBins1] = {-1, -1, -1, -1, -1, -1, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, -1, -1, -1, -1, -1};
+ const Int_t binBeginTPCITS[maxBins] = {-1, -1, -1, -1, -1, -1, 4, 4, 4, 4, 4, 4, 4, 4, -1, -1, -1, -1, -1, -1};
if (fAnalysisMode & AliPWG0Helper::kSPD)
{
- const Int_t binBeginSPD[maxBins1] = {19, 18, 17, 15, 14, 12, 10, 9, 8, 7, 6, 6, 5, 5, 5, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4};
binBegin = binBeginSPD;
}
else if (fAnalysisMode & AliPWG0Helper::kTPC)
{
- const Int_t binBeginTPC[maxBins1] = {-1, -1, -1, -1, -1, -1, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, -1, -1, -1, -1, -1};
binBegin = binBeginTPC;
}
else if (fAnalysisMode & AliPWG0Helper::kTPCITS)
{
- const Int_t binBeginTPCITS[maxBins] = {-1, -1, -1, -1, -1, -1, 4, 4, 4, 4, 4, 4, 4, 4, -1, -1, -1, -1, -1, -1};
binBegin = binBeginTPCITS;
}
fPtDist = dynamic_cast<TH1F*> (gDirectory->Get(fPtDist->GetName()));
- if (dynamic_cast<TNamed*> (gDirectory->Get("fTag")))
- fTag = (dynamic_cast<TNamed*> (gDirectory->Get("fTag")))->GetTitle();
+ if (static_cast<TNamed*> (gDirectory->Get("fTag")))
+ fTag = (static_cast<TNamed*> (gDirectory->Get("fTag")))->GetTitle();
- if (dynamic_cast<TParameter<Int_t>*> (gDirectory->Get("fAnalysisMode")))
- fAnalysisMode = (AliPWG0Helper::AnalysisMode) (dynamic_cast<TParameter<Int_t>*> (gDirectory->Get("fAnalysisMode")))->GetVal();
+ if (static_cast<TParameter<Int_t>*> (gDirectory->Get("fAnalysisMode")))
+ fAnalysisMode = (AliPWG0Helper::AnalysisMode) (static_cast<TParameter<Int_t>*> (gDirectory->Get("fAnalysisMode")))->GetVal();
gDirectory->cd("../");
}
fChipsFired = dynamic_cast<TH2F*> (file->Get("fChipsFired"));
fClusterZL1 = dynamic_cast<TH1F*> (file->Get("fClusterZL1"));
fClusterZL2 = dynamic_cast<TH1F*> (file->Get("fClusterZL2"));
+
+ if (!fPrimaryL1 || !fPrimaryL2 || !fChipsFired || !fClusterZL1 || !fClusterZL2)
+ return;
#define MULT 1001, -0.5, 1000.5
#define BINNING_LAYER1 401, -0.5, 400.5
Double_t windowsPerSecond = 1.0 / windowSize;
TH1* triggerEffHist = (TH1*) GetTriggerEfficiency(fMvsL, cut)->Clone("triggerEff");
- Float_t* triggerEff = new Float_t[max];
+ Float_t* triggerEff = new Float_t[max*2];
for (Int_t mult = 0; mult < max; mult++)
triggerEff[mult] = triggerEffHist->GetBinContent(mult+1);
Int_t cut = cuts[currentCut];
TH1* triggerEffHist = (TH1*) GetTriggerEfficiency(fMvsL, cut)->Clone("triggerEff");
- Float_t* triggerEff = new Float_t[max];
+ Float_t* triggerEff = new Float_t[max*2];
for (Int_t mult = 0; mult < max; mult++)
triggerEff[mult] = triggerEffHist->GetBinContent(mult+1);
Printf("Cut at %d", cut);
TH1* triggerEffHist = (TH1*) GetTriggerEfficiency(fMvsL, cut)->Clone("triggerEff");
- Float_t* triggerEff = new Float_t[max];
+ Float_t* triggerEff = new Float_t[max*2];
for (Int_t mult = 0; mult < max; mult++)
triggerEff[mult] = triggerEffHist->GetBinContent(mult+1);
return 0;
TH1* xSection;
- xSection = dynamic_cast<TH1*> (gFile->Get("xSection2Ex"));
+ xSection = static_cast<TH1*> (gFile->Get("xSection2Ex"));
TGraph* result = new TGraph;
for (Int_t i = 0; i < kCorrHists; ++i)
{
canvas3->cd(i+1);
- TH3* hist = dynamic_cast<TH3*> (fCorrelation[i]->Clone());
+ TH3* hist = static_cast<TH3*> (fCorrelation[i]->Clone());
for (Int_t y=1; y<=hist->GetYaxis()->GetNbins(); ++y)
{
for (Int_t z=1; z<=hist->GetZaxis()->GetNbins(); ++z)
// plot (MC - Unfolded) / error (MC)
canvas1->cd(3);
- TH1* diffMCUnfolded2 = dynamic_cast<TH1*> (proj->Clone("diffMCUnfolded2"));
+ TH1* diffMCUnfolded2 = static_cast<TH1*> (proj->Clone("diffMCUnfolded2"));
diffMCUnfolded2->Add(esdCorrected, -1);
Int_t ndfQual[kQualityRegions];
Float_t regPar = 0;
- Int_t correlationID = inputRange + ((fullPhaseSpace == kFALSE) ? 0 : 4);
- Int_t mcTarget = ((fullPhaseSpace == kFALSE) ? inputRange : 4);
+ Int_t correlationID = inputRange; // + ((fullPhaseSpace == kFALSE) ? 0 : 4);
+ Int_t mcTarget = inputRange; //((fullPhaseSpace == kFALSE) ? inputRange : 4);
SetupCurrentHists(inputRange, fullPhaseSpace, eventType);
//normalize ESD
// correct spectrum using a simple Gaussian approach, that is model-dependent
//
- Int_t correlationID = inputRange + ((fullPhaseSpace == kFALSE) ? 0 : 4);
- Int_t mcTarget = ((fullPhaseSpace == kFALSE) ? inputRange : 4);
+ Int_t correlationID = inputRange; // + ((fullPhaseSpace == kFALSE) ? 0 : 4);
+ Int_t mcTarget = inputRange; //((fullPhaseSpace == kFALSE) ? inputRange : 4);
SetupCurrentHists(inputRange, fullPhaseSpace, kTrVtx);
//normalize ESD
}
}
- TH2F* target = dynamic_cast<TH2F*> (fMultiplicityESD[0]->Clone(Form("%s_measured", inputMC->GetName())));
+ TH2F* target = static_cast<TH2F*> (fMultiplicityESD[0]->Clone(Form("%s_measured", inputMC->GetName())));
target->Reset();
for (Int_t meas=1; meas<=corr->GetNbinsY(); ++meas)