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

Inheritance diagram for SoftLepton:

Classes
struct	TrackCompare

Public Types
using	LeptonIds = std::map< unsigned int, float >

using	Leptons = std::map< edm::RefToBase< reco::Track >, LeptonIds, TrackCompare >

Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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

Public Member Functions
	SoftLepton (const edm::ParameterSet &iConfig)

reco::SoftLeptonTagInfo	tag (const edm::RefToBase< reco::Jet > &jet, const reco::TrackRefVector &tracks, const Leptons &leptons, const reco::Vertex &primaryVertex, const TransientTrackBuilder &builder) const

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () 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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
GlobalVector	refineJetAxis (const edm::RefToBase< reco::Jet > &jet, const reco::TrackRefVector &tracks, const edm::RefToBase< reco::Track > &exclude=edm::RefToBase< reco::Track >()) const

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Static Protected Member Functions
static double	boostedPPar (const math::XYZVector &vector, const math::XYZVector &axis)

static double	relativeEta (const math::XYZVector &vector, const math::XYZVector &axis)

Private Member Functions
void	produce (edm::StreamID, edm::Event &event, const edm::EventSetup &setup) const final

Private Attributes
const double	m_chi2Cut

const double	m_deltaRCut

const edm::InputTag	m_jets

const edm::InputTag	m_leptonCands

const edm::InputTag	m_leptonId

const edm::InputTag	m_leptons

const muon::SelectionType	m_muonSelection

const unsigned int	m_refineJetAxis

const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord >	token_builder

const edm::EDGetTokenT< edm::View< reco::Electron > >	token_electrons

const edm::EDGetTokenT< edm::View< reco::GsfElectron > >	token_gsfElectrons

const edm::EDGetTokenT< edm::View< reco::Jet > >	token_jets

const edm::EDGetTokenT< reco::JetTracksAssociationCollection >	token_jtas

const edm::EDGetTokenT< edm::ValueMap< float > >	token_leptonCands

const edm::EDGetTokenT< edm::ValueMap< float > >	token_leptonId

const edm::EDGetTokenT< edm::View< reco::Muon > >	token_muons

const edm::EDGetTokenT< reco::PFCandidateCollection >	token_pfElectrons

const edm::EDGetTokenT< reco::VertexCollection >	token_primaryVertex

const edm::EDPutTokenT< reco::SoftLeptonTagInfoCollection >	token_put

const edm::EDGetTokenT< edm::View< reco::Track > >	token_tracks

Static Private Attributes
static const reco::Vertex	s_nominalBeamSpot

Detailed Description

Description: CMSSW EDProducer for soft lepton b tagging.

Implementation: The actual tagging is performed by SoftLeptonAlgorithm.

Definition at line 65 of file SoftLepton.cc.

Member Typedef Documentation

◆ LeptonIds

using SoftLepton::LeptonIds = std::map<unsigned int, float>

Definition at line 77 of file SoftLepton.cc.

◆ Leptons

using SoftLepton::Leptons = std::map<edm::RefToBase<reco::Track>, LeptonIds, TrackCompare>

Definition at line 78 of file SoftLepton.cc.

Constructor & Destructor Documentation

◆ SoftLepton()

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

explicit

