Better to do
```
vector<double> tfit( nptfit, 0.);
```
and get rid of the loop with ```tfit.push_back``` since the above will do the same work more efficiently (since the size of the vector will be properly allocated right at the beginning).