Your problem is that f and tpas are both range variables. Range variables are NOT vectors. They are iterators (i.e. loop variables). That means TF1, for example, is calculated for every value of f AND every value of tpas. You need to make f and tpas vectors (as a good rule of thumb, only use range variable to index vectors. Use vectors for everything else).
