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

Inheritance diagram for JetVertexChecker:

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

	~JetVertexChecker () override

Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

	~EDFilterBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

const edm::EDGetTokenT< reco::BeamSpot >	m_beamSpot

const double	m_cutMinPt

const double	m_cutMinPtRatio

const bool	m_doFilter

const double	m_maxChi2

const double	m_maxETA

const int32_t	m_maxNjets

const int32_t	m_maxNjetsOutput

const double	m_maxTrackPt

const bool	m_newMethod

const double	m_pvErr_x

const double	m_pvErr_y

const double	m_pvErr_z

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 88 of file JetVertexChecker.cc.

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

JetVertexChecker::~JetVertexChecker ( )

override

Definition at line 110 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 181 of file JetVertexChecker.cc.

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

                                                                              {
    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 
    desc.add<double> ("maxETA",         2.4   );
    desc.add<double> ("pvErr_x",        0.0015);
    desc.add<double> ("pvErr_y",        0.0015);
    desc.add<double> ("pvErr_z",        1.5   );
    descriptions.add("jetVertexChecker",desc);
 }

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

overrideprivate

Definition at line 125 of file JetVertexChecker.cc.

References funct::abs(), ecalDrivenElectronSeedsParameters_cff::beamSpot, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::begin(), edm::RefVector< C, T, F >::begin(), reco::TrackBase::chi2(), MillePedeFileConverter_cfg::e, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::end(), edm::RefVector< C, T, F >::end(), stringResolutionProvider_cfi::et, edm::Event::getByToken(), mps_fire::i, m_associator, m_beamSpot, m_cutMinPt, m_cutMinPtRatio, m_doFilter, m_maxChi2, m_maxETA, m_maxNjets, m_maxTrackPt, m_newMethod, m_pvErr_x, m_pvErr_y, m_pvErr_z, min(), reco::TrackBase::momentum(), eostools::move(), AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::pt(), edm::Event::put(), mps_fire::result, l1t::tracks, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

 {  
    using namespace edm;
    Handle<reco::JetTracksAssociationCollection> jetTracksAssociation;
    iEvent.getByToken(m_associator, jetTracksAssociation);
    auto pOut = std::make_unique<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(), et = jetTracksAssociation->end();
        it != et && i < m_maxNjets; it++, i++) {
      if(std::abs(it->first->eta()) < m_maxETA)
      {
       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(std::move(pOut));
 
    edm::Handle<reco::BeamSpot> beamSpot;
    iEvent.getByToken(m_beamSpot,beamSpot);
  
    reco::Vertex::Error e;
    e(0, 0) = m_pvErr_x * m_pvErr_x;
    e(1, 1) = m_pvErr_y * m_pvErr_y;
    e(2, 2) = m_pvErr_z * m_pvErr_z;
    reco::Vertex::Point p(beamSpot->x0(), beamSpot->y0(), beamSpot->z0());
    reco::Vertex thePV(p, e, 0, 0, 0);
    auto pOut2 = std::make_unique<reco::VertexCollection>();
    pOut2->push_back(thePV);
    iEvent.put(std::move(pOut2));
 
    if(m_doFilter) return result;
    else 
    return true;
 }

Member Data Documentation

const edm::EDGetTokenT<reco::JetTracksAssociationCollection> JetVertexChecker::m_associator

private

Definition at line 58 of file JetVertexChecker.cc.

Referenced by filter().

const edm::EDGetTokenT<reco::BeamSpot> JetVertexChecker::m_beamSpot

private

Definition at line 59 of file JetVertexChecker.cc.

Referenced by filter().

const double JetVertexChecker::m_cutMinPt

private

Definition at line 61 of file JetVertexChecker.cc.

Referenced by filter().

const double JetVertexChecker::m_cutMinPtRatio

private

Definition at line 62 of file JetVertexChecker.cc.

Referenced by filter().

const bool JetVertexChecker::m_doFilter

private

Definition at line 60 of file JetVertexChecker.cc.

Referenced by filter().

const double JetVertexChecker::m_maxChi2

private

Definition at line 64 of file JetVertexChecker.cc.

Referenced by filter().

const double JetVertexChecker::m_maxETA

private

Definition at line 70 of file JetVertexChecker.cc.

Referenced by filter().

const int32_t JetVertexChecker::m_maxNjets

private

Definition at line 65 of file JetVertexChecker.cc.

Referenced by filter().

const int32_t JetVertexChecker::m_maxNjetsOutput

private

Definition at line 66 of file JetVertexChecker.cc.

const double JetVertexChecker::m_maxTrackPt

private

Definition at line 63 of file JetVertexChecker.cc.

Referenced by filter().

const bool JetVertexChecker::m_newMethod

private

Definition at line 68 of file JetVertexChecker.cc.

Referenced by filter().

const double JetVertexChecker::m_pvErr_x

private

Definition at line 71 of file JetVertexChecker.cc.

Referenced by filter().

const double JetVertexChecker::m_pvErr_y

private

Definition at line 72 of file JetVertexChecker.cc.

Referenced by filter().

const double JetVertexChecker::m_pvErr_z

private

Definition at line 73 of file JetVertexChecker.cc.

Referenced by filter().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation