from matplotlib.patches import PathPatch
from matplotlib.path import Path

vertices_mu1 = []
codes_mu1= []
codes_mu1 += [Path.MOVETO] + [Path.LINETO]*3 + [Path.CLOSEPOLY]
vertices_mu1 += [(15, 10), (30,10), (30,0.1),(15,0.1),(0, 0)]  
path_mu1 = Path(vertices_mu1, codes_mu1)
pathpatch_mu1 = PathPatch(path_mu1, facecolor='None', edgecolor='black',linewidth=3)

vertices_mu2 = []
codes_mu2 = []
codes_mu2 += [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_mu2 += [(35,6),  (45, 6),(50, 10), (50,20),(37,20),(37,10),(0, 0)] 
path_mu2 = Path(vertices_mu2, codes_mu2)
pathpatch_mu2 = PathPatch(path_mu2, facecolor='None', edgecolor='red',linewidth=3)

vertices_mu3 = []
codes_mu3 = []
codes_mu3 += [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_mu3 += [(37,0.2),  (47, 0.2), (50, 0.1), (50,0.08),(37,0.08),(37,0.1),(0, 0)] 
path_mu3 = Path(vertices_mu3, codes_mu3)
pathpatch_mu3 = PathPatch(path_mu3, facecolor='None', edgecolor='red',linewidth=3)

vertices_pr1 = []
codes_pr1 = []
codes_pr1+= [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_pr1 += [(65,3),  (80, 9),(80, 20), (110,20),(110,9),(95,3),(0, 0)] 
path_pr1 = Path(vertices_pr1, codes_pr1)
pathpatch_pr1 = PathPatch(path_pr1, facecolor='None', edgecolor='red',linewidth=3)

vertices_pr3 = []
codes_pr3 = []
codes_pr3 += [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_pr3 += [(65,0.4),  (80, 0.2),(80, 0.1), (110,0.1),(110,0.2),(95,0.4),(0, 0)] 
path_pr3 = Path(vertices_pr3, codes_pr3)
pathpatch_pr3 = PathPatch(path_pr3, facecolor='None', edgecolor='red',linewidth=3)
# alphas >700 MeV
vertices_al1 = []
codes_a1 = []
codes_a1 += [Path.MOVETO] + [Path.LINETO]*3 + [Path.CLOSEPOLY]
vertices_al1 += [(70, 1.6), (180, 1.6), (180, 0.5), (70,0.5),(0, 0)]  
path_a1 = Path(vertices_al1, codes_a1)
pathpatch_a1 = PathPatch(path_a1, facecolor='None', edgecolor='black',linewidth=3)

# alphas 100 MeV -> ~550 MeV forwards
vertices_a2 = []
codes_a2 = []
codes_a2 += [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_a2 += [(180, 1.6), (355, 1.6), (460, 3.5), (460,14),(320, 14), (320,3.5),(0,0)]  
path_a2 = Path(vertices_a2, codes_a2)
pathpatch_a2 = PathPatch(path_a2, facecolor='None', edgecolor='black',linewidth=3)
# alphas 100 MeV -> ~550 MeV backwards
vertices_a3 = []
codes_a3 = []
codes_a3 += [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_a3 += [ (355, 0.5), (470, 0.2), (470,0.07),(350, 0.07), (350,0.2), (180,0.5),(0,0)]  
path_a3 = Path(vertices_a3, codes_a3)
pathpatch_a3 = PathPatch(path_a3, facecolor='None', edgecolor='black',linewidth=3)
#alphas ~550 MeV -700 MeV fowards and backwards
vertices_a4 = []
codes_a4 = []
codes_a4 += [Path.MOVETO] + [Path.LINETO]*3 + [Path.CLOSEPOLY]
vertices_a4 += [ (180, 1.6), (355, 1.6), (355,0.5),(180, 0.5),(0,0)]  
path_a4 = Path(vertices_a4, codes_a4)
pathpatch_a4 = PathPatch(path_a4, facecolor='None', edgecolor='black',linewidth=3)

# Boxes CNO
# 4000-openend
vertices_cno1 = []
codes_cno1 = []
codes_cno1 += [Path.MOVETO] + [Path.LINETO]*3 + [Path.CLOSEPOLY]
vertices_cno1 += [(700, 1.5), (700, 0.7), (1300,0.7),(1300, 1.5),(0,0)]  
path_cno1 = Path(vertices_cno1, codes_cno1)
pathpatch_cno1 = PathPatch(path_cno1, facecolor='None', edgecolor='black',linewidth=3)
# 3000-4000
vertices_cno4 = []
codes_cno4 = []
codes_cno4 += [Path.MOVETO] + [Path.LINETO]*3 + [Path.CLOSEPOLY]
vertices_cno4 += [(1300, 1.5), (2700, 1.5), (2700,0.7),(1300, 0.7),(0,0)]  
path_cno4 = Path(vertices_cno4, codes_cno4)
pathpatch_cno4 = PathPatch(path_cno4, facecolor='None', edgecolor='black',linewidth=3)

# 2400-3500
vertices_cno3 = []
codes_cno3 = []
codes_cno3 += [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_cno3 += [(1000, 0.7), (1700, 0.3), (1700,0.05),(4500, 0.05),(4500, 0.3),(2800, 0.7),(0,0)]  
path_cno3 = Path(vertices_cno3, codes_cno3)
pathpatch_cno3 = PathPatch(path_cno3, facecolor='None', edgecolor='black',linewidth=3)
#3000-4000
vertices_cno2= []
codes_cno2= []
codes_cno2 += [Path.MOVETO] + [Path.LINETO]*5 + [Path.CLOSEPOLY]
vertices_cno2 += [(1000, 1.5), (1700, 4.0), (1700,16),(4500, 16),(4500, 4.0),(2800,1.5),(0,0)]  
path_cno2 = Path(vertices_cno2, codes_cno2)
pathpatch_cno2 = PathPatch(path_cno2, facecolor='None', edgecolor='black',linewidth=3)