Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:34:49

0001 #include <cmath>
0002 #include <cassert>
0003 
0004 #include "RecoJets/FFTJetAlgorithms/interface/JetToPeakDistance.h"
0005 
0006 namespace fftjetcms {
0007   JetToPeakDistance::JetToPeakDistance(const double etaToPhiBandwidthRatio)
0008       : etaBw_(sqrt(etaToPhiBandwidthRatio)), phiBw_(1.0 / etaBw_) {
0009     assert(etaToPhiBandwidthRatio > 0.0);
0010   }
0011 
0012   double JetToPeakDistance::operator()(const fftjet::RecombinedJet<VectorLike>& j1, const fftjet::Peak& peak) const {
0013     if (peak.membershipFactor() <= 0.0)
0014       // This peak essentially does not exist...
0015       return 2.0e300;
0016 
0017     const double deta = (j1.vec().Eta() - peak.eta()) / etaBw_;
0018     double dphi = j1.vec().Phi() - peak.phi();
0019     if (dphi > M_PI)
0020       dphi -= (2.0 * M_PI);
0021     else if (dphi < -M_PI)
0022       dphi += (2.0 * M_PI);
0023     dphi /= phiBw_;
0024     return sqrt(deta * deta + dphi * dphi);
0025   }
0026 }  // namespace fftjetcms