JetFinder.h

Go to the documentation of this file.

#ifndef JetFinder_h
#define JetFinder_h 1

#include "lcfiplus.h"
#include "TObject.h"

using namespace std;
using namespace lcfiplus;

namespace lcfiplus {

struct JetConfig {
  string algo;
  //string algoY;
  int nJet;
  double Ycut;
  int useBeamJets;
  double rParameter;     // for Kt
  double alphaParameter; // for Durham
  double betaParameter;  // for Valencia
  double coneR;
  double epsCut;
  string coreAlgo;
  double coreThreshold;
  double distCut;
  int nIteration;
  double YaddVV;
  double YaddVL;
  double YaddLL;
  int useMuonID;
  int muonIDExternal;
  double muonIDMinEnergy;
  double muonIDMinD0Sig;
  double muonIDMinZ0Sig;
  double muonIDMaxDist;
  double muonIDMinProb;

  JetConfig() : 
    algo("DurhamVertex"),
      nJet(0),
      Ycut(1),
      useBeamJets(0),
      rParameter(1.0),
      alphaParameter(1.0),
      betaParameter(1.0),
      coneR(0),
      epsCut(0), 
      coreThreshold(0),
      distCut(0),
      nIteration(5),
      YaddVV(100), 
      YaddVL(100), 
      YaddLL(100),
      useMuonID(1),
      muonIDExternal(1),
      muonIDMinEnergy(5.),
      muonIDMinD0Sig(5.),
      muonIDMinZ0Sig(5.),
      muonIDMaxDist(5.),
      muonIDMinProb(.5)
      {}
};

Jet* convertJetVertex(const Jet* jet);

class JetFinder {
 public:
  static double funcDurham(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);
  static double funcJade(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);
  static double funcJadeE(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);
  static double funcDurhamCheat(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);

  static double funcKt(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);
  static double funcValencia(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);

  static double funcVertex(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg, double(*func)(Jet&, Jet&,double,JetConfig&));

  // Vertex functions using above funcVertex()
  // using template functions are more smart, but we maintain this for readability of non-experts
  static double funcDurhamVertex(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);
  static double funcKtVertex(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);
  static double funcValenciaVertex(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);

  static double funcDurhamBeamDistance(Jet& jet, double Evis2, JetConfig& cfg);
  static double funcKtBeamDistance(Jet& jet, double Evis2, JetConfig& cfg);
  static double funcValenciaBeamDistance(Jet& jet, double Evis2, JetConfig& cfg);

  JetFinder(const JetConfig& cfg);
  ~JetFinder() {};
  void Configure(const JetConfig& cfg);
  vector<Jet*> run(TrackVec& tracks);
  vector<Jet*> run(TrackVec& tracks, NeutralVec& neutrals, double* pymin = 0, int ynjetmax = 10);
  vector<Jet*> run(TrackVec& tracks, NeutralVec& neutrals, VertexVec& vertices, double* pymin = 0, int ynjetmax = 10);
  vector<Jet*> prerun(TrackVec& tracks, NeutralVec& neutrals, VertexVec& vertices, int* pnVertexJets = 0);
  vector<Jet*> run(vector<Jet*> input, double* pymin = 0, int ynjetmax = 10);

 private:
  double (*_Yfunc)(Jet& jet1, Jet& jet2, double Evis2, JetConfig& cfg);
  double (*_YfuncBeam)(Jet& jet1, double Evis2, JetConfig& cfg);
  JetConfig _cfg;
};

/*
class CheatedJetFinder {
  public:
    CheatedJetFinder(const JetConfig& cfg);
    ~CheatedJetFinder() {};
    vector<Jet*> run(Event* event);
    vector<MCParticle*> originalPartons( MCParticleVec & mcps );

  private:
    double (*_Yfunc)(Jet& jet1, Jet& jet2, double Evis2);
    JetConfig _cfg;
};
*/
}

#endif