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) const
	~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() [1/2]

ConeIsolationAlgorithm::ConeIsolationAlgorithm

(

const edm::ParameterSet &

parameters

)

Definition at line 8 of file ConeIsolationAlgorithm.cc.

                                                                             {
  //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() [2/2]

ConeIsolationAlgorithm::ConeIsolationAlgorithm ( )

default

~ConeIsolationAlgorithm()

ConeIsolationAlgorithm::~ConeIsolationAlgorithm ( )

inline

Definition at line 26 of file ConeIsolationAlgorithm.h.

{}

Member Function Documentation

fillDescription()

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

static

Definition at line 31 of file ConeIsolationAlgorithm.cc.

References submitPVResolutionJobs::desc.

Referenced by ConeIsolation::fillDescriptions().

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

tag()

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

Definition at line 51 of file ConeIsolationAlgorithm.cc.

References PDWG_TauSkim_cff::discriminator, reco::IsolatedTauTagInfo::discriminator(), PVValHelper::dxy, PVValHelper::dz, createfilelist::int, HLT_2023v12_cff::jetTracks, reco::IsolatedTauTagInfo::leadingSignalTrack(), SiStripPI::max, SiStripPI::min, TrackCollections2monitor_cff::normalizedChi2, DiDispStaMuonMonitor_cfi::pt, edm::RefVector< C, T, F >::push_back(), reco::IsolatedTauTagInfo::setDiscriminator(), reco::IsolatedTauTagInfo::setLeadingTrack(), HLT_2023v12_cff::track, and pwdgSkimBPark_cfi::tracks.

Referenced by ConeIsolation::produce().

                                                                                    {
  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

dZ_vertex

double ConeIsolationAlgorithm::dZ_vertex

private

Definition at line 45 of file ConeIsolationAlgorithm.h.

isolation_cone

double ConeIsolationAlgorithm::isolation_cone

private

Definition at line 42 of file ConeIsolationAlgorithm.h.

m_cutMaxChiSquared

double ConeIsolationAlgorithm::m_cutMaxChiSquared

private

Definition at line 39 of file ConeIsolationAlgorithm.h.

m_cutMaxTIP

double ConeIsolationAlgorithm::m_cutMaxTIP

private

Definition at line 37 of file ConeIsolationAlgorithm.h.

m_cutMinPt

double ConeIsolationAlgorithm::m_cutMinPt

private

Definition at line 38 of file ConeIsolationAlgorithm.h.

m_cutPixelHits

int ConeIsolationAlgorithm::m_cutPixelHits

private

Definition at line 35 of file ConeIsolationAlgorithm.h.

m_cutTotalHits

int ConeIsolationAlgorithm::m_cutTotalHits

private

Definition at line 36 of file ConeIsolationAlgorithm.h.

matching_cone

double ConeIsolationAlgorithm::matching_cone

private

Definition at line 40 of file ConeIsolationAlgorithm.h.

n_tracks_isolation_ring

int ConeIsolationAlgorithm::n_tracks_isolation_ring

private

Definition at line 46 of file ConeIsolationAlgorithm.h.

pt_min_isolation

double ConeIsolationAlgorithm::pt_min_isolation

private

Definition at line 43 of file ConeIsolationAlgorithm.h.

pt_min_leadTrack

double ConeIsolationAlgorithm::pt_min_leadTrack

private

Definition at line 44 of file ConeIsolationAlgorithm.h.

signal_cone

double ConeIsolationAlgorithm::signal_cone

private

Definition at line 41 of file ConeIsolationAlgorithm.h.

useFixedSizeCone

bool ConeIsolationAlgorithm::useFixedSizeCone

private

Definition at line 48 of file ConeIsolationAlgorithm.h.

useVertexConstrain_

bool ConeIsolationAlgorithm::useVertexConstrain_

private

Definition at line 47 of file ConeIsolationAlgorithm.h.

variableConeParameter

double ConeIsolationAlgorithm::variableConeParameter

private

Definition at line 49 of file ConeIsolationAlgorithm.h.

variableMaxCone

double ConeIsolationAlgorithm::variableMaxCone

private

Definition at line 50 of file ConeIsolationAlgorithm.h.

variableMinCone

double ConeIsolationAlgorithm::variableMinCone

private

Definition at line 51 of file ConeIsolationAlgorithm.h.