|
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
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
|
#ifndef _DATAFORMATS_METRECO_HCALNOISEHPD_H_
#define _DATAFORMATS_METRECO_HCALNOISEHPD_H_
//
// HcalNoiseHPD.h
//
// description: Container class of HPD information to study anomalous noise in the HCAL.
// The information for HcalNoiseHPD's are filled in RecoMET/METProducers/HcalNoiseInfoProducer,
// but the idnumber is managed by DataFormats/METReco/HcalNoiseRBXArray.
// Provides relevant digi, rechit, and calotower information.
//
// author: J.P. Chou, Brown
//
#include <algorithm>
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/HcalDigi/interface/HBHEDataFrame.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "DataFormats/HcalRecHit/interface/HBHERecHitAuxSetter.h"
#include "DataFormats/HcalRecHit/interface/CaloRecHitAuxSetter.h"
#include "DataFormats/Common/interface/RefVector.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/Common/interface/RefProd.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
#include "DataFormats/METReco/interface/HcalCaloFlagLabels.h"
namespace reco {
//
// forward declarations
//
class HcalNoiseHPD;
class RefHBHERecHitEnergyComparison;
//
// typdefs
//
typedef std::vector<HcalNoiseHPD> HcalNoiseHPDCollection;
//
// RefHBHERecHitEnergyComparison is a class functor to compare energies between Ref<HBHERecHitCollection>
//
class RefHBHERecHitEnergyComparison {
public:
bool operator()(const edm::Ref<HBHERecHitCollection>& x, const edm::Ref<HBHERecHitCollection>& y) const {
return x->energy() > y->energy();
}
};
//
// class definition
//
class HcalNoiseHPD {
friend class HcalNoiseInfoProducer; // allows this class the fill the HPDs with info
friend class HcalNoiseRBXArray; // allows this class to manage the idnumber
public:
// constructor
HcalNoiseHPD();
// destructor
virtual ~HcalNoiseHPD();
//
// Detector ID accessors
//
// unique integer specifier for the hpd [0,NUM_HPDS-1]
// correlates roughly with the detector phi slice
int idnumber(void) const;
//
// Digi accessors
//
// pedestal subtracted fC information for the highest energy pixel in the HPD by timeslice
const std::vector<float> bigCharge(void) const;
float bigChargeTotal(void) const;
float bigChargeHighest2TS(unsigned int firstts = 4) const;
float bigChargeHighest3TS(unsigned int firstts = 4) const;
// same as above but the integral over the 5 highest energy pixels in the HPD
const std::vector<float> big5Charge(void) const;
float big5ChargeTotal(void) const;
float big5ChargeHighest2TS(unsigned int firstts = 4) const;
float big5ChargeHighest3TS(unsigned int firstts = 4) const;
// total number of adc zeros
int totalZeros(void) const;
// largest number of adc zeros in a digi in the HPD
int maxZeros(void) const;
//
// RecHit accessors
//
// returns a reference to a vector of references to the rechits
const edm::RefVector<HBHERecHitCollection> recHits(void) const;
// integral of rechit energies in the HPD with E>threshold (default is 1.5 GeV)
float recHitEnergy(float threshold = 1.5) const;
float recHitEnergyFailR45(float threshold = 1.5) const;
// minimum and maximum time for rechits with E>threshold (default is 10.0 GeV)
float minRecHitTime(float threshold = 10.0) const;
float maxRecHitTime(float threshold = 10.0) const;
// number of rechits with E>threshold (default is 1.5 GeV)
int numRecHits(float threshold = 1.5) const;
int numRecHitsFailR45(float threshold = 1.5) const;
//
// CaloTower accessors
//
// returns a reference to a vector of references to the calotowers
const edm::RefVector<CaloTowerCollection> caloTowers(void) const;
// calotower properties integrated over the entire HPD
double caloTowerHadE(void) const;
double caloTowerEmE(void) const;
double caloTowerTotalE(void) const;
double caloTowerEmFraction(void) const;
private:
// unique id number specifying the HPD
int idnumber_;
// digi data members
int totalZeros_;
int maxZeros_;
std::vector<float> bigCharge_;
std::vector<float> big5Charge_;
// a vector of references to rechits
edm::RefVector<HBHERecHitCollection> rechits_;
// a transient set of rechits for sorting purposes
// at some point before storing, these get transfered to the RefVector rechits_
std::set<edm::Ref<HBHERecHitCollection>, RefHBHERecHitEnergyComparison> refrechitset_;
// a vector of references to calotowers
edm::RefVector<CaloTowerCollection> calotowers_;
};
} // namespace reco
#endif
|