#include <DAFTrackProducerAlgorithm.h>

Inheritance diagram for DAFTrackProducerAlgorithm:

Public Types
using	AlgoProductCollection = typename Base::AlgoProductCollection

using	Base = AlgoProductTraits< reco::Track >

using	TrackCollection = typename Base::TrackCollection

using	TrajAnnealingCollection = std::vector< TrajAnnealing >

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

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

using	TrackView = edm::View< reco::Track >

Public Member Functions
	DAFTrackProducerAlgorithm (const edm::ParameterSet &conf)

void	runWithCandidate (const TrackingGeometry , const MagneticField , const TrajTrackAssociationCollection &, const MeasurementTrackerEvent measTk, const TrajectoryFitter , const TransientTrackingRecHitBuilder , const MultiRecHitCollector measurementTracker, const SiTrackerMultiRecHitUpdator *, const reco::BeamSpot &, AlgoProductCollection &, TrajAnnealingCollection &, bool, AlgoProductCollection &, AlgoProductCollection &) const
	Run the Final Fit taking TrackCandidates as input. More...

	~DAFTrackProducerAlgorithm ()

Private Member Functions
bool	buildTrack (const Trajectory, AlgoProductCollection &algoResults, float, const reco::BeamSpot &, const reco::TrackRef *) const
	Construct Tracks to be put in the event. More...

float	calculateNdof (const Trajectory vtraj) const

int	checkHits (Trajectory iInitTraj, const Trajectory iFinalTraj) const

std::pair< TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface >	collectHits (const Trajectory vtraj, const MultiRecHitCollector measurementCollector, const MeasurementTrackerEvent measTk) const

int	countingGoodHits (const Trajectory traj) const

void	filter (const TrajectoryFitter fitter, std::vector< Trajectory > &input, int minhits, std::vector< Trajectory > &output, const TransientTrackingRecHitBuilder builder) const

Trajectory	fit (const std::pair< TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface > &hits, const TrajectoryFitter *theFitter, Trajectory vtraj) const
	accomplishes the fitting-smoothing step for each annealing value More...

void	PrintHit (const TrackingRecHit *const &hit, TrajectoryStateOnSurface &tsos) const

std::pair< TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface >	updateHits (const Trajectory vtraj, const SiTrackerMultiRecHitUpdator updator, const MeasurementTrackerEvent theMTE, double annealing) const

Private Attributes
edm::ParameterSet	conf_

int	minHits_

Detailed Description

All is needed to run the deterministic annealing algorithm. Ported from ORCA.

Author: tropiano, genta \review in May 2014 by brondolin

Definition at line 34 of file DAFTrackProducerAlgorithm.h.

Member Typedef Documentation

◆ AlgoProductCollection

using DAFTrackProducerAlgorithm::AlgoProductCollection = typename Base::AlgoProductCollection

Definition at line 38 of file DAFTrackProducerAlgorithm.h.

◆ Base

using DAFTrackProducerAlgorithm::Base = AlgoProductTraits<reco::Track>

Definition at line 36 of file DAFTrackProducerAlgorithm.h.

◆ TrackCollection

using DAFTrackProducerAlgorithm::TrackCollection = typename Base::TrackCollection

Definition at line 37 of file DAFTrackProducerAlgorithm.h.

◆ TrajAnnealingCollection

using DAFTrackProducerAlgorithm::TrajAnnealingCollection = std::vector<TrajAnnealing>

Definition at line 40 of file DAFTrackProducerAlgorithm.h.

Constructor & Destructor Documentation

◆ DAFTrackProducerAlgorithm()

DAFTrackProducerAlgorithm::DAFTrackProducerAlgorithm ( const edm::ParameterSet & conf )

Definition at line 27 of file DAFTrackProducerAlgorithm.cc.

28 : conf_(conf), minHits_(conf.getParameter<int>("MinHits")) {}

◆ ~DAFTrackProducerAlgorithm()

DAFTrackProducerAlgorithm::~DAFTrackProducerAlgorithm ( )

inline

Definition at line 44 of file DAFTrackProducerAlgorithm.h.

44 {}

Member Function Documentation

◆ buildTrack()

bool DAFTrackProducerAlgorithm::buildTrack	(	const Trajectory	vtraj,
		AlgoProductCollection &	algoResults,
		float	ndof,
		const reco::BeamSpot &	bs,
		const reco::TrackRef *	BeforeDAFTrack
	)		const

