cms::CosmicTrackFinder Class Reference

#include <CosmicTrackFinder.h>

Inheritance diagram for cms::CosmicTrackFinder:

Public Member Functions
	CosmicTrackFinder (const edm::ParameterSet &conf)
void	produce (edm::Event &e, const edm::EventSetup &c) override
	~CosmicTrackFinder () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Types
typedef TrajectoryStateOnSurface	TSOS

Private Attributes
edm::ParameterSet	conf_
CosmicTrajectoryBuilder	cosmicTrajectoryBuilder_
CRackTrajectoryBuilder	crackTrajectoryBuilder_
std::string	geometry
edm::EDGetTokenT < SiStripMatchedRecHit2DCollection >	matchedrecHitsToken_
edm::EDGetTokenT < SiPixelRecHitCollection >	pixelRecHitsToken_
edm::EDGetTokenT < SiStripRecHit2DCollection >	rphirecHitsToken_
edm::EDGetTokenT < TrajectorySeedCollection >	seedToken_
edm::EDGetTokenT < SiStripRecHit2DCollection >	stereorecHitsToken_
bool	trinevents
bool	useHitsSplitting_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 69 of file CosmicTrackFinder.h.

Member Typedef Documentation

typedef TrajectoryStateOnSurface cms::CosmicTrackFinder::TSOS

private

Definition at line 70 of file CosmicTrackFinder.h.

Constructor & Destructor Documentation

cms::CosmicTrackFinder::CosmicTrackFinder ( const edm::ParameterSet &  conf )

explicit

Definition at line 28 of file CosmicTrackFinder.cc.

References geometry, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), matchedrecHitsToken_, pixelRecHitsToken_, rphirecHitsToken_, seedToken_, stereorecHitsToken_, AlCaHLTBitMon_QueryRunRegistry::string, and useHitsSplitting_.

      : cosmicTrajectoryBuilder_(conf, consumesCollector()), crackTrajectoryBuilder_(conf, consumesCollector()) {
    geometry = conf.getUntrackedParameter<std::string>("GeometricStructure", "STANDARD");
    useHitsSplitting_ = conf.getParameter<bool>("useHitsSplitting");
    matchedrecHitsToken_ =
        consumes<SiStripMatchedRecHit2DCollection>(conf.getParameter<edm::InputTag>("matchedRecHits"));
    rphirecHitsToken_ = consumes<SiStripRecHit2DCollection>(conf.getParameter<edm::InputTag>("rphirecHits"));
    stereorecHitsToken_ = consumes<SiStripRecHit2DCollection>(conf.getParameter<edm::InputTag>("stereorecHits"));
    pixelRecHitsToken_ = consumes<SiPixelRecHitCollection>(conf.getParameter<edm::InputTag>("pixelRecHits"));
    seedToken_ = consumes<TrajectorySeedCollection>(conf.getParameter<edm::InputTag>("cosmicSeeds"));

    produces<TrackCandidateCollection>();
  }

cms::CosmicTrackFinder::~CosmicTrackFinder ( )

override

Definition at line 43 of file CosmicTrackFinder.cc.

{}

Member Function Documentation

void cms::CosmicTrackFinder::produce ( edm::Event &  e, const edm::EventSetup &  c  )

override

.....

Definition at line 46 of file CosmicTrackFinder.cc.

