HitCandidate.h

#ifndef HIT_CANDIDATE_H
#define HIT_CANDIDATE_H

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

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

//stl
#include <vector>

namespace marlintpc
{
class HitCandidate
{
    public:
        HitCandidate(const gear::TPCParameters & tpcParams, double vDrift,
                     EVENT::TrackerPulse * pulse = 0);

        void addPulse(EVENT::TrackerPulse*);
        void addDeadPad();

        void setBoundaryFlag();


        const EVENT::TrackerPulse* getMaximumPulse() const
        {
            return _maxPulse;
        };

        double getCharge();

        double getMaxPulseHeight() const;

        const std::vector<EVENT::TrackerPulse*>& getPulses() const
        {
            return _thePulses;
        };

        unsigned int getWidth() const;

        bool   isAtBoundary() const
        {
            return _isAtBoundary;
        };

        double getLocalFixedCoordinate();

        double getLocalFreeCoordinate();

        double getZCoordinate();

        double getLocalFixedVariance();

        double getLocalFreeVariance();

        double getLocalSkewness();

        double getZVariance();

        lcio::DoubleVec getLocalCoordinates();

        lcio::DoubleVec getGlobalCoordinates();

        lcio::DoubleVec getGlobalCartesianCoordinates();

        lcio::FloatVec getGlobalCovarianceMatrix();

        lcio::FloatVec getGlobalCartecianCovarianceMatrix();

        IMPL::TrackerHitImpl * createTrackerHit();

    protected:
        std::vector<EVENT::TrackerPulse*> _thePulses;
        EVENT::TrackerPulse* _maxPulse;
        double _maxPulseCharge;
        int _numberOfDeadPads;
        bool _isAtBoundary;
        double _vDrift;
        double _zAnode;

        const gear::TPCParameters & _tpcParameters;
        const gear::TPCModule * _tpcModule;

        // the parameters calculated by calculateHitProperties

        lcio::DoubleVec _localCoordinates;
        lcio::DoubleVec _globalCoordinates;
        lcio::FloatVec _globalCovariance;

        double _charge;
        double _localFixedVariance;
        double _localFreeVariance;
        double _localSkewness; 
        double _zVariance;

        void calculateHitProperties();

        bool _hitParametersAreValid; 
};

}

#endif //HIT_CANDIDATE_H