top::Sonnenschein Class Reference

An AnalysisTop friendly wrapper for the Sonnenschein reconstruction that is implemented in SonnenscheinEngine. More...

#include <Sonnenschein.h>

Inheritance diagram for top::Sonnenschein:

Collaboration diagram for top::Sonnenschein:

Public Member Functions
	Sonnenschein ()
virtual	~Sonnenschein ()
virtual 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
bool	tWb (const top::Event &event, int start, TLorentzVector &t, TLorentzVector &W, TLorentzVector &b, TLorentzVector &l, TLorentzVector &nu) const

Private Attributes
top::SonnenscheinEngine	m_sonn

Detailed Description

An AnalysisTop friendly wrapper for the Sonnenschein reconstruction that is implemented in SonnenscheinEngine.

Definition at line 20 of file Sonnenschein.h.

Constructor & Destructor Documentation

Sonnenschein()

top::Sonnenschein::Sonnenschein

(

)

Definition at line 14 of file Sonnenschein.cxx.

                             {
  }

~Sonnenschein()

top::Sonnenschein::~Sonnenschein ( )

virtual

Definition at line 17 of file Sonnenschein.cxx.

                              {
  }

Member Function Documentation

apply()

bool top::Sonnenschein::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 71 of file Sonnenschein.cxx.

                                                      {
    TLorentzVector t;
    TLorentzVector Wp;
    TLorentzVector b;
    TLorentzVector lp;
    TLorentzVector nu;
    bool cando1 = tWb(event, 6, t, Wp, b, lp, nu);
 
    TLorentzVector tbar;
    TLorentzVector Wm;
    TLorentzVector bbar;
    TLorentzVector lm;
    TLorentzVector nubar;
    bool cando2 = tWb(event, -6, tbar, Wm, bbar, lm, nubar);
 
    //isn't dilepton at truth level, so return
    if (!cando1 || !cando2) return true;
 
    const double met_ex = nu.Px() + nubar.Px();
    const double met_ey = nu.Py() + nubar.Py();
 
    //solve
    std::vector<std::pair<TLorentzVector, TLorentzVector> > allSolutions = m_sonn.solve(lp, b, t.M(),
                                                                                        Wp.M(), lm, bbar,
                                                                                        tbar.M(), Wm.M(), met_ex,
                                                                                        met_ey);
 
    std::cout << "EVENT" << std::endl;
    unsigned int counter = 0;
    for (const auto& possible : allSolutions) {
      std::cout << "Possible solution number " << counter << "\n";
      std::cout << "    first neutrino : " << possible.first << "\n";
      std::cout << "    second neutrino: " << possible.second << "\n";
      ++counter;
    }
 
    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::Sonnenschein::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 27 of file Sonnenschein.h.

{return "RECO:SONNENSCHEIN";}

tWb()

bool top::Sonnenschein::tWb ( const top::Event &  event, int  start, TLorentzVector &  t, TLorentzVector &  W, TLorentzVector &  b, TLorentzVector &  l, TLorentzVector &  nu  ) const

private

Definition at line 20 of file Sonnenschein.cxx.

                                                                      {
    bool hasT = false;
    bool hasB = false;
    bool hasL = false;
    bool hasNu = false;
 
    for (const xAOD::TruthParticle* particle : *event.m_truth) {
      if (particle->pdgId() == start) {
        //std::cout << *particle << std::endl;
        t = particle->p4();
        hasT = true;
 
        if (particle->hasDecayVtx()) {
          const xAOD::TruthVertex* vtx = particle->decayVtx();
          for (size_t i = 0; i < vtx->nOutgoingParticles(); ++i) {
            const xAOD::TruthParticle* tDecayProduct = vtx->outgoingParticle(i);
            //std::cout << "    " << *tDecayProduct << std::endl;
 
            if (abs(tDecayProduct->pdgId()) == 24 && tDecayProduct->hasDecayVtx()) {
              W = tDecayProduct->p4();
 
              const xAOD::TruthVertex* wDecayVtx = tDecayProduct->decayVtx();
              for (size_t j = 0; j < wDecayVtx->nOutgoingParticles(); ++j) {
                const xAOD::TruthParticle* wDecayProduct = wDecayVtx->outgoingParticle(j);
                //std::cout << "       " << *wDecayProduct << std::endl;
 
                if (abs(wDecayProduct->pdgId()) == 11 || abs(wDecayProduct->pdgId()) == 13 ||
                    abs(wDecayProduct->pdgId()) == 15) {
                  l = wDecayProduct->p4();
                  hasL = true;
                }
 
                if (abs(wDecayProduct->pdgId()) == 12 || abs(wDecayProduct->pdgId()) == 14 ||
                    abs(wDecayProduct->pdgId()) == 16) {
                  nu = wDecayProduct->p4();
                  hasNu = true;
                }
              }
            } else if (abs(tDecayProduct->pdgId()) == 5) {
              b = tDecayProduct->p4();
              hasB = true;
            }
          }
        }
      }
    }
 
    return hasT && hasB && hasL && hasNu;
  }

Member Data Documentation

m_sonn

top::SonnenscheinEngine top::Sonnenschein::m_sonn

private

Definition at line 32 of file Sonnenschein.h.

---

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

• Sonnenschein.h
• Sonnenschein.cxx