#include <GsfTrackProducerBase.h>

Inheritance diagram for GsfTrackProducerBase:

Public Member Functions
	GsfTrackProducerBase (bool trajectoryInEvent, bool split)
	Constructor. More...

virtual void	putInEvt (edm::Event &, const Propagator prop, const MeasurementTracker measTk, std::unique_ptr< TrackingRecHitCollection > &, std::unique_ptr< reco::GsfTrackCollection > &, std::unique_ptr< reco::TrackExtraCollection > &, std::unique_ptr< reco::GsfTrackExtraCollection > &, std::unique_ptr< std::vector< Trajectory > > &, AlgoProductCollection &, TransientTrackingRecHitBuilder const , const reco::BeamSpot &, const TrackerTopology ttopo)
	Put produced collections in the event. More...

Public Member Functions inherited from TrackProducerBase< reco::GsfTrack >
const edm::ParameterSet &	getConf () const

virtual void	getFromES (const edm::EventSetup &, edm::ESHandle< TrackerGeometry > &, edm::ESHandle< MagneticField > &, edm::ESHandle< TrajectoryFitter > &, edm::ESHandle< Propagator > &, edm::ESHandle< MeasurementTracker > &, edm::ESHandle< TransientTrackingRecHitBuilder > &)
	Get needed services from the Event Setup. More...

virtual void	getFromEvt (edm::Event &, edm::Handle< TrackCandidateCollection > &, reco::BeamSpot &)
	Get TrackCandidateCollection from the Event (needed by TrackProducer) More...

virtual void	getFromEvt (edm::Event &, edm::Handle< TrackCollection > &, reco::BeamSpot &)
	Get TrackCollection from the Event (needed by TrackRefitter) More...

virtual void	produce (edm::Event &, const edm::EventSetup &)=0
	Method where the procduction take place. To be implemented in concrete classes. More...

void	setAlias (std::string alias)
	set the aliases of produced collections More...

void	setClusterRemovalInfo (const edm::InputTag &clusterRemovalInfo)
	Sets the information on cluster removal, and turns it on. More...

void	setConf (const edm::ParameterSet &conf)
	Set parameter set. More...

void	setSecondHitPattern (Trajectory traj, reco::GsfTrack &track, const Propagator prop, const MeasurementTrackerEvent measTk, const TrackerTopology ttopo)

void	setSrc (const edm::EDGetToken &src, const edm::EDGetTokenT< reco::BeamSpot > &bsSrc, const edm::EDGetTokenT< MeasurementTrackerEvent > &mteSrc)
	set label of source collection More...

	TrackProducerBase (bool trajectoryInEvent=false)
	Constructor. More...

virtual	~TrackProducerBase ()
	Destructor. More...

Protected Member Functions
void	fillMode (reco::GsfTrack &track, const TrajectoryStateOnSurface innertsos, const Propagator &gsfProp, const TransverseImpactPointExtrapolator &tipExtrapolator, TrajectoryStateClosestToBeamLineBuilder &tscblBuilder, const reco::BeamSpot &bs) const

void	fillStates (TrajectoryStateOnSurface tsos, std::vector< reco::GsfComponent5D > &states) const

Private Member Functions
bool	computeModeAtTM (const TrajectoryMeasurement &tm, reco::GsfTrackExtra::Point &position, reco::GsfTrackExtra::Vector &momentum, Measurement1D &deltaP) const
	position, momentum and estimated deltaP at an intermediate measurement (true if successful) More...

void	localParametersFromQpMode (const TrajectoryStateOnSurface tsos, AlgebraicVector5 &parameters, AlgebraicSymMatrix55 &covariance) const
	local parameters rescaled with q/p from mode More...

Private Attributes
bool	useSplitting

Additional Inherited Members
Public Types inherited from TrackProducerBase< reco::GsfTrack >
using	AlgoProductCollection = typename Base::AlgoProductCollection

using	Base = AlgoProductTraits< reco::GsfTrack >

using	TrackCollection = typename Base::TrackCollection

Public Types inherited from AlgoProductTraits< reco::GsfTrack >
using	AlgoProductCollection = std::vector< AlgoProduct >

