FastJetProcessor.h

/*
 * FastJetProcessor.h
 *
 *  Created on: 25.05.2010
 *      Author: Lars Weuste (MPP Munich) - weuste@mpp.mpg.de
 *                      iterative inclusive algorithm based on design by Marco Battaglia (CERN) - Marco.Battaglia@cern.ch
 */

#ifndef FASTJETPROCESSOR_H_
#define FASTJETPROCESSOR_H_

#include "EClusterMode.h"

#include "marlin/Processor.h"
#include "marlin/VerbosityLevels.h"
#include "EVENT/LCCollection.h"
#include "EVENT/ReconstructedParticle.h"
#include "IMPL/LCCollectionVec.h"
#include "LCIOSTLTypes.h"

#include <fastjet/PseudoJet.hh>
#include <fastjet/JetDefinition.hh>

#include <vector>
#include <string>

 //Forward Declaration
 class FastJetUtil;
typedef std::vector< fastjet::PseudoJet > PseudoJetList;



class FastJetProcessor : marlin::Processor {
public:
  FastJetProcessor();
  virtual ~FastJetProcessor();


  virtual Processor* newProcessor() {
    return new FastJetProcessor();
  }

  virtual void init();

  virtual void processRunHeader(LCRunHeader* ){};

  virtual void processEvent(LCEvent * evt);

  virtual void check(LCEvent * ) {};

  virtual void end();

  friend class FastJetUtil;

private:

  // the LC Collection names for input/output
  std::string   _lcParticleInName;
  std::string   _lcParticleOutName;
  std::string   _lcJetOutName;

  EVENT::ReconstructedParticle* getRecPar(fastjet::PseudoJet& fj, const PseudoJetList& constituents);

  LCCollection* _reconstructedPars;

  int   _statsFoundJets;
  int _statsNrEvents;
  int _statsNrSkippedEmptyEvents;
  int _statsNrSkippedFixedNrJets;
  int _statsNrSkippedMaxIterations;
  bool _storeParticlesInJets;

  FastJetUtil* _fju;

private:
  FastJetProcessor(const FastJetProcessor& rhs);
  FastJetProcessor& operator=(const FastJetProcessor&) { return *this; }

};

std::ostream& operator<<(std::ostream& out, EClusterMode& m);

#endif /* FASTJETPROCESSOR_H_ */