#include <iostream>
#include "TLorentzVector.h"

void test()
{
   double pt = 5366.49744762, eta = 0.305951957105, phi = 1.619407720288;
   double pt2 = pt;
   
   TLorentzVector a,b;
   
   a.SetPtEtaPhiM( pt, eta, phi, 0 );
   b = a;
   
   
   if( pt < pt2 )
     std::cout << "pt < pt2" << std::endl;
   if( pt > pt2 )
     std::cout << "pt > pt2" << std::endl;
   if( pt == pt2 )
     std::cout << "pt == pt2" << std::endl;
     
   if( a.Pt() < b.Pt() )
     std::cout << "a.Pt() < b.Pt()" << std::endl;
   if( a.Pt() > b.Pt() )
     std::cout << "a.Pt() > b.Pt()" << std::endl;
   if( a.Pt() == b.Pt() )
     std::cout << "a.Pt() == b.Pt()" << std::endl;

   if( a.Pt() < a.Pt() )
     std::cout << "a.Pt() < a.Pt()" << std::endl;
   if( a.Pt() > a.Pt() )
     std::cout << "a.Pt() > a.Pt()" << std::endl;
   if( a.Pt() == a.Pt() )
     std::cout << "a.Pt() == a.Pt()" << std::endl;

}