Classes
class	Base_Constrainer
	Base class for $\chi^{2}$ constrained fitter. More...

class	Base_Constrainer_Args
	Hold on to parameters for the Base_Constrainer class. More...

class	Chisq_Constrainer
	Minimize a $\chi^{2}$ subject to a set of constraints. Based on the SQUAW algorithm. More...

class	Chisq_Constrainer_Args
	Hold on to parameters for the Chisq_Constrainer class. More...

class	Constrained_Top
	Do a constrained kinematic fitting for a $t\bar{t}\to\ell + \rm{jets}$ event. More...

class	Constrained_Top_Args
	Hold on to parameters for the Constrained_Top class. More...

class	Constrained_Z
	Do a constrained kinematic fitting for a $Z\to\ell^{+}\ell^{-} + \rm{jets}$ event. More...

class	Constrained_Z_Args
	Hold on to parameters for the Constrained_Z class. More...

class	Constraint
	Represent a mass constraint equation. Mass constraints come in two varieties, either saying that the sum of a set of labels should equal a constant: More...

class	Constraint_Calculator
	Abstract base class for evaluating constraints. Users derive from this class and implement the eval() method. More...

class	Constraint_Intermed
	Abstract base classes for describing one side of a mass constraint. More...

class	Constraint_Intermed_Constant
	Concrete class for one side of mass constraint equation of the type: . More...

class	Constraint_Intermed_Labels
	Concrete class for one side of mass constraint equation of the type: . More...

class	Defaults
	Define an interface for getting parameter settings. More...

class	Defaults_Text
	A lightweight implementation of the Defaults interface that uses simple ASCII text files. More...

class	Defaults_Textrep
	The internal representation for a Defaults_Text object. More...

class	EtaDepResElement
	Represent a resolution and an $\eta$ range in which the resolution is valid. More...

class	EtaDepResolution
	Hold on to $\eta$ -dependent resolution. This class acts as a function object and returns Vector_Resolution as a function of $\eta$ . It does not assume symmetry between $+\eta$ and $-\eta$ . The interpretation of $\eta$ as physics $\eta$ or detector $\eta$ is left to users's implementation. More...

class	FE_Obj
	Represent a single object in a Fourvec_Event, this is just a dumb data container. Each object in a Fourvec_Event has the following attributes: More...

class	Fit_Result
	Hold the result of one kinematic fit. More...

class	Fit_Result_Vec
	Holds pointers to a set of Fit_Result objects, resulting from different jet permutation with some consistent selection. The results are ordered by increasing $\chi^{2}$ values. A maximum length for the list of Fit_Result objects may be specified; when new results, those with the largest $\chi^{2}$ value fall off the end. More...

class	Fit_Results
	Holds set(s) of results from more than one kinematic fits. More...

class	Fourvec_Constrainer
	Do a kinematic fit for a set of four-momenta, given a set of mass constraints. More...

class	Fourvec_Constrainer_Args
	Hold on to parameters for the Fourvec_Constrainer class. More...

class	Fourvec_Constraint_Calculator
	Concrete realization of the Constraint_Calculator class. Evaluate constraints at the point described by x and y (well-measured and poorly-measured variables, respectively). The results should be stored in F. Bx and By should be set to the gradients of F with respect to x and y, respectively. More...

class	Fourvec_Event
	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...

class	Gentop_Args
	Hold on to parameters for the toy event generator. More...

class	JetTranslatorBase
	Template class of function object to translate jet physics object to HitFit's Lepjets_Event_Jet object. Users need to write an implementation of a template specialization of this class for their jet physics object class. Then users combine this header file and their implementation for their analysis code. With this approach, it is possible to use HitFit for different jet physics object class in different experiments. More...

class	Lepjets_Event
	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...

class	Lepjets_Event_Jet
	A class to represent a jet in an instance of Lepjets_Event class. The class is derived from the Lepjets_Event_Lep class. In addition to the information stored in Lepjets_Event_Lep class, this class holds the following information: More...

class	Lepjets_Event_Lep
	Represent a lepton in an instance of Lepjets_Event class. This class hold the following information: More...

class	LeptonTranslatorBase
	Template class of function object to translate lepton physics object to HitFit's Lepjets_Event_Lep object. Users need to write an implementation of a template specialization of this class for their lepton physics object class. Then users combine this header file and their implementation for their analysis code. With this approach, it is possible to use HitFit for different lepton physics object class in different experiments. More...

class	METTranslatorBase
	Template class of function object to translate missing transverse energy physics object to HitFit's Fourvec object. Users need to write an implementation of a template specialization of this class for their missing transverse energy physics object class, Then users combine this header file and their implementation for their analysis code. With this approach, it is possible to use HitFit for different missing transverse energy physics object class indifferent experiments. More...

class	Objpair
	Represent a pair of objects in Pair_Table. More...

class	Pair_Table
	A lookup table to speed up constraint evaluation using Fourvec_Constrainer. More...

class	Resolution

class	Row_Vector
	Row-vector class. CLHEP doesn't have a row-vector class, so HitFit uses its own. This is only a simple wrapper around Matrix that tries to constrain the shape to a row vector. It will raise an assertion if you try to assign to it something that isn't a row vector. More...

class	RunHitFit
	Template class of experiment-independent interface to HitFit. This class is intended to be used inside the programming environment of a specific experiment, where each type of physics objects has its own class/type. For using HitFit with generic four-vector classes, user can't use this class and have to use the Top_Fit class directly. The reason is: this class is designed assuming electron and muon are represented by different object type, a situation which is guaranteed to happen in any experiments. The class contains some static integer constants to limit the maximum amount of jets in an event before fitting. See the description of those constants for details. The numbers of permutations for $t\bar{t}$ and $t\bar{t}H$ as a function of the number of jets $N_{\mathrm{jet}}$ in the event for a few values of are. More...

class	Top_Decaykin
	A class to hold functions to calculate kinematic quantities of interest in $t\bar{t} \to \ell + 4 \mathrm{jets}$ events. This class has no state, only static member functions. More...

