cov=cov,
        tau=tau,
    )

MASKS_SOME_PLACES = np.datetime64('2021-12-23')
MASKS_ALL_PLACES = np.datetime64('2022-01-02')
END_MASKS = np.datetime64('2022-03-04')

fig1 = plt.figure(figsize=(10, 6))
ax1 = plt.axes()

ax1.fill_betweenx(
    [-10, 10],
    [MASKS_SOME_PLACES, MASKS_SOME_PLACES],
    [MASKS_ALL_PLACES, MASKS_ALL_PLACES],
    color=whiten("green", 0.25),
    linewidth=0,
    label="Masks in high-risk environments only",
)

ax1.fill_betweenx(
    [-10, 10],
    [MASKS_ALL_PLACES, MASKS_ALL_PLACES],
    [END_MASKS, END_MASKS],
    color=whiten("green", 0.5),
    linewidth=0,
    label="Indoor masks",
)

ax1.fill_betweenx(
    [-10, 10],
        R_eff=R[-1],
        cov=cov,
        tau=tau,
    )

MASKS = np.datetime64('2021-12-31')
LOCKOUT = np.datetime64('2022-01-06')

fig1 = plt.figure(figsize=(10, 6))
ax1 = plt.axes()

ax1.fill_betweenx(
    [-10, 10],
    [MASKS, MASKS],
    [LOCKOUT, LOCKOUT],
    color=whiten("green", 0.5),
    linewidth=0,
    label="Indoor mask mandate",
)

ax1.fill_betweenx(
    [-10, 10],
    [LOCKOUT, LOCKOUT],
    [END_PLOT, END_PLOT],
    color=whiten("yellow", 0.5),
    linewidth=0,
    label="Unvaccinated lockout",
)

ax1.fill_between(
    dates[1:] + 1,

LOCKDOWN = np.datetime64('2021-08-13')
END_LOCKDOWN = np.datetime64('2021-10-15')
FURTHER_EASING = np.datetime64('2021-10-22')
EASING_95  = np.datetime64('2021-11-12')

fig1 = plt.figure(figsize=(10, 6))
ax1 = plt.axes()


ax1.fill_betweenx(
    [-10, 10],
    [LOCKDOWN, LOCKDOWN],
    [END_LOCKDOWN, END_LOCKDOWN],
    color=whiten("red", 0.45),
    linewidth=0,
    label="Lockdown",
)

ax1.fill_betweenx(
    [-10, 10],
    [END_LOCKDOWN, END_LOCKDOWN],
    [FURTHER_EASING, FURTHER_EASING],
    color=whiten("yellow", 0.5),
    linewidth=0,
    label="Easing",
)

ax1.fill_betweenx(
    [-10, 10],
END_LOCKDOWN = np.datetime64('2021-10-11')
EASING_80 = np.datetime64('2021-10-18')
END_MASKS = np.datetime64('2021-12-15')
MASKS_AGAIN = np.datetime64('2021-12-24')
DENSITY_LIMITS = np.datetime64('2021-12-27')
END_DENSITY_LIMITS = np.datetime64('2022-02-18')
END_MASKS_AGAIN = np.datetime64('2022-02-25')

fig1 = plt.figure(figsize=(10, 6))
ax1 = plt.axes()

ax1.fill_betweenx(
    [-10, 10],
    [MASKS, MASKS],
    [LGA_LOCKDOWN, LGA_LOCKDOWN],
    color=whiten("yellow", 0.5),
    linewidth=0,
    label="Density limits/70% easing",
)

ax1.fill_betweenx(
    [-10, 10],
    [LGA_LOCKDOWN, LGA_LOCKDOWN],
    [LOCKDOWN, LOCKDOWN],
    color=whiten("yellow", 0.5),
    edgecolor=whiten("orange", 0.5),
    linewidth=0,
    hatch="//////",
    label="East Sydney LGA lockdown",
)
ALERT_LEVEL_1 = np.datetime64('2021-06-29')
ALERT_LEVEL_4 = np.datetime64('2021-08-17')
ALERT_LEVEL_3 = np.datetime64('2021-09-22')
STEP_1 = np.datetime64('2021-10-06')
STEP_2 = np.datetime64('2021-11-10')
TRAFFIC_LIGHT_SYSTEM = np.datetime64('2021-12-03')
# END_LOCKDOWN = all_dates[-1] + 28

fig1 = plt.figure(figsize=(10, 6))
ax1 = plt.axes()

ax1.fill_betweenx(
    [-10, 10],
    [ALERT_LEVEL_4, ALERT_LEVEL_4],
    [ALERT_LEVEL_3, ALERT_LEVEL_3],
    color=whiten("#F27824", 0.5),
    # alpha=0.45,
    linewidth=0,
    label="Alert Level 4",
)

ax1.fill_betweenx(
    [-10, 10],
    [ALERT_LEVEL_3, ALERT_LEVEL_3],
    [STEP_1, STEP_1],
    color=whiten("#F6AE2F", 0.5),
    # alpha=0.45,
    linewidth=0,
    label="Alert Level 3",
)