private

Construct Tracks to be put in the event.

Definition at line 229 of file DAFTrackProducerAlgorithm.cc.

                                                                                      {
   LogDebug("DAFTrackProducerAlgorithm") << " BUILDER " << std::endl;
   ;
   TrajectoryStateOnSurface innertsos;
  
   if (vtraj.isValid()) {
     std::unique_ptr<Trajectory> theTraj(new Trajectory(vtraj));
  
     if (vtraj.direction() == alongMomentum) {
       //if (theTraj->direction() == oppositeToMomentum) {
       innertsos = vtraj.firstMeasurement().updatedState();
     } else {
       innertsos = vtraj.lastMeasurement().updatedState();
     }
  
     TSCBLBuilderNoMaterial tscblBuilder;
     TrajectoryStateClosestToBeamLine tscbl = tscblBuilder(*(innertsos.freeState()), bs);
  
     if (tscbl.isValid() == false)
       return false;
  
     GlobalPoint v = tscbl.trackStateAtPCA().position();
     math::XYZPoint pos(v.x(), v.y(), v.z());
     GlobalVector p = tscbl.trackStateAtPCA().momentum();
     math::XYZVector mom(p.x(), p.y(), p.z());
  
     //    LogDebug("TrackProducer") <<v<<p<<std::endl;
  
     auto algo = (*BeforeDAFTrack)->algo();
     std::unique_ptr<reco::Track> theTrack(new reco::Track(vtraj.chiSquared(),
                                                           ndof,  //in the DAF the ndof is not-integer
                                                           pos,
                                                           mom,
                                                           tscbl.trackStateAtPCA().charge(),
                                                           tscbl.trackStateAtPCA().curvilinearError(),
                                                           algo));
     theTrack->setQualityMask((*BeforeDAFTrack)->qualityMask());
  
     AlgoProduct aProduct{theTraj.release(), theTrack.release(), vtraj.direction(), 0};
     algoResults.push_back(aProduct);
  
     return true;
   } else {
     LogDebug("DAFTrackProducerAlgorithm") << " BUILDER NOT POSSIBLE: traj is not valid" << std::endl;
     ;
     return false;
   }
 }

References alongMomentum, cms::cuda::bs, FreeTrajectoryState::charge(), Trajectory::chiSquared(), FreeTrajectoryState::curvilinearError(), Trajectory::direction(), Trajectory::firstMeasurement(), TrajectoryStateOnSurface::freeState(), TrajectoryStateClosestToBeamLine::isValid(), Trajectory::isValid(), Trajectory::lastMeasurement(), LogDebug, FreeTrajectoryState::momentum(), ndof, AlCaHLTBitMon_ParallelJobs::p, FreeTrajectoryState::position(), TrajectoryStateClosestToBeamLine::trackStateAtPCA(), TrajectoryMeasurement::updatedState(), and findQualityFiles::v.

Referenced by runWithCandidate().

◆ calculateNdof()

float DAFTrackProducerAlgorithm::calculateNdof ( const Trajectory vtraj ) const

private

Definition at line 376 of file DAFTrackProducerAlgorithm.cc.

                                                                            {
   if (!vtraj.isValid())
     return 0;
   float ndof = 0;
   const std::vector<TrajectoryMeasurement>& meas = vtraj.measurements();
   for (std::vector<TrajectoryMeasurement>::const_iterator iter = meas.begin(); iter != meas.end(); iter++) {
     if (iter->recHit()->isValid()) {
       SiTrackerMultiRecHit const& mHit = dynamic_cast<SiTrackerMultiRecHit const&>(*iter->recHit());
       std::vector<const TrackingRecHit*> components = mHit.recHits();
       int iComp = 0;
       for (std::vector<const TrackingRecHit*>::const_iterator iter2 = components.begin(); iter2 != components.end();
            iter2++, iComp++) {
         if ((*iter2)->isValid())
           ndof += ((*iter2)->dimension()) * mHit.weight(iComp);
       }
     }
   }
  
   return ndof - 5;
 }

References makeMuonMisalignmentScenario::components, Trajectory::isValid(), Trajectory::measurements(), ndof, SiTrackerMultiRecHit::recHits(), and SiTrackerMultiRecHit::weight().

