FinalPartons.cc

Go to the documentation of this file.

// -*- C++ -*-

#include "Rivet/Projections/FinalPartons.hh"
#include "Rivet/Tools/RivetHepMC.hh"

namespace Rivet {


    bool FinalPartons::accept(const Particle& p) const {

        // reject if *not* a parton
        if (!isParton(p))
            return false;

        // accept partons if they end on a standard hadronization vertex
        if (p.genParticle()->end_vertex() != NULL &&
                p.genParticle()->end_vertex()->id() == 5)
          return true;

        // reject if p has a parton child.
        foreach (const Particle& c, p.children()) {
            if (isParton(c))
                return false;
        }

        // reject if from a hadron or tau decay.
        if (p.fromDecay())
            return false;

        return _cuts->accept(p);
    }


    void FinalPartons::project(const Event& e) {
        _theParticles.clear();

        foreach (const GenParticle* gp, Rivet::particles(e.genEvent())) {
            if (!gp) continue;

            const Particle p(gp);
            if (accept(p)) _theParticles.push_back(p);
        }
    }


}