fjcore.hh File Reference

#include <cstdlib>
#include <iostream>
#include <string>
#include <list>
#include <vector>
#include <valarray>
#include <cassert>
#include <cmath>
#include <limits>
#include <memory>
#include <map>
#include <set>

Include dependency graph for fjcore.hh:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	integral_type< T, _t >
struct	__inheritance_helper< B, D >
struct	IsBaseAndDerived< B, D >
struct	IsBaseAndDerived< B, D >::Host
class	SharedPtr< T >
class	SharedPtr< T >::__SharedCountingPtr
class	LimitedWarning
class	Error
class	InternalError
class	PseudoJetStructureBase
class	PseudoJet
class	PseudoJet::UserInfoBase
class	PseudoJet::InexistentUserInfo
class	IndexedSortHelper
class	FunctionOfPseudoJet< TOut >
class	SelectorWorker
class	Selector
class	Selector::InvalidWorker
class	Selector::InvalidArea
class	JetDefinition
class	JetDefinition::Recombiner
class	JetDefinition::DefaultRecombiner
class	JetDefinition::Plugin
class	CompositeJetStructure
class	ClusterSequenceStructure
class	ClusterSequence
class	ClusterSequence::Extras
struct	ClusterSequence::history_element
class	ClusterSequence::_Parabola
class	ClusterSequence::_Line
struct	ClusterSequence::BriefJet
class	ClusterSequence::TiledJet
struct	ClusterSequence::Tile
struct	ClusterSequence::EEBriefJet
class	_NoInfo
class	NNInfo< I >
class	NNInfo< _NoInfo >
class	NNBase< I >
class	NNH< BJ, I >
class	NNH< BJ, I >::NNBJ

Macros
#define	__FJCORE__
#define	__FJCORE_DROP_CGAL
#define	_INCLUDE_FJCORE_CONFIG_AUTO_H   1
#define	FJCORE_HAVE_DLFCN_H   1
#define	FJCORE_HAVE_INTTYPES_H   1
#define	FJCORE_HAVE_LIBM   1
#define	FJCORE_HAVE_MEMORY_H   1
#define	FJCORE_HAVE_STDINT_H   1
#define	FJCORE_HAVE_STDLIB_H   1
#define	FJCORE_HAVE_STRINGS_H   1
#define	FJCORE_HAVE_STRING_H   1
#define	FJCORE_HAVE_SYS_STAT_H   1
#define	FJCORE_HAVE_SYS_TYPES_H   1
#define	FJCORE_HAVE_UNISTD_H   1
#define	FJCORE_LT_OBJDIR   ".libs/"
#define	FJCORE_PACKAGE   "fastjet"
#define	FJCORE_PACKAGE_BUGREPORT   ""
#define	FJCORE_PACKAGE_NAME   "FastJet"
#define	FJCORE_PACKAGE_STRING   "FastJet 3.3.2"
#define	FJCORE_PACKAGE_TARNAME   "fastjet"
#define	FJCORE_PACKAGE_URL   ""
#define	FJCORE_PACKAGE_VERSION   "3.3.2"
#define	FJCORE_STDC_HEADERS   1
#define	FJCORE_VERSION   "3.3.2"
#define	FJCORE_VERSION_MAJOR   3
#define	FJCORE_VERSION_MINOR   3
#define	FJCORE_VERSION_NUMBER   30302
#define	FJCORE_VERSION_PATCHLEVEL   2
#define	__FJCORE_CONFIG_H__
#define	__FJCORE_FASTJET_BASE_HH__
#define	FJCORE_BEGIN_NAMESPACE   namespace fjcore {
#define	FJCORE_END_NAMESPACE   }
#define	FJCORE_OVERRIDE
#define	__FJCORE_NUMCONSTS__
#define	__FJCORE_INTERNAL_IS_BASE_HH__
#define	__FJCORE_FJCORE_DEPRECATED_HH__
#define	FJCORE_DEPRECATED
#define	FJCORE_DEPRECATED_MSG(message)
#define	__FJCORE_SHARED_PTR_HH__
#define	__FJCORE_LIMITEDWARNING_HH__
#define	__FJCORE_ERROR_HH__
#define	__FJCORE_PSEUDOJET_STRUCTURE_BASE_HH__
#define	__FJCORE_PSEUDOJET_HH__
#define	__FJCORE_FUNCTION_OF_PSEUDOJET_HH__
#define	__FJCORE_SELECTOR_HH__
#define	__FJCORE_JETDEFINITION_HH__
#define	__FJCORE_COMPOSITEJET_STRUCTURE_HH__
#define	__FJCORE_CLUSTER_SEQUENCE_STRUCTURE_HH__
#define	__FJCORE_CLUSTERSEQUENCE_HH__
#define	__FJCORE_NNBASE_HH__
#define	__FJCORE_NNH_HH__

Typedefs
typedef integral_type< bool, true >	true_type
typedef integral_type< bool, false >	false_type
typedef char(&	__yes_type)[1]
typedef char(&	__no_type)[2]
typedef JetAlgorithm	JetFinder

