VHbbDataFormats/src/QjetsPlugin.C

#include "QjetsPlugin.h"

QjetsPlugin::QjetsPlugin(double zcut, double dcut_fctr, double exp_min, double exp_max, double rigidity)
  : _zcut(zcut), _dcut_fctr(dcut_fctr), _exp_min(exp_min), _exp_max(exp_max), _rigidity(rigidity), _rand_seed_set(false)
{
}

void QjetsPlugin::SetRandSeed(unsigned int seed){
  _rand_seed_set = true;
  _seed = seed;
}

double QjetsPlugin::R()const{
  return 0.;
}

string QjetsPlugin::description() const{
  string desc("Qjets pruning plugin");
  return desc;
}

void QjetsPlugin::run_clustering(fastjet::ClusterSequence & cs) const{
  Qjets qjets(_zcut, _dcut_fctr, _exp_min, _exp_max, _rigidity);
  if(_rand_seed_set)
    qjets.SetRandSeed(_seed);
  
  qjets.Cluster(cs);
}
Revision:	1.1.2.1
Committed:	Tue Nov 27 05:04:18 2012 UTC (12 years, 5 months ago) by ntran
Content type:	text/plain
Branch:	hbbsubstructDevPostHCP
CVS Tags:	hbbsubstructDev_6, hbbsubstructDev_5
Changes since 1.1:	+28 -0 lines
Log Message:	major update, including new proposal for FJ3 inputs to step 2
#	User	Rev	Content
1	ntran	1.1.2.1	#include "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			}