MinMassesSqCombJigsaw.cc

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/>.
 
#include "RestFrames/RestFrame.hh"
#include "RestFrames/MinMassesSqCombJigsaw.hh"
 
namespace RestFrames {
 
  MinMassesSqCombJigsaw::MinMassesSqCombJigsaw(const std::string& sname, 
                           const std::string& stitle,
                           int N_comb, int N_mass) 
    : CombinatoricJigsaw(sname, stitle, N_comb, N_mass),
      m_Ncomb(N_comb), m_Nmass(N_mass) {}
  
  MinMassesSqCombJigsaw::~MinMassesSqCombJigsaw() {}
 
  bool MinMassesSqCombJigsaw::EvaluateMetric(double& metric) const {
    double sum = 0.;
 
    if(!IsTransverse()){
      
      for(int i = 0; i < m_Nmass; i++)
    sum += GetDependancyStates(i).GetFourVector().M2();
      
    } else { // set transverse
      
      for(int i = 0; i < m_Nmass; i++){
    TLorentzVector CM;
    int N = GetDependancyStates(i).GetN();
    for(int j = 0; j < N; j++){
      TLorentzVector v = GetDependancyStates(i)[j].GetFourVector();
      TVector3 p = v.Vect() - v.Vect().Dot(GetTransverseAxis())*GetTransverseAxis();
      v.SetVectM(p, v.M());
      CM += v;
    }
    sum += CM.M2();
      }
 
    }
    
    metric = sum;
    return true;
  }
 
}