Enumerations
enum	Strategy {   N2MHTLazy9AntiKtSeparateGhosts = -10, N2MHTLazy9 = -7, N2MHTLazy25 = -6, N2MHTLazy9Alt = -5,   N2MinHeapTiled = -4, N2Tiled = -3, N2PoorTiled = -2, N2Plain = -1,   N3Dumb = 0, Best = 1, NlnN = 2, NlnN3pi = 3,   NlnN4pi = 4, NlnNCam4pi = 14, NlnNCam2pi2R = 13, NlnNCam = 12,   BestFJ30 = 21, plugin_strategy = 999 }
enum	JetAlgorithm {   kt_algorithm =0, cambridge_algorithm =1, antikt_algorithm =2, genkt_algorithm =3,   cambridge_for_passive_algorithm =11, genkt_for_passive_algorithm =13, ee_kt_algorithm =50, ee_genkt_algorithm =53,   plugin_algorithm = 99, undefined_jet_algorithm = 999 }
enum	RecombinationScheme {   E_scheme =0, pt_scheme =1, pt2_scheme =2, Et_scheme =3,   Et2_scheme =4, BIpt_scheme =5, BIpt2_scheme =6, WTA_pt_scheme =7,   WTA_modp_scheme =8, external_scheme = 99 }

Functions
template<class B , class D >
B *	cast_if_derived (D *d)
template<class T , class U >
bool	operator== (SharedPtr< T > const &t, SharedPtr< U > const &u)
template<class T , class U >
bool	operator!= (SharedPtr< T > const &t, SharedPtr< U > const &u)
template<class T , class U >
bool	operator< (SharedPtr< T > const &t, SharedPtr< U > const &u)
template<class T >
void	swap (SharedPtr< T > &a, SharedPtr< T > &b)
template<class T >
T *	get_pointer (SharedPtr< T > const &t)
PseudoJet	operator+ (const PseudoJet &, const PseudoJet &)
PseudoJet	operator- (const PseudoJet &, const PseudoJet &)
PseudoJet	operator* (double, const PseudoJet &)
PseudoJet	operator* (const PseudoJet &, double)
PseudoJet	operator/ (const PseudoJet &, double)
bool	operator== (const PseudoJet &, const PseudoJet &)
bool	operator!= (const PseudoJet &a, const PseudoJet &b)
bool	operator== (const PseudoJet &jet, const double val)
bool	operator== (const double val, const PseudoJet &jet)
bool	operator!= (const PseudoJet &a, const double val)
bool	operator!= (const double val, const PseudoJet &a)
double	dot_product (const PseudoJet &a, const PseudoJet &b)
double	cos_theta (const PseudoJet &a, const PseudoJet &b)
double	theta (const PseudoJet &a, const PseudoJet &b)
bool	have_same_momentum (const PseudoJet &, const PseudoJet &)
PseudoJet	PtYPhiM (double pt, double y, double phi, double m=0.0)
std::vector< PseudoJet >	sorted_by_pt (const std::vector< PseudoJet > &jets)
std::vector< PseudoJet >	sorted_by_rapidity (const std::vector< PseudoJet > &jets)
std::vector< PseudoJet >	sorted_by_E (const std::vector< PseudoJet > &jets)
std::vector< PseudoJet >	sorted_by_pz (const std::vector< PseudoJet > &jets)
void	sort_indices (std::vector< int > &indices, const std::vector< double > &values)
template<class T >
std::vector< T >	objects_sorted_by_values (const std::vector< T > &objects, const std::vector< double > &values)
PseudoJet	join (const std::vector< PseudoJet > &pieces)
PseudoJet	join (const PseudoJet &j1)
PseudoJet	join (const PseudoJet &j1, const PseudoJet &j2)
PseudoJet	join (const PseudoJet &j1, const PseudoJet &j2, const PseudoJet &j3)
PseudoJet	join (const PseudoJet &j1, const PseudoJet &j2, const PseudoJet &j3, const PseudoJet &j4)
Selector	SelectorIdentity ()
Selector	operator! (const Selector &s)
Selector	operator|| (const Selector &s1, const Selector &s2)
Selector	operator&& (const Selector &s1, const Selector &s2)
Selector	operator* (const Selector &s1, const Selector &s2)
Selector	SelectorPtMin (double ptmin)
Selector	SelectorPtMax (double ptmax)
Selector	SelectorPtRange (double ptmin, double ptmax)
Selector	SelectorEtMin (double Etmin)
Selector	SelectorEtMax (double Etmax)
Selector	SelectorEtRange (double Etmin, double Etmax)
Selector	SelectorEMin (double Emin)
Selector	SelectorEMax (double Emax)
Selector	SelectorERange (double Emin, double Emax)
Selector	SelectorMassMin (double Mmin)
Selector	SelectorMassMax (double Mmax)
Selector	SelectorMassRange (double Mmin, double Mmax)
Selector	SelectorRapMin (double rapmin)
Selector	SelectorRapMax (double rapmax)
Selector	SelectorRapRange (double rapmin, double rapmax)
Selector	SelectorAbsRapMin (double absrapmin)
Selector	SelectorAbsRapMax (double absrapmax)
Selector	SelectorAbsRapRange (double absrapmin, double absrapmax)
Selector	SelectorEtaMin (double etamin)
Selector	SelectorEtaMax (double etamax)
Selector	SelectorEtaRange (double etamin, double etamax)
Selector	SelectorAbsEtaMin (double absetamin)
Selector	SelectorAbsEtaMax (double absetamax)
Selector	SelectorAbsEtaRange (double absetamin, double absetamax)
Selector	SelectorPhiRange (double phimin, double phimax)
Selector	SelectorRapPhiRange (double rapmin, double rapmax, double phimin, double phimax)
Selector	SelectorNHardest (unsigned int n)
Selector	SelectorCircle (const double radius)
Selector	SelectorDoughnut (const double radius_in, const double radius_out)
Selector	SelectorStrip (const double half_width)
Selector	SelectorRectangle (const double half_rap_width, const double half_phi_width)
Selector	SelectorPtFractionMin (double fraction)
Selector	SelectorIsZero ()
FJCORE_BEGIN_NAMESPACE std::string	fastjet_version_string ()
PseudoJet	join (const std::vector< PseudoJet > &pieces, const JetDefinition::Recombiner &recombiner)
PseudoJet	join (const PseudoJet &j1, const JetDefinition::Recombiner &recombiner)
PseudoJet	join (const PseudoJet &j1, const PseudoJet &j2, const JetDefinition::Recombiner &recombiner)
PseudoJet	join (const PseudoJet &j1, const PseudoJet &j2, const PseudoJet &j3, const JetDefinition::Recombiner &recombiner)
PseudoJet	join (const PseudoJet &j1, const PseudoJet &j2, const PseudoJet &j3, const PseudoJet &j4, const JetDefinition::Recombiner &recombiner)

