top::Chi2LJets Class Reference

Not complete. More...

#include <Chi2LJets.h>

Inheritance diagram for top::Chi2LJets:

Collaboration diagram for top::Chi2LJets:

Public Member Functions
	Chi2LJets ()
virtual	~Chi2LJets ()
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...

Static Public Attributes
static const double	input_topmass = 172500
	The top mass that we're using as the baseline - used in chi2 function. More...
static const double	input_wmass = 80300
	The W mass that we're using as the baseline - used in chi2 function. More...
static const double	input_topsigma = 13000
	The top resolution that we're using as the baseline - used in chi2 function. More...
static const double	input_wsigma = 7000
	The W resolution that we're using as the baseline - used in chi2 function. More...

Private Attributes
TMinuit *	min
	So we can do a minimisation. More...
FitInfo	params

Detailed Description

Not complete.

Work in progress

Definition at line 44 of file Chi2LJets.h.

Constructor & Destructor Documentation

Chi2LJets()

top::Chi2LJets::Chi2LJets

(

)

Definition at line 16 of file Chi2LJets.cxx.

                       {
    min = new TMinuit(2);
    min->SetPrintLevel(-1);
    min->SetObjectFit(&params);
    min->SetFCN(fcn);
  }

~Chi2LJets()

top::Chi2LJets::~Chi2LJets ( )

virtual

Definition at line 23 of file Chi2LJets.cxx.

                        {
    delete min;
  }

Member Function Documentation

apply()

bool top::Chi2LJets::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.

The other jet in the W

The hadronic b quark

Implements top::EventSelectorBase.