Referenced by runWithCandidate().

◆ checkHits()

int DAFTrackProducerAlgorithm::checkHits	(	Trajectory	iInitTraj,
		const Trajectory	iFinalTraj
	)		const

private

Definition at line 397 of file DAFTrackProducerAlgorithm.cc.

                                                                                                 {
   std::vector<TrajectoryMeasurement> initmeasurements = iInitTraj.measurements();
   std::vector<TrajectoryMeasurement> finalmeasurements = iFinalTraj.measurements();
   std::vector<TrajectoryMeasurement>::iterator jmeas;
   int nSame = 0;
   int ihit = 0;
  
   if (initmeasurements.empty() || finalmeasurements.empty()) {
     LogDebug("DAFTrackProducerAlgorithm") << "Initial or Final Trajectory empty.";
     return 0;
   }
  
   if (initmeasurements.size() != finalmeasurements.size()) {
     LogDebug("DAFTrackProducerAlgorithm")
         << "Initial Trajectory size(" << initmeasurements.size() << " hits) "
         << "is different to final traj size (" << finalmeasurements.size() << ")! No checkHits possible! ";
     return 0;
   }
  
   for (jmeas = initmeasurements.begin(); jmeas != initmeasurements.end(); jmeas++) {
     const TrackingRecHit* initHit = jmeas->recHit()->hit();
     if (!initHit->isValid() && ihit == 0)
       continue;
  
     if (initHit->isValid()) {
       TrajectoryMeasurement imeas = finalmeasurements.at(ihit);
       const TrackingRecHit* finalHit = imeas.recHit()->hit();
       const TrackingRecHit* MaxWeightHit = nullptr;
       float maxweight = 0;
  
       const SiTrackerMultiRecHit* mrh = dynamic_cast<const SiTrackerMultiRecHit*>(finalHit);
       if (mrh) {
         std::vector<const TrackingRecHit*> components = mrh->recHits();
         std::vector<const TrackingRecHit*>::const_iterator icomp;
         int hitcounter = 0;
  
         for (icomp = components.begin(); icomp != components.end(); icomp++) {
           if ((*icomp)->isValid()) {
             double weight = mrh->weight(hitcounter);
             if (weight > maxweight) {
               MaxWeightHit = *icomp;
               maxweight = weight;
             }
           }
           hitcounter++;
         }
         if (!MaxWeightHit)
           continue;
  
         auto myref1 = reinterpret_cast<const BaseTrackerRecHit*>(initHit)->firstClusterRef();
         auto myref2 = reinterpret_cast<const BaseTrackerRecHit*>(MaxWeightHit)->firstClusterRef();
  
         if (myref1 == myref2) {
           nSame++;
         } else {
           LogDebug("DAFTrackProducerAlgorithm") << "diverso hit!" << std::endl;
           TrajectoryStateOnSurface dummState;
           LogTrace("DAFTrackProducerAlgorithm") << "  This hit was:\n ";
           PrintHit(initHit, dummState);
           LogTrace("DAFTrackProducerAlgorithm") << "  instead now is:\n ";
           PrintHit(MaxWeightHit, dummState);
         }
       }
     } else {
       TrajectoryMeasurement imeas = finalmeasurements.at(ihit);
       const TrackingRecHit* finalHit = imeas.recHit()->hit();
       if (!finalHit->isValid()) {
         nSame++;
       }
     }
  
     ihit++;
   }
  
   return nSame;
 }

References makeMuonMisalignmentScenario::components, TrackingRecHit::hit(), TrackingRecHit::isValid(), LogDebug, LogTrace, Trajectory::measurements(), PrintHit(), TrajectoryMeasurement::recHit(), SiTrackerMultiRecHit::recHits(), SiTrackerMultiRecHit::weight(), and mps_merge::weight.

Referenced by runWithCandidate().

◆ collectHits()

std::pair< TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface > DAFTrackProducerAlgorithm::collectHits	(	const Trajectory	vtraj,
		const MultiRecHitCollector *	measurementCollector,
		const MeasurementTrackerEvent *	measTk
	)		const

private

