TrackDetectorAssociator.h

Go to the documentation of this file.

#ifndef TrackAssociator_TrackDetectorAssociator_h
#define TrackAssociator_TrackDetectorAssociator_h 1
 
// -*- C++ -*-
//
// Package:    TrackAssociator
// Class:      TrackDetectorAssociator
//
/*
 
 Description: main class of tools to associate a track to calorimeter and muon detectors
 
*/
//
// Original Author:  Dmytro Kovalskyi
//         Created:  Fri Apr 21 10:59:41 PDT 2006
//
//
 
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESWatcher.h"
#include "DataFormats/Common/interface/OrphanHandle.h"
 
#include "DataFormats/TrackReco/interface/TrackBase.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
 
#include "TrackingTools/TrackAssociator/interface/DetIdAssociator.h"
#include "TrackingTools/TrackAssociator/interface/TrackDetMatchInfo.h"
#include "TrackingTools/TrackAssociator/interface/CachedTrajectory.h"
 
#include "DataFormats/CaloTowers/interface/CaloTower.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackingRecHit/interface/RecSegment.h"
 
#include "SimDataFormats/Track/interface/SimTrack.h"
#include "SimDataFormats/Vertex/interface/SimVertex.h"
 
#include "TrackingTools/TrackAssociator/interface/TrackAssociatorParameters.h"
 
#include "FWCore/Utilities/interface/Visibility.h"
 
class TrackDetectorAssociator {
public:
  TrackDetectorAssociator();
  ~TrackDetectorAssociator();
 
  typedef TrackAssociatorParameters AssociatorParameters;
  enum Direction { Any, InsideOut, OutsideIn };
 
  TrackDetMatchInfo associate(const edm::Event&,
                              const edm::EventSetup&,
                              const FreeTrajectoryState&,
                              const AssociatorParameters&);
  TrackDetMatchInfo associate(const edm::Event& iEvent,
                              const edm::EventSetup& iSetup,
                              const AssociatorParameters& parameters,
                              const FreeTrajectoryState* innerState,
                              const FreeTrajectoryState* outerState = nullptr);
  TrackDetMatchInfo associate(const edm::Event&,
                              const edm::EventSetup&,
                              const reco::Track&,
                              const AssociatorParameters&,
                              Direction direction = Any);
  TrackDetMatchInfo associate(
      const edm::Event&, const edm::EventSetup&, const SimTrack&, const SimVertex&, const AssociatorParameters&);
  TrackDetMatchInfo associate(const edm::Event&,
                              const edm::EventSetup&,
                              const GlobalVector&,
                              const GlobalPoint&,
                              const int,
                              const AssociatorParameters&);
  const CachedTrajectory& getCachedTrajector() const { return cachedTrajectory_; }
 
  void setPropagator(const Propagator*);
 
  void useDefaultPropagator();
 
  static FreeTrajectoryState getFreeTrajectoryState(const edm::EventSetup&, const reco::Track&);
  static FreeTrajectoryState getFreeTrajectoryState(const edm::EventSetup&, const SimTrack&, const SimVertex&);
  static FreeTrajectoryState getFreeTrajectoryState(const edm::EventSetup&,
                                                    const GlobalVector&,
                                                    const GlobalPoint&,
                                                    const int);
 
  static bool crossedIP(const reco::Track& track);
 
private:
  DetIdAssociator::MapRange getMapRange(const std::pair<float, float>& delta, const float dR) dso_internal;
 
  void fillEcal(const edm::Event&, TrackDetMatchInfo&, const AssociatorParameters&) dso_internal;
 
  void fillCaloTowers(const edm::Event&, TrackDetMatchInfo&, const AssociatorParameters&) dso_internal;
 
  void fillHcal(const edm::Event&, TrackDetMatchInfo&, const AssociatorParameters&) dso_internal;
 
  void fillHO(const edm::Event&, TrackDetMatchInfo&, const AssociatorParameters&) dso_internal;
 
  void fillPreshower(const edm::Event& iEvent, TrackDetMatchInfo& info, const AssociatorParameters&) dso_internal;
 
  void fillMuon(const edm::Event&, TrackDetMatchInfo&, const AssociatorParameters&) dso_internal;
 
  void fillCaloTruth(const edm::Event&, TrackDetMatchInfo&, const AssociatorParameters&) dso_internal;
 
  bool addTAMuonSegmentMatch(TAMuonChamberMatch&, const RecSegment*, const AssociatorParameters&) dso_internal;
 
  void getTAMuonChamberMatches(std::vector<TAMuonChamberMatch>& matches,
                               const AssociatorParameters& parameters) dso_internal;
 
  void init(const edm::EventSetup&) dso_internal;
 
  math::XYZPoint getPoint(const GlobalPoint& point) dso_internal {
    return math::XYZPoint(point.x(), point.y(), point.z());
  }
 
  math::XYZPoint getPoint(const LocalPoint& point) dso_internal {
    return math::XYZPoint(point.x(), point.y(), point.z());
  }
 
  math::XYZVector getVector(const GlobalVector& vec) dso_internal { return math::XYZVector(vec.x(), vec.y(), vec.z()); }
 
  math::XYZVector getVector(const LocalVector& vec) dso_internal { return math::XYZVector(vec.x(), vec.y(), vec.z()); }
 
  const Propagator* ivProp_;
  Propagator* defProp_;
  CachedTrajectory cachedTrajectory_;
  bool useDefaultPropagator_;
 
  edm::ESHandle<DetIdAssociator> ecalDetIdAssociator_;
  edm::ESHandle<DetIdAssociator> hcalDetIdAssociator_;
  edm::ESHandle<DetIdAssociator> hoDetIdAssociator_;
  edm::ESHandle<DetIdAssociator> caloDetIdAssociator_;
  edm::ESHandle<DetIdAssociator> muonDetIdAssociator_;
  edm::ESHandle<DetIdAssociator> preshowerDetIdAssociator_;
 
  edm::ESHandle<CaloGeometry> theCaloGeometry_;
  edm::ESHandle<GlobalTrackingGeometry> theTrackingGeometry_;
 
  edm::ESWatcher<IdealMagneticFieldRecord> theMagneticFeildWatcher_;
};
#endif