Definition at line 165 of file SoftLepton.cc.

     : m_jets(iConfig.getParameter<edm::InputTag>("jets")),
       token_jtas(mayConsume<reco::JetTracksAssociationCollection>(m_jets)),
       token_jets(mayConsume<edm::View<reco::Jet> >(m_jets)),
       token_primaryVertex(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("primaryVertex"))),
       m_leptons(iConfig.getParameter<edm::InputTag>("leptons")),
       token_gsfElectrons(mayConsume<GsfElectronView>(m_leptons)),
       token_electrons(mayConsume<ElectronView>(m_leptons)),
       token_pfElectrons(mayConsume<reco::PFCandidateCollection>(m_leptons)),
       token_muons(mayConsume<MuonView>(m_leptons)),
       token_tracks(mayConsume<edm::View<reco::Track> >(m_leptons)),
       m_leptonCands(iConfig.getParameter<edm::InputTag>("leptonCands")),
       token_leptonCands(mayConsume<edm::ValueMap<float> >(m_leptonCands)),
       m_leptonId(iConfig.getParameter<edm::InputTag>("leptonId")),
       token_leptonId(mayConsume<edm::ValueMap<float> >(m_leptonId)),
       token_builder(esConsumes(edm::ESInputTag("", "TransientTrackBuilder"))),
       token_put(produces()),
       m_refineJetAxis(iConfig.getParameter<unsigned int>("refineJetAxis")),
       m_deltaRCut(iConfig.getParameter<double>("leptonDeltaRCut")),
       m_chi2Cut(iConfig.getParameter<double>("leptonChi2Cut")),
       m_muonSelection((muon::SelectionType)iConfig.getParameter<unsigned int>("muonSelection")) {}

Member Function Documentation

◆ boostedPPar()

double SoftLepton::boostedPPar	(	const math::XYZVector &	vector,
		const math::XYZVector &	axis
	)

staticprotected

Definition at line 499 of file SoftLepton.cc.

References metsig::jet, and trackerHitRTTI::vector.

Referenced by tag().

                                                                                      {
   static const double lepton_mass = 0.00;  // assume a massless (ultrarelativistic) lepton
   static const double jet_mass = 5.279;    // use B±/B0 mass as the jet rest mass [PDG 2007 updates]
   ROOT::Math::LorentzVector<ROOT::Math::PxPyPzM4D<double> > lepton(
       vector.Dot(axis) / axis.r(), Perp(vector, axis), 0., lepton_mass);
   ROOT::Math::LorentzVector<ROOT::Math::PxPyPzM4D<double> > jet(axis.r(), 0., 0., jet_mass);
   ROOT::Math::BoostX boost(-jet.Beta());
   return boost(lepton).x();
 }

◆ fillDescriptions()

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

static

Definition at line 510 of file SoftLepton.cc.

References edm::ConfigurationDescriptions::addDefault(), submitPVResolutionJobs::desc, and HLT_2022v15_cff::InputTag.

                                                                             {
   edm::ParameterSetDescription desc;
   desc.add<unsigned int>("muonSelection", 1);
   desc.add<edm::InputTag>("leptons", edm::InputTag("muons"));
   desc.add<edm::InputTag>("primaryVertex", edm::InputTag("offlinePrimaryVertices"));
   desc.add<edm::InputTag>("leptonCands", edm::InputTag());
   desc.add<edm::InputTag>("leptonId", edm::InputTag());
   desc.add<unsigned int>("refineJetAxis", 0);
   desc.add<edm::InputTag>("jets", edm::InputTag("ak4PFJetsCHS"));
   desc.add<double>("leptonDeltaRCut", 0.4);
   desc.add<double>("leptonChi2Cut", 9999.0);
   descriptions.addDefault(desc);
 }

◆ produce()

void SoftLepton::produce	(	edm::StreamID	,
		edm::Event &	event,
		const edm::EventSetup &	setup
	)		const

finalprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 188 of file SoftLepton.cc.

References edm::View< T >::begin(), HPSPFTauProducerPuppi_cfi::electron, edm::View< T >::end(), Exception, edm::ValueMap< T >::get(), mps_fire::i, l1ctLayer2EG_cff::id, HLT_2022v15_cff::InputTag, muon::isGoodMuon(), edm::HandleBase::isValid(), PDWG_EXODelayedJetMET_cff::jets, HLT_2022v15_cff::leptonId, HLT_2022v15_cff::leptons, m_jets, m_leptonCands, m_leptonId, m_leptons, m_muonSelection, eostools::move(), HLT_2022v15_cff::muon, edm::errors::NotFound, edm::View< T >::refAt(), mps_fire::result, s_nominalBeamSpot, singleTopDQM_cfi::setup, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::size(), findQualityFiles::size, tag(), token_builder, token_electrons, token_gsfElectrons, token_jets, token_jtas, token_leptonCands, token_leptonId, token_muons, token_pfElectrons, token_primaryVertex, token_put, token_tracks, pwdgSkimBPark_cfi::tracks, and bphysicsOniaDQM_cfi::vertex.

                                                                                        {
   // grab a TransientTrack helper from the Event Setup
   auto const &transientTrackBuilder = setup.getData(token_builder);
 
   // input objects
 
   // input jets (and possibly tracks)
   ProductID jets_id;
   std::vector<edm::RefToBase<reco::Jet> > jets;
   std::vector<reco::TrackRefVector> tracks;
   do {
     {
       // look for a JetTracksAssociationCollection
       edm::Handle<reco::JetTracksAssociationCollection> h_jtas = event.getHandle(token_jtas);
       if (h_jtas.isValid()) {
         unsigned int size = h_jtas->size();
         jets.resize(size);
         tracks.resize(size);
         for (unsigned int i = 0; i < size; ++i) {
           jets[i] = (*h_jtas)[i].first;
           tracks[i] = (*h_jtas)[i].second;
         }
         break;
       }
     }
     {  // else...
       // look for a View<Jet>
       edm::Handle<edm::View<reco::Jet> > h_jets = event.getHandle(token_jets);
       if (h_jets.isValid()) {
         unsigned int size = h_jets->size();
         jets.resize(size);
         tracks.resize(size);
         for (unsigned int i = 0; i < h_jets->size(); i++)
           jets[i] = h_jets->refAt(i);
         break;
       }
     }
     {  // else...
       throw edm::Exception(edm::errors::NotFound)
           << "Object " << m_jets
           << " of type among (\"reco::JetTracksAssociationCollection\", \"edm::View<reco::Jet>\") not found";
     }
   } while (false);
 
   // input primary vetex
   reco::Vertex vertex;
   Handle<reco::VertexCollection> h_primaryVertex = event.getHandle(token_primaryVertex);
   if (h_primaryVertex.isValid() and not h_primaryVertex->empty())
     vertex = h_primaryVertex->front();
   else
     // fall back to nominal beam spot
     vertex = s_nominalBeamSpot;
 
   // input leptons (can be of different types)
   Leptons leptons;
 
   Handle<edm::ValueMap<float> > h_leptonCands;
   bool haveLeptonCands = !(m_leptonCands == edm::InputTag());
   if (haveLeptonCands)
     h_leptonCands = event.getHandle(token_leptonCands);
 
   // try to access the input collection as a collection of GsfElectrons, Muons or Tracks
 
   unsigned int leptonId = SoftLeptonProperties::Quality::leptonId;
   do {
     {
       // look for View<GsfElectron>
       Handle<GsfElectronView> h_electrons = event.getHandle(token_gsfElectrons);
 
       if (h_electrons.isValid()) {
         leptonId = SoftLeptonProperties::Quality::egammaElectronId;
         for (GsfElectronView::const_iterator electron = h_electrons->begin(); electron != h_electrons->end();
              ++electron) {
           LeptonIds &id = leptons[reco::TrackBaseRef(electron->gsfTrack())];
           id[SoftLeptonProperties::Quality::pfElectronId] = electron->mva_e_pi();
           if (haveLeptonCands)
             id[SoftLeptonProperties::Quality::btagElectronCands] =
                 (*h_leptonCands)[h_electrons->refAt(electron - h_electrons->begin())];
         }
         break;
       }
     }
     {  // else
       // look for View<Electron>
       // FIXME: is this obsolete?
       Handle<ElectronView> h_electrons = event.getHandle(token_electrons);
       if (h_electrons.isValid()) {
         leptonId = SoftLeptonProperties::Quality::egammaElectronId;
         for (ElectronView::const_iterator electron = h_electrons->begin(); electron != h_electrons->end(); ++electron) {
           LeptonIds &id = leptons[reco::TrackBaseRef(electron->track())];
           if (haveLeptonCands)
             id[SoftLeptonProperties::Quality::btagElectronCands] =
                 (*h_leptonCands)[h_electrons->refAt(electron - h_electrons->begin())];
         }
         break;
       }
     }
     {  // else
       // look for PFElectrons
       // FIXME: is this obsolete?
       Handle<reco::PFCandidateCollection> h_electrons = event.getHandle(token_pfElectrons);
       if (h_electrons.isValid()) {
         leptonId = SoftLeptonProperties::Quality::egammaElectronId;
         for (reco::PFCandidateCollection::const_iterator electron = h_electrons->begin();
              electron != h_electrons->end();
              ++electron) {
           LeptonIds *id;
           if (electron->gsfTrackRef().isNonnull())
             id = &leptons[reco::TrackBaseRef(electron->gsfTrackRef())];
           else if (electron->trackRef().isNonnull())
             id = &leptons[reco::TrackBaseRef(electron->trackRef())];
           else
             continue;
           (*id)[SoftLeptonProperties::Quality::pfElectronId] = electron->mva_e_pi();
           if (haveLeptonCands)
             (*id)[SoftLeptonProperties::Quality::btagElectronCands] =
                 (*h_leptonCands)[reco::PFCandidateRef(h_electrons, electron - h_electrons->begin())];
         }
         break;
       }
     }
     {  // else
       // look for View<Muon>
       Handle<MuonView> h_muons = event.getHandle(token_muons);
       if (h_muons.isValid()) {
         for (MuonView::const_iterator muon = h_muons->begin(); muon != h_muons->end(); ++muon) {
           // FIXME -> turn this selection into a muonCands input?
           if (muon::isGoodMuon(*muon, m_muonSelection)) {
             LeptonIds *id;
             if (muon->globalTrack().isNonnull())
               id = &leptons[reco::TrackBaseRef(muon->globalTrack())];
             else if (muon->innerTrack().isNonnull())
               id = &leptons[reco::TrackBaseRef(muon->innerTrack())];
             else if (muon->outerTrack().isNonnull())
               // does this makes sense ?
               id = &leptons[reco::TrackBaseRef(muon->outerTrack())];
             else
               continue;
             if (haveLeptonCands)
               (*id)[SoftLeptonProperties::Quality::btagMuonCands] =
                   (*h_leptonCands)[h_muons->refAt(muon - h_muons->begin())];
           }
         }
         break;
       }
     }
     {  // else
       // look for edm::View<Track>
       Handle<edm::View<reco::Track> > h_tracks = event.getHandle(token_tracks);
       if (h_tracks.isValid()) {
         for (unsigned int i = 0; i < h_tracks->size(); i++) {
           LeptonIds &id = leptons[h_tracks->refAt(i)];
           if (haveLeptonCands)
             id[SoftLeptonProperties::Quality::btagLeptonCands] = (*h_leptonCands)[h_tracks->refAt(i)];
         }
         break;
       }
     }
     {  // else
       throw edm::Exception(edm::errors::NotFound) << "Object " << m_leptons
                                                   << " of type among (\"edm::View<reco::GsfElectron>\", "
                                                      "\"edm::View<reco::Muon>\", \"edm::View<reco::Track>\") !found";
     }
   } while (false);
 
   if (!(m_leptonId == edm::InputTag())) {
     edm::ValueMap<float> const &h_leptonId = event.get(token_leptonId);
 
     for (Leptons::iterator lepton = leptons.begin(); lepton != leptons.end(); ++lepton)
       lepton->second[leptonId] = h_leptonId[lepton->first];
   }
 
   // output collections
   reco::SoftLeptonTagInfoCollection outputCollection;
   for (unsigned int i = 0; i < jets.size(); ++i) {
     reco::SoftLeptonTagInfo result = tag(jets[i], tracks[i], leptons, vertex, transientTrackBuilder);
     outputCollection.push_back(result);
   }
   event.emplace(token_put, std::move(outputCollection));
 }

