HEPTopTaggerV2 Class Reference

#include <HEPTopTaggerWrapperV2.h>

Inheritance diagram for HEPTopTaggerV2:

Public Types
typedef HEPTopTaggerV2Structure	StructureType

Public Member Functions
std::string	description () const override
	returns a textual description of the tagger More...
	HEPTopTaggerV2 (bool DoOptimalR, bool DoQjets, double minSubjetPt, double minCandPt, double subjetMass, double muCut, double filtR, int filtN, int mode, double minCandMass, double maxCandMass, double massRatioWidth, double minM23Cut, double minM13Cut, double maxM13Cut, bool optRrejectMin)
PseudoJet	result (const PseudoJet &jet) const override
void	set_rng (CLHEP::HepRandomEngine *engine)

Private Attributes
bool	DoOptimalR_
bool	DoQjets_
CLHEP::HepRandomEngine *	engine_
int	filtN_
double	filtR_
double	massRatioWidth_
double	maxCandMass_
double	maxM13Cut_
double	minCandMass_
double	minCandPt_
double	minM13Cut_
double	minM23Cut_
double	minSubjetPt_
int	mode_
double	muCut_
bool	optRrejectMin_
double	subjetMass_

Detailed Description

Definition at line 49 of file HEPTopTaggerWrapperV2.h.

Member Typedef Documentation

typedef HEPTopTaggerV2Structure HEPTopTaggerV2::StructureType

Definition at line 97 of file HEPTopTaggerWrapperV2.h.

Constructor & Destructor Documentation

HEPTopTaggerV2::HEPTopTaggerV2 ( bool  DoOptimalR, bool  DoQjets, double  minSubjetPt, double  minCandPt, double  subjetMass, double  muCut, double  filtR, int  filtN, int  mode, double  minCandMass, double  maxCandMass, double  massRatioWidth, double  minM23Cut, double  minM13Cut, double  maxM13Cut, bool  optRrejectMin  )

inline

Definition at line 51 of file HEPTopTaggerWrapperV2.h.

References description(), metsig::jet, result(), and AlCaHLTBitMon_QueryRunRegistry::string.

                            : DoOptimalR_(DoOptimalR),
    DoQjets_(DoQjets),
    minSubjetPt_(minSubjetPt),
    minCandPt_(minCandPt),
    subjetMass_(subjetMass),
    muCut_(muCut),
    filtR_(filtR),
    filtN_(filtN),
    mode_(mode),
    minCandMass_(minCandMass),
    maxCandMass_(maxCandMass),
    massRatioWidth_(massRatioWidth),
    minM23Cut_(minM23Cut),
    minM13Cut_(minM13Cut),
    maxM13Cut_(maxM13Cut),
    optRrejectMin_(optRrejectMin),
    engine_(nullptr)
  {}

Member Function Documentation

std::string HEPTopTaggerV2::description ( ) const

inlineoverride

returns a textual description of the tagger

Definition at line 277 of file HEPTopTaggerWrapperV2.h.

References filtN_, filtR_, massRatioWidth_, maxCandMass_, maxM13Cut_, minCandMass_, minCandPt_, minM13Cut_, minM23Cut_, minSubjetPt_, mode_, muCut_, and subjetMass_.

Referenced by HEPTopTaggerV2().

                                                  { 

  std::ostringstream oss;
  oss << "HEPTopTaggerV2 with: "
      << "minSubjetPt = " << minSubjetPt_ 
      << "minCandPt = " << minCandPt_ 
      << "subjetMass = " << subjetMass_ 
      << "muCut = " << muCut_ 
      << "filtR = " << filtR_ 
      << "filtN = " << filtN_     
      << "mode = " << mode_ 
      << "minCandMass = " << minCandMass_ 
      << "maxCandMass = " << maxCandMass_ 
      << "massRatioWidth = " << massRatioWidth_ 
      << "minM23Cut = " << minM23Cut_ 
      << "minM13Cut = " << minM13Cut_ 
      << "maxM13Cut = " << maxM13Cut_ << std::endl;
  return oss.str();
}

PseudoJet HEPTopTaggerV2::result ( const PseudoJet &  jet ) const

override

runs the tagger on the given jet and returns the tagged PseudoJet if successful, or a PseudoJet==0 otherwise (standard access is through operator()).

Parameters

jet	the PseudoJet to tag

Definition at line 65 of file HEPTopTaggerWrapperV2.cc.

