DecayTreeFitter Namespace Reference

Classes
class	ChiSquare
class	Configuration
class	Constraint
class	DecayChain
class	ErrCode
class	FitParams
class	Fitter
class	InteractionPoint
class	InternalParticle
class	KalmanCalculator
class	Line
class	MergedConstraint
class	MissingParticle
class	ParticleBase
class	Projection
class	RecoComposite
class	RecoParticle
class	RecoPhoton
class	RecoResonance
class	RecoTrack
class	RecoTrackStateProvider
class	Resonance
class	State

Functions
std::ostream &	operator<< (std::ostream &os, const Line &rhs)
bool	closestPointParams (const Line &line0, const Line &line1, double &mu0, double &mu1)
bool	sortByType (const ParticleBase *lhs, const ParticleBase *rhs)
bool	compTrkTransverseMomentum (const RecoTrack *lhs, const RecoTrack *rhs)
std::ostream &	operator<< (std::ostream &str, const State &obj)

Function Documentation

closestPointParams()

bool DecayTreeFitter::closestPointParams ( const Line &  line0, const Line &  line1, double &  mu0, double &  mu1  )

inline

Definition at line 58 of file LineTool.h.

References DecayTreeFitter::Line::beginPoint(), DecayTreeFitter::Line::direction(), and fabs().

{
  // lhs:

  bool OK = true;

  // the matrix:
  const double a00 = line0.direction().Mag2();
  const double a10 = line0.direction().Dot(line1.direction());
  const double a01 = -a10;
  const double a11 = -line1.direction().Mag2();

  // the inverse determinant:
  const double det = (a00 * a11 - a01 * a10); // det = -sin^2(angle(line0.dir,line1.dir))
  if (std::fabs(det) < 1e-10) {
    OK = false; // parallel
  } else {

    const double detinv = 1.0 / det;

    // rhs:
    const TVector3 p1_p0 = line1.beginPoint() - line0.beginPoint();

    const double b0 = p1_p0.Dot(line0.direction());
    const double b1 = p1_p0.Dot(line1.direction());

    // get the Kramer solutions:

    mu0 = (b0 * a11 - b1 * a01) * detinv;
    mu1 = (a00 * b1 - a10 * b0) * detinv;
  }

  return OK;
}

compTrkTransverseMomentum()

bool DecayTreeFitter::compTrkTransverseMomentum ( const RecoTrack *  lhs, const RecoTrack *  rhs  )

inline

Definition at line 36 of file SortTool.h.

References DecayTreeFitter::ParticleBase::particle(), and RhoCandidate::Pt().

{
  return lhs->particle()->Pt() > rhs->particle()->Pt();
}

operator<<() [1/2]

std::ostream& DecayTreeFitter::operator<< ( std::ostream &  os, const Line &  rhs  )

inline

Definition at line 51 of file LineTool.h.

References DecayTreeFitter::Line::fillStream().

{
  return rhs.fillStream(os);
}

operator<<() [2/2]

std::ostream& DecayTreeFitter::operator<< ( std::ostream &  str, const State &  obj  )

inline

Definition at line 177 of file State.h.

References DecayTreeFitter::State::fillStream().

{
  return obj.fillStream(str);
}

sortByType()

bool DecayTreeFitter::sortByType ( const ParticleBase *  lhs, const ParticleBase *  rhs  )

inline

Definition at line 23 of file SortTool.h.

References DecayTreeFitter::ParticleBase::kRecoTrack, RhoCandidate::NDaughters(), DecayTreeFitter::ParticleBase::nFinalChargedCandidates(), DecayTreeFitter::ParticleBase::particle(), RhoCandidate::Pt(), and DecayTreeFitter::ParticleBase::type().

{
  int lhstype = lhs->type();
  int rhstype = rhs->type();
  bool rc = false;
  if (lhstype == rhstype && lhstype == ParticleBase::kRecoTrack)
    rc = lhs->particle()->Pt() > rhs->particle()->Pt();
  else if (lhs->particle()->NDaughters() > 0 && rhs->particle()->NDaughters() > 0)
    rc = lhs->nFinalChargedCandidates() > rhs->nFinalChargedCandidates();
  else
    rc = lhstype < rhstype;
  return rc;
}