From ade8a56f86005a646a0d0e329e695892b187e6e5 Mon Sep 17 00:00:00 2001 From: rpreghen Date: Wed, 14 Mar 2012 13:14:26 +0000 Subject: [PATCH 1/1] spectra-utils update --- PWGLF/SPECTRA/PiKaPr/COMBINED/SpectraUtils.C | 43 ++++++++++++++++++++ 1 file changed, 43 insertions(+) diff --git a/PWGLF/SPECTRA/PiKaPr/COMBINED/SpectraUtils.C b/PWGLF/SPECTRA/PiKaPr/COMBINED/SpectraUtils.C index 95de1f93a4e..ef527616267 100644 --- a/PWGLF/SPECTRA/PiKaPr/COMBINED/SpectraUtils.C +++ b/PWGLF/SPECTRA/PiKaPr/COMBINED/SpectraUtils.C @@ -603,6 +603,25 @@ Convert_dNdy_1over2pipt_dNdeta(TH1 *hin, Double_t mass, Double_t eta = 0.8) return hout; } +TH1 * +Convert_dNdy_1over2pipt_dNdy(TH1 *hin) +{ + + TH1 *hout = hin->Clone("hout"); + hout->Reset(); + Double_t pt, mt, conv, val, vale; + for (Int_t ibin = 0; ibin < hin->GetNbinsX(); ibin++) { + pt = hin->GetBinCenter(ibin + 1); + val = hin->GetBinContent(ibin + 1); + vale = hin->GetBinError(ibin + 1); + val /= (2. * TMath::Pi() * pt); + vale /= (2. * TMath::Pi() * pt); + hout->SetBinContent(ibin + 1, val); + hout->SetBinError(ibin + 1, vale); + } + return hout; +} + TH1 * Convert_dNdy_dNdeta(TH1 *hin, Double_t mass, Double_t eta = 0.8) { @@ -623,6 +642,30 @@ Convert_dNdy_dNdeta(TH1 *hin, Double_t mass, Double_t eta = 0.8) return hout; } +TGraph * +Convert_dNdy_dNdeta(TGraph *hin, Double_t mass, Double_t eta = 0.8) +{ + + TGraph *hout = hin->Clone("hout"); + // hout->Reset(); + Double_t pt, mt, conv, val, valelo, valehi; + for (Int_t ibin = 0; ibin < hin->GetN(); ibin++) { + pt = hin->GetX()[ibin]; + conv = eta2y(pt, mass, eta) / eta; + val = hin->GetY()[ibin]; + valelo = hin->GetEYlow()[ibin]; + valehi = hin->GetEYhigh()[ibin]; + val *= conv; + valelo *= conv; + valehi *= conv; + hout->GetX()[ibin] = pt; + hout->GetY()[ibin] = val; + hout->GetEYlow()[ibin] = valelo; + hout->GetEYhigh()[ibin] = valehi; + } + return hout; +} + TH1 * SummedId_1over2pipt_dNdeta(const Char_t *filename, Int_t icent) { -- 2.39.3