top::PTMaxReco Class Reference

Not complete. More...

#include <PTMaxReco.h>

Inheritance diagram for top::PTMaxReco:

Collaboration diagram for top::PTMaxReco:

Public Member Functions
	PTMaxReco ()
virtual	~PTMaxReco ()
bool	apply (const top::Event &) const override
	This does stuff based on the information in an event. More...
std::string	name () const override
	A human readable name. More...
virtual bool	applyParticleLevel (const top::ParticleLevelEvent &) const
	This does stuff based on the information in a particle level event. More...

Private Member Functions
TLorentzVector	neutrinoCandidate (const TLorentzVector &lep, const xAOD::MissingET &met, bool dealWithNegative_nu) const

Private Attributes
double	m_wmass
double	m_topmass

Detailed Description

Not complete.

Work in progress

Definition at line 18 of file PTMaxReco.h.

Constructor & Destructor Documentation

PTMaxReco()

top::PTMaxReco::PTMaxReco

(

)

Definition at line 10 of file PTMaxReco.cxx.

                       : m_wmass(80400.), m_topmass(172500.) {
  }

~PTMaxReco()

top::PTMaxReco::~PTMaxReco ( )

virtual

Definition at line 13 of file PTMaxReco.cxx.

                        {
  }

Member Function Documentation

apply()

bool top::PTMaxReco::apply ( const top::Event &  ) const

overridevirtual

This does stuff based on the information in an event.

The idea is that you implement this to return either true or false based on the information held within top::Event. If this returns true then the event is kept. If it returns false then the event is removed.

Parameters

top::Event

The current event.

Returns

true if the event should be kept, false otherwise.

Implements top::EventSelectorBase.

Definition at line 16 of file PTMaxReco.cxx.

                                                   {
    //get the hadronic top
 
    //jets that make-up the hadronic top quark
    const xAOD::Jet* jet1 = nullptr;
    const xAOD::Jet* jet2 = nullptr;
    const xAOD::Jet* jet3 = nullptr;
 
    double maxptsq = 0.;
 
    for (xAOD::JetContainer::const_iterator j1 = event.m_jets.begin(); j1 != event.m_jets.end(); ++j1) {
      for (xAOD::JetContainer::const_iterator j2 = j1 + 1; j2 != event.m_jets.end(); ++j2) {
        for (xAOD::JetContainer::const_iterator j3 = j2 + 1; j3 != event.m_jets.end(); ++j3) {
          const double px = (*j1)->px() + (*j2)->px() + (*j3)->px();
          const double py = (*j1)->py() + (*j2)->py() + (*j3)->py();
 
          //pt2 to avoid the slow sqrt
          const double ptsq = px * px + py * py;
 
          if (ptsq > maxptsq) {
            jet1 = (*j1);
            jet2 = (*j2);
            jet3 = (*j3);
          }
        }
      }
    }
 
    //get the leptonic top
 
    //e+jets
    TLorentzVector lepton;
    if (event.m_electrons.size() == 1) lepton = event.m_electrons.at(0)->p4();
 
    //mu+jets
    if (event.m_muons.size() == 1) lepton = event.m_muons.at(0)->p4();
 
    TLorentzVector nu = neutrinoCandidate(lepton, *event.m_met, true);
 
    const xAOD::Jet* lepb = nullptr;
    double topdiff = 1000000.;
    for (xAOD::JetContainer::const_iterator j1 = event.m_jets.begin(); j1 != event.m_jets.end(); ++j1) {
      //ignore the three jets that make-up the hadronic top quark
      if (*j1 == jet1 || *j1 == jet2 || *j1 == jet3) continue;
 
      //pick the jet which makes a leptonic top closest to the true top mass
      const double tempdiff = fabs((lepton + nu + (*j1)->p4()).M());
 
      if ((tempdiff - m_topmass) < topdiff) {
        lepb = *j1;
        topdiff = tempdiff;
      }
    }
 
    if (lepb) std::cout << lepb->pt() << std::endl;
 
    return true;
  }

