TtFullHadSignalSel.cc File Reference

#include "TopQuarkAnalysis/TopEventSelection/interface/TtFullHadSignalSel.h"
#include "PhysicsTools/CandUtils/interface/EventShapeVariables.h"
#include "PhysicsTools/CandUtils/interface/Thrust.h"
#include "DataFormats/Math/interface/deltaR.h"
#include "TLorentzVector.h"
#include "TVector3.h"

Go to the source code of this file.

Functions
std::vector< math::XYZVector >	makeVecForEventShape (std::vector< pat::Jet > jets, bool only6Jets=true, ROOT::Math::Boost boost=ROOT::Math::Boost(0., 0., 0.))

Function Documentation

std::vector<math::XYZVector> makeVecForEventShape	(	std::vector< pat::Jet >	jets,
		bool	only6Jets = true,
		ROOT::Math::Boost	boost = ROOT::Math::Boost(0., 0., 0.)
	)

Definition at line 10 of file TtFullHadSignalSel.cc.

References metsig::jet, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by TtFullHadSignalSel::TtFullHadSignalSel().

                                                                                                   {
  std::vector<math::XYZVector> p;
  bool doBoost = (boost == ROOT::Math::Boost(0., 0., 0.)) ? false : true;
  for (std::vector<pat::Jet>::const_iterator jet = jets.begin(); jet != jets.end(); ++jet) {
    math::XYZVector Vjet;
    if (doBoost) {
      ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > Ljet(jet->px(), jet->py(), jet->pz(), jet->energy());
      Vjet = math::XYZVector(boost(Ljet).Px(), boost(Ljet).Py(), boost(Ljet).Pz());
    } else
      Vjet = math::XYZVector(jet->px(), jet->py(), jet->pz());
    p.push_back(Vjet);
    if (only6Jets && jet - jets.begin() == 5)
      break;
  }
  return p;
}