References alongMomentum, cms::cuda::assert(), gen::C, cosmicTrajectoryBuilder_, gather_cfg::cout, crackTrajectoryBuilder_, edm::OwnVector< T, P >::front(), TrackingRecHit::geographicalId(), geometry, edm::Event::getByToken(), CosmicTrajectoryBuilder::hitBuilder(), CRackTrajectoryBuilder::hitBuilder(), TrajectoryStateOnSurface::isValid(), cmsLHEtoEOSManager::l, TrajectoryStateOnSurface::localParameters(), LogDebug, TrajectoryStateOnSurface::magneticField(), matchedrecHitsToken_, eostools::move(), convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), beampixel_dqm_sourceclient-live_cfg::pixelHitCollection, pixelRecHitsToken_, edm::Event::put(), DetId::rawId(), TrackingRecHit::rawId(), HLT_FULL_cff::recHits, edm::OwnVector< T, P >::reverse(), rphirecHitsToken_, CosmicTrajectoryBuilder::run(), CRackTrajectoryBuilder::run(), fileCollector::seed, seedToken_, stereorecHitsToken_, TrajectoryStateOnSurface::surface(), and useHitsSplitting_.

                                                                      {
    using namespace std;

    // retrieve seeds
    edm::Handle<TrajectorySeedCollection> seed;
    e.getByToken(seedToken_, seed);

    //retrieve PixelRecHits
    static const SiPixelRecHitCollection s_empty;
    const SiPixelRecHitCollection* pixelHitCollection = &s_empty;
    edm::Handle<SiPixelRecHitCollection> pixelHits;
    if (geometry != "MTCC" && (geometry != "CRACK")) {
      if (e.getByToken(pixelRecHitsToken_, pixelHits)) {
        pixelHitCollection = pixelHits.product();
      } else {
        Labels l;
        labelsForToken(pixelRecHitsToken_, l);
        edm::LogWarning("CosmicTrackFinder") << "Collection SiPixelRecHitCollection with InputTag " << l.module
                                             << " cannot be found, using empty collection of same type.";
      }
    }

    //retrieve StripRecHits
    edm::Handle<SiStripMatchedRecHit2DCollection> matchedrecHits;
    e.getByToken(matchedrecHitsToken_, matchedrecHits);
    edm::Handle<SiStripRecHit2DCollection> rphirecHits;
    e.getByToken(rphirecHitsToken_, rphirecHits);
    edm::Handle<SiStripRecHit2DCollection> stereorecHits;
    e.getByToken(stereorecHitsToken_, stereorecHits);

    // Step B: create empty output collection
    auto output = std::make_unique<TrackCandidateCollection>();

    edm::LogVerbatim("CosmicTrackFinder") << "========== Cosmic Track Finder Info ==========";
    edm::LogVerbatim("CosmicTrackFinder") << " Numbers of Seeds " << (*seed).size();
    if (!(*seed).empty()) {
      std::vector<Trajectory> trajoutput;

      const TransientTrackingRecHitBuilder* hitBuilder = nullptr;
      if (geometry != "CRACK") {
        cosmicTrajectoryBuilder_.run(
            *seed, *stereorecHits, *rphirecHits, *matchedrecHits, *pixelHitCollection, es, e, trajoutput);
        hitBuilder = cosmicTrajectoryBuilder_.hitBuilder();
      } else {
        crackTrajectoryBuilder_.run(
            *seed, *stereorecHits, *rphirecHits, *matchedrecHits, *pixelHitCollection, es, e, trajoutput);
        hitBuilder = crackTrajectoryBuilder_.hitBuilder();
      }
      assert(hitBuilder);
      Traj2TrackHits t2t(hitBuilder, true);

      edm::LogVerbatim("CosmicTrackFinder") << " Numbers of Temp Trajectories " << trajoutput.size();
      edm::LogVerbatim("CosmicTrackFinder") << "========== END Info ==========";
      if (!trajoutput.empty()) {
        // crazyness...
        std::vector<Trajectory*> tmpTraj;
        std::vector<Trajectory>::iterator itr;
        for (itr = trajoutput.begin(); itr != trajoutput.end(); itr++)
          tmpTraj.push_back(&(*itr));

        //The best track is selected
        //FOR MTCC the criteria are:
        //1)# of layers,2) # of Hits,3)Chi2
        if (geometry == "MTCC")
          stable_sort(tmpTraj.begin(), tmpTraj.end(), CompareTrajLay());
        else
          stable_sort(tmpTraj.begin(), tmpTraj.end(), CompareTrajChi());


        const Trajectory theTraj = *(*tmpTraj.begin());
        bool seedplus = (theTraj.seed().direction() == alongMomentum);

        // std::cout << "CosmicTrackFinder " <<"Reconstruction " <<  (seedplus ? "along" : "opposite to") << " momentum" << std::endl;
        LogDebug("CosmicTrackFinder") << "Reconstruction " << (seedplus ? "along" : "opposite to") << " momentum";

        /*
    // === the convention is to save always final tracks with hits sorted *along* momentum
        // --- this is NOT what was necessaraly happening before and not consistent with what done in standard CKF (this is a candidate not a track)    
    */
        edm::OwnVector<TrackingRecHit> recHits;
        if (theTraj.direction() == alongMomentum)
          std::cout << "cosmic: along momentum... " << std::endl;
        t2t(theTraj, recHits, useHitsSplitting_);
        recHits.reverse();  // according to original code

        /*
    Trajectory::RecHitContainer thits;
    //it->recHitsV(thits);
        theTraj.recHitsV(thits,useHitsSplitting_);
    edm::OwnVector<TrackingRecHit> recHits;
    recHits.reserve(thits.size());
    // reverse hit order
    for (Trajectory::RecHitContainer::const_iterator hitIt = thits.end()-1;
         hitIt >= thits.begin(); hitIt--) {
      recHits.push_back( (**hitIt).hit()->clone());
    }
       */

        TSOS firstState;
        unsigned int firstId;

        // assume not along momentum....
        firstState = theTraj.lastMeasurement().updatedState();
        firstId = theTraj.lastMeasurement().recHitR().rawId();
        //firstId = recHits.front().rawId();

        /*
    cout << "firstState y, z: " << firstState.globalPosition().y() 
         << " , " << firstState.globalPosition().z() <<  endl;

    */

        AlgebraicSymMatrix55 C = AlgebraicMatrixID();
        TSOS startingState(
            firstState.localParameters(), LocalTrajectoryError(C), firstState.surface(), firstState.magneticField());

        // protection againt invalid initial states
        if (!firstState.isValid()) {
          edm::LogWarning("CosmicTrackFinder") << "invalid innerState, will not make TrackCandidate";
          edm::OrphanHandle<TrackCandidateCollection> rTrackCand = e.put(std::move(output));
          return;
        }

        // FIXME in case of slitting this can happen see CkfTrackCandidateMakerBase
        if (firstId != recHits.front().rawId()) {
          edm::LogWarning("CosmicTrackFinder") << "Mismatch in DetID of first hit: firstID= " << firstId
                                               << "   DetId= " << recHits.front().geographicalId().rawId();
          edm::OrphanHandle<TrackCandidateCollection> rTrackCand = e.put(std::move(output));
          return;
        }

        PTrajectoryStateOnDet const& state = trajectoryStateTransform::persistentState(startingState, firstId);

        output->push_back(TrackCandidate(recHits, theTraj.seed(), state, theTraj.seedRef()));
      }
    }
    edm::OrphanHandle<TrackCandidateCollection> rTrackCand = e.put(std::move(output));
  }