References external::HEPTopTaggerV2::do_optimalR(), external::HEPTopTaggerV2::do_qjets(), external::HEPTopTaggerV2::f_rec(), objects.autophobj::float, external::HEPTopTaggerV2::is_masscut_passed(), external::HEPTopTaggerV2::is_tagged(), join(), external::HEPTopTaggerV2::nsub_filtered(), external::HEPTopTaggerV2::nsub_unfiltered(), external::HEPTopTaggerV2::pruned_mass(), external::HEPTopTaggerV2::pt_for_Ropt_calc(), external::HEPTopTaggerV2::q_weight(), R_min_expected_function(), mps_fire::result, external::HEPTopTaggerV2::Ropt(), external::HEPTopTaggerV2::Ropt_calc(), external::HEPTopTaggerV2::run(), alignCSCRings::s, external::HEPTopTaggerV2::set_filtering_minpt_subjet(), external::HEPTopTaggerV2::set_filtering_n(), external::HEPTopTaggerV2::set_filtering_R(), external::HEPTopTaggerV2::set_mass_drop_threshold(), external::HEPTopTaggerV2::set_mass_ratio_cut(), external::HEPTopTaggerV2::set_mass_ratio_range(), external::HEPTopTaggerV2::set_max_subjet_mass(), external::HEPTopTaggerV2::set_mode(), external::HEPTopTaggerV2::set_optimalR_calc_fun(), external::HEPTopTaggerV2::set_optimalR_reject_minimum(), external::HEPTopTaggerV2::set_qjets(), external::HEPTopTaggerV2::set_qjets_rng(), external::HEPTopTaggerV2::set_top_mass_range(), external::HEPTopTaggerV2::set_top_minpt(), mathSSE::sqrt(), external::HEPTopTaggerV2::t(), external::HEPTopTaggerV2::top_subjets(), and external::HEPTopTaggerV2::unfiltered_mass().

Referenced by HEPTopTaggerV2().

                                                           {

  // make sure that there is a "regular" cluster sequence associated
  // with the jet. Note that we also check it is valid (to avoid a
  // more criptic error later on)
  if (!jet.has_valid_cluster_sequence()){
    throw Error("HEPTopTagger can only be applied on jets having an associated (and valid) ClusterSequence");
  }

  external::HEPTopTaggerV2 tagger(jet);
  
  external::HEPTopTaggerV2 best_tagger;

  // translate the massRatioWidth (which should be the half-width given in %) 
  // to values useful for the A-shape cuts
  double mw_over_mt = 80.4/172.3;
  double ratio_min = mw_over_mt * (100.-massRatioWidth_)/100.;
  double ratio_max = mw_over_mt * (100.+massRatioWidth_)/100.;
  
  // Unclustering, Filtering & Subjet Settings
  tagger.set_max_subjet_mass(subjetMass_);
  tagger.set_mass_drop_threshold(muCut_);
  tagger.set_filtering_R(filtR_);
  tagger.set_filtering_n(filtN_);
  tagger.set_filtering_minpt_subjet(minSubjetPt_); 

  // Optimal R
  tagger.do_optimalR(DoOptimalR_);
  tagger.set_optimalR_reject_minimum(optRrejectMin_);

  // How to select among candidates
  tagger.set_mode((external::Mode)mode_);
  
  // Requirements to accept a candidate
  tagger.set_top_minpt(minCandPt_); 
  tagger.set_top_mass_range(minCandMass_, maxCandMass_); 
  tagger.set_mass_ratio_cut(minM23Cut_, minM13Cut_, maxM13Cut_);
  tagger.set_mass_ratio_range(ratio_min, ratio_max);

  // Set function to calculate R_min_expected
  tagger.set_optimalR_calc_fun(R_min_expected_function);

  
  double qweight  = -1;
  double qepsilon = -1;
  double qsigmaM  = -1;

  if (DoQjets_){
    
    int niter(100);
    double q_zcut(0.1);
    double q_dcut_fctr(0.5);
    double q_exp_min(0.);
    double q_exp_max(0.);
    double q_rigidity(0.1);
    double q_truncation_fctr(0.0);
           
    double weight_q1 = -1.;
    double m_sum = 0.;
    double m2_sum = 0.;
    int qtags = 0;

    tagger.set_qjets(q_zcut,
             q_dcut_fctr,
             q_exp_min,
             q_exp_max,
             q_rigidity,
             q_truncation_fctr);
    tagger.set_qjets_rng(engine_);    
    tagger.do_qjets(true);
    tagger.run();

    for (int iq = 0; iq < niter; iq++) {
      tagger.run();
      if (tagger.is_tagged()) {
    qtags++;
    m_sum += tagger.t().m();
    m2_sum += tagger.t().m() * tagger.t().m();
    if (tagger.q_weight() > weight_q1) {
      best_tagger = tagger;
      weight_q1=tagger.q_weight();     
    }        
      }
    }
    
    tagger = best_tagger;
    qweight = weight_q1;
    qepsilon = float(qtags)/float(niter);
    
    // calculate width of tagged mass distribution if we have at least one candidate
    if (qtags > 0){
      double mean_m = m_sum / qtags;
      double mean_m2 = m2_sum / qtags;
      qsigmaM = sqrt(mean_m2 - mean_m*mean_m);
    }
  }
  else{
    tagger.run();
  }

  // Requires:
  //   - top mass window
  //   - mass ratio cuts
  //   - minimal candidate pT
  // If this is not intended: use loose top mass and ratio windows
  if (!tagger.is_tagged())
    return PseudoJet();
  
  // create the result and its structure
  const JetDefinition::Recombiner *rec
    = jet.associated_cluster_sequence()->jet_def().recombiner();

  const vector<PseudoJet>& subjets = tagger.top_subjets();
  assert(subjets.size() == 3);

  PseudoJet non_W = subjets[0];
  PseudoJet W1 = subjets[1];
  PseudoJet W2 = subjets[2];
  PseudoJet W = join(subjets[1], subjets[2], *rec);

  PseudoJet result = join<HEPTopTaggerV2Structure>( W1, W2, non_W, *rec);
  HEPTopTaggerV2Structure *s = (HEPTopTaggerV2Structure*) result.structure_non_const_ptr();

  s->_fj_mass  = jet.m();
  s->_fj_pt    = jet.perp();
  s->_fj_eta   = jet.eta();
  s->_fj_phi   = jet.phi();

  s->_top_mass = tagger.t().m();
  s->_pruned_mass = tagger.pruned_mass();
  s->_unfiltered_mass = tagger.unfiltered_mass();
  s->_fRec = tagger.f_rec();
  s->_mass_ratio_passed = tagger.is_masscut_passed();

  if (DoOptimalR_){
    s->_tau1Unfiltered = tagger.nsub_unfiltered(1);
    s->_tau2Unfiltered = tagger.nsub_unfiltered(2);
    s->_tau3Unfiltered = tagger.nsub_unfiltered(3);
    s->_tau1Filtered = tagger.nsub_filtered(1);
    s->_tau2Filtered = tagger.nsub_filtered(2);
    s->_tau3Filtered = tagger.nsub_filtered(3);
  }

  s->_qweight = qweight;
  s->_qepsilon = qepsilon;
  s->_qsigmaM = qsigmaM;


  if (DoOptimalR_){
    s->_ropt = tagger.Ropt();
    s->_roptCalc = tagger.Ropt_calc();
    s->_ptForRoptCalc = tagger.pt_for_Ropt_calc();
  }
  else {
    s->_ropt = -1;
    s->_roptCalc = -1;
    s->_ptForRoptCalc = -1;
  }

  // Removed selectors as all cuts are applied in HTT

  return result;
}

