reco::IPTagInfo< Container, Base > Class Template Reference

#include <IPTagInfo.h>

Inheritance diagram for reco::IPTagInfo< Container, Base >:

Public Types
typedef Base	base_class
typedef Container	input_container

Public Member Functions
const GlobalVector &	axis () const
virtual IPTagInfo *	clone (void) const
	clone More...
const TrackRef &	ghostTrack () const
virtual bool	hasProbabilities () const
const std::vector < btag::TrackIPData > &	impactParameterData () const
	IPTagInfo (const std::vector< btag::TrackIPData > &ipData, const std::vector< float > &prob2d, const std::vector< float > &prob3d, const Container &selected, const Base &base, const edm::Ref< VertexCollection > &pv, const GlobalVector &axis, const TrackRef &ghostTrack)
	IPTagInfo ()
const edm::Ref < VertexCollection > &	primaryVertex () const
const std::vector< float > &	probabilities (int ip) const
const Container &	selected () const
const Track *	selectedTrack (size_t i) const
const Container &	selectedTracks () const
Container	sorted (const std::vector< size_t > &indexes) const
std::vector< size_t >	sortedIndexes (btag::SortCriteria mode=reco::btag::IP3DSig) const
std::vector< size_t >	sortedIndexesWithCut (float cut, btag::SortCriteria mode=reco::btag::IP3DSig) const
Container	sortedTracks (const std::vector< size_t > &indexes) const
virtual TaggingVariableList	taggingVariables (void) const
std::vector< bool >	variableJTA (const btag::variableJTAParameters &params) const
virtual	~IPTagInfo ()

Static Public Member Functions
static short	Class_Version ()
static bool	passVariableJTA (const btag::variableJTAParameters &params, double jetpt, double trackpt, double jettrackdr)

Private Attributes
GlobalVector	m_axis
std::vector< btag::TrackIPData >	m_data
TrackRef	m_ghostTrack
std::vector< float >	m_prob2d
std::vector< float >	m_prob3d
edm::Ref< VertexCollection >	m_pv
Container	m_selected

Detailed Description

template<class Container, class Base>
class reco::IPTagInfo< Container, Base >

Definition at line 51 of file IPTagInfo.h.

Member Typedef Documentation

template<class Container, class Base>

typedef Base reco::IPTagInfo< Container, Base >::base_class

Definition at line 54 of file IPTagInfo.h.

template<class Container, class Base>

typedef Container reco::IPTagInfo< Container, Base >::input_container

Definition at line 53 of file IPTagInfo.h.

Constructor & Destructor Documentation

template<class Container, class Base>

reco::IPTagInfo< Container, Base >::IPTagInfo ( const std::vector< btag::TrackIPData > &  ipData, const std::vector< float > &  prob2d, const std::vector< float > &  prob3d, const Container &  selected, const Base &  base, const edm::Ref< VertexCollection > &  pv, const GlobalVector &  axis, const TrackRef &  ghostTrack  )

inline

Definition at line 56 of file IPTagInfo.h.

                                 :
      Base(base), m_data(ipData),  m_prob2d(prob2d),
      m_prob3d(prob3d), m_selected(selected), m_pv(pv),
      m_axis(axis), m_ghostTrack(ghostTrack) {}

template<class Container, class Base>

reco::IPTagInfo< Container, Base >::IPTagInfo ( )

inline

Definition at line 69 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::clone().

{}

template<class Container, class Base>

virtual reco::IPTagInfo< Container, Base >::~IPTagInfo ( )

inlinevirtual

Definition at line 71 of file IPTagInfo.h.

{}

Member Function Documentation

template<class Container, class Base>

const GlobalVector& reco::IPTagInfo< Container, Base >::axis ( ) const

inline

Definition at line 135 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_axis.

{ return m_axis; }

template<class Container, class Base>

static short reco::IPTagInfo< Container, Base >::Class_Version ( )

inlinestatic

Definition at line 141 of file IPTagInfo.h.

:
  std::vector<btag::TrackIPData> m_data;

template<class Container, class Base>

virtual IPTagInfo* reco::IPTagInfo< Container, Base >::clone ( void  ) const

inlinevirtual

clone

Definition at line 74 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::IPTagInfo().

  { return new IPTagInfo(*this); }

template<class Container, class Base>

const TrackRef& reco::IPTagInfo< Container, Base >::ghostTrack ( ) const

inline

Definition at line 136 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_ghostTrack.

{ return m_ghostTrack; }

template<class Container, class Base>

