PhotoelectricFittingProcessor.h

#ifdef USE_LCCD

#ifndef PHOTOELECTRICFITTINGPROCESSOR_H
#define PHOTOELECTRICFITTINGPROCESSOR_H 1

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

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

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

/* NOTE this uses the C++ version of Minuit - not the Fortran version */
#include <Minuit2/FCNBase.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 PhotoelectricFittingProcessor : public marlin::Processor
{

    public:

        PhotoelectricFittingProcessor() ;

        virtual void init() ;

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

        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;

        double _sigma;

        AIDA::ITuple *_photodot_fit_results_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::list< CLHEP::Hep2Vector > _photodot_centroids;

        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;

        class PhotodotLikelihood : public ROOT::Minuit2::FCNBase
        {

            public:

                PhotodotLikelihood(const std::set< EVENT::TrackerPulse* > &photodot_data);

                virtual double operator()(const std::vector<double>& x) const;

                virtual double Up() const;

            private:

                const std::set< EVENT::TrackerPulse* > &_photodot_data;
        };

};

}

#endif // PHOTOELECTRICFITTINGPROCESSOR_H
#endif // USE_LCCD