hitfit::Constrained_Top Class Reference

Do a constrained kinematic fitting for a event. More...

#include <Constrained_Top.h>

Public Member Functions
double	constrain (Lepjets_Event &ev, double &mt, double &sigmt, Column_Vector &pullx, Column_Vector &pully)
	Do a constrained fit of events. Returns the top mass and its error in mt and sigmt, and the pull quantities in pullx and pully. Returns the , this will be negative if the fit failed to converge. More...
	Constrained_Top (const Constrained_Top_Args &args, double lepw_mass, double hadw_mass, double top_mass)
	Constructor, create an instance of the Constrained_Top object from the arguments object and the mass constraints. More...

Private Attributes
const Constrained_Top_Args	_args
Fourvec_Constrainer	_constrainer

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

Detailed Description

Do a constrained kinematic fitting for a event.

Definition at line 128 of file Constrained_Top.h.

Constructor & Destructor Documentation

hitfit::Constrained_Top::Constrained_Top	(	const Constrained_Top_Args &	args,
		double	lepw_mass,
		double	hadw_mass,
		double	top_mass
	)

Constructor, create an instance of the Constrained_Top object from the arguments object and the mass constraints.

Parameters

args	Argument for this instance of Constrained_Top object.
lepw_mass	The mass to which the leptonic W should be constrained. If this parameter is set to 0, the constraint is skied.
hadw_mass	The mass to which the hadronic W should be constrained. If this parameter is set to 0, the constraint is skied.
top_mass	The mass to which the top quarks should be constrained. If this parameter is set to 0, the constraints is skied.

Definition at line 106 of file Constrained_Top.cc.

References _constrainer, hitfit::Fourvec_Event::add(), hitfit::Fourvec_Constrainer::add_constraint(), hitfit::Constrained_Top_Args::equal_side(), hitfit::Lepjets_Event_Lep::eta_sigma(), ev, hitfit::hadb_label, hitfit::hadw1_label, hitfit::hadw2_label, hitfit::higgs_label, hitfit::Resolution::inverse(), hitfit::isr_label, hitfit::Lepjets_Event::jet(), gen::k, hitfit::Lepjets_Event::kt(), hitfit::Lepjets_Event::kt_res(), hitfit::Lepjets_Event::lep(), hitfit::lepb_label, hitfit::lepton_label, ResonanceBuilder::mass, hitfit::Fourvec_Constrainer::mass_constraint(), hitfit::Lepjets_Event::met(), hitfit::Lepjets_Event::njets(), hitfit::Lepjets_Event::nleps(), hitfit::Fourvec_Event::nu(), hitfit::nu_label, hgcalPlots::obj, hitfit::Fourvec_Event::obj(), hitfit::Lepjets_Event_Lep::p(), hitfit::FE_Obj::p, hitfit::Vector_Resolution::p_res(), hitfit::Lepjets_Event_Lep::p_sigma(), hitfit::Lepjets_Event_Lep::phi_sigma(), hitfit::Lepjets_Event_Lep::res(), hitfit::Fourvec_Event::set_kt_error(), hitfit::Fourvec_Event::set_nu_p(), hitfit::Resolution::sigma(), and hitfit::Lepjets_Event_Lep::type().

  : _args (args),
    _constrainer (args.fourvec_constrainer_args())
{
  char buf[256];
  if (lepw_mass > 0) {
    sprintf (buf, "(%d %d) = %f", nu_label, lepton_label, lepw_mass);
    _constrainer.add_constraint (buf);
  }

  if (hadw_mass > 0) {
    sprintf (buf, "(%d %d) = %f", hadw1_label, hadw2_label, hadw_mass);
    _constrainer.add_constraint (buf);
  }

  if (args.equal_side()) {
      sprintf (buf, "(%d %d %d) = (%d %d %d)",
               nu_label, lepton_label, lepb_label,
               hadw1_label, hadw2_label, hadb_label);
      _constrainer.add_constraint (buf);
  }

  if (top_mass > 0) {
    sprintf (buf, "(%d %d %d) = %f",
             hadw1_label, hadw2_label, hadb_label, top_mass);
    _constrainer.add_constraint (buf);
  }
  else {
    sprintf (buf, "(%d %d %d) = 0", hadw1_label, hadw2_label, hadb_label);
    _constrainer.mass_constraint (buf);
  }
}

Member Function Documentation

double hitfit::Constrained_Top::constrain	(	Lepjets_Event &	ev,
		double &	mt,
		double &	sigmt,
		Column_Vector &	pullx,
		Column_Vector &	pully
	)

Do a constrained fit of events. Returns the top mass and its error in mt and sigmt, and the pull quantities in pullx and pully. Returns the , this will be negative if the fit failed to converge.

Parameters

ev	The event to be fitted (input), and the event after fitting (output).
mt	The fitted top mass.
sigmt	The uncertainty on the fitted top mass.
pullx	Pull quantities for the well-measured variables.
pully	Pull quantities for the poorly-measured variables.

Input:

• Lepjets_Event ev.

Output:

• Lepjets_Event ev.
• double mt.
• double sigmt.
• Column_Vector pullx.
• Column_Vector pully.

Return:

The of the fit. Return a negative value if the fit didn't converge.

Definition at line 290 of file Constrained_Top.cc.

References _args, _constrainer, hitfit::Constrained_Top_Args::bmass(), and hitfit::Fourvec_Constrainer::constrain().

Referenced by hitfit::Top_Fit::fit_one_perm().

{
  Fourvec_Event fe;
  do_import (ev, _args.bmass (), fe);
  double chisq = _constrainer.constrain (fe, mt, sigmt, pullx, pully);
  do_export (fe, ev);

  return chisq;
}

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  s, const Constrained_Top &  ct  )

friend

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

Parameters

s	The output stream to which to write.
ct	The instance of Constrained_Top to be printed.

Definition at line 333 of file Constrained_Top.cc.

{
  return s << ct._constrainer;
}

Member Data Documentation

const Constrained_Top_Args hitfit::Constrained_Top::_args

private

Parameter settings for the constrainer.

Definition at line 208 of file Constrained_Top.h.

Referenced by Vispa.Main.Application.Application::_readCommandLineAttributes(), and constrain().

Fourvec_Constrainer hitfit::Constrained_Top::_constrainer

private

The guy that actually does the work.

Definition at line 214 of file Constrained_Top.h.

Referenced by constrain(), Constrained_Top(), and hitfit::operator<<().