virtual bool reco::IPTagInfo< Container, Base >::hasProbabilities ( ) const

inlinevirtual

Check if probability information is globally available impact parameters in the collection

Even if true for some tracks it is possible that a -1 probability is returned if some problem occured

Definition at line 85 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_data, and reco::IPTagInfo< Container, Base >::m_prob3d.

Referenced by IPTagPlotter< Container, Base >::analyzeTag(), and reco::IPTagInfo< Container, Base >::sortedIndexesWithCut().

  { return m_data.size() == m_prob3d.size(); }

template<class Container, class Base>

const std::vector<btag::TrackIPData>& reco::IPTagInfo< Container, Base >::impactParameterData ( ) const

inline

Vectors of TrackIPData orderd as the selectedTracks()

Definition at line 91 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_data.

Referenced by IPTagPlotter< Container, Base >::analyzeTag(), CandidateBoostedDoubleSecondaryVertexComputer::discriminator(), TemplatedJetProbabilityComputer< Container, Base >::discriminator(), TemplatedJetBProbabilityComputer< Container, Base >::discriminator(), GhostTrackComputer::operator()(), PromptTrackCountingComputer::orderedSignificances(), TemplatedTrackCountingComputer< Container, Base >::orderedSignificances(), and GhostTrackComputer::threshTrack().

  { return m_data; }

template<class Container , class Base >

bool reco::IPTagInfo< Container, Base >::passVariableJTA ( const btag::variableJTAParameters &  params, double  jetpt, double  trackpt, double  jettrackdr  )

static

Definition at line 303 of file IPTagInfo.h.

References reco::btag::variableJTAParameters::a_dR, reco::btag::variableJTAParameters::a_pT, reco::btag::variableJTAParameters::b_dR, reco::btag::variableJTAParameters::b_pT, reco::btag::variableJTAParameters::max_pT, reco::btag::variableJTAParameters::max_pT_dRcut, reco::btag::variableJTAParameters::max_pT_trackPTcut, reco::btag::variableJTAParameters::min_pT, and reco::btag::variableJTAParameters::min_pT_dRcut.

Referenced by reco::TrackSelector::operator()(), and reco::IPTagInfo< Container, Base >::variableJTA().

                                                                                                                                                                              {

  bool pass = false;

  // intermediate pt range (between min_pT and max_pT), apply variable JTA !
  if ( jetpT > params.min_pT && jetpT < params.max_pT ) {
    double deltaRfunction_highpt = -jetpT * params.a_dR + params.b_dR;
    double ptfunction_highpt = jetpT * params.a_pT + params.b_pT;
    
    if (jettrackdR < deltaRfunction_highpt
    &&
    trackpT > ptfunction_highpt) 
      pass = true;
    
    //  cout << "IPTagInfo: passVariableJTA: dR and TrackpT " << jettrackdR << " " << trackpT << endl;
    
    //high pt range, apply fixed default cuts
  }else if (jetpT > params.max_pT ) {
    if (jettrackdR < params.max_pT_dRcut
    &&
    trackpT > params.max_pT_trackPTcut)
      pass = true;

    // low pt range, apply fixed default cuts
  }else {
    if (jettrackdR < params.min_pT_dRcut)
      pass = true;
  }
  
  return pass;
}

template<class Container, class Base>

const edm::Ref<VertexCollection>& reco::IPTagInfo< Container, Base >::primaryVertex ( ) const

inline

Definition at line 133 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_pv.

Referenced by IPTagPlotter< Container, Base >::analyzeTag(), CandidateBoostedDoubleSecondaryVertexComputer::discriminator(), TemplatedJetProbabilityComputer< Container, Base >::discriminator(), TemplatedJetBProbabilityComputer< Container, Base >::discriminator(), GhostTrackComputer::operator()(), PromptTrackCountingComputer::orderedSignificances(), TemplatedTrackCountingComputer< Container, Base >::orderedSignificances(), and reco::IPTagInfo< Container, Base >::taggingVariables().

{ return m_pv; }

template<class Container, class Base>

const std::vector<float>& reco::IPTagInfo< Container, Base >::probabilities ( int  ip ) const

inline

Definition at line 103 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_prob2d, and reco::IPTagInfo< Container, Base >::m_prob3d.

Referenced by IPTagPlotter< Container, Base >::analyzeTag(), TemplatedJetProbabilityComputer< Container, Base >::discriminator(), and TemplatedJetBProbabilityComputer< Container, Base >::discriminator().

{return (ip==0)?m_prob3d:m_prob2d; }

template<class Container, class Base>

