def METHOD_NAME():
    # start dev
    logger.info('cells_filter_columns: %s' % par["cells_filter_columns"])
    logger.info('min_cells_filter_columns: %s' % par["min_cells_filter_columns"])
    # end dev
    census_connection = connect_census(
        par["input_database"],
        par["cellxgene_release"]
        ) 
    query_data = get_anndata(
        census_connection,
        par["cell_query"],
        par["species"]
        )
    query_data.obs = add_cellcensus_metadata_obs(
        census_connection,
        query_data
        )
    census_connection.close()
    del census_connection
    if par["cells_filter_columns"]:
        if not par["min_cells_filter_columns"]:
            raise NotImplementedError(
            "You specified cells_filter_columns, thus add min_cells_filter_columns!"
            )
        query_data = cellcensus_cell_filter(
            query_data,
            par["cells_filter_columns"],
            par["min_cells_filter_columns"]
            )
    mdata = mu.MuData(
        {par["modality"]: query_data}
        )
    mdata["rna"].var_names = mdata["rna"].var["feature_id"]
    mdata["rna"].var["gene_symbol"] = mdata["rna"].var["feature_name"]
    write_mudata(
        mdata,
        par["output"],
        par["output_compression"]
        )