void HEPTopTaggerV2::set_rng ( CLHEP::HepRandomEngine *  engine )

inline

Definition at line 94 of file HEPTopTaggerWrapperV2.h.

References engine_.

{ engine_ = engine;}

Member Data Documentation

bool HEPTopTaggerV2::DoOptimalR_

private

Definition at line 100 of file HEPTopTaggerWrapperV2.h.

bool HEPTopTaggerV2::DoQjets_

private

Definition at line 101 of file HEPTopTaggerWrapperV2.h.

CLHEP::HepRandomEngine* HEPTopTaggerV2::engine_

private

Definition at line 134 of file HEPTopTaggerWrapperV2.h.

Referenced by set_rng().

int HEPTopTaggerV2::filtN_

private

Definition at line 110 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::filtR_

private

Definition at line 109 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::massRatioWidth_

private

Definition at line 125 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::maxCandMass_

private

Definition at line 123 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::maxM13Cut_

private

Definition at line 128 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::minCandMass_

private

Definition at line 122 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::minCandPt_

private

Definition at line 104 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::minM13Cut_

private

Definition at line 127 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::minM23Cut_

private

Definition at line 126 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::minSubjetPt_

private

Definition at line 103 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

int HEPTopTaggerV2::mode_

private

Definition at line 119 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

double HEPTopTaggerV2::muCut_

private

Definition at line 107 of file HEPTopTaggerWrapperV2.h.

Referenced by description().

bool HEPTopTaggerV2::optRrejectMin_

private

Definition at line 130 of file HEPTopTaggerWrapperV2.h.

double HEPTopTaggerV2::subjetMass_

private

Definition at line 106 of file HEPTopTaggerWrapperV2.h.

Referenced by description().