InOutConversionTrackFinder Class Reference

#include <InOutConversionTrackFinder.h>

Inheritance diagram for InOutConversionTrackFinder:

Public Member Functions
	InOutConversionTrackFinder (const edm::ParameterSet &config, const BaseCkfTrajectoryBuilder *trajectoryBuilder)
std::vector< Trajectory >	tracks (const TrajectorySeedCollection &seeds, TrackCandidateCollection &candidate) const override
	~InOutConversionTrackFinder () override
Public Member Functions inherited from ConversionTrackFinder
	ConversionTrackFinder (const edm::ParameterSet &config, const BaseCkfTrajectoryBuilder *trajectoryBuilder)
void	setEventSetup (const edm::EventSetup &es)
	Initialize EventSetup objects at each event. More...
virtual	~ConversionTrackFinder ()

Private Attributes
RedundantSeedCleaner *	theSeedCleaner_
TrajectoryCleanerBySharedHits *	theTrajectoryCleaner_

Additional Inherited Members
Protected Attributes inherited from ConversionTrackFinder
const BaseCkfTrajectoryBuilder *	theCkfTrajectoryBuilder_
std::unique_ptr< TransientInitialStateEstimator >	theInitialState_
const MeasurementTracker *	theMeasurementTracker_
std::string	theMeasurementTrackerName_
const MagneticField *	theMF_
edm::ESHandle< Propagator >	thePropagator_
const TrackerGeometry *	theTrackerGeom_
KFUpdator *	theUpdator_
bool	useSplitHits_

Detailed Description

Author

Nancy Marinelli, U. of Notre Dame, US

Definition at line 28 of file InOutConversionTrackFinder.h.

Constructor & Destructor Documentation

InOutConversionTrackFinder::InOutConversionTrackFinder	(	const edm::ParameterSet &	config,
		const BaseCkfTrajectoryBuilder *	trajectoryBuilder
	)

Definition at line 20 of file InOutConversionTrackFinder.cc.

References Exception, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theSeedCleaner_, and theTrajectoryCleaner_.

    : ConversionTrackFinder(conf, trajectoryBuilder) {
  theTrajectoryCleaner_ = new TrajectoryCleanerBySharedHits(conf);

  // get the seed cleaner
  std::string cleaner = conf.getParameter<std::string>("InOutRedundantSeedCleaner");
  if (cleaner == "CachingSeedCleanerBySharedInput") {
    theSeedCleaner_ = new CachingSeedCleanerBySharedInput();
  } else if (cleaner == "none") {
    theSeedCleaner_ = nullptr;
  } else {
    throw cms::Exception("InOutRedundantSeedCleaner not found, please use CachingSeedCleanerBySharedInput or none",
                         cleaner);
  }
}

InOutConversionTrackFinder::~InOutConversionTrackFinder ( )

override

Definition at line 37 of file InOutConversionTrackFinder.cc.

References theSeedCleaner_, and theTrajectoryCleaner_.

                                                        {
  delete theTrajectoryCleaner_;
  if (theSeedCleaner_)
    delete theSeedCleaner_;
}

Member Function Documentation

std::vector< Trajectory > InOutConversionTrackFinder::tracks ( const TrajectorySeedCollection &  seeds, TrackCandidateCollection &  candidate  ) const

overridevirtual

Implements ConversionTrackFinder.

Definition at line 43 of file InOutConversionTrackFinder.cc.

