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

Inheritance diagram for ConeIsolation:

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

void	produce (edm::Event &, const edm::EventSetup &) override

	~ConeIsolation () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Private Attributes
edm::InputTag	beamSpotTag

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken

edm::InputTag	jetTrackTag

edm::EDGetTokenT< reco::JetTracksAssociationCollection >	jetTrackToken

ConeIsolationAlgorithm *	m_algo

bool	usingBeamSpot

edm::InputTag	vertexTag

edm::EDGetTokenT< reco::VertexCollection >	vertexToken

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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>="">

EDProducer of the tagged TauJet with the ConeIsolationAlgorithm. It returns two collections: base collection is the JetTag, and extended Collection which is the IsolatedTauTagInfo. The method implemented in the IsolatedTauTagInfo class are used to compute the discriminator variable. A trick is used to link the IsolatedTauTagInfo to a smart reference to the JetTag.

Author: Simone Gennai

Definition at line 24 of file ConeIsolation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 53 of file ConeIsolation.cc.

References AlignmentProducer_cff::beamSpotTag, edm::ParameterSet::getParameter(), and electronIdMVABased_cfi::vertexTag.

 {
   jetTrackTag   = iConfig.getParameter<InputTag>("JetTrackSrc");
   jetTrackToken = consumes<reco::JetTracksAssociationCollection>(jetTrackTag);
   vertexTag     = iConfig.getParameter<InputTag>("vertexSrc");
   vertexToken   = consumes<reco::VertexCollection>(vertexTag);
   beamSpotTag   = iConfig.getParameter<edm::InputTag>("BeamSpotProducer");
   beamSpotToken = consumes<reco::BeamSpot>(beamSpotTag);
   usingBeamSpot = iConfig.getParameter<bool>("useBeamSpot"); //If false the OfflinePrimaryVertex will be used.
 
   m_algo = new ConeIsolationAlgorithm(iConfig);
   
   produces<reco::JetTagCollection>(); 
   produces<reco::IsolatedTauTagInfoCollection>();       
 
 }

ConeIsolation::~ConeIsolation ( )

override

Definition at line 71 of file ConeIsolation.cc.

 {
   delete m_algo;
 }

Member Function Documentation

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

static

Definition at line 83 of file ConeIsolation.cc.

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

                                                                           {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("JetTrackSrc",edm::InputTag("ic5JetTracksAssociatorAtVertex"));
   desc.add<edm::InputTag>("vertexSrc",edm::InputTag("offlinePrimaryVertices"));
   desc.add<edm::InputTag>("BeamSpotProducer",edm::InputTag("offlineBeamSpot"));
   desc.add<bool>("useBeamSpot",false);
   ConeIsolationAlgorithm::fillDescription(desc);
   descriptions.add("coneIsolationTauJetTags", desc);
 }

void ConeIsolation::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 95 of file ConeIsolation.cc.

References edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::begin(), reco::Vertex::chi2(), MillePedeFileConverter_cfg::e, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::empty(), reco::Vertex::error(), edm::Event::getByToken(), mps_fire::i, iEvent, edm::makeRefToBaseProdFrom(), eostools::move(), reco::Vertex::ndof(), AlCaHLTBitMon_ParallelJobs::p, parseEventContent::prod, edm::Handle< T >::product(), edm::Event::put(), edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::size(), reco::Vertex::tracksSize(), electrons_cff::vertices, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::Vertex::z().

 {
    using namespace edm;
    //Get jets with tracks
    Handle<reco::JetTracksAssociationCollection> jetTracksAssociation;
    iEvent.getByToken(jetTrackToken,jetTracksAssociation);
 
    std::unique_ptr<reco::JetTagCollection>             tagCollection;
    std::unique_ptr<reco::IsolatedTauTagInfoCollection> extCollection( new reco::IsolatedTauTagInfoCollection() );
    if (not jetTracksAssociation->empty()) {
      RefToBaseProd<reco::Jet> prod( edm::makeRefToBaseProdFrom(jetTracksAssociation->begin()->first, iEvent) );
      tagCollection.reset( new reco::JetTagCollection(prod) );
    } else {
      tagCollection.reset( new reco::JetTagCollection() );
    }
 
    Vertex::Error e;
    e(0,0)=1;
    e(1,1)=1;
    e(2,2)=1;
    Vertex::Point p(0,0,-1000);
    Vertex myPVtmp(p,e);//Fake vertex to be used in case no vertex is found
    Vertex myPV;
 
    //Get pixel vertices
    Handle<reco::VertexCollection> vertices;
    iEvent.getByToken(vertexToken,vertices);
    const reco::VertexCollection vertCollection = *(vertices.product());
    //Check if there is the PV!!!!
    if(vertCollection.begin() != vertCollection.end())
      myPVtmp = *(vertCollection.begin());
 
    //In case the beam spot is used, the Z of the vertex still comes from the PV, while the (x,y) is taken from the beamspot
    reco::BeamSpot vertexBeamSpot;
    edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
    
    if(usingBeamSpot)
      {
        //Create a new vertex with the information on x0 and Y0 from the beamspot, to be used in HLT.
        iEvent.getByToken(beamSpotToken,recoBeamSpotHandle);
        vertexBeamSpot = *recoBeamSpotHandle;
        Vertex::Point bspoint(vertexBeamSpot.x0(),vertexBeamSpot.y0(),myPVtmp.z());
        Vertex combinedVertex = Vertex(bspoint,myPVtmp.error(),myPVtmp.chi2(),myPVtmp.ndof(),myPVtmp.tracksSize());
        myPV = combinedVertex;
      }else{
        myPV = myPVtmp;
      }
    
    for (unsigned int i = 0; i < jetTracksAssociation->size(); ++i)
      {
        pair<float,IsolatedTauTagInfo> myPair =m_algo->tag(edm::Ref<JetTracksAssociationCollection>(jetTracksAssociation,i),myPV); 
        tagCollection->setValue(i, myPair.first);
        extCollection->push_back(myPair.second);
      }
 
    iEvent.put(std::move(extCollection));
    iEvent.put(std::move(tagCollection));
 
 }

Member Data Documentation

edm::InputTag ConeIsolation::beamSpotTag

private

Definition at line 36 of file ConeIsolation.h.

edm::EDGetTokenT<reco::BeamSpot> ConeIsolation::beamSpotToken

private

Definition at line 37 of file ConeIsolation.h.

edm::InputTag ConeIsolation::jetTrackTag

private

Definition at line 32 of file ConeIsolation.h.

edm::EDGetTokenT<reco::JetTracksAssociationCollection> ConeIsolation::jetTrackToken

private

Definition at line 33 of file ConeIsolation.h.

ConeIsolationAlgorithm* ConeIsolation::m_algo

private

Definition at line 39 of file ConeIsolation.h.

bool ConeIsolation::usingBeamSpot

private

Definition at line 38 of file ConeIsolation.h.

edm::InputTag ConeIsolation::vertexTag

private

Definition at line 34 of file ConeIsolation.h.

edm::EDGetTokenT<reco::VertexCollection> ConeIsolation::vertexToken

private

Definition at line 35 of file ConeIsolation.h.

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation