LOGO

RestFrames  v1.0.1
RestFrames HEP Event Analysis Software Library
MinMassDiffCombJigsaw.cc
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 #include "RestFrames/RestFrame.hh"
32 
33 namespace RestFrames {
34 
36  const std::string& stitle,
37  int N_comb, int N_mass)
38  : CombinatoricJigsaw(sname, stitle, N_comb, N_mass),
39  m_Ncomb(N_comb), m_Nmass(N_mass) {}
40 
41  MinMassDiffCombJigsaw::~MinMassDiffCombJigsaw() {}
42 
43  bool MinMassDiffCombJigsaw::EvaluateMetric(double& metric) const {
44  double diff = 0.;
45  double m1, m2;
46  TLorentzVector P1, P2;
47 
48  for(int i = 0; i < m_Nmass-1; i++){
49  if(!IsTransverse()){
50  int N = GetDependancyStates(i).GetN();
51  P1.SetPxPyPzE(0.,0.,0.,0.);
52  for(int j = 0; j < N; j++){
53  TLorentzVector v = GetDependancyStates(i)[j].GetFourVector();
54  TVector3 p = v.Vect() - v.Vect().Dot(GetTransverseAxis())*GetTransverseAxis();
55  v.SetVectM(p, v.M());
56  P1 += v;
57  }
58  } else {
59  P1 = GetDependancyStates(i).GetFourVector();
60  }
61  m1 = P1.M();
62  for(int j = i+1; i < m_Nmass; j++){
63  if(!IsTransverse()){
64  int N = GetDependancyStates(j).GetN();
65  P2.SetPxPyPzE(0.,0.,0.,0.);
66  for(int k = 0; k < N; k++){
67  TLorentzVector v = GetDependancyStates(j)[k].GetFourVector();
68  TVector3 p = v.Vect() - v.Vect().Dot(GetTransverseAxis())*GetTransverseAxis();
69  v.SetVectM(p, v.M());
70  P2 += v;
71  }
72  } else {
73  P2 = GetDependancyStates(j).GetFourVector();
74  }
75  m2 = P2.M();
76  diff += (m1-m2)*(m1-m2);
77  }
78  }
79 
80  metric = diff;
81  return true;
82  }
83 
84 }
85 
MinMassDiffCombJigsaw.hh
RestFrames::MinMassDiffCombJigsaw::MinMassDiffCombJigsaw
MinMassDiffCombJigsaw(const std::string &sname, const std::string &stitle, int N_comb, int N_mass)
Standard constructor.
Definition: MinMassDiffCombJigsaw.cc:35
RestFrames::CombinatoricJigsaw
Definition: CombinatoricJigsaw.hh:40
RestFrame.hh