TopJetAnalyzer Class Reference

#include <TopJetAnalyzer.h>

Inheritance diagram for TopJetAnalyzer:

Public Member Functions
	TopJetAnalyzer (const edm::ParameterSet &)
	~TopJetAnalyzer ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()

Private Attributes
TH1F *	en_
TH1F *	eta_
edm::InputTag	input_
TH1F *	mult_
TH1F *	phi_
TH1F *	pt_
bool	verbose_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 13 of file TopJetAnalyzer.h.

Constructor & Destructor Documentation

TopJetAnalyzer::TopJetAnalyzer ( const edm::ParameterSet &  cfg )

explicit

Definition at line 5 of file TopJetAnalyzer.cc.

References en_, eta_, TFileDirectory::make(), mult_, phi_, and pt_.

                                                        :
  input_  (cfg.getParameter<edm::InputTag>("input"  )),
  verbose_(cfg.getParameter<bool>         ("verbose"))
{
  edm::Service<TFileService> fs;
  
  mult_ = fs->make<TH1F>("mult", "multiplicity (jets)", 30,  0 ,   30);
  en_   = fs->make<TH1F>("en"  , "energy (jets)",       60,  0., 300.);
  pt_   = fs->make<TH1F>("pt"  , "pt (jets)",           60,  0., 300.);
  eta_  = fs->make<TH1F>("eta" , "eta (jets)",          30, -3.,   3.);
  phi_  = fs->make<TH1F>("phi" , "phi (jets)",          40, -4.,   4.);
}

TopJetAnalyzer::~TopJetAnalyzer

(

)

Definition at line 18 of file TopJetAnalyzer.cc.

{
}

Member Function Documentation

void TopJetAnalyzer::analyze ( const edm::Event &  evt, const edm::EventSetup &  setup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 23 of file TopJetAnalyzer.cc.

References gather_cfg::cout, en_, eta_, edm::Event::getByLabel(), i, input_, metsig::jet, fwrapper::jets, mult_, phi_, pt_, and verbose_.

{
  edm::Handle<std::vector<pat::Jet> > jets;
  evt.getByLabel(input_, jets); 

  // fill histograms
  
  mult_->Fill( jets->size() );
  for(std::vector<pat::Jet>::const_iterator jet=jets->begin(); jet!=jets->end(); ++jet){
    pt_ ->Fill( jet->pt()     );
    en_ ->Fill( jet->energy() );
    eta_->Fill( jet->eta()    );
    phi_->Fill( jet->phi()    );
  }

  // produce printout if desired

  if( jets->size()<1 || !verbose_ )
    return;

  int lineWidth = 75;
  if( jets->begin()->isCaloJet() )
    lineWidth = 100;
  else if( jets->begin()->isPFJet() )
    lineWidth = 120;

  std::cout << std::setfill('=') << std::setw(lineWidth) << "\n" << std::setfill(' ');
  std::cout << std::setw( 5) << "jet :"
            << std::setw(11) << "pt :"
            << std::setw( 9) << "eta :"
            << std::setw( 9) << "phi :"
        << std::setw(11) << "TCHE :"
            << std::setw(11) << "TCHP :"
            << std::setw( 9) << "SSVHE :"
            << std::setw( 9) << "SSVHP :";
  if( jets->begin()->isCaloJet() ) {
    std::cout << std::setw( 8) << "emf :"
          << std::setw(10) << "n90Hits :"
          << std::setw( 7) << "fHPD";
  }
  if( jets->begin()->isPFJet() ) {
    std::cout << std::setw(9) << "chf : "
          << std::setw(8) << "nhf : "
          << std::setw(8) << "cef : "
          << std::setw(8) << "nef : "
          << std::setw(6) << "nCh : "
          << std::setw(6) << "nConst";
  }
  std::cout << std::endl
        << std::setfill('-') << std::setw(lineWidth) << "\n" << std::setfill(' ');
  unsigned i=0;
  for(std::vector<pat::Jet>::const_iterator jet=jets->begin(); jet!=jets->end(); ++jet){
    std::cout << std::setw(3) << i << " : " << std::setprecision(3) << std::fixed
          << std::setw(8) << jet->pt() << " : "
          << std::setw(6) << jet->eta() << " : "
          << std::setw(6) << jet->phi() << " : "
          << std::setw(8) << jet->bDiscriminator("trackCountingHighEffBJetTags") << " : "
          << std::setw(8) << jet->bDiscriminator("trackCountingHighPurBJetTags") << " : "
          << std::setw(6) << jet->bDiscriminator("simpleSecondaryVertexHighEffBJetTags") << " : "
          << std::setw(6) << jet->bDiscriminator("simpleSecondaryVertexHighPurBJetTags") << " : ";
    if( jet->isCaloJet() ) {
      std::cout << std::setw(5) << jet->emEnergyFraction() << " : "
        << std::setw(7) << jet->jetID().n90Hits << " : "
        << std::setw(6) << jet->jetID().fHPD;
    }
    if( jet->isPFJet() ) {
      std::cout << std::setw(5) << jet->chargedHadronEnergyFraction() << " : "
        << std::setw(5) << jet->neutralHadronEnergyFraction() << " : "
        << std::setw(5) << jet->chargedEmEnergyFraction() << " : "
        << std::setw(5) << jet->neutralEmEnergyFraction() << " : "
        << std::setw(3) << jet->chargedMultiplicity() << " : "
        << std::setw(6) << jet->nConstituents();
    }
    std::cout << std::endl;
    i++;
  }
  std::cout << std::setfill('=') << std::setw(lineWidth) << "\n" << std::setfill(' ');
}

void TopJetAnalyzer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 102 of file TopJetAnalyzer.cc.

{
}

void TopJetAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 106 of file TopJetAnalyzer.cc.

{
}

Member Data Documentation

TH1F* TopJetAnalyzer::en_

private

Definition at line 30 of file TopJetAnalyzer.h.

Referenced by analyze(), and TopJetAnalyzer().

TH1F* TopJetAnalyzer::eta_

private

Definition at line 32 of file TopJetAnalyzer.h.

Referenced by analyze(), and TopJetAnalyzer().

edm::InputTag TopJetAnalyzer::input_

private

Definition at line 26 of file TopJetAnalyzer.h.

Referenced by analyze().

TH1F* TopJetAnalyzer::mult_

private

Definition at line 29 of file TopJetAnalyzer.h.

Referenced by analyze(), and TopJetAnalyzer().

TH1F* TopJetAnalyzer::phi_

private

Definition at line 33 of file TopJetAnalyzer.h.

Referenced by analyze(), and TopJetAnalyzer().

TH1F* TopJetAnalyzer::pt_

private

Definition at line 31 of file TopJetAnalyzer.h.

Referenced by analyze(), and TopJetAnalyzer().

bool TopJetAnalyzer::verbose_

private

Definition at line 27 of file TopJetAnalyzer.h.

Referenced by analyze().