JetSubStructureUtils Namespace Reference

Classes
class	Angularity
class	BoostedXbbTag
class	Charge
class	Dipolarity
class	EnergyCorrelator
class	EnergyCorrelatorGeneralized
class	FoxWolfram
struct	jet_distance
class	JetDistanceCompare
class	KtDeltaR
class	KtMassDrop
class	KtSplittingScale
class	Multiplicities
class	Nsubjettiness
class	PlanarFlow
class	Pull
class	Qjets
class	QjetsPlugin
class	Qw
class	SoftDropObservables
class	SphericityTensor
class	SubjetFinder
class	SubstructureCalculator
class	Thrust
class	Volatility
class	ZCut

Functions
std::vector< fastjet::PseudoJet >	boostToCenterOfMass (const fastjet::PseudoJet &jet, std::vector< fastjet::PseudoJet > constituents)

Function Documentation

boostToCenterOfMass()

vector< fastjet::PseudoJet > JetSubStructureUtils::boostToCenterOfMass	(	const fastjet::PseudoJet &	jet,
		std::vector< fastjet::PseudoJet >	constituents )

Definition at line 9 of file BoostToCenterOfMass.cxx.

                                                                                                {
  vector<fastjet::PseudoJet> clusters;
  if(jet.e() < 1e-20) { // FPE
    return clusters;
  }
 
  double bx = jet.px()/jet.e();
  double by = jet.py()/jet.e();
  double bz = jet.pz()/jet.e();
 
  if(bx*bx + by*by + bz*bz >= 1) { // Faster than light
    return clusters;
  }
 
  for(unsigned int i1=0; i1 < constit_pseudojets.size(); i1++) {
    TLorentzVector v;
    v.SetPxPyPzE(constit_pseudojets.at(i1).px(), constit_pseudojets.at(i1).py(),constit_pseudojets.at(i1).pz(),constit_pseudojets.at(i1).e());
    v.Boost(-bx,-by,-bz);
    fastjet::PseudoJet v2(v.Px(), v.Py(), v.Pz(), v.E());
    clusters.push_back(v2);
  }
 
  return clusters;
}