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 $k_{T}$ four-momentum. More...

double	kt_x_error () const
	Return the x uncertainty in $k_{T}$ . More...

double	kt_xy_covar () const
	Return the xy covariance in $k_{T}$ . More...

double	kt_y_error () const
	Return the y uncertainty in $k_{T}$ . 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 $k_{T}$ . 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 $x$ momentum, that will be added into the $k_{T}$ 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.

106 : _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.

   {
     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;
   }

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().

◆ 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.

   {
     return _has_neutrino;
   }

References _has_neutrino.

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

◆ kt()

const Fourvec & hitfit::Fourvec_Event::kt ( ) const

Access the $k_{T}$ four-momentum.

Definition at line 163 of file Fourvec_Event.cc.

   {
     return _kt;
   }

References _kt.

◆ kt_x_error()

double hitfit::Fourvec_Event::kt_x_error ( ) const

Return the x uncertainty in $k_{T}$ .

Definition at line 179 of file Fourvec_Event.cc.

   {
     return _kt_x_error;
   }

References _kt_x_error.

Referenced by set_kt_error().

◆ kt_xy_covar()

double hitfit::Fourvec_Event::kt_xy_covar ( ) const

Return the xy covariance in $k_{T}$ .

Definition at line 201 of file Fourvec_Event.cc.

   {
     return _kt_xy_covar;
   }

References _kt_xy_covar.

Referenced by set_kt_error().

◆ kt_y_error()

double hitfit::Fourvec_Event::kt_y_error ( ) const

Return the y uncertainty in $k_{T}$ .

Definition at line 190 of file Fourvec_Event.cc.

   {
     return _kt_y_error;
   }

References _kt_y_error.

Referenced by set_kt_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.

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

References _has_neutrino, and _objs.

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

◆ 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.

   {
     return _objs.size();
   }

References _objs.

◆ nu()

const Fourvec & hitfit::Fourvec_Event::nu ( ) const

Access the neutrino four-momentum.

Definition at line 154 of file Fourvec_Event.cc.

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

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

◆ 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.

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

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

Referenced by add().

◆ 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 $k_{T}$ .

Parameters

kt_x_error	The uncertainty in the component of $k_{T}$ .
kt_y_error	The uncertainty in the component of $k_{T}$ .
kt_xy_covar	The covariance between the and component of $k_{T}$ .

Definition at line 310 of file Fourvec_Event.cc.

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

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

◆ set_nu_p()

void hitfit::Fourvec_Event::set_nu_p ( const Fourvec & p )

Set the neutrino four-momentum to $p$ . 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.

   {
     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;
   }

References _has_neutrino, _kt, _objs, hitfit::nu_label, and AlCaHLTBitMon_ParallelJobs::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 $p$ .

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.

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

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

◆ set_x_p()

void hitfit::Fourvec_Event::set_x_p ( const Fourvec & p )

Set the four-momentum of the $x$ object.

Parameters

p	The new four-momentum of the object.

Definition at line 297 of file Fourvec_Event.cc.

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

References _kt, _x, and AlCaHLTBitMon_ParallelJobs::p.

◆ x()

const Fourvec & hitfit::Fourvec_Event::x ( ) const

Access the $x$ four-momentum.

Definition at line 171 of file Fourvec_Event.cc.

   {
     return _x;
   }

References _x.

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

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 $k_{T}$ , this should always be equal to the sum of all the object momenta, including $x$ .

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 $x-$ component of $k_{T}$ .

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 $x-$ and $y-$ component of $k_{T}$ .

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 $y-$ component of $k_{T}$ .

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 $x$ object.

Definition at line 351 of file Fourvec_Event.h.

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

Public Member Functions

Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ Fourvec_Event()

Member Function Documentation

◆ add()

◆ has_neutrino()

◆ kt()

◆ kt_x_error()

◆ kt_xy_covar()

◆ kt_y_error()

◆ nobjs()

◆ nobjs_all()

◆ nu()

◆ obj()

◆ set_kt_error()

◆ set_nu_p()

◆ set_obj_p()

◆ set_x_p()

◆ x()

Friends And Related Function Documentation

◆ operator<<

Member Data Documentation

◆ _has_neutrino

◆ _kt

◆ _kt_x_error

◆ _kt_xy_covar

◆ _kt_y_error

◆ _objs

◆ _x