LOGO

RestFrames  v1.0.1
RestFrames HEP Event Analysis Software Library
SetRapidityInvJigsaw.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 
31 #include "RestFrames/RestFrame.hh"
32 
33 namespace RestFrames {
34 
35  SetRapidityInvJigsaw::SetRapidityInvJigsaw(const std::string& sname,
36  const std::string& stitle) :
37  InvisibleJigsaw(sname, stitle, 1, 1)
38  {
39  m_Axis = RestFrame::GetAxis();
40  }
41 
42  SetRapidityInvJigsaw::SetRapidityInvJigsaw() : InvisibleJigsaw() {}
43 
44  SetRapidityInvJigsaw::~SetRapidityInvJigsaw() {}
45 
46  void SetRapidityInvJigsaw::SetAxis(const TVector3& axis){
47  m_Axis = axis;
48  }
49 
50  bool SetRapidityInvJigsaw::AnalyzeEvent(){
51  if(!IsSoundMind())
52  return SetSpirit(false);
53 
54  TLorentzVector inv_P = GetParentState().GetFourVector();
55  TLorentzVector vis_P = GetDependancyStates(0).GetFourVector();
56 
57  TVector3 boostZ = vis_P.BoostVector();
58  boostZ = boostZ.Dot(m_Axis.Unit())*m_Axis.Unit();
59 
60  inv_P.SetZ(0.0);
61  inv_P.Boost(boostZ);
62 
63  GetChildState(0).SetFourVector(inv_P);
64 
65  return SetSpirit(true);
66  }
67 
68 }