FastJetFilter.h

#ifndef __FastJetFilter_h__
#define __FastJetFilter_h__

#include <vector>

#include "StarGenerator/FILT/StarFilterMaker.h"

#ifndef __CINT__
#include <fastjet/PseudoJet.hh>
#include <fastjet/ClusterSequence.hh>
#include <fastjet/JetDefinition.hh>
#include <fastjet/Selector.hh>
#endif

class FastJetFilter : public StarFilterMaker
{
public:
  FastJetFilter( const char* name = "partfilt" );
  virtual ~FastJetFilter(){ /* nada */ }

  int Init();

  int Filter( StarGenEvent *event = 0 );

  void AddTrigger( int _pdgid, double _ptmn=0, double _ptmx=-1, double _etamn=-1, double _etamx=1, int _pdgidParent=0 );

private:
protected:

  struct Trigger_t {
    int    pdgid;
    double ptmn;
    double ptmx;
    double etamn;
    double etamx;
    int    pdgid_parent;
    Trigger_t() :
      pdgid(0), 
      ptmn (0.), 
      ptmx (0.), 
      etamn(0.), 
      etamx(0.),
      pdgid_parent(0)
    { };
    Trigger_t( int _pdgid, double _ptmn, double _ptmx, double _etamn, double _etamx, int _pdgid_parent ) : 
      pdgid(_pdgid), 
      ptmn (_ptmn), 
      ptmx (_ptmx), 
      etamn(_etamn), 
      etamx(_etamx),
      pdgid_parent(_pdgid_parent)
    {
    };
  };
  vector<Trigger_t> mTriggers;

#ifndef __CINT__
  fastjet::Strategy            strategy = fastjet::Best;
  fastjet::RecombinationScheme recombScheme = fastjet::E_scheme; //Change as you need
  fastjet::JetDefinition*      jetdefinition = nullptr;
  fastjet::JetAlgorithm        algorithm;
#endif

#if 1 // we dont really need, but this triggers cons to create dictionary
  ClassDef(FastJetFilter,0);
#endif

};

#endif