|
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
|
#include "EventFilter/GctRawToDigi/src/GctUnpackCollections.h"
GctUnpackCollections::GctUnpackCollections(edm::Event& event)
: m_event(event),
m_gctFibres(new L1GctFibreCollection()), // GCT input collections
m_rctEm(new L1CaloEmCollection()),
m_rctCalo(new L1CaloRegionCollection()),
m_gctInternEm(new L1GctInternEmCandCollection()), // GCT internal collections
m_gctInternJets(new L1GctInternJetDataCollection()),
m_gctInternEtSums(new L1GctInternEtSumCollection()),
m_gctInternHFData(new L1GctInternHFDataCollection()),
m_gctInternHtMiss(new L1GctInternHtMissCollection()),
m_gctIsoEm(new L1GctEmCandCollection()), // GCT output collections
m_gctNonIsoEm(new L1GctEmCandCollection()),
m_gctCenJets(new L1GctJetCandCollection()),
m_gctForJets(new L1GctJetCandCollection()),
m_gctTauJets(new L1GctJetCandCollection()),
m_gctHfBitCounts(new L1GctHFBitCountsCollection()),
m_gctHfRingEtSums(new L1GctHFRingEtSumsCollection()),
m_gctEtTot(new L1GctEtTotalCollection()),
m_gctEtHad(new L1GctEtHadCollection()),
m_gctEtMiss(new L1GctEtMissCollection()),
m_gctHtMiss(new L1GctHtMissCollection()),
m_gctJetCounts(new L1GctJetCountsCollection()), // Deprecated (empty collection still needed by GT)
m_errors(new L1TriggerErrorCollection()) // Misc
{
m_gctIsoEm->reserve(4);
m_gctCenJets->reserve(4);
m_gctForJets->reserve(4);
m_gctTauJets->reserve(4);
// ** DON'T RESERVE SPACE IN VECTORS FOR DEBUG UNPACK ITEMS! **
}
GctUnpackCollections::~GctUnpackCollections() {
// GCT input collections
m_event.put(std::move(m_gctFibres));
m_event.put(std::move(m_rctEm));
m_event.put(std::move(m_rctCalo));
// GCT internal collections
m_event.put(std::move(m_gctInternEm));
m_event.put(std::move(m_gctInternJets));
m_event.put(std::move(m_gctInternEtSums));
m_event.put(std::move(m_gctInternHFData));
m_event.put(std::move(m_gctInternHtMiss));
// GCT output collections
m_event.put(std::move(m_gctIsoEm), "isoEm");
m_event.put(std::move(m_gctNonIsoEm), "nonIsoEm");
m_event.put(std::move(m_gctCenJets), "cenJets");
m_event.put(std::move(m_gctForJets), "forJets");
m_event.put(std::move(m_gctTauJets), "tauJets");
m_event.put(std::move(m_gctHfBitCounts));
m_event.put(std::move(m_gctHfRingEtSums));
m_event.put(std::move(m_gctEtTot));
m_event.put(std::move(m_gctEtHad));
m_event.put(std::move(m_gctEtMiss));
m_event.put(std::move(m_gctHtMiss));
m_event.put(std::move(m_gctJetCounts)); // Deprecated (empty collection still needed by GT)
// Misc
m_event.put(std::move(m_errors));
}
std::ostream& operator<<(std::ostream& os, const GctUnpackCollections& rhs) {
// GCT input collections
os << "Read " << rhs.gctFibres()->size() << " GCT raw fibre data\n"
<< "Read " << rhs.rctEm()->size() << " RCT EM candidates\n"
<< "Read " << rhs.rctCalo()->size()
<< " RCT Calo Regions\n"
// GCT internal collections
<< "Read " << rhs.gctInternEm()->size() << " GCT intermediate EM candidates\n"
<< "Read " << rhs.gctInternJets()->size() << " GCT intermediate jet candidates\n"
<< "Read " << rhs.gctInternEtSums()->size() << " GCT intermediate et sums\n"
<< "Read " << rhs.gctInternHFData()->size() << " GCT intermediate HF data\n"
<< "Read " << rhs.gctInternHtMiss()->size()
<< " GCT intermediate Missing Ht\n"
// GCT output collections
<< "Read " << rhs.gctIsoEm()->size() << " GCT iso EM candidates\n"
<< "Read " << rhs.gctNonIsoEm()->size() << " GCT non-iso EM candidates\n"
<< "Read " << rhs.gctCenJets()->size() << " GCT central jet candidates\n"
<< "Read " << rhs.gctForJets()->size() << " GCT forward jet candidates\n"
<< "Read " << rhs.gctTauJets()->size() << " GCT tau jet candidates\n"
<< "Read " << rhs.gctHfBitCounts()->size() << " GCT HF ring bit counts\n"
<< "Read " << rhs.gctHfRingEtSums()->size() << " GCT HF ring et sums\n"
<< "Read " << rhs.gctEtTot()->size() << " GCT total et\n"
<< "Read " << rhs.gctEtHad()->size() << " GCT ht\n"
<< "Read " << rhs.gctEtMiss()->size() << " GCT met\n"
<< "Read " << rhs.gctHtMiss()->size() << " GCT mht";
// Any point in putting in an m_errors()->size()? Not sure.
return os;
}
|