JetVertexChecker Class Reference

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

Inheritance diagram for JetVertexChecker:

Public Member Functions
	JetVertexChecker (const edm::ParameterSet &)
	~JetVertexChecker ()
Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default
Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDFilterBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	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
std::vector< ConsumesInfo >	consumesInfo () const
	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
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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::stream::EDFilterBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual 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
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 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

(

)

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 HLT_25ns14e33_v1_cff::InputTag.

                                                                             {
   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  )

overrideprivatevirtual

Implements edm::stream::EDFilterBase.

Definition at line 125 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_maxETA, m_maxNjets, m_maxTrackPt, m_newMethod, m_pvErr_x, m_pvErr_y, m_pvErr_z, 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(), 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(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);
   std::auto_ptr<reco::VertexCollection> pOut2(new reco::VertexCollection);
   pOut2->push_back(thePV);
   iEvent.put(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().