PhotoelectricPulseExtractionProcessor.h

#ifndef PHOTOELECTRICPULSEEXTRACTIONPROCESSOR_H
#define PHOTOELECTRICPULSEEXTRACTIONPROCESSOR_H 1

#include "PhotoelectricXMLHelper.h"

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

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

// LCIO
#include <lcio.h>
#include <EVENT/TrackerPulse.h>

//ROOT
#include "TH2D.h"
#include "TH1D.h"

//AIDA
namespace AIDA
{

class ITuple;
}

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

typedef std::map<int, std::set<EVENT::TrackerPulse*> > photodotCluster;

namespace marlintpc
{

class PhotoelectricPulseExtractionProcessor : public marlin::Processor
{

    public:

        PhotoelectricPulseExtractionProcessor() ;

        virtual void init() ;

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

        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;

        std::string _output_photoelectric_tracker_pulses_collection_name_dots;

        std::string _output_photoelectric_tracker_pulses_collection_name_lines;

        float _arrival_time_width;

        bool _output_is_transient;

        std::string _projected_photogeometry_file_name;

        float _photopulse_interval_start;

        float _photopulse_interval_end;

        std::set< gear::GlobalPadIndex > _global_pad_indices_dots;

        std::set< gear::GlobalPadIndex > _global_pad_indices_lines;

        PhotoelectricXMLFile* _photodotFile; 

        /* @{ */
        std::string _distortion_file_name;

        photodotCluster _photodot_pulse_sets;

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

        int _runNumber; 

        std::string _detectorName; 

        std::vector<std::vector<marlintpc::Photodot>::const_iterator> _dots; 

        double _sigma; 
        
        bool _fitSigma; 

        TH2D * _numberOfPadsWithPulsesPerCluster;
        TH2D * _numberOfPulsesPerCluster;
        AIDA::ITuple *_pulse_tuple;
        AIDA::ITuple *_cluster_tuple;
        /* @} */
};

}

#endif // PHOTOELECTRICPULSEEXTRACTIONPROCESSOR_H