#include <RecoBTag/SoftLepton/plugin/SoftLepton.h>

Inheritance diagram for SoftLepton:

Classes
struct	TrackCompare

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

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

Public Types inherited from edm::stream::EDProducer<>
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::EDProducerBase
typedef EDProducerAdaptorBase	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

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

	~SoftLepton ()

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

Public Member Functions inherited from edm::stream::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

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)

virtual	~ProducerBase () noexcept(false)

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

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::stream::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::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)

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
virtual void	produce (edm::Event &event, const edm::EventSetup &setup)

Private Attributes
double	m_chi2Cut

double	m_deltaRCut

const edm::InputTag	m_jets

const edm::InputTag	m_leptonCands

const edm::InputTag	m_leptonId

const edm::InputTag	m_leptons

muon::SelectionType	m_muonSelection

unsigned int	m_refineJetAxis

const TransientTrackBuilder *	m_transientTrackBuilder

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::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:

Description: CMSSW EDProducer wrapper for sot lepton b tagging.

Implementation: The actual tagging is performed by SoftLeptonAlgorithm.

Definition at line 51 of file SoftLepton.h.

Member Typedef Documentation

typedef std::map<unsigned int, float> SoftLepton::LeptonIds

Definition at line 65 of file SoftLepton.h.

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

Definition at line 66 of file SoftLepton.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 66 of file SoftLepton.cc.

                                                       :
   m_jets(              iConfig.getParameter<edm::InputTag>( "jets" ) ),
   token_jtas(          mayConsume<reco::JetTracksAssociationCollection>( iConfig.getParameter<edm::InputTag>( "jets" ) ) ),
   token_jets(          mayConsume<edm::View<reco::Jet> >( iConfig.getParameter<edm::InputTag>( "jets" ) ) ),
   token_primaryVertex( consumes<reco::VertexCollection>( iConfig.getParameter<edm::InputTag>( "primaryVertex" ) ) ),
   m_leptons(           iConfig.getParameter<edm::InputTag>( "leptons" ) ),
   token_gsfElectrons(  mayConsume<GsfElectronView>( iConfig.getParameter<edm::InputTag>( "leptons" ) ) ),
   token_electrons(     mayConsume<ElectronView>( iConfig.getParameter<edm::InputTag>( "leptons" ) ) ),
   token_pfElectrons(   mayConsume<reco::PFCandidateCollection>( iConfig.getParameter<edm::InputTag>( "leptons" ) ) ),
   token_muons(         mayConsume<MuonView>( iConfig.getParameter<edm::InputTag>( "leptons" ) ) ),
   token_tracks(        mayConsume<edm::View<reco::Track> >( iConfig.getParameter<edm::InputTag>( "leptons" ) ) ),
   m_leptonCands(       iConfig.exists("leptonCands") ? iConfig.getParameter<edm::InputTag>( "leptonCands" ) : edm::InputTag() ),
   token_leptonCands(   mayConsume<edm::ValueMap<float> >( iConfig.exists("leptonCands") ? iConfig.getParameter<edm::InputTag>( "leptonCands" ) : edm::InputTag() ) ),
   m_leptonId(          iConfig.exists("leptonId") ? iConfig.getParameter<edm::InputTag>( "leptonId" ) : edm::InputTag() ),
   token_leptonId(      mayConsume<edm::ValueMap<float> >( iConfig.exists("leptonId") ? iConfig.getParameter<edm::InputTag>( "leptonId" ) : edm::InputTag() ) ),
   m_transientTrackBuilder( 0 ),
   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" ) )
 {
   produces<reco::SoftLeptonTagInfoCollection>();
 }

SoftLepton::~SoftLepton ( void )

Definition at line 91 of file SoftLepton.cc.

91 {

92 }

Member Function Documentation

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

staticprotected

Definition at line 398 of file SoftLepton.cc.

References metsig::jet.

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();
 }

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

static

Definition at line 409 of file SoftLepton.cc.

References edm::ParameterSetDescription::add(), and edm::ConfigurationDescriptions::addDefault().

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

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

privatevirtual

Definition at line 96 of file SoftLepton.cc.

