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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 69 of file CosmicTrackFinder.h.

Member Typedef Documentation

TSOS

typedef TrajectoryStateOnSurface cms::CosmicTrackFinder::TSOS

private

Definition at line 70 of file CosmicTrackFinder.h.

Constructor & Destructor Documentation

CosmicTrackFinder()

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

explicit

Definition at line 28 of file CosmicTrackFinder.cc.

      : cosmicTrajectoryBuilder_(conf), crackTrajectoryBuilder_(conf) {
    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>();
  }

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

~CosmicTrackFinder()

cms::CosmicTrackFinder::~CosmicTrackFinder ( )

override

Definition at line 43 of file CosmicTrackFinder.cc.

{}

Member Function Documentation

produce()

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

override

.....

Definition at line 46 of file CosmicTrackFinder.cc.

                                                                      {
    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::ESHandle<TrackerGeometry> tracker;
    es.get<TrackerDigiGeometryRecord>().get(tracker);
    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));
  }

References alongMomentum, cms::cuda::assert(), gen::C, cosmicTrajectoryBuilder_, gather_cfg::cout, crackTrajectoryBuilder_, MillePedeFileConverter_cfg::e, edm::EventSetup::get(), get, CosmicTrajectoryBuilder::hitBuilder(), CRackTrajectoryBuilder::hitBuilder(), TrajectoryStateOnSurface::isValid(), cmsLHEtoEOSManager::l, TrajectoryStateOnSurface::localParameters(), LogDebug, TrajectoryStateOnSurface::magneticField(), AlignmentTrackSelector_cfi::matchedrecHits, matchedrecHitsToken_, eostools::move(), convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), beampixel_dqm_sourceclient-live_cfg::pixelHitCollection, pixelRecHitsToken_, edm::Handle< T >::product(), FastTrackerRecHitMaskProducer_cfi::recHits, AlignmentTrackSelector_cfi::rphirecHits, rphirecHitsToken_, CosmicTrajectoryBuilder::run(), CRackTrajectoryBuilder::run(), SurveyInfoScenario_cff::seed, seedToken_, convBremSeeds_cfi::stereorecHits, stereorecHitsToken_, TrajectoryStateOnSurface::surface(), PbPb_ZMuSkimMuonDPG_cff::tracker, and useHitsSplitting_.

Member Data Documentation

conf_

edm::ParameterSet cms::CosmicTrackFinder::conf_

private

Definition at line 82 of file CosmicTrackFinder.h.

cosmicTrajectoryBuilder_

CosmicTrajectoryBuilder cms::CosmicTrackFinder::cosmicTrajectoryBuilder_

private

Definition at line 80 of file CosmicTrackFinder.h.

Referenced by produce().

crackTrajectoryBuilder_

CRackTrajectoryBuilder cms::CosmicTrackFinder::crackTrajectoryBuilder_

private

Definition at line 81 of file CosmicTrackFinder.h.

Referenced by produce().

geometry

std::string cms::CosmicTrackFinder::geometry

private

Definition at line 83 of file CosmicTrackFinder.h.

matchedrecHitsToken_

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

private

Definition at line 86 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

pixelRecHitsToken_

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

private

Definition at line 89 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

rphirecHitsToken_

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

private

Definition at line 87 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

seedToken_

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

private

Definition at line 90 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

stereorecHitsToken_

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

private

Definition at line 88 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().

trinevents

bool cms::CosmicTrackFinder::trinevents

private

Definition at line 84 of file CosmicTrackFinder.h.

useHitsSplitting_

bool cms::CosmicTrackFinder::useHitsSplitting_

private

Definition at line 85 of file CosmicTrackFinder.h.

Referenced by CosmicTrackFinder(), and produce().