hitfit::Fourvec_Event Class Reference

Represent an event for kinematic fitting as a collection of four-momenta. Each object is represented as an instance of FE_Obj. There may be an object for a neutrino. If that is the case, it is always at the end of the object list. It is not included in the count returned by nobjs(). But is is included in nobjs_all(). More...

#include <Fourvec_Event.h>

Public Member Functions
void	add (const FE_Obj &obj)
	Add an object to the event. The object should not be a neutrino, use the method set_nu_p for that. More...
	Fourvec_Event ()
	Default constructor. More...
bool	has_neutrino () const
	Return TRUE is this event contains a neutrino, otherwise returns FALSE. More...
const Fourvec &	kt () const
	Access the four-momentum. More...
double	kt_x_error () const
	Return the x uncertainty in . More...
double	kt_xy_covar () const
	Return the xy covariance in . More...
double	kt_y_error () const
	Return the y uncertainty in . More...
int	nobjs () const
	Return the number of objects in the event not including any neutrinos. More...
int	nobjs_all () const
	Return the number of objects in the event including any neutrinos. More...
const Fourvec &	nu () const
	Access the neutrino four-momentum. More...
const FE_Obj &	obj (std::vector< FE_Obj >::size_type i) const
	Access object at index i, with the convention that the index starts at 0. More...
void	set_kt_error (double kt_x_error, double kt_y_error, double kt_xy_covar)
	Set the uncertainties on . More...
void	set_nu_p (const Fourvec &p)
	Set the neutrino four-momentum to . This method adds a neutrino if there wasn't already one. More...
void	set_obj_p (std::vector< FE_Obj >::size_type i, const Fourvec &p)
	Set the four-momentum of object at index i to . More...
void	set_x_p (const Fourvec &p)
	Set the four-momentum of the object. More...
const Fourvec &	x () const
	Access the four-momentum. More...

Private Attributes
bool	_has_neutrino
Fourvec	_kt
double	_kt_x_error
double	_kt_xy_covar
double	_kt_y_error
std::vector< FE_Obj >	_objs
Fourvec	_x

Friends
std::ostream &	operator<< (std::ostream &s, const Fourvec_Event &fe)
	Output stream operator, print the content of this Fourvec_Event object to an output stream. More...

Detailed Description

Represent an event for kinematic fitting as a collection of four-momenta. Each object is represented as an instance of FE_Obj. There may be an object for a neutrino. If that is the case, it is always at the end of the object list. It is not included in the count returned by nobjs(). But is is included in nobjs_all().

We can also record the other momentum, that will be added into the sum. This can be used to store a missing transverse energy that is not attributed to a neutrino but is instead due to mismeasurement. Typically this will be set to zero in events that have a neutrino, and to the measured missing transverse energy in events that do not.

Definition at line 193 of file Fourvec_Event.h.

Constructor & Destructor Documentation

Fourvec_Event()

hitfit::Fourvec_Event::Fourvec_Event

(

)

Default constructor.

Definition at line 102 of file Fourvec_Event.cc.

      : _kt_x_error(0), _kt_y_error(0), _kt_xy_covar(0), _has_neutrino(false) {}

Member Function Documentation

add()

void hitfit::Fourvec_Event::add

(

const FE_Obj &

obj

)

Add an object to the event. The object should not be a neutrino, use the method set_nu_p for that.

Parameters

obj	The FE_Obj to add.

Definition at line 240 of file Fourvec_Event.cc.

References _has_neutrino, _kt, _objs, cms::cuda::assert(), hitfit::FE_Obj::label, hitfit::nu_label, obj(), and hitfit::FE_Obj::p.

Referenced by counter.Counter::register().

  {
    assert(obj.label != nu_label);

    // Add to the end of the list, but before any neutrino.
    if (_has_neutrino) {
      assert(!_objs.empty());
      _objs.insert(_objs.begin() + _objs.size() - 1, obj);
    } else
      _objs.push_back(obj);

    // Maintain kt.
    _kt += obj.p;
  }

has_neutrino()

