TrackFitterSimpleMinimizerProcessor.h

#ifndef TRACK_FITTER_SIMPLE_MINIMIZER_PROCESSOR_H
#define TRACK_FITTER_SIMPLE_MINIMIZER_PROCESSOR_H 1

#include "marlin/Processor.h"
//#include "lcio.h"
//#include <string>
#include <vector>
#include "TrackFitterFactory.h"

//#include <EVENT/TrackerHit.h>

//namespace AIDA {
//
//  class IHistogram1D;
//  class IHistogram2D;
//
//
//}


namespace marlintpc
{

class TrackFitterSimpleMinimizerProcessor : public marlin::Processor
{

    public:

        TrackFitterSimpleMinimizerProcessor();


        ~TrackFitterSimpleMinimizerProcessor();


        Processor * newProcessor();


        void init();


        void processRunHeader(lcio::LCRunHeader *run);


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


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


        virtual void end();


    protected:

        std::string _input_seed_tracks_collection_name;

        std::string _output_tracks_collection_name;

        bool _writeOutput; //< Flag wheter to write the output to storage

        int _maxIterations; //< Maximum number of iterations for the outlier rejection

        bool _rejectOutliers; //< Flag whether to reject outliers

        float _cutValue; //< Cut outliers wich are _cutValue*RMS away from the fitted track.

        EVENT::LCParameters * _trackCollectionParameters;

        std::vector<double> calculateMeanAndRMS( EVENT::Track *track, TrackFitterBase *fitter);
};

}// namespace marlintoc

#endif // TRACKSEEDERPROCESSOR_H