ConeIsolationAlgorithm Class Reference

#include <ConeIsolationAlgorithm.h>

Public Member Functions
	ConeIsolationAlgorithm (const edm::ParameterSet &parameters)
	ConeIsolationAlgorithm ()
std::pair< float,reco::IsolatedTauTagInfo >	tag (const reco::JetTracksAssociationRef &jetTracks, const reco::Vertex &pv)
	~ConeIsolationAlgorithm ()

Static Public Member Functions
static void	fillDescription (edm::ParameterSetDescription &desc)

Private Attributes
double	dZ_vertex
double	isolation_cone
double	m_cutMaxChiSquared
double	m_cutMaxTIP
double	m_cutMinPt
int	m_cutPixelHits
int	m_cutTotalHits
double	matching_cone
int	n_tracks_isolation_ring
double	pt_min_isolation
double	pt_min_leadTrack
double	signal_cone
bool	useFixedSizeCone
bool	useVertexConstrain_
double	variableConeParameter
double	variableMaxCone
double	variableMinCone

Detailed Description

Definition at line 17 of file ConeIsolationAlgorithm.h.

Constructor & Destructor Documentation

ConeIsolationAlgorithm::ConeIsolationAlgorithm

(

const edm::ParameterSet &

parameters

)

Definition at line 8 of file ConeIsolationAlgorithm.cc.

References edm::ParameterSet::getParameter().

{
  //FIXME: use unsigned int where needed
  m_cutPixelHits     = parameters.getParameter<int>("MinimumNumberOfPixelHits");    // not used
  m_cutTotalHits     = parameters.getParameter<int>("MinimumNumberOfHits");
  m_cutMaxTIP        = parameters.getParameter<double>("MaximumTransverseImpactParameter");
  m_cutMinPt         = parameters.getParameter<double>("MinimumTransverseMomentum");
  m_cutMaxChiSquared = parameters.getParameter<double>("MaximumChiSquared");
  dZ_vertex          = parameters.getParameter<double>("DeltaZetTrackVertex");      //  to be modified
  useVertexConstrain_ = parameters.getParameter<bool>("useVertex");

  matching_cone      = parameters.getParameter<double>("MatchingCone");
  signal_cone        = parameters.getParameter<double>("SignalCone");
  isolation_cone     = parameters.getParameter<double>("IsolationCone"); 
  pt_min_isolation   = parameters.getParameter<double>("MinimumTransverseMomentumInIsolationRing"); 
  pt_min_leadTrack   = parameters.getParameter<double>("MinimumTransverseMomentumLeadingTrack"); 
  n_tracks_isolation_ring = parameters.getParameter<int>("MaximumNumberOfTracksIsolationRing"); 
  
  useFixedSizeCone = parameters.getParameter<bool>("UseFixedSizeCone"); 
  variableConeParameter = parameters.getParameter<double>("VariableConeParameter");
  variableMaxCone = parameters.getParameter<double>("VariableMaxCone");
  variableMinCone = parameters.getParameter<double>("VariableMinCone");
}

ConeIsolationAlgorithm::ConeIsolationAlgorithm ( )

default

ConeIsolationAlgorithm::~ConeIsolationAlgorithm ( )

inline

Definition at line 27 of file ConeIsolationAlgorithm.h.

References fillDescription(), impactParameterTagInfos_cfi::jetTracks, MetAnalyzer::pv(), and tag().

{}

Member Function Documentation

void ConeIsolationAlgorithm::fillDescription ( edm::ParameterSetDescription &  desc )

static

Definition at line 32 of file ConeIsolationAlgorithm.cc.

References edm::ParameterSetDescription::add().

Referenced by ConeIsolation::fillDescriptions(), and ~ConeIsolationAlgorithm().

                                                                            {
  desc.add<int>("MinimumNumberOfPixelHits",2);
  desc.add<int>("MinimumNumberOfHits",8);
  desc.add<double>("MaximumTransverseImpactParameter",0.03);
  desc.add<double>("MinimumTransverseMomentum",1.0);
  desc.add<double>("MaximumChiSquared",100.0);
  desc.add<double>("DeltaZetTrackVertex",0.2);
  desc.add<bool>("useVertex",true);
  desc.add<double>("MatchingCone",0.1);
  desc.add<double>("SignalCone",0.07);
  desc.add<double>("IsolationCone",0.45);
  desc.add<double>("MinimumTransverseMomentumInIsolationRing",0.0);
  desc.add<double>("MinimumTransverseMomentumLeadingTrack",6.0);
  desc.add<int>("MaximumNumberOfTracksIsolationRing",0);
  desc.add<bool>("UseFixedSizeCone",true);
  desc.add<double>("VariableConeParameter",3.5);
  desc.add<double>("VariableMaxCone",0.17);
  desc.add<double>("VariableMinCone",0.05);
}

pair< float, IsolatedTauTagInfo > ConeIsolationAlgorithm::tag	(	const reco::JetTracksAssociationRef &	jetTracks,
		const reco::Vertex &	pv
	)

