LOGO

RestFrames  v1.0.1
RestFrames HEP Event Analysis Software Library
MaxProbBreitWignerInvJigsaw.hh
Go to the documentation of this file.
1 // RestFrames: particle physics event analysis library
3 // --------------------------------------------------------------------
4 // Copyright (c) 2014-2016, Christopher Rogan
14 // This file is part of RestFrames.
15 //
16 // RestFrames is free software; you can redistribute it and/or modify
17 // it under the terms of the GNU General Public License as published by
18 // the Free Software Foundation; either version 2 of the License, or
19 // (at your option) any later version.
20 //
21 // RestFrames is distributed in the hope that it will be useful,
22 // but WITHOUT ANY WARRANTY; without even the implied warranty of
23 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
24 // GNU General Public License for more details.
25 //
26 // You should have received a copy of the GNU General Public License
27 // along with RestFrames. If not, see <http://www.gnu.org/licenses/>.
29 
30 #ifndef MaxProbBreitWignerInvJigsaw_HH
31 #define MaxProbBreitWignerInvJigsaw_HH
32 
33 #include "Math/Minimizer.h"
34 #include "Math/Functor.h"
35 
37 
38 namespace RestFrames {
39 
41  public:
42  MaxProbBreitWignerInvJigsaw(const std::string& sname,
43  const std::string& stitle,
44  int N_vis_inv_pair);
46 
47  virtual std::string GetLabel() const { return "Max Prob Breit-Wigner"; }
48 
49  virtual void SetPoleMass(double mass, int i = 0);
50  virtual void SetWidth(double width, int i = 0);
51 
52  virtual double GetMinimumMass() const;
53 
54  virtual bool AnalyzeEvent();
55 
56  private:
57  const int m_Npair;
58  mutable std::vector<TLorentzVector> m_Pvis;
59  mutable std::vector<TLorentzVector> m_Pinv;
60  std::vector<double> m_Minv;
61 
62  std::vector<double> m_Mass;
63  std::vector<double> m_Width;
64 
65  bool m_2D;
66  TVector3 m_X;
67  TVector3 m_Y;
68  TVector3 m_Z;
69 
70  double GetPScale(double Minv);
71  void ApplyOptimalRotation();
72 
73  ROOT::Math::Minimizer* m_minimizer;
74  ROOT::Math::Functor* m_functor;
75 
76  double EvaluateMetric(const double* param);
77 
78  };
79 
80 }
81 
82 #endif