|
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
|
// CaloJet.cc
// Fedor Ratnikov UMd
#include <sstream>
#include "FWCore/Utilities/interface/Exception.h"
#include "DataFormats/RecoCandidate/interface/RecoCaloTowerCandidate.h"
//Own header file
#include "DataFormats/JetReco/interface/CaloJet.h"
using namespace reco;
CaloJet::CaloJet(const LorentzVector& fP4, const Point& fVertex, const Specific& fSpecific)
: Jet(fP4, fVertex), m_specific(fSpecific) {}
CaloJet::CaloJet(const LorentzVector& fP4,
const Point& fVertex,
const Specific& fSpecific,
const Jet::Constituents& fConstituents)
: Jet(fP4, fVertex, fConstituents), m_specific(fSpecific) {}
CaloJet::CaloJet(const LorentzVector& fP4, const Specific& fSpecific, const Jet::Constituents& fConstituents)
: Jet(fP4, Point(0, 0, 0), fConstituents), m_specific(fSpecific) {}
/// Physics Eta (use reference Z and jet kinematics only)
//float CaloJet::physicsEtaQuick (float fZVertex) const {
// return Jet::physicsEta (fZVertex, eta());
//}
/// Physics Eta (loop over constituents)
//float CaloJet::physicsEtaDetailed (float fZVertex) const {
// Jet::LorentzVector correctedMomentum;
// std::vector<const Candidate*> towers = getJetConstituentsQuick ();
// for (unsigned i = 0; i < towers.size(); ++i) {
// const Candidate* c = towers[i];
// double etaRef = Jet::physicsEta (fZVertex, c->eta());
// math::PtEtaPhiMLorentzVectorD p4 (c->p()/cosh(etaRef), etaRef, c->phi(), c->mass());
// correctedMomentum += p4;
// }
// return correctedMomentum.eta();
//}
/// Physics p4 (use jet Z and kinematics only)
//deprecated, with 3d vertex correction not clear anymore!
//CaloJet::LorentzVector CaloJet::physicsP4 (float fZVertex) const {
// double physicsEta = Jet::physicsEta (fZVertex, eta());
// math::PtEtaPhiMLorentzVectorD p4 (p()/cosh(physicsEta), physicsEta, phi(), mass());
// return CaloJet::LorentzVector (p4);
//}
CaloJet::LorentzVector CaloJet::physicsP4(const Particle::Point& vertex) const {
return Jet::physicsP4(vertex, *this, this->vertex());
}
CaloJet::LorentzVector CaloJet::detectorP4() const { return Jet::detectorP4(this->vertex(), *this); }
CaloTowerPtr CaloJet::getCaloConstituent(unsigned fIndex) const {
Constituent dau = daughterPtr(fIndex);
if (dau.isNonnull() && dau.isAvailable()) {
const CaloTower* towerCandidate = dynamic_cast<const CaloTower*>(dau.get());
if (towerCandidate) {
// return towerCandidate;
// 086 Ptr(ProductID const& productID, T const* item, key_type item_key) :
return edm::Ptr<CaloTower>(dau.id(), towerCandidate, dau.key());
} else {
throw cms::Exception("Invalid Constituent") << "CaloJet constituent is not of CaloTowere type";
}
}
else {
return CaloTowerPtr();
}
}
std::vector<CaloTowerPtr> CaloJet::getCaloConstituents() const {
std::vector<CaloTowerPtr> result;
for (unsigned i = 0; i < numberOfDaughters(); i++)
result.push_back(getCaloConstituent(i));
return result;
}
CaloJet* CaloJet::clone() const { return new CaloJet(*this); }
bool CaloJet::overlap(const Candidate&) const { return false; }
std::string CaloJet::print() const {
std::ostringstream out;
out << Jet::print() // generic jet info
<< " CaloJet specific:" << std::endl
<< " energy fractions em/had: " << emEnergyFraction() << '/' << energyFractionHadronic() << std::endl
<< " em energy in EB/EE/HF: " << emEnergyInEB() << '/' << emEnergyInEE() << '/' << emEnergyInHF()
<< std::endl
<< " had energy in HB/HO/HE/HF: " << hadEnergyInHB() << '/' << hadEnergyInHO() << '/' << hadEnergyInHE()
<< '/' << hadEnergyInHF() << std::endl
<< " constituent towers area: " << towersArea() << std::endl;
out << " Towers:" << std::endl;
std::vector<CaloTowerPtr> towers = getCaloConstituents();
for (unsigned i = 0; i < towers.size(); i++) {
if (towers[i].get()) {
out << " #" << i << " " << *(towers[i]) << std::endl;
} else {
out << " #" << i << " tower is not available in the event" << std::endl;
}
}
return out.str();
}
//----------------------------------------------------------
// here are methods extracting information from constituents
//----------------------------------------------------------
std::vector<CaloTowerDetId> CaloJet::getTowerIndices() const {
std::vector<CaloTowerDetId> result;
std::vector<CaloTowerPtr> towers = getCaloConstituents();
for (unsigned i = 0; i < towers.size(); ++i) {
result.push_back(towers[i]->id());
}
return result;
}
|