SatoruJetFinderProcessor.h

/* -*- Mode: C++; tab-width: 2; indent-tabs-mode: nil; c-basic-offset: 2 -*- */
/*
** This file is part of the MarlinReco Project.
** Forming part of the SubPackage: SatoruJetFinder.
**
** For the latest version download from Web CVS:
** http://www-zeuthen.desy.de/lc-cgi-bin/cvsweb.cgi/marlinreco/?cvsroot=MarlinReco
**
** $Id$
**
**
*/ 

#ifndef SatoruJetFinderProcessor_h
#define SatoruJetFinderProcessor_h 1

#include "marlin/Processor.h"

#include "lcio.h"
#include "IMPL/LCCollectionVec.h" 


using namespace lcio ;


struct SatoruPartonArray{
    int NumberOfPartons;
    float Momentum[1200];
    int PointerParticleToJets[300];
}; 
struct SatoruJetsArray{
    int NumberOfJets;
    float Momentum[80];
};

extern "C" {
    void syjkrn_( const char* GlobalModus_,
                 // Primary jet finding
                 int &NJetRequested_,float &Threshold_,
                 int &PrimaryJetFindingMode_,float* YCut_,
                 //  Reassosiation after first jet finding
                 int &MergingMode_,float &MergingThreshold_,
                 // Second Jet finding mode
                 int &SecondJetFindingMode_,
                 // input array--> array of 
                 // PPartons(DimensionOfInputArray(4-6),NumberOfPartons)
                 int &NumberOfPartons_,
                 int &DimensionOfInputArray_,float *PPartons_,
                 // Output 
                 int &MaximalNumberOfJets_,
                 int &NJetsFound_,int *PointerParicleToJet_,
                 int &DimensionOfOutputArray_,float *PJets_,
                 float &YMinus_,float &YPlus_,
                 // error flag
                 int &IError_,int GlobalModusLenght_)
        ;}


namespace marlin
{


  class SatoruJetFinderProcessor : public Processor {

  
  public:
    
    virtual Processor*  newProcessor() { return new SatoruJetFinderProcessor ; }
    
  
    SatoruJetFinderProcessor() ;
    
    virtual void init() ;
    
    virtual void processRunHeader( LCRunHeader* run ) ;

    virtual void processEvent( LCEvent * evt ) ; 
    

    virtual void end() ;



  private:

    std::string _inputCollection;
    std::string _outputCollection;
    std::string _successTag;
    bool        _writeTag;
    std::string _jetFindingMode;
    SatoruPartonArray _partonsWorkArray;
    SatoruJetsArray _jetsWorkArray;
    //    int _nJetVal;
    std::string _globalMode;
    int _nJetRequested;
    float _threshold;
    int _primaryJetFindingMode;
    float _yCutParam;
    float _rConeParam;
    float _epsConeParam;
    float _yCut[2];
    int _mergingMode;
    float _mergingThreshold;
    int _secondJetFindingMode;
    int _debug;
    float _YMinus;
    float _YPlus;
    void putPartons(LCEvent * evt);
    void callSatoru(LCEvent * evt);
    void goSatoru(LCEvent * evt,LCCollection *JetsCol);                     
    void getJets(LCEvent * evt,LCCollection *JetsCol);

    /* debug routines */
    void writePartons();
    void writeJets();
    void writeParameters();

    // true constants for mode selection
    struct MD 
    {
      const static int JADE_E0   = 1;
      const static int JADE_P    = 2;
      const static int JADE_P0   = 3;
      const static int JADE_E    = 4;
      const static int DURHAM    = 5;
      const static int GENEVA    = 6;
      const static int CAMBRIDGE = 7;
      const static int LUCLUS_1  = 8;
      const static int LUCLUS_2  = 9;
      const static int LUCLUS_3  = 10;
      const static int LUCLUS_4  = 11;
      const static int CONE      = 12;
    };

  } ;



} // namespace marlin
#endif