Added plot script again, and upper limit is now reading from chi2 thresholds inside main