hitfit::Lepjets_Event Class Reference

Represent a simple event consisting of lepton(s) and jet(s). An instance of this class holds a list of leptons (as represented by the Lepjets_Event_Lep class) and jets (as represented by Lepjets_Event_Jet class). Also recorded are: More...

#include <Lepjets_Event.h>

Public Member Functions
void	add_jet (const Lepjets_Event_Jet &jet)
	Add a new jet to the event. More...
void	add_lep (const Lepjets_Event_Lep &lep)
	Add a new lepton to the event. More...
int	cut_jets (double pt_cut, double eta_cut)
	Remove jets which fail transverse momentum and pseudorapidity cut. More...
int	cut_leps (double pt_cut, double eta_cut)
	Remove leptons which fail transverse momentum and pseudorapidity cut. More...
double &	dlb ()
	Return a reference to the value of low-bias (LB) discriminant (Irrelevant for non-D0 experiment). More...
double	dlb () const
	Return the value of low-bias (LB) discriminant (Irrelevant for non-D0 experiment). More...
double &	dnn ()
	Return a reference to the value of neural network (NN) discriminant (Irrelevant for non-D0 experiment). More...
double	dnn () const
	Return a the value of neural network (NN) discriminant (Irrelevant for non-D0 experiment). More...
std::ostream &	dump (std::ostream &s, bool full=false) const
	Print the content of this object. More...
const int &	evnum () const
	Return a constant reference to the event number. More...
int &	evnum ()
	Return a reference to the event number. More...
bool	isMC () const
	Return the Monte Carlo flag. More...
Lepjets_Event_Jet &	jet (std::vector< Lepjets_Event_Jet >::size_type i)
	Return a reference to jet at index position i. More...
const Lepjets_Event_Jet &	jet (std::vector< Lepjets_Event_Jet >::size_type i) const
	Return a constant reference to jet at index position i. More...
std::string	jet_permutation () const
	Return a string representing the jet permutation. The following notation is used for each type of jet: More...
std::vector< int >	jet_types () const
	Return the jet types in the event. More...
Fourvec	kt () const
	Return the sum of all objects' four-momentum and missing transverse energy. More...
Resolution &	kt_res ()
	Return a reference to the resolution. More...
const Resolution &	kt_res () const
	Return a const reference to the resolution. More...
Lepjets_Event_Lep &	lep (std::vector< Lepjets_Event_Lep >::size_type i)
	Return a reference to lepton at index position i. More...
const Lepjets_Event_Lep &	lep (std::vector< Lepjets_Event_Lep >::size_type i) const
	Return a constant reference to lepton at index position i. More...
	Lepjets_Event (int runnum, int evnum)
	Constructor. More...
Fourvec &	met ()
	Return a reference to the missing transverse energy. More...
const Fourvec &	met () const
	Return a constant reference to the missing transverse energy. More...
std::vector< Lepjets_Event_Jet >::size_type	njets () const
	Return the number of jets in the event. More...
std::vector< Lepjets_Event_Lep >::size_type	nleps () const
	Return the number of leptons in the event. More...
const int &	runnum () const
	Return a constant reference to the run number. More...
int &	runnum ()
	Return a reference to the run number. More...
bool	set_jet_types (const std::vector< int > &)
	Set the jet types in the event. More...
void	setMC (bool isMC)
	Set the Monte Carlo flag. More...
void	smear (CLHEP::HepRandomEngine &engine, bool smear_dir=false)
	Smear the objects in the event according to their resolutions. More...
void	sort ()
	Sort objects in the event according to their transverse momentum . More...
Fourvec	sum (int type) const
	Return the sum of all objects' four-momentum which have a particular type. More...
void	trimjets (std::vector< Lepjets_Event_Jet >::size_type n)
	Remove all but the first n jets. More...
double	zvertex () const
	Return the value of z-vertex. More...
double &	zvertex ()
	Return a reference to the value of z-vertex. More...

Private Attributes
double	_dlb
double	_dnn
int	_evnum
bool	_isMC
std::vector< Lepjets_Event_Jet >	_jets
Resolution	_kt_res
std::vector< Lepjets_Event_Lep >	_leps
Fourvec	_met
int	_runnum
double	_zvertex

