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