def load_MACHO(macho_id):
    """lightcurve of 2 bands (R, B) from the MACHO survey.

    Notes
    -----

    The files are gathered from the original FATS project tutorial:
    https://github.com/isadoranun/tsfeat

    """
    tarfname = "{}.tar.bz2".format(macho_id)
    tarpath = os.path.join(DATA_PATH, tarfname)

    rpath = "{}.R.mjd".format(macho_id)
    bpath = "{}.B.mjd".format(macho_id)
    with tarfile.open(tarpath, mode="r:bz2") as tf:
        rlc = np.loadtxt(tf.extractfile(rpath))
        blc = np.loadtxt(tf.extractfile(bpath))

    bands = ("R", "B")
    data = {
        "R": {
            "time": rlc[:, 0],
            "magnitude": rlc[:, 1],
            "error": rlc[:, 2]},
        "B": {
            "time": blc[:, 0],
            "magnitude": blc[:, 1],
            "error": blc[:, 2]}
    }
    descr = ("The files are gathered from the original FATS project "
             "tutorial: https://github.com/isadoranun/tsfeat")

    return Data(
        id=macho_id, metadata=None, ds_name="MACHO",
        description=descr, bands=bands, data=data)