Detailed Description

Represent a simple event consisting of lepton(s) and jet(s). An instance of this class holds a list of leptons (as represented by the Lepjets_Event_Lep class) and jets (as represented by Lepjets_Event_Jet class). Also recorded are:

• Missing transverse energy.
• vertex (Irrelevant for non-D0 experiment)
• Run and event number (Irrelevant for non-D0 experiment)

Definition at line 62 of file Lepjets_Event.h.

Constructor & Destructor Documentation

hitfit::Lepjets_Event::Lepjets_Event	(	int	runnum,
		int	evnum
	)

Constructor.

Parameters

runnum	The run number.
evnum	The event number.

Definition at line 46 of file Lepjets_Event.cc.

      : _zvertex(0), _isMC(false), _runnum(runnum), _evnum(evnum), _dlb(-1), _dnn(-1) {}

Member Function Documentation

void hitfit::Lepjets_Event::add_jet

(

const Lepjets_Event_Jet &

jet

)

Add a new jet to the event.

Parameters

jet	The jet to be added.

Definition at line 361 of file Lepjets_Event.cc.

References _jets.

Referenced by hitfit::RunHitFit< AElectron, AMuon, AJet, AMet >::FitAllPermutation(), hitfit::gentop(), and hitfit::gentth().

  {
    _jets.push_back(jet);
  }

void hitfit::Lepjets_Event::add_lep

(

const Lepjets_Event_Lep &

lep

)

Add a new lepton to the event.

Parameters

lep	The lepton to be added.

Definition at line 350 of file Lepjets_Event.cc.

References _leps.

Referenced by hitfit::RunHitFit< AElectron, AMuon, AJet, AMet >::AddLepton(), hitfit::gentop(), and hitfit::gentth().

  {
    _leps.push_back(lep);
  }

int hitfit::Lepjets_Event::cut_jets	(	double	pt_cut,
		double	eta_cut
	)

Remove jets which fail transverse momentum and pseudorapidity cut.

Parameters

pt_cut	Remove jets which have transverse momentum less than this value, in GeV.
eta_cut	Remove jetss which have absolute pseudorapidity more than this value.

Definition at line 481 of file Lepjets_Event.cc.

References _jets.

  {
    _jets.erase(remove_if(_jets.begin(), _jets.end(), Lepjets_Event_Cutter(pt_cut, eta_cut)), _jets.end());
    return _jets.size();
  }

int hitfit::Lepjets_Event::cut_leps	(	double	pt_cut,
		double	eta_cut
	)

Remove leptons which fail transverse momentum and pseudorapidity cut.

Parameters

pt_cut	Remove leptons which have transverse momentum less than this value, in GeV.
eta_cut	Remove leptons which have absolute pseudorapidity more than this value.

Definition at line 465 of file Lepjets_Event.cc.

References _leps.

  {
    _leps.erase(remove_if(_leps.begin(), _leps.end(), Lepjets_Event_Cutter(pt_cut, eta_cut)), _leps.end());
    return _leps.size();
  }

double & hitfit::Lepjets_Event::dlb

(

)

Return a reference to the value of low-bias (LB) discriminant (Irrelevant for non-D0 experiment).

Definition at line 281 of file Lepjets_Event.cc.

References _dlb.

  {
    return _dlb;
  }

double hitfit::Lepjets_Event::dlb

(

)

const

Return the value of low-bias (LB) discriminant (Irrelevant for non-D0 experiment).

Definition at line 270 of file Lepjets_Event.cc.

References _dlb.

  {
    return _dlb;
  }

double & hitfit::Lepjets_Event::dnn

(

)

Return a reference to the value of neural network (NN) discriminant (Irrelevant for non-D0 experiment).

Definition at line 303 of file Lepjets_Event.cc.

References _dnn.

  {
    return _dnn;
  }

double hitfit::Lepjets_Event::dnn

(

)

const

Return a the value of neural network (NN) discriminant (Irrelevant for non-D0 experiment).

Definition at line 292 of file Lepjets_Event.cc.

