aboutsummaryrefslogtreecommitdiff
# Copyright (C) University of Tennessee Health Science Center, Memphis, TN.
#
# This program is free software: you can redistribute it and/or modify it
# under the terms of the GNU Affero General Public License
# as published by the Free Software Foundation, either version 3 of the
# License, or (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
# See the GNU Affero General Public License for more details.
#
# This program is available from Source Forge: at GeneNetwork Project
# (sourceforge.net/projects/genenetwork/).
#
# Contact Dr. Robert W. Williams at rwilliams@uthsc.edu
#
#
# This module is used by GeneNetwork project (www.genenetwork.org)
from gn2.wqflask.database import database_connection
from gn2.utility import helper_functions
from gn2.utility.tools import get_setting


class SendToGeneWeaver:
    def __init__(self, start_vars):
        trait_db_list = [trait.strip()
                         for trait in start_vars['trait_list'].split(',')]
        helper_functions.get_trait_db_obs(self, trait_db_list)

        self.chip_name = test_chip(self.trait_list)
        self.wrong_input = "False"
        if self.chip_name == "mixed" or self.chip_name == "not_microarray" or '_NA' in self.chip_name:
            self.wrong_input = "True"
        else:
            species = self.trait_list[0][1].group.species
            if species == "rat":
                species_name = "Rattus norvegicus"
            elif species == "human":
                species_name = "Homo sapiens"
            elif species == "mouse":
                species_name = "Mus musculus"
            else:
                species_name = ""

            trait_name_list = get_trait_name_list(self.trait_list)

            self.hidden_vars = {
                'client': "genenetwork",
                'species': species_name,
                'idtype': self.chip_name,
                'list': ",".join(trait_name_list),
            }


def get_trait_name_list(trait_list):
    name_list = []
    for trait_db in trait_list:
        name_list.append(trait_db[0].name)

    return name_list


def test_chip(trait_list):
    final_chip_name = ""
    with database_connection(get_setting("SQL_URI")) as conn, conn.cursor() as cursor:
        for trait_db in trait_list:
            dataset = trait_db[1]
            cursor.execute(
                "SELECT GeneChip.GO_tree_value "
                "FROM GeneChip, ProbeFreeze, ProbeSetFreeze "
                "WHERE GeneChip.Id = ProbeFreeze.ChipId "
                "AND ProbeSetFreeze.ProbeFreezeId = ProbeFreeze.Id "
                "AND ProbeSetFreeze.Name = %s",
                (dataset.name,)
            )

            if result := cursor.fetchone():
                chip_name = result[0]
                if chip_name:
                    if chip_name != final_chip_name:
                        if final_chip_name:
                            return "mixed"
                        else:
                            final_chip_name = chip_name
                    else:
                        pass
                else:
                    cursor.execute(
                        "SELECT GeneChip.Name "
                        "FROM GeneChip, ProbeFreeze, ProbeSetFreeze "
                        "WHERE GeneChip.Id = ProbeFreeze.ChipId "
                        "AND ProbeSetFreeze.ProbeFreezeId = ProbeFreeze.Id "
                        "AND ProbeSetFreeze.Name = %s",
                        (dataset.name,)
                    )
                    chip_name = f'{cursor.fetchone()[0]}_NA'
                    return chip_name
            else:
                cursor.execute(
                    "SELECT GeneChip.Name FROM GeneChip, "
                    "ProbeFreeze, ProbeSetFreeze WHERE "
                    "GeneChip.Id = ProbeFreeze.ChipId "
                    "AND ProbeSetFreeze.ProbeFreezeId = ProbeFreeze.Id "
                    "AND ProbeSetFreeze.Name = %s",
                    (dataset.name,)
                )
                if result := cursor.fetchone():
                    chip_name = f'{result[0]}_NA'
                    return chip_name
                return "not_microarray"

    return chip_name