30 #ifndef __FASTJET_COMPOSITEJET_STRUCTURE_HH__ 31 #define __FASTJET_COMPOSITEJET_STRUCTURE_HH__ 33 #include <fastjet/PseudoJet.hh> 34 #include <fastjet/PseudoJetStructureBase.hh> 37 #include <fastjet/JetDefinition.hh> 39 FASTJET_BEGIN_NAMESPACE
60 if (_area_4vector_ptr)
delete _area_4vector_ptr;
64 virtual std::string description()
const;
70 virtual bool has_constituents()
const;
78 virtual std::vector<PseudoJet> constituents(
const PseudoJet &jet)
const;
87 virtual std::vector<PseudoJet> pieces(
const PseudoJet &jet)
const;
93 virtual bool has_area()
const;
96 virtual double area(
const PseudoJet &reference)
const;
102 virtual double area_error(
const PseudoJet &reference)
const;
110 virtual bool is_pure_ghost(
const PseudoJet &reference)
const;
115 void set_area_information(
PseudoJet *area_4vector_ptr){
116 _area_4vector_ptr = area_4vector_ptr;
137 template<
typename T>
PseudoJet join(
const std::vector<PseudoJet> & pieces){
139 for (
unsigned int i=0; i<pieces.size(); i++){
144 T *cj_struct =
new T(pieces);
154 return join<T>(std::vector<PseudoJet>(1,j1));
160 std::vector<PseudoJet> pieces;
161 pieces.push_back(j1);
162 pieces.push_back(j2);
163 return join<T>(pieces);
170 std::vector<PseudoJet> pieces;
171 pieces.push_back(j1);
172 pieces.push_back(j2);
173 pieces.push_back(j3);
174 return join<T>(pieces);
181 std::vector<PseudoJet> pieces;
182 pieces.push_back(j1);
183 pieces.push_back(j2);
184 pieces.push_back(j3);
185 pieces.push_back(j4);
186 return join<T>(pieces);
200 template<
typename T>
PseudoJet join(
const std::vector<PseudoJet> & pieces,
203 if (pieces.size()>0){
205 for (
unsigned int i=1; i<pieces.size(); i++){
210 T *cj_struct =
new T(pieces, &recombiner);
220 return join<T>(std::vector<PseudoJet>(1,j1), recombiner);
227 std::vector<PseudoJet> pieces;
228 pieces.push_back(j1);
229 pieces.push_back(j2);
230 return join<T>(pieces, recombiner);
238 std::vector<PseudoJet> pieces;
239 pieces.push_back(j1);
240 pieces.push_back(j2);
241 pieces.push_back(j3);
242 return join<T>(pieces, recombiner);
250 std::vector<PseudoJet> pieces;
251 pieces.push_back(j1);
252 pieces.push_back(j2);
253 pieces.push_back(j3);
254 pieces.push_back(j4);
255 return join<T>(pieces, recombiner);
259 FASTJET_END_NAMESPACE
261 #endif // __FASTJET_MERGEDJET_STRUCTURE_HH__ virtual bool has_pieces(const PseudoJet &) const
true if it has pieces (always the case)
virtual ~CompositeJetStructure()
default dtor
The structure for a jet made of pieces.
Contains any information related to the clustering that should be directly accessible to PseudoJet...
std::vector< PseudoJet > _pieces
the pieces building the jet
PseudoJet * _area_4vector_ptr
pointer to the 4-vector jet area
an implementation of C++0x shared pointers (or boost's)
void set_structure_shared_ptr(const SharedPtr< PseudoJetStructureBase > &structure)
set the associated structure
CompositeJetStructure()
default ctor
An abstract base class that will provide the recombination scheme facilities and/or allow a user to e...
Class to contain pseudojets, including minimal information of use to jet-clustering routines...
void plus_equal(PseudoJet &pa, const PseudoJet &pb) const
pa += pb in the given recombination scheme.