CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
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 9 of file ConeIsolationAlgorithm.cc.

References edm::ParameterSet::getParameter().

10 {
11  //FIXME: use unsigned int where needed
12  m_cutPixelHits = parameters.getParameter<int>("MinimumNumberOfPixelHits"); // not used
13  m_cutTotalHits = parameters.getParameter<int>("MinimumNumberOfHits");
14  m_cutMaxTIP = parameters.getParameter<double>("MaximumTransverseImpactParameter");
15  m_cutMinPt = parameters.getParameter<double>("MinimumTransverseMomentum");
16  m_cutMaxChiSquared = parameters.getParameter<double>("MaximumChiSquared");
17  dZ_vertex = parameters.getParameter<double>("DeltaZetTrackVertex"); // to be modified
18  useVertexConstrain_ = parameters.getParameter<bool>("useVertex");
19 
20  matching_cone = parameters.getParameter<double>("MatchingCone");
21  signal_cone = parameters.getParameter<double>("SignalCone");
22  isolation_cone = parameters.getParameter<double>("IsolationCone");
23  pt_min_isolation = parameters.getParameter<double>("MinimumTransverseMomentumInIsolationRing");
24  pt_min_leadTrack = parameters.getParameter<double>("MinimumTransverseMomentumLeadingTrack");
25  n_tracks_isolation_ring = parameters.getParameter<int>("MaximumNumberOfTracksIsolationRing");
26 
27  useFixedSizeCone = parameters.getParameter<bool>("UseFixedSizeCone");
28  variableConeParameter = parameters.getParameter<double>("VariableConeParameter");
29  variableMaxCone = parameters.getParameter<double>("VariableMaxCone");
30  variableMinCone = parameters.getParameter<double>("VariableMinCone");
31 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
ConeIsolationAlgorithm::ConeIsolationAlgorithm ( void  )

Definition at line 6 of file ConeIsolationAlgorithm.cc.

7 { }
ConeIsolationAlgorithm::~ConeIsolationAlgorithm ( )
inline

Definition at line 27 of file ConeIsolationAlgorithm.h.

27 {}

Member Function Documentation

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

Definition at line 33 of file ConeIsolationAlgorithm.cc.

References edm::ParameterSetDescription::add().

Referenced by ConeIsolation::fillDescriptions().

33  {
34  desc.add<int>("MinimumNumberOfPixelHits",2);
35  desc.add<int>("MinimumNumberOfHits",8);
36  desc.add<double>("MaximumTransverseImpactParameter",0.03);
37  desc.add<double>("MinimumTransverseMomentum",1.0);
38  desc.add<double>("MaximumChiSquared",100.0);
39  desc.add<double>("DeltaZetTrackVertex",0.2);
40  desc.add<bool>("useVertex",true);
41  desc.add<double>("MatchingCone",0.1);
42  desc.add<double>("SignalCone",0.07);
43  desc.add<double>("IsolationCone",0.45);
44  desc.add<double>("MinimumTransverseMomentumInIsolationRing",0.0);
45  desc.add<double>("MinimumTransverseMomentumLeadingTrack",6.0);
46  desc.add<int>("MaximumNumberOfTracksIsolationRing",0);
47  desc.add<bool>("UseFixedSizeCone",true);
48  desc.add<double>("VariableConeParameter",3.5);
49  desc.add<double>("VariableMaxCone",0.17);
50  desc.add<double>("VariableMinCone",0.05);
51 }
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:89
pair< float, IsolatedTauTagInfo > ConeIsolationAlgorithm::tag ( const reco::JetTracksAssociationRef &  jetTracks,
const reco::Vertex &  pv 
)

Definition at line 53 of file ConeIsolationAlgorithm.cc.

References edm::RefVector< C, T, F >::begin(), reco::IsolatedTauTagInfo::discriminator(), edm::RefVector< C, T, F >::end(), reco::IsolatedTauTagInfo::leadingSignalTrack(), bookConverter::max, min(), reco::Vertex::position(), edm::RefVector< C, T, F >::push_back(), reco::IsolatedTauTagInfo::setDiscriminator(), reco::IsolatedTauTagInfo::setLeadingTrack(), testEve_cfg::tracks, and reco::Vertex::z().

Referenced by Inspector.Inspector::SetTag().

54 {
55  const edm::RefVector<reco::TrackCollection> & tracks = jetTracks->second;
56  edm::RefVector<reco::TrackCollection> myTracks;
57 
58  // Selection of the Tracks
59  float z_pv = pv.z();
60  for(edm::RefVector<reco::TrackCollection>::const_iterator it = tracks.begin(); it!= tracks.end(); ++it)
61  {
62  if ( (*it)->pt() > m_cutMinPt &&
63  (*it)->normalizedChi2() < m_cutMaxChiSquared &&
64  fabs((*it)->dxy(pv.position())) < m_cutMaxTIP &&
65  (*it)->recHitsSize() >= (unsigned int) m_cutTotalHits &&
66  (*it)->hitPattern().numberOfValidPixelHits() >= m_cutPixelHits )
67  {
68  if (useVertexConstrain_ && z_pv > -500.) {
69  if (fabs((*it)->dz(pv.position())) < dZ_vertex)
70  myTracks.push_back(*it);
71  } else
72  myTracks.push_back(*it);
73  }
74  }
75  IsolatedTauTagInfo resultExtended(myTracks,jetTracks);
76 
77  double r_sigCone = signal_cone;
78  double energyJet = jetTracks->first->energy();
79  if (not useFixedSizeCone) {
80  r_sigCone = std::min(variableMaxCone, variableConeParameter / energyJet);
81  r_sigCone = std::max((double)r_sigCone, variableMinCone);
82  }
83 
84  // now I can use it for the discriminator;
85  math::XYZVector jetDir(jetTracks->first->px(), jetTracks->first->py(), jetTracks->first->pz());
86  float discriminator = 0.;
87  if (useVertexConstrain_) {
88  // In this case all the selected tracks comes from the same vertex, so no need to pass the dZ_vertex requirement to the discriminator
89  const TrackRef myLeadTk = resultExtended.leadingSignalTrack( matching_cone, pt_min_leadTrack);
90  resultExtended.setLeadingTrack(myLeadTk);
91  discriminator = resultExtended.discriminator(jetDir, matching_cone, r_sigCone, isolation_cone, pt_min_leadTrack, pt_min_isolation, n_tracks_isolation_ring);
92  resultExtended.setDiscriminator(discriminator);
93  } else {
94  // In this case the dZ_vertex is used to associate the tracks to the Z_imp parameter of the Leading Track
95  const TrackRef myLeadTk = resultExtended.leadingSignalTrack( matching_cone, pt_min_leadTrack);
96  resultExtended.setLeadingTrack(myLeadTk);
97  discriminator = resultExtended.discriminator(jetDir, matching_cone, r_sigCone, isolation_cone, pt_min_leadTrack, pt_min_isolation, n_tracks_isolation_ring, dZ_vertex);
98  resultExtended.setDiscriminator(discriminator);
99  }
100 
101  return std::make_pair(discriminator, resultExtended);
102 }
edm::RefVector::end
const_iterator end() const
Termination of iteration.
Definition: RefVector.h:249
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:106
edm::RefVector::begin
const_iterator begin() const
Initialize an iterator over the RefVector.
Definition: RefVector.h:244
reco::Vertex::z
double z() const
y coordinate
Definition: Vertex.h:112
min
T min(T a, T b)
Definition: MathUtil.h:58
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
edm::RefVector::push_back
void push_back(value_type const &ref)
Add a Ref&lt;C, T&gt; to the RefVector.
Definition: RefVector.h:64

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.

Generated for CMSSW Reference Manual by doxygen 1.8.5