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

Inheritance diagram for PixelJetPuId:

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

virtual	~PixelJetPuId ()

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	produce (edm::StreamID sid, edm::Event &, const edm::EventSetup &) const override

Private Attributes
edm::EDGetTokenT< edm::View < reco::Jet > >	generaljetsToken

edm::EDGetTokenT< edm::View < reco::CaloJet > >	jetsToken

bool	m_fwjets

edm::InputTag	m_jets

double	m_MaxTrackChi2

double	m_MaxTrackDistanceToJet

double	m_minet_fwjets

double	m_mineta_fwjets

double	m_MinGoodJetTrackPt

double	m_MinGoodJetTrackPtRatio

double	m_MinTrackPt

edm::InputTag	m_primaryVertex

edm::InputTag	m_tracks

edm::EDGetTokenT < reco::VertexCollection >	primaryVertexToken

edm::EDGetTokenT< std::vector < reco::Track > >	tracksToken

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	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: The PixelJetPuId module select all the pixel tracks compatible with a jet. If the sum of the tracks momentum is under a threshold the jet is tagged as "PUjets".

Implementation: [Notes on implementation]

Definition at line 56 of file PixelJetPuId.cc.

Constructor & Destructor Documentation

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

Definition at line 91 of file PixelJetPuId.cc.

References generaljetsToken, edm::ParameterSet::getParameter(), jetsToken, m_fwjets, m_jets, m_MaxTrackChi2, m_MaxTrackDistanceToJet, m_minet_fwjets, m_mineta_fwjets, m_MinGoodJetTrackPt, m_MinGoodJetTrackPtRatio, m_MinTrackPt, m_primaryVertex, m_tracks, primaryVertexToken, and tracksToken.

 {
     //InputTag
     m_tracks           = iConfig.getParameter<edm::InputTag>("tracks");
     tracksToken        = consumes<std::vector<reco::Track> >(m_tracks);
     m_jets             = iConfig.getParameter<edm::InputTag>("jets");
     jetsToken          = consumes<edm::View<reco::CaloJet> >(m_jets);
     generaljetsToken   = consumes<edm::View<reco::Jet> >(m_jets);
     m_primaryVertex    = iConfig.getParameter<edm::InputTag>("primaryVertex");
     primaryVertexToken = consumes<reco::VertexCollection>(m_primaryVertex);
     
     //Tracks Selection
     m_MinTrackPt             = iConfig.getParameter<double>("MinTrackPt");
     m_MaxTrackDistanceToJet  = iConfig.getParameter<double>("MaxTrackDistanceToJet");
     m_MaxTrackChi2           = iConfig.getParameter<double>("MaxTrackChi2");
     
     //A jet is defined as a signal jet if Sum(trackPt) > minPt or Sum(comp.trackPt)/CaloJetPt > minPtRatio
     m_MinGoodJetTrackPt      = iConfig.getParameter<double>("MinGoodJetTrackPt");
     m_MinGoodJetTrackPtRatio = iConfig.getParameter<double>("MinGoodJetTrackPtRatio");
     
     m_fwjets                 = iConfig.getParameter<bool>("UseForwardJetsAsNoPU");
     m_mineta_fwjets          = iConfig.getParameter<double>("MinEtaForwardJets");
     m_minet_fwjets           = iConfig.getParameter<double>("MinEtForwardJets");
     
     produces<std::vector<reco::CaloJet> >(); 
     produces<std::vector<reco::CaloJet> >("PUjets"); 
     produces<reco::JetTagCollection>();
 }

PixelJetPuId::~PixelJetPuId ( )

virtual

Definition at line 121 of file PixelJetPuId.cc.

121 {}

Member Function Documentation

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

static

Definition at line 125 of file PixelJetPuId.cc.

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

                                                                          {
     edm::ParameterSetDescription desc;
     desc.add<edm::InputTag> ("jets",edm::InputTag("hltCaloJetL1FastJetCorrected"));
     desc.add<edm::InputTag> ("tracks",edm::InputTag("hltPixelTracksNoPU"));
     desc.add<edm::InputTag> ("primaryVertex",edm::InputTag("hltFastPVPixelVertices"));
     desc.add<double>("MinGoodJetTrackPtRatio",0.045);
     desc.add<double>("MinGoodJetTrackPt",1.8);
     desc.add<double>("MaxTrackDistanceToJet",0.04);
     desc.add<double>("MinTrackPt",0.6);
     desc.add<double>("MaxTrackChi2",20.);
     desc.add<bool>("UseForwardJetsAsNoPU",true);
     desc.add<double>("MinEtaForwardJets",2.4);
     desc.add<double>("MinEtForwardJets",40.);
     descriptions.add("pixelJetPuId",desc);
 }

