ConversionTrackPairFinder.h

Go to the documentation of this file.

#ifndef ConversionTrackPairFinder_H
#define ConversionTrackPairFinder_H
 
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackExtra.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrackFwd.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrackExtra.h"
 
//
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
#include "DataFormats/Common/interface/ValueMap.h"
 
//
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"
#include "TrackingTools/TransientTrack/interface/TrackTransientTrack.h"
#include "DataFormats/EgammaTrackReco/interface/TrackCaloClusterAssociation.h"
 
struct CompareTwoTracks {
  bool operator()(const reco::TransientTrack& a, const reco::TransientTrack& b) const {
    return a.impactPointState().globalMomentum().perp() > b.impactPointState().globalMomentum().perp();
  }
};
 
struct CompareTwoTracksVectors {
  bool operator()(const std::vector<reco::TransientTrack>& a, const std::vector<reco::TransientTrack>& b) const {
    return (a[0].impactPointState().globalMomentum().perp() > b[0].impactPointState().globalMomentum().perp());
  }
};
 
class ConversionTrackPairFinder {
public:
  ConversionTrackPairFinder();
 
  ~ConversionTrackPairFinder();
 
  std::map<std::vector<reco::TransientTrack>, reco::CaloClusterPtr, CompareTwoTracksVectors> run(
      const std::vector<reco::TransientTrack>& outIn,
      const edm::Handle<reco::TrackCollection>& outInTrkHandle,
      const edm::Handle<reco::TrackCaloClusterPtrAssociation>& outInTrackSCAssH,
      const std::vector<reco::TransientTrack>& inOut,
      const edm::Handle<reco::TrackCollection>& inOutTrkHandle,
      const edm::Handle<reco::TrackCaloClusterPtrAssociation>& inOutTrackSCAssH);
 
private:
  class ByNumOfHits {
  public:
    bool operator()(reco::TransientTrack const& a, reco::TransientTrack const& b) {
      if (a.numberOfValidHits() == b.numberOfValidHits()) {
        return a.normalizedChi2() < b.normalizedChi2();
      } else {
        return a.numberOfValidHits() > b.numberOfValidHits();
      }
    }
  };
};
 
#endif  // ConversionTrackPairFinder_H