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