#!/usr/bin/env python
# coding: utf-8

# In[ ]:


import numpy as np
import matplotlib.pyplot as plt


# In[ ]:


from solcore import material, si


# In[ ]:


from solcore.solar_cell import SolarCell, Layer, Junction
from solcore.solar_cell_solver import solar_cell_solver, prepare_solar_cell
from solcore.state import State


# In[ ]:


GaInP = material("GaInP")(In=0.5)
GaAs = material("GaAs")()
Ge = material("Ge")()


# In[ ]:


optical_struct = SolarCell(
    [
        Layer(material=GaInP, width=si("5000nm")),
        Junction(
            [
                Layer(material=GaAs, width=si("200nm")),
                Layer(material=GaAs, width=si("5um")),
            ],
            kind="DA",
        ),
        Layer(material=Ge, width=si("50um")),
    ]
)


# In[ ]:


wl = np.linspace(250, 1700, 300) * 1e-9


# In[ ]:


options = State()
options.wavelength = wl
options.optics_method = "TMM"
options.no_back_reflection = False
options.BL_correction = True
options.recalculate_absorption = True
options.positions = [1e-8, 1e-9, 1e-8, 1e-7]
options.theta = 0


# In[ ]:


c_list = [
    ["c", "c", "c", "c"],
    ["c", "c", "c", "i"],
    ["c", "i", "i", "c"],
    ["i", "i", "i", "i"],
]


# In[ ]:


titles = [
    "All coherent",
    "Bottom Ge layer explicity incoherent",
    "Both layers of GaAs junction incoherent",
    "All layers incoherent",
]


# In[ ]:


for i1, cl in enumerate(c_list):
    plt.figure(i1)
    options.coherency_list = cl
    solar_cell_solver(optical_struct, "optics", options)
    plt.plot(wl * 1e9, optical_struct[0].layer_absorption)
    plt.plot(wl * 1e9, optical_struct[1].layer_absorption)
    plt.plot(wl * 1e9, optical_struct[2].layer_absorption)
    plt.plot(wl * 1e9, optical_struct.reflected, "--")
    plt.plot(wl * 1e9, optical_struct.transmitted, "--")
    plt.plot(
        wl * 1e9,
        optical_struct[0].layer_absorption
        + optical_struct[1].layer_absorption
        + optical_struct[2].layer_absorption
        + optical_struct.reflected
        + optical_struct.transmitted,
    )
    plt.legend(["GaInP", "GaAs", "Ge", "R", "T", "R+A+T"], loc="upper left")
    plt.title(titles[i1])
    plt.xlabel("Wavelength (nm)")
    plt.show()