Constraint.h

Go to the documentation of this file.

//
//
// 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::unique_ptr<Constraint_Intermed> _lhs;

    std::unique_ptr<Constraint_Intermed> _rhs;
  };

}  // namespace hitfit

#endif  // not HITFIT_CONSTRAINT_H