class	Top_Fit
	Handle and fit jet permutations of an event. This is the primary interface between user's Lepjets_Event and HitFit kinematic fitting algorithm. More...

class	Top_Fit_Args
	Hold on to parameters for the Top_Fit class. More...

class	Vector_Resolution

Typedefs
typedef CLHEP::HepVector	Column_Vector

typedef CLHEP::HepDiagMatrix	Diagonal_Matrix

typedef CLHEP::HepLorentzVector	Fourvec
	Typedef for a HepLorentzVector. More...

typedef CLHEP::HepMatrix	Matrix

typedef CLHEP::Hep3Vector	Threevec
	Typedef for a Hep3Vector. More...

Enumerations
enum	Jet_Labels { isr_label = 0, lepb_label = 11, hadb_label = 12, hadw1_label = 13, hadw2_label = 14, higgs_label = 15, unknown_label = 20 }

enum	Kt_Offsets { x_offs = 0, y_offs = 1 }

enum	Lepton_Labels { lepton_label = 1, electron_label = 2, muon_label = 3 }

enum	Lists { all_list = 0, noperm_list = 1, semicorrect_list = 2, limited_isr_list = 3, topfour_list = 4, btag_list = 5, htag_list = 6, n_lists = 7 }

enum	Offsets { p_offs = 0, phi_offs = 1, eta_offs = 2 }

enum	Unmeasured_Variables { nu_z = 1 }

Functions
void	adjust_e_for_mass (Fourvec &v, double mass)
	Adjust the energy component of four-vector v (leaving the three-vector part unchanged) so that the four-vector has mass as specified in the argument. More...

void	adjust_p_for_mass (Fourvec &v, double mass)
	Adjust the three-vector part of v, leaving the energy unchanged,. More...

void	clear (CLHEP::HepGenMatrix &m)
	Helper function: Reset all elements of a matrix to 0. More...

double	delta_r (const Fourvec &a, const Fourvec &b)
	Find the distance between two four-vectors in the two-dimensional space $\eta-\phi$ . More...

double	deteta (const Fourvec &v, double zvert)
	NOT USED ANYMORE: Get the detector $\eta$ (D0-specific), requires z-vertex. More...

double	eta_to_theta (double eta)
	Convert pseudorapidity to polar angle. More...

Lepjets_Event	gentop (const Gentop_Args &args, CLHEP::HepRandomEngine &engine)
	Generate a $t\bar{t} \to \ell + 4~\mathrm{jets}$ event. More...

Lepjets_Event	gentth (const Gentop_Args &args, CLHEP::HepRandomEngine &engine)
	Generate a $t\bar{t}H \to \ell + b\bar{b} + 4~\mathrm{jets}$ event. More...

char	jetTypeChar (int type)
	Helper function: Translate jet type code from integer to char. The following notation is used for each type of jet: More...

std::string	jetTypeString (int type)
	Helper function: Translate jet type code from integer to char. The following notation is used for each type of jet: More...

template<class T >
std::string	jetTypeString (std::vector< T > jet_types)
	Helper function: Translate jet type code from a list of numbers to a string. More...

std::unique_ptr< Constraint_Intermed >	make_constraint_intermed (std::string s)

unique_ptr< Constraint_Intermed >	make_constraint_intermed (string s)

bool	operator< (const EtaDepResElement &a, const EtaDepResElement &b)
	Comparison operator, compare two EtaDepResElement instances based on their respective valid $\eta$ ranges. More...

bool	operator< (const Fit_Result &a, const Fit_Result &b)
	Sort fit results based on their $\chi^{2}$ . More...

std::ostream &	operator<< (std::ostream &s, const Base_Constrainer &f)
	Output stream operator, print the content of this Base_Constrainer to an output stream. More...

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

std::ostream &	operator<< (std::ostream &s, const Constrained_Z &cz)
	Output stream operator, print the content of this Constrained_Z to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Constraint &c)
	Output stream operator, print the content of this Constraint to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Defaults_Text &def)

std::ostream &	operator<< (std::ostream &s, const EtaDepResElement &e)

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

std::ostream &	operator<< (std::ostream &s, const Fit_Result &res)
	Output stream operator, print the content of this Fit_Result to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Fit_Result_Vec &resvec)
	Output stream operator, print the content of this Fit_Result_Vec to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Fit_Results &res)
	Output stream operator, print the content of this Fit_Results to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Fourvec_Constrainer &c)
	Output stream operator, print the content of this Fourvec_Constrainer to an output stream. More...

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

std::ostream &	operator<< (std::ostream &s, const hitfit::Constraint_Intermed &ci)
	Output stream operator, print the content of this Constraint_Intermed to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Lepjets_Event &ev)
	Output stream operator, print the content of this Lepjets_Event to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Lepjets_Event_Jet &l)
	Output stream operator, print the content of this Lepjets_Event_Jet to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Lepjets_Event_Lep &l)
	Output stream operator, print the content of this Lepjets_Event_Lep to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Objpair &o)
	Output stream operator, print the content of this Objpair to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Pair_Table &p)
	Output stream operator, print the content of this Pair_Table to an output stream. More...

std::ostream &	operator<< (std::ostream &s, const Resolution &r)
	Output stream operator, print the content of this Resolution to an output stream. More...

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

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

double	phidiff (double phi)
	Normalized difference in azimuthal angles to a range between $[-\pi \dot \pi]$ . More...

void	roteta (Fourvec &v, double eta)
	Rotate four-vector v through a polar angle such that the four-vector pseudorapidity changes by a desired value. More...

void	rottheta (Fourvec &v, double theta)
	Rotate four-vector v through a polar angle. More...

double	scalar (const CLHEP::HepGenMatrix &m)
	Return the $1 \times 1$ matrix as a scalar. Raise an assertion if the matris is not $1 \times 1$ . More...

