#include<bits/stdc++.h>
using namespace std;
const int MAXN=3e5+5,MOD=1e9+7;
#define ll long long
int n,m,k,siz[MAXN],son[MAXN],dep[MAXN],fa[MAXN];
int dfn[MAXN],rnk[MAXN],top[MAXN],cnt;
int s[MAXN],lc[MAXN];
ll fac[MAXN],ans;
vector<int> ve[MAXN];
void dfs(int o,int pa)
{
	fa[o]=pa;
	siz[o]=1;dep[o]=dep[fa[o]]+1;
	for(auto u:ve[o])
	{
		if(u==pa)continue;
		dfs(u,o);siz[o]+=siz[u];
		if(siz[u]>siz[son[o]])son[o]=u;
	}
}
void dfs1(int o,int tp)
{
	top[o]=tp;dfn[o]=++cnt;rnk[cnt]=o;
	if(son[o])dfs1(son[o],tp);
	for(auto u:ve[o])
	{
		if(u==son[o]||u==fa[o])continue;
		dfs1(u,u);
	}
}
int lca(int x,int y)
{
	while(top[x]^top[y])
	{
		if(dep[top[x]]<dep[top[y]])swap(x,y);
		x=fa[top[x]];
	}
	if(dep[x]>dep[y])swap(x,y);
	return x;
}
ll qpow(ll a,ll b)
{
	ll ans=1;
	while(b)
	{
		if(b&1)ans*=a,ans%=MOD;
		a*=a,a%=MOD;b>>=1;
	}
	return ans;
}
ll cbn(ll x,ll y)
{
	if(y>x)return 0;
	return fac[x]*qpow(fac[y]*fac[x-y]%MOD,MOD-2)%MOD;
}
void calc(int o)
{
	for(auto u:ve[o])
	{
		if(u==fa[o])continue;
		calc(u);s[o]+=s[u];
	}
	ans+=cbn(s[o],k)-cbn(s[o]-lc[o],k);ans%=MOD;
}
int main()
{
	int T;cin>>T;fac[0]=1;
	for(int i=1;i<=3e5;i++)fac[i]=fac[i-1]*i%MOD;
	while(T--)
	{
		cin>>n>>m>>k;
		for(int i=1;i<n;i++)
		{
			int u,v;
			scanf("%d%d",&u,&v);
			ve[u].push_back(v);
			ve[v].push_back(u);
		}
		dfs(1,1);
		dfs1(1,1);fa[1]=0;
		for(int i=1;i<=m;i++)
		{
			int u,v,c;
			scanf("%d%d",&u,&v);c=lca(u,v);
			s[u]++,s[v]++;s[c]--,s[fa[c]]--;
			lc[c]++;
		}
		calc(1);
		cout<<(ans+MOD)%MOD<<endl;
		for(int i=1;i<=n;i++)ve[i].clear(),s[i]=lc[i]=0;
		for(int i=1;i<=n;i++)
			son[i]=0;cnt=ans=0;
	}
	return 0;
}