InvMassFinalState.hh

// -*- C++ -*-
#ifndef RIVET_InvMassFinalState_HH
#define RIVET_InvMassFinalState_HH
 
#include "Rivet/Projections/FinalState.hh"
 
namespace Rivet {
 

  class InvMassFinalState : public FinalState {
  public:

    InvMassFinalState(const FinalState& fsp,
                      const std::pair<PdgId, PdgId>& idpair, // pair of decay products
                      double minmass,                        // min inv mass
                      double maxmass,                        // max inv mass
                      double masstarget = -1.0);
 

    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
        double masstarget = -1.0);
 

    InvMassFinalState(const std::pair<PdgId, PdgId>& idpair, // pair of decay products
                      double minmass,                        // min inv mass
                      double maxmass,                        // max inv mass
                      double masstarget = -1.0);
    InvMassFinalState(
        const std::vector<std::pair<PdgId, PdgId>>& idpairs, // vector of pairs of decay products
        double minmass,                                      // min inv mass
        double maxmass,                                      // max inv mass
        double masstarget = -1.0);
 

    RIVET_DEFAULT_PROJ_CLONE(InvMassFinalState);

    using Projection::operator=;
 

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

    void useTransverseMass(bool usetrans = true) {
      _useTransverseMass = usetrans;
    }

    void calc(const Particles& inparticles);
 

    void project(const Event& e);

    CmpState compare(const Projection& p) const;
 
 
  protected:

    vector<PdgIdPair> _decayids;

    vector<pair<Particle, Particle>> _particlePairs;

    double _minmass;

    double _maxmass;

    double _masstarget;

    bool _useTransverseMass;
  };
 
 
}
 
#endif