bool hitfit::Fourvec_Event::has_neutrino

(

)

const

Return TRUE is this event contains a neutrino, otherwise returns FALSE.

Definition at line 108 of file Fourvec_Event.cc.

References _has_neutrino.

Referenced by hitfit::Fourvec_Constraint_Calculator::eval().

  {
    return _has_neutrino;
  }

kt()

const Fourvec & hitfit::Fourvec_Event::kt

(

)

const

Access the four-momentum.

Definition at line 163 of file Fourvec_Event.cc.

References _kt.

  {
    return _kt;
  }

kt_x_error()

double hitfit::Fourvec_Event::kt_x_error

(

)

const

Return the x uncertainty in .

Definition at line 179 of file Fourvec_Event.cc.

References _kt_x_error.

Referenced by set_kt_error().

  {
    return _kt_x_error;
  }

kt_xy_covar()

double hitfit::Fourvec_Event::kt_xy_covar

(

)

const

Return the xy covariance in .

Definition at line 201 of file Fourvec_Event.cc.

References _kt_xy_covar.

Referenced by set_kt_error().

  {
    return _kt_xy_covar;
  }

kt_y_error()

double hitfit::Fourvec_Event::kt_y_error

(

)

const

Return the y uncertainty in .

Definition at line 190 of file Fourvec_Event.cc.

References _kt_y_error.

Referenced by set_kt_error().

  {
    return _kt_y_error;
  }

nobjs()

int hitfit::Fourvec_Event::nobjs

(

)

const

Return the number of objects in the event not including any neutrinos.

Definition at line 119 of file Fourvec_Event.cc.

References _has_neutrino, and _objs.

Referenced by hitfit::Fourvec_Constraint_Calculator::eval().

  {
    return _objs.size() - (_has_neutrino ? 1 : 0);
  }

nobjs_all()

int hitfit::Fourvec_Event::nobjs_all

(

)

const

Return the number of objects in the event including any neutrinos.

Definition at line 131 of file Fourvec_Event.cc.

References _objs.

  {
    return _objs.size();
  }

nu()

const Fourvec & hitfit::Fourvec_Event::nu

(

)

const

Access the neutrino four-momentum.

Definition at line 154 of file Fourvec_Event.cc.

References _has_neutrino, _objs, and cms::cuda::assert().

  {
    assert(_has_neutrino);
    return _objs.back().p;
  }

obj()

const FE_Obj & hitfit::Fourvec_Event::obj

(

std::vector< FE_Obj >::size_type

i

)

const

Access object at index i, with the convention that the index starts at 0.

Parameters

i	The index of the desired object.

Definition at line 142 of file Fourvec_Event.cc.

References cms::cuda::assert(), and mps_fire::i.

Referenced by add().

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

set_kt_error()

void hitfit::Fourvec_Event::set_kt_error	(	double	kt_x_error,
		double	kt_y_error,
		double	kt_xy_covar
	)

Set the uncertainties on .

Parameters

kt_x_error	The uncertainty in the component of .
kt_y_error	The uncertainty in the component of .
kt_xy_covar	The covariance between the and component of .

Definition at line 310 of file Fourvec_Event.cc.

References _kt_x_error, _kt_xy_covar, _kt_y_error, kt_x_error(), kt_xy_covar(), and kt_y_error().

  {
    _kt_x_error = kt_x_error;
    _kt_y_error = kt_y_error;
    _kt_xy_covar = kt_xy_covar;
  }

set_nu_p()

void hitfit::Fourvec_Event::set_nu_p

(

const Fourvec &

p

)

Set the neutrino four-momentum to . This method adds a neutrino if there wasn't already one.

Parameters

p	The new four-momentum of the neutrino.

Definition at line 262 of file Fourvec_Event.cc.

References _has_neutrino, _kt, _objs, hitfit::nu_label, and AlCaHLTBitMon_ParallelJobs::p.

  {
    if (_has_neutrino) {
      _kt -= _objs.back().p;
      _objs.back().p = p;
    } else {
      _has_neutrino = true;
      _objs.push_back(FE_Obj(p, 0, nu_label, 0, 0, 0, false));
    }

    _kt += p;
  }