Variables
FJCORE_BEGIN_NAMESPACE const double	pi = 3.141592653589793238462643383279502884197
const double	twopi = 6.283185307179586476925286766559005768394
const double	pisq = 9.869604401089358618834490999876151135314
const double	zeta2 = 1.644934066848226436472415166646025189219
const double	zeta3 = 1.202056903159594285399738161511449990765
const double	eulergamma = 0.577215664901532860606512090082402431042
const double	ln2 = 0.693147180559945309417232121458176568076
FJCORE_BEGIN_NAMESPACE const double	MaxRap = 1e5
const double	pseudojet_invalid_phi = -100.0
const double	pseudojet_invalid_rap = -1e200
const JetAlgorithm	aachen_algorithm = cambridge_algorithm
const JetAlgorithm	cambridge_aachen_algorithm = cambridge_algorithm

Macro Definition Documentation

__FJCORE__

#define __FJCORE__

__FJCORE_CLUSTER_SEQUENCE_STRUCTURE_HH__

#define __FJCORE_CLUSTER_SEQUENCE_STRUCTURE_HH__

__FJCORE_CLUSTERSEQUENCE_HH__

#define __FJCORE_CLUSTERSEQUENCE_HH__

__FJCORE_COMPOSITEJET_STRUCTURE_HH__

#define __FJCORE_COMPOSITEJET_STRUCTURE_HH__

__FJCORE_CONFIG_H__

#define __FJCORE_CONFIG_H__

__FJCORE_DROP_CGAL

#define __FJCORE_DROP_CGAL

__FJCORE_ERROR_HH__

#define __FJCORE_ERROR_HH__

__FJCORE_FASTJET_BASE_HH__

#define __FJCORE_FASTJET_BASE_HH__

__FJCORE_FJCORE_DEPRECATED_HH__

#define __FJCORE_FJCORE_DEPRECATED_HH__

__FJCORE_FUNCTION_OF_PSEUDOJET_HH__

#define __FJCORE_FUNCTION_OF_PSEUDOJET_HH__

__FJCORE_INTERNAL_IS_BASE_HH__

#define __FJCORE_INTERNAL_IS_BASE_HH__

__FJCORE_JETDEFINITION_HH__

#define __FJCORE_JETDEFINITION_HH__

__FJCORE_LIMITEDWARNING_HH__

#define __FJCORE_LIMITEDWARNING_HH__

__FJCORE_NNBASE_HH__

#define __FJCORE_NNBASE_HH__

__FJCORE_NNH_HH__

#define __FJCORE_NNH_HH__

__FJCORE_NUMCONSTS__

#define __FJCORE_NUMCONSTS__

__FJCORE_PSEUDOJET_HH__

#define __FJCORE_PSEUDOJET_HH__

__FJCORE_PSEUDOJET_STRUCTURE_BASE_HH__

#define __FJCORE_PSEUDOJET_STRUCTURE_BASE_HH__

__FJCORE_SELECTOR_HH__

#define __FJCORE_SELECTOR_HH__

__FJCORE_SHARED_PTR_HH__

#define __FJCORE_SHARED_PTR_HH__

_INCLUDE_FJCORE_CONFIG_AUTO_H

#define _INCLUDE_FJCORE_CONFIG_AUTO_H   1

FJCORE_BEGIN_NAMESPACE

#define FJCORE_BEGIN_NAMESPACE   namespace fjcore {

FJCORE_DEPRECATED

#define FJCORE_DEPRECATED

FJCORE_DEPRECATED_MSG

#define FJCORE_DEPRECATED_MSG

(

message

)

FJCORE_END_NAMESPACE

#define FJCORE_END_NAMESPACE   }

FJCORE_HAVE_DLFCN_H

