def METHOD_NAME():
    parfile_ell1k = StringIO(
        """
        # Based on J0636+5128. Only for testing. Relativistic 
        # effects may not be physically consistent.
        PSRJ            PSRTEST
        ELONG           96.363147518                  1.400e-08
        ELAT            28.24309900                   3.000e-08
        DM              11.108159
        PEPOCH          57277.00
        F0              348.5592316999902             9.000e-13
        F1              -4.1895E-16                   8.000e-20
        POSEPOCH        57277.00
        DMEPOCH         56307.00
        BINARY          ELL1k
        A1              0.0898636            1         6.000e-08
        TASC            57277.01594431       1         1.100e-07
        EPS1            1.4E-6               1         1.000e-05
        EPS2            1.7E-5               1         9.000e-06
        OMDOT           0                    1
        LNEDOT          0                    1
        CLK             TT(BIPM2017)
        EPHEM           DE436
        RM              -7                            1.000e+00
        PX              1.4                           3.000e-01
        FB0             0.00017391195942     1         4.000e-14
        FB1             -7.7E-20             1         3.000e-21
        PMELONG         3.33                          4.000e-02
        PMELAT          -1.44                         1.100e-01
        UNITS           TDB
        """
    )
    model_ell1k = get_model(parfile_ell1k)
    parfile_ell1 = StringIO(
        """
        # Based on J0636+5128. Only for testing. Relativistic 
        # effects may not be physically consistent.
        PSRJ            PSRTEST
        ELONG           96.363147518                  1.400e-08
        ELAT            28.24309900                   3.000e-08
        DM              11.108159
        PEPOCH          57277.00
        F0              348.5592316999902             9.000e-13
        F1              -4.1895E-16                   8.000e-20
        POSEPOCH        57277.00
        DMEPOCH         56307.00
        BINARY          ELL1
        A1              0.0898636            1         6.000e-08
        TASC            57277.01594431       1         1.100e-07
        EPS1            1.4E-6               1         1.000e-05
        EPS2            1.7E-5               1         9.000e-06
        CLK             TT(BIPM2017)
        EPHEM           DE436
        RM              -7                            1.000e+00
        PX              1.4                           3.000e-01
        FB0             0.00017391195942     1         4.000e-14
        FB1             -7.7E-20             1         3.000e-21
        PMELONG         3.33                          4.000e-02
        PMELAT          -1.44                         1.100e-01
        UNITS           TDB
        """
    )
    model_ell1 = get_model(parfile_ell1)
    fake_toas = make_fake_toas_uniform(
        startMJD=50000, endMJD=51000, ntoas=100, model=model_ell1k, add_noise=True
    )
    return model_ell1, model_ell1k, fake_toas