def METHOD_NAME(self):
    """
    Compare calculated completenesses for multiple targets with constrain orbits set to true
    """
    with RedirectStreams(stdout=self.dev_null):
        TL = TargetList(
            ntargs=100, constrainOrbits=True, **copy.deepcopy(self.spec)
        )
        mode = list(
            filter(
                lambda mode: mode["detectionMode"] == True,
                TL.OpticalSystem.observingModes,
            )
        )[0]
        IWA = mode["IWA"]
        OWA = mode["OWA"]
        rrange = TL.PlanetPopulation.rrange
        maxd = (rrange[1] / np.tan(IWA)).to(u.pc).value
        mind = (rrange[0] / np.tan(OWA)).to(u.pc).value
        # want distances to span from outer edge below IWA to inner edge above OWA
        TL.dist = (
            np.logspace(np.log10(mind / 10.0), np.log10(maxd * 10.0), TL.nStars)
            * u.pc
        )
    Brown = EXOSIMS.Completeness.BrownCompleteness.BrownCompleteness(
        constrainOrbits=True, **copy.deepcopy(self.spec)
    )
    Garrett = EXOSIMS.Completeness.GarrettCompleteness.GarrettCompleteness(
        constrainOrbits=True, **copy.deepcopy(self.spec)
    )
    cBrown = Brown.target_completeness(TL)
    cGarrett = Garrett.target_completeness(TL)
    np.testing.assert_allclose(cGarrett, cBrown, rtol=0.1, atol=1e-6)