29 #include "fastjet/D0RunIIConePlugin.hh" 30 #include "fastjet/ClusterSequence.hh" 31 #include "fastjet/Error.hh" 36 #include "ILConeAlgorithm.hpp" 37 #include "HepEntity.h" 39 FASTJET_BEGIN_NAMESPACE
44 const double D0RunIIConePlugin::_DEFAULT_split_ratio = 0.5 ;
45 const double D0RunIIConePlugin::_DEFAULT_far_ratio = 0.5 ;
46 const double D0RunIIConePlugin::_DEFAULT_Et_min_ratio = 0.5 ;
47 const bool D0RunIIConePlugin::_DEFAULT_kill_duplicate = true ;
48 const double D0RunIIConePlugin::_DEFAULT_duplicate_dR = 0.005;
49 const double D0RunIIConePlugin::_DEFAULT_duplicate_dPT = 0.01 ;
50 const double D0RunIIConePlugin::_DEFAULT_search_factor = 1.0 ;
51 const double D0RunIIConePlugin::_DEFAULT_pT_min_leading_protojet = 0. ;
52 const double D0RunIIConePlugin::_DEFAULT_pT_min_second_protojet = 0. ;
53 const int D0RunIIConePlugin::_DEFAULT_merge_max = 10000;
54 const double D0RunIIConePlugin::_DEFAULT_pT_min_nomerge = 0. ;
56 bool D0RunIIConePlugin::_first_time =
true;
58 string D0RunIIConePlugin::description ()
const {
61 desc <<
"D0 Run II Improved Legacy (midpoint) cone jet algorithm, with ";
62 desc <<
"cone_radius = " << cone_radius () <<
", " 63 <<
"min_jet_Et = " << min_jet_Et () <<
", " 64 <<
"split_ratio = " << split_ratio ();
75 vector<HepEntity> entities(clust_seq.
jets().size());
76 list<const HepEntity * > ensemble;
77 for (
unsigned i = 0; i < clust_seq.
jets().size(); i++) {
78 entities[i].Fill(clust_seq.
jets()[i].E(),
79 clust_seq.
jets()[i].px(),
80 clust_seq.
jets()[i].py(),
81 clust_seq.
jets()[i].pz(),
84 if (abs(entities[i].pz) < entities[i].E) {
85 ensemble.push_back(& (entities[i]));
90 ILConeAlgorithm<HepEntity>
100 pT_min_leading_protojet(),
101 pT_min_second_protojet(),
106 float Item_ET_Threshold = 0.;
107 list<HepEntity> jets;
108 ilegac.makeClusters(jets, ensemble, Item_ET_Threshold);
112 for(
int i = ilegac.ilcv.size()-1; i >= 0; i--) {
114 std::list<const HepEntity*> tlist = ilegac.ilcv[i].LItems();
115 std::list<const HepEntity*>::iterator tk;
126 int jet_k = (*tk)->index;
129 for (; tk != tlist.end(); tk++) {
131 int jet_j = (*tk)->index;
138 double d_iB = clust_seq.
jets()[jet_k].perp2();
145 void D0RunIIConePlugin::_print_banner(ostream *ostr)
const{
146 if (! _first_time)
return;
152 (*ostr) <<
"#--------------------------------------------------------------------------" << endl;
153 (*ostr) <<
"# You are running the D0 Run II Cone plugin for FastJet " << endl;
154 (*ostr) <<
"# Original code by the D0 collaboration, provided by Lars Sonnenschein; " << endl;
155 (*ostr) <<
"# interface by FastJet authors " << endl;
156 (*ostr) <<
"# If you use this plugin, please cite " << endl;
157 (*ostr) <<
"# G. C. Blazey et al., hep-ex/0005012 " << endl;
158 (*ostr) <<
"# V. M. Abazov et al. [D0 Collaboration], arXiv:1110.3771 [hep-ex] " << endl;
159 (*ostr) <<
"# in addition to the usual FastJet reference. " << endl;
160 (*ostr) <<
"#--------------------------------------------------------------------------" << endl;
166 FASTJET_END_NAMESPACE
void plugin_record_ij_recombination(int jet_i, int jet_j, double dij, int &newjet_k)
record the fact that there has been a recombination between jets()[jet_i] and jets()[jet_k], with the specified dij, and return the index (newjet_k) allocated to the new jet, whose momentum is assumed to be the 4-vector sum of that of jet_i and jet_j
void plugin_record_iB_recombination(int jet_i, double diB)
record the fact that there has been a recombination between jets()[jet_i] and the beam...
static std::ostream * fastjet_banner_stream()
returns a pointer to the stream to be used to print banners (cout by default).
const std::vector< PseudoJet > & jets() const
allow the user to access the internally stored _jets() array, which contains both the initial particl...