#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::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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
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::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Private Attributes
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::EDProducer
typedef EDProducer	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: 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 54 of file PixelJetPuId.cc.

Constructor & Destructor Documentation

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

Definition at line 89 of file PixelJetPuId.cc.

References 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);
   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"); 
 }

PixelJetPuId::~PixelJetPuId ( )

virtual

Definition at line 117 of file PixelJetPuId.cc.

117 {}

Member Function Documentation

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

static

Definition at line 121 of file PixelJetPuId.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_25ns14e33_v1_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::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 142 of file PixelJetPuId.cc.

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

 {
   using namespace edm;
   std::auto_ptr<std::vector<reco::CaloJet> > pOut(new std::vector<reco::CaloJet> );
   std::auto_ptr<std::vector<reco::CaloJet> > pOut_PUjets(new std::vector<reco::CaloJet> );
   
    //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);
    
   //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);
 
   //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;
     
     //loop on tracks
     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 
       {
         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
           }
       }
       }
     }
   iEvent.put(pOut);
   iEvent.put(pOut_PUjets,"PUjets");
   
 }