Rivet/Tools/RivetFastJet.hh

0001 #ifndef RIVET_RIVETFASTJET_HH
0002 #define RIVET_RIVETFASTJET_HH
0003
0004 #include "Rivet/Config/RivetCommon.hh"
0005 #include "Rivet/Math/Vectors.hh"
0006
0007 #include "fastjet/JetDefinition.hh"
0008 #include "fastjet/AreaDefinition.hh"
0009 #include "fastjet/ClusterSequence.hh"
0010 #include "fastjet/ClusterSequenceArea.hh"
0011 #include "fastjet/PseudoJet.hh"
0012 #include "fastjet/tools/Filter.hh"
0013 #include "fastjet/tools/Recluster.hh"
0014
0015 namespace fastjet {
0016   namespace contrib { }
0017 }
0018
0019 namespace Rivet {
0020
0021   namespace fjcontrib = fastjet::contrib;
0022
0023   /// Unscoped awareness of FastJet's PseudoJet
0024   using fastjet::PseudoJet;
0025   using fastjet::ClusterSequence;
0026   using fastjet::JetDefinition;
0027
0028   /// Typedef for a collection of PseudoJet objects.
0029   /// @todo Make into an explicit container with conversion to Jets and FourMomenta?
0030   typedef std::vector<PseudoJet> PseudoJets;
0031
0032
0033   /// Make a 3-momentum vector from a FastJet pseudojet
0034   inline Vector3 momentum3(const fastjet::PseudoJet& pj) {
0035     return Vector3(pj.px(), pj.py(), pj.pz());
0036   }
0037
0038   /// Make a 4-momentum vector from a FastJet pseudojet
0039   inline FourMomentum momentum(const fastjet::PseudoJet& pj) {
0040     return FourMomentum(pj.E(), pj.px(), pj.py(), pj.pz());
0041   }
0042
0043
0044 }
0045
0046 #endif