| 1 | #ifndef _JETCLU_ALGORITHM_HH_
|
|---|
| 2 | #define _JETCLU_ALGORITHM_HH_
|
|---|
| 3 |
|
|---|
| 4 | //----------------------------------------------------------------------
|
|---|
| 5 | // This file distributed with FastJet has been obtained from
|
|---|
| 6 | // http://www.pa.msu.edu/~huston/Les_Houches_2005/JetClu+Midpoint-StandAlone.tgz
|
|---|
| 7 | //
|
|---|
| 8 | // Permission to distribute it with FastJet has been granted by Joey
|
|---|
| 9 | // Huston (see the COPYING file in the main FastJet directory for
|
|---|
| 10 | // details).
|
|---|
| 11 | // Changes from the original file are listed below.
|
|---|
| 12 | //----------------------------------------------------------------------
|
|---|
| 13 |
|
|---|
| 14 | // History of changes compared to the original JetCluAlgorithm.hh file
|
|---|
| 15 | //
|
|---|
| 16 | // 2009-01-17 Gregory Soyez <soyez@fastjet.fr>
|
|---|
| 17 | //
|
|---|
| 18 | // * put the code in the fastjet::cdf namespace
|
|---|
| 19 | //
|
|---|
| 20 | // 2006-09-24 Gavin Salam <salam@lpthe.jussieu.fr>
|
|---|
| 21 | //
|
|---|
| 22 | // * added JetClu+MidPoint to FastJet
|
|---|
| 23 |
|
|---|
| 24 | #include "PhysicsTower.hh"
|
|---|
| 25 | #include "Cluster.hh"
|
|---|
| 26 | #include <vector>
|
|---|
| 27 |
|
|---|
| 28 | #include <fastjet/internal/base.hh>
|
|---|
| 29 |
|
|---|
| 30 | FASTJET_BEGIN_NAMESPACE
|
|---|
| 31 |
|
|---|
| 32 | namespace cdf{
|
|---|
| 33 |
|
|---|
| 34 | class JetCluAlgorithm
|
|---|
| 35 | {
|
|---|
| 36 | private:
|
|---|
| 37 | double _seedThreshold;
|
|---|
| 38 | double _coneRadius;
|
|---|
| 39 | int _adjacencyCut;
|
|---|
| 40 | int _maxIterations;
|
|---|
| 41 | int _iratch;
|
|---|
| 42 | double _overlapThreshold;
|
|---|
| 43 |
|
|---|
| 44 | public:
|
|---|
| 45 | JetCluAlgorithm():
|
|---|
| 46 | _seedThreshold(1),
|
|---|
| 47 | _coneRadius(0.7),
|
|---|
| 48 | _adjacencyCut(2),
|
|---|
| 49 | _maxIterations(100),
|
|---|
| 50 | _iratch(1),
|
|---|
| 51 | _overlapThreshold(0.75)
|
|---|
| 52 | {}
|
|---|
| 53 | JetCluAlgorithm(double st, double cr, int ac, int mi, int ir, double ot):
|
|---|
| 54 | _seedThreshold(st),
|
|---|
| 55 | _coneRadius(cr),
|
|---|
| 56 | _adjacencyCut(ac),
|
|---|
| 57 | _maxIterations(mi),
|
|---|
| 58 | _iratch(ir),
|
|---|
| 59 | _overlapThreshold(ot)
|
|---|
| 60 | {}
|
|---|
| 61 | void makeSeedTowers(std::vector<PhysicsTower>& towers, std::vector<Cluster>& seedTowers);
|
|---|
| 62 | void buildPreClusters(std::vector<Cluster>& seedTowers, std::vector<PhysicsTower>& towers, std::vector<Cluster>& preClusters);
|
|---|
| 63 | void findStableCones(std::vector<Cluster>& preClusters, std::vector<PhysicsTower>& towers, std::vector<Cluster>& stableCones);
|
|---|
| 64 | void splitAndMerge(std::vector<Cluster>& stableCones, std::vector<Cluster>& jets);
|
|---|
| 65 | void run(std::vector<PhysicsTower>& towers, std::vector<Cluster>& jets);
|
|---|
| 66 | };
|
|---|
| 67 |
|
|---|
| 68 | } // namespace cdf
|
|---|
| 69 |
|
|---|
| 70 | FASTJET_END_NAMESPACE
|
|---|
| 71 |
|
|---|
| 72 | #endif
|
|---|
