CATopJetHelperUser Class Reference

#include <HLTCATopTagFilter.h>

Public Member Functions
	CATopJetHelperUser (double TopMass)
reco::CATopJetProperties	operator() (reco::Jet const &ihardJet) const

Protected Attributes
double	TopMass_

Detailed Description

Definition at line 27 of file HLTCATopTagFilter.h.

Constructor & Destructor Documentation

CATopJetHelperUser()

CATopJetHelperUser::CATopJetHelperUser ( double  TopMass )

inline

Definition at line 29 of file HLTCATopTagFilter.h.

: TopMass_(TopMass) {}

Member Function Documentation

operator()()

reco::CATopJetProperties CATopJetHelperUser::operator() ( reco::Jet const &  ihardJet ) const

inline

Definition at line 69 of file HLTCATopTagFilter.h.

References reco::Jet::getJetConstituents(), reco::LeafCandidate::mass(), reco::CATopJetProperties::minMass, reco::CATopJetProperties::nSubJets, jetsAK4_CHS_cff::sort, and reco::CATopJetProperties::topMass.

                                                                                          {
  reco::CATopJetProperties properties;
  // Get subjets
  reco::Jet::Constituents subjets = ihardJet.getJetConstituents();
  properties.nSubJets = subjets.size();  // number of subjets
  properties.topMass = ihardJet.mass();  // jet mass
  properties.minMass = 999999.;          // minimum mass pairing

  // Require at least three subjets in all cases, if not, untagged
  if (properties.nSubJets >= 3) {
    // Take the highest 3 pt subjets for cuts
    sort(subjets.begin(), subjets.end(), GreaterByPtCandPtrUser());

    // Now look at the subjets that were formed
    for (int isub = 0; isub < 2; ++isub) {
      // Get this subjet
      reco::Jet::Constituent icandJet = subjets[isub];

      // Now look at the "other" subjets than this one, form the minimum invariant mass
      // pairing, as well as the "closest" combination to the W mass
      for (int jsub = isub + 1; jsub < 3; ++jsub) {
        // Get the second subjet
        reco::Jet::Constituent jcandJet = subjets[jsub];

        reco::Candidate::LorentzVector wCand = icandJet->p4() + jcandJet->p4();

        // Get the candidate mass
        double imw = wCand.mass();

        // Find the minimum mass pairing.
        if (fabs(imw) < properties.minMass) {
          properties.minMass = imw;
        }
      }  // end second loop over subjets
    }    // end first loop over subjets
  }      // endif 3 subjets

  if (properties.minMass == 999999)
    properties.minMass = -1;

  return properties;
}

Member Data Documentation

TopMass_

double CATopJetHelperUser::TopMass_

protected

Definition at line 34 of file HLTCATopTagFilter.h.