FastTrackingUtilities.h

Go to the documentation of this file.

#ifndef FASTSIMULATION_TRACKING_FASTTRACKINGHELPER_H
#define FASTSIMULATION_TRACKING_FASTTRACKINGHELPER_H

#include "DataFormats/TrackerRecHit2D/interface/FastTrackerRecHitCollection.h"
#include "DataFormats/TrackerRecHit2D/interface/FastTrackerRecHit.h"

namespace fastTrackingUtilities {
    
    template<class T> inline void setRecHitCombinationIndex(edm::OwnVector<T> & recHits,int32_t icomb){
    for(auto & recHit : recHits){
        if(!trackerHitRTTI::isFast(recHit)){
        throw cms::Exception("fastTrackingHelpers::setRecHitCombinationIndex: one of hits in OwnVector is non-fastsim");
        }
        static_cast<FastTrackerRecHit &>(recHit).setRecHitCombinationIndex(icomb);
    }
    }

    // get recHitCombination for TrackCandidate and TrajectorySeed
    template<class T> int32_t getRecHitCombinationIndex(const T & object){
    // seed must have at least one hit
    if(object.recHits().first == object.recHits().second){
        throw cms::Exception("fastTrackingHelpers::getRecHitCombinationIndex") << "  given object has 0 hits" << std::endl;
    }

    const TrackingRecHit & recHit = *object.recHits().first;
    if(!trackerHitRTTI::isFast(recHit)){
        throw cms::Exception("fastTrackingHelpers::setRecHitCombinationIndex") << "  one of hits in OwnVector is non-fastsim" << std::endl;
    }
    // cast and return combination index
    return static_cast<const FastTrackerRecHit &>(recHit).recHitCombinationIndex();
    // return combination index of first hit
    }
    
    
    inline bool hitIsMasked(const FastTrackerRecHit * hit,const std::vector<bool> * hitMasks){
    for(unsigned int i = 0;i<hit->nIds();i++){
        if(!(*hitMasks)[hit->id(i)]){
        return false;
        }
    }
    return true;
    }

    inline double hitLocalError(const TrackingRecHit * hit)
    { 
    
    double xx = hit->localPositionError().xx();
    double yy = hit->localPositionError().yy();
    double xy = hit->localPositionError().xy();
    double delta = std::sqrt((xx-yy)*(xx-yy)+4.*xy*xy);
    return 0.5 * (xx+yy-delta);
    }
    
    
}

#endif