#define FJCORE_HAVE_DLFCN_H   1

FJCORE_HAVE_INTTYPES_H

#define FJCORE_HAVE_INTTYPES_H   1

FJCORE_HAVE_LIBM

#define FJCORE_HAVE_LIBM   1

FJCORE_HAVE_MEMORY_H

#define FJCORE_HAVE_MEMORY_H   1

FJCORE_HAVE_STDINT_H

#define FJCORE_HAVE_STDINT_H   1

FJCORE_HAVE_STDLIB_H

#define FJCORE_HAVE_STDLIB_H   1

FJCORE_HAVE_STRING_H

#define FJCORE_HAVE_STRING_H   1

FJCORE_HAVE_STRINGS_H

#define FJCORE_HAVE_STRINGS_H   1

FJCORE_HAVE_SYS_STAT_H

#define FJCORE_HAVE_SYS_STAT_H   1

FJCORE_HAVE_SYS_TYPES_H

#define FJCORE_HAVE_SYS_TYPES_H   1

FJCORE_HAVE_UNISTD_H

#define FJCORE_HAVE_UNISTD_H   1

FJCORE_LT_OBJDIR

#define FJCORE_LT_OBJDIR   ".libs/"

FJCORE_OVERRIDE

#define FJCORE_OVERRIDE

FJCORE_PACKAGE

#define FJCORE_PACKAGE   "fastjet"

FJCORE_PACKAGE_BUGREPORT

#define FJCORE_PACKAGE_BUGREPORT   ""

FJCORE_PACKAGE_NAME

#define FJCORE_PACKAGE_NAME   "FastJet"

FJCORE_PACKAGE_STRING

#define FJCORE_PACKAGE_STRING   "FastJet 3.3.2"

FJCORE_PACKAGE_TARNAME

#define FJCORE_PACKAGE_TARNAME   "fastjet"

FJCORE_PACKAGE_URL

#define FJCORE_PACKAGE_URL   ""

FJCORE_PACKAGE_VERSION

#define FJCORE_PACKAGE_VERSION   "3.3.2"

FJCORE_STDC_HEADERS

#define FJCORE_STDC_HEADERS   1

FJCORE_VERSION

#define FJCORE_VERSION   "3.3.2"

FJCORE_VERSION_MAJOR

#define FJCORE_VERSION_MAJOR   3

FJCORE_VERSION_MINOR

#define FJCORE_VERSION_MINOR   3

FJCORE_VERSION_NUMBER

#define FJCORE_VERSION_NUMBER   30302

FJCORE_VERSION_PATCHLEVEL

#define FJCORE_VERSION_PATCHLEVEL   2

Typedef Documentation

__no_type

typedef char(& __no_type)[2]

__yes_type

typedef char(& __yes_type)[1]

false_type

typedef integral_type<bool, false> false_type

the bool 'false' value promoted to a type

JetFinder

typedef JetAlgorithm JetFinder

true_type

typedef integral_type<bool, true> true_type

the bool 'true' value promoted to a type

Enumeration Type Documentation

JetAlgorithm

enum JetAlgorithm

Enumerator
kt_algorithm
cambridge_algorithm
antikt_algorithm
genkt_algorithm
cambridge_for_passive_algorithm
genkt_for_passive_algorithm
ee_kt_algorithm
ee_genkt_algorithm
plugin_algorithm
undefined_jet_algorithm

                  {
  kt_algorithm=0,
  cambridge_algorithm=1,
  antikt_algorithm=2,
  genkt_algorithm=3,
  cambridge_for_passive_algorithm=11,
  genkt_for_passive_algorithm=13,
  ee_kt_algorithm=50,
  ee_genkt_algorithm=53,
  plugin_algorithm = 99,
  undefined_jet_algorithm = 999
};

RecombinationScheme

enum RecombinationScheme

Enumerator
E_scheme
pt_scheme
pt2_scheme
Et_scheme
Et2_scheme
BIpt_scheme
BIpt2_scheme
WTA_pt_scheme
WTA_modp_scheme
external_scheme

                         {
  E_scheme=0,
  pt_scheme=1,
  pt2_scheme=2,
  Et_scheme=3,
  Et2_scheme=4,
  BIpt_scheme=5,
  BIpt2_scheme=6,
  WTA_pt_scheme=7,
  WTA_modp_scheme=8,
  external_scheme = 99
};

Strategy

enum Strategy

Enumerator
N2MHTLazy9AntiKtSeparateGhosts
N2MHTLazy9
N2MHTLazy25
N2MHTLazy9Alt
N2MinHeapTiled
N2Tiled
N2PoorTiled
N2Plain
N3Dumb
Best
NlnN
NlnN3pi
NlnN4pi
NlnNCam4pi
NlnNCam2pi2R
NlnNCam
BestFJ30
plugin_strategy

              {
  N2MHTLazy9AntiKtSeparateGhosts   = -10,
  N2MHTLazy9   = -7,
  N2MHTLazy25   = -6,
  N2MHTLazy9Alt   = -5,
  N2MinHeapTiled   = -4,
  N2Tiled     = -3,
  N2PoorTiled = -2,
  N2Plain     = -1,
  N3Dumb      =  0,
  Best        =  1,
  NlnN        =  2,
  NlnN3pi     =  3,
  NlnN4pi     =  4,
  NlnNCam4pi   = 14,
  NlnNCam2pi2R = 13,
  NlnNCam      = 12, // 2piMultD
  BestFJ30     =  21,
  plugin_strategy = 999
};