void PixelJetPuId::produce	(	edm::StreamID	sid,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 146 of file PixelJetPuId.cc.

References reco::deltaR2(), generaljetsToken, edm::EventSetup::get(), edm::Event::getByToken(), i, fwrapper::jets, jetsToken, IPTools::jetTrackDistance(), findQualityFiles::jj, m_fwjets, m_MaxTrackChi2, m_MaxTrackDistanceToJet, m_minet_fwjets, m_mineta_fwjets, m_MinGoodJetTrackPt, m_MinGoodJetTrackPtRatio, m_MinTrackPt, edm::makeRefToBaseProdFrom(), eostools::move(), HLT_25ns10e33_v2_cff::primaryVertex, primaryVertexToken, edm::Event::put(), MetAnalyzer::pv(), edm::second(), testEve_cfg::tracks, and tracksToken.

 {
     using namespace edm;
     std::unique_ptr<std::vector<reco::CaloJet> > pOut(new std::vector<reco::CaloJet> );
     std::unique_ptr<std::vector<reco::CaloJet> > pOut_PUjets(new std::vector<reco::CaloJet> );
     std::unique_ptr<reco::JetTagCollection > pOut_jetTagCollection(new reco::JetTagCollection );
     
     //get tracks
     Handle<std::vector<reco::Track> > tracks;
     iEvent.getByToken(tracksToken, tracks);
     unsigned int tsize = tracks->size();
     float teta[tsize], tphi[tsize];
     unsigned int i=0;
     for (auto const & tr : *tracks) { teta[i]=tr.eta(); tphi[i]=tr.phi();++i;}
     
     //get jets
     Handle<edm::View<reco::CaloJet> > jets;
     iEvent.getByToken(jetsToken, jets);
 
     Handle<edm::View<reco::Jet> > generaljets;
     iEvent.getByToken(generaljetsToken, generaljets);
     
     //get primary vertices
     Handle<reco::VertexCollection> primaryVertex;
     iEvent.getByToken(primaryVertexToken, primaryVertex);
     
     //get Transient Track Builder
     edm::ESHandle<TransientTrackBuilder> builder;
     iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", builder);
     
     //init JetTagCollection
     if(generaljets.product()->size()>0)
     {
         edm::RefToBase<reco::Jet> jj = edm::RefToBase<reco::Jet>(generaljets,0);
         pOut_jetTagCollection.reset(new reco::JetTagCollection(edm::makeRefToBaseProdFrom(jj, iEvent)));
     }
      
     //loop on trackIPTagInfos
     if(primaryVertex->size()>0)
     {
         const reco::Vertex* pv = &*primaryVertex->begin();
         //loop on jets
         for(edm::View<reco::CaloJet>::const_iterator itJet = jets->begin(); itJet != jets->end(); itJet++ ) {
             
             math::XYZVector jetMomentum = itJet->momentum();
             GlobalVector direction(jetMomentum.x(), jetMomentum.y(), jetMomentum.z());
             
             math::XYZVector trMomentum;
             
                 if(fabs(itJet->eta())>m_mineta_fwjets)
                 {
                     if((m_fwjets) && (itJet->et()>m_minet_fwjets))
                         pOut->push_back(*itJet);// fill forward jet as signal jet
                 }
                 else 
                 {
                     //loop on tracks
                     std::vector<reco::Track>::const_iterator itTrack = tracks->begin();
                     for (unsigned int i=0; i<tsize; ++i) {
                         float deltaR2=reco::deltaR2(itJet->eta(),itJet->phi(), teta[i],tphi[i]);
                         if(deltaR2<0.25) {
                             reco::TransientTrack transientTrack = builder->build(*itTrack);
                             float jetTrackDistance = -((IPTools::jetTrackDistance(transientTrack, direction, *pv)).second).value();
                             
                             //select the tracks compabible with the jet
                             if(( itTrack->pt() > m_MinTrackPt) && ( itTrack->normalizedChi2() < m_MaxTrackChi2) && (jetTrackDistance<m_MaxTrackDistanceToJet))
                             {
                                 trMomentum += itTrack->momentum(); //calculate the Sum(trackPt)
                             }
                         }
                         itTrack++;
                     }
                     //if Sum(comp.trackPt)/CaloJetPt > minPtRatio or Sum(trackPt) > minPt  the jet is a signal jet
                     if(trMomentum.rho()/jetMomentum.rho() > m_MinGoodJetTrackPtRatio || trMomentum.rho() > m_MinGoodJetTrackPt ) 
                     {
                         pOut->push_back(*itJet);        // fill it as signal jet
                     }
                     else//else it is a PUjet
                     {
                         pOut_PUjets->push_back(*itJet); // fill it as PUjets
                     }
                 }
             RefToBase<reco::Jet> jRef(generaljets, itJet-jets->begin());
             (*pOut_jetTagCollection)[jRef] = trMomentum.rho();  // fill jetTagCollection
         }
     }
     iEvent.put(std::move(pOut));
     iEvent.put(std::move(pOut_PUjets),"PUjets");
     iEvent.put(std::move(pOut_jetTagCollection));
 }