KinematicConstrainedVertexUpdatorT< nTrk, nConstraint > Class Template Reference

#include <KinematicConstrainedVertexUpdatorT.h>

Public Member Functions
	KinematicConstrainedVertexUpdatorT ()
RefCountedKinematicVertex	update (const ROOT::Math::SVector< double, 3+7 *nTrk > &inState, ROOT::Math::SMatrix< double, 3+7 *nTrk, 3+7 *nTrk, ROOT::Math::MatRepSym< double, 3+7 *nTrk > > &inCov, std::vector< KinematicState > &lStates, const GlobalPoint &lPoint, GlobalVector const &fieldValue, MultiTrackKinematicConstraintT< nTrk, nConstraint > *cs)
	~KinematicConstrainedVertexUpdatorT ()

Private Attributes
ROOT::Math::SVector< double, 3+7 *nTrk >	delta_alpha
ROOT::Math::SVector< double, 3+7 *nTrk >	finPar
ROOT::Math::SMatrix< double, nConstraint+4, 3+7 *nTrk >	g
ROOT::Math::SVector< double, nConstraint+4 >	lambda
ROOT::Math::SMatrix< double, 7, 7, ROOT::Math::MatRepSym< double, 7 > >	nCovariance = ROOT::Math::SMatrixNoInit()
ROOT::Math::SMatrix< double, 3, 3, ROOT::Math::MatRepSym< double, 3 > >	pCov = ROOT::Math::SMatrixNoInit()
ROOT::Math::SMatrix< double, nConstraint+4, nConstraint+4, ROOT::Math::MatRepSym< double, nConstraint+4 > >	v_g_sym = ROOT::Math::SMatrixNoInit()
ROOT::Math::SVector< double, nConstraint+4 >	val
VertexKinematicConstraintT	vConstraint
KinematicVertexFactory	vFactory

Detailed Description

template<int nTrk, int nConstraint>
class KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >

Class caching the math part for KinematicConstrainedVertexFitter

Definition at line 32 of file KinematicConstrainedVertexUpdatorT.h.

Constructor & Destructor Documentation

template<int nTrk, int nConstraint>

KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::KinematicConstrainedVertexUpdatorT

(

)

Default constructor and destructor

Definition at line 74 of file KinematicConstrainedVertexUpdatorT.h.

{}

template<int nTrk, int nConstraint>

KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::~KinematicConstrainedVertexUpdatorT

(

)

Definition at line 78 of file KinematicConstrainedVertexUpdatorT.h.

{}

Member Function Documentation

template<int nTrk, int nConstraint>

RefCountedKinematicVertex KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::update	(	const ROOT::Math::SVector< double, 3+7 *nTrk > &	inState,
		ROOT::Math::SMatrix< double, 3+7 *nTrk, 3+7 *nTrk, ROOT::Math::MatRepSym< double, 3+7 *nTrk > > &	inCov,
		std::vector< KinematicState > &	lStates,
		const GlobalPoint &	lPoint,
		GlobalVector const &	fieldValue,
		MultiTrackKinematicConstraintT< nTrk, nConstraint > *	cs
	)

Method updating the states. Takes a vector of full parameters: (x,y,z,particle_1,...,particle_n), corresponding linearization point: vector of states and GlobalPoint, and constraint to be applied during the vertex fit. Returns refitted vector of 7n+3 parameters and corresponding covariance matrix, where n - number of tracks.

Definition at line 83 of file KinematicConstrainedVertexUpdatorT.h.

References KineDebug3::count(), g, mps_fire::i, MultiTrackKinematicConstraintBaseT::init(), invertPosDefMatrix(), LogDebug, MultiTrackKinematicConstraintT< NTRK, DIM >::parametersDerivative(), MultiTrackKinematicConstraintT< NTRK, DIM >::positionDerivative(), parseEventContent::prod, heppy_batch::val, MultiTrackKinematicConstraintT< NTRK, DIM >::value(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
  KineDebug3::count();

  int vSize = lStates.size();

  assert( nConstraint==0 || cs!=0);
  assert(vSize == nConstraint);

  const MagneticField* field=lStates.front().magneticField();
  
  delta_alpha=inPar;
  delta_alpha(0)-=lPoint.x();
  delta_alpha(1)-=lPoint.y();
  delta_alpha(2)-=lPoint.z();
  int cst=3;
  for(std::vector<KinematicState>::const_iterator i = lStates.begin(); i != lStates.end(); i++) 
    for ( int j=0; j<7; j++) {
      delta_alpha(cst)-=i->kinematicParameters()(j);
      cst++;
    }
  
  // cout<<"delta_alpha"<<delta_alpha<<endl;
  //resulting matrix of derivatives and vector of values.
  //their size  depends of number of tracks to analyze and number of
  //additional constraints to apply 
  
  if( nConstraint !=0) {
    cs->init(lStates, lPoint, fieldValue);
    val.Place_at(cs->value(),0);
    g.Place_at(cs->positionDerivative(),0,0);
    g.Place_at(cs->parametersDerivative(),0,3);
  }

  vConstraint.init(lStates, lPoint, fieldValue);
  val.Place_at(vConstraint.value(),nConstraint);
  g.Place_at(vConstraint.positionDerivative(),nConstraint, 0);
  g.Place_at(vConstraint.parametersDerivative(),nConstraint, 3);
  
  

  
  //debug code   
  v_g_sym = ROOT::Math::Similarity(g,inCov);
  
  // bool ifl1 = v_g_sym.Invert();
  bool ifl1 = invertPosDefMatrix(v_g_sym);
  if(!ifl1) {
    edm::LogWarning("KinematicConstrainedVertexUpdatorFailed")<< "invert failed\n"
                                  << v_g_sym;
    LogDebug("KinematicConstrainedVertexFitter3")
      << "Fit failed: unable to invert SYM gain matrix\n";
    return  RefCountedKinematicVertex();    
  }
  
  // delta alpha is now valid!
  //full math case now!
  val += g*delta_alpha;
  lambda = v_g_sym *val;
  
  //final parameters  
  finPar = inPar -  inCov * (ROOT::Math::Transpose(g) * lambda);
  
  //refitted covariance 
  ROOT::Math::SMatrix<double,3+7*nTrk,3+7*nTrk,ROOT::Math::MatRepSym<double,3+7*nTrk> > prod = ROOT::Math::SimilarityT(g,v_g_sym);
  ROOT::Math::SMatrix<double,3+7*nTrk,3+7*nTrk,ROOT::Math::MatRepSym<double,3+7*nTrk> > prod1;
  ROOT::Math::AssignSym::Evaluate(prod1, inCov * prod * inCov); 
  // ROOT::Math::AssignSym::Evaluate(prod, prod1 * inCov); 
  inCov -= prod1;
  
  pCov = inCov.template Sub< ROOT::Math::SMatrix<double,3,3,ROOT::Math::MatRepSym<double,3> > >(0,0);
  
  // chi2
  double chi  = ROOT::Math::Dot(lambda,val); //??
  
  //this is ndf without significant prior
  //vertex so -3 factor exists here 
  float ndf = 2*vSize - 3;
  ndf += nConstraint;
 
  
  //making resulting vertex 
  GlobalPoint vPos (finPar(0),finPar(1),finPar(2)); 
  VertexState st(vPos,GlobalError(pCov));
  RefCountedKinematicVertex rVtx = vFactory.vertex(st,chi,ndf);
  
  //making refitted states of Kinematic Particles
  AlgebraicVector7 newPar; 
  int i_int = 0;
  for(std::vector<KinematicState>::iterator i_st=lStates.begin(); i_st != lStates.end(); i_st++)
    {
      for(int i =0; i<7; i++)
    {newPar(i) = finPar(3 + i_int*7 + i);}
      
   nCovariance = inCov.template Sub<ROOT::Math::SMatrix<double, 7,7,ROOT::Math::MatRepSym<double,7> > >(3 + i_int*7, 3 + i_int*7);
   TrackCharge chl = i_st->particleCharge();
   KinematicParameters nrPar(newPar);
   KinematicParametersError nrEr(nCovariance);
   KinematicState newState(nrPar,nrEr,chl, field);
   (*i_st) = newState;
   i_int++;
  }
  return rVtx;  
}

Member Data Documentation

template<int nTrk, int nConstraint>

ROOT::Math::SVector<double,3+7*nTrk> KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::delta_alpha

private

Definition at line 60 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

ROOT::Math::SVector<double, 3+7*nTrk> KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::finPar

private

Definition at line 63 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

ROOT::Math::SMatrix<double,nConstraint+4,3+7*nTrk> KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::g

private

Definition at line 61 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

ROOT::Math::SVector<double, nConstraint+4> KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::lambda

private

Definition at line 64 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

ROOT::Math::SMatrix<double,7,7,ROOT::Math::MatRepSym<double,7> > KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::nCovariance = ROOT::Math::SMatrixNoInit()

private

Definition at line 66 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

ROOT::Math::SMatrix<double,3,3,ROOT::Math::MatRepSym<double,3> > KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::pCov = ROOT::Math::SMatrixNoInit()

private

Definition at line 65 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

ROOT::Math::SMatrix<double,nConstraint+4,nConstraint+4,ROOT::Math::MatRepSym<double,nConstraint+4> > KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::v_g_sym = ROOT::Math::SMatrixNoInit()

private

Definition at line 67 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

ROOT::Math::SVector<double,nConstraint+4> KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::val

private

Definition at line 62 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

VertexKinematicConstraintT KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::vConstraint

private

Definition at line 59 of file KinematicConstrainedVertexUpdatorT.h.

template<int nTrk, int nConstraint>

KinematicVertexFactory KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::vFactory

private

Definition at line 58 of file KinematicConstrainedVertexUpdatorT.h.