#include <PromptTrackCountingComputer.h>

Inheritance diagram for PromptTrackCountingComputer:

Public Types
using	Tokens = void

Public Member Functions
float	discriminator (const TagInfoHelper &ti) const override

	PromptTrackCountingComputer (const edm::ParameterSet &parameters)

Public Member Functions inherited from JetTagComputer
const std::vector< std::string > &	getInputLabels () const

virtual void	initialize (const JetTagComputerRecord &)

	JetTagComputer ()

	JetTagComputer (const edm::ParameterSet &configuration)

float	operator() (const reco::BaseTagInfo &info) const

float	operator() (const TagInfoHelper &helper) const

void	setupDone ()

virtual	~JetTagComputer ()

Protected Member Functions
std::multiset< float >	orderedSignificances (const reco::TrackIPTagInfo &tkip) const

Protected Member Functions inherited from JetTagComputer
virtual float	discriminator (const reco::BaseTagInfo &) const

void	uses (unsigned int id, const std::string &label)

void	uses (const std::string &label)

Protected Attributes
double	m_cutMaxDecayLen

double	m_cutMaxDistToAxis

double	m_deltaR

double	m_deltaRmin

int	m_ipType

int	m_nthTrack

reco::TrackBase::TrackQuality	m_trackQuality

bool	m_useAllQualities

double	maxImpactParameter

double	maxImpactParameterSig

Detailed Description

Definition at line 16 of file PromptTrackCountingComputer.h.

Member Typedef Documentation

◆ Tokens

using PromptTrackCountingComputer::Tokens = void

Definition at line 18 of file PromptTrackCountingComputer.h.

Constructor & Destructor Documentation

◆ PromptTrackCountingComputer()

PromptTrackCountingComputer::PromptTrackCountingComputer ( const edm::ParameterSet & parameters )

inline

Definition at line 20 of file PromptTrackCountingComputer.h.

References m_cutMaxDecayLen, m_cutMaxDistToAxis, m_deltaR, m_deltaRmin, m_ipType, m_nthTrack, m_trackQuality, m_useAllQualities, maxImpactParameter, maxImpactParameterSig, reco::TrackBase::qualityByName(), AlCaHLTBitMon_QueryRunRegistry::string, and JetTagComputer::uses().

                                                                  {
     m_nthTrack = parameters.getParameter<int>("nthTrack");
     m_ipType = parameters.getParameter<int>("impactParameterType");
     // Maximum and minimum allowed deltaR respectively.
     m_deltaR = parameters.getParameter<double>("deltaR");
     m_deltaRmin = parameters.getParameter<double>("deltaRmin");
     maxImpactParameter = parameters.getParameter<double>("maxImpactParameter");
     maxImpactParameterSig = parameters.getParameter<double>("maxImpactParameterSig");
     m_cutMaxDecayLen = parameters.getParameter<double>("maximumDecayLength");          //used
     m_cutMaxDistToAxis = parameters.getParameter<double>("maximumDistanceToJetAxis");  //used
     //
     // access track quality class; "any" takes everything
     //
     std::string trackQualityType = parameters.getParameter<std::string>("trackQualityClass");  //used
     m_trackQuality = reco::TrackBase::qualityByName(trackQualityType);
     m_useAllQualities = false;
     if (trackQualityType == "any" || trackQualityType == "Any" || trackQualityType == "ANY")
       m_useAllQualities = true;
 
     uses("ipTagInfos");
   }

Member Function Documentation

◆ discriminator()

float PromptTrackCountingComputer::discriminator ( const TagInfoHelper & ti ) const

inlineoverridevirtual

Reimplemented from JetTagComputer.

Definition at line 42 of file PromptTrackCountingComputer.h.

References JetTagComputer::TagInfoHelper::get(), maxImpactParameterSig, and orderedSignificances().

                                                               {
     const reco::TrackIPTagInfo& tkip = ti.get<reco::TrackIPTagInfo>();
     std::multiset<float> significances = orderedSignificances(tkip);
     std::multiset<float>::iterator sig;
     unsigned int nPromptTrk = 0;
     for (sig = significances.begin(); sig != significances.end(); sig++) {
       if (fabs(*sig) < maxImpactParameterSig)
         nPromptTrk++;
       //       edm::LogDebug("") << "Track "<< nPromptTrk << " sig=" << *sig;
     }
     return double(nPromptTrk);
   }

◆ orderedSignificances()

std::multiset<float> PromptTrackCountingComputer::orderedSignificances ( const reco::TrackIPTagInfo & tkip ) const

inlineprotected

Definition at line 56 of file PromptTrackCountingComputer.h.

References electronAnalyzer_cfi::DeltaR, mps_fire::i, reco::IPTagInfo< Container, Base >::impactParameterData(), edm::Ref< C, T, F >::isNull(), m_cutMaxDecayLen, m_cutMaxDistToAxis, m_deltaR, m_deltaRmin, m_ipType, m_trackQuality, m_useAllQualities, maxImpactParameter, reco::IPTagInfo< Container, Base >::primaryVertex(), MetAnalyzer::pv(), quality, reco::IPTagInfo< Container, Base >::selectedTracks(), pwdgSkimBPark_cfi::tracks, and relativeConstraints::value.

Referenced by discriminator().

                                                                                 {
     const std::vector<reco::btag::TrackIPData>& impactParameters((tkip.impactParameterData()));
     const edm::RefVector<reco::TrackCollection>& tracks(tkip.selectedTracks());
     std::multiset<float> significances;
     int i = 0;
     if (tkip.primaryVertex().isNull()) {
       return std::multiset<float>();
     }
 
     GlobalPoint pv(tkip.primaryVertex()->position().x(),
                    tkip.primaryVertex()->position().y(),
                    tkip.primaryVertex()->position().z());
 
     for (std::vector<reco::btag::TrackIPData>::const_iterator it = impactParameters.begin();
          it != impactParameters.end();
          ++it, i++) {
       if (fabs(impactParameters[i].distanceToJetAxis.value()) < m_cutMaxDistToAxis &&  // distance to JetAxis
           (impactParameters[i].closestToJetAxis - pv).mag() < m_cutMaxDecayLen &&      // max decay len
           (m_useAllQualities == true || (*tracks[i]).quality(m_trackQuality))          // use selected track qualities
       ) {
         if ((m_deltaR <= 0 ||
              ROOT::Math::VectorUtil::DeltaR((*tkip.jet()).p4().Vect(), (*tracks[i]).momentum()) < m_deltaR) &&
             (m_deltaRmin <= 0 ||
              ROOT::Math::VectorUtil::DeltaR((*tkip.jet()).p4().Vect(), (*tracks[i]).momentum()) > m_deltaRmin)) {
           if (fabs(((m_ipType == 0) ? it->ip3d : it->ip2d).value()) < maxImpactParameter) {
             significances.insert(((m_ipType == 0) ? it->ip3d : it->ip2d).significance());
           }
         }
       }
     }
 
     return significances;
   }