def METHOD_NAME(device, q0_name, times, MC=None):
    if MC is None:
        MC = qc.station.components['MC']
    q0 = device.qubits()[q0_name]
    cal_points = 4
    cal_pts = times[-1] + \
        np.arange(1, 1+cal_points)*(times[1]-times[0])
    times = np.concatenate((times, cal_pts))
    operation_dict = device.get_operation_dict()
    AWG = q0.AWG
    busT1swf = awg_swf.awg_seq_swf(
        fsqs.BusT1,
        parameter_name='times',
        unit='s',
        AWG=q0.AWG,
        fluxing_channels=[q0.fluxing_channel()],
        awg_seq_func_kwargs={'operation_dict': operation_dict,
                             'q0': q0_name,
                             'distortion_dict': q0.dist_dict()})
    MC.set_sweep_function(busT1swf)
    MC.set_sweep_points(times)
    MC.set_detector_function(q0.int_avg_det_rot)
    MC.run('BusT1_{}'.format(q0.name))
    ma.T1_Analysis(label='BusT1')