Fork me on GitHub

source: git/external/fastjet/contribs/Nsubjettiness/Nsubjettiness.cc@ 239e1d0

ImprovedOutputFile Timing dual_readout llp 3.3.3pre11
Last change on this file since 239e1d0 was 973b92a, checked in by Pavel Demin <pavel.demin@…>, 9 years ago

update FastJet library to 3.1.3 and Nsubjettiness library to 2.2.1

  • Property mode set to 100644
File size: 2.0 KB
Line 
1// Nsubjettiness Package
2// Questions/Comments? jthaler@jthaler.net
3//
4// Copyright (c) 2011-14
5// Jesse Thaler, Ken Van Tilburg, Christopher K. Vermilion, and TJ Wilkason
6//
7// $Id: Nsubjettiness.cc 821 2015-06-15 18:50:53Z jthaler $
8//----------------------------------------------------------------------
9// This file is part of FastJet contrib.
10//
11// It is free software; you can redistribute it and/or modify it under
12// the terms of the GNU General Public License as published by the
13// Free Software Foundation; either version 2 of the License, or (at
14// your option) any later version.
15//
16// It is distributed in the hope that it will be useful, but WITHOUT
17// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
18// or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public
19// License for more details.
20//
21// You should have received a copy of the GNU General Public License
22// along with this code. If not, see <http://www.gnu.org/licenses/>.
23//----------------------------------------------------------------------
24
25#include "Nsubjettiness.hh"
26
27FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh
28
29namespace contrib {
30
31LimitedWarning Nsubjettiness::_old_constructor_warning;
32
33
34//result returns tau_N with normalization dependent on what is specified in constructor
35double Nsubjettiness::result(const PseudoJet& jet) const {
36 std::vector<fastjet::PseudoJet> particles = jet.constituents();
37 return _njettinessFinder.getTau(_N, particles);
38}
39
40TauComponents Nsubjettiness::component_result(const PseudoJet& jet) const {
41 std::vector<fastjet::PseudoJet> particles = jet.constituents();
42 return _njettinessFinder.getTauComponents(_N, particles);
43}
44
45//ratio result uses Nsubjettiness result to find the ratio tau_N/tau_M, where N and M are specified by user
46double NsubjettinessRatio::result(const PseudoJet& jet) const {
47 double numerator = _nsub_numerator.result(jet);
48 double denominator = _nsub_denominator.result(jet);
49 return numerator/denominator;
50}
51
52} // namespace contrib
53
54FASTJET_END_NAMESPACE
Note: See TracBrowser for help on using the repository browser.