Function Documentation

cast_if_derived()

B* cast_if_derived

(

D *

d

)

                        {
  return IsBaseAndDerived<B,D>::value ? (B*)(d) : 0;
}

cos_theta()

double cos_theta ( const PseudoJet &  a, const PseudoJet &  b  )

inline

                                                                  {
  double dot_3d = a.px()*b.px() + a.py()*b.py() + a.pz()*b.pz();
  return std::min(1.0, std::max(-1.0, dot_3d/sqrt(a.modp2()*b.modp2())));
}

dot_product()

double dot_product ( const PseudoJet &  a, const PseudoJet &  b  )

inline

                                                                    {
  return a.E()*b.E() - a.px()*b.px() - a.py()*b.py() - a.pz()*b.pz();
}

fastjet_version_string()

FJCORE_BEGIN_NAMESPACE std::string fastjet_version_string

(

)

                                {
  return "FastJet version "+string(fastjet_version)+" [fjcore]";
}

get_pointer()

T* get_pointer ( SharedPtr< T > const &  t )

inline

                                             {
  return t.get();
}

have_same_momentum()

bool have_same_momentum	(	const PseudoJet &	,
		const PseudoJet &
	)

                                                                        {
  return jeta.px() == jetb.px()
    &&   jeta.py() == jetb.py()
    &&   jeta.pz() == jetb.pz()
    &&   jeta.E()  == jetb.E();
}

join() [1/10]

PseudoJet join

(

const PseudoJet &

j1

)

                                                         {
  return join<T>(std::vector<PseudoJet>(1,j1));
}

join() [2/10]

PseudoJet join	(	const PseudoJet &	j1,
		const JetDefinition::Recombiner &	recombiner
	)

                                                               {
  return join<T>(std::vector<PseudoJet>(1,j1), recombiner);
}

join() [3/10]

PseudoJet join	(	const PseudoJet &	j1,
		const PseudoJet &	j2
	)

                                                                               {
  std::vector<PseudoJet> pieces;
  pieces.push_back(j1);
  pieces.push_back(j2);
  return join<T>(pieces);
}

join() [4/10]

PseudoJet join	(	const PseudoJet &	j1,
		const PseudoJet &	j2,
		const JetDefinition::Recombiner &	recombiner
	)

                                                               {
  std::vector<PseudoJet> pieces;
  pieces.reserve(2);
  pieces.push_back(j1);
  pieces.push_back(j2);
  return join<T>(pieces, recombiner);
}

join() [5/10]

PseudoJet join	(	const PseudoJet &	j1,
		const PseudoJet &	j2,
		const PseudoJet &	j3
	)

                                         {
  std::vector<PseudoJet> pieces;
  pieces.push_back(j1);
  pieces.push_back(j2);
  pieces.push_back(j3);
  return join<T>(pieces);
}

join() [6/10]

PseudoJet join	(	const PseudoJet &	j1,
		const PseudoJet &	j2,
		const PseudoJet &	j3,
		const JetDefinition::Recombiner &	recombiner
	)

                                                               {
  std::vector<PseudoJet> pieces;
  pieces.reserve(3);
  pieces.push_back(j1);
  pieces.push_back(j2);
  pieces.push_back(j3);
  return join<T>(pieces, recombiner);
}

join() [7/10]

PseudoJet join	(	const PseudoJet &	j1,
		const PseudoJet &	j2,
		const PseudoJet &	j3,
		const PseudoJet &	j4
	)

                                                               {
  std::vector<PseudoJet> pieces;
  pieces.push_back(j1);
  pieces.push_back(j2);
  pieces.push_back(j3);
  pieces.push_back(j4);
  return join<T>(pieces);
}

join() [8/10]

PseudoJet join	(	const PseudoJet &	j1,
		const PseudoJet &	j2,
		const PseudoJet &	j3,
		const PseudoJet &	j4,
		const JetDefinition::Recombiner &	recombiner
	)

                                                               {
  std::vector<PseudoJet> pieces;
  pieces.reserve(4);
  pieces.push_back(j1);
  pieces.push_back(j2);
  pieces.push_back(j3);
  pieces.push_back(j4);
  return join<T>(pieces, recombiner);
}

join() [9/10]

PseudoJet join

(

const std::vector< PseudoJet > &

pieces

)

                                                                        {
  PseudoJet result(0.0,0.0,0.0,0.0);
  for (unsigned int i=0; i<pieces.size(); i++){
    const PseudoJet it = pieces[i];
    result += it;
  }
  T *cj_struct = new T(pieces);
  result.set_structure_shared_ptr(SharedPtr<PseudoJetStructureBase>(cj_struct));
  return result;
}

join() [10/10]

PseudoJet join	(	const std::vector< PseudoJet > &	pieces,
		const JetDefinition::Recombiner &	recombiner
	)

                                                               {
  PseudoJet result;
  if (pieces.size()>0){
    result = pieces[0];
    for (unsigned int i=1; i<pieces.size(); i++){
      recombiner.plus_equal(result, pieces[i]);
    }
  }
  T *cj_struct = new T(pieces, &recombiner);
  result.set_structure_shared_ptr(SharedPtr<PseudoJetStructureBase>(cj_struct));
  return result;
}

