spectra-utils update

diff --git a/PWGLF/SPECTRA/PiKaPr/COMBINED/SpectraUtils.C b/PWGLF/SPECTRA/PiKaPr/COMBINED/SpectraUtils.C
index 95de1f93a4edfdb346df17d6a336b688fc3baf43..ef527616267504f0408900cdc5cf770739cc74e3 100644 (file)
--- 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)
 {