DAFTrackProducerAlgorithm Class Reference

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

Definition at line 34 of file DAFTrackProducerAlgorithm.h.

Member Typedef Documentation

using DAFTrackProducerAlgorithm::AlgoProductCollection = typename Base::AlgoProductCollection

Definition at line 38 of file DAFTrackProducerAlgorithm.h.

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

Definition at line 36 of file DAFTrackProducerAlgorithm.h.

using DAFTrackProducerAlgorithm::TrackCollection = typename Base::TrackCollection

Definition at line 37 of file DAFTrackProducerAlgorithm.h.

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

Definition at line 40 of file DAFTrackProducerAlgorithm.h.

Constructor & Destructor Documentation

DAFTrackProducerAlgorithm::DAFTrackProducerAlgorithm

(

const edm::ParameterSet &

conf

)

Definition at line 27 of file DAFTrackProducerAlgorithm.cc.

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

DAFTrackProducerAlgorithm::~DAFTrackProducerAlgorithm ( )

inline

Definition at line 44 of file DAFTrackProducerAlgorithm.h.

References ALCARECOTkAlBeamHalo_cff::filter, trackingPlots::fit, hfClusterShapes_cfi::hits, input, HLTSiStripMonitoring_cff::measurementTracker, convertSQLitetoXML_cfg::output, and HLT_2018_cff::updator.

{}

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

float DAFTrackProducerAlgorithm::calculateNdof ( const Trajectory  vtraj ) const

private

Definition at line 380 of file DAFTrackProducerAlgorithm.cc.

References makeMuonMisalignmentScenario::components, watchdog::const, 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 401 of file DAFTrackProducerAlgorithm.cc.

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().

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

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.

References hfClusterShapes_cfi::hits, 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 285 of file DAFTrackProducerAlgorithm.cc.

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

Referenced by runWithCandidate().

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

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

private

Definition at line 315 of file DAFTrackProducerAlgorithm.cc.

References makeMuonMisalignmentScenario::components, watchdog::const, MillePedeFileConverter_cfg::e, TrajectoryFitter::fit(), edm::OwnVector< T, P >::front(), AnalysisDataFormats_SUSYBSMObjects::hc, hfClusterShapes_cfi::hits, 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 214 of file DAFTrackProducerAlgorithm.cc.

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

Referenced by trackingPlots.Iteration::modules(), and 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 478 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 30 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(), hfClusterShapes_cfi::hits, 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 172 of file DAFTrackProducerAlgorithm.cc.

References heavyFlavorValidationHarvestingSequence_cff::combiner, hfClusterShapes_cfi::hits, 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().