diff --git a/CHANGELOG.md b/CHANGELOG.md
index 38038d909e..f60655c8eb 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -208,6 +208,7 @@ it in future.
   - Update Python code in shipDigiReco.py to use modern constructors
   - Bump ClassDef versions to 2 for schema evolution
   - Add TrackInfo to RNTuple I/O test suite
+* Fixed the weird cut on P in the script that produces P and P-Pt histograms for neutrinos
 
 ### Removed
 
diff --git a/macro/makeDecay.py b/macro/makeDecay.py
index e0d712f552..a27d24e96b 100644
--- a/macro/makeDecay.py
+++ b/macro/makeDecay.py
@@ -142,8 +142,8 @@
          if par.id()<0: idhnu+=1000
          pt2=par.px()**2+par.py()**2
          ptot=ROOT.TMath.Sqrt(pt2+par.pz()**2)
-         l10ptot=min(max(ROOT.TMath.Log10(ptot),-0.3),1.69999)
-         l10pt=min(max(ROOT.TMath.Log10(ROOT.TMath.Sqrt(pt2)),-2.),0.4999)
+         l10ptot= max(ROOT.TMath.Log10(ptot),-0.3)
+         l10pt= max(ROOT.TMath.Log10(ROOT.TMath.Sqrt(pt2)),-2.)
          h[str(idhnu)].Fill(ptot,wspill)
          h[str(idhnu+100)].Fill(l10ptot,l10pt,wspill)
          h[str(idhnu+200)].Fill(l10ptot,l10pt,wspill)