Page MenuHomeHEPForge
Files F8310297

WinnerTakeAllRecombiner.cc
No OneTemporary
Actions

WinnerTakeAllRecombiner.cc
View Options

// Nsubjettiness Package
// Questions/Comments? jthaler@jthaler.net
//
// Copyright (c) 2011-14
// Jesse Thaler, Ken Van Tilburg, Christopher K. Vermilion, and TJ Wilkason
//
// $Id: WinnerTakeAllRecombiner.cc 811 2015-05-27 04:31:47Z tjwilk $
//----------------------------------------------------------------------
// This file is part of FastJet contrib.
//
// It 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.
//
// It 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 this code. If not, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------
#include "WinnerTakeAllRecombiner.hh"
FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh
namespace contrib{
std::string GeneralEtSchemeRecombiner::description() const {
return "General Et-scheme recombination";
}
// recombine pa and pb according to a generalized Et-scheme parameterized by the power delta
void GeneralEtSchemeRecombiner::recombine(const fastjet::PseudoJet & pa, const fastjet::PseudoJet & pb, fastjet::PseudoJet & pab) const {
// Define new weights for recombination according to delta
// definition of ratio done so that we do not encounter issues about numbers being too large for huge values of delta
double ratio;
if (abs(_delta - 1.0) < std::numeric_limits<double>::epsilon()) ratio = pb.perp()/pa.perp(); // save computation time of pow()
else ratio = pow(pb.perp()/pa.perp(), _delta);
double weighta = 1.0/(1.0 + ratio);
double weightb = 1.0/(1.0 + 1.0/ratio);
double perp_ab = pa.perp() + pb.perp();
// reweight the phi and rap sums according to the weights above
if (perp_ab != 0.0) {
double y_ab = (weighta * pa.rap() + weightb * pb.rap());
double phi_a = pa.phi(), phi_b = pb.phi();
if (phi_a - phi_b > pi) phi_b += twopi;
if (phi_a - phi_b < -pi) phi_b -= twopi;
double phi_ab = (weighta * phi_a + weightb * phi_b);
pab.reset_PtYPhiM(perp_ab, y_ab, phi_ab);
}
else {
pab.reset(0.0,0.0,0.0,0.0);
}
}
std::string WinnerTakeAllRecombiner::description() const {
return "Winner-Take-All recombination";
}
// recombine pa and pb by creating pab with energy of the sum of particle energies in the direction of the harder particle
// updated recombiner to use more general form of a metric equal to E*(pT/E)^(alpha), which reduces to pT*cosh(rap)^(1-alpha)
// alpha is specified by the user. The default is alpha = 1, which is the typical behavior. alpha = 2 provides a metric which more
// favors central jets
void WinnerTakeAllRecombiner::recombine(const fastjet::PseudoJet & pa, const fastjet::PseudoJet & pb, fastjet::PseudoJet & pab) const {
double a_pt = pa.perp(), b_pt = pb.perp(), a_rap = pa.rap(), b_rap = pb.rap();
// special case of alpha = 1, everything is just pt (made separate so that pow function isn't called)
if (_alpha == 1.0) {
if (a_pt >= b_pt) {
pab.reset_PtYPhiM(a_pt + b_pt, a_rap, pa.phi());
}
else if (b_pt > a_pt) {
pab.reset_PtYPhiM(a_pt + b_pt, b_rap, pb.phi());
}
}
// every other case uses additional cosh(rap) term
else {
double a_metric = a_pt*pow(cosh(a_rap), 1.0-_alpha);
double b_metric = b_pt*pow(cosh(b_rap), 1.0-_alpha);
if (a_metric >= b_metric) {
double new_pt = a_pt + b_pt*pow(cosh(b_rap)/cosh(a_rap), 1.0-_alpha);
pab.reset_PtYPhiM(new_pt, a_rap, pa.phi());
}
if (b_metric > a_metric) {
double new_pt = b_pt + a_pt*pow(cosh(a_rap)/cosh(b_rap), 1.0-_alpha);
pab.reset_PtYPhiM(new_pt, b_rap, pb.phi());
}
}
}
} //namespace contrib
FASTJET_END_NAMESPACE

File Metadata

Mime Type
text/x-c
Expires
Sat, Dec 21, 6:00 PM (9 h, 54 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4023729
Default Alt Text
WinnerTakeAllRecombiner.cc (4 KB)

Event Timeline