|
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
|
#include "DataFormats/TauReco/interface/PFRecoTauChargedHadron.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
namespace reco {
PFRecoTauChargedHadron::PFRecoTauChargedHadron() : CompositePtrCandidate(), algo_(kUndefined) {}
PFRecoTauChargedHadron::PFRecoTauChargedHadron(PFRecoTauChargedHadronAlgorithm algo, Charge q)
: CompositePtrCandidate(), algo_(algo) {
if (q > 0.)
this->setPdgId(+211);
else if (q < 0.)
this->setPdgId(-211);
}
PFRecoTauChargedHadron::PFRecoTauChargedHadron(Charge q,
const LorentzVector& p4,
const Point& vtx,
int status,
bool integerCharge,
PFRecoTauChargedHadronAlgorithm algo)
: CompositePtrCandidate(q, p4, vtx, 211, status, integerCharge), algo_(algo) {
if (q > 0.)
this->setPdgId(+211);
else if (q < 0.)
this->setPdgId(-211);
}
PFRecoTauChargedHadron::PFRecoTauChargedHadron(const Candidate& c, PFRecoTauChargedHadronAlgorithm algo)
: CompositePtrCandidate(c), algo_(algo) {
if (c.charge() > 0.)
this->setPdgId(+211);
else if (c.charge() < 0.)
this->setPdgId(-211);
}
PFRecoTauChargedHadron::~PFRecoTauChargedHadron() {}
const CandidatePtr& PFRecoTauChargedHadron::getChargedPFCandidate() const { return chargedPFCandidate_; }
const PFRecoTauChargedHadron::TrackPtr& PFRecoTauChargedHadron::getTrack() const { return track_; }
const CandidatePtr& PFRecoTauChargedHadron::getLostTrackCandidate() const { return lostTrackCandidate_; }
const std::vector<CandidatePtr>& PFRecoTauChargedHadron::getNeutralPFCandidates() const {
return neutralPFCandidates_;
}
const math::XYZPointF& PFRecoTauChargedHadron::positionAtECALEntrance() const { return positionAtECALEntrance_; }
PFRecoTauChargedHadron::PFRecoTauChargedHadronAlgorithm PFRecoTauChargedHadron::algo() const { return algo_; }
bool PFRecoTauChargedHadron::algoIs(PFRecoTauChargedHadron::PFRecoTauChargedHadronAlgorithm algo) const {
return (algo_ == algo);
}
void PFRecoTauChargedHadron::print(std::ostream& stream) const {
stream << " Pt = " << this->pt() << ", eta = " << this->eta() << ", phi = " << this->phi()
<< " (mass = " << this->mass() << ")" << std::endl;
stream << " charge = " << this->charge() << " (pdgId = " << this->pdgId() << ")" << std::endl;
stream << "charged PFCandidate";
if (chargedPFCandidate_.isNonnull()) {
stream << " (" << chargedPFCandidate_.id() << ":" << chargedPFCandidate_.key() << "):"
<< " Pt = " << chargedPFCandidate_->pt() << ", eta = " << chargedPFCandidate_->eta()
<< ", phi = " << chargedPFCandidate_->phi() << " (pdgId = " << chargedPFCandidate_->pdgId() << ")"
<< std::endl;
} else {
stream << ": N/A" << std::endl;
}
stream << "reco::Track: ";
if (track_.isNonnull()) {
stream << "Pt = " << track_->pt() << " +/- " << track_->ptError() << ", eta = " << track_->eta()
<< ", phi = " << track_->phi() << std::endl;
} else if (lostTrackCandidate_.isNonnull()) {
stream << "(lostTrackCandidate: " << lostTrackCandidate_.id() << ":" << lostTrackCandidate_.key() << "):"
<< " Pt = " << lostTrackCandidate_->pt() << ", eta = " << lostTrackCandidate_->eta()
<< ", phi = " << lostTrackCandidate_->phi() << std::endl;
} else {
stream << "N/A" << std::endl;
}
stream << "neutral PFCandidates:";
if (!neutralPFCandidates_.empty()) {
stream << std::endl;
int idx = 0;
for (std::vector<CandidatePtr>::const_iterator neutralPFCandidate = neutralPFCandidates_.begin();
neutralPFCandidate != neutralPFCandidates_.end();
++neutralPFCandidate) {
stream << " #" << idx << " (" << neutralPFCandidate->id() << ":" << neutralPFCandidate->key() << "):"
<< " Pt = " << (*neutralPFCandidate)->pt() << ", eta = " << (*neutralPFCandidate)->eta()
<< ", phi = " << (*neutralPFCandidate)->phi() << " (pdgId = " << (*neutralPFCandidate)->pdgId() << ")"
<< std::endl;
++idx;
}
} else {
stream << " ";
stream << "N/A" << std::endl;
}
stream << "position@ECAL entrance: x = " << this->positionAtECALEntrance().x()
<< ", y = " << this->positionAtECALEntrance().y() << ", z = " << this->positionAtECALEntrance().z()
<< " (eta = " << this->positionAtECALEntrance().eta() << ", phi = " << this->positionAtECALEntrance().phi()
<< ")" << std::endl;
std::string algo_string = "undefined";
if (algo_ == kChargedPFCandidate)
algo_string = "chargedPFCandidate";
else if (algo_ == kTrack)
algo_string = "Track";
else if (algo_ == kPFNeutralHadron)
algo_string = "PFNeutralHadron";
stream << "algo = " << algo_string << std::endl;
}
std::ostream& operator<<(std::ostream& stream, const reco::PFRecoTauChargedHadron& c) {
c.print(stream);
return stream;
}
} // namespace reco
|