MaxProbBreitWignerInvJigsaw.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 MaxProbBreitWignerInvJigsaw_HH
#define MaxProbBreitWignerInvJigsaw_HH

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

#include "RestFrames/InvisibleJigsaw.hh"

namespace RestFrames {

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

    virtual std::string GetLabel() const { return "Max Prob Breit-Wigner"; }

    virtual void SetPoleMass(double mass, int i = 0);
    virtual void SetWidth(double width, int i = 0);

    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;
    
    std::vector<double> m_Mass;
    std::vector<double> m_Width;

    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