MinMassesSqInvJigsaw.hh

//   RestFrames: particle physics event analysis library
//   --------------------------------------------------------------------
//   Copyright (c) 2014-2016, Christopher Rogan
//   This file is part of RestFrames.
//
//   RestFrames is free software; you can redistribute it and/or modify
//   it under the terms of the GNU General Public License as published by
//   the Free Software Foundation; either version 2 of the License, or
//   (at your option) any later version.
// 
//   RestFrames is distributed in the hope that it will be useful,
//   but WITHOUT ANY WARRANTY; without even the implied warranty of
//   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//   GNU General Public License for more details.
// 
//   You should have received a copy of the GNU General Public License
//   along with RestFrames. If not, see <http://www.gnu.org/licenses/>.

#ifndef MinMassesSqInvJigsaw_HH
#define MinMassesSqInvJigsaw_HH

#include "RestFrames/InvisibleJigsaw.hh"

namespace RestFrames {

  class MinMassesSqInvJigsaw : public InvisibleJigsaw {
  public:
    MinMassesSqInvJigsaw(const std::string& sname, 
               const std::string& stitle,
               int N_vis_inv_pair);
    MinMassesSqInvJigsaw();
    ~MinMassesSqInvJigsaw();

    virtual std::string GetLabel() const { return "Min Masses Sq."; }

    virtual double GetMinimumMass() const;
    
    virtual bool AnalyzeEvent();

    static MinMassesSqInvJigsaw& Empty();


  private:
    const int m_Npair;
    mutable std::vector<TLorentzVector> m_Pvis;
    mutable std::vector<TLorentzVector> m_Pinv;
    std::vector<double>        m_Minv;

    double GetPScale(double Minv);
    void ApplyOptimalRotation();

    static MinMassesSqInvJigsaw m_Empty;
  };

}

#endif