float ReconstructionMCMatching::DrH(Interpretation& i){
    float dr=999;
    if(higgs.Pt()>0.01 && i.Higgs().Pt()>0.01){
	dr=higgs.DeltaR(i.Higgs());
    }
    return dr;
}
float ReconstructionQuality::TTH_ME(Interpretation& i){
  if(i.HasTag("TTH_ME")) return i.GetTag("TTH_ME");
  float tag=1e4*me.GetTTHMEsq(i.TopHad(),i.TopLep(),i.Higgs());
  i.SetTag("TTH_ME",tag);
  return tag;
}