References RedundantSeedCleaner::add(), alongMomentum, TrajectoryCleanerBySharedHits::clean(), TrackingRecHit::det(), RedundantSeedCleaner::done(), edm::OwnVector< T, P >::front(), RedundantSeedCleaner::good(), BaseCkfTrajectoryBuilder::hitBuilder(), mps_fire::i, RedundantSeedCleaner::init(), TrajectoryStateOnSurface::isValid(), LogDebug, trajectoryStateTransform::persistentState(), precomputed_value_sort(), Propagator::propagate(), TrackingRecHit::rawId(), FastTrackerRecHitMaskProducer_cfi::recHits, mps_fire::result, GeomDet::surface(), ConversionTrackFinder::theCkfTrajectoryBuilder_, ConversionTrackFinder::theInitialState_, ConversionTrackFinder::thePropagator_, theSeedCleaner_, theTrajectoryCleaner_, TrajectoryBuilder::trajectories(), and ConversionTrackFinder::useSplitHits_.

                                                                                                     {
  //  std::cout << " InOutConversionTrackFinder::tracks getting " <<  inOutSeeds.size() << " In-Out seeds " << "\n";

  std::vector<Trajectory> tmpO;
  tmpO.erase(tmpO.begin(), tmpO.end());

  std::vector<Trajectory> result;
  result.erase(result.begin(), result.end());

  std::vector<Trajectory> rawResult;
  if (theSeedCleaner_)
    theSeedCleaner_->init(&rawResult);

  // Loop over the seeds
  int goodSeed = 0;
  for (TrajectorySeedCollection::const_iterator iSeed = inOutSeeds.begin(); iSeed != inOutSeeds.end(); iSeed++) {
    if (!theSeedCleaner_ || theSeedCleaner_->good(&(*iSeed))) {
      goodSeed++;

      LogDebug("InOutConversionTrackFinder")
          << " InOutConversionTrackFinder::tracks hits in the seed " << iSeed->nHits() << "\n";
      LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder::tracks seed starting state position  "
                                             << iSeed->startingState().parameters().position() << " momentum "
                                             << iSeed->startingState().parameters().momentum() << " charge "
                                             << iSeed->startingState().parameters().charge() << "\n";
      LogDebug("InOutConversionTrackFinder")
          << " InOutConversionTrackFinder::tracks seed  starting state para, vector  "
          << iSeed->startingState().parameters().vector() << "\n";

      std::vector<Trajectory> theTmpTrajectories;

      theTmpTrajectories = theCkfTrajectoryBuilder_->trajectories(*iSeed);

      LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder::track returned "
                                             << theTmpTrajectories.size() << " trajectories for this seed "
                                             << "\n";

      theTrajectoryCleaner_->clean(theTmpTrajectories);

      for (std::vector<Trajectory>::const_iterator it = theTmpTrajectories.begin(); it != theTmpTrajectories.end();
           it++) {
        if (it->isValid()) {
          rawResult.push_back(*it);
          if (theSeedCleaner_)
            theSeedCleaner_->add(&(*it));
        }
      }
    }
  }  // end loop over the seeds

  LogDebug("InOutConversionTrackFinder") << "InOutConversionTrackFinder::track Good seeds " << goodSeed << "\n";
  LogDebug("InOutConversionTrackFinder") << "InOutConversionTrackFinder::track rawResult size after cleaning "
                                         << rawResult.size() << "\n";

  if (theSeedCleaner_)
    theSeedCleaner_->done();

  std::vector<Trajectory> unsmoothedResult;
  LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder::track  Start second cleaning "
                                         << "\n";
  theTrajectoryCleaner_->clean(rawResult);
  LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder::track rawResult size after cleaning "
                                         << rawResult.size() << "\n";

  int tra = 0;
  for (std::vector<Trajectory>::const_iterator itraw = rawResult.begin(); itraw != rawResult.end(); itraw++) {
    tra++;
    LogDebug("InOutConversionTrackFinder")
        << " InOutConversionTrackFinder looping of rawResult after cleaning " << tra << "\n";
    if ((*itraw).isValid()) {
      // unsmoothedResult.push_back( *itraw);
      tmpO.push_back(*itraw);
      LogDebug("InOutConversionTrackFinder")
          << " InOutConversionTrackFinder::track rawResult num of valid recHits per trajectory " << (*itraw).foundHits()
          << "\n";
    }
  }

  LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder  tmpO size " << tmpO.size()
                                         << " before sorting "
                                         << "\n";
  //  for (std::vector<Trajectory>::const_iterator it =tmpO.begin(); it != tmpO.end(); it++) {
  // LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder  tmpO num of hits " << (*it).foundHits() << " before ordering " << "\n";
  //}

  precomputed_value_sort(tmpO.begin(), tmpO.end(), ExtractNumOfHits());

  LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder  tmpO after sorting "
                                         << "\n";
  //  for (std::vector<Trajectory>::const_iterator it =tmpO.begin(); it != tmpO.end(); it++) {
  // LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder  tmpO  num of hits " << (*it).foundHits() << "\n";
  // }

  for (int i = tmpO.size() - 1; i >= 0; i--) {
    unsmoothedResult.push_back(tmpO[i]);
  }
  LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder  unsmoothedResult size  "
                                         << unsmoothedResult.size() << "\n";

  //  for (std::vector<Trajectory>::const_iterator it =  unsmoothedResult.begin(); it !=  unsmoothedResult.end(); it++) {
  // LogDebug("InOutConversionTrackFinder") << " InOutConversionTrackFinder  unsmoothedResult  after reordering " <<(*it).foundHits() <<  "\n";
  //  }

  // Convert to TrackCandidates and fill in the output_p
  Traj2TrackHits t2t(theCkfTrajectoryBuilder_->hitBuilder(), true);
  for (std::vector<Trajectory>::const_iterator it = unsmoothedResult.begin(); it != unsmoothedResult.end(); it++) {
    edm::OwnVector<TrackingRecHit> recHits;
    if (it->direction() != alongMomentum)
      LogDebug("InOutConversionTrackFinder") << "InOutConv not along momentum... " << std::endl;

    t2t(*it, recHits, useSplitHits_);

    std::pair<TrajectoryStateOnSurface, const GeomDet*> initState = theInitialState_->innerState(*it);

    // temporary protection againt invalid initial states
    if ((!initState.first.isValid()) | (initState.second == nullptr)) {
      LogDebug("InOutConversionTrackFinder") << "invalid innerState, will not make TrackCandidate" << std::endl;
      continue;
    }

    PTrajectoryStateOnDet state;
    if (useSplitHits_ && (initState.second != recHits.front().det()) && recHits.front().det()) {
      TrajectoryStateOnSurface propagated =
          thePropagator_->propagate(initState.first, recHits.front().det()->surface());
      if (!propagated.isValid())
        continue;
      state = trajectoryStateTransform::persistentState(propagated, recHits.front().rawId());
    } else
      state = trajectoryStateTransform::persistentState(initState.first, initState.second->geographicalId().rawId());

    LogDebug("InOutConversionTrackFinder") << "  InOutConversionTrackFinder::track Making the result: seed position "
                                           << it->seed().startingState().parameters().position() << " seed momentum "
                                           << it->seed().startingState().parameters().momentum() << " charge "
                                           << it->seed().startingState().parameters().charge() << "\n";
    LogDebug("InOutConversionTrackFinder")
        << "  InOutConversionTrackFinder::track TSOS charge  " << initState.first.charge() << "\n";

    LogDebug("InOutConversionTrackFinder")
        << " InOutConversionTrackFinder::track  PTrajectoryStateOnDet* state position  "
        << state.parameters().position() << " momentum " << state.parameters().momentum() << " charge "
        << state.parameters().charge() << "\n";

    result.push_back(*it);
    output_p.push_back(TrackCandidate(recHits, it->seed(), state));
  }
  // assert(result.size()==output_p.size());
  LogDebug("InOutConversionTrackFinder") << "  InOutConversionTrackFinder::track Returning " << result.size()
                                         << " valid In Out Trajectories "
                                         << "\n";
  return result;
}

Member Data Documentation

RedundantSeedCleaner* InOutConversionTrackFinder::theSeedCleaner_

private

Definition at line 38 of file InOutConversionTrackFinder.h.

Referenced by InOutConversionTrackFinder(), tracks(), and ~InOutConversionTrackFinder().

TrajectoryCleanerBySharedHits* InOutConversionTrackFinder::theTrajectoryCleaner_

private

Definition at line 37 of file InOutConversionTrackFinder.h.

Referenced by InOutConversionTrackFinder(), tracks(), and ~InOutConversionTrackFinder().