StFastJetPars.h

// -*- mode: c++ -*-
//
// Pibero Djawotho <pibero@tamu.edu>
// Texas A&M University
// 31 Aug 2011
//
// $Id: StFastJetPars.h,v 1.4 2016/01/06 22:00:17 gdwebb Exp $
//
// $Log: StFastJetPars.h,v $
// Revision 1.4  2016/01/06 22:00:17  gdwebb
// This is code to implement the off axis cone underlying event analysis.
//
// Revision 1.3  2012/03/10 23:18:28  pibero
// Added destructor for StCDFMidPointPlugin
//
// Revision 1.2  2012/03/10 23:09:53  pibero
// Addeed support for fastjet plugins
//
// Revision 1.1  2011/08/31 17:57:36  pibero
// Support for FastJet
//
// http://www.lpthe.jussieu.fr/~salam/fastjet/
// http://www.star.bnl.gov/HyperNews-star/protected/get/starsoft/8521.html
//
//

#ifndef ST_FAST_JET_PARS_H
#define ST_FAST_JET_PARS_H

#include "StJetPars.h"
#include "StFastJetAreaPars.h"

class StPlugin : public TObject {
public:
  StPlugin(void* i) : mImpl(i) {}
  virtual ~StPlugin() {}
  void* impl() const { return mImpl; }
protected:
  void* mImpl;
  ClassDef(StPlugin,0)
};

class StCDFMidPointPlugin : public StPlugin {
public:
  StCDFMidPointPlugin(double coneRadius,
                      double overlapTreshold,
                      double seedThreshold,
                      double coneAreaFraction);
  ~StCDFMidPointPlugin();
  ClassDef(StCDFMidPointPlugin,0)
};

class StFastJetPars : public StJetPars {
public:
  static const int kt_algorithm;
  static const int cambridge_algorithm;
  static const int antikt_algorithm;
  static const int genkt_algorithm;
  static const int cambridge_for_passive_algorithm;
  static const int genkt_for_passive_algorithm;
  //.................................................................
  static const int ee_kt_algorithm;
  static const int ee_genkt_algorithm;
  //.................................................................
  static const int plugin_algorithm;

  static const int E_scheme;
  static const int pt_scheme;
  static const int pt2_scheme;
  static const int Et_scheme;
  static const int Et2_scheme;
  static const int BIpt_scheme;
  static const int BIpt2_scheme;
  static const int external_scheme;

  static const int N2MinHeapTiled;
  static const int N2Tiled;
  static const int N2PoorTiled;
  static const int N2Plain;
  static const int N3Dumb;
  static const int Best;
  static const int NlnN;
  static const int NlnN3pi;
  static const int NlnN4pi;
  static const int NlnNCam4pi;
  static const int NlnNCam2pi2R;
  static const int NlnNCam;
  static const int plugin_strategy;

  StFastJetPars()
    : mJetAlgorithm(kt_algorithm)
    , mRparam(1.0)
    , mRecombScheme(E_scheme)
    , mStrategy(Best)
    , mPtMin(5.0)
    , mPlugin(0)
    , mJetAreaFlag(false)
  {
  }

  StJetFinder* constructJetFinder();

  int    jetAlgorithm       () const { return mJetAlgorithm; }
  double Rparam             () const { return mRparam;       }
  int    recombinationScheme() const { return mRecombScheme; }
  int    strategy           () const { return mStrategy;     }
  double ptMin              () const { return mPtMin;        }
  StFastJetAreaPars* jetArea () const { return mJetArea;      }
  bool jetAreaFlag          () const { return mJetAreaFlag;  }

  void* plugin() const { return mPlugin->impl(); }

  void setJetAlgorithm       (int    jetAlgorithm) { mJetAlgorithm = jetAlgorithm; }
  void setRparam             (double Rparam      ) { mRparam       = Rparam;       }
  void setRecombinationScheme(int    recombScheme) { mRecombScheme = recombScheme; }
  void setStrategy           (int    strategy    ) { mStrategy     = strategy;     }
  void setPtMin              (double ptmin       ) { mPtMin        = ptmin;        }
  void setPlugin(StPlugin* plugin) { mPlugin = plugin; }
  void setJetArea (StFastJetAreaPars* area_pars){ mJetArea = area_pars; mJetAreaFlag = true;};
private:
  int    mJetAlgorithm;
  double mRparam;
  int    mRecombScheme;
  int    mStrategy;
  double mPtMin;
  StPlugin* mPlugin;
  StFastJetAreaPars *mJetArea;
  bool mJetAreaFlag;

  ClassDef(StFastJetPars,0)
};

#endif  // ST_FAST_JET_PARS_H