MomentumHelpers.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
 
 
#ifndef COLUMNAR_CORE_MOMENTUM_HELPERS_H
#define COLUMNAR_CORE_MOMENTUM_HELPERS_H
 
#include <ColumnarCore/MomentumAccessors.h>
#include <FourMomUtils/xAODP4Helpers.h>
 
namespace columnar
{
 
  inline
  double deltaPhi( double phiA, double phiB )
  {
    return  -remainder( -phiA + phiB, 2*M_PI );
  }
 
 
 
  template<ContainerIdConcept CI1,ContainerIdConcept CI2,typename CM>
  double deltaRapidity( const MomentumAccessors<CI1,CM>& momAcc1, ObjectId<CI1,CM> p1, const MomentumAccessors<CI2,CM>& momAcc2, ObjectId<CI2,CM> p2, bool useRapidity=true )
  {
    if (useRapidity)
      return momAcc1.rapidity(p1) - momAcc2.rapidity(p2);
    else
      return momAcc1.eta(p1) - momAcc2.eta(p2);
  }
 
 
 
  template<ContainerIdConcept CI1,ContainerIdConcept CI2,typename CM>
  bool isInDeltaR( const MomentumAccessors<CI1,CM>& momAcc1, ObjectId<CI1,CM> p1, const MomentumAccessors<CI2,CM>& momAcc2, ObjectId<CI2,CM> p2,
                    double dR, bool useRapidity=true )
  {
    const double dPhi = std::abs( xAOD::P4Helpers::deltaPhi(momAcc1.phi(p1),momAcc2.phi(p2)) ); // in [-pi,pi)
    if ( CxxUtils::fpcompare::greater(dPhi,dR) ) return false;        // <==
    const double dRapidity = std::abs( deltaRapidity(momAcc1,p1,momAcc2,p2,useRapidity) );
    if ( CxxUtils::fpcompare::greater(dRapidity,dR) ) return false;        // <==
    const double deltaR2 = dRapidity*dRapidity + dPhi*dPhi;
    if ( CxxUtils::fpcompare::greater(deltaR2,dR*dR) ) return false;  // <==
    return true;
  }
}
 
#endif