PhotoelectricProjectionFinder.h

#ifdef USE_LCCD

#ifndef PHOTOELECTRICPROJECTIONFINDER_H
#define PHOTOELECTRICPROJECTIONFINDER_H 1

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

// LCIO
#include <lcio.h>

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

// stl
#include <map>
#include <set>
#include <string>
#include <vector>

namespace CLHEP
{

class Hep2Vector;
}

namespace EVENT
{

class TrackerPulse;
}

namespace marlintpc
{

class PhotoelectricProjectionFinder : public marlin::Processor
{

    public:

        PhotoelectricProjectionFinder() ;

        virtual void init() ;

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

        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_tracker_pulses_collection_name;

        float _max_separation;

        float _density_threshold;

        std::string _projected_photogeometry_file_name;

        std::string _photogeometry_file_name;

        std::map< gear::GlobalPadIndex, int > _counts_by_pad;

        int _total_event_count;

        float _photopulse_interval_start;

        float _photopulse_interval_end;

        float _event_fraction_threshold;

        float _maximum_dot_displacement;

        std::string _distortion_file_name;

        int _runNumber; 

        std::string _detectorName; 

        void create_points(std::vector<CLHEP::Hep2Vector*> &points, std::map< CLHEP::Hep2Vector*, gear::GlobalPadIndex> &pad_by_point);

        void find_neighbours(const std::vector<CLHEP::Hep2Vector*> &points, const float max_seperation, std::map< CLHEP::Hep2Vector*, std::set<CLHEP::Hep2Vector*> > &neighbours_by_point);

        void calculate_point_clusters(const std::vector< CLHEP::Hep2Vector* > &points, const std::map< CLHEP::Hep2Vector*, std::set< CLHEP::Hep2Vector* > > &neighbours_by_point, std::set< std::set< CLHEP::Hep2Vector* > > &point_clusters);

        void classify_point_clusters(const std::set< std::set< CLHEP::Hep2Vector* > > &point_clusters, const std::map< CLHEP::Hep2Vector*, gear::GlobalPadIndex > pad_by_point, std::set< const std::set< CLHEP::Hep2Vector*>* > &photodots, std::set< const std::set< CLHEP::Hep2Vector* >* > &photolines);

        void store_results(const std::set< const std::set< CLHEP::Hep2Vector* >* > &photodots, const std::set< const std::set< CLHEP::Hep2Vector* >* > &photolines, std::map< CLHEP::Hep2Vector*, gear::GlobalPadIndex > pad_by_point, bool storeLCIO);

        void produce_tuples(const std::set< const std::set< CLHEP::Hep2Vector* >* > &photodots, const std::set< const std::set< CLHEP::Hep2Vector* >* > &photolines);
};

}

#endif // PHOTOELECTRICPROJECTIONFINDER_H
#endif // USE_LCCD