Loading...
Searching...
No Matches
TopTaggerBase.cc
//FJSTARTHEADER
// $Id$
//
// Copyright (c) 2005-2021, Matteo Cacciari, Gavin P. Salam and Gregory Soyez
//
//----------------------------------------------------------------------
// This file is part of FastJet.
//
// FastJet is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// The algorithms that underlie FastJet have required considerable
// development. They are described in the original FastJet paper,
// hep-ph/0512210 and in the manual, arXiv:1111.6097. If you use
// FastJet as part of work towards a scientific publication, please
// quote the version you use and include a citation to the manual and
// optionally also to hep-ph/0512210.
//
// FastJet is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with FastJet. If not, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------
//FJENDHEADER
#include "fastjet/tools/TopTaggerBase.hh"
FASTJET_BEGIN_NAMESPACE
using namespace std;
// compute the W helicity angle
//
// The helicity angle is a standard observable in top decays, used to
// determine the Lorentz structure of the top- W coupling [13]. It is
// defined as the angle, measured in the rest frame of the
// reconstructed W, between the reconstructed top's flight direction
// and one of the W decay products. Normally, it is studied in
// semi-leptonic top decays, where the charge of the lepton uniquely
// identifies these decay products. In hadronic top decays there is an
// ambiguity which we resolve by choosing the lower pT subjet, as
// measured in the lab frame.
//
// The jet passed to this function is expected to already have
// the structure of a top, including a functional "W()" call;
// the W must be made of two pieces.
double TopTaggerBase::_cos_theta_W(const PseudoJet & res) const{
// the two jets of interest: top and lower-pt prong of W
const PseudoJet & W = res.structure_of<TopTaggerBase>().W();
vector<PseudoJet> W_pieces = W.pieces();
assert(W_pieces.size() == 2);
//assert(W_pieces[0].perp2() >= W_pieces[1].perp2());
//PseudoJet W2 = W_pieces[1];
// extract the softer of the two W pieces.
PseudoJet W2 = (W_pieces[0].perp2() < W_pieces[1].perp2())
? W_pieces[0]
: W_pieces[1];
PseudoJet top = res;
// transform these jets into jets in the rest frame of the W
W2.unboost(W);
top.unboost(W);
return (W2.px()*top.px() + W2.py()*top.py() + W2.pz()*top.pz())/
sqrt(W2.modp2() * top.modp2());
}
FASTJET_END_NAMESPACE
fastjet::PseudoJet
Class to contain pseudojets, including minimal information of use to jet-clustering routines.
Definition PseudoJet.hh:68
fastjet::PseudoJet::structure_of
const TransformerType::StructureType & structure_of() const
this is a helper to access any structure created by a Transformer (that is, of type Transformer::Stru...
Definition PseudoJet.hh:1122
fastjet::PseudoJet::modp2
double modp2() const
return the squared 3-vector modulus = px^2+py^2+pz^2
Definition PseudoJet.hh:178
fastjet::PseudoJet::pieces
virtual std::vector< PseudoJet > pieces() const
retrieve the pieces that make up the jet.
Definition PseudoJet.cc:792
fastjet::PseudoJet::unboost
PseudoJet & unboost(const PseudoJet &prest)
transform this jet (given in lab) into a jet in the rest frame of prest
Definition PseudoJet.cc:414
fastjet::TopTaggerBase::_cos_theta_W
double _cos_theta_W(const PseudoJet &result) const
computes the W helicity angle
Definition TopTaggerBase.cc:52