const Container& reco::IPTagInfo< Container, Base >::selected ( ) const

inline

Return the vector of tracks for which the IP information is available Quality cuts are applied to reject fake tracks

Definition at line 98 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_selected.

{ return m_selected; }

template<class Container, class Base>

const Track* reco::IPTagInfo< Container, Base >::selectedTrack ( size_t  i ) const

inline

Definition at line 138 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_selected, and reco::btag::toTrack().

Referenced by reco::IPTagInfo< Container, Base >::taggingVariables(), and reco::IPTagInfo< Container, Base >::variableJTA().

{return reco::btag::toTrack(m_selected[i]);}

template<class Container, class Base>

const Container& reco::IPTagInfo< Container, Base >::selectedTracks ( ) const

inline

Definition at line 101 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::m_selected.

Referenced by CandidateBoostedDoubleSecondaryVertexComputer::discriminator(), TemplatedJetProbabilityComputer< Container, Base >::discriminator(), TemplatedJetBProbabilityComputer< Container, Base >::discriminator(), GhostTrackComputer::operator()(), PromptTrackCountingComputer::orderedSignificances(), TemplatedTrackCountingComputer< Container, Base >::orderedSignificances(), and GhostTrackComputer::threshTrack().

{ return m_selected; }

template<class Container , class Base >

Container reco::IPTagInfo< Container, Base >::sorted

(

const std::vector< size_t > &

indexes

)

const

Definition at line 205 of file IPTagInfo.h.

References i, and reco::IPTagInfo< Container, Base >::m_selected.

Referenced by reco::IPTagInfo< Container, Base >::sortedTracks().

{
 Container tr;
 for(size_t i =0 ; i < indexes.size(); i++) tr.push_back(m_selected[indexes[i]]);
 return tr;
}

template<class Container , class Base >

std::vector< size_t > reco::IPTagInfo< Container, Base >::sortedIndexes

(

btag::SortCriteria

mode = reco::btag::IP3DSig

)

const

Return the list of track index sorted by mode

Definition at line 238 of file IPTagInfo.h.

References BTagSF::btag, GOODCOLL_filter_cfg::cut, reco::btag::Prob2D, reco::btag::Prob3D, and reco::IPTagInfo< Container, Base >::sortedIndexesWithCut().

Referenced by IPTagPlotter< Container, Base >::analyzeTag(), CandidateBoostedDoubleSecondaryVertexComputer::discriminator(), GhostTrackComputer::operator()(), reco::IPTagInfo< Container, Base >::taggingVariables(), and GhostTrackComputer::threshTrack().

{
 using namespace reco::btag;
 float cut=-1e99;
 if((mode == Prob3D || mode == Prob2D)) cut=1e99;
 return sortedIndexesWithCut(cut,mode);
}

template<class Container , class Base >

std::vector< size_t > reco::IPTagInfo< Container, Base >::sortedIndexesWithCut	(	float	cut,
		btag::SortCriteria	mode = reco::btag::IP3DSig
	)		const

Return the list of track index sorted by mode A cut can is specified to select only tracks with IP value or significance > cut or probability < cut (according to the specified mode)

Definition at line 246 of file IPTagInfo.h.

References BTagSF::btag, reco::IPTagInfo< Container, Base >::hasProbabilities(), i, reco::btag::IP2DSig, reco::btag::IP2DValue, reco::btag::IP3DSig, reco::btag::IP3DValue, reco::IPTagInfo< Container, Base >::m_data, reco::IPTagInfo< Container, Base >::m_prob2d, reco::IPTagInfo< Container, Base >::m_prob3d, reco::IPTagInfo< Container, Base >::m_selected, reco::btag::Prob2D, reco::btag::Prob3D, and query::result.

Referenced by reco::IPTagInfo< Container, Base >::sortedIndexes().