applyParticleLevel()

virtual bool top::EventSelectorBase::applyParticleLevel ( const top::ParticleLevelEvent &  ) const

inlinevirtualinherited

This does stuff based on the information in a particle level event.

The idea is that you implement this to return either true or false, based on the information held within the top::ParticleLevelEvent. If this function returns true, then the event is kept, otherwise it is removed. The function has a default implementation (which returns true) because it is expected that many EventSelector objects do not operate on ParticleLevelEvent objects.

Parameters

top::ParticleLevelEvent	the current particle level event.
true	if the event should be kept (i.e. it passed the selector criteria), false otherwise.

Reimplemented in top::JetNGhostSelector, top::PrintEventSelector, top::PseudoTopRecoRun, top::NElectronNMuonTightSelector, top::NElectronNMuonSelector, top::NFwdElectronSelector, top::HTSelector, top::OSLeptonTightSelector, top::MLLSelector, top::MWTSelector, top::NElectronTightSelector, top::NFwdElectronTightSelector, top::NMuonTightSelector, top::OSLeptonSelector, top::METMWTSelector, top::METSelector, top::MLLWindow, top::NElectronSelector, top::NJetSelector, top::NMuonSelector, top::NPhotonSelector, top::NSoftMuonSelector, top::NTauSelector, top::SSLeptonTightSelector, top::SSLeptonSelector, top::ParticleLevelSelector, top::RecoLevelSelector, top::NVarRCJetSelector, top::NLargeJetSelector, and top::NRCJetSelector.

Definition at line 73 of file EventSelectorBase.h.

{return true;}

name()

std::string top::PTMaxReco::name ( ) const

inlineoverridevirtual

A human readable name.

Mostly used for printing the cut and value to the screen. Must be implemented for each tool.

Implements top::EventSelectorBase.

Definition at line 25 of file PTMaxReco.h.

{return "RECO:PTMAX";}

neutrinoCandidate()

TLorentzVector top::PTMaxReco::neutrinoCandidate ( const TLorentzVector &  lep, const xAOD::MissingET &  met, bool  dealWithNegative_nu  ) const

private

Definition at line 75 of file PTMaxReco.cxx.

                                                                              {
    const double px = met.mpx();
    const double py = met.mpy();
 
    // solve quadratic to find neutrino four vector
    double alpha = pow(m_wmass, 2) + pow((px + lep.Px()), 2) + pow((py + lep.Py()), 2) - pow(lep.E(), 2);
    double beta = 0.5 * (alpha - pow(met.met(), 2) + pow(lep.Pz(), 2));
    double gamma = -(beta * beta - (pow(lep.E(), 2) * pow(met.met(), 2))) / (pow(lep.E(), 2) - pow(lep.Pz(), 2));
    double lambda = 2. * beta * lep.Pz() / (pow(lep.E(), 2) - pow(lep.Pz(), 2));
    double delta = pow(lambda, 2) - 4 * gamma;
 
    if (delta < 0) { // ignore non real solutions
      if (dealWithNegative_nu) delta = 0;
      else return TLorentzVector();
    }
 
    delta = sqrt(delta);
 
    double pz_pos = (lambda + delta) / 2.;
    double pz_neg = (lambda - delta) / 2.;
 
    double pz = 0.0;
 
    if (fabs(pz_pos) > fabs(pz_neg)) pz = pz_neg;
    else pz = pz_pos;
 
    double e = sqrt(px * px + py * py + pz * pz);
 
    TLorentzVector neutrino;
    neutrino.SetPxPyPzE(px, py, pz, e);
    return neutrino;
  }

Member Data Documentation

m_topmass

double top::PTMaxReco::m_topmass

private

Definition at line 31 of file PTMaxReco.h.

m_wmass

double top::PTMaxReco::m_wmass

private

Definition at line 30 of file PTMaxReco.h.

---

The documentation for this class was generated from the following files:

• PTMaxReco.h
• PTMaxReco.cxx