objects_sorted_by_values()

std::vector<T> objects_sorted_by_values	(	const std::vector< T > &	objects,
		const std::vector< double > &	values
	)

                                                            {
  if (objects.size() != values.size()){
    throw Error("fjcore::objects_sorted_by_values(...): the size of the 'objects' vector must match the size of the 'values' vector");
  }
  std::vector<int> indices(values.size());
  for (size_t i = 0; i < indices.size(); i++) {indices[i] = i;}
  sort_indices(indices, values);
  std::vector<T> objects_sorted(objects.size());
  for (size_t i = 0; i < indices.size(); i++) {
    objects_sorted[i] = objects[indices[i]];
  }
  return objects_sorted;
}

operator!()

Selector operator!

(

const Selector &

s

)

                                       {
  return Selector(new SW_Not(s));
}

operator!=() [1/4]

bool operator!= ( const double  val, const PseudoJet &  a  )

inline

{return !(a==val);}

operator!=() [2/4]

bool operator!= ( const PseudoJet &  a, const double  val  )

inline

{return !(a==val);}

operator!=() [3/4]

bool operator!= ( const PseudoJet &  a, const PseudoJet &  b  )

inline

{return !(a==b);}

operator!=() [4/4]

bool operator!= ( SharedPtr< T > const &  t, SharedPtr< U > const &  u  )

inline

                                                                      {
  return t.get() != u.get();
}

operator&&()

Selector operator&&	(	const Selector &	s1,
		const Selector &	s2
	)

                                                              {
  return Selector(new SW_And(s1,s2));
}

operator*() [1/3]

PseudoJet operator*	(	const PseudoJet &	,
		double
	)

                                                          {
  return coeff*jet;
}

operator*() [2/3]

Selector operator*	(	const Selector &	s1,
		const Selector &	s2
	)

                                                             {
  return Selector(new SW_Mult(s1,s2));
}

operator*() [3/3]

PseudoJet operator*	(	double	,
		const PseudoJet &
	)

                                                          {
  jet._ensure_valid_rap_phi();
  PseudoJet coeff_times_jet(jet);
  coeff_times_jet *= coeff;
  return coeff_times_jet;
}

operator+()

PseudoJet operator+	(	const PseudoJet &	,
		const PseudoJet &
	)

                                                                     {
  return PseudoJet(jet1.px()+jet2.px(),
           jet1.py()+jet2.py(),
           jet1.pz()+jet2.pz(),
           jet1.E() +jet2.E()  );
}

operator-()

PseudoJet operator-	(	const PseudoJet &	,
		const PseudoJet &
	)

                                                                     {
  return PseudoJet(jet1.px()-jet2.px(),
           jet1.py()-jet2.py(),
           jet1.pz()-jet2.pz(),
           jet1.E() -jet2.E()  );
}

operator/()

PseudoJet operator/	(	const PseudoJet &	,
		double
	)

                                                          {
  return (1.0/coeff)*jet;
}

operator<()

bool operator< ( SharedPtr< T > const &  t, SharedPtr< U > const &  u  )

inline

                                                                     {
  return t.get() < u.get();
}

operator==() [1/4]

bool operator== ( const double  val, const PseudoJet &  jet  )

inline

{return jet == val;}

operator==() [2/4]

bool operator==	(	const PseudoJet &	,
		const PseudoJet &
	)

                                                          {
  if (a.px() != b.px()) return false;
  if (a.py() != b.py()) return false;
  if (a.pz() != b.pz()) return false;
  if (a.E () != b.E ()) return false;
  if (a.user_index()    != b.user_index()) return false;
  if (a.cluster_hist_index() != b.cluster_hist_index()) return false;
  if (a.user_info_ptr() != b.user_info_ptr()) return false;
  if (a.structure_ptr() != b.structure_ptr()) return false;
  return true;
}

operator==() [3/4]

bool operator==	(	const PseudoJet &	jet,
		const double	val
	)

                                                         {
  if (val != 0)
    throw Error("comparing a PseudoJet with a non-zero constant (double) is not allowed.");
  return (jet.px() == 0 && jet.py() == 0 &&
      jet.pz() == 0 && jet.E() == 0);
}

operator==() [4/4]

bool operator== ( SharedPtr< T > const &  t, SharedPtr< U > const &  u  )

inline

                                                                      {
  return t.get() == u.get();
}

operator||()

Selector operator||	(	const Selector &	s1,
		const Selector &	s2
	)

                                                               {
  return Selector(new SW_Or(s1,s2));
}

PtYPhiM()

PseudoJet PtYPhiM	(	double	pt,
		double	y,
		double	phi,
		double	m = 0.0
	)

                                                             {
  assert(phi < 2*twopi && phi > -twopi);
  double ptm = (m == 0) ? pt : sqrt(pt*pt+m*m);
  double exprap = exp(y);
  double pminus = ptm/exprap;
  double pplus  = ptm*exprap;
  double px = pt*cos(phi);
  double py = pt*sin(phi);
  PseudoJet mom(px,py,0.5*(pplus-pminus),0.5*(pplus+pminus));
  mom.set_cached_rap_phi(y,phi);
  return mom;
}

