|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
#ifndef DataFormats_METReco_interface_PhiWedge_h
#define DataFormats_METReco_interface_PhiWedge_h
/*
[class]: PhiWedge
[authors]: R. Remington, The University of Florida
[description]: Simple class analogous to CaloTower but in z-direction. Stores basic information related to Hcal and Ecal rechits within constant 5-degree phi windows. The idea will be to match these reconstructed phi-wedges with csc tracks for BeamHalo identification.
[date]: October 15, 2009
*/
#include <vector>
#include "TMath.h"
namespace reco {
class PhiWedge {
public:
// Constructors
constexpr PhiWedge()
: energy_{0.f},
iphi_{0},
constituents_{0},
min_time_{0.f},
max_time_{0.f},
PlusZOriginConfidence_{0.f},
OverlappingCSCTracks_{0},
OverlappingCSCSegments_{0},
OverlappingCSCRecHits_{0},
OverlappingCSCHaloTriggers_{0} {}
constexpr PhiWedge(float E, int iphi, int constituents)
: energy_{E},
iphi_{iphi},
constituents_{constituents},
min_time_{0.f},
max_time_{0.f},
PlusZOriginConfidence_{0.f},
OverlappingCSCTracks_{0},
OverlappingCSCSegments_{0},
OverlappingCSCRecHits_{0},
OverlappingCSCHaloTriggers_{0} {}
constexpr PhiWedge(float E, int iphi, int constituents, float min_time, float max_time)
: energy_{E},
iphi_{iphi},
constituents_{constituents},
min_time_{min_time},
max_time_{max_time},
PlusZOriginConfidence_{0.f},
OverlappingCSCTracks_{0},
OverlappingCSCSegments_{0},
OverlappingCSCRecHits_{0},
OverlappingCSCHaloTriggers_{0} {}
constexpr PhiWedge(PhiWedge const &) = default;
constexpr PhiWedge(PhiWedge &&) = default;
constexpr PhiWedge &operator=(PhiWedge const &) = default;
constexpr PhiWedge &operator=(PhiWedge &&) = default;
// Destructor
constexpr ~PhiWedge() = default;
// Energy sum of all rechits above threshold in this 5-degree window
float Energy() const { return energy_; }
// Number of rechits above threshold in this 5-degree window
int NumberOfConstituents() const { return constituents_; }
// iPhi value of this 5-degree window
int iPhi() const { return iphi_; }
// Global phi lower bound of this 5-degree window (between 0 and 2Pi)
float PhiLow() const { return 2. * TMath::Pi() * (float)((iphi_ * 5) - (5.)); }
// Global phi upper bound of this 5-degree window (between 0 and 2Pi
float PhiHigh() const { return 2. * TMath::Pi() * (float)((iphi_ * 5)); }
// Get Min/Max Time
float MinTime() const { return min_time_; }
float MaxTime() const { return max_time_; }
// Get halo direction confidence based on time ordering of the rechits ( must be within range of -1 to + 1 )
// Direction is calculated by counting the number of pair-wise time-ascending rechits from -Z to +Z and then normalizing this count by number of pair-wise combinations
// If all pair-wise combinations are consistent with a common z-direction, then this value will be plus or minus 1 exactly. Otherwise it will be somewhere in between.
float ZDirectionConfidence() const { return (1. - PlusZOriginConfidence_) * 2. - 1.; }
float PlusZDirectionConfidence() const { return 1. - PlusZOriginConfidence_; }
float MinusZDirectionConfidence() const { return PlusZOriginConfidence_; }
// Get halo origin confidence based on time ordering of the rechits
float PlusZOriginConfidence() const { return PlusZOriginConfidence_; }
float MinusZOriginConfidence() const { return 1. - PlusZOriginConfidence_; }
// To be filled later or removed
int OverlappingCSCTracks() const { return OverlappingCSCTracks_; }
int OverlappingCSCSegments() const { return OverlappingCSCSegments_; }
int OverlappingCSCRecHits() const { return OverlappingCSCRecHits_; }
int OverlappingCSCHaloTriggers() const { return OverlappingCSCHaloTriggers_; }
// Setters
void SetOverlappingCSCTracks(int x) { OverlappingCSCTracks_ = x; }
void SetOverlappingCSCSegments(int x) { OverlappingCSCSegments_ = x; }
void SetOverlappingCSCRecHits(int x) { OverlappingCSCRecHits_ = x; }
void SetOverlappingCSCHaloTriggers(int x) { OverlappingCSCHaloTriggers_ = x; }
void SetMinMaxTime(float min, float max) {
min_time_ = min;
max_time_ = max;
}
void SetPlusZOriginConfidence(float x) { PlusZOriginConfidence_ = x; }
private:
float energy_;
int iphi_;
int constituents_;
float min_time_;
float max_time_;
float PlusZOriginConfidence_;
int OverlappingCSCTracks_;
int OverlappingCSCSegments_;
int OverlappingCSCRecHits_;
int OverlappingCSCHaloTriggers_;
};
using PhiWedgeCollection = std::vector<PhiWedge>;
} // namespace reco
#endif // DataFormats_METReco_interface_PhiWedge_h
|