31 #include <fastjet/tools/TopTaggerBase.hh>
33 FASTJET_BEGIN_NAMESPACE
52 double TopTaggerBase::_cos_theta_W(
const PseudoJet & res)
const{
55 vector<PseudoJet> W_pieces = W.
pieces();
56 assert(W_pieces.size() == 2);
60 PseudoJet W2 = (W_pieces[0].perp2() < W_pieces[1].perp2())
69 return (W2.px()*top.px() + W2.py()*top.py() + W2.pz()*top.pz())/
70 sqrt(W2.modp2() * top.
modp2());