| 1 |
#include "UHHAnalysis/NtupleWriter/interface/QjetsPlugin.h"
|
| 2 |
|
| 3 |
QjetsPlugin::QjetsPlugin(double zcut, double dcut_fctr, double exp_min, double exp_max, double rigidity)
|
| 4 |
: _zcut(zcut), _dcut_fctr(dcut_fctr), _exp_min(exp_min), _exp_max(exp_max), _rigidity(rigidity), _rand_seed_set(false)
|
| 5 |
{
|
| 6 |
}
|
| 7 |
|
| 8 |
void QjetsPlugin::SetRandSeed(unsigned int seed){
|
| 9 |
_rand_seed_set = true;
|
| 10 |
_seed = seed;
|
| 11 |
}
|
| 12 |
|
| 13 |
double QjetsPlugin::R()const{
|
| 14 |
return 0.;
|
| 15 |
}
|
| 16 |
|
| 17 |
string QjetsPlugin::description() const{
|
| 18 |
string desc("Qjets pruning plugin");
|
| 19 |
return desc;
|
| 20 |
}
|
| 21 |
|
| 22 |
void QjetsPlugin::run_clustering(fastjet::ClusterSequence & cs) const{
|
| 23 |
Qjets qjets(_zcut, _dcut_fctr, _exp_min, _exp_max, _rigidity);
|
| 24 |
if(_rand_seed_set)
|
| 25 |
qjets.SetRandSeed(_seed);
|
| 26 |
|
| 27 |
qjets.Cluster(cs);
|
| 28 |
}
|
