#include <RecoBTag/JetVertexChecker/src/JetVertexChecker.cc>

Inheritance diagram for JetVertexChecker:

Public Member Functions
	JetVertexChecker (const edm::ParameterSet &)

	~JetVertexChecker ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Private Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::EDGetTokenT < reco::JetTracksAssociationCollection >	m_associator

edm::EDGetTokenT< reco::BeamSpot >	m_beamSpot

double	m_cutMinPt

double	m_cutMinPtRatio

bool	m_doFilter

double	m_maxChi2

int32_t	m_maxNjets

int32_t	m_maxNjetsOutput

double	m_maxTrackPt

bool	m_newMethod

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 47 of file JetVertexChecker.cc.

Constructor & Destructor Documentation

JetVertexChecker::JetVertexChecker ( const edm::ParameterSet & iConfig )

explicit

Definition at line 84 of file JetVertexChecker.cc.

References edm::ParameterSet::getParameter(), m_associator, m_beamSpot, m_cutMinPt, m_cutMinPtRatio, m_doFilter, m_maxChi2, m_maxNjets, m_maxNjetsOutput, m_maxTrackPt, and m_newMethod.

 {
    //now do what ever initialization is needed
   m_beamSpot                = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpot"));
   m_associator              = consumes<reco::JetTracksAssociationCollection>(iConfig.getParameter<edm::InputTag>("jetTracks"));
   m_doFilter                = iConfig.getParameter<bool>("doFilter");
   m_cutMinPt                = iConfig.getParameter<double>("minPt");
   m_cutMinPtRatio           = iConfig.getParameter<double>("minPtRatio");
   m_maxNjets                = iConfig.getParameter<int32_t>("maxNJetsToCheck");
   m_maxNjetsOutput          = iConfig.getParameter<int32_t>("maxNjetsOutput");
   m_newMethod                = iConfig.getParameter<bool>("newMethod");
   m_maxTrackPt           = iConfig.getParameter<double>("maxTrackPt");
   m_maxChi2           = iConfig.getParameter<double>("maxChi2");
   produces<std::vector<reco::CaloJet> >(); 
   produces<reco::VertexCollection >(); 
 }

JetVertexChecker::~JetVertexChecker ( )

Definition at line 102 of file JetVertexChecker.cc.

 {
  
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void JetVertexChecker::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 173 of file JetVertexChecker.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                              {
    edm::ParameterSetDescription desc;
    desc.add<edm::InputTag> ("beamSpot",edm::InputTag("hltOnlineBeamSpot"));
    desc.add<edm::InputTag> ("jetTracks",edm::InputTag("hltFastPVJetTracksAssociator"));
    desc.add<double> ("minPtRatio",0.1);
    desc.add<double> ("minPt",0.0);
    desc.add<bool> ("doFilter",false);
    desc.add<int> ("maxNJetsToCheck",2);
    desc.add<int> ("maxNjetsOutput",2);
    desc.add<double> ("maxChi2",20.0);
    desc.add<double> ("maxTrackPt",20.0);
    desc.add<bool> ("newMethod",false);      // <---- newMethod 
    descriptions.add("jetVertexChecker",desc);
 }

bool JetVertexChecker::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDFilter.

Definition at line 117 of file JetVertexChecker.cc.

References funct::abs(), SiPixelRawToDigiRegional_cfi::beamSpot, edm::RefVector< C, T, F >::begin(), reco::TrackBase::chi2(), alignCSCRings::e, edm::RefVector< C, T, F >::end(), edm::Event::getByToken(), i, m_associator, m_beamSpot, m_cutMinPt, m_cutMinPtRatio, m_doFilter, m_maxChi2, m_maxNjets, m_maxTrackPt, m_newMethod, min(), reco::TrackBase::momentum(), AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::pt(), edm::Event::put(), query::result, and testEve_cfg::tracks.

 {  
    using namespace edm;
    Handle<reco::JetTracksAssociationCollection> jetTracksAssociation;
    iEvent.getByToken(m_associator, jetTracksAssociation);
    std::auto_ptr<std::vector<reco::CaloJet> > pOut(new std::vector<reco::CaloJet> );
 
    bool result=true;
    int i = 0;
    float calopt=0;
    float trkpt=0;
    //limit to first two jets
    for(reco::JetTracksAssociationCollection::const_iterator it = jetTracksAssociation->begin();
        it != jetTracksAssociation->end() && i < m_maxNjets; it++, i++) {
      if(std::abs(it->first->eta()) < 2.4)
      {
       reco::TrackRefVector tracks = it->second;
       math::XYZVector jetMomentum = it->first->momentum();
       math::XYZVector trMomentum;
       for(reco::TrackRefVector::const_iterator itTrack = tracks.begin(); itTrack != tracks.end(); ++itTrack) 
       {
          const reco::Track& iTrack = **itTrack;
          if(m_newMethod && iTrack.chi2()>m_maxChi2) continue;
              trMomentum += iTrack.momentum();
          if(m_newMethod) trkpt += std::min(m_maxTrackPt,( iTrack.pt()));
          else trkpt += iTrack.pt();
       }
       calopt += jetMomentum.rho();
       if(trMomentum.rho()/jetMomentum.rho() < m_cutMinPtRatio || trMomentum.rho() < m_cutMinPt) 
       {
         pOut->push_back(* dynamic_cast<const reco::CaloJet *>(&(*it->first)));
       }
      }
     }
    iEvent.put(pOut);
 
    edm::Handle<reco::BeamSpot> beamSpot;
    iEvent.getByToken(m_beamSpot,beamSpot);
  
    reco::Vertex::Error e;
    e(0, 0) = 0.0015 * 0.0015;
    e(1, 1) = 0.0015 * 0.0015;
    e(2, 2) = 1.5 * 1.5;
    reco::Vertex::Point p(beamSpot->x0(), beamSpot->y0(), beamSpot->z0());
    reco::Vertex thePV(p, e, 0, 0, 0);
    std::auto_ptr<reco::VertexCollection> pOut2(new reco::VertexCollection);
    pOut2->push_back(thePV);
    iEvent.put(pOut2);
 
    if(m_doFilter) return result;
    else 
    return true;
 }