InvMassFinalState.hh

Go to the documentation of this file.

// -*- C++ -*-
#ifndef RIVET_InvMassFinalState_HH
#define RIVET_InvMassFinalState_HH

#include "Rivet/Projections/FinalState.hh"

namespace Rivet {


  /// Identify particles which can be paired to make an invariant mass within a given mass window
  class InvMassFinalState : public FinalState {
  public:
 
    // Constructor for a single inv-mass pair
    InvMassFinalState(const FinalState& fsp,
                      const std::pair<PdgId, PdgId>& idpair, // pair of decay products
                      double minmass, // min inv mass
                      double maxmass); // max inv mass


    InvMassFinalState(const FinalState& fsp,
                      const std::vector<std::pair<PdgId, PdgId> >& idpairs,  // vector of pairs of decay products
                      double minmass, // min inv mass
                      double maxmass); // max inv mass
 
 
    /// Clone on the heap.
    virtual const Projection* clone() const {
        return new InvMassFinalState(*this);
    }


  public:

    /// Constituent pairs
    const std::vector<std::pair<Particle, Particle> >& particlePairs() const;
        

  protected:
 
    /// Apply the projection on the supplied event.
    void project(const Event& e);
 
    /// Compare projections.
    int compare(const Projection& p) const;


  private:
 
    /// Handy typedef for a pair of PID codes
    typedef pair<PdgId,PdgId> PidPair;

    /// IDs of the decay products
    std::vector<PidPair> _decayids;

    /// Constituent pairs
    std::vector<std::pair<Particle, Particle> > _particlePairs;

    /// Min inv mass
    double _minmass;

    /// Max inv mass
    double _maxmass;
 
  };


}


#endif