def METHOD_NAME(
    self,
    sage_with_planar_monovalent_carbonyl,
    carbonyl_planar,
    distance_,
    theta,
):
    sage_with_planar_monovalent_carbonyl["VirtualSites"].parameters[
        0
    ].distance = Quantity(
        distance_,
        unit.angstrom,
    )
    sage_with_planar_monovalent_carbonyl["VirtualSites"].parameters[
        0
    ].inPlaneAngle = Quantity(
        theta,
        unit.degree,
    )
    out = sage_with_planar_monovalent_carbonyl.create_interchange(
        carbonyl_planar.to_topology(),
    )
    positions = get_positions_with_virtual_sites(out).to(unit.angstrom)
    distance = numpy.linalg.norm((positions[-1] - positions[0]).m_as(unit.angstrom))
    assert distance == pytest.approx(distance_)