Definition at line 52 of file ConeIsolationAlgorithm.cc.

References reco::IsolatedTauTagInfo::discriminator(), mvaPFMET_leptons_cfi::discriminator, particleFlowDisplacedVertexCandidate_cfi::dxy, allConversions_cfi::dz, createfilelist::int, reco::IsolatedTauTagInfo::leadingSignalTrack(), hpstanc_transforms::max, min(), reco::Vertex::position(), EnergyCorrector::pt, edm::RefVector< C, T, F >::push_back(), reco::IsolatedTauTagInfo::setDiscriminator(), reco::IsolatedTauTagInfo::setLeadingTrack(), HiIsolationCommonParameters_cff::track, l1t::tracks, and reco::Vertex::z().

Referenced by ~ConeIsolationAlgorithm().

{
  const edm::RefVector<reco::TrackCollection> & tracks = jetTracks->second;
  edm::RefVector<reco::TrackCollection> myTracks;

  // Selection of the Tracks
  float z_pv = pv.z();
  for(auto && track : tracks)
  {
    if ( (track)->pt()                                  >  m_cutMinPt                     &&
         (track)->normalizedChi2()                      <  m_cutMaxChiSquared             &&
         fabs((track)->dxy(pv.position()))                            <  m_cutMaxTIP                    &&
         (track)->recHitsSize()                         >= (unsigned int) m_cutTotalHits  &&
         (track)->hitPattern().numberOfValidPixelHits() >= m_cutPixelHits ) 
    {
      if (useVertexConstrain_ && z_pv > -500.) {
        if (fabs((track)->dz(pv.position())) < dZ_vertex)
          myTracks.push_back(track);
      } else
        myTracks.push_back(track);
    }
  }
  IsolatedTauTagInfo resultExtended(myTracks,jetTracks);

  double r_sigCone = signal_cone;
  double energyJet = jetTracks->first->energy();
  if (not useFixedSizeCone) {
    r_sigCone = std::min(variableMaxCone, variableConeParameter / energyJet);
    r_sigCone = std::max((double)r_sigCone, variableMinCone);
  }

  // now I can use it for the discriminator;
  math::XYZVector jetDir(jetTracks->first->px(), jetTracks->first->py(), jetTracks->first->pz());   
  float discriminator = 0.;
  if (useVertexConstrain_) {
    // In this case all the selected tracks comes from the same vertex, so no need to pass the dZ_vertex requirement to the discriminator 
   const TrackRef myLeadTk = resultExtended.leadingSignalTrack( matching_cone, pt_min_leadTrack);
      resultExtended.setLeadingTrack(myLeadTk);
    discriminator = resultExtended.discriminator(jetDir, matching_cone, r_sigCone, isolation_cone, pt_min_leadTrack, pt_min_isolation,  n_tracks_isolation_ring); 
    resultExtended.setDiscriminator(discriminator);
  } else {
    // In this case the dZ_vertex is used to associate the tracks to the Z_imp parameter of the Leading Track
 const TrackRef myLeadTk = resultExtended.leadingSignalTrack( matching_cone, pt_min_leadTrack);
   resultExtended.setLeadingTrack(myLeadTk);    
    discriminator = resultExtended.discriminator(jetDir, matching_cone, r_sigCone, isolation_cone, pt_min_leadTrack, pt_min_isolation,  n_tracks_isolation_ring, dZ_vertex); 
    resultExtended.setDiscriminator(discriminator);
  }

  return std::make_pair(discriminator, resultExtended); 
}

Member Data Documentation

double ConeIsolationAlgorithm::dZ_vertex

private

Definition at line 46 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::isolation_cone

private

Definition at line 43 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::m_cutMaxChiSquared

private

Definition at line 40 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::m_cutMaxTIP

private

Definition at line 38 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::m_cutMinPt

private

Definition at line 39 of file ConeIsolationAlgorithm.h.

int ConeIsolationAlgorithm::m_cutPixelHits

private

Definition at line 36 of file ConeIsolationAlgorithm.h.

int ConeIsolationAlgorithm::m_cutTotalHits

private

Definition at line 37 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::matching_cone

private

Definition at line 41 of file ConeIsolationAlgorithm.h.

int ConeIsolationAlgorithm::n_tracks_isolation_ring

private

Definition at line 47 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::pt_min_isolation

private

Definition at line 44 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::pt_min_leadTrack

private

Definition at line 45 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::signal_cone

private

Definition at line 42 of file ConeIsolationAlgorithm.h.

bool ConeIsolationAlgorithm::useFixedSizeCone

private

Definition at line 49 of file ConeIsolationAlgorithm.h.

bool ConeIsolationAlgorithm::useVertexConstrain_

private

Definition at line 48 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::variableConeParameter

private

Definition at line 50 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::variableMaxCone

private

Definition at line 51 of file ConeIsolationAlgorithm.h.

double ConeIsolationAlgorithm::variableMinCone

private

Definition at line 52 of file ConeIsolationAlgorithm.h.