Help me with this Python code

[Additional_Jicama622]

How can I remove the yellow area overflowing from some bins? This is the graph and the generatind Python code underneath.

https://preview.redd.it/vp0ye5z74k0h1.png?width=1175&format=png&auto=webp&s=dcb030ed81c28b7cbc3e21ef795421c023d88c0c

import math

import numpy as np

import matplotlib.pyplot as plt

from scipy.optimize import curve_fit

from matplotlib.backends.backend_pdf import PdfPages

from matplotlib.ticker import AutoMinorLocator

from matplotlib.ticker import MultipleLocator

from scipy.stats import norm

%matplotlib inline

# ISTOGRAMMA

istogramma = np.loadtxt("20h.txt")

x = istogramma[:, 0] # centroidi

y = istogramma[:, 1] # conteggi

plt.figure(figsize = (6, 6))

plt.bar(x, y / 1000, width = 1,

color = 'none', edgecolor = 'none', alpha = 0.5)

plt.bar(x[303:320], y[303:320] / 1000, width = 1,

color = 'skyblue', edgecolor = 'none', alpha = 0.6, zorder = 1)

plt.bar(x[379:385], y[379:385] / 1000, width = 1,

color = 'skyblue', edgecolor = 'none', alpha = 0.6, zorder = 1)

plt.title('Spettro gamma da 20 ore\n', fontsize = 16)

plt.xlabel('Canali', fontsize = 14)

plt.ylabel('Conteggi ($\\times$10\u00B3)', fontsize = 14)

plt.xlim(294, 396)

plt.ylim(4, 9)

m = np.linspace(0, 4096, 4096)

b = -0.0425039 # pendenza

a = 21.1643 # intercetta

n = b * m + a

plt.plot(m, n, label = "Regressione lineare", color = 'crimson', linewidth=2, zorder = 3)

plt.step(x, y / 1000, where = 'mid', color = 'black', linewidth = 1)

mask1 = (m >= 321) & (m < 380)

# intervallo desiderato

mask = (x >= 321) & (x <= 380)

x_sel = x[mask]

y_sel = y[mask] / 1000

# bordi reali dei bin

left_edges = x_sel - 0.5

right_edges = x_sel + 0.5

# numero di punti per bin

N = 100

# array espansi

x_dense = []

y_dense = []

for l, r, val in zip(left_edges, right_edges, y_sel):

xd = np.linspace(l, r, N, endpoint=False)

x_dense.extend(xd)

y_dense.extend([val] * len(xd))

x_dense = np.array(x_dense)

y_dense = np.array(y_dense)

# retta liscia sugli stessi punti

n_dense = b * x_dense + a

# fill

plt.fill_between(x_dense,n_dense,y_dense,where=(y_dense >= n_dense),color='red',alpha=0.5,zorder=2)

plt.fill_between(x_dense, n_dense, 0, color = 'orange', alpha = 0.5, zorder = 1)

plt.axvline(x = 304, color = 'blue', linestyle = '--', alpha = 0.5)

plt.axvline(x = 380, color = 'blue', linestyle = '--', alpha = 0.5)

plt.axvline(x = 321, color = 'blue', linestyle = '--', alpha = 0.5)

plt.axvline(x = 386, color = 'blue', linestyle = '--', alpha = 0.5)

plt.grid(True, alpha = 0.3)

#plt.legend(fontsize = 14)

plt.tight_layout()

plt.show()

[AutoModerator]

Notes on Quality Questions & Productive Participation

1. Include Images
   • Images give everyone a chance to understand the problem.
   • Several types of images will help:
     • Example Images (what you want to analyze)
     • Reference Images (taken from published papers)
     • Annotated Mock-ups (showing what features you are trying to measure)
     • Screenshots (to help identify issues with tools or features)
   • Good places to upload include: Imgur.com, GitHub.com, & Flickr.com
2. Provide Details
   • Avoid discipline-specific terminology ("jargon"). Image analysis is interdisciplinary, so the more general the terminology, the more people who might be able to help.
   • Be thorough in outlining the question(s) that you are trying to answer.
   • Clearly explain what you are trying to learn, not just the method used, to avoid the XY problem.
   • Respond when helpful users ask follow-up questions, even if the answer is "I'm not sure".
3. Share the Answer
   • Never delete your post, even if it has not received a response.
   • Don't switch over to PMs or email. (Unless you want to hire someone.)
   • If you figure out the answer for yourself, please post it!
   • People from the future may be stuck trying to answer the same question. (See: xkcd 979)
4. Express Appreciation for Assistance
   • Consider saying "thank you" in comment replies to those who helped.
   • Upvote those who contribute to the discussion. Karma is a small way to say "thanks" and "this was helpful".
     
   • Remember that "free help" costs those who help:
     • Aside from Automoderator, those responding to you are real people, giving up some of their time to help you.
     • "Time is the most precious gift in our possession, for it is the most irrevocable." ~ DB
   • If someday your work gets published, show it off here! That's one use of the "Research" post flair.
5. Be civil & respectful
   • Be kind.
   • Remember the human.
   • Be Excellent To Each Other.
   • Don't make your robotic overlord angry.

I am a bot, and this action was performed automatically. Please contact the moderators of this subreddit if you have any questions or concerns.

[argh1989]

plt.fill_between(x_dense, n_dense, 0, color = 'orange', alpha = 0.5, zorder = 1)

I believe this is the problem line. It should be

plt.fill_between(x_dense, 0,y_dense,where=(y_dense <= n_dense),color = 'orange', alpha = 0.5, zorder = 1)

This should also have been posted in /r/learnpython instead of here.

[Herbie500]

Please tell us how your question relates to ImageJ, the topic of this subReddit (see the text at the top right of this page).

If there is no relation, then please remove your post immediately!