1 | //---------------------------------------------------------------------- |
---|
2 | // This file distributed with FastJet has been obtained from SpartyJet |
---|
3 | // v2.20.0 by Pierre-Antoine Delsart, Kurtis L. Geerlings, Joey |
---|
4 | // Huston, Brian T. Martin and Chris Vermilion |
---|
5 | // For details, see http://www.pa.msu.edu/~huston/SpartyJet/ |
---|
6 | // http://projects.hepforge.org/spartyjet/ |
---|
7 | // |
---|
8 | // Changes from the original file are listed below. |
---|
9 | //---------------------------------------------------------------------- |
---|
10 | |
---|
11 | // History of changes from the original JetDistance.hh file in |
---|
12 | // SpartyJet v2.20 |
---|
13 | // |
---|
14 | // 2009-01-15 Gregory Soyez <soyez@fastjet.fr> |
---|
15 | // |
---|
16 | // * put the code in the fastjet::atlas namespace |
---|
17 | |
---|
18 | #ifndef JETUTIL_JETDISTANCES_H |
---|
19 | #define JETUTIL_JETDISTANCES_H |
---|
20 | |
---|
21 | |
---|
22 | #include "Jet.hh" |
---|
23 | |
---|
24 | #include <cmath> |
---|
25 | |
---|
26 | #include <fastjet/internal/base.hh> |
---|
27 | |
---|
28 | FASTJET_BEGIN_NAMESPACE |
---|
29 | |
---|
30 | namespace atlas { |
---|
31 | |
---|
32 | struct JetDistances { |
---|
33 | |
---|
34 | |
---|
35 | // distance in eta |
---|
36 | inline static double deltaEta(const Jet& jet1, const Jet& jet2) |
---|
37 | { return jet1.eta() - jet2.eta(); } |
---|
38 | inline static double deltaEta(const Jet* jet1, |
---|
39 | const Jet* jet2) |
---|
40 | { return jet1->eta() - jet2->eta(); } |
---|
41 | inline static double deltaEta(const double eta1, |
---|
42 | const double eta2) |
---|
43 | { return (eta1 - eta2); } |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | inline static double deltaPhi(const Jet& jet1, |
---|
48 | const Jet& jet2) |
---|
49 | { |
---|
50 | return fixedPhi(jet1.phi() - jet2.phi()); |
---|
51 | } |
---|
52 | inline static double deltaPhi(const Jet* jet1, |
---|
53 | const Jet* jet2) |
---|
54 | { |
---|
55 | return fixedPhi(jet1->phi()-jet2->phi()); |
---|
56 | } |
---|
57 | inline static double deltaPhi(const double phi1, |
---|
58 | const double phi2) |
---|
59 | { |
---|
60 | return fixedPhi( phi1 - phi2 ); |
---|
61 | } |
---|
62 | |
---|
63 | |
---|
64 | // distance in (eta,phi) |
---|
65 | inline static double deltaR(const Jet& jet1, |
---|
66 | const Jet& jet2) |
---|
67 | { return sqrt( deltaEta(jet1,jet2) * deltaEta(jet1,jet2) + |
---|
68 | deltaPhi(jet1,jet2) * deltaPhi(jet1,jet2) ); } |
---|
69 | inline static double deltaR(const Jet* jet1, |
---|
70 | const Jet* jet2) |
---|
71 | { return sqrt( deltaEta(jet1,jet2) * deltaEta(jet1,jet2) + |
---|
72 | deltaPhi(jet1,jet2) * deltaPhi(jet1,jet2) ); } |
---|
73 | |
---|
74 | inline static double deltaR(const double eta1, const double phi1, |
---|
75 | const double eta2, const double phi2 ) |
---|
76 | { return sqrt( deltaEta(eta1,eta2) * deltaEta(eta1,eta2) + |
---|
77 | deltaPhi(phi1,phi2) * deltaPhi(phi1,phi2) ); } |
---|
78 | |
---|
79 | |
---|
80 | |
---|
81 | // phi convention enforcement |
---|
82 | inline static double fixedPhi(double aPhi) |
---|
83 | { |
---|
84 | while ( aPhi < -M_PI ) |
---|
85 | { |
---|
86 | aPhi += 2.*M_PI; |
---|
87 | } |
---|
88 | while ( aPhi > M_PI ) |
---|
89 | { |
---|
90 | aPhi -= 2.*M_PI; |
---|
91 | } |
---|
92 | return aPhi; |
---|
93 | } |
---|
94 | |
---|
95 | }; |
---|
96 | |
---|
97 | } // namespace atlas |
---|
98 | |
---|
99 | FASTJET_END_NAMESPACE |
---|
100 | #endif |
---|