def lorentzian(x,dc,a,gamma,centre):
    """Single Lorentzian function. 
    
    Parameters:
    dc: float
        baseline
    gamma: float
           linewidth 
    centre: float
            peak centre
    x: array
       input parameter
    
    returns a Lorentzian function
    """
    lorentzian = dc + a*gamma**2 / ((x-centre)**2+gamma**2)
    return lorentzian