◆ refineJetAxis()

GlobalVector SoftLepton::refineJetAxis	(	const edm::RefToBase< reco::Jet > &	jet,
		const reco::TrackRefVector &	tracks,
		const edm::RefToBase< reco::Track > &	exclude = `edm::RefToBase<reco::Track>()`
	)		const

protected

Definition at line 415 of file SoftLepton.cc.

References AXIS_CALORIMETRIC_NOLEPTON, AXIS_CHARGED_AVERAGE, AXIS_CHARGED_AVERAGE_NOLEPTON, AXIS_CHARGED_SUM, AXIS_CHARGED_SUM_NOLEPTON, edm::RefToBase< T >::isNonnull(), metsig::jet, M_PI, m_refineJetAxis, reco::TrackBase::momentum(), or, perp(), HLT_2022v15_cff::track, and pwdgSkimBPark_cfi::tracks.

Referenced by tag().

         {
   math::XYZVector axis = jet->momentum();
 
   if (m_refineJetAxis == AXIS_CHARGED_AVERAGE or m_refineJetAxis == AXIS_CHARGED_AVERAGE_NOLEPTON) {
     double sum_pT = 0.;
     double sum_eta_by_pT = 0.;
     double sum_phi_by_pT = 0.;
 
     double perp;
     double phi_rel;
     double eta_rel;
 
     // refine jet eta and phi with charged tracks measurements, if available
     for (reco::TrackRefVector::const_iterator track_it = tracks.begin(); track_it != tracks.end(); ++track_it) {
       const reco::Track &track = **track_it;
 
       perp = track.pt();
       eta_rel = (double)track.eta() - axis.eta();
       phi_rel = (double)track.phi() - axis.phi();
       while (phi_rel < -M_PI)
         phi_rel += 2 * M_PI;
       while (phi_rel > M_PI)
         phi_rel -= 2 * M_PI;
 
       sum_pT += perp;
       sum_phi_by_pT += perp * phi_rel;
       sum_eta_by_pT += perp * eta_rel;
     }
 
     // "remove" excluded track
     if (m_refineJetAxis == AXIS_CHARGED_AVERAGE_NOLEPTON and exclude.isNonnull()) {
       const reco::Track &track = *exclude;
 
       perp = track.pt();
       eta_rel = (double)track.eta() - axis.eta();
       phi_rel = (double)track.phi() - axis.phi();
       while (phi_rel < -M_PI)
         phi_rel += 2 * M_PI;
       while (phi_rel > M_PI)
         phi_rel -= 2 * M_PI;
 
       sum_pT -= perp;
       sum_phi_by_pT -= perp * phi_rel;
       sum_eta_by_pT -= perp * eta_rel;
     }
 
     if (sum_pT > 1.)  // avoid the case of only the lepton-track with small rounding errors
       axis =
           math::RhoEtaPhiVector(axis.rho(), axis.eta() + sum_eta_by_pT / sum_pT, axis.phi() + sum_phi_by_pT / sum_pT);
 
   } else if (m_refineJetAxis == AXIS_CHARGED_SUM or m_refineJetAxis == AXIS_CHARGED_SUM_NOLEPTON) {
     math::XYZVector sum;
 
     // recalculate the jet direction as the sum of charget tracks momenta
     for (reco::TrackRefVector::const_iterator track_it = tracks.begin(); track_it != tracks.end(); ++track_it) {
       const reco::Track &track = **track_it;
       sum += track.momentum();
     }
 
     // "remove" excluded track
     if (m_refineJetAxis == AXIS_CHARGED_SUM_NOLEPTON and exclude.isNonnull()) {
       const reco::Track &track = *exclude;
       sum -= track.momentum();
     }
 
     if (sum.R() > 1.)  // avoid the case of only the lepton-track with small rounding errors
       axis = sum;
   } else if (m_refineJetAxis == AXIS_CALORIMETRIC_NOLEPTON) {
     axis -= exclude->momentum();
   }
 
   return GlobalVector(axis.x(), axis.y(), axis.z());
 }