{
 std::multimap<float,size_t> sortedIdx;
 size_t nSelectedTracks = m_selected.size();
 std::vector<size_t> result;
 using namespace reco::btag;
 
//check if probabilities are available
 if((mode == Prob3D || mode == Prob2D) && ! hasProbabilities()) 
  {
   return result;
  }

 for(size_t i=0;i<nSelectedTracks;i++) 
  {
     float sortingKey;
     switch(mode)
     {
      case IP3DSig:
           sortingKey=m_data[i].ip3d.significance();
           break;
      case IP2DSig:
           sortingKey=m_data[i].ip2d.significance();
           break;
      case IP3DValue:
           sortingKey=m_data[i].ip3d.value();
           break;
      case IP2DValue:
           sortingKey=m_data[i].ip2d.value();
           break;
      case Prob3D:
           sortingKey=m_prob3d[i];
           break;
      case Prob2D:
           sortingKey=m_prob2d[i];
           break;

      default:
       sortingKey=i;
     }   
     sortedIdx.insert(std::pair<float,size_t>(sortingKey,i));
  }

//Descending: 
if(mode == IP3DSig || mode == IP2DSig ||mode ==  IP3DValue || mode == IP2DValue)
 { 
   for(std::multimap<float,size_t>::reverse_iterator it = sortedIdx.rbegin(); it!=sortedIdx.rend(); it++)
    if(it->first >= cut) result.push_back(it->second);
 } else
//Ascending:
 {
  for(std::multimap<float,size_t>::iterator it = sortedIdx.begin(); it!=sortedIdx.end(); it++)
    if(it->first <= cut) result.push_back(it->second);
 }
 return result;
}

template<class Container, class Base>

Container reco::IPTagInfo< Container, Base >::sortedTracks ( const std::vector< size_t > &  indexes ) const

inline

Definition at line 129 of file IPTagInfo.h.

References reco::IPTagInfo< Container, Base >::sorted().

Referenced by IPTagPlotter< Container, Base >::analyzeTag().

{return sorted(indexes);}

template<class Container , class Base >

TaggingVariableList reco::IPTagInfo< Container, Base >::taggingVariables ( void  ) const

virtual

Definition at line 158 of file IPTagInfo.h.

References data, HLT_25ns14e33_v1_cff::DeltaR, reco::btau::etaRel(), reco::TaggingVariableList::finalize(), reco::TaggingVariableList::insert(), edm::Ref< C, T, F >::isNonnull(), metsig::jet, reco::IPTagInfo< Container, Base >::primaryVertex(), MetAnalyzer::pv(), reco::IPTagInfo< Container, Base >::selectedTrack(), reco::IPTagInfo< Container, Base >::sortedIndexes(), mathSSE::sqrt(), reco::btau::trackChi2, reco::btau::trackDecayLenVal, reco::btau::trackDeltaR, reco::btau::trackEta, reco::btau::trackEtaRel, reco::btau::trackGhostTrackDistSig, reco::btau::trackGhostTrackDistVal, reco::btau::trackGhostTrackWeight, reco::btau::trackJetDistSig, reco::btau::trackJetDistVal, reco::btau::trackMomentum, reco::btau::trackNPixelHits, reco::btau::trackNTotalHits, reco::btau::trackPPar, reco::btau::trackPParRatio, reco::btau::trackPtRatio, reco::btau::trackPtRel, reco::btau::trackSip2dSig, reco::btau::trackSip2dVal, reco::btau::trackSip3dSig, reco::btau::trackSip3dVal, and x().

                                                                                                                 {
  TaggingVariableList vars;

  math::XYZVector jetDir = Base::jet()->momentum().Unit();
  bool havePv = primaryVertex().isNonnull();
  GlobalPoint pv;
  if (havePv)
    pv = GlobalPoint(primaryVertex()->x(),
                     primaryVertex()->y(),
                     primaryVertex()->z());

  std::vector<size_t> indexes = sortedIndexes(); // use default criterium
  for(std::vector<size_t>::const_iterator it = indexes.begin();
      it != indexes.end(); ++it)
   {
     using namespace ROOT::Math;
     const Track * track = selectedTrack(*it);
     const btag::TrackIPData *data = &m_data[*it];
     math::XYZVector trackMom = track->momentum();
     double trackMag = std::sqrt(trackMom.Mag2());

     vars.insert(btau::trackMomentum, trackMag, true);
     vars.insert(btau::trackEta, trackMom.Eta(), true);
     vars.insert(btau::trackEtaRel, reco::btau::etaRel(jetDir, trackMom), true);
     vars.insert(btau::trackPtRel, VectorUtil::Perp(trackMom, jetDir), true);
     vars.insert(btau::trackPPar, jetDir.Dot(trackMom), true);
     vars.insert(btau::trackDeltaR, VectorUtil::DeltaR(trackMom, jetDir), true);
     vars.insert(btau::trackPtRatio, VectorUtil::Perp(trackMom, jetDir) / trackMag, true);
     vars.insert(btau::trackPParRatio, jetDir.Dot(trackMom) / trackMag, true);
     vars.insert(btau::trackSip3dVal, data->ip3d.value(), true);
     vars.insert(btau::trackSip3dSig, data->ip3d.significance(), true);
     vars.insert(btau::trackSip2dVal, data->ip2d.value(), true);
     vars.insert(btau::trackSip2dSig, data->ip2d.significance(), true);
     vars.insert(btau::trackDecayLenVal, havePv ? (data->closestToJetAxis - pv).mag() : -1.0, true);
     vars.insert(btau::trackJetDistVal, data->distanceToJetAxis.value(), true);
     vars.insert(btau::trackJetDistSig, data->distanceToJetAxis.significance(), true);
     vars.insert(btau::trackGhostTrackDistVal, data->distanceToGhostTrack.value(), true);
     vars.insert(btau::trackGhostTrackDistSig, data->distanceToGhostTrack.significance(), true);
     vars.insert(btau::trackGhostTrackWeight, data->ghostTrackWeight, true);
     vars.insert(btau::trackChi2, track->normalizedChi2(), true);
     vars.insert(btau::trackNTotalHits, track->hitPattern().numberOfValidHits(), true);
     vars.insert(btau::trackNPixelHits, track->hitPattern().numberOfValidPixelHits(), true);
   } 
  vars.finalize();
  return vars;
}