Definition at line 137 of file DAFTrackProducerAlgorithm.cc.

                                                  {
   LogDebug("DAFTrackProducerAlgorithm") << "Calling DAFTrackProducerAlgorithm::collectHits";
  
   //getting the traj measurements from the MeasurementCollector
   int nHits = 0;
   TransientTrackingRecHit::RecHitContainer hits;
   std::vector<TrajectoryMeasurement> collectedmeas = measurementCollector->recHits(vtraj, measTk);
  
   //if the MeasurementCollector is empty, make an "empty" pair
   //else taking the collected measured hits and building the pair
   if (collectedmeas.empty())
     return std::make_pair(TransientTrackingRecHit::RecHitContainer(), TrajectoryStateOnSurface());
  
   for (std::vector<TrajectoryMeasurement>::const_iterator iter = collectedmeas.begin(); iter != collectedmeas.end();
        iter++) {
     nHits++;
     hits.push_back(iter->recHit());
   }
  
   //TrajectoryStateWithArbitraryError() == Creates a TrajectoryState with the same parameters
   //     as the input one, but with "infinite" errors, i.e. errors so big that they don't
   //     bias a fit starting with this state.
   //return std::make_pair(hits,TrajectoryStateWithArbitraryError()(collectedmeas.front().predictedState()));
  
   // I do not have to rescale the error because it is already rescaled in the fit code
   TrajectoryStateOnSurface initialStateFromTrack = collectedmeas.front().predictedState();
  
   LogDebug("DAFTrackProducerAlgorithm")
       << "Pair (hits, TSOS)  with TSOS predicted(collectedmeas.front().predictedState())";
   return std::make_pair(hits, initialStateFromTrack);
 }

References hfClusterShapes_cfi::hits, LogDebug, and MultiRecHitCollector::recHits().

Referenced by runWithCandidate().

◆ countingGoodHits()

int DAFTrackProducerAlgorithm::countingGoodHits ( const Trajectory traj ) const

private

Definition at line 282 of file DAFTrackProducerAlgorithm.cc.

                                                                            {
   int ngoodhits = 0;
   std::vector<TrajectoryMeasurement> vtm = traj.measurements();
  
   for (std::vector<TrajectoryMeasurement>::const_iterator tm = vtm.begin(); tm != vtm.end(); tm++) {
     //if the rechit is valid
     if (tm->recHit()->isValid()) {
       SiTrackerMultiRecHit const& mHit = dynamic_cast<SiTrackerMultiRecHit const&>(*tm->recHit());
       std::vector<const TrackingRecHit*> components = mHit.recHits();
  
       int iComp = 0;
  
       for (std::vector<const TrackingRecHit*>::const_iterator iter = components.begin(); iter != components.end();
            iter++, iComp++) {
         //if there is at least one component with weight higher than 1e-6 then the hit is not an outlier
         if (mHit.weight(iComp) > 1e-6) {
           ngoodhits++;
           iComp++;
           break;
         }
       }
     }
   }
  
   LogDebug("DAFTrackProducerAlgorithm") << "Original number of valid hits " << traj.foundHits()
                                         << " -> hit with good weight (>1e-6) are " << ngoodhits;
   return ngoodhits;
 }

References makeMuonMisalignmentScenario::components, MillePedeFileConverter_cfg::e, Trajectory::foundHits(), LogDebug, Trajectory::measurements(), SiTrackerMultiRecHit::recHits(), and SiTrackerMultiRecHit::weight().

Referenced by runWithCandidate().

◆ filter()

void DAFTrackProducerAlgorithm::filter	(	const TrajectoryFitter *	fitter,
		std::vector< Trajectory > &	input,
		int	minhits,
		std::vector< Trajectory > &	output,
		const TransientTrackingRecHitBuilder *	builder
	)		const

private

Definition at line 312 of file DAFTrackProducerAlgorithm.cc.

                                                                                             {
   if (input.empty())
     return;
  
   int ngoodhits = 0;
   std::vector<TrajectoryMeasurement> vtm = input[0].measurements();
   TransientTrackingRecHit::RecHitContainer hits;
  
   //count the number of non-outlier and non-invalid hits
   for (std::vector<TrajectoryMeasurement>::reverse_iterator tm = vtm.rbegin(); tm != vtm.rend(); tm++) {
     //if the rechit is valid
     if (tm->recHit()->isValid()) {
       SiTrackerMultiRecHit const& mHit = dynamic_cast<SiTrackerMultiRecHit const&>(*tm->recHit());
       std::vector<const TrackingRecHit*> components = mHit.recHits();
       int iComp = 0;
       bool isGood = false;
       for (std::vector<const TrackingRecHit*>::const_iterator iter = components.begin(); iter != components.end();
            iter++, iComp++) {
         //if there is at least one component with weight higher than 1e-6 then the hit is not an outlier
         if (mHit.weight(iComp) > 1e-6) {
           ngoodhits++;
           iComp++;
           isGood = true;
           break;
         }
       }
       if (isGood) {
         TkClonerImpl hc = static_cast<TkTransientTrackingRecHitBuilder const*>(builder)->cloner();
         auto tempHit = hc.makeShared(tm->recHit(), tm->updatedState());
         hits.push_back(tempHit);
       } else
         hits.push_back(std::make_shared<InvalidTrackingRecHit>(*tm->recHit()->det(), TrackingRecHit::missing));
     } else {
       TkClonerImpl hc = static_cast<TkTransientTrackingRecHitBuilder const*>(builder)->cloner();
       auto tempHit = hc.makeShared(tm->recHit(), tm->updatedState());
       hits.push_back(tempHit);
     }
   }
  
   LogDebug("DAFTrackProducerAlgorithm") << "Original number of valid hits " << input[0].foundHits()
                                         << "; after filtering " << ngoodhits;
   if (ngoodhits > input[0].foundHits())
     edm::LogError("DAFTrackProducerAlgorithm")
         << "Something wrong: the number of good hits from DAFTrackProducerAlgorithm::filter " << ngoodhits
         << " is higher than the original one " << input[0].foundHits();
  
   if (ngoodhits < minhits)
     return;
  
   TrajectoryStateOnSurface curstartingTSOS = input.front().lastMeasurement().updatedState();
   LogDebug("DAFTrackProducerAlgorithm") << "starting tsos for final refitting " << curstartingTSOS;
   //curstartingTSOS.rescaleError(100);
  
   output = fitter->fit(TrajectorySeed({}, {}, input.front().seed().direction()),
                        hits,
                        TrajectoryStateWithArbitraryError()(curstartingTSOS));
  
   LogDebug("DAFTrackProducerAlgorithm") << "After filtering " << output.size() << " trajectories";
 }

References makeMuonMisalignmentScenario::components, MillePedeFileConverter_cfg::e, TrajectoryFitter::fit(), hfClusterShapes_cfi::hits, input, LogDebug, TkClonerImpl::makeShared(), TrackingRecHit::missing, convertSQLitetoXML_cfg::output, SiTrackerMultiRecHit::recHits(), and SiTrackerMultiRecHit::weight().

◆ fit()

Trajectory DAFTrackProducerAlgorithm::fit	(	const std::pair< TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface > &	hits,
		const TrajectoryFitter *	theFitter,
		Trajectory	vtraj
	)		const

private

accomplishes the fitting-smoothing step for each annealing value

Definition at line 214 of file DAFTrackProducerAlgorithm.cc.

                             {
   //creating a new trajectory starting from the direction of the seed of the input one and the hits
   Trajectory newVec = theFitter->fitOne(TrajectorySeed({}, {}, vtraj.seed().direction()), hits.first, hits.second);
  
   if (newVec.isValid())
     return newVec;
   else {
     LogDebug("DAFTrackProducerAlgorithm") << "Fit no valid.";
     return Trajectory();
   }
 }

References TrajectorySeed::direction(), TrajectoryFitter::fitOne(), hfClusterShapes_cfi::hits, Trajectory::isValid(), LogDebug, and Trajectory::seed().

Referenced by trackingPlots.Iteration::modules(), and runWithCandidate().

◆ PrintHit()

void DAFTrackProducerAlgorithm::PrintHit	(	const TrackingRecHit *const &	hit,
		TrajectoryStateOnSurface &	tsos
	)		const

private

Definition at line 474 of file DAFTrackProducerAlgorithm.cc.

                                                                                                                {
   if (hit->isValid()) {
     LogTrace("DAFTrackProducerAlgorithm")
         << "  Valid Hit with DetId " << hit->geographicalId().rawId() << " and dim:" << hit->dimension()
         << " local position " << hit->localPosition() << " global position " << hit->globalPosition() << " and r "
         << hit->globalPosition().perp();
     if (tsos.isValid())
       LogTrace("DAFTrackProducerAlgorithm") << "  TSOS combtsos " << tsos.localPosition();
  
   } else {
     LogTrace("DAFTrackProducerAlgorithm") << "  Invalid Hit with DetId " << hit->geographicalId().rawId();
   }
 }

References TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localPosition(), and LogTrace.

Referenced by checkHits(), and updateHits().

◆ runWithCandidate()

void DAFTrackProducerAlgorithm::runWithCandidate	(	const TrackingGeometry *	theG,
		const MagneticField *	theMF,
		const TrajTrackAssociationCollection &	TTmap,
		const MeasurementTrackerEvent *	measTk,
		const TrajectoryFitter *	theFitter,
		const TransientTrackingRecHitBuilder *	builder,
		const MultiRecHitCollector *	measurementTracker,
		const SiTrackerMultiRecHitUpdator *	updator,
		const reco::BeamSpot &	bs,
		AlgoProductCollection &	algoResults,
		TrajAnnealingCollection &	trajann,
		bool	TrajAnnSaving_,
		AlgoProductCollection &	algoResultsBeforeDAF,
		AlgoProductCollection &	algoResultsAfterDAF
	)		const

Run the Final Fit taking TrackCandidates as input.

Definition at line 30 of file DAFTrackProducerAlgorithm.cc.

                                                                                                    {
   LogDebug("DAFTrackProducerAlgorithm") << "Size of map: " << TTmap.size() << "\n";
  
   int cont = 0;
   int nTracksChanged = 0;
  
   for (TrajTrackAssociationCollection::const_iterator itTTmap = TTmap.begin(); itTTmap != TTmap.end(); itTTmap++) {
     const edm::Ref<std::vector<Trajectory> > BeforeDAFTraj = itTTmap->key;
     std::vector<TrajectoryMeasurement> BeforeDAFTrajMeas = BeforeDAFTraj->measurements();
     const reco::TrackRef BeforeDAFTrack = itTTmap->val;
  
     float ndof = 0;
     Trajectory CurrentTraj;
  
     if (BeforeDAFTraj->isValid()) {
       LogDebug("DAFTrackProducerAlgorithm") << "The trajectory #" << cont + 1 << " is valid. \n";
  
       //getting the MultiRecHit collection and the trajectory with a first fit-smooth round
       std::pair<TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface> hits =
           collectHits(*BeforeDAFTraj, measurementCollector, &*measTk);
  
       //new initial fit
       CurrentTraj = fit(hits, theFitter, *BeforeDAFTraj);
  
       //starting the annealing program
       for (std::vector<double>::const_iterator ian = updator->getAnnealingProgram().begin();
            ian != updator->getAnnealingProgram().end();
            ian++) {
         if (CurrentTraj.isValid()) {
           LogDebug("DAFTrackProducerAlgorithm") << "Seed direction is " << CurrentTraj.seed().direction()
                                                 << ".Traj direction is " << CurrentTraj.direction() << std::endl;
  
           //updating MultiRecHits and fit-smooth again
           std::pair<TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface> curiterationhits =
               updateHits(CurrentTraj, updator, &*measTk, *ian);
           if (curiterationhits.first.size() < 3) {
             LogDebug("DAFTrackProducerAlgorithm")
                 << "Rejecting trajectory with " << curiterationhits.first.size() << " hits" << std::endl;
             CurrentTraj = Trajectory();
             break;
           }
  
           CurrentTraj = fit(curiterationhits, theFitter, CurrentTraj);
  
           //saving trajectory for each annealing cycle ...
           if (TrajAnnSaving_) {
             TrajAnnealing temp(CurrentTraj, *ian);
             trajann.push_back(temp);
           }
  
           LogDebug("DAFTrackProducerAlgorithm") << "done annealing value " << (*ian);
  
         } else
           break;
       }  //end annealing cycle
  
       int percOfHitsUnchangedAfterDAF =
           (1. * checkHits(*BeforeDAFTraj, CurrentTraj) / (1. * BeforeDAFTrajMeas.size())) * 100.;
       LogDebug("DAFTrackProducerAlgorithm")
           << "Ended annealing program with " << percOfHitsUnchangedAfterDAF << " unchanged." << std::endl;
  
       //computing the ndof keeping into account the weights
       ndof = calculateNdof(CurrentTraj);
  
       //checking if the trajectory has the minimum number of valid hits ( weight (>1e-6) )
       //in order to remove tracks with too many outliers.
       int goodHits = countingGoodHits(CurrentTraj);
  
       if (goodHits >= minHits_) {
         bool ok = buildTrack(CurrentTraj, algoResults, ndof, bs, &BeforeDAFTrack);
         // or filtered?
         if (ok)
           cont++;
  
         //saving tracks before and after DAF
         if ((100. - percOfHitsUnchangedAfterDAF) > 0.) {
           bool okBefore = buildTrack(*BeforeDAFTraj, algoResultsBeforeDAF, ndof, bs, &BeforeDAFTrack);
           bool okAfter = buildTrack(CurrentTraj, algoResultsAfterDAF, ndof, bs, &BeforeDAFTrack);
           if (okBefore && okAfter)
             nTracksChanged++;
         }
       } else {
         LogDebug("DAFTrackProducerAlgorithm") << "Rejecting trajectory with " << CurrentTraj.foundHits() << " hits";
       }
     }  //end run on track collection
     else {
       LogDebug("DAFTrackProducerAlgorithm") << "Rejecting empty trajectory" << std::endl;
     }
   }  //end run on track collection
  
   LogDebug("DAFTrackProducerAlgorithm") << "Number of Tracks found:   " << cont << "\n";
   LogDebug("DAFTrackProducerAlgorithm") << "Number of Tracks changed: " << nTracksChanged << "\n";
 }

