Fork me on GitHub

source: git/external/fastjet/plugins/ATLASCone/JetDistances.hh@ 282f591

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 282f591 was d7d2da3, checked in by pavel <pavel@…>, 12 years ago

move branches/ModularDelphes to trunk

  • Property mode set to 100644
File size: 2.5 KB
RevLine 
[d7d2da3]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
28FASTJET_BEGIN_NAMESPACE
29
30namespace atlas {
31
32struct 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
99FASTJET_END_NAMESPACE
100#endif
Note: See TracBrowser for help on using the repository browser.