set_obj_p()

void hitfit::Fourvec_Event::set_obj_p	(	std::vector< FE_Obj >::size_type	i,
		const Fourvec &	p
	)

Set the four-momentum of object at index i to .

Parameters

i	The position index of the object to change.
p	The new four-momentum of object at index i.

Definition at line 282 of file Fourvec_Event.cc.

References _kt, _objs, cms::cuda::assert(), mps_fire::i, and AlCaHLTBitMon_ParallelJobs::p.

  {
    assert(i < _objs.size());
    _kt -= _objs[i].p;
    _objs[i].p = p;
    _kt += p;
  }

set_x_p()

void hitfit::Fourvec_Event::set_x_p

(

const Fourvec &

p

)

Set the four-momentum of the object.

Parameters

p	The new four-momentum of the object.

Definition at line 297 of file Fourvec_Event.cc.

References _kt, _x, and AlCaHLTBitMon_ParallelJobs::p.

  {
    _kt -= _x;
    _x = p;
    _kt += p;
  }

x()

const Fourvec & hitfit::Fourvec_Event::x

(

)

const

Access the four-momentum.

Definition at line 171 of file Fourvec_Event.cc.

References _x.

Referenced by svgfig.Curve.Sample::__repr__(), svgfig.Ellipse::__repr__(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().

  {
    return _x;
  }

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  s, const Fourvec_Event &  fe  )

friend

Output stream operator, print the content of this Fourvec_Event object to an output stream.

Parameters

s	The stream to which to write.
fe	The instance of Fourvec_Event to be printed.

Definition at line 220 of file Fourvec_Event.cc.

  {
    s << "kt: (" << fe._kt.x() << ", " << fe._kt.y() << "); "
      << " error: " << fe._kt_x_error << " " << fe._kt_y_error << " " << fe._kt_xy_covar << "\n";
    s << "x: " << fe._x << "\n";
    for (unsigned i = 0; i < fe._objs.size(); i++)
      s << i + 1 << ": " << fe._objs[i];
    return s;
  }

Member Data Documentation

_has_neutrino

bool hitfit::Fourvec_Event::_has_neutrino

private

Flag that a neutrino has been added to the event.

Definition at line 372 of file Fourvec_Event.h.

Referenced by add(), has_neutrino(), nobjs(), nu(), and set_nu_p().

_kt

Fourvec hitfit::Fourvec_Event::_kt

private

Cached , this should always be equal to the sum of all the object momenta, including .

Definition at line 345 of file Fourvec_Event.h.

Referenced by add(), kt(), hitfit::operator<<(), set_nu_p(), set_obj_p(), and set_x_p().

_kt_x_error

double hitfit::Fourvec_Event::_kt_x_error

private

The uncertainty in the component of .

Definition at line 356 of file Fourvec_Event.h.

Referenced by kt_x_error(), hitfit::operator<<(), and set_kt_error().

_kt_xy_covar

double hitfit::Fourvec_Event::_kt_xy_covar

private

The covariance between the and component of .

Definition at line 366 of file Fourvec_Event.h.

Referenced by kt_xy_covar(), hitfit::operator<<(), and set_kt_error().

_kt_y_error

double hitfit::Fourvec_Event::_kt_y_error

private

The uncertainty in the component of .

Definition at line 360 of file Fourvec_Event.h.

Referenced by kt_y_error(), hitfit::operator<<(), and set_kt_error().

_objs

std::vector<FE_Obj> hitfit::Fourvec_Event::_objs

private

The list of contained objects in the event.

Definition at line 337 of file Fourvec_Event.h.

Referenced by add(), nobjs(), nobjs_all(), nu(), hitfit::operator<<(), set_nu_p(), and set_obj_p().

_x

Fourvec hitfit::Fourvec_Event::_x

private

Four-momentum of the object.

Definition at line 351 of file Fourvec_Event.h.

Referenced by plotting.PlotText::Draw(), hitfit::operator<<(), set_x_p(), and x().