QCDSmearSequence.h

#ifndef FastJetStability_QCDSmearSequence_h
#define FastJetStability_QCDSmearSequence_h

#include <vector>
#include <cassert>

#include "npstat/nm/Triple.hh" 

#include "AbsSmearingAgent.h"
#include "QCDSmearSequenceCallback.h"
#include "StandardModel.h"

namespace stab {

    class QCDSmearSequence: public AbsSmearingAgent
    {
    public:
        explicit QCDSmearSequence(const AlphaStrong& alphaS);

        // This object will not assume the ownership of the callbacks.
        QCDSmearSequence(const AlphaStrong& alphaS,
                            const std::vector<QCDSmearSequenceCallback*> cb);

    // initialize
    inline virtual void initialize(const Event &event, const StabilityEventInfo &info) { m_event = event; }

        virtual void smear(Event& event);
        
        // The following function returns the number of callback
        // invocations which returned "true"
    inline unsigned nCallBacks() const { return m_callBackCount; }

        // Add a callback. This object will not assume the callback ownership.
        inline void addCallback(QCDSmearSequenceCallback* ptr)
            {assert(ptr); callbacks_.push_back(ptr);}


    private:
        typedef npstat::Triple<double,unsigned,unsigned> ProbTriple;
    
        QCDSmearSequence();

        AlphaStrong alphaS_;
        std::vector<QCDSmearSequenceCallback*> callbacks_;

        mutable std::vector<ProbTriple> probs_;
        mutable std::vector<double> ptValues_;
        mutable std::vector<int> used_;

        // number of callbac invocations returning true from previous smear 
    unsigned m_callBackCount;

    Event m_event;

    };


} // end stab namespace

#endif