float ReconstructionMCMatching::DrTopHad(Interpretation& i){
    float dr=999;
    if(top_had.Pt()>0.01 && i.TopHad().Pt()>0.01){
	dr=top_had.DeltaR(i.TopHad());
    }
    return dr;
}
float ReconstructionQuality::TTHBB_ME(Interpretation& i){
  if(i.HasTag("TTHBB_ME")) return i.GetTag("TTHBB_ME");
  float tag=me.GetTTHBBMEsq(i.TopHad(),i.TopLep(),i.B1(),i.B2());
  i.SetTag("TTHBB_ME",tag);
  return tag;
}