HitSplitter.h

#ifndef MARLINTPC_HIT_SPLITTER_H
#define MARLINTPC_HIT_SPLITTER_H

//stl
#include <set>

// lcio
#include <IMPL/LCCollectionVec.h>
#include <EVENT/TrackerHit.h>
#include <EVENT/TrackerPulse.h>
#include <EVENT/LCIO.h>

namespace marlintpc
{
class HitSplitter
{
    public:
        HitSplitter(IMPL::LCCollectionVec * SplitHitsCollection,
                    IMPL::LCCollectionVec * PulsesInSplitHitsCollection,
                    float vDrift);

        EVENT::TrackerHitVec splitHit(EVENT::TrackerHit const * const);

    protected:

        IMPL::LCCollectionVec * _splitHitsCollection;
        IMPL::LCCollectionVec * _pulsesInSplitHitsCollection;
        float _vDrift; 

        class SortablePulse
        {
            public:
                bool operator<(const SortablePulse & p) const;
                int operator-(const SortablePulse & p) const;

                SortablePulse(EVENT::TrackerPulse * pulse);
                EVENT::TrackerPulse * getTrackerPulse() const;
            protected:
                EVENT::TrackerPulse * _pulse;
        };

        typedef std::pair < std::multiset<SortablePulse>::iterator,
        std::multiset<SortablePulse>::iterator > IteratorPair;

        struct Distance
        {
            float distance;
            std::multiset<SortablePulse>::iterator firstPulse;
            std::multiset<SortablePulse>::iterator secondPulse;

            Distance(float d, std::multiset<SortablePulse>::iterator first,
                     std::multiset<SortablePulse>::iterator  second);

            bool operator<(Distance const &right) const
            {
                return distance < right.distance;
            }
        };

};
}// namespace marlintpc

#endif // MARLINTPC_HIT_SPLITTER_H