[PWGLF] using deltay instead of y in normalization histo

Author: alcaliva

PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

@@ -258,8 +258,8 @@ struct AntinucleiInJets {
       registryMC.add("antiproton_gen_ue", "antiproton_gen_ue", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
 
       // Normalization histogram
-      registryMC.add("antiproton_y_phi_jet", "antiproton_y_phi_jet", HistType::kTH2F, {{2000, -1.0, 1.0, "#it{y}"}, {2000, 0.0, TwoPI, "#Delta#phi"}});
-      registryMC.add("antiproton_y_phi_ue", "antiproton_y_phi_ue", HistType::kTH2F, {{2000, -1.0, 1.0, "#it{y}"}, {2000, 0.0, TwoPI, "#Delta#phi"}});
+      registryMC.add("antiproton_deltay_deltaphi_jet", "antiproton_deltay_deltaphi_jet", HistType::kTH2F, {{2000, -1.0, 1.0, "#Delta#it{y}"}, {2000, 0.0, 2.0, "#Delta#phi"}});
+      registryMC.add("antiproton_deltay_deltaphi_ue", "antiproton_deltay_deltaphi_ue", HistType::kTH2F, {{2000, -1.0, 1.0, "#Delta#it{y}"}, {2000, 0.0, 2.0, "#Delta#phi"}});
     }
 
     // Reconstructed antiproton spectra in jets and UE (MC-matched) with TPC/TOF PID
@@ -1464,12 +1464,7 @@ struct AntinucleiInJets {
             continue;
 
           // Fill normalization histogram
-          double px = particle.px();
-          double py = particle.py();
-          double pz = particle.pz();
-          double energy = std::sqrt(MassProton * MassProton + px * px + py * py + pz * pz);
-          double y = 0.5 * std::log((energy + pz) / (energy - pz));
-          registryMC.fill(HIST("antiproton_y_phi_jet"), y, getDeltaPhi(particle.phi(), jet.phi()));
+          registryMC.fill(HIST("antiproton_deltay_deltaphi_jet"), particle.eta() - jet.eta(), getDeltaPhi(particle.phi(), jet.phi()));
 
           // Fill histogram for generated antiprotons
           registryMC.fill(HIST("antiproton_gen_jet"), particle.pt());
@@ -1504,17 +1499,11 @@ struct AntinucleiInJets {
             continue;
 
           // Fill normalization histogram
-          double px = protonVec.Px();
-          double py = protonVec.Py();
-          double pz = protonVec.Pz();
-          double energy = std::sqrt(MassProton * MassProton + px * px + py * py + pz * pz);
-          double y = 0.5 * std::log((energy + pz) / (energy - pz));
-
           if (deltaRUe1 < maxConeRadius) {
-            registryMC.fill(HIST("antiproton_y_phi_ue"), y, getDeltaPhi(protonVec.Phi(), ueAxis1.Phi()));
+            registryMC.fill(HIST("antiproton_deltay_deltaphi_ue"), protonVec.Eta() - ueAxis1.Eta(), getDeltaPhi(protonVec.Phi(), ueAxis1.Phi()));
           }
           if (deltaRUe2 < maxConeRadius) {
-            registryMC.fill(HIST("antiproton_y_phi_ue"), y, getDeltaPhi(protonVec.Phi(), ueAxis2.Phi()));
+            registryMC.fill(HIST("antiproton_deltay_deltaphi_ue"), protonVec.Eta() - ueAxis2.Eta(), getDeltaPhi(protonVec.Phi(), ueAxis2.Phi()));
           }
 
           // Fill histogram for antiprotons in the UE