◆ relativeEta()

double SoftLepton::relativeEta	(	const math::XYZVector &	vector,
		const math::XYZVector &	axis
	)

staticprotected

Definition at line 492 of file SoftLepton.cc.

References dot(), dqm-mbProfile::log, mag(), and trackerHitRTTI::vector.

Referenced by tag().

                                                                                      {
   double mag = vector.r() * axis.r();
   double dot = vector.Dot(axis);
   return -log((mag - dot) / (mag + dot)) / 2;
 }

◆ tag()

reco::SoftLeptonTagInfo SoftLepton::tag	(	const edm::RefToBase< reco::Jet > &	jet,
		const reco::TrackRefVector &	tracks,
		const Leptons &	leptons,
		const reco::Vertex &	primaryVertex,
		const TransientTrackBuilder &	builder
	)		const

Definition at line 370 of file SoftLepton.cc.

Referenced by produce().

                                                                                                   {
   reco::SoftLeptonTagInfo info;
   info.setJetRef(jet);
 
   for (Leptons::const_iterator lepton = leptons.begin(); lepton != leptons.end(); ++lepton) {
     const math::XYZVector &lepton_momentum = lepton->first->momentum();
     if (m_chi2Cut > 0.0 && lepton->first->normalizedChi2() > m_chi2Cut)
       continue;
 
     const GlobalVector jetAxis = refineJetAxis(jet, tracks, lepton->first);
     const math::XYZVector axis(jetAxis.x(), jetAxis.y(), jetAxis.z());
     float deltaR = Geom::deltaR(lepton_momentum, axis);
     if (deltaR > m_deltaRCut)
       continue;
 
     reco::SoftLeptonProperties properties;
 
     reco::TransientTrack transientTrack = transientTrackBuilder.build(*lepton->first);
     Measurement1D ip2d = IPTools::signedTransverseImpactParameter(transientTrack, jetAxis, primaryVertex).second;
     Measurement1D ip3d = IPTools::signedImpactParameter3D(transientTrack, jetAxis, primaryVertex).second;
     properties.sip2dsig = ip2d.significance();
     properties.sip3dsig = ip3d.significance();
     properties.sip2d = ip2d.value();
     properties.sip3d = ip3d.value();
     properties.deltaR = deltaR;
     properties.ptRel = Perp(lepton_momentum, axis);
     properties.p0Par = boostedPPar(lepton_momentum, axis);
     properties.etaRel = relativeEta(lepton_momentum, axis);
     properties.ratio = lepton_momentum.R() / axis.R();
     properties.ratioRel = lepton_momentum.Dot(axis) / axis.Mag2();
 
     for (LeptonIds::const_iterator iter = lepton->second.begin(); iter != lepton->second.end(); ++iter)
       properties.setQuality(static_cast<SoftLeptonProperties::Quality::Generic>(iter->first), iter->second);
 
     info.insert(lepton->first, properties);
   }
 
   return info;
 }

