TopoFinder.h

#ifndef TOPO_FINDER_H
#define TOPO_FINDER_H 1

#include <map>
#include <list>

#include <lcio.h>
#include <EVENT/TrackerHit.h>
#include <EVENT/TrackerPulse.h>
#include <gear/TPCParameters.h>
#include <gear/PadRowLayout2D.h>
#include <gearimpl/GlobalPadIndex.h>

namespace marlintpc
{

class TopoFinder
{

    protected:

        TopoFinder();

//    void setPulseMap(std::multimap<int, EVENT::TrackerPulse*> *map);

        EVENT::TrackerHit* findHit(
            std::multimap<gear::GlobalPadIndex, EVENT::TrackerPulse*>::iterator currentPulseIter);

        std::list<EVENT::TrackerPulse *> findNeighbouringPulses(
            gear::GlobalPadIndex currentPadIndex,
            EVENT::TrackerPulse *  currentPulse,
            const gear::TPCModule &module);

        EVENT::TrackerHitVec *findHitsInNeighbouringRows(EVENT::TrackerHit* hit, int moduleID);

        EVENT::TrackerHitVec *findHitsInRange(gear::GlobalPadIndex index1, gear::GlobalPadIndex index2,
                                              float minADCIndex, float maxADCIndex);

        int _maxEmptyPads;   
        int _minHitSize; 
        float _maxTimeSpread;
        int _minTrackHits; 
        int _maxSkipRows;    

        float _vDrift; 

        //  /// The pad layout as defined in GEAR
        //    int _padLayoutType; ///< Pad layout type (PadRowLayout2D::CARTESIAN or PadRowLayout2D::POLAR)
        //    gear::PadRowLayout2D const *_padLayout;

        float _zAnode;

        int _negativeZDirection;

        std::multimap<gear::GlobalPadIndex, EVENT::TrackerPulse*> *_pulsesByIndex;

//    /**
//     * A nested helper class for sorting gear::GlobalPadIndex.
//     * They are sorted by moduleID, and within moduleID in padIndex
//     */
//    class  sortGlobalIndex {
//    public:
//      bool operator()(gear::GlobalPadIndex index1, gear::GlobalPadIndex index2) {
//  if (index1.getModuleID() == index2.getModuleID() )
//    { // the modules IDs are equal, distiguish using pad index
//      return ( index1.getPadIndex() < index2.getPadIndex() );
//    }
//  else
//    { // the modules IDs are equal, distiguish using pad index
//      return ( index1.getModuleID() < index2.getModuleID() );
//    }
//      }
//    };
        unsigned int _recursionCounter;

    public:
        static void addTrackerHitVecs(EVENT::TrackerHitVec **vectorToAddTo, EVENT::TrackerHitVec **vectorToBeAdded);

        std::string getRevision() ;

};

}// namespace marlintpc

#endif //TOPO_FINDER_H