JetSplitMergeTool.hh

#ifndef JETREC_JETSPLITMERGETOOL_H
#define JETREC_JETSPLITMERGETOOL_H

//----------------------------------------------------------------------
// This file distributed with FastJet has been obtained from SpartyJet
// v2.20.0 by Pierre-Antoine Delsart, Kurtis L. Geerlings, Joey
// Huston, Brian T. Martin and Chris Vermilion
// For details, see http://www.pa.msu.edu/~huston/SpartyJet/
//                  http://projects.hepforge.org/spartyjet/
//
// Changes from the original file are listed below.
//----------------------------------------------------------------------

//******************************************************************************
// Filename :  JetSplitMergeTool   
// Author   :  Ambreesh Gupta
// Created  :  November, 2001
//
// DESCRIPTION:
//
// Cone algorithms generally produce overlapping jets. These jets are then
// passed through a Split Merge algorithm to create seperated jets. This
// is an implementation of one such algorithm.
// 
// The steps of the algorithm are:
//          - <to be written>
//
// PROPERTIES (JobOption Parameters):
//
//    ProtoJetContainerLoc  string "Default"   Key for ProtoJet input.
//    JetContainerLoc       string "Default"   Key for Jet list to output.
//    OverlapFraction       double 0.5
//
// HISTORY
//   05Nov01 agupta.  First version
//   15Aug02 agupta.  Fix phi-wrapping. 
//*****************************************************************************


// History of changes from the original JetSplitMergeTool.hh file in
// SpartyJet v2.20
// 
// 2009-01-15  Gregory Soyez  <soyez@fastjet.fr>
// 
//        * put the code in the fastjet::atlas namespace
// 
// 2009-02-14  Gregory Soyez  <soyez@fastjet.fr>
// 
//        * imported into FastJet
//        * removed the message logs
//        * replaced StatusCode by int


//class StoreGateSvc;

//Includes
#include "Jet.hh"

#include <vector>
#include <list>

//#include "JetCore/CustomMessage.hh"

#include <fastjet/internal/base.hh>

FASTJET_BEGIN_NAMESPACE

namespace atlas { 

//typedef int StatusCode ;

class JetSplitMergeTool 
{
public:
  typedef Jet::jet_list_t jetcollection_t;
    //Constructors
    JetSplitMergeTool();

    //Destructor
    virtual ~JetSplitMergeTool();

    virtual int execute( jetcollection_t* theJets );

    //Handle negative energy ProtoJets:
    double etaTrue(Jet::constit_vect_t::iterator pj);
    double phiTrue(Jet::constit_vect_t::iterator pj);

    //private:
    void                                     split_merge();   
    double                                             m_f;
    jetcollection_t                               m_preJet;
    jetcollection_t                                  m_jet;
    jetcollection_t*                              m_jetVec;

    int m_ctr;
    int m_dctr;

  //    Message m_log;
};

}  // namespace atlas

FASTJET_END_NAMESPACE

#endif // JETREC_JETSPLITMERGETOOL_H