References _dnn.

  {
    return _dnn;
  }

std::ostream & hitfit::Lepjets_Event::dump	(	std::ostream &	s,
		bool	full = false
	)		const

Print the content of this object.

Parameters

s	The output stream to which to write
full	If TRUE, print all information about this instance of Lepjets_Event. If FALSE, print partial information about this instance of Lepjets_Event.

Definition at line 510 of file Lepjets_Event.cc.

References _evnum, _jets, _leps, _met, _runnum, _zvertex, full, mps_fire::i, and alignCSCRings::s.

Referenced by hitfit::operator<<().

  {
    s << "Run: " << _runnum << "  Event: " << _evnum << "\n";
    s << "Leptons:\n";
    for (std::vector<Lepjets_Event_Lep>::size_type i = 0; i < _leps.size(); i++) {
      s << "  ";
      _leps[i].dump(s, full);
      s << "\n";
    }
    s << "Jets:\n";
    for (std::vector<Lepjets_Event_Jet>::size_type i = 0; i < _jets.size(); i++) {
      s << "  ";
      _jets[i].dump(s, full);
      s << "\n";
    }
    s << "Missing Et: " << _met << "\n";
    if (_zvertex != 0)
      s << "z-vertex: " << _zvertex << "\n";
    return s;
  }

const int & hitfit::Lepjets_Event::evnum

(

)

const

Return a constant reference to the event number.

Definition at line 78 of file Lepjets_Event.cc.

References _evnum.

  {
    return _evnum;
  }

int & hitfit::Lepjets_Event::evnum

(

)

Return a reference to the event number.

Definition at line 89 of file Lepjets_Event.cc.

References _evnum.

  {
    return _evnum;
  }

bool hitfit::Lepjets_Event::isMC

(

)

const

Return the Monte Carlo flag.

Definition at line 248 of file Lepjets_Event.cc.

References _isMC.

Referenced by setMC().

  {
    return _isMC;
  }

Lepjets_Event_Jet & hitfit::Lepjets_Event::jet

(

std::vector< Lepjets_Event_Jet >::size_type

i

)

Return a reference to jet at index position i.

Parameters

i	The jet index position.

Definition at line 137 of file Lepjets_Event.cc.

References _jets, and mps_fire::i.

Referenced by hitfit::Constrained_Top::Constrained_Top(), hitfit::Constrained_Z::Constrained_Z(), hitfit::Top_Fit_Args::constrainer_args(), jet_permutation(), jet_types(), hitfit::Gentop_Args::kt_res_str(), TtSemiLepHitFitProducer< LeptonCollection >::produce(), and set_jet_types().

  {
    assert(i < _jets.size());
    return _jets[i];
  }

const Lepjets_Event_Jet & hitfit::Lepjets_Event::jet

(

std::vector< Lepjets_Event_Jet >::size_type

i

)

const

Return a constant reference to jet at index position i.

Parameters

i	The jet index position.

Definition at line 167 of file Lepjets_Event.cc.

References _jets, and mps_fire::i.

  {
    assert(i < _jets.size());
    return _jets[i];
  }

std::string hitfit::Lepjets_Event::jet_permutation

(

)

const

Return a string representing the jet permutation. The following notation is used for each type of jet:

• g ISR/gluon.
• b leptonic quark.
• B hadronic quark.
• w hadronic jet from boson.
• H jet from Higgs boson.
• ? Unknown.

Definition at line 541 of file Lepjets_Event.cc.

References _jets, jet(), hitfit::jetTypeString(), and AlCaHLTBitMon_QueryRunRegistry::string.

  {
    std::string permutation;
    for (size_t jet = 0; jet != _jets.size(); ++jet) {
      permutation = permutation + hitfit::jetTypeString(_jets[jet].type());
    }
    return permutation;
  }

std::vector< int > hitfit::Lepjets_Event::jet_types

(

)

const

Return the jet types in the event.

Definition at line 408 of file Lepjets_Event.cc.

References _jets, jet(), and runTheMatrix::ret.