double	theta_to_eta (double theta)
	Convert polar angle to pseudorapidity. More...

Variables
const int	nu_label = -1

Typedef Documentation

◆ Column_Vector

typedef CLHEP::HepVector hitfit::Column_Vector

Definition at line 63 of file matutil.h.

◆ Diagonal_Matrix

typedef CLHEP::HepDiagMatrix hitfit::Diagonal_Matrix

Definition at line 64 of file matutil.h.

◆ Fourvec

typedef CLHEP::HepLorentzVector hitfit::Fourvec

Typedef for a HepLorentzVector.

Definition at line 55 of file fourvec.h.

◆ Matrix

typedef CLHEP::HepMatrix hitfit::Matrix

Definition at line 62 of file matutil.h.

◆ Threevec

typedef CLHEP::Hep3Vector hitfit::Threevec

Typedef for a Hep3Vector.

Definition at line 60 of file fourvec.h.

Enumeration Type Documentation

◆ Jet_Labels

enum hitfit::Jet_Labels

Possible types of jet in an instance of Lepjets_Event class.

Enumerator
isr_label
lepb_label
hadb_label
hadw1_label
hadw2_label
higgs_label
unknown_label

Definition at line 54 of file Lepjets_Event_Jet.h.

                   {
     isr_label = 0,
     lepb_label = 11,
     hadb_label = 12,
     hadw1_label = 13,
     hadw2_label = 14,
     higgs_label = 15,
     unknown_label = 20
   };

◆ Kt_Offsets

enum hitfit::Kt_Offsets

Offset indices for the components of missing transverse energy (or $k_{T}$ ) variables.

Enumerator
x_offs
y_offs

Definition at line 133 of file Fourvec_Constrainer.cc.

133 { x_offs = 0, y_offs = 1 } Kt_Offsets;

◆ Lepton_Labels

enum hitfit::Lepton_Labels

Possible types of lepton in an instance of Lepjets_Event class.

Enumerator
lepton_label
electron_label
muon_label

Definition at line 51 of file Lepjets_Event_Lep.h.

                      {
     lepton_label = 1,  // generic lepton
     electron_label = 2,
     muon_label = 3
   };

◆ Lists

enum hitfit::Lists

Indices for the various results lists we store in Fit_Results.

Enumerator
all_list
noperm_list
semicorrect_list
limited_isr_list
topfour_list
btag_list
htag_list
n_lists

Definition at line 53 of file Top_Fit.h.

              {
     all_list = 0,          // All events.
     noperm_list = 1,       // All jet assignments are correct.
     semicorrect_list = 2,  // Jets assigned to the correct top.
     limited_isr_list = 3,  // Top three jets are not ISR.
     topfour_list = 4,      // Top four jets are not ISR, any other jets are ISR.
     btag_list = 5,         // All tagged jets were assigned as b-jets.
     htag_list = 6,         // All tagged jets were assigned as b-jets or higgs.
     n_lists = 7
   };

◆ Offsets

enum hitfit::Offsets

Offset indices for the component of four-momentum variables.

Enumerator
p_offs
phi_offs
eta_offs

Definition at line 127 of file Fourvec_Constrainer.cc.

127 { p_offs = 0, phi_offs = 1, eta_offs = 2 } Offsets;

◆ Unmeasured_Variables

enum hitfit::Unmeasured_Variables

If there is a neutrino, then it is at index 1 of the vector of poorly-measured variables.

Enumerator
nu_z

Definition at line 143 of file Fourvec_Constrainer.cc.

143 { nu_z = 1 } Unmeasured_Variables;

Function Documentation

◆ adjust_e_for_mass()

void hitfit::adjust_e_for_mass	(	Fourvec &	v,
		double	mass
	)

Adjust the energy component of four-vector v (leaving the three-vector part unchanged) so that the four-vector has mass as specified in the argument.

Parameters

v	The four-vector to scale.
mass	The desired-mass of the four-vector.

Output:: v The scaled four-vector.

Definition at line 119 of file fourvec.cc.

   {
     v.setE(sqrt(v.vect().mag2() + mass * mass));
   }

References EgHLTOffHistBins_cfi::mass, mathSSE::sqrt(), and findQualityFiles::v.

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

◆ adjust_p_for_mass()

void hitfit::adjust_p_for_mass	(	Fourvec &	v,
		double	mass
	)

Adjust the three-vector part of v, leaving the energy unchanged,.

so that the four-vector has mass as specified in the argument.

Parameters

v	The four-vector to scale.
mass	The desired mass of the four-vector.

Output:: v The scaled four-vector.

Definition at line 95 of file fourvec.cc.

   {
     CLHEP::Hep3Vector vect = v.vect();
     double old_p2 = vect.mag2();
     if (old_p2 == 0)
       return;
     double new_p2 = v.e() * v.e() - mass * mass;
     if (new_p2 < 0)
       new_p2 = 0;
     vect *= sqrt(new_p2 / old_p2);
     v.setVect(vect);
   }

References EgHLTOffHistBins_cfi::mass, mathSSE::sqrt(), and findQualityFiles::v.

◆ clear()

void hitfit::clear ( CLHEP::HepGenMatrix & m )

Helper function: Reset all elements of a matrix to 0.

Parameters

m	The matrix to reset.

Definition at line 151 of file matutil.cc.

   {
     int nrow = m.num_row();
     int ncol = m.num_col();
     for (int i = 1; i <= nrow; i++)
       for (int j = 1; j <= ncol; j++)
         m(i, j) = 0;
   }

References mps_fire::i, dqmiolumiharvest::j, and visualization-live-secondInstance_cfg::m.

◆ delta_r()

double hitfit::delta_r	(	const Fourvec &	a,
		const Fourvec &	b
	)

Find the distance between two four-vectors in the two-dimensional space $\eta-\phi$ .

Definition at line 238 of file fourvec.cc.

   {
     double deta = a.pseudoRapidity() - b.pseudoRapidity();
     double dphi = phidiff(a.phi() - b.phi());
     return sqrt(deta * deta + dphi * dphi);
   }

