#include <KinematicParticleVertexFitter.h>

Public Member Functions
RefCountedKinematicTree	fit (const std::vector< RefCountedKinematicParticle > &particles) const

	KinematicParticleVertexFitter ()

	KinematicParticleVertexFitter (const edm::ParameterSet &pSet)

	~KinematicParticleVertexFitter ()

Private Member Functions
edm::ParameterSet	defaultParameters () const

void	setup (const edm::ParameterSet &pSet)

Private Attributes
VertexFitter< 6 > *	fitter

LinearizationPointFinder *	pointFinder

VertexTrackFactory< 6 > *	vFactory

Detailed Description

Class creating a kinematic particle out of set of daughter particles. Daughter particles are supposed to come from a common vertex. Arbitrary VerexFitter can be used to fit the common vertex and refit the daughter particles with the knowledge of vertex. The Kinematic Vertex is also created and the resulting KinematicParticle points on it.

Kirill Prokofiev, December 2002

Definition at line 25 of file KinematicParticleVertexFitter.h.

Constructor & Destructor Documentation

◆ KinematicParticleVertexFitter() [1/2]

KinematicParticleVertexFitter::KinematicParticleVertexFitter ( )

Constructor with LMSLinearizationPointFinder used as default.

Definition at line 15 of file KinematicParticleVertexFitter.cc.

                                                              {
   edm::ParameterSet pSet = defaultParameters();
   setup(pSet);
 }

References defaultParameters(), and setup().

◆ KinematicParticleVertexFitter() [2/2]

KinematicParticleVertexFitter::KinematicParticleVertexFitter ( const edm::ParameterSet & pSet )

Definition at line 20 of file KinematicParticleVertexFitter.cc.

20 { setup(pSet); }

References setup().

◆ ~KinematicParticleVertexFitter()

KinematicParticleVertexFitter::~KinematicParticleVertexFitter ( )

Definition at line 34 of file KinematicParticleVertexFitter.cc.

                                                               {
   delete vFactory;
   delete pointFinder;
   delete fitter;
 }

References fitter, pointFinder, and vFactory.

Member Function Documentation

◆ defaultParameters()

edm::ParameterSet KinematicParticleVertexFitter::defaultParameters ( ) const

private

Definition at line 40 of file KinematicParticleVertexFitter.cc.

                                                                        {
   edm::ParameterSet pSet;
   pSet.addParameter<double>("maxDistance", 0.01);
   pSet.addParameter<int>("maxNbrOfIterations", 100);  //10
   return pSet;
 }

References edm::ParameterSet::addParameter().

Referenced by KinematicParticleVertexFitter().

◆ fit()

RefCountedKinematicTree KinematicParticleVertexFitter::fit ( const std::vector< RefCountedKinematicParticle > & particles ) const

Fit method taking set of particles, fitting them to the common vertex and creating tree out of them. Input particles can belong to kinmaticTrees. In such a case it should be TOP particle of corresponding tree.

Definition at line 47 of file KinematicParticleVertexFitter.cc.

                                                                    {
   typedef ReferenceCountingPointer<VertexTrack<6> > RefCountedVertexTrack;
   //sorting the input
   if (particles.size() < 2)
     throw VertexException("KinematicParticleVertexFitter::input states are less than 2");
   InputSort iSort;
   std::pair<std::vector<RefCountedKinematicParticle>, std::vector<FreeTrajectoryState> > input = iSort.sort(particles);
   std::vector<RefCountedKinematicParticle> &newPart = input.first;
   std::vector<FreeTrajectoryState> &freeStates = input.second;
  
   GlobalPoint linPoint = pointFinder->getLinearizationPoint(freeStates);
  
   // cout<<"Linearization point found"<<endl;
  
   //making initial veretx seed with lin point as position and a fake error
   AlgebraicSymMatrix33 we;
   we(0, 0) = we(1, 1) = we(2, 2) = 10000.;
   GlobalError error(we);
   VertexState state(linPoint, error);
  
   //vector of Vertex Tracks to fit
   std::vector<RefCountedVertexTrack> ttf;
   TrackKinematicStatePropagator propagator_;
   for (auto const &i : newPart) {
     if (!(i)->currentState().isValid() || !propagator_.willPropagateToTheTransversePCA((i)->currentState(), linPoint)) {
       // std::cout << "Here's the bad state." << std::endl;
       return ReferenceCountingPointer<KinematicTree>(new KinematicTree());  // return invalid vertex
     }
     ttf.push_back(vFactory->vertexTrack((i)->particleLinearizedTrackState(linPoint), state, 1.));
   }
  
   // //debugging code to check neutrals:
   //  for(std::vector<RefCountedVertexTrack>::const_iterator i = ttf.begin(); i!=ttf.end(); i++)
   //  {
   // //   cout<<"predicted state momentum error"<<(*i)->linearizedTrack()->predictedStateMomentumError()<<endl;
   // //  cout<<"Momentum jacobian"<<(*i)->linearizedTrack()->momentumJacobian() <<endl;
   //  //  cout<<"predicted state momentum "<<(*i)->linearizedTrack()->predictedStateMomentum()<<endl;
   // //   cout<<"constant term"<<(*i)->linearizedTrack()->constantTerm()<<endl;
   //
   //  }
  
   CachingVertex<6> vtx = fitter->vertex(ttf);
   if (!vtx.isValid()) {
     LogDebug("RecoVertex/KinematicParticleVertexFitter") << "Fitted position is invalid. Returned Tree is invalid\n";
     return ReferenceCountingPointer<KinematicTree>(new KinematicTree());  // return invalid vertex
   }
   FinalTreeBuilder tBuilder;
   return tBuilder.buildTree(vtx, newPart);
 }

References FinalTreeBuilder::buildTree(), relativeConstraints::error, fitter, LinearizationPointFinder::getLinearizationPoint(), mps_fire::i, input, sistrip::SpyUtilities::isValid(), LogDebug, ecalTrigSettings_cff::particles, pointFinder, InputSort::sort(), VertexFitter< N >::vertex(), VertexTrackFactory< N >::vertexTrack(), vFactory, extraflags_cff::vtx, and TrackKinematicStatePropagator::willPropagateToTheTransversePCA().

Referenced by KineExample::analyze(), trackingPlots.Iteration::modules(), and PFTau3ProngReco::produce().

◆ setup()

void KinematicParticleVertexFitter::setup ( const edm::ParameterSet & pSet )

private

Definition at line 22 of file KinematicParticleVertexFitter.cc.

                                                                      {
   pointFinder = new DefaultLinearizationPointFinder();
   vFactory = new VertexTrackFactory<6>();
  
   KalmanVertexTrackUpdator<6> vtu;
   KalmanSmoothedVertexChi2Estimator<6> vse;
   KalmanTrackToTrackCovCalculator<6> covCalc;
   SequentialVertexSmoother<6> smoother(vtu, vse, covCalc);
   fitter = new SequentialVertexFitter<6>(
       pSet, *pointFinder, KalmanVertexUpdator<6>(), smoother, ParticleKinematicLinearizedTrackStateFactory());
 }