SetRapidityInvJigsaw.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/SetRapidityInvJigsaw.hh"
#include "RestFrames/RestFrame.hh"
 
namespace RestFrames {
 
  SetRapidityInvJigsaw::SetRapidityInvJigsaw(const std::string& sname, 
                         const std::string& stitle) : 
    InvisibleJigsaw(sname, stitle, 1, 1)
  {
    m_Axis = RestFrame::GetAxis();
    m_InvMassDependancy = true;
  }
 
  SetRapidityInvJigsaw::SetRapidityInvJigsaw() : InvisibleJigsaw() {}
 
  SetRapidityInvJigsaw::~SetRapidityInvJigsaw() {}
 
  void SetRapidityInvJigsaw::SetAxis(const TVector3& axis){
    m_Axis = axis;
  }
 
  bool SetRapidityInvJigsaw::AnalyzeEvent(){
    if(!IsSoundMind()) 
      return SetSpirit(false);
    
    TLorentzVector inv_P = GetParentState().GetFourVector();
    TLorentzVector vis_P = GetDependancyStates(0).GetFourVector();
 
    TVector3 v_inv_P = inv_P.Vect();
    v_inv_P = v_inv_P - v_inv_P.Dot(m_Axis.Unit())*m_Axis.Unit();
    inv_P.SetVectM(v_inv_P, inv_P.M());
 
    if(vis_P.E() > 1e-8){
      TVector3 boostZ = vis_P.BoostVector();
      boostZ = boostZ.Dot(m_Axis.Unit())*m_Axis.Unit();
      if(!(boostZ.Mag() >= 1.))
    inv_P.Boost(boostZ);
    }
 
    GetChildState(0).SetFourVector(inv_P);
    
    return SetSpirit(true);
  }
 
}