/data/refman/pasoursint/CMSSW_5_3_10_patch1/src/TopQuarkAnalysis/TopHitFit/interface/Constraint.h

Go to the documentation of this file.

//
// $Id: Constraint.h,v 1.1 2011/05/26 09:46:53 mseidel Exp $
//
// File: hitfit/private/Constraint.h
// Purpose: Represent a mass constraint equation.
// Created: Jul, 2000, sss, based on run 1 mass analysis code.
//
// This class represents a mass constraint equation.
// Mass constraints come in two varieties, either saying that the sum
// of a set of labels should equal a constant:
//
//     (1 + 2) = 80
//
// or that two such sums should equal each other:
//
//     (1 + 2) = (3 + 4)
//
// We represent such a constraint equation by two Constraint_Intermed
// instances, each of which represents one side of the equation.
//
// CMSSW File      : interface/Constraint.h
// Original Author : Scott Stuart Snyder <snyder@bnl.gov> for D0
// Imported to CMSSW by Haryo Sumowidagdo <Suharyo.Sumowidagdo@cern.ch>


#ifndef HITFIT_CONSTRAINT_H
#define HITFIT_CONSTRAINT_H


#include <memory>
#include <string>
#include "TopQuarkAnalysis/TopHitFit/interface/Constraint_Intermed.h"


namespace hitfit {


class Fourvec_Event;


class Constraint
//
// Purpose: Represent a mass constraint equation.
//
{
public:
  // Constructor, destructor.  S is the string to parse describing
  // the constraint.

  Constraint (std::string s);

  Constraint (const Constraint& c);


  ~Constraint () {}

  // Assignment.
  Constraint& operator= (const Constraint& c);

  // See if this guy references both labels ILABEL and JLABEL
  // on a single side of the constraint equation.

  int has_labels (int ilabel, int jlabel) const;

  // Evaluate the mass constraint, using the data in EV.
  // Return m(lhs)^2/2 - m(rhs)^2/2.
  double sum_mass_terms (const Fourvec_Event& ev) const;

  // Print this object.
  friend std::ostream& operator<< (std::ostream& s, const Constraint& c);


private:
  // The two sides of the constraint.

  std::auto_ptr<Constraint_Intermed> _lhs;

  std::auto_ptr<Constraint_Intermed> _rhs;
};


} // namespace hitfit


#endif // not HITFIT_CONSTRAINT_H