// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/RivetAIDA.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FastJets.hh"

namespace Rivet {

  class JETS_TEMPLATE : public Analysis {
  public:

    JETS_TEMPLATE()
      : Analysis("JETS_TEMPLATE")
    {
      setBeams(PROTON, PROTON);
      setNeedsCrossSection(false);
    }

  public:

    void init() {

      const FinalState fsj (-5.2, 5.2, 0.0 * GeV);
        addProjection (FastJets (fsj, FastJets::ANTIKT, 0.5), "jetlabel");

    }

    void analyze(const Event& event) {

      const Jets& jets = applyProjection < FastJets > (event, "jetlabel").jetsByPt ();
	getLog() << Log::DEBUG 
		 << "[jetSize, LeadJetpT] = ["
		 << jets.size() << ", "
		 << jets[0].momentum().pT() << "]"
		 << endl;
    }

    void finalize() {
     
    }

  private:


  };

  AnalysisBuilder<JETS_TEMPLATE> plugin_JETS_TEMPLATE;

}