MinMassDiffCombJigsaw.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/MinMassDiffCombJigsaw.hh"
 
namespace RestFrames {
 
  MinMassDiffCombJigsaw::MinMassDiffCombJigsaw(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) {}
  
  MinMassDiffCombJigsaw::~MinMassDiffCombJigsaw() {}
 
  bool MinMassDiffCombJigsaw::EvaluateMetric(double& metric) const {
    double diff = 0.;
    double m1, m2;
    TLorentzVector P1, P2;
    
    for(int i = 0; i < m_Nmass-1; i++){
      if(!IsTransverse()){
    int N = GetDependancyStates(i).GetN();
    P1.SetPxPyPzE(0.,0.,0.,0.);
    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());
      P1 += v;
    }
      } else {
    P1 = GetDependancyStates(i).GetFourVector();
      }
      m1 = P1.M();
      for(int j = i+1; i < m_Nmass; j++){
    if(!IsTransverse()){
      int N = GetDependancyStates(j).GetN();
      P2.SetPxPyPzE(0.,0.,0.,0.);
      for(int k = 0; k < N; k++){
        TLorentzVector v = GetDependancyStates(j)[k].GetFourVector();
        TVector3 p = v.Vect() - v.Vect().Dot(GetTransverseAxis())*GetTransverseAxis();
        v.SetVectM(p, v.M());
        P2 += v;
      }
    } else {
      P2 = GetDependancyStates(j).GetFourVector();
    }
    m2 = P2.M();
    diff += (m1-m2)*(m1-m2);
      }
    }
    
    metric = diff;
    return true;
  }
 
}