Referenced by hitfit::Fit_Result::jet_types().

  {
    std::vector<int> ret;
    for (std::vector<Lepjets_Event_Jet>::size_type ijet = 0; ijet != _jets.size(); ijet++) {
      ret.push_back(jet(ijet).type());
    }
    return ret;
  }

Fourvec hitfit::Lepjets_Event::kt

(

)

const

Return the sum of all objects' four-momentum and missing transverse energy.

Definition at line 335 of file Lepjets_Event.cc.

References _jets, _leps, _met, mps_fire::i, AlCaHLTBitMon_ParallelJobs::p, and findQualityFiles::v.

Referenced by hitfit::Constrained_Top::Constrained_Top(), hitfit::Constrained_Z::Constrained_Z(), and smear().

  {
    Fourvec v = _met;
    for (std::vector<Lepjets_Event_Lep>::size_type i = 0; i < _leps.size(); i++)
      v += _leps[i].p();
    for (std::vector<Lepjets_Event_Jet>::size_type i = 0; i < _jets.size(); i++)
      v += _jets[i].p();
    return v;
  }

Resolution & hitfit::Lepjets_Event::kt_res

(

)

Return a reference to the resolution.

Definition at line 204 of file Lepjets_Event.cc.

References _kt_res.

Referenced by hitfit::Constrained_Top::Constrained_Top(), hitfit::Constrained_Z::Constrained_Z(), hitfit::gentop(), hitfit::gentth(), hitfit::RunHitFit< AElectron, AMuon, AJet, AMet >::SetKtResolution(), and hitfit::RunHitFit< AElectron, AMuon, AJet, AMet >::SetMet().

  {
    return _kt_res;
  }

const Resolution & hitfit::Lepjets_Event::kt_res

(

)

const

Return a const reference to the resolution.

Definition at line 215 of file Lepjets_Event.cc.

References _kt_res.

  {
    return _kt_res;
  }

Lepjets_Event_Lep & hitfit::Lepjets_Event::lep

(

std::vector< Lepjets_Event_Lep >::size_type

i

)

Return a reference to lepton at index position i.

Parameters

i	The lepton index position.

Definition at line 122 of file Lepjets_Event.cc.

References _leps, and mps_fire::i.

Referenced by hitfit::Constrained_Top::Constrained_Top(), hitfit::Constrained_Z::Constrained_Z(), and TtSemiLepHitFitProducer< LeptonCollection >::produce().

  {
    assert(i < _leps.size());
    return _leps[i];
  }

const Lepjets_Event_Lep & hitfit::Lepjets_Event::lep

(

std::vector< Lepjets_Event_Lep >::size_type

i

)

const

Return a constant reference to lepton at index position i.

Parameters

i	The lepton index position.

Definition at line 152 of file Lepjets_Event.cc.

References _leps, and mps_fire::i.

  {
    assert(i < _leps.size());
    return _leps[i];
  }

Fourvec & hitfit::Lepjets_Event::met

(

)

Return a reference to the missing transverse energy.

Definition at line 182 of file Lepjets_Event.cc.

References _met.

Referenced by hitfit::Constrained_Top::Constrained_Top(), hitfit::Constrained_Z::Constrained_Z(), hitfit::Top_Fit::fit_one_perm(), hitfit::gentop(), hitfit::gentth(), hitfit::Top_Decaykin::lepw(), TtSemiLepHitFitProducer< LeptonCollection >::produce(), hitfit::RunHitFit< AElectron, AMuon, AJet, AMet >::SetMet(), hitfit::Top_Decaykin::solve_nu(), and hitfit::Top_Decaykin::solve_nu_tmass().

  {
    return _met;
  }

const Fourvec & hitfit::Lepjets_Event::met

(

)

const

Return a constant reference to the missing transverse energy.

Definition at line 193 of file Lepjets_Event.cc.

References _met.

  {
    return _met;
  }

std::vector< Lepjets_Event_Jet >::size_type hitfit::Lepjets_Event::njets

(

)

const

Return the number of jets in the event.

Definition at line 111 of file Lepjets_Event.cc.

References _jets.

Referenced by hitfit::Constrained_Top::Constrained_Top(), hitfit::Constrained_Z::Constrained_Z(), hitfit::Top_Fit_Args::constrainer_args(), hitfit::Top_Fit::fit(), hitfit::Gentop_Args::kt_res_str(), and set_jet_types().

  {
    return _jets.size();
  }

std::vector< Lepjets_Event_Lep >::size_type hitfit::Lepjets_Event::nleps

(

)

const

Return the number of leptons in the event.

Definition at line 100 of file Lepjets_Event.cc.

References _leps.

Referenced by hitfit::Constrained_Top::Constrained_Top(), and hitfit::Constrained_Z::Constrained_Z().

  {
    return _leps.size();
  }

const int & hitfit::Lepjets_Event::runnum

(

)

const

Return a constant reference to the run number.

Definition at line 56 of file Lepjets_Event.cc.

References _runnum.

  {
    return _runnum;
  }

int & hitfit::Lepjets_Event::runnum

(

)

Return a reference to the run number.

Definition at line 67 of file Lepjets_Event.cc.

References _runnum.

  {
    return _runnum;
  }

bool hitfit::Lepjets_Event::set_jet_types

(

const std::vector< int > &

_jet_types

)

Set the jet types in the event.

Definition at line 420 of file Lepjets_Event.cc.

References _eta_cut, _jets, _pt_cut, hitfit::hadw1_label, hitfit::hadw2_label, mps_fire::i, jet(), cmsLHEtoEOSManager::l, njets(), OrderedSet::t, and hitfit::Lepjets_Event_Lep::type().

Referenced by hitfit::RunHitFit< AElectron, AMuon, AJet, AMet >::FitAllPermutation().

  {
    if (_jets.size() != _jet_types.size()) {
      return false;
    }
    bool saw_hadw1 = false;
    for (std::vector<Lepjets_Event_Jet>::size_type i = 0; i < njets(); i++) {
      int t = _jet_types[i];
      if (t == hadw1_label) {
        if (saw_hadw1)
          t = hadw2_label;
        saw_hadw1 = true;
      }
      jet(i).type() = t;
    }
    return true;
  }

void hitfit::Lepjets_Event::setMC

(

bool

isMC

)

Set the Monte Carlo flag.

Definition at line 259 of file Lepjets_Event.cc.

References _isMC, and isMC().

  {
    _isMC = isMC;
  }

void hitfit::Lepjets_Event::smear	(	CLHEP::HepRandomEngine &	engine,
		bool	smear_dir = false
	)

Smear the objects in the event according to their resolutions.

Parameters

engine	The underlying random number generator.
smear_dir	If TRUE, also smear the object's direction. If FALSE, then only smear the magnitude of three-momentum.

Definition at line 372 of file Lepjets_Event.cc.

References _jets, _kt_res, _leps, _met, mps_fire::i, kt(), hitfit::Resolution::pick(), X, and DOFs::Y.

Referenced by hitfit::gentop(), and hitfit::gentth().

  {
    Fourvec before, after;
    for (std::vector<Lepjets_Event_Lep>::size_type i = 0; i < _leps.size(); i++) {
      before += _leps[i].p();
      _leps[i].smear(engine, smear_dir);
      after += _leps[i].p();
    }
    for (std::vector<Lepjets_Event_Jet>::size_type i = 0; i < _jets.size(); i++) {
      before += _jets[i].p();
      _jets[i].smear(engine, smear_dir);
      after += _jets[i].p();
    }

    Fourvec kt = _met + before;
    kt(Fourvec::X) = _kt_res.pick(kt(Fourvec::X), kt(Fourvec::X), engine);
    kt(Fourvec::Y) = _kt_res.pick(kt(Fourvec::Y), kt(Fourvec::Y), engine);
    _met = kt - after;
  }

void hitfit::Lepjets_Event::sort

(

)

Sort objects in the event according to their transverse momentum .

Definition at line 399 of file Lepjets_Event.cc.

References _jets, and _leps.

  {
    std::stable_sort(_leps.begin(), _leps.end(), std::not_fn(std::less<Lepjets_Event_Lep>()));
    std::stable_sort(_jets.begin(), _jets.end(), std::not_fn(std::less<Lepjets_Event_Lep>()));
  }

Fourvec hitfit::Lepjets_Event::sum

(

int

type

)

const

Return the sum of all objects' four-momentum which have a particular type.

Parameters

type	The type code of the objects to be summed up.

Definition at line 314 of file Lepjets_Event.cc.

References _jets, _leps, mps_fire::i, MillePedeFileConverter_cfg::out, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by hitfit::Top_Fit_Args::constrainer_args(), hitfit::Top_Decaykin::hadt(), hitfit::Top_Decaykin::hadw(), hitfit::Top_Decaykin::hadw1(), hitfit::Top_Decaykin::hadw2(), hitfit::Top_Decaykin::lept(), and hitfit::Top_Decaykin::solve_nu_tmass().

  {
    Fourvec out;
    for (std::vector<Lepjets_Event_Lep>::size_type i = 0; i < _leps.size(); i++)
      if (_leps[i].type() == type)
        out += _leps[i].p();
    for (std::vector<Lepjets_Event_Jet>::size_type i = 0; i < _jets.size(); i++)
      if (_jets[i].type() == type)
        out += _jets[i].p();
    return out;
  }

void hitfit::Lepjets_Event::trimjets

(

std::vector< Lepjets_Event_Jet >::size_type

n

)

Remove all but the first n jets.

Parameters

n	The number of jets to keep.

Definition at line 497 of file Lepjets_Event.cc.

References _jets, and dqmiodumpmetadata::n.

  {
    if (n >= _jets.size())
      return;
    _jets.erase(_jets.begin() + n, _jets.end());
  }

double hitfit::Lepjets_Event::zvertex

(

)

const

Return the value of z-vertex.

Definition at line 226 of file Lepjets_Event.cc.

References _zvertex.

  {
    return _zvertex;
  }

double & hitfit::Lepjets_Event::zvertex

(

)

Return a reference to the value of z-vertex.

Definition at line 237 of file Lepjets_Event.cc.

References _zvertex.

  {
    return _zvertex;
  }

Member Data Documentation

double hitfit::Lepjets_Event::_dlb

private

The low-bias (LB) discriminant.

Definition at line 378 of file Lepjets_Event.h.

Referenced by dlb().

double hitfit::Lepjets_Event::_dnn

private

The neural network (NN) discriminant.

Definition at line 383 of file Lepjets_Event.h.

Referenced by dnn().

int hitfit::Lepjets_Event::_evnum

private

The event number.

Definition at line 373 of file Lepjets_Event.h.

Referenced by dump(), and evnum().

bool hitfit::Lepjets_Event::_isMC

private

The Monte Calro flag.

Definition at line 363 of file Lepjets_Event.h.

Referenced by isMC(), and setMC().

std::vector<Lepjets_Event_Jet> hitfit::Lepjets_Event::_jets

private

The list of jets in the event.

Definition at line 342 of file Lepjets_Event.h.

Referenced by add_jet(), cut_jets(), dump(), jet(), jet_permutation(), jet_types(), kt(), njets(), set_jet_types(), smear(), sort(), sum(), and trimjets().

Resolution hitfit::Lepjets_Event::_kt_res

private

The resolution.

Definition at line 353 of file Lepjets_Event.h.

Referenced by kt_res(), and smear().

std::vector<Lepjets_Event_Lep> hitfit::Lepjets_Event::_leps

private

The list of leptons in the event.

Definition at line 337 of file Lepjets_Event.h.

Referenced by add_lep(), cut_leps(), dump(), kt(), lep(), nleps(), smear(), sort(), and sum().

Fourvec hitfit::Lepjets_Event::_met

private

Missing transverse energy.

Definition at line 348 of file Lepjets_Event.h.

Referenced by dump(), kt(), met(), and smear().

int hitfit::Lepjets_Event::_runnum

private

The run number.

Definition at line 368 of file Lepjets_Event.h.

Referenced by dump(), and runnum().

double hitfit::Lepjets_Event::_zvertex

private

The vertex of the event.

Definition at line 358 of file Lepjets_Event.h.

Referenced by dump(), and zvertex().