#include <TopJetAnalyzer.h>
Definition at line 15 of file TopJetAnalyzer.h.
Definition at line 4 of file TopJetAnalyzer.cc.
References en_, eta_, TFileService::make(), mult_, phi_, and pt_.
10 mult_ = fs->
make<TH1F>(
"mult",
"multiplicity (jets)", 30, 0 , 30);
11 en_ = fs->
make<TH1F>(
"en" ,
"energy (jets)", 60, 0., 300.);
12 pt_ = fs->make<TH1F>(
"pt" ,
"pt (jets)", 60, 0., 300.);
13 eta_ = fs->make<TH1F>(
"eta" ,
"eta (jets)", 30, -3., 3.);
14 phi_ = fs->make<TH1F>(
"phi" ,
"phi (jets)", 40, -4., 4.);
T getParameter(std::string const &) const
T * make(const Args &...args) const
make new ROOT object
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
edm::EDGetTokenT< std::vector< pat::Jet > > inputToken_
TopJetAnalyzer::~TopJetAnalyzer |
( |
| ) |
|
Implements edm::EDAnalyzer.
Definition at line 22 of file TopJetAnalyzer.cc.
References gather_cfg::cout, en_, eta_, edm::Event::getByToken(), i, inputToken_, metsig::jet, fwrapper::jets, mult_, phi_, pt_, and verbose_.
29 mult_->Fill( jets->size() );
30 for(std::vector<pat::Jet>::const_iterator
jet=jets->begin();
jet!=jets->end(); ++
jet){
32 en_ ->Fill(
jet->energy() );
43 if( jets->begin()->isCaloJet() )
45 else if( jets->begin()->isPFJet() )
48 std::cout << std::setfill(
'=') << std::setw(lineWidth) <<
"\n" << std::setfill(
' ');
50 << std::setw(11) <<
"pt :"
51 << std::setw( 9) <<
"eta :"
52 << std::setw( 9) <<
"phi :"
53 << std::setw(11) <<
"TCHE :"
54 << std::setw(11) <<
"TCHP :"
55 << std::setw( 9) <<
"SSVHE :"
56 << std::setw( 9) <<
"SSVHP :";
57 if( jets->begin()->isCaloJet() ) {
59 << std::setw(10) <<
"n90Hits :"
60 << std::setw( 7) <<
"fHPD";
62 if( jets->begin()->isPFJet() ) {
64 << std::setw(8) <<
"nhf : "
65 << std::setw(8) <<
"cef : "
66 << std::setw(8) <<
"nef : "
67 << std::setw(6) <<
"nCh : "
68 << std::setw(6) <<
"nConst";
71 << std::setfill(
'-') << std::setw(lineWidth) <<
"\n" << std::setfill(
' ');
73 for(std::vector<pat::Jet>::const_iterator
jet=jets->begin();
jet!=jets->end(); ++
jet){
74 std::cout << std::setw(3) << i <<
" : " << std::setprecision(3) << std::fixed
75 << std::setw(8) <<
jet->pt() <<
" : "
76 << std::setw(6) <<
jet->eta() <<
" : "
77 << std::setw(6) <<
jet->phi() <<
" : "
78 << std::setw(8) <<
jet->bDiscriminator(
"trackCountingHighEffBJetTags") <<
" : "
79 << std::setw(8) <<
jet->bDiscriminator(
"trackCountingHighPurBJetTags") <<
" : "
80 << std::setw(6) <<
jet->bDiscriminator(
"simpleSecondaryVertexHighEffBJetTags") <<
" : "
81 << std::setw(6) <<
jet->bDiscriminator(
"simpleSecondaryVertexHighPurBJetTags") <<
" : ";
82 if(
jet->isCaloJet() ) {
83 std::cout << std::setw(5) <<
jet->emEnergyFraction() <<
" : "
84 << std::setw(7) <<
jet->jetID().n90Hits <<
" : "
85 << std::setw(6) <<
jet->jetID().fHPD;
87 if(
jet->isPFJet() ) {
88 std::cout << std::setw(5) <<
jet->chargedHadronEnergyFraction() <<
" : "
89 << std::setw(5) <<
jet->neutralHadronEnergyFraction() <<
" : "
90 << std::setw(5) <<
jet->chargedEmEnergyFraction() <<
" : "
91 << std::setw(5) <<
jet->neutralEmEnergyFraction() <<
" : "
92 << std::setw(3) <<
jet->chargedMultiplicity() <<
" : "
93 << std::setw(6) <<
jet->nConstituents();
98 std::cout << std::setfill(
'=') << std::setw(lineWidth) <<
"\n" << std::setfill(
' ');
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::EDGetTokenT< std::vector< pat::Jet > > inputToken_
void TopJetAnalyzer::beginJob |
( |
void |
| ) |
|
|
privatevirtual |
void TopJetAnalyzer::endJob |
( |
void |
| ) |
|
|
privatevirtual |
TH1F* TopJetAnalyzer::en_ |
|
private |
TH1F* TopJetAnalyzer::eta_ |
|
private |
TH1F* TopJetAnalyzer::mult_ |
|
private |
TH1F* TopJetAnalyzer::phi_ |
|
private |
TH1F* TopJetAnalyzer::pt_ |
|
private |
bool TopJetAnalyzer::verbose_ |
|
private |