#include<bits/stdc++.h>
using namespace std;
#define FOR(i,s,t) for(int i=(s),_t=(t); i<=_t; ++i)
#define DOR(i,s,t) for(int i=(s),_t=(t); i>=_t; --i)
#define int long long 
typedef long long ll;
const int N=1e5+50;
const int Mod=1e9+7;
ll Fast(ll x,int b) {
    ll t=1;
    for(; b; b>>=1,x=x*x%Mod) if(b&1) t=t*x%Mod;
    return t;
}
vector<pair<int,int>> E[N];
int n,K,L;
int sz[N],mx[N];
int mark[N];
vector<int> vec;
void dfs1(int x,int fr) {
    vec.push_back(x);
    sz[x]=mx[x]=0;
    for(auto [y,v]: E[x]) if(y!=fr && !mark[y]) {
        dfs1(y,x);
        sz[x]+=sz[y];
        mx[x]=max(mx[x],sz[y]);
    }
}
vector<ll> f;
vector<pair<int,ll>> g;
void dfs2(int x,int fr,ll d) {
    g.push_back({x,d});
    f.push_back(d);
    for(auto [y,v]: E[x]) if(y!=fr && !mark[y]) {
        dfs2(y,x,d+v);
    }
}
int cnt[N],val[N];
void divide(int x) {
    vec.clear();
    dfs1(x,0);
    int rt=x;
    for(auto y: vec) {
        if(max(sz[x]-sz[y],mx[y])<max(sz[x]-sz[rt],mx[rt])) {
            rt=y;
        }
    }
    mark[rt]=1;

    g.clear();
    f.clear();
    dfs2(rt,0,0);
    sort(f.begin(),f.end());
    for(auto [x,d]: g) {
        cnt[x]+=upper_bound(f.begin(),f.end(),L-d)-f.begin();
    }
    for(auto [y,w]: E[rt]) if(!mark[y]) {
        g.clear();
        f.clear();
        dfs2(y,rt,w);
        sort(f.begin(),f.end());
        for(auto [x,d]: g) {
            cnt[x]-=upper_bound(f.begin(),f.end(),L-d)-f.begin();
        }
    }
    for(auto [y,w]: E[rt]) if(!mark[y]) divide(y);
}
ll Fac[N],Inv[N];
ll C(int n,int m) {
    if(n<m) return 0;
    return Fac[n]*Inv[m]%Mod*Inv[n-m]%Mod;
}
ll dis[N],Fa[N][21];
void dfs3(int x,int fr) {
    Fa[x][0]=fr;
    for(auto [y,w]: E[x]) if(y!=fr) {
        dis[y]=dis[x]+w;
        dfs3(y,x);
    }
}
signed main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cin >> n >> K >> L;
    FOR(i,2,n) {
        int u,v,w;
        cin >> u >> v >> w;
        E[u].push_back({v,w});
        E[v].push_back({u,w});
    }
    divide(1);
    dfs3(1,0);
    FOR(j,1,20) FOR(i,1,n) Fa[i][j]=Fa[Fa[i][j-1]][j-1];
    FOR(i,1,n) {
        int x=i;
        DOR(j,20,0) if(Fa[x][j] && dis[i]-dis[Fa[x][j]]<=L) x=Fa[x][j];
        ++val[x];
    }
    Fac[0]=1;
    FOR(i,1,n) Fac[i]=Fac[i-1]*i%Mod;
    Inv[n]=Fast(Fac[n],Mod-2);
    DOR(i,n,1) Inv[i-1]=Inv[i]*i%Mod;

    ll res=0;
    FOR(i,1,n) {
        // cerr << cnt[i] << ' ' << val[i] << endl;
        res=(res+C(cnt[i],K)-C(cnt[i]-val[i],K))%Mod;
    }
    if(res<0) res+=Mod;
    res=res*Fac[K]%Mod;
    cout << res << '\n';
    return 0;
}