Definition at line 27 of file Chi2LJets.cxx.

                                                   {
    //set lepton
    if (event.m_electrons.size() == 1) params.lepton = event.m_electrons.at(0)->p4();
 
    if (event.m_muons.size() == 1) params.lepton = event.m_muons.at(0)->p4();
 
    //and met
    params.neutrino_px = event.m_met->mpx();
    params.neutrino_py = event.m_met->mpy();
 
    //Output
    double chi2min = 1e200;
    double topmass = 0.;
    //double wmass = 0.;
 
    //First loop is for the leptonic b quark
    for (xAOD::JetContainer::const_iterator j1 = event.m_jets.begin(); j1 != event.m_jets.end(); ++j1) {
      params.leptonic_b = (*j1)->p4();
 
      //Second loop is for one of the jets that makes a W
      for (xAOD::JetContainer::const_iterator j2 = event.m_jets.begin(); j2 != event.m_jets.end(); ++j2) {
        if (j1 == j2) continue;
 
        for (xAOD::JetContainer::const_iterator j3 = event.m_jets.begin(); j3 != event.m_jets.end(); ++j3) {
          if (j1 == j3 || j2 == j3) continue;
 
          params.hadronic_w_mass = ((*j2)->p4() + (*j3)->p4()).M();
 
          for (xAOD::JetContainer::const_iterator j4 = event.m_jets.begin(); j4 != event.m_jets.end(); ++j4) {
            if (j1 == j4 || j2 == j4 || j3 == j4) continue;
 
            params.hadronic_t_mass = ((*j1)->p4() + (*j2)->p4() + (*j3)->p4()).M();
 
            //Set input parameters to something reasonable (a guess)
            min->DefineParameter(0, "top_mass", 172500., 10000., 0, 500000);         //starting guess, step, min, max
            min->DefineParameter(1, "neutrino_pz", 10000, 10000, -1000000, 1000000);
 
            //Perform the minimisation
            int status = min->Migrad();
 
            double chi2 = 0;
            double edm = 0;
            double errdef = 0;
            int nvpar = 0;
            int nparx = 0;
            int icstat = 0;
 
            min->mnstat(chi2, edm, errdef, nvpar, nparx, icstat);
 
            if (status == 0 && chi2 < chi2min) {
              double error = 0;
              chi2min = chi2;
              min->GetParameter(0, topmass, error);
              //wmass = params.hadronic_w_mass;
            } //min chi2
          } //j4 = hadronic b
        } //j3 = j3 + j2 = hadronic W
      } //j2
    } //j1 = leptonic b
 
 
 
 
    //
    //  577     MiniElectronContainer electrons = event.electrons();
    //  578     if (electrons.size() == 1)
    //  579         params.lepton.SetPtEtaPhiE(electrons[0]->pt(), electrons[0]->eta(), electrons[0]->phi(),
    // electrons[0]->E());
    //  580
    //  581     MiniMuonContainer muons = event.muons();
    //  582     if (muons.size() == 1)
    //  583         params.lepton.SetPtEtaPhiE(muons[0]->pt(), muons[0]->eta(), muons[0]->phi(), muons[0]->E());
    //  584
    //  585     //set met
    //  586     const MiniMET* met = event.MET();
    //  587     params.neutrino_px = met->et()*cos(met->phi());
    //  588     params.neutrino_py = met->et()*sin(met->phi());
    //  589
    //  590     //Out
    //  591     double chi2, edm, errdef;
    //  592     int nvpar, nparx, icstat;
    //  593
    //  594     double chi2min = 1e200;
    //  595     double topmass = 0.;
    //  596     double wmass = 0.;
    //  597
    //  598     //TObjects have a significant creation overhead so move them out of the loop
    //  599     TLorentzVector a;
    //  600     TLorentzVector b;
    //  601     TLorentzVector c;
    //  602
    //  603     const MiniJetContainer jets = event.jets();
    //  604
    //  605     //checks
    //  606     assert((electrons.size() == 1 && muons.empty()) || (muons.size() == 1 && electrons.empty()));
    //  607     assert(jets.size() >= 4);
    //  608
    //  609     //First loop is for the leptonic b quark
    //  610     for (MiniJetContainer::const_iterator j1 = jets.begin(); j1 !=jets.end(); ++j1) {
    //  611
    //  612         params.leptonic_b.SetPtEtaPhiE((*j1)->pt(), (*j1)->eta(), (*j1)->phi(), (*j1)->E());
    //  613
    //  614         //Second loop is for one of the jets that makes a W
    //  615         for (MiniJetContainer::const_iterator j2 = jets.begin(); j2 != jets.end(); ++j2) {
    //  616             if (j1 == j2)
    //  617                 continue;
    //  618
    //  619             a.SetPtEtaPhiE((*j2)->pt(), (*j2)->eta(), (*j2)->phi(), (*j2)->E());
    //  620
    //  621             ///The other jet in the W
    //  622             for (MiniJetContainer::const_iterator j3 = jets.begin(); j3 != jets.end(); ++j3) {
    //  623                 if (j1 == j3 || j2 == j3)
    //  624                     continue;
    //  625
    //  626                 b.SetPtEtaPhiE((*j3)->pt(), (*j3)->eta(), (*j3)->phi(), (*j3)->E());
    //  627                 params.hadronic_w_mass = (a + b).M();
    //  628
    //  629                 ///The hadronic b quark
    //  630                 for (MiniJetContainer::const_iterator j4 = jets.begin(); j4 != jets.end(); ++j4) {
    //  631                     if (j1 == j4 || j2 == j4 || j3 == j4)
    //  632                         continue;
    //  633
    //  634                     c.SetPtEtaPhiE((*j4)->pt(), (*j4)->eta(), (*j4)->phi(), (*j4)->E());
    //  635                     params.hadronic_t_mass = (a + b + c).M();
    //  636
    //  637                     //Set input parameters to something reasonable (a guess)
    //  638                     min->DefineParameter(0, "top_mass", 172500., 10000., 0, 500000);         //starting
    // guess, step, min, max
    //  639                     min->DefineParameter(1, "neutrino_pz", 10000, 10000, -1000000, 1000000);
    //  640
    //  641                     //Perform the minimisation
    //  642                     int status = min->Migrad();
    //  643                     min->mnstat(chi2, edm, errdef, nvpar, nparx, icstat);
    //  644
    //  645                     if (status == 0 && chi2 < chi2min) {
    //  646                         double error;
    //  647                         chi2min = chi2;
    //  648                         min->GetParameter(0, topmass, error);
    //  649                         wmass = params.hadronic_w_mass;
    //  650                     } //min chi2
    //  651
    //  652                 } //j4 = hadronic b
    //  653             } //j3 = j3 + j2 = hadronic W
    //  654         } //j2
    //  655     } //j1 = leptonic b
 
    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::Chi2LJets::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 51 of file Chi2LJets.h.

{return "RECO:CHI2LJETS";}

Member Data Documentation

input_topmass

const double top::Chi2LJets::input_topmass = 172500

static

The top mass that we're using as the baseline - used in chi2 function.

Definition at line 61 of file Chi2LJets.h.

input_topsigma

const double top::Chi2LJets::input_topsigma = 13000

static

The top resolution that we're using as the baseline - used in chi2 function.

Definition at line 67 of file Chi2LJets.h.

input_wmass

const double top::Chi2LJets::input_wmass = 80300

static

The W mass that we're using as the baseline - used in chi2 function.

Definition at line 64 of file Chi2LJets.h.

input_wsigma

const double top::Chi2LJets::input_wsigma = 7000

static

The W resolution that we're using as the baseline - used in chi2 function.

Definition at line 70 of file Chi2LJets.h.

min

TMinuit* top::Chi2LJets::min

private

So we can do a minimisation.

Definition at line 54 of file Chi2LJets.h.

params

FitInfo top::Chi2LJets::params

mutableprivate

Definition at line 56 of file Chi2LJets.h.

---

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

• Chi2LJets.h
• Chi2LJets.cxx