|
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
|
#include "DataFormats/L1TParticleFlow/interface/HPSPFTau.h"
#include "FWCore/Utilities/interface/Exception.h"
// default constructor
l1t::HPSPFTau::HPSPFTau()
: tauType_(kUndefined),
sumChargedIso_(0.),
sumNeutralIso_(0.),
sumCombinedIso_(0.),
sumChargedIsoPileup_(0.),
passTightIso_(false),
passMediumIso_(false),
passLooseIso_(false),
passVLooseIso_(false),
passTightRelIso_(false),
passMediumRelIso_(false),
passLooseRelIso_(false),
passVLooseRelIso_(false) {}
// destructor
l1t::HPSPFTau::~HPSPFTau() {}
// print to stream
ostream& operator<<(ostream& os, const l1t::HPSPFTau& l1PFTau) {
os << "pT = " << l1PFTau.pt() << ", eta = " << l1PFTau.eta() << ", phi = " << l1PFTau.phi()
<< " (type = " << l1PFTau.tauType() << ")" << std::endl;
os << "lead. ChargedPFCand:" << std::endl;
if (l1PFTau.leadChargedPFCand().isNonnull()) {
printPFCand(os, *l1PFTau.leadChargedPFCand(), l1PFTau.primaryVertex());
} else {
os << " N/A" << std::endl;
}
os << "seed:";
if (l1PFTau.isChargedPFCandSeeded()) {
os << " chargedPFCand";
} else if (l1PFTau.isJetSeeded()) {
os << " CaloJet";
} else {
cms::Exception ex("InconsistentTau");
ex.addContext("Calling HPSPFTau::operator <<");
ex.addAdditionalInfo("This tau is not seed by either a chargedPFCand or a PFJet!");
throw ex;
}
os << std::endl;
os << "signalPFCands:" << std::endl;
for (const auto& l1PFCand : l1PFTau.signalAllL1PFCandidates()) {
printPFCand(os, *l1PFCand, l1PFTau.primaryVertex());
}
os << "stripPFCands:" << std::endl;
for (const auto& l1PFCand : l1PFTau.stripAllL1PFCandidates()) {
printPFCand(os, *l1PFCand, l1PFTau.primaryVertex());
}
os << "strip pT = " << l1PFTau.stripP4().pt() << std::endl;
os << "isolationPFCands:" << std::endl;
for (const auto& l1PFCand : l1PFTau.isoAllL1PFCandidates()) {
printPFCand(os, *l1PFCand, l1PFTau.primaryVertex());
}
os << "isolation pT-sum: charged = " << l1PFTau.sumChargedIso() << ", neutral = " << l1PFTau.sumNeutralIso()
<< " (charged from pileup = " << l1PFTau.sumChargedIsoPileup() << ")" << std::endl;
return os;
}
void printPFCand(ostream& os, const l1t::PFCandidate& l1PFCand, const l1t::VertexWordRef& primaryVertex) {
float primaryVertexZ = (primaryVertex.isNonnull()) ? primaryVertex->z0() : 0.;
printPFCand(os, l1PFCand, primaryVertexZ);
}
void printPFCand(ostream& os, const l1t::PFCandidate& l1PFCand, float primaryVertexZ) {
std::string typeString;
if (l1PFCand.id() == l1t::PFCandidate::ChargedHadron)
typeString = "PFChargedHadron";
else if (l1PFCand.id() == l1t::PFCandidate::Electron)
typeString = "PFElectron";
else if (l1PFCand.id() == l1t::PFCandidate::NeutralHadron)
typeString = "PFNeutralHadron";
else if (l1PFCand.id() == l1t::PFCandidate::Photon)
typeString = "PFPhoton";
else if (l1PFCand.id() == l1t::PFCandidate::Muon)
typeString = "PFMuon";
else
typeString = "N/A";
os << " " << typeString << " with pT = " << l1PFCand.pt() << ", eta = " << l1PFCand.eta()
<< ", phi = " << l1PFCand.phi() << ","
<< " mass = " << l1PFCand.mass() << ", charge = " << l1PFCand.charge();
if (l1PFCand.charge() != 0 && primaryVertexZ != 0.) {
os << " (dz = " << std::fabs(l1PFCand.pfTrack()->vertex().z() - primaryVertexZ) << ")";
}
os << std::endl;
}
|