QCDSmearSequenceCallback.h

#ifndef FastJetStability_QCDSmearSequenceCallback_h
#define FastJetStability_QCDSmearSequenceCallback_h

#include "JetStabilityDefs.h"

namespace stab {

   struct QCDSmearSequenceCallback
    {
        inline virtual ~QCDSmearSequenceCallback() {}

        //
        // The callback arguments are as follows:
        //
        // prob               -- g -> gg splitting probability from QCD
        //
        // combinationNumber  -- Sequential number in pairwise combinations
        //                       of particles. These combinations are arranged
        //                       in the order of decreasing probability.
        //                       Once considered, the combinations containing
        //                       the softer particle are no longes taken
        //                       into consideration (as if the softer
        //                       particle is discarded)
        //
        // hardNumber         -- Number of the harder particle in the
        //                       event
        //
        // softNumber         -- Number of the softer particle in the
        //                       event (for each event, all "softNumber"
        //                       arguments should be unique).
        //
        // hard, soft         -- Four-momenta of the harder and softer
        //                       particles, respectively
        //
        // Note that the callback can actually modify "hard" and "soft"
        // arguments. The probabilities given by "prob" will not be
        // recalculated in this case.
        //
        // The callback should return "true" if the combination sequence
        // should be continued and "false" if it should be stopped.
        //
        virtual bool operator()(double prob, unsigned combinationNumber,
                                unsigned hardNumber, unsigned softNumber,
                                PseudoJet& hard, PseudoJet& soft) = 0;
    };


} // end stab namespace

#endif