Member Data Documentation

edm::ParameterSet cms::CosmicTrackFinder::conf_

private

Definition at line 82 of file CosmicTrackFinder.h.

CosmicTrajectoryBuilder cms::CosmicTrackFinder::cosmicTrajectoryBuilder_

private

Definition at line 80 of file CosmicTrackFinder.h.

Referenced by produce().

CRackTrajectoryBuilder cms::CosmicTrackFinder::crackTrajectoryBuilder_

private

Definition at line 81 of file CosmicTrackFinder.h.

Referenced by produce().

std::string cms::CosmicTrackFinder::geometry

private

Definition at line 83 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> cms::CosmicTrackFinder::matchedrecHitsToken_

private

Definition at line 86 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<SiPixelRecHitCollection> cms::CosmicTrackFinder::pixelRecHitsToken_

private

Definition at line 89 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<SiStripRecHit2DCollection> cms::CosmicTrackFinder::rphirecHitsToken_

private

Definition at line 87 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<TrajectorySeedCollection> cms::CosmicTrackFinder::seedToken_

private

Definition at line 90 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

edm::EDGetTokenT<SiStripRecHit2DCollection> cms::CosmicTrackFinder::stereorecHitsToken_

private

Definition at line 88 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

bool cms::CosmicTrackFinder::trinevents

private

Definition at line 84 of file CosmicTrackFinder.h.

bool cms::CosmicTrackFinder::useHitsSplitting_

private

Definition at line 85 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().