References cms::cuda::bs, buildTrack(), calculateNdof(), checkHits(), collectHits(), edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, generateEDF::cont, countingGoodHits(), TrajectorySeed::direction(), Trajectory::direction(), fit(), Trajectory::foundHits(), hfClusterShapes_cfi::hits, Trajectory::isValid(), edm::Ref< C, T, F >::key(), LogDebug, minHits_, ndof, convertSQLiteXML::ok, Trajectory::seed(), groupFilesInBlocks::temp, updateHits(), and HLT_FULL_cff::updator.

Referenced by DAFTrackProducer::produce().

◆ updateHits()

std::pair< TransientTrackingRecHit::RecHitContainer, TrajectoryStateOnSurface > DAFTrackProducerAlgorithm::updateHits	(	const Trajectory	vtraj,
		const SiTrackerMultiRecHitUpdator *	updator,
		const MeasurementTrackerEvent *	theMTE,
		double	annealing
	)		const

private

Definition at line 172 of file DAFTrackProducerAlgorithm.cc.

                             {
   LogDebug("DAFTrackProducerAlgorithm") << "Calling DAFTrackProducerAlgorithm::updateHits";
   TransientTrackingRecHit::RecHitContainer hits;
   std::vector<TrajectoryMeasurement> vmeas = vtraj.measurements();
   std::vector<TrajectoryMeasurement>::reverse_iterator imeas;
   unsigned int hitcounter = 1;
  
   //I run inversely on the trajectory obtained and update the state
   for (imeas = vmeas.rbegin(); imeas != vmeas.rend(); imeas++, hitcounter++) {
     DetId id = imeas->recHit()->geographicalId();
     MeasurementDetWithData measDet = theMTE->idToDet(id);
  
     TrajectoryStateCombiner combiner;
     TrajectoryStateOnSurface combtsos;
     if (hitcounter == vmeas.size())
       combtsos = imeas->predictedState();  //fwd
     else if (hitcounter == 1)
       combtsos = imeas->backwardPredictedState();  //bwd
     else
       combtsos = combiner(imeas->backwardPredictedState(), imeas->predictedState());
  
     PrintHit(&*imeas->recHit(), combtsos);
     if (imeas->recHit()->isValid()) {
       TransientTrackingRecHit::RecHitPointer updated = updator->update(imeas->recHit(), combtsos, measDet, annealing);
       hits.push_back(updated);
     } else {
       hits.push_back(imeas->recHit());
     }
   }
  
   TrajectoryStateOnSurface updatedStateFromTrack = vmeas.back().predictedState();
  
   //return std::make_pair(hits,TrajectoryStateWithArbitraryError()(vmeas.back().updatedState()));
   LogDebug("DAFTrackProducerAlgorithm") << "Pair (hits, TSOS)  with TSOS predicted (vmeas.back().predictedState())";
  
   return std::make_pair(hits, updatedStateFromTrack);
 }