Member Data Documentation

◆ m_chi2Cut

const double SoftLepton::m_chi2Cut

private

Definition at line 124 of file SoftLepton.cc.

Referenced by tag().

◆ m_deltaRCut

const double SoftLepton::m_deltaRCut

private

Definition at line 123 of file SoftLepton.cc.

Referenced by tag().

◆ m_jets

const edm::InputTag SoftLepton::m_jets

private

Definition at line 102 of file SoftLepton.cc.

Referenced by produce().

◆ m_leptonCands

const edm::InputTag SoftLepton::m_leptonCands

private

Definition at line 112 of file SoftLepton.cc.

Referenced by produce().

◆ m_leptonId

const edm::InputTag SoftLepton::m_leptonId

private

Definition at line 114 of file SoftLepton.cc.

Referenced by produce().

◆ m_leptons

const edm::InputTag SoftLepton::m_leptons

private

Definition at line 106 of file SoftLepton.cc.

Referenced by produce().

◆ m_muonSelection

const muon::SelectionType SoftLepton::m_muonSelection

private

Definition at line 127 of file SoftLepton.cc.

Referenced by produce().

◆ m_refineJetAxis

const unsigned int SoftLepton::m_refineJetAxis

private

Definition at line 122 of file SoftLepton.cc.

Referenced by refineJetAxis().

