#include <DAFTrackProducerAlgorithm.h>

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 Types
typedef std::pair< Trajectory , std::pair< reco::Track , PropagationDirection > >	AlgoProduct

typedef std::vector< AlgoProduct >	AlgoProductCollection

typedef std::vector < TrajAnnealing >	TrajAnnealingCollection

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 in May 2014 by brondolin

Definition at line 32 of file DAFTrackProducerAlgorithm.h.

Member Typedef Documentation

typedef std::pair<Trajectory*, std::pair<reco::Track*,PropagationDirection> > DAFTrackProducerAlgorithm::AlgoProduct

private

Definition at line 34 of file DAFTrackProducerAlgorithm.h.

typedef std::vector< AlgoProduct > DAFTrackProducerAlgorithm::AlgoProductCollection

private

Definition at line 35 of file DAFTrackProducerAlgorithm.h.

typedef std::vector<TrajAnnealing> DAFTrackProducerAlgorithm::TrajAnnealingCollection

private

Definition at line 36 of file DAFTrackProducerAlgorithm.h.

Constructor & Destructor Documentation

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

Definition at line 28 of file DAFTrackProducerAlgorithm.cc.

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

DAFTrackProducerAlgorithm::~DAFTrackProducerAlgorithm ( )

inline

Definition at line 41 of file DAFTrackProducerAlgorithm.h.

41 {}

Member Function Documentation

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 239 of file DAFTrackProducerAlgorithm.cc.

References ecalcalib_dqm_sourceclient-live_cfg::algo, alongMomentum, 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(), findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by runWithCandidate().

 {
   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(), std::make_pair(theTrack.release(), vtraj.direction()));
     algoResults.push_back(aProduct);
 
     return true;
   } 
   else {
     LogDebug("DAFTrackProducerAlgorithm") <<" BUILDER NOT POSSIBLE: traj is not valid" << std::endl;;
     return false;
   }
 }

float DAFTrackProducerAlgorithm::calculateNdof ( const Trajectory vtraj ) const

private

Definition at line 378 of file DAFTrackProducerAlgorithm.cc.

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

Referenced by runWithCandidate().

 {
 
   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;
 
 }

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

private

Definition at line 402 of file DAFTrackProducerAlgorithm.cc.

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

Referenced by runWithCandidate().

                                                                                                  {
 
   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=0;
       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++;
         }
 
         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;
 }

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

private

Definition at line 139 of file DAFTrackProducerAlgorithm.cc.

References LogDebug, and MultiRecHitCollector::recHits().

Referenced by runWithCandidate().

 {
 
   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);
 
 }

int DAFTrackProducerAlgorithm::countingGoodHits ( const Trajectory traj ) const

private

Definition at line 289 of file DAFTrackProducerAlgorithm.cc.

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

Referenced by runWithCandidate().

                                                                           {
 
   int ngoodhits = 0;
   Trajectory myTraj = traj;
   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;
 
 }

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

private

Definition at line 321 of file DAFTrackProducerAlgorithm.cc.

References makeMuonMisalignmentScenario::components, alignCSCRings::e, TrajectoryFitter::fit(), edm::OwnVector< T, P >::front(), AnalysisDataFormats_SUSYBSMObjects::hc, LogDebug, TkClonerImpl::makeShared(), TrackingRecHit::missing, SiTrackerMultiRecHit::recHits(), and SiTrackerMultiRecHit::weight().

 {
   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(PTrajectoryStateOnDet(),
                                       BasicTrajectorySeed::recHitContainer(),
                                       input.front().seed().direction()),
                                   hits,
                                   TrajectoryStateWithArbitraryError()(curstartingTSOS));
 
   LogDebug("DAFTrackProducerAlgorithm") << "After filtering " << output.size() << " trajectories";
 
 }

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 220 of file DAFTrackProducerAlgorithm.cc.

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

Referenced by runWithCandidate().

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

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

private

Definition at line 481 of file DAFTrackProducerAlgorithm.cc.

References TrackingRecHit::dimension(), TrackingRecHit::geographicalId(), TrackingRecHit::globalPosition(), TrajectoryStateOnSurface::isValid(), TrackingRecHit::isValid(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPosition(), LogTrace, PV3DBase< T, PVType, FrameType >::perp(), and DetId::rawId().

Referenced by checkHits(), and updateHits().

 {
   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();
   }
 
 }

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 33 of file DAFTrackProducerAlgorithm.cc.

References edm::AssociationMap< Tag >::begin(), buildTrack(), calculateNdof(), checkHits(), collectHits(), edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, generateEDF::cont, countingGoodHits(), TrajectorySeed::direction(), Trajectory::direction(), edm::AssociationMap< Tag >::end(), fit(), Trajectory::foundHits(), SiTrackerMultiRecHitUpdator::getAnnealingProgram(), Trajectory::isValid(), edm::Ref< C, T, F >::key(), LogDebug, minHits_, ndof, convertSQLiteXML::ok, Trajectory::seed(), edm::AssociationMap< Tag >::size(), groupFilesInBlocks::temp, and updateHits().

Referenced by DAFTrackProducer::produce().

 {
   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";
 
 }

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

private

Definition at line 179 of file DAFTrackProducerAlgorithm.cc.

References heavyFlavorValidationHarvestingSequence_cff::combiner, MeasurementTrackerEvent::idToDet(), LogDebug, Trajectory::measurements(), PrintHit(), and SiTrackerMultiRecHitUpdator::update().

Referenced by runWithCandidate().

 {
   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);
 }

Member Data Documentation

edm::ParameterSet DAFTrackProducerAlgorithm::conf_

private

Definition at line 102 of file DAFTrackProducerAlgorithm.h.

int DAFTrackProducerAlgorithm::minHits_

private

Definition at line 103 of file DAFTrackProducerAlgorithm.h.

Referenced by runWithCandidate().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation