StConePars.h

// -*- mode: c++;-*-
// $Id: StConePars.h,v 1.8 2008/07/06 01:13:43 tai Exp $
#ifndef STCONEPARS_H
#define STCONEPARS_H

#include "StJetPars.h"

#include <cmath>

class StConePars : public StJetPars {

public:

  StConePars()
    : mNeta(105)
    , mNphi(120)
    , mEtaMin(-3.0)
    , mEtaMax(3.0)
    , mPhiMin(-atan(1.0)*4.0)
    , mPhiMax(atan(1.0)*4.0)
    , mR(0.7)
    , mSeedEtMin(0.5)
    , mAssocEtMin(0.1)
    , mSplitFraction(0.5)
    , mDoMinimization(true)
    , mAddMidpoints(true)
    , mRequireStableMidpoints(true)
    , mDoSplitMerge(true)
    , mDebug(false)
  { }
    

  virtual StJetFinder* constructJetFinder();

  void setGridSpacing(int nEta, double etaMin, double etaMax,
                      int nPhi, double phiMin, double phiMax)
  {
    mNeta = nEta; mEtaMin = etaMin; mEtaMax = etaMax;
    mNphi = nPhi; mPhiMin = phiMin; mPhiMax = phiMax;
  }
    
  void setSeedEtMin(double v) { mSeedEtMin = v; }

  void setAssocEtMin(double v) { mAssocEtMin = v; }

  void setSplitFraction(double v) { mSplitFraction = v; }

  void setPerformMinimization(bool v) {mDoMinimization=v;}

  void setAddMidpoints(bool v) {mAddMidpoints=v;}

  void setDoSplitMerge(bool v) {mDoSplitMerge=v;}

  void setRequireStableMidpoints(bool v) {mRequireStableMidpoints=v;}

  void setConeRadius(double v) {mR = v;}

  void setDebug(bool v) {mDebug = v;}

  int    Neta()   const { return mNeta; }
  int    Nphi()   const { return mNphi; } 
  double EtaMin() const { return mEtaMin; }
  double EtaMax() const { return mEtaMax; }
  double PhiMin() const { return mPhiMin; }
  double PhiMax() const { return mPhiMax; }

  double coneRadius() const {return mR;}

  double seedEtMin() const { return mSeedEtMin; }
  double assocEtMin() const { return mAssocEtMin; }

  double splitFraction() const { return mSplitFraction; }

  bool performMinimization() const {return mDoMinimization;}

  bool addMidpoints() const {return mAddMidpoints;}

  bool requiredStableMidpoints() const {return mRequireStableMidpoints;}

  bool doSplitMerge() const {return mDoSplitMerge;}

  bool debug() const {return mDebug;}
    
  double phiWidth() const { return (PhiMax() - PhiMin())/static_cast<double>(Nphi()); }
  double etaWidth() const { return (EtaMax() - EtaMin())/static_cast<double>(Neta()); }
  int    deltaPhi() const { return static_cast<int>(floor(coneRadius()/phiWidth())) + 1; }
  int    deltaEta() const { return static_cast<int>(floor(coneRadius()/etaWidth())) + 1; }


private:

  int mNeta;
  int mNphi;
  double mEtaMin;
  double mEtaMax;
  double mPhiMin;
  double mPhiMax;
  double mR;
  double mSeedEtMin;
  double mAssocEtMin;
  double mSplitFraction;
  bool mDoMinimization;
  bool mAddMidpoints;
  bool mRequireStableMidpoints;
  bool mDoSplitMerge;
  bool mDebug;

  ClassDef(StConePars,1)
};

#endif // STCONEPARS_H