Page MenuHomeHEPForge
Files F10227126

SoftDrop.cc
No OneTemporary
Actions

SoftDrop.cc
View Options

// $Id: SoftDrop.cc 1059 2017-09-07 20:48:00Z gsoyez $
//
// Copyright (c) 2014-, Gregory Soyez, Jesse. Thaler
// based on arXiv:1402.2657 by Andrew J. Larkoski, Simone Marzani,
// Gregory Soyez, Jesse Thaler
//
//----------------------------------------------------------------------
// 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 "SoftDrop.hh"
#include <cmath>
#include <sstream>
using namespace std;
FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh
namespace contrib{
//----------------------------------------------------------------------
// TODO:
//
// - implement reclustering (at the moment we assume it's C/A as for mMDT
//
// - what is returned if no substructure is found? [and for negativeve
// pt, m2 or other situations where mMDT currentlty returns an empty
// PseudoJet]
//
// GS.: mMDT seeks substructure so it makes sense for it to return
// an empry PseudoJet when no substructure is found (+ it is the
// original behaviour). For SoftDrop, in grooming mode (beta>0), if
// would make sense to return a jet with a single particle. In
// tagging mode (beta<0), the situation is less clear. At the level
// of the implementation, having a virtual function could work
// (with a bit of care to cover the -ve pt or m2 cases)
//
// - Do we include Andrew and Simone in the "contrib-author" list
// since they are on the paper? Do we include Gavin in the author's
// list since he started this contrib?
//
//----------------------------------------------------------------------
//----------------------------------------------------------------------
double SoftDrop::symmetry_cut_fn(const PseudoJet & p1,
const PseudoJet & p2,
void *optional_R0sqr_ptr) const{
double R0sqr = (optional_R0sqr_ptr == 0) ? _R0sqr : *((double*) optional_R0sqr_ptr);
return _symmetry_cut * pow(squared_geometric_distance(p1,p2)/R0sqr, 0.5*_beta);
}
//----------------------------------------------------------------------
string SoftDrop::symmetry_cut_description() const {
ostringstream ostr;
ostr << _symmetry_cut << " (theta/" << sqrt(_R0sqr) << ")^" << _beta << " [SoftDrop]";
return ostr.str();
}
} // namespace contrib
FASTJET_END_NAMESPACE

File Metadata

Mime Type
text/x-c
Expires
Thu, Apr 3, 8:01 PM (2 h, 49 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4724294
Default Alt Text
SoftDrop.cc (2 KB)

Event Timeline