PhotoelectricLikelihoodSeederProcessor.h

#ifndef PHOTOELECTRICLIKELIHOODSEEDERPROCESSOR_H
#define PHOTOELECTRICLIKELIHOODSEEDERPROCESSOR_H 1

// Marlin
#include <marlin/Processor.h>

// GEAR
#include <gear/TPCParameters.h>
#include <gearimpl/GlobalPadIndex.h>

// LCIO
#include <EVENT/TrackerHit.h>

// CLHEP
#include <CLHEP/Vector/TwoVector.h>

// stl
#include <string>
#include <set>
#include <list>

namespace EVENT
{

class TrackerPulse;
}

namespace AIDA
{

class ITuple;
}

namespace marlintpc
{

class PhotoelectricLikelihoodSeederProcessor : public marlin::Processor
{

    public:

        PhotoelectricLikelihoodSeederProcessor() ;

        virtual void init() ;

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

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

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

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

        virtual void end();

    private:

        std::string _input_photoelectric_tracker_pulses_collection_name;

        std::string _output_photodot_collection_name;

        std::string _output_photoline_collection_name;

        bool _output_is_transient;

        std::string _projected_photogeometry_file_name;

        bool _write_seed_tuple;

        double _sigma;


        AIDA::ITuple *_seed_position_tuple;

        std::map< gear::GlobalPadIndex, std::list< std::set< EVENT::TrackerPulse* > >::iterator > _photodot_pulse_set_by_global_index;

        std::list< std::set< EVENT::TrackerPulse* > > _photodot_pulse_sets;

        std::map< gear::GlobalPadIndex, std::list< std::set< EVENT::TrackerPulse* > >::iterator > _photoline_pulse_set_by_global_index;

        std::list< std::set< EVENT::TrackerPulse* > > _photoline_pulse_sets;
};

}

#endif // PHOTOELECTRICLIKELIHOODSEEDERPROCESSOR_H