template<class Container , class Base >

std::vector< bool > reco::IPTagInfo< Container, Base >::variableJTA

(

const btag::variableJTAParameters &

params

)

const

variable jet-to track association: returns vector of bool, indicating for each track whether it passed the variable JTA.

Definition at line 212 of file IPTagInfo.h.

References HLT_25ns14e33_v1_cff::DeltaR, i, metsig::jet, reco::IPTagInfo< Container, Base >::m_selected, reco::TrackBase::momentum(), reco::IPTagInfo< Container, Base >::passVariableJTA(), reco::TrackBase::pt(), query::result, and reco::IPTagInfo< Container, Base >::selectedTrack().

Referenced by TemplatedJetProbabilityComputer< Container, Base >::discriminator(), TemplatedJetBProbabilityComputer< Container, Base >::discriminator(), and TemplatedTrackCountingComputer< Container, Base >::orderedSignificances().

                                                                                                                                         {
  
  std::vector<bool> result;

  //Jet parameters
  double jetpT = Base::jet()->pt();
  math::XYZVector jetDir = Base::jet()->momentum().Unit();

  for(size_t  i = 0 ; i<  m_selected.size();  i++) {
    
    //Track parameters
    const Track * track =  selectedTrack(i);    
    double trackpT = track->pt();
    math::XYZVector trackMom = track->momentum();

    // do the math in passVariableJTA
    result.push_back(passVariableJTA( params, jetpT, trackpT, ROOT::Math::VectorUtil::DeltaR(trackMom, jetDir)));

  }  
  
  return result;
}

Member Data Documentation

template<class Container, class Base>

GlobalVector reco::IPTagInfo< Container, Base >::m_axis

private

Definition at line 149 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::axis().

template<class Container, class Base>

std::vector<btag::TrackIPData> reco::IPTagInfo< Container, Base >::m_data

private

Definition at line 144 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::hasProbabilities(), reco::IPTagInfo< Container, Base >::impactParameterData(), and reco::IPTagInfo< Container, Base >::sortedIndexesWithCut().

template<class Container, class Base>

TrackRef reco::IPTagInfo< Container, Base >::m_ghostTrack

private

Definition at line 150 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::ghostTrack().

template<class Container, class Base>

std::vector<float> reco::IPTagInfo< Container, Base >::m_prob2d

private

Definition at line 145 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::probabilities(), and reco::IPTagInfo< Container, Base >::sortedIndexesWithCut().

template<class Container, class Base>

std::vector<float> reco::IPTagInfo< Container, Base >::m_prob3d

private

Definition at line 146 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::hasProbabilities(), reco::IPTagInfo< Container, Base >::probabilities(), and reco::IPTagInfo< Container, Base >::sortedIndexesWithCut().

template<class Container, class Base>

edm::Ref<VertexCollection> reco::IPTagInfo< Container, Base >::m_pv

private

Definition at line 148 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::primaryVertex().

template<class Container, class Base>

Container reco::IPTagInfo< Container, Base >::m_selected

private

Definition at line 147 of file IPTagInfo.h.

Referenced by reco::IPTagInfo< Container, Base >::selected(), reco::IPTagInfo< Container, Base >::selectedTrack(), reco::IPTagInfo< Container, Base >::selectedTracks(), reco::IPTagInfo< Container, Base >::sorted(), reco::IPTagInfo< Container, Base >::sortedIndexesWithCut(), and reco::IPTagInfo< Container, Base >::variableJTA().