SelectorAbsEtaMax()

Selector SelectorAbsEtaMax

(

double

absetamax

)

select objects with |eta| <= absetamax

                                             {
  return Selector(new SW_QuantityMax<QuantityAbsEta>(absetamax));
}

SelectorAbsEtaMin()

Selector SelectorAbsEtaMin

(

double

absetamin

)

select objects with |eta| >= absetamin

                                             {
  return Selector(new SW_QuantityMin<QuantityAbsEta>(absetamin));
}

SelectorAbsEtaRange()

Selector SelectorAbsEtaRange	(	double	absetamin,
		double	absetamax
	)

select objects with absetamin <= |eta| <= absetamax

                                                                 {
  return Selector(new SW_QuantityRange<QuantityAbsEta>(absetamin, absetamax));
}

SelectorAbsRapMax()

Selector SelectorAbsRapMax

(

double

absrapmax

)

select objects with |rap| <= absrapmax

                                             {
  return Selector(new SW_AbsRapMax(absrapmax));
}

SelectorAbsRapMin()

Selector SelectorAbsRapMin

(

double

absrapmin

)

select objects with |rap| >= absrapmin

                                             {
  return Selector(new SW_QuantityMin<QuantityAbsRap>(absrapmin));
}

SelectorAbsRapRange()

Selector SelectorAbsRapRange	(	double	absrapmin,
		double	absrapmax
	)

select objects with absrapmin <= |rap| <= absrapmax

                                                           {
  return Selector(new SW_AbsRapRange(rapmin, rapmax));
}

SelectorCircle()

Selector SelectorCircle

(

const double

radius

)

                                             {
  return Selector(new SW_Circle(radius));
}

SelectorDoughnut()

Selector SelectorDoughnut	(	const double	radius_in,
		const double	radius_out
	)

                                                                           {
  return Selector(new SW_Doughnut(radius_in, radius_out));
}

SelectorEMax()

Selector SelectorEMax

(

double

Emax

)

select objects with E <= Emax

                                   {
  return Selector(new SW_QuantityMax<QuantityE>(Emax));
}

SelectorEMin()

Selector SelectorEMin

(

double

Emin

)

select objects with E >= Emin

                                   {
  return Selector(new SW_QuantityMin<QuantityE>(Emin));
}

SelectorERange()

Selector SelectorERange	(	double	Emin,
		double	Emax
	)

select objects with Emin <= E <= Emax

                                                  {
  return Selector(new SW_QuantityRange<QuantityE>(Emin, Emax));
}

SelectorEtaMax()

Selector SelectorEtaMax

(

double

etamax

)

select objects with eta <= etamax

                                       {
  return Selector(new SW_QuantityMax<QuantityEta>(etamax));
}

SelectorEtaMin()

Selector SelectorEtaMin

(

double

etamin

)

select objects with eta >= etamin

                                       {
  return Selector(new SW_QuantityMin<QuantityEta>(etamin));
}

SelectorEtaRange()

Selector SelectorEtaRange	(	double	etamin,
		double	etamax
	)

select objects with etamin <= eta <= etamax

                                                        {
  return Selector(new SW_QuantityRange<QuantityEta>(etamin, etamax));
}

SelectorEtMax()

Selector SelectorEtMax

(

double

Etmax

)

select objects with Et <= Etmax

                                     {
  return Selector(new SW_QuantityMax<QuantityEt2>(Etmax));
}

SelectorEtMin()

Selector SelectorEtMin

(

double

Etmin

)

select objects with Et >= Etmin

                                     {
  return Selector(new SW_QuantityMin<QuantityEt2>(Etmin));
}

SelectorEtRange()

Selector SelectorEtRange	(	double	Etmin,
		double	Etmax
	)

select objects with Etmin <= Et <= Etmax

                                                     {
  return Selector(new SW_QuantityRange<QuantityEt2>(Etmin, Etmax));
}

SelectorIdentity()

Selector SelectorIdentity

(

)

                            {
  return Selector(new SW_Identity);
}

SelectorIsZero()

Selector SelectorIsZero

(

)

                         {
  return Selector(new SW_IsZero());
}

SelectorMassMax()

Selector SelectorMassMax

(

double

Mmax

)

select objects with Mass <= Mmax

                                      {
  return Selector(new SW_QuantityMax<QuantityM2>(mmax));
}

SelectorMassMin()

Selector SelectorMassMin

(

double

Mmin

)

select objects with Mass >= Mmin

                                      {
  return Selector(new SW_QuantityMin<QuantityM2>(mmin));
}

SelectorMassRange()

Selector SelectorMassRange	(	double	Mmin,
		double	Mmax
	)

select objects with Mmin <= Mass <= Mmax

                                                     {
  return Selector(new SW_QuantityRange<QuantityM2>(mmin, mmax));
}

SelectorNHardest()

Selector SelectorNHardest

(

unsigned int

n

)

                                          {
  return Selector(new SW_NHardest(n));
}

SelectorPhiRange()

Selector SelectorPhiRange	(	double	phimin,
		double	phimax
	)

select objects with phimin <= phi <= phimax

                                                        {
  return Selector(new SW_PhiRange(phimin, phimax));
}

SelectorPtFractionMin()

Selector SelectorPtFractionMin

(

double

fraction

)

                                               {
  return Selector(new SW_PtFractionMin(fraction));
}

SelectorPtMax()

Selector SelectorPtMax

(

double

ptmax

)

select objects with pt <= ptmax

                                     {
  return Selector(new SW_QuantityMax<QuantityPt2>(ptmax));
}

SelectorPtMin()

Selector SelectorPtMin

(

double

ptmin

)

select objects with pt >= ptmin

                                     {
  return Selector(new SW_QuantityMin<QuantityPt2>(ptmin));
}

SelectorPtRange()

Selector SelectorPtRange	(	double	ptmin,
		double	ptmax
	)

select objects with ptmin <= pt <= ptmax

                                                     {
  return Selector(new SW_QuantityRange<QuantityPt2>(ptmin, ptmax));
}

SelectorRapMax()

Selector SelectorRapMax

(

double

rapmax

)

select objects with rap <= rapmax

                                       {
  return Selector(new SW_RapMax(rapmax));
}

SelectorRapMin()

Selector SelectorRapMin

(

double

rapmin

)

select objects with rap >= rapmin

                                       {
  return Selector(new SW_RapMin(rapmin));
}

SelectorRapPhiRange()

Selector SelectorRapPhiRange	(	double	rapmin,
		double	rapmax,
		double	phimin,
		double	phimax
	)

                                                                                         {
  return Selector(new SW_RapPhiRange(rapmin, rapmax, phimin, phimax));
}

SelectorRapRange()

Selector SelectorRapRange	(	double	rapmin,
		double	rapmax
	)

select objects with rapmin <= rap <= rapmax

                                                        {
  return Selector(new SW_RapRange(rapmin, rapmax));
}

SelectorRectangle()

Selector SelectorRectangle	(	const double	half_rap_width,
		const double	half_phi_width
	)

                                                                                     {
  return Selector(new SW_Rectangle(half_rap_width, half_phi_width));
}

SelectorStrip()

Selector SelectorStrip

(

const double

half_width

)

                                                {
  return Selector(new SW_Strip(half_width));
}

sort_indices()

void sort_indices	(	std::vector< int > &	indices,
		const std::vector< double > &	values
	)

                                            {
  IndexedSortHelper index_sort_helper(&values);
  sort(indices.begin(), indices.end(), index_sort_helper);
}

sorted_by_E()

std::vector<PseudoJet> sorted_by_E

(

const std::vector< PseudoJet > &

jets

)

                                                              {
  vector<double> energies(jets.size());
  for (size_t i = 0; i < jets.size(); i++) {energies[i] = -jets[i].E();}
  return objects_sorted_by_values(jets, energies);
}

sorted_by_pt()

std::vector<PseudoJet> sorted_by_pt

(

const std::vector< PseudoJet > &

jets

)

                                                               {
  vector<double> minus_kt2(jets.size());
  for (size_t i = 0; i < jets.size(); i++) {minus_kt2[i] = -jets[i].kt2();}
  return objects_sorted_by_values(jets, minus_kt2);
}

sorted_by_pz()

std::vector<PseudoJet> sorted_by_pz

(

const std::vector< PseudoJet > &

jets

)

                                                               {
  vector<double> pz(jets.size());
  for (size_t i = 0; i < jets.size(); i++) {pz[i] = jets[i].pz();}
  return objects_sorted_by_values(jets, pz);
}

sorted_by_rapidity()

std::vector<PseudoJet> sorted_by_rapidity

(

const std::vector< PseudoJet > &

jets

)

                                                                     {
  vector<double> rapidities(jets.size());
  for (size_t i = 0; i < jets.size(); i++) {rapidities[i] = jets[i].rap();}
  return objects_sorted_by_values(jets, rapidities);
}

swap()

void swap ( SharedPtr< T > &  a, SharedPtr< T > &  b  )

inline

                                                    {
  return a.swap(b);
}

theta()

double theta ( const PseudoJet &  a, const PseudoJet &  b  )

inline

                                                              {
  return acos(cos_theta(a,b));
}

Variable Documentation

aachen_algorithm

const JetAlgorithm aachen_algorithm = cambridge_algorithm

cambridge_aachen_algorithm

const JetAlgorithm cambridge_aachen_algorithm = cambridge_algorithm

eulergamma

const double eulergamma = 0.577215664901532860606512090082402431042

ln2

const double ln2 = 0.693147180559945309417232121458176568076

MaxRap

FJCORE_BEGIN_NAMESPACE const double MaxRap = 1e5

pi

FJCORE_BEGIN_NAMESPACE const double pi = 3.141592653589793238462643383279502884197

pisq

const double pisq = 9.869604401089358618834490999876151135314

pseudojet_invalid_phi

const double pseudojet_invalid_phi = -100.0

pseudojet_invalid_rap

const double pseudojet_invalid_rap = -1e200

twopi

const double twopi = 6.283185307179586476925286766559005768394

zeta2

const double zeta2 = 1.644934066848226436472415166646025189219

zeta3

const double zeta3 = 1.202056903159594285399738161511449990765