JiveXML::TrackParticleRetrieverHelpers Namespace Reference

Functions
void	getPolylineFromTrackParticle (const Rec::TrackParticle *track, DataVect &polylineX, DataVect &polylineY, DataVect &polylineZ, DataVect &numPolyline)
	Get polyline hits if available.

Function Documentation

getPolylineFromTrackParticle()

void JiveXML::TrackParticleRetrieverHelpers::getPolylineFromTrackParticle	(	const Rec::TrackParticle *	track,
		DataVect &	polylineX,
		DataVect &	polylineY,
		DataVect &	polylineZ,
		DataVect &	numPolyline )

Get polyline hits if available.

Modified from TrackRetriever

Definition at line 46 of file TrackParticleRetriever.cxx.

                                                                                               {
 
     DataVect myPositionX;
     DataVect myPositionY;
     DataVect myPositionZ;
 
     int numPoly = 0 ;
 
// Migration Dec'13: Not yet done here ! Need equivalent code
// to access spacepoints. jpt 13Dec13
 
     unsigned int i = 0;
     const std::vector<const Trk::TrackParameters*>&  trackpars =track->trackParameters();
     
     if (!trackpars.empty() && trackpars.size() == 3 ) { //polylines with less than 3 entries not useful 
 
       bool needresorting = trackpars.at(0)!=track->perigee();//Needed since TrackParticles are (at the moment)
                                                              //created with the first parameter put last
       for ( i=0; i<trackpars.size(); i++){
     const Amg::Vector3D& pos = trackpars.at(i)->position();
             myPositionX.push_back(DataType(pos.x()/10.));
             myPositionY.push_back(DataType(pos.y()/10.));
             myPositionZ.push_back(DataType(pos.z()/10.));
       //std::cout << " Polyline # " << i << " x: " << pos.x() << std::endl;
       ++numPoly;
       } 
       if (needresorting){ // swapping vector around. Is this too naive ? 
     polylineX.push_back( myPositionX.at( trackpars.size()-1 ) );
         polylineY.push_back( myPositionY.at( trackpars.size()-1 ) );
     polylineZ.push_back( myPositionZ.at( trackpars.size()-1 ) );
         polylineX.push_back( myPositionX.at( trackpars.size()-3 ) );
         polylineY.push_back( myPositionY.at( trackpars.size()-3 ) );
         polylineZ.push_back( myPositionZ.at( trackpars.size()-3 ) );
         polylineX.push_back( myPositionX.at( trackpars.size()-2 ) );
         polylineY.push_back( myPositionY.at( trackpars.size()-2 ) );
         polylineZ.push_back( myPositionZ.at( trackpars.size()-2 ) );
       }
     }
 
     // std::cout << " numPolyline" << numPoly << std::endl;
 
     //Store counter as well
     numPolyline.push_back(DataType(numPoly)); 
  }