MinMassDiffInvJigsaw.hh

Go to the documentation of this file.

//   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 MinMassDiffInvJigsaw_HH
#define MinMassDiffInvJigsaw_HH

#include "Math/Minimizer.h"
#include "Math/Functor.h"

#include "RestFrames/InvisibleJigsaw.hh"

namespace RestFrames {

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

    virtual std::string GetLabel() const { return "Min #Sigma M(vis , inv)_{i}^{2}"; }

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

  private:
    const int m_Npair;
    mutable std::vector<TLorentzVector> m_Pvis;
    mutable std::vector<TLorentzVector> m_Pinv;
    std::vector<double> m_Minv;
    
    bool m_2D;
    TVector3 m_X;
    TVector3 m_Y;
    TVector3 m_Z;

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

    ROOT::Math::Minimizer* m_minimizer;
    ROOT::Math::Functor* m_functor;

    double EvaluateMetric(const double* param);

  };

}

#endif