test.cc

#include <fastjet/ClusterSequenceArea.hh>
#include "fastjet/tools/BackgroundEstimator.hh"
#include <iostream>

using namespace fastjet;
using namespace std;

class Test{
public:
  Test(const BackgroundEstimator & bge) : _bge(bge){}
  double get_rho(){ return _bge.rho();}
private:
  BackgroundEstimator _bge;
};

int main(){
  vector<PseudoJet> input;
  input.push_back(PtYPhiM(10.0,0.0,1.0,0.0));
  ClusterSequenceArea csa(input, JetDefinition(kt_algorithm, 0.5),
                          AreaDefinition(active_area_explicit_ghosts, GhostedAreaSpec(5.0)));
  Selector sel = SelectorAbsRapMax(4.0);

  BackgroundEstimator bge(csa, sel);

  Test test(bge);
  cout << test.get_rho() << endl;

  return 0;
}