◆ s_nominalBeamSpot

const reco::Vertex SoftLepton::s_nominalBeamSpot

staticprivate

Definition at line 130 of file SoftLepton.cc.

Referenced by produce().

◆ token_builder

const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> SoftLepton::token_builder

private

Definition at line 118 of file SoftLepton.cc.

Referenced by produce().

◆ token_electrons

const edm::EDGetTokenT<edm::View<reco::Electron> > SoftLepton::token_electrons

private

Definition at line 108 of file SoftLepton.cc.

Referenced by produce().

◆ token_gsfElectrons

const edm::EDGetTokenT<edm::View<reco::GsfElectron> > SoftLepton::token_gsfElectrons

private

Definition at line 107 of file SoftLepton.cc.

Referenced by produce().

◆ token_jets

const edm::EDGetTokenT<edm::View<reco::Jet> > SoftLepton::token_jets

private

Definition at line 104 of file SoftLepton.cc.

Referenced by produce().

◆ token_jtas

const edm::EDGetTokenT<reco::JetTracksAssociationCollection> SoftLepton::token_jtas

private

Definition at line 103 of file SoftLepton.cc.

Referenced by produce().

◆ token_leptonCands

const edm::EDGetTokenT<edm::ValueMap<float> > SoftLepton::token_leptonCands

private

Definition at line 113 of file SoftLepton.cc.

Referenced by produce().

◆ token_leptonId

const edm::EDGetTokenT<edm::ValueMap<float> > SoftLepton::token_leptonId

private

Definition at line 115 of file SoftLepton.cc.

Referenced by produce().

◆ token_muons

const edm::EDGetTokenT<edm::View<reco::Muon> > SoftLepton::token_muons

private

Definition at line 110 of file SoftLepton.cc.

Referenced by produce().

◆ token_pfElectrons

const edm::EDGetTokenT<reco::PFCandidateCollection> SoftLepton::token_pfElectrons

private

Definition at line 109 of file SoftLepton.cc.

Referenced by produce().

◆ token_primaryVertex

const edm::EDGetTokenT<reco::VertexCollection> SoftLepton::token_primaryVertex

private

Definition at line 105 of file SoftLepton.cc.

Referenced by produce().

◆ token_put

const edm::EDPutTokenT<reco::SoftLeptonTagInfoCollection> SoftLepton::token_put

private

Definition at line 120 of file SoftLepton.cc.

Referenced by produce().

◆ token_tracks

const edm::EDGetTokenT<edm::View<reco::Track> > SoftLepton::token_tracks

private

Definition at line 111 of file SoftLepton.cc.

Referenced by produce().

Classes

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Static Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ LeptonIds

◆ Leptons

Constructor & Destructor Documentation

◆ SoftLepton()

Member Function Documentation

◆ boostedPPar()

◆ fillDescriptions()

◆ produce()

◆ refineJetAxis()

◆ relativeEta()

◆ tag()

Member Data Documentation

◆ m_chi2Cut

◆ m_deltaRCut

◆ m_jets

◆ m_leptonCands

◆ m_leptonId

◆ m_leptons

◆ m_muonSelection

◆ m_refineJetAxis

◆ s_nominalBeamSpot

◆ token_builder

◆ token_electrons

◆ token_gsfElectrons

◆ token_jets

◆ token_jtas

◆ token_leptonCands

◆ token_leptonId

◆ token_muons

◆ token_pfElectrons

◆ token_primaryVertex

◆ token_put

◆ token_tracks