using	TrackCollection = std::vector< reco::GsfTrack >

Protected Attributes inherited from TrackProducerBase< reco::GsfTrack >
std::string	alias_

edm::EDGetTokenT< reco::BeamSpot >	bsSrc_

edm::InputTag	clusterRemovalInfo_

edm::EDGetTokenT < MeasurementTrackerEvent >	mteSrc_

bool	rekeyClusterRefs_

edm::OrphanHandle < TrackCollection >	rTracks_

edm::ESHandle< NavigationSchool >	theSchool

bool	trajectoryInEvent_

Detailed Description

Produce Tracks from TrackCandidates

Author: cerati

Definition at line 26 of file GsfTrackProducerBase.h.

Constructor & Destructor Documentation

GsfTrackProducerBase::GsfTrackProducerBase	(	bool	trajectoryInEvent,
		bool	split
	)

inlineexplicit

Constructor.

Definition at line 30 of file GsfTrackProducerBase.h.

                                                                     :
   TrackProducerBase<reco::GsfTrack>(trajectoryInEvent),
     useSplitting(split){}

Member Function Documentation

bool GsfTrackProducerBase::computeModeAtTM	(	const TrajectoryMeasurement &	tm,
		reco::GsfTrackExtra::Point &	position,
		reco::GsfTrackExtra::Vector &	momentum,
		Measurement1D &	deltaP
	)		const

private

position, momentum and estimated deltaP at an intermediate measurement (true if successful)

Definition at line 379 of file GsfTrackProducerBase.cc.

