+//_________________________________________________________________________________________________________
+void AliHFPtSpectrum::CalculateCorrectedSpectrumNoFeedDown(){
+ //
+ // Compute the corrected spectrum with no feed-down correction
+ //
+
+
+ // declare the output histograms
+ fhYieldCorr = new TH1D("fhYieldCorr","corrected yield (no feed-down corr)",fnPtBins,fPtBinLimits);
+ fhYieldCorrMax = new TH1D("fhYieldCorrMax","max corrected yield (no feed-down corr)",fnPtBins,fPtBinLimits);
+ fhYieldCorrMin = new TH1D("fhYieldCorrMin","min corrected yield (no feed-down corr)",fnPtBins,fPtBinLimits);
+ fgYieldCorr = new TGraphAsymmErrors(fnPtBins+1);
+
+ //
+ // Do the calculation
+ //
+ for (Int_t ibin=1; ibin<=fnPtBins; ibin++) {
+ Double_t value = fhRECpt->GetBinContent(ibin) /fhRECpt->GetBinWidth(ibin);
+ Double_t errvalue = fhRECpt->GetBinError(ibin) /fhRECpt->GetBinWidth(ibin);
+ fhYieldCorr->SetBinContent(ibin,value);
+ fhYieldCorr->SetBinError(ibin,errvalue);
+ fhYieldCorrMax->SetBinContent(ibin,value);
+ fhYieldCorrMax->SetBinError(ibin,errvalue);
+ fhYieldCorrMin->SetBinContent(ibin,value);
+ fhYieldCorrMin->SetBinError(ibin,errvalue);
+ }
+
+ fAsymUncertainties=kFALSE;
+}
+