docs: def getIcis_geno(self,geneID,cis_window=50E3):
          """ if eqtl==True it returns a bool vec for cis """
          assert self.eqtl == True, 'Only for eqtl data'
          index = self.geneID==geneID
          [_chrom,_gene_start,_gene_end] = self.gene_pos[index][0,:]
          Icis = (self.genoChrom==_chrom)*(self.genoPos>=_gene_start-cis_window)*(self.genoPos<=_gene_end+cis_window)
          return Icis