TrackFinderRectangularProcessor.h

#ifndef TRACK_FINDER_RECTANGULAR_PROCESSOR_H
#define TRACK_FINDER_RECTANGULAR_PROCESSOR_H

// C++
#include <string>
#include <vector>

// LCIO
#include <lcio.h>

// Marlin
#include <marlin/Processor.h>
#include <EVENT/LCEvent.h>
#include <EVENT/LCRunHeader.h>
#include <IMPL/TrackerHitImpl.h>
#include <EVENT/Track.h>
#include <gear/PadRowLayout2D.h>

namespace marlintpc
{

class TrackFinderRectangularProcessor : public marlin::Processor
{

    public:

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


        TrackFinderRectangularProcessor() ;

        virtual void init() ;

        virtual void processRunHeader(lcio::LCRunHeader* run) ;

        virtual void processEvent(lcio::LCEvent * evt) ;


        virtual void check(lcio::LCEvent * evt) ;


        virtual void end() ;


    protected:

        std::string _inputColName ;    
        std::string _outputColName ;   

        int _outputIsTransient;

        gear::PadRowLayout2D const * _padLayout;


        int _rowHeight;
        int _minTrackHits; 
        int _maxSkipRows;    

        float _deltaX, _deltaY, _deltaZ;

        int _nRun ;
        int _nEvt ;

        int _verbosityLevel;

        bool *isInTrack; 
        bool *isInTrackCandidate;  
        unsigned int nHits; 

        void zeroFlags(bool *IsIn);

        void addTrackCandidateFlags();


        void createNewIsInTrackFlags();

        IMPL::TrackerHitImpl extrapolateNextHit(double y);

        struct
        {
            float a, b, c, d;
        } straightLine;

        void doRegression(EVENT::Track *trackCandidate);
} ;

} // namespace martintpc

#endif//TRACK_FINDER_LINEAR_ROW_BASED_H