#include <bits/stdc++.h>

#define x first
#define y second

using namespace std;

typedef double db;
typedef pair<int, int> PII;

const int N = 110;
const db PI = acos(-1);

int n, sx, sy, d, t, order[N];
PII cv[N];
db rads[N], alpha[N], rr[2];

void rotate90(int &x, int &y) {
	int t = y;
	y = x, x = -t;
}
db rad(int x, int y) {
	return atan2(y, x);
}
bool check() {
	for(int i = 1; i <= n; i ++) {
		if(cv[i].x <= 0 && cv[i].y <= 0) return false;
	}
	return true;
}
db dist(int x, int y) {
	return sqrt(x*x + y*y);
}
db shift(db r) {
	db ret = r + 3./4 * PI;
	ret = ret - int(ret/2/PI)*2*PI;
	if(ret >= 2*PI) ret -= 2*PI;
	if(ret < 0) ret += 2*PI;
	return ret;
}
int main() {
	ios::sync_with_stdio(0); cin.tie(0);
	cout << fixed << setprecision(12);
	cin >> n >> sx >> sy >> d >> t;
	int cnt = 1.*t/2/PI;
	db res = t - 1.*cnt*2*PI;
	for(int i = 1; i <= n; i ++) cin >> cv[i].x >> cv[i].y;
	for(int u = 0; u < 3; u ++) {
		if(check()) break;
		rotate90(sx, sy);
		for(int i = 1; i <= n; i ++) rotate90(cv[i].x, cv[i].y);
	}
	rads[0] = rad(sx, sy);
	iota(order+1, order+n+1, 1);
	for(int i = 1; i <= n; i ++) rads[i] = rad(cv[i].x, cv[i].y);
	for(int i = 1; i <= n; i ++) alpha[i] = asin(1.*d/dist(cv[i].x, cv[i].y));
	sort(order+1, order+n+1, [](auto x, auto y) {
		return shift(rads[x]-alpha[x]) < shift(rads[y]-alpha[y]);
	});
	rr[0] = shift(rads[order[1]] - alpha[order[1]]);
	sort(order+1, order+n+1, [](auto x, auto y) {
		return shift(rads[x]+alpha[x]) < shift(rads[y]+alpha[y]);
	});
	rr[1] = shift(rads[order[n]] + alpha[order[n]]);
	// cerr << rr[0] << ' ' << rr[1] << '\n';
	db ans = 1.*cnt * (rr[1]-rr[0]);
	// cerr << ans << '\n';
	rads[0] = shift(rads[0]);
	if(rads[0] >= rr[0] && rads[0] <= rr[1]) {
		ans += min(rr[1]-rads[0], res);
	}
	else {
		if(rads[0] > rr[1] && res >= 2*PI-rads[0] + rr[0]) {
			res -= 2*PI-rads[0] + rr[0];
			ans += min(rr[1]-rr[0], res);
		}
		if(rads[0] < rr[0] && res >= rr[0]-rads[0]) {
			res -= rr[0]-rads[0];
			ans += min(rr[1]-rr[0], res);
		}
	}
	cout << ans << '\n';
	return 0;
}