QCDSmear.h

#ifndef FastJetStability_QCDSmear_h
#define FastJetStability_QCDSmear_h

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

namespace stab {

class QCDSmear: public QCDSmearSequenceCallback {

public:
        inline explicit QCDSmear(const double cut,
                                      const AlphaStrong& alphaS,
                                      const bool doUniformPhaseSpace = true)
            : cutoff_(cut),
              alphaS_(alphaS),
              uniformPhaseSpace_(doUniformPhaseSpace) {}

        inline virtual ~QCDSmear() {}

        inline double getCutoff() const {return cutoff_;}
        inline const AlphaStrong& getAlphaS() const {return alphaS_;}

        bool operator()(double prob, unsigned recoNum,
                        unsigned ihard, unsigned isoft,
                        PseudoJet& hard, PseudoJet& soft);

private:
        virtual bool passesCut(double prob, unsigned recoNum,
                               unsigned ihard, unsigned isoft,
                               const PseudoJet& hard,
                               const PseudoJet& soft) const;
        double cutoff_;
        AlphaStrong alphaS_;
        bool uniformPhaseSpace_;

}; // end QCDSmear class declaration

} // end stab namespace

#endif