Skip to content
Snippets Groups Projects

[Run2_2017] Additional information needed for the HH+MET Synchronization Exercise

Merged Alexx Perloff requested to merge github/fork/aperloff/Run2_2017_Synchronization into Run2_2017

Created by: aperloff

This PR adds some additional information needed for the HH+MET synchronization exercise. In particular it adds the following variables to the PFCandidates:

  1. The PF Relative Isolation (DR=0.3) for only the charged component
  2. The PF Relative Isolation (DR=0.3) for all components
  3. The dxy impact parameter These branches are only added if the 'debugtracks' option is set to True.

Adding these branches necessitated adding a new isolation calculator to the SUSYIsolation class. There is one change (L89 in Utils/interface/get_isolation_activity.h) that needs some special review. This as added based upon https://github.com/cms-sw/cmssw/blob/56d3b7f00bac17989db6efb3efaf8e1ba31f2dfe/PhysicsTools/PatAlgos/plugins/PATIsolatedTrackProducer.cc#L554, which is the producer that makes the 'isolatedTrack' collection for MiniAOD. I think the change is well motivated.

This PR has been unit tested and compared to 37 events (50 PF candidates) from NanoAOD. The comparisons show that the new isolation values match up to 4 decimal places.

Merge request reports

Loading
Loading

Activity

Filter activity
  • Approvals
  • Assignees & reviewers
  • Comments (from bots)
  • Comments (from users)
  • Commits & branches
  • Edits
  • Labels
  • Lock status
  • Mentions
  • Merge request status
  • Tracking
127 127 return iso;
128 128 }
129
129
130 void GetPFIsolation(edm::Handle<edm::View<pat::PackedCandidate> > pfcands, const unsigned pfcand_idx,
131 double& pf_iso_all, double& pf_iso_chg, //return values
132 double r_pfiso=0.3, double ptThresh=0.0) const
133 {
134 pf_iso_all = 0.0;
135 pf_iso_chg = 0.0;
136 const auto& ptcl = pfcands->at(pfcand_idx);
137 if (ptcl.p4().Pt()<5.) {
138 pf_iso_all = 99999.;
139 pf_iso_chg = 99999.;
140 return;
141 }
  • Kevin Pedro
  • Alexx Perloff @aperloff started a thread on commit 438a17f1
  • 127 127 return iso;
    128 128 }
    129
    129
    130 void GetPFIsolation(edm::Handle<edm::View<pat::PackedCandidate> > pfcands, const unsigned pfcand_idx,
    131 double& pf_iso_all, double& pf_iso_chg, //return values
    132 double r_pfiso=0.3, double ptThresh=0.0) const
    133 {
    134 pf_iso_all = 0.0;
    135 pf_iso_chg = 0.0;
    136 const auto& ptcl = pfcands->at(pfcand_idx);
    137 if (ptcl.p4().Pt()<5.) {
    138 pf_iso_all = 99999.;
    139 pf_iso_chg = 99999.;
    140 return;
    141 }
    Please register or sign in to reply
    Loading