References TrajectoryMeasurement::backwardPredictedState(), CollinearFitAtTM2::deltaP(), TrajectoryMeasurement::forwardPredictedState(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), localParametersFromQpMode(), LocalTrajectoryError::matrix(), TrajectoryMeasurement::recHit(), mps_fire::result, TrajectoryMeasurement::updatedState(), LocalTrajectoryParameters::vector(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by putInEvt().

 {  
   //
   // states
   //
   TrajectoryStateOnSurface fwdState = tm.forwardPredictedState();
   TrajectoryStateOnSurface bwdState = tm.backwardPredictedState();
   TrajectoryStateOnSurface upState  = tm.updatedState();
 
 
   if ( !fwdState.isValid() || !bwdState.isValid() || !upState.isValid() ) {
     return false;
   }
   //
   // position from mean, momentum from mode (in cartesian coordinates)
   //  following PF code
   //
   GlobalPoint pos = upState.globalPosition();
   position = reco::GsfTrackExtra::Point(pos.x(),pos.y(),pos.z());
   MultiTrajectoryStateMode mts;
   GlobalVector mom;
   bool result = mts.momentumFromModeCartesian(upState,mom);
   if ( !result ) {
 //     std::cout << "momentumFromModeCartesian failed" << std::endl;
     return false;
   }
   momentum = reco::GsfTrackExtra::Vector(mom.x(),mom.y(),mom.z());
   //
   // calculation from deltaP from fit to forward & backward predictions
   //  (momentum from mode) and hit
   //
   // prepare input parameter vectors and covariance matrices
   AlgebraicVector5 fwdPars = fwdState.localParameters().vector();
   AlgebraicSymMatrix55 fwdCov = fwdState.localError().matrix();
   localParametersFromQpMode(fwdState,fwdPars,fwdCov);
   AlgebraicVector5 bwdPars = bwdState.localParameters().vector();
   AlgebraicSymMatrix55 bwdCov = bwdState.localError().matrix();
   localParametersFromQpMode(bwdState,bwdPars,bwdCov);
   LocalPoint hitPos(0.,0.,0.);
   LocalError hitErr(-1.,-1.,-1.);
   if ( tm.recHit()->isValid() ) {
     hitPos = tm.recHit()->localPosition();
     hitErr = tm.recHit()->localPositionError();
   }    
   CollinearFitAtTM2 collinearFit(fwdPars,fwdCov,bwdPars,bwdCov,hitPos,hitErr);
   deltaP = collinearFit.deltaP();
 
   return true;
 }

void GsfTrackProducerBase::fillMode	(	reco::GsfTrack &	track,
		const TrajectoryStateOnSurface	innertsos,
		const Propagator &	gsfProp,
		const TransverseImpactPointExtrapolator &	tipExtrapolator,
		TrajectoryStateClosestToBeamLineBuilder &	tscblBuilder,
		const reco::BeamSpot &	bs
	)		const

protected

Definition at line 282 of file GsfTrackProducerBase.cc.

References FreeTrajectoryState::charge(), FreeTrajectoryState::curvilinearError(), reco::GsfTrack::dimensionMode, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), TransverseImpactPointExtrapolator::extrapolate(), CurvilinearTrajectoryError::matrix(), SingleGaussianState1D::mean(), GaussianSumUtilities1D::mean(), GaussianSumUtilities1D::mode(), GaussianSumUtilities1D::modeIsValid(), FreeTrajectoryState::momentum(), MultiGaussianStateTransform::multiState1D(), reco::BeamSpot::position(), Propagator::propagate(), reco::GsfTrack::setMode(), mathSSE::sqrt(), TrajectoryStateOnSurface::surface(), SingleGaussianState1D::variance(), GaussianSumUtilities1D::variance(), reco::TrackBase::vz(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by putInEvt().

 {
   // Get transverse impact parameter plane (from mean). This is a first approximation;
   // the mode is then extrapolated to the
   // final position closest to the beamline.
   GlobalPoint bsPos(bs.position().x()+(track.vz()-bs.position().z())*bs.dxdz(),
             bs.position().y()+(track.vz()-bs.position().z())*bs.dydz(),
             track.vz());
   TrajectoryStateOnSurface vtxTsos = tipExtrapolator.extrapolate(innertsos,bsPos);
   if ( !vtxTsos.isValid() )  vtxTsos = innertsos;
   // extrapolate mixture
   vtxTsos = gsfProp.propagate(innertsos,vtxTsos.surface());
   if ( !vtxTsos.isValid() )  return;              // failed (GsfTrack keeps mode = mean)
   // extract mode
   // build perigee parameters (for covariance to be stored)
   AlgebraicVector5 modeParameters;
   AlgebraicSymMatrix55 modeCovariance;
   // set parameters and variances for "mode" state (local parameters)
   for ( unsigned int iv=0; iv<5; ++iv ) {
     MultiGaussianState1D state1D = MultiGaussianStateTransform::multiState1D(vtxTsos,iv);
     GaussianSumUtilities1D utils(state1D);
     modeParameters(iv) = utils.mode().mean();
     modeCovariance(iv,iv) = utils.mode().variance();
     if ( !utils.modeIsValid() ) {
       // if mode calculation fails: use mean
       modeParameters(iv) = utils.mean();
       modeCovariance(iv,iv) = utils.variance();
     }
   }
   // complete covariance matrix
   // approximation: use correlations from mean
   const AlgebraicSymMatrix55& meanCovariance(vtxTsos.localError().matrix());
   for ( unsigned int iv1=0; iv1<5; ++iv1 ) {
     for ( unsigned int iv2=0; iv2<iv1; ++iv2 ) {
       double cov12 = meanCovariance(iv1,iv2) * 
     sqrt(modeCovariance(iv1,iv1)/meanCovariance(iv1,iv1)*
          modeCovariance(iv2,iv2)/meanCovariance(iv2,iv2));
       modeCovariance(iv1,iv2) = modeCovariance(iv2,iv1) = cov12;
     }
   }
   TrajectoryStateOnSurface modeTsos(LocalTrajectoryParameters(modeParameters,
                                   vtxTsos.localParameters().pzSign()),
                     LocalTrajectoryError(modeCovariance),
                     vtxTsos.surface(),
                     vtxTsos.magneticField(),
                     vtxTsos.surfaceSide());
   TrajectoryStateClosestToBeamLine tscbl = tscblBuilder(*modeTsos.freeState(),bs);
   if ( !tscbl.isValid() )  return;            // failed (GsfTrack keeps mode = mean)
   //
   // extract state at PCA and create momentum vector and covariance matrix
   //
   FreeTrajectoryState fts = tscbl.trackStateAtPCA();
   GlobalVector tscblMom = fts.momentum();
   reco::GsfTrack::Vector mom(tscblMom.x(),tscblMom.y(),tscblMom.z());
   reco::GsfTrack::CovarianceMatrixMode cov;
   const AlgebraicSymMatrix55& tscblCov = fts.curvilinearError().matrix();
   for ( unsigned int iv1=0; iv1<reco::GsfTrack::dimensionMode; ++iv1 ) {
     for ( unsigned int iv2=0; iv2<reco::GsfTrack::dimensionMode; ++iv2 ) {
       cov(iv1,iv2) = tscblCov(iv1,iv2);
     }
   } 
   track.setMode(fts.charge(),mom,cov);
 }

void GsfTrackProducerBase::fillStates	(	TrajectoryStateOnSurface	tsos,
		std::vector< reco::GsfComponent5D > &	states
	)		const

protected

Definition at line 269 of file GsfTrackProducerBase.cc.

References makeMuonMisalignmentScenario::components, TrajectoryStateOnSurface::components(), and i.

Referenced by putInEvt().

 {
   reco::GsfComponent5D::ParameterVector pLocS;
   reco::GsfComponent5D::CovarianceMatrix cLocS;
   std::vector<TrajectoryStateOnSurface> components(tsos.components());
   for ( std::vector<TrajectoryStateOnSurface>::const_iterator i=components.begin();
     i!=components.end(); ++i ) {
     states.push_back(reco::GsfComponent5D(i->weight(),i->localParameters().vector(),i->localError().matrix()));
   }
 }

void GsfTrackProducerBase::localParametersFromQpMode	(	const TrajectoryStateOnSurface	tsos,
		AlgebraicVector5 &	parameters,
		AlgebraicSymMatrix55 &	covariance
	)		const

private

local parameters rescaled with q/p from mode

Definition at line 351 of file GsfTrackProducerBase.cc.

References i, TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LocalTrajectoryError::matrix(), SingleGaussianState1D::mean(), GaussianSumUtilities1D::mode(), GaussianSumUtilities1D::modeIsValid(), MultiGaussianStateTransform::multiState1D(), HLT_25ns10e33_v2_cff::parameters, mathSSE::sqrt(), SingleGaussianState1D::variance(), and LocalTrajectoryParameters::vector().

Referenced by computeModeAtTM().

 {
   //
   // parameters and errors from combined state
   //
   parameters = tsos.localParameters().vector();
   covariance = tsos.localError().matrix();
   //
   // mode for parameter 0 (q/p)
   //
   MultiGaussianState1D qpState(MultiGaussianStateTransform::multiState1D(tsos,0));
   GaussianSumUtilities1D qpGS(qpState);
   if ( !qpGS.modeIsValid() )  return;
   double qp = qpGS.mode().mean();
   double varQp = qpGS.mode().variance();
   //
   // replace q/p value and variance, rescale correlation terms
   //   (heuristic procedure - alternative would be mode in 5D ...)
   //
   double VarQpRatio = sqrt(varQp/covariance(0,0));
   parameters(0) = qp;
   covariance(0,0) = varQp;
   for ( int i=1; i<5; ++i )  covariance(i,0) *= VarQpRatio;
 }

void GsfTrackProducerBase::putInEvt	(	edm::Event &	evt,
		const Propagator *	prop,
		const MeasurementTracker *	measTk,
		std::unique_ptr< TrackingRecHitCollection > &	selHits,
		std::unique_ptr< reco::GsfTrackCollection > &	selTracks,
		std::unique_ptr< reco::TrackExtraCollection > &	selTrackExtras,
		std::unique_ptr< reco::GsfTrackExtraCollection > &	selGsfTrackExtras,
		std::unique_ptr< std::vector< Trajectory > > &	selTrajectories,
		AlgoProductCollection &	algoResults,
		TransientTrackingRecHitBuilder const *	hitBuilder,
		const reco::BeamSpot &	bs,
		const TrackerTopology *	ttopo
	)

virtual

Put produced collections in the event.

Definition at line 29 of file GsfTrackProducerBase.cc.

References alongMomentum, AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, reco::TrackBase::appendHitPattern(), assert(), TrajectoryStateOnSurface::components(), computeModeAtTM(), TrajectoryStateOnSurface::curvilinearError(), Trajectory::direction(), fillMode(), fillStates(), Trajectory::firstMeasurement(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), TrajectoryStateOnSurface::globalParameters(), i, customizeTrackingMonitorSeedNumber::idx, edm::ESHandleBase::isValid(), TrajectoryStateOnSurface::isValid(), Trajectory::lastMeasurement(), TrajectoryStateOnSurface::localParameters(), LogDebug, LogTrace, TrajectoryStateOnSurface::magneticField(), Trajectory::measurements(), GlobalTrajectoryParameters::momentum(), eostools::move(), TrackProducerBase< reco::GsfTrack >::mteSrc_, AlCaHLTBitMon_ParallelJobs::p, GlobalTrajectoryParameters::position(), position, edm::Event::put(), LocalTrajectoryParameters::pzSign(), TrajectoryMeasurement::recHit(), TrackProducerBase< reco::GsfTrack >::rTracks_, Trajectory::seedRef(), reco::Track::setExtra(), reco::GsfTrack::setGsfExtra(), reco::TrackExtraBase::setHits(), TrackProducerBase< reco::GsfTrack >::setSecondHitPattern(), DetId::subdetId(), lumiQTWidget::t, SiStripDetId::TEC, TrackProducerBase< reco::GsfTrack >::theSchool, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, DetId::Tracker, TrackProducerBase< reco::GsfTrack >::trajectoryInEvent_, TrajectoryMeasurement::updatedState(), useSplitting, findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by GsfTrackProducer::produce(), and GsfTrackRefitter::produce().

 {
 
   TrackingRecHitRefProd rHits = evt.getRefBeforePut<TrackingRecHitCollection>();
   reco::TrackExtraRefProd rTrackExtras = evt.getRefBeforePut<reco::TrackExtraCollection>();
   reco::GsfTrackExtraRefProd rGsfTrackExtras = evt.getRefBeforePut<reco::GsfTrackExtraCollection>();
   reco::GsfTrackRefProd rTracks = evt.getRefBeforePut<reco::GsfTrackCollection>();
 
   edm::Ref<reco::TrackExtraCollection>::key_type idx = 0;
   edm::Ref<reco::TrackExtraCollection>::key_type hidx = 0;
   edm::Ref<reco::GsfTrackExtraCollection>::key_type idxGsf = 0;
   edm::Ref<reco::GsfTrackCollection>::key_type iTkRef = 0;
   edm::Ref< std::vector<Trajectory> >::key_type iTjRef = 0;
   std::map<unsigned int, unsigned int> tjTkMap;
 
   TSCBLBuilderNoMaterial tscblBuilder;
   
   for(auto & i : algoResults){
     Trajectory * theTraj = i.trajectory;
     if(trajectoryInEvent_) {
       selTrajectories->push_back(*theTraj);
       iTjRef++;
     }
 
     reco::GsfTrack * theTrack = i.track;
     PropagationDirection seedDir = i.pDir;  
     
     LogDebug("TrackProducer") << "In GsfTrackProducerBase::putInEvt - seedDir=" << seedDir;
 
     reco::GsfTrack t = * theTrack;
     selTracks->push_back( t );
     iTkRef++;
 
     // Store indices in local map (starts at 0)
     if(trajectoryInEvent_) tjTkMap[iTjRef-1] = iTkRef-1;
 
     //sets the outermost and innermost TSOSs
     TrajectoryStateOnSurface outertsos;
     TrajectoryStateOnSurface innertsos;
     unsigned int innerId, outerId;
     
     // ---  NOTA BENE: the convention is to sort hits and measurements "along the momentum".
     // This is consistent with innermost and outermost labels only for tracks from LHC collision
     if (theTraj->direction() == alongMomentum) {
       outertsos = theTraj->lastMeasurement().updatedState();
       innertsos = theTraj->firstMeasurement().updatedState();
       outerId = theTraj->lastMeasurement().recHit()->geographicalId().rawId();
       innerId = theTraj->firstMeasurement().recHit()->geographicalId().rawId();
     } else { 
       outertsos = theTraj->firstMeasurement().updatedState();
       innertsos = theTraj->lastMeasurement().updatedState();
       outerId = theTraj->firstMeasurement().recHit()->geographicalId().rawId();
       innerId = theTraj->lastMeasurement().recHit()->geographicalId().rawId();
     }
     //build the TrackExtra
     GlobalPoint v = outertsos.globalParameters().position();
     GlobalVector p = outertsos.globalParameters().momentum();
     math::XYZVector outmom( p.x(), p.y(), p.z() );
     math::XYZPoint  outpos( v.x(), v.y(), v.z() );
     v = innertsos.globalParameters().position();
     p = innertsos.globalParameters().momentum();
     math::XYZVector inmom( p.x(), p.y(), p.z() );
     math::XYZPoint  inpos( v.x(), v.y(), v.z() );
 
     reco::TrackExtraRef teref= reco::TrackExtraRef ( rTrackExtras, idx ++ );
     reco::GsfTrack & track = selTracks->back();
     track.setExtra( teref );
     
     //======= I want to set the second hitPattern here =============
     if (theSchool.isValid())
       {
         edm::Handle<MeasurementTrackerEvent> mte;
         evt.getByToken(mteSrc_, mte);
     // NavigationSetter setter( *theSchool );
     setSecondHitPattern(theTraj,track,prop,&*mte, ttopo);
       }
     //==============================================================
     
     selTrackExtras->push_back( reco::TrackExtra (outpos, outmom, true, inpos, inmom, true,
                          outertsos.curvilinearError(), outerId,
                          innertsos.curvilinearError(), innerId,
                          seedDir, theTraj->seedRef()));
 
 
     reco::TrackExtra & tx = selTrackExtras->back();
 
     // ---  NOTA BENE: the convention is to sort hits and measurements "along the momentum".
     // This is consistent with innermost and outermost labels only for tracks from LHC collisions
     Traj2TrackHits t2t(hitBuilder,false);
     auto ih = selHits->size();
     assert(ih==hidx);
     t2t(*theTraj,*selHits,useSplitting);
     auto ie = selHits->size();
     unsigned int nHitsAdded = 0;
     for (;ih<ie; ++ih) {
       auto const & hit = (*selHits)[ih];
       track.appendHitPattern(hit, *ttopo);
       ++nHitsAdded;
     }
     tx.setHits(rHits, hidx, nHitsAdded);
     hidx += nHitsAdded;
 
     /*
     TrajectoryFitter::RecHitContainer transHits; theTraj->recHitsV(transHits,useSplitting);
     // ---  NOTA BENE: the convention is to sort hits and measurements "along the momentum".
     // This is consistent with innermost and outermost labels only for tracks from LHC collisions
     if (theTraj->direction() == alongMomentum) {
         for(TrajectoryFitter::RecHitContainer::const_iterator j = transHits.begin(); 
                 j != transHits.end(); j++) {
             if ((**j).hit() != 0){
                 TrackingRecHit *hit = (**j).hit()->clone();
                 track.appendHitPattern(*hit, *ttopo);
                 selHits->push_back(hit);
                 tx.add(TrackingRecHitRef(rHits, hidx++));
             }
         }
     }else{
         for(TrajectoryFitter::RecHitContainer::const_iterator j = transHits.end() - 1;
                 j != transHits.begin() - 1; --j) {
             if ((**j).hit() != 0){
                 TrackingRecHit *hit = (**j).hit()->clone();
                 track.appendHitPattern(*hit, *ttopo);
                 selHits->push_back(hit);
                 tx.add(TrackingRecHitRef(rHits, hidx++));
             }
         }
     }
     */
     // ----
 
     std::vector<reco::GsfTangent> tangents;
     const Trajectory::DataContainer& measurements = theTraj->measurements();
     if ( measurements.size()>2 ) {
       tangents.reserve(measurements.size()-2);
       Trajectory::DataContainer::const_iterator ibegin,iend;
       int increment(0);
       if (theTraj->direction() == alongMomentum) {
     ibegin = measurements.begin() + 1;
     iend = measurements.end() - 1;
     increment = 1;
       }
       else {
     ibegin = measurements.end() - 2;
     iend = measurements.begin();
     increment = -1;
       }
       math::XYZPoint position;
       math::XYZVector momentum;
       Measurement1D deltaP;
       // only measurements on "mono" detectors
       for ( Trajectory::DataContainer::const_iterator i=ibegin;
         i!=iend; i+=increment ) {
     if ( i->recHit().get() ) {
       DetId detId(i->recHit()->geographicalId());
       if ( detId.det()==DetId::Tracker ) {
         int subdetId = detId.subdetId();
         if ( subdetId==SiStripDetId::TIB || subdetId==SiStripDetId::TID || 
          subdetId==SiStripDetId::TOB || subdetId==SiStripDetId::TEC ) {
           if ( SiStripDetId(detId).stereo() )  continue;        
         }
       }
     }
     bool valid = computeModeAtTM(*i,position,momentum,deltaP);
     if ( valid ) {
       tangents.push_back(reco::GsfTangent(position,momentum,deltaP));
     }
       }
     }
     
 
     //build the GsfTrackExtra
     std::vector<reco::GsfComponent5D> outerStates;
     outerStates.reserve(outertsos.components().size());
     fillStates(outertsos,outerStates);
     std::vector<reco::GsfComponent5D> innerStates;
     innerStates.reserve(innertsos.components().size());
     fillStates(innertsos,innerStates);
     
 
     reco::GsfTrackExtraRef terefGsf = reco::GsfTrackExtraRef ( rGsfTrackExtras, idxGsf ++ );
     track.setGsfExtra( terefGsf );
     selGsfTrackExtras->push_back( reco::GsfTrackExtra (outerStates, outertsos.localParameters().pzSign(),
                                innerStates, innertsos.localParameters().pzSign(),
                                tangents));
 
     if ( innertsos.isValid() ) {
       GsfPropagatorAdapter gsfProp(AnalyticalPropagator(innertsos.magneticField(),anyDirection));
       TransverseImpactPointExtrapolator tipExtrapolator(gsfProp);
       fillMode(track,innertsos,gsfProp,tipExtrapolator,tscblBuilder,bs);
     }
 
     delete theTrack;
     delete theTraj;
   }
 
   LogTrace("TrackingRegressionTest") << "========== TrackProducer Info ===================";
   LogTrace("TrackingRegressionTest") << "number of finalGsfTracks: " << selTracks->size();
   for (reco::GsfTrackCollection::const_iterator it = selTracks->begin(); it != selTracks->end(); it++) {
     LogTrace("TrackingRegressionTest") << "track's n valid and invalid hit, chi2, pt : " 
                        << it->found() << " , " 
                        << it->lost()  <<" , " 
                        << it->normalizedChi2() << " , "
                        << it->pt() << " , "
                        << it->eta() ;
   }
   LogTrace("TrackingRegressionTest") << "=================================================";
   
 
   rTracks_ = evt.put(std::move(selTracks) );
   evt.put(std::move(selTrackExtras) );
   evt.put(std::move(selGsfTrackExtras) );
   evt.put(std::move(selHits) );
 
   if(trajectoryInEvent_) {
     edm::OrphanHandle<std::vector<Trajectory> > rTrajs = evt.put(std::move(selTrajectories));
 
     // Now Create traj<->tracks association map
     std::unique_ptr<TrajGsfTrackAssociationCollection> trajTrackMap( new TrajGsfTrackAssociationCollection(rTrajs, rTracks_) );
     for ( std::map<unsigned int, unsigned int>::iterator i = tjTkMap.begin(); 
       i != tjTkMap.end(); i++ ) {
       edm::Ref<std::vector<Trajectory> > trajRef( rTrajs, (*i).first );
       edm::Ref<reco::GsfTrackCollection>    tkRef( rTracks_, (*i).second );
       trajTrackMap->insert( edm::Ref<std::vector<Trajectory> >( rTrajs, (*i).first ),
                 edm::Ref<reco::GsfTrackCollection>( rTracks_, (*i).second ) );
     }
     evt.put( std::move(trajTrackMap) );
   }
 }

Member Data Documentation

bool GsfTrackProducerBase::useSplitting

private

Definition at line 68 of file GsfTrackProducerBase.h.

Referenced by putInEvt().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation