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 66 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 50 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 397 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 385 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 533 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 514 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 311 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 299 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 335 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 323 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 566 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 92 of file Lepjets_Event.cc.

References _evnum.

{
  return _evnum;
}

int & hitfit::Lepjets_Event::evnum

(

)

Return a reference to the event number.

Definition at line 104 of file Lepjets_Event.cc.

References _evnum.

{
  return _evnum;
}

bool hitfit::Lepjets_Event::isMC

(

)

const

Return the Monte Carlo flag.

Definition at line 276 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 156 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 188 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 598 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 447 of file Lepjets_Event.cc.

References _jets, and jet().

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 369 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 228 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 240 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 140 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 172 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 204 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 216 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 128 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 116 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 68 of file Lepjets_Event.cc.

References _runnum.

{
  return _runnum;
}

int & hitfit::Lepjets_Event::runnum

(

)

Return a reference to the run number.

Definition at line 80 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 462 of file Lepjets_Event.cc.

References _eta_cut, _jets, _pt_cut, hitfit::hadw1_label, hitfit::hadw2_label, mps_fire::i, jet(), checklumidiff::l, njets(), protons_cff::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 288 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 409 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 437 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 347 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 552 of file Lepjets_Event.cc.

References _jets, and gen::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 252 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 264 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 383 of file Lepjets_Event.h.

Referenced by dlb().

double hitfit::Lepjets_Event::_dnn

private

The neural network (NN) discriminant.

Definition at line 388 of file Lepjets_Event.h.

Referenced by dnn().

int hitfit::Lepjets_Event::_evnum

private

The event number.

Definition at line 378 of file Lepjets_Event.h.

Referenced by dump(), and evnum().

bool hitfit::Lepjets_Event::_isMC

private

The Monte Calro flag.

Definition at line 368 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 347 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 358 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 342 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 353 of file Lepjets_Event.h.

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

int hitfit::Lepjets_Event::_runnum

private

The run number.

Definition at line 373 of file Lepjets_Event.h.

Referenced by dump(), and runnum().

double hitfit::Lepjets_Event::_zvertex

private

The vertex of the event.

Definition at line 363 of file Lepjets_Event.h.

Referenced by dump(), and zvertex().