References edm::View< T >::begin(), metsig::electron, edm::View< T >::end(), Exception, edm::EventSetup::get(), mps_fire::i, triggerObjects_cff::id, muon::isGoodMuon(), edm::HandleBase::isValid(), fwrapper::jets, softMuonTagInfos_cfi::leptonId, AK4GenJetFlavourInfos_cfi::leptons, m_jets, m_leptonCands, m_leptonId, m_leptons, m_muonSelection, m_transientTrackBuilder, eostools::move(), metsig::muon, edm::errors::NotFound, edm::ESHandle< T >::product(), edm::View< T >::refAt(), mps_fire::result, s_nominalBeamSpot, edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::size(), findQualityFiles::size, tag(), token_electrons, token_gsfElectrons, token_jets, token_jtas, token_leptonCands, token_leptonId, token_muons, token_pfElectrons, token_primaryVertex, token_tracks, and l1t::tracks.

                                                                  {
 
   // grab a TransientTrack helper from the Event Setup
   edm::ESHandle<TransientTrackBuilder> builder;
   setup.get<TransientTrackRecord>().get( "TransientTrackBuilder", builder );
   m_transientTrackBuilder = builder.product();
 
   // 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.getByToken(token_jtas, h_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.getByToken(token_jets, h_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.getByToken(token_primaryVertex, h_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)
     event.getByToken(token_leptonCands, h_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.getByToken(token_gsfElectrons, h_electrons);
 
     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.getByToken(token_electrons, h_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.getByToken(token_pfElectrons, h_electrons);
     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.getByToken(token_muons, h_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.getByToken(token_tracks, h_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())) {
     Handle<edm::ValueMap<float> > h_leptonId;
     event.getByToken(token_leptonId, h_leptonId);
 
     for (Leptons::iterator lepton = leptons.begin(); lepton != leptons.end(); ++lepton)
       lepton->second[leptonId] = (*h_leptonId)[lepton->first];
   }
 
   // output collections
   auto outputCollection = std::make_unique<reco::SoftLeptonTagInfoCollection>();
   for (unsigned int i = 0; i < jets.size(); ++i) {
     reco::SoftLeptonTagInfo result = tag( jets[i], tracks[i], leptons, vertex );
     outputCollection->push_back( result );
   }
   event.put(std::move(outputCollection));
 }

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 315 of file SoftLepton.cc.

References AXIS_CALORIMETRIC_NOLEPTON, AXIS_CHARGED_AVERAGE, AXIS_CHARGED_AVERAGE_NOLEPTON, AXIS_CHARGED_SUM, AXIS_CHARGED_SUM_NOLEPTON, edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::end(), reco::TrackBase::eta(), edm::RefToBase< T >::isNonnull(), M_PI, m_refineJetAxis, reco::LeafCandidate::momentum(), reco::TrackBase::momentum(), or, perp(), reco::TrackBase::phi(), reco::TrackBase::pt(), and HiIsolationCommonParameters_cff::track.

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());
 }

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

staticprotected

Definition at line 391 of file SoftLepton.cc.

References dot(), cmsBatch::log, and mag().

Referenced by tag().

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

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

Definition at line 268 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 = m_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

double SoftLepton::m_chi2Cut

private

Definition at line 120 of file SoftLepton.h.

Referenced by tag().

double SoftLepton::m_deltaRCut

private

Definition at line 119 of file SoftLepton.h.

Referenced by tag().

const edm::InputTag SoftLepton::m_jets

private

Definition at line 99 of file SoftLepton.h.

Referenced by produce().

const edm::InputTag SoftLepton::m_leptonCands

private

Definition at line 109 of file SoftLepton.h.

Referenced by produce().

const edm::InputTag SoftLepton::m_leptonId

private

Definition at line 111 of file SoftLepton.h.

Referenced by produce().

const edm::InputTag SoftLepton::m_leptons

private

Definition at line 103 of file SoftLepton.h.

Referenced by produce().

muon::SelectionType SoftLepton::m_muonSelection

private

Definition at line 123 of file SoftLepton.h.

Referenced by produce().

unsigned int SoftLepton::m_refineJetAxis

private

Definition at line 118 of file SoftLepton.h.

Referenced by refineJetAxis().

const TransientTrackBuilder* SoftLepton::m_transientTrackBuilder

private

Definition at line 115 of file SoftLepton.h.

Referenced by produce(), and tag().

const reco::Vertex SoftLepton::s_nominalBeamSpot

staticprivate

Definition at line 126 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 105 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 104 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 101 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 100 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 110 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 112 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 107 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 106 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 102 of file SoftLepton.h.

Referenced by produce().

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

private

Definition at line 108 of file SoftLepton.h.

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

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation