---

AlVertexFitter.cpp

Go to the documentation of this file.

#include "AlVertexFitter.h"

//
// Base class for vertex fitters
//
// 11/01  K. Huettmann
//
   
bool AlVertexFitter::doFit(AlephCollection<QvecBase*>& track_collection,
                           AlUserTrack& motherTrack, AlUsrVertex& endVertex) const {
 
  vector<Helix> vTracks;
  Helix tempT;

  // get geometrical track data
  AlephCollection<QvecBase*>::iterator Iqv;
  for (Iqv = track_collection.begin(); Iqv < track_collection.end(); Iqv++) {
    tempT.setTrackParameters(*Iqv);
    vTracks.push_back(tempT);
  }    

  // call the vertex fit routine using a specific algorithm 
  HepLorentzVector track4V; 
  Hep3Vector vertex3V;
  HepSymMatrix covTrack;
  HepSymMatrix covVertex;
  HepMatrix corrVertexTrack;
  float chi2; 
  bool failed = _algo->getResult(vTracks, track4V, covTrack,  
                                 vertex3V, covVertex, corrVertexTrack, chi2);
  
  if (! failed ) { 
    if (chi2 < 0.) {
      failed = ! failed;
    } else {
    // set vertex information
    endVertex.setA4V(vertex3V.x(), vertex3V.y(), vertex3V.z(), 0.);
    endVertex.setEM(covVertex);  
    endVertex.setCHIF((double) chi2);  
    
    // set mother track information and vertex-track relationships
    motherTrack.setA4V(track4V);                
    motherTrack.setEndVertex(& endVertex);
    motherTrack.setEM_kin(covTrack); 
    motherTrack.setCorrMatrix(corrVertexTrack);

    for (Iqv = track_collection.begin(); Iqv < track_collection.end(); Iqv++) 
      motherTrack.addDaughter(*Iqv);
    }
  }

  return failed;
  
}

---