Assuming you mean the most abundant known splice variant of each gene you could do the following. Download an existing RNA-seq data set from GEO for HeLa such as: GSM759888. Since a BAM is provided, you can skip the alignment step and go straight to estimating isoform expression levels.
If you want to do your own alignments (say to take advantage of a new aligner version), convert the BAM back to FASTQs using Picard SamToFastq. Then run Tophat to get new alignments.
Next get a GTF file representing all transcripts of all genes for human build hg19/build37. For example, you can get one for all Ensembl transcripts here: Homo_sapiens.GRCh37.64.gtf.gz
Use this GTF to run Cufflinks with the -G option. Once it finishes you should have an isoforms.fpkm file. This contains the FPKM expression estimate for each transcript (each of which should also be marked with a gene ID). With that data in hand it will be trivial to identify the most highly expressed transcript for each gene that has multiple known transcripts.
If you do not want to be limited to the known or Ensembl predicted alternative transcripts of each gene the problem becomes more complicated. You could attempt to merge transcript GTFs from multiple sources to get a more comprehensive representation of existing transcript annotations. Or you could run Cufflinks in de novo mode, then figure out which correspond to the same loci (using Cuffcompare perhaps), and select the one with the highest abundance for each locus.
Useful tips here on how to do the analysis but it seems that this has not been pre-calculated. I might put this forward for an MSc project.