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
double	fastsymmatrixaccess (double *m, int row, int col)
double	symmatrixaccess (double *m, int row, int col)
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

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

inline

Definition at line 67 of file LineTool.h.

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

Referenced by DecayTreeFitter::InternalParticle::initPar1().

  {
    // 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 ;
  }

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

inline

Definition at line 37 of file SortTool.h.

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

Referenced by DecayTreeFitter::InternalParticle::initPar1().

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

double DecayTreeFitter::fastsymmatrixaccess ( double *  m, int  row, int  col  )

inline

Definition at line 20 of file KalmanCalculator.cxx.

Referenced by symmatrixaccess().

  {
    return *(m+(row*(row-1))/2+(col-1));
  }

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

inline

Definition at line 58 of file LineTool.h.

  {
    return rhs.fillStream(os);
  }

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

inline

Definition at line 174 of file State.h.

References DecayTreeFitter::State::fillStream().

  {
    return obj.fillStream(str);
  }

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

inline

Definition at line 24 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 ;
  }

double DecayTreeFitter::symmatrixaccess ( double *  m, int  row, int  col  )

inline

Definition at line 25 of file KalmanCalculator.cxx.

References fastsymmatrixaccess().

  {
    return (row>=col? fastsymmatrixaccess(m,row,col) : fastsymmatrixaccess(m,col,row)) ;
  }