|
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
|
// JPTJet.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/JPTJet.h"
using namespace reco;
JPTJet::JPTJet(const LorentzVector& fP4,
const Point& fVertex,
const Specific& fSpecific,
const Jet::Constituents& fConstituents)
: Jet(fP4, fVertex), mspecific(fSpecific) {}
JPTJet::JPTJet(const LorentzVector& fP4, const Specific& fSpecific, const Jet::Constituents& fConstituents)
: Jet(fP4, Point(0, 0, 0)), mspecific(fSpecific) {}
JPTJet* JPTJet::clone() const { return new JPTJet(*this); }
bool JPTJet::overlap(const Candidate&) const { return false; }
void JPTJet::printJet() const {
std::cout << " Raw Calo jet " << getCaloJetRef()->et() << " " << getCaloJetRef()->eta() << " "
<< getCaloJetRef()->phi() << " JPTJet specific:" << std::endl
<< " charged multiplicity: " << chargedMultiplicity() << std::endl;
std::cout << " JPTCandidate constituents:" << std::endl;
std::cout << " Number of pions: " << getPionsInVertexInCalo().size() + getPionsInVertexOutCalo().size() << std::endl;
std::cout << " Number of muons: " << getMuonsInVertexInCalo().size() + getMuonsInVertexOutCalo().size() << std::endl;
std::cout << " Number of Electrons: " << getElecsInVertexInCalo().size() + getElecsInVertexOutCalo().size()
<< std::endl;
}
std::string JPTJet::print() const {
std::ostringstream out;
out << Jet::print() // generic jet info
<< " JPTJet specific:" << std::endl
<< " charged: " << chargedMultiplicity() << std::endl;
out << " JPTCandidate constituents:" << std::endl;
out << " Number of pions: " << getPionsInVertexInCalo().size() + getPionsInVertexOutCalo().size() << std::endl;
out << " Number of muons: " << getMuonsInVertexInCalo().size() + getMuonsInVertexOutCalo().size() << std::endl;
out << " Number of Electrons: " << getElecsInVertexInCalo().size() + getElecsInVertexOutCalo().size() << std::endl;
return out.str();
}
|