References a, b, phidiff(), and mathSSE::sqrt().

Referenced by HGCalCLUEAlgoT< TILE >::calculateDistanceToHigher(), HGCalCLUEAlgoT< TILE >::calculateLocalDensity(), TrackingMaterialPlotter::fill_gradient(), DD4hep_TrackingMaterialPlotter::fill_gradient(), DD4hep_ListGroups::fillGradient(), ListGroups::fillGradient(), HGCalCLUEAlgoT< TILE >::findAndAssignClusters(), and HGCalCLUEAlgoT< TILE >::makeClusters().

◆ deteta()

double hitfit::deteta	(	const Fourvec &	v,
		double	zvert
	)

NOT USED ANYMORE: Get the detector $\eta$ (D0-specific), requires z-vertex.

Parameters

v	The four-vector on which to operate.
zvert	z-vertex of the event.

Definition at line 205 of file fourvec.cc.

   {
     return theta_to_eta(cal_th(v.theta(), zvert));
   }

References theta_to_eta(), and findQualityFiles::v.

Referenced by HcalNumberingFromPS::unitID().

◆ eta_to_theta()

double hitfit::eta_to_theta ( double eta )

Convert pseudorapidity to polar angle.

Parameters

eta	The value of pseudorapidity to convert.

Definition at line 177 of file fourvec.cc.

   {
     return 2 * atan(exp(-eta));
   }

References PVValHelper::eta, and JetChargeProducer_cfi::exp.

Referenced by roteta().

◆ gentop()

Lepjets_Event hitfit::gentop	(	const Gentop_Args &	args,
		CLHEP::HepRandomEngine &	engine
	)

Generate a $t\bar{t} \to \ell + 4~\mathrm{jets}$ event.

Parameters

args	The parameter settings for this event.
engine	The underlying random number generator.

Definition at line 513 of file gentop.cc.

   {
     CLHEP::RandBreitWigner rbw(engine);
     CLHEP::RandGauss rgauss(engine);
  
     // Get the t decay momentum in the ttbar rest frame.
     Threevec p = rand_pt(args.t_pt_mean(), engine);
  
     // Make the t/tbar vectors.
     Fourvec lept = make_massive(p, args.mt(), args.sigma_mt(), engine);
     Fourvec hadt = make_massive(-p, args.mt(), args.sigma_mt(), engine);
  
     // Boost the rest frame.
     Fourvec boost = rand_boost(args, engine);
     lept.boost(boost.boostVector());
     hadt.boost(boost.boostVector());
  
     // Decay t -> b W, leptonic side.
     Fourvec lepb, lepw;
     double mlb = rgauss.fire(args.mb(), args.sigma_mb());
     double mlw = rbw.fire(args.mw(), args.sigma_mw());
     decay(lept, mlb, mlw, engine, lepb, lepw);
  
     // Decay t -> b W, hadronic side.
     Fourvec hadb, hadw;
     double mhb = rgauss.fire(args.mb(), args.sigma_mb());
     double mhw = rbw.fire(args.mw(), args.sigma_mw());
     decay(hadt, mhb, mhw, engine, hadb, hadw);
  
     // Decay W -> l nu.
     Fourvec lep, nu;
     decay(lepw, 0, 0, engine, lep, nu);
  
     // Decay W -> qqbar.
     Fourvec q1, q2;
     decay(hadw, 0, 0, engine, q1, q2);
  
     // Fill in the event.
     Lepjets_Event ev(0, ++next_evnum);
     Vector_Resolution lep_res(args.muon() ? args.muo_res_str() : args.ele_res_str());
     Vector_Resolution jet_res(args.jet_res_str());
     Resolution kt_res = (args.kt_res_str());
  
     ev.add_lep(Lepjets_Event_Lep(lep, args.muon() ? muon_label : electron_label, lep_res));
  
     ev.add_jet(Lepjets_Event_Jet(lepb, lepb_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(hadb, hadb_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(q1, hadw1_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(q2, hadw2_label, jet_res));
  
     ev.met() = nu;
     ev.kt_res() = kt_res;
  
     // Simulate SVX tagging.
     tagsim(args, ev, engine);
  
     // Smear the event, if requested.
     if (args.smear())
       ev.smear(engine, args.smear_dir());
  
     // Done!
     return ev;
   }

References writedatasetfile::args, PA_ZEESkim_cff::decay, electron_label, ev, hadb_label, hadw1_label, hadw2_label, lepb_label, muon_label, AlCaHLTBitMon_ParallelJobs::p, q1, and q2.

◆ gentth()

Lepjets_Event hitfit::gentth	(	const Gentop_Args &	args,
		CLHEP::HepRandomEngine &	engine
	)

Generate a $t\bar{t}H \to \ell + b\bar{b} + 4~\mathrm{jets}$ event.

Parameters

args	The parameter settings for this event.
engine	The underlying random number generator.

Definition at line 587 of file gentop.cc.

   {
     CLHEP::RandBreitWigner rbw(engine);
     CLHEP::RandGauss rgauss(engine);
  
     // Generate three-vectors for two tops.
     Threevec p_t1 = rand_pt(args.t_pt_mean(), engine);
     Threevec p_t2 = rand_pt(args.t_pt_mean(), engine);
  
     // Conserve momentum.
     Threevec p_h = -(p_t1 + p_t2);
  
     // Construct the 4-vectors.
     Fourvec lept = make_massive(p_t1, args.mt(), args.sigma_mt(), engine);
     Fourvec hadt = make_massive(p_t2, args.mt(), args.sigma_mt(), engine);
     Fourvec higgs = make_massive(p_h, args.mh(), args.sigma_mh(), engine);
  
     // Boost the rest frame.
     Fourvec boost = rand_boost(args, engine);
     lept.boost(boost.boostVector());
     hadt.boost(boost.boostVector());
     higgs.boost(boost.boostVector());
  
     // Decay t -> b W, leptonic side.
     Fourvec lepb, lepw;
     decay(lept, rgauss.fire(args.mb(), args.sigma_mb()), rbw.fire(args.mw(), args.sigma_mw()), engine, lepb, lepw);
  
     // Decay t -> b W, hadronic side.
     Fourvec hadb, hadw;
     decay(hadt, rgauss.fire(args.mb(), args.sigma_mb()), rbw.fire(args.mw(), args.sigma_mw()), engine, hadb, hadw);
  
     // Decay W -> l nu.
     Fourvec lep, nu;
     decay(lepw, 0, 0, engine, lep, nu);
  
     // Decay W -> qqbar.
     Fourvec q1, q2;
     decay(hadw, 0, 0, engine, q1, q2);
  
     // Decay H -> bbbar.
     Fourvec hb1, hb2;
     decay(higgs, rgauss.fire(args.mb(), args.sigma_mb()), rgauss.fire(args.mb(), args.sigma_mb()), engine, hb1, hb2);
  
     // Fill in the event.
     Lepjets_Event ev(0, ++next_evnum);
     Vector_Resolution lep_res(args.muon() ? args.muo_res_str() : args.ele_res_str());
     Vector_Resolution jet_res(args.jet_res_str());
     Resolution kt_res = (args.kt_res_str());
  
     ev.add_lep(Lepjets_Event_Lep(lep, args.muon() ? muon_label : electron_label, lep_res));
  
     ev.add_jet(Lepjets_Event_Jet(lepb, lepb_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(hadb, hadb_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(q1, hadw1_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(q2, hadw2_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(hb1, higgs_label, jet_res));
     ev.add_jet(Lepjets_Event_Jet(hb2, higgs_label, jet_res));
  
     ev.met() = nu;
     ev.kt_res() = kt_res;
  
     // Simulate SVX tagging.
     tagsim(args, ev, engine);
  
     // Smear the event, if requested.
     if (args.smear())
       ev.smear(engine, args.smear_dir());
  
     // Done!
     return ev;
   }

References writedatasetfile::args, PA_ZEESkim_cff::decay, electron_label, ev, hadb_label, hadw1_label, hadw2_label, higgs_label, lepb_label, muon_label, q1, and q2.

◆ jetTypeChar()

char hitfit::jetTypeChar ( int type )

Helper function: Translate jet type code from integer to char. 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.

Parameters

type	The jet type code

Definition at line 224 of file Lepjets_Event_Jet.cc.

   {
     switch (j) {
       case hitfit::isr_label:
         return 'g';
       case hitfit::lepb_label:
         return 'b';
       case hitfit::hadb_label:
         return 'B';
       case hitfit::hadw1_label:
         return 'W';
       case hitfit::hadw2_label:
         return 'W';
       case hitfit::higgs_label:
         return 'h';
       case hitfit::unknown_label:
         return '?';
       default:
         return '?';
     }
  
     return '?';
   }

References hadb_label, hadw1_label, hadw2_label, higgs_label, isr_label, dqmiolumiharvest::j, lepb_label, and unknown_label.

Referenced by jetTypeString().

◆ jetTypeString() [1/2]

std::string hitfit::jetTypeString ( int type )

Helper function: Translate jet type code from integer to char. 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.

Parameters

type	The jet type code

Definition at line 257 of file Lepjets_Event_Jet.cc.

   {
     switch (j) {
       case hitfit::isr_label:
         return std::string("g");
       case hitfit::lepb_label:
         return std::string("b");
       case hitfit::hadb_label:
         return std::string("B");
       case hitfit::hadw1_label:
         return std::string("W");
       case hitfit::hadw2_label:
         return std::string("W");
       case hitfit::higgs_label:
         return std::string("h");
       case hitfit::unknown_label:
         return std::string("?");
       default:
         return std::string("?");
     }
  
     return std::string("?");
   }

References hadb_label, hadw1_label, hadw2_label, higgs_label, isr_label, dqmiolumiharvest::j, lepb_label, AlCaHLTBitMon_QueryRunRegistry::string, and unknown_label.

Referenced by hitfit::Lepjets_Event::jet_permutation().

◆ jetTypeString() [2/2]

template<class T >

std::string hitfit::jetTypeString ( std::vector< T > jet_types )

Helper function: Translate jet type code from a list of numbers to a string.

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

Parameters

jet_types The jet type codes in vector form.

Definition at line 248 of file Lepjets_Event_Jet.h.

                                                   {
     std::ostringstream ret;
  
     for (size_t j = 0; j != jet_types.size(); ++j) {
       ret << jetTypeChar((int)(jet_types[j]));
     }
  
     return ret.str();
   }

References dqmiolumiharvest::j, jetTypeChar(), and runTheMatrix::ret.

◆ make_constraint_intermed() [1/2]

std::unique_ptr<Constraint_Intermed> hitfit::make_constraint_intermed ( std::string s )

Helper function to parse input string s and construct the appropriate Constraint_Intermed instance. Returns null if the input string cannot be interpreted as a mass constraint.

The string should be either a numeric constant like

80.4

or a list of integers in parenthesis like

(1 2 4)

Leading spaces are ignored, as is text in a leading <> construction.

Parameters

s	The string to parse which contains information about the mass constraint.

Return:: Valid pointer to the appropriate type of constraint object, i.e. a Constraint_Intermed_Constant or a Constraint_Intermed_Labels.
NULL pointer if the input string cannot be interpreted as a mass constraint.

Referenced by hitfit::Constraint::Constraint().

◆ make_constraint_intermed() [2/2]

unique_ptr<Constraint_Intermed> hitfit::make_constraint_intermed ( string s )

Definition at line 284 of file Constraint_Intermed.cc.

   {
     // Skip leading spaces, `=', '< ... >'.
     string::size_type i = 0;
     while (i < s.size() && s[i] == ' ')
       ++i;
     if (s[i] == '=')
       ++i;
     while (i < s.size() && s[i] == ' ')
       ++i;
     if (i < s.size() && s[i] == '<') {
       i = s.find('>', i);
       if (i == string::npos)
         return unique_ptr<Constraint_Intermed>();
       ++i;
     }
     while (i < s.size() && s[i] == ' ')
       ++i;
  
     // Fail if there's nothing left.
     if (i == s.size())
       return unique_ptr<Constraint_Intermed>();
  
     if (s[i] == '(') {
       // List of labels.
       // Make a Constraint_Intermed_Labels instance.
       vector<int> labels;
       ++i;
       while (i < s.size()) {
         while (i < s.size() && s[i] == ' ')
           ++i;
         if (i < s.size() && s[i] == ')')
           break;
         if (i < s.size())
           labels.push_back(atoi(s.c_str() + i));
         while (i < s.size() && s[i] != ' ' && s[i] != ')')
           ++i;
       }
       return unique_ptr<Constraint_Intermed>(new Constraint_Intermed_Labels(labels));
     } else {
       // Make a Constraint_Intermed_Constant instance.
       return unique_ptr<Constraint_Intermed>(new Constraint_Intermed_Constant(atof(s.c_str() + i)));
     }
   }

References mps_fire::i, SummaryClient_cfi::labels, and alignCSCRings::s.

◆ operator<() [1/2]

bool hitfit::operator<	(	const EtaDepResElement &	a,
		const EtaDepResElement &	b
	)

Comparison operator, compare two EtaDepResElement instances based on their respective valid $\eta$ ranges.

Parameters

a	The first instance of EtaDepResElement to be compared.
b	The second instance of EtaDepResElement to be compared.

Return:: TRUE if a's upper limit is less than b's lower limit.
FALSE all other cases.

Definition at line 79 of file EtaDepResElement.cc.

                                                                        {
     if (a.IsOverlap(b)) {
       return false;
     }
     return !(a._EtaMax > b._EtaMin);
   }

References a, and b.

◆ operator<() [2/2]

bool hitfit::operator<	(	const Fit_Result &	a,
		const Fit_Result &	b
	)

Sort fit results based on their $\chi^{2}$ .

Parameters

a	The first instance of Fit_Result to compare.
b	The second instance of Fit_Result to compare.

Return:: TRUE if the first instance has smaller absolute value of $\chi^{2}$ than b.
FALSE all other cases.

Definition at line 174 of file Fit_Result.cc.

   {
     return abs(a._chisq) < abs(b._chisq);
   }

References a, funct::abs(), and b.

◆ operator<<() [1/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Base_Constrainer &	f
	)

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

Parameters

s	The output stream to which to write.
f	The instance of Base_Constrainer to be printed.

Definition at line 191 of file Base_Constrainer.cc.

   {
     return f.print(s);
   }

References f, and alignCSCRings::s.

◆ operator<<() [2/21]

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

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 298 of file Constrained_Top.cc.

   {
     return s << ct._constrainer;
   }

References hitfit::Constrained_Top::_constrainer, and alignCSCRings::s.

◆ operator<<() [3/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Constrained_Z &	cz
	)

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

Parameters

s	The output stream to which to wrire.
cz	The instance of Constrained_Z to be printed.

Definition at line 236 of file Constrained_Z.cc.

   {
     return s << cz._constrainer;
   }

References hitfit::Constrained_Z::_constrainer, and alignCSCRings::s.

◆ operator<<() [4/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Constraint &	c
	)

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

Parameters

s	The stream to which to write.
c	The instance of Constraint to be printed.

Definition at line 149 of file Constraint.cc.

   {
     s << *c._lhs.get() << " = " << *c._rhs.get();
     return s;
   }

References c, and alignCSCRings::s.

◆ operator<<() [5/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Defaults_Text &	def
	)

Output stream operator. Print out all parameters' names and their values.

Parameters

s	The output stream to write.
def	The instance to print.

Return:: The output stream s

Definition at line 430 of file Defaults_Text.cc.

   {
     for (std::map<std::string, std::string>::const_iterator it = def._rep->_map.begin(); it != def._rep->_map.end();
          it++) {
       s << "[" << it->first << "] = [" << it->second << "]\n";
     }
  
     return s;
   }

References spu::def(), and alignCSCRings::s.

◆ operator<<() [6/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const EtaDepResElement &	e
	)

Parameters

s	The output stream to write to.
e	The instance of EtaDepResElement to be printed.

Return:: The output stream s

Definition at line 110 of file EtaDepResElement.cc.

                                                                    {
     s << "(" << e._EtaMin << " to " << e._EtaMax << ")"
       << " / " << e.GetResolution();
     return s;
   }

References MillePedeFileConverter_cfg::e, and alignCSCRings::s.

◆ operator<<() [7/21]

std::ostream & hitfit::operator<<	(	std::ostream &	s,
		const FE_Obj &	o
	)

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

Parameters

s	The output stream to which to write.
o	The instance of FE_Obj to be printed.

Definition at line 80 of file Fourvec_Event.cc.

   {
     s << o.p << " - " << o.mass << " - " << o.label << "\n";
     s << "    errors: " << o.p_error << " " << o.phi_error << " " << o.eta_error;
     if (o.muon_p)
       s << " (mu)";
     s << "\n";
     return s;
   }

References EcalTangentSkim_cfg::o, and alignCSCRings::s.

◆ operator<<() [8/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Fit_Result &	res
	)

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

Parameters

s	The output stream to which to write.
res	The instance of Fit_Result to be printed.

Definition at line 199 of file Fit_Result.cc.

   {
     s << "chisq: " << res._chisq << "\n";
     s << "umwhad: " << res._umwhad << "\n";
     s << "utmass: " << res._utmass << "\n";
     s << "mt: " << res._mt << "\n";
     s << "sigmt: " << res._sigmt << "\n";
     s << res._ev;
     s << "pullx: " << res._pullx.T();
     s << "pully: " << res._pully.T();
     return s;
   }

References alignCSCRings::s.

◆ operator<<() [9/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Fit_Result_Vec &	resvec
	)

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

Parameters

s	The output stream to which to write.
resvec	The instance of Fit_Result_Vec to be printed.

Definition at line 159 of file Fit_Result_Vec.cc.

   {
     for (std::vector<std::shared_ptr<Fit_Result>>::size_type i = 0; i < resvec._v.size(); i++)
       s << "Entry " << i << "\n" << *resvec._v[i];
     return s;
   }

References hitfit::Fit_Result_Vec::_v, mps_fire::i, alignCSCRings::s, and trackerHitRTTI::vector.

◆ operator<<() [10/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Fit_Results &	res
	)

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

Parameters

s	The output stream to which to write.
res	The instance Fit_Results to be printed.

Definition at line 104 of file Fit_Results.cc.

   {
     for (std::vector<Fit_Result_Vec>::size_type i = 0; i < res._v.size(); i++)
       s << "List " << i << "\n" << res._v[i];
     return s;
   }

References mps_fire::i, and alignCSCRings::s.

◆ operator<<() [11/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Fourvec_Constrainer &	c
	)

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

Parameters

s	The output stream to which to write.
c	The instance of Fourvec_Constrainer to be printed.

Definition at line 220 of file Fourvec_Constrainer.cc.

   {
     s << "Constraints: (e_com = " << c._args.e_com() << ") ";
     if (c._args.use_e())
       s << "(E)";
     s << "\n";
  
     for (std::vector<Constraint>::size_type i = 0; i < c._constraints.size(); i++)
       s << "  " << c._constraints[i] << "\n";
  
     if (!c._mass_constraint.empty()) {
       s << "Mass constraint:\n";
       s << c._mass_constraint[0] << "\n";
     }
     return s;
   }

References c, mps_fire::i, if(), and alignCSCRings::s.

◆ operator<<() [12/21]

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

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

References hitfit::Fourvec_Event::_kt, hitfit::Fourvec_Event::_kt_x_error, hitfit::Fourvec_Event::_kt_xy_covar, hitfit::Fourvec_Event::_kt_y_error, hitfit::Fourvec_Event::_objs, hitfit::Fourvec_Event::_x, mps_fire::i, and alignCSCRings::s.

◆ operator<<() [13/21]

std::ostream & hitfit::operator<<	(	std::ostream &	s,
		const hitfit::Constraint_Intermed &	ci
	)

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

Parameters

s	The output stream to write.
ci	The instance of Constraint_Intermed to be printed.

Return:: The stream s.

Definition at line 268 of file Constraint_Intermed.cc.

   {
     ci.print(s);
     return s;
   }

References hitfit::Constraint_Intermed::print(), and alignCSCRings::s.

◆ operator<<() [14/21]

std::ostream & hitfit::operator<<	(	std::ostream &	s,
		const Lepjets_Event &	ev
	)

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

Parameters

s	The output stream to which to write.
ev	The instance of Lepjets_Event to be printed.

Definition at line 566 of file Lepjets_Event.cc.

   {
     return ev.dump(s);
   }

References ev, and alignCSCRings::s.

◆ operator<<() [15/21]

std::ostream & hitfit::operator<<	(	std::ostream &	s,
		const Lepjets_Event_Jet &	l
	)

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

Parameters

s	The stream to which to write.
l	The instance of Lepjets_Event_Jet to be printed.

Definition at line 209 of file Lepjets_Event_Jet.cc.

   {
     return l.dump(s);
   }

References cmsLHEtoEOSManager::l, and alignCSCRings::s.

◆ operator<<() [16/21]

std::ostream & hitfit::operator<<	(	std::ostream &	s,
		const Lepjets_Event_Lep &	l
	)

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

Parameters

s	The stream to which to write.
l	The instance of Lepjets_Event_Lep to be printed.

Definition at line 187 of file Lepjets_Event_Lep.cc.

   {
     return l.dump(s);
   }

References cmsLHEtoEOSManager::l, and alignCSCRings::s.

◆ operator<<() [17/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Objpair &	o
	)

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

Parameters

s	The stream to which to write.
o	The instance of Objpair to be printed.

Definition at line 74 of file Objpair.cc.

   {
     s << o._i << " " << o._j;
     for (unsigned k = 0; k < o._for_constraint.size(); ++k)
       s << " " << static_cast<int>(o._for_constraint[k]);
     return s;
   }

References dqmdumpme::k, EcalTangentSkim_cfg::o, and alignCSCRings::s.

◆ operator<<() [18/21]

std::ostream & hitfit::operator<<	(	std::ostream &	s,
		const Pair_Table &	p
	)

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

Parameters

s	The stream to which to write.
p	The instance of Pair_Table to be printed.

Definition at line 115 of file Pair_Table.cc.

   {
     for (int i = 0; i < p.npairs(); i++)
       s << " " << p.get_pair(i) << "\n";
     return s;
   }

References mps_fire::i, AlCaHLTBitMon_ParallelJobs::p, and alignCSCRings::s.

◆ operator<<() [19/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Resolution &	r
	)

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

Parameters

s	The stream to which to write.
r	The instance of Resolution to be printed.

Definition at line 211 of file Resolution.cc.

   {
     if (r._inverse)
       s << "-";
     s << r._constant_sigma << "," << r._resolution_sigma << "," << r._noise_sigma;
     return s;
   }

References alignCSCRings::r, and alignCSCRings::s.

◆ operator<<() [20/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Top_Fit &	fitter
	)

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

Parameters

s	The output stream to which to write.
fitter	The instance of Top_Fit to be printed.

Definition at line 548 of file Top_Fit.cc.

   {
     return s << fitter._constrainer;
   }

References hitfit::Top_Fit::_constrainer, and alignCSCRings::s.

◆ operator<<() [21/21]

std::ostream& hitfit::operator<<	(	std::ostream &	s,
		const Vector_Resolution &	r
	)

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

Parameters

s	The stream to which to write.
r	The instance of Vector_Resolution to be printed.

Definition at line 346 of file Vector_Resolution.cc.

   {
     s << r._p_res << "/ " << r._eta_res << "/ " << r._phi_res;
     if (r._use_et)
       s << "/et";
     s << "\n";
     return s;
   }

References alignCSCRings::r, and alignCSCRings::s.

◆ phidiff()

double hitfit::phidiff ( double phi )

Normalized difference in azimuthal angles to a range between $[-\pi \dot \pi]$ .

Parameters

phi	The azimuthal to be normalized.

Definition at line 220 of file fourvec.cc.

   {
     while (phi < -M_PI)
       phi += 2 * M_PI;
     while (phi > M_PI)
       phi -= 2 * M_PI;
     return phi;
   }

References M_PI.

Referenced by ResolutionCreator::analyze(), delta_r(), MuonMesh::isClusteredWith(), MuonMesh::isDuplicateOf(), ForwardDiskSectorBuilderFromWedges::operator()(), PtAssignment::process(), L1MuDTMuonSorter::runCOL(), and L1MuBMMuonSorter::runCOL().

◆ roteta()

void hitfit::roteta	(	Fourvec &	v,
		double	eta
	)

Rotate four-vector v through a polar angle such that the four-vector pseudorapidity changes by a desired value.

Parameters

v	The four-vector to rotate.
eta	The desired change in the pseudorapidity.

Output:: v The rotated four-vector.

Definition at line 156 of file fourvec.cc.

   {
     double theta1 = v.theta();
     double eta1 = theta_to_eta(theta1);
     double eta2 = eta1 + eta;
     double theta2 = eta_to_theta(eta2);
  
     rottheta(v, theta1 - theta2);
   }

References PVValHelper::eta, HLT_FULL_cff::eta1, HLT_FULL_cff::eta2, eta_to_theta(), rottheta(), theta_to_eta(), and findQualityFiles::v.

◆ rottheta()

void hitfit::rottheta	(	Fourvec &	v,
		double	theta
	)

Rotate four-vector v through a polar angle.

Parameters

v	The four-vector to rotate.
theta	The rotation angle.

Output:: v The rotated vector.

Definition at line 135 of file fourvec.cc.

   {
     double s = sin(theta), c = cos(theta);
     double old_pt = v.perp();
     double new_pt = old_pt * c - v.z() * s;
     v.setZ(old_pt * s + v.z() * c);
  
     v.setX(v.x() * new_pt / old_pt);
     v.setY(v.y() * new_pt / old_pt);
   }

References c, funct::cos(), alignCSCRings::s, funct::sin(), theta(), and findQualityFiles::v.

Referenced by roteta().

◆ scalar()

double hitfit::scalar ( const CLHEP::HepGenMatrix & m )

Return the $1 \times 1$ matrix as a scalar. Raise an assertion if the matris is not $1 \times 1$ .

Parameters

m	The matrix to convert, must be $1 \times 1$ .

Definition at line 166 of file matutil.cc.

   {
     assert(m.num_row() == 1 && m.num_col() == 1);
     return m(1, 1);
   }

References cms::cuda::assert(), and visualization-live-secondInstance_cfg::m.

Referenced by PhotonValidator::analyze(), PhotonAnalyzer::analyze(), hitfit::Chisq_Constrainer::fit(), conddblib::listObject(), Traj2TrackHits::split(), and brokenline::translateKarimaki().

◆ theta_to_eta()

double hitfit::theta_to_eta ( double theta )

Convert polar angle to pseudorapidity.

Parameters

theta The polar angle to convert.

Definition at line 191 of file fourvec.cc.

   {
     return -log(tan(theta / 2));
   }

References dqm-mbProfile::log, funct::tan(), and theta().

Referenced by deteta(), and roteta().

Variable Documentation

◆ nu_label

const int hitfit::nu_label = -1

A special label used for neutrino.

Definition at line 178 of file Fourvec_Event.h.

Referenced by hitfit::Fourvec_Event::add(), hitfit::Constrained_Top::Constrained_Top(), and hitfit::Fourvec_Event::set_nu_p().

Classes

Typedefs

Enumerations

Functions

Variables

Typedef Documentation

◆ Column_Vector

◆ Diagonal_Matrix

◆ Fourvec

◆ Matrix

◆ Threevec

Enumeration Type Documentation

◆ Jet_Labels

◆ Kt_Offsets

◆ Lepton_Labels

◆ Lists

◆ Offsets

◆ Unmeasured_Variables

Function Documentation

◆ adjust_e_for_mass()

◆ adjust_p_for_mass()

◆ clear()

◆ delta_r()

◆ deteta()

◆ eta_to_theta()

◆ gentop()

◆ gentth()

◆ jetTypeChar()

◆ jetTypeString() [1/2]

◆ jetTypeString() [2/2]

◆ make_constraint_intermed() [1/2]

◆ make_constraint_intermed() [2/2]

◆ operator<() [1/2]

◆ operator<() [2/2]

◆ operator<<() [1/21]

◆ operator<<() [2/21]

◆ operator<<() [3/21]

◆ operator<<() [4/21]

◆ operator<<() [5/21]

◆ operator<<() [6/21]

◆ operator<<() [7/21]

◆ operator<<() [8/21]

◆ operator<<() [9/21]

◆ operator<<() [10/21]

◆ operator<<() [11/21]

◆ operator<<() [12/21]

◆ operator<<() [13/21]

◆ operator<<() [14/21]

◆ operator<<() [15/21]

◆ operator<<() [16/21]

◆ operator<<() [17/21]

◆ operator<<() [18/21]

◆ operator<<() [19/21]

◆ operator<<() [20/21]

◆ operator<<() [21/21]

◆ phidiff()

◆ roteta()

◆ rottheta()

◆ scalar()

◆ theta_to_eta()

Variable Documentation

◆ nu_label