TrackExtrapolator.h

Go to the documentation of this file.

#ifndef RecoJets_JetAssociationProducers_interface_TrackExtrapolator_h
#define RecoJets_JetAssociationProducers_interface_TrackExtrapolator_h

// -*- C++ -*-
//
// Package:    TrackExtrapolator
// Class:      TrackExtrapolator
//
//
// Original Author:  Salvatore Rappoccio (salvatore.rappoccio@cern.ch)
//         Created:  Mon Feb 22 11:54:41 CET 2010
//
// Revision by: John Paul Chou (chou@hep.brown.edu)
//              Modified algorithm to extrapolate correctly to the endcap front face.
//
//

// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/JetReco/interface/TrackExtrapolation.h"
#include "DataFormats/GeometrySurface/interface/Cylinder.h"
#include "DataFormats/GeometrySurface/interface/Plane.h"
#include "DataFormats/Math/interface/deltaR.h"
#include "DataFormats/Math/interface/Vector3D.h"
#include "MagneticField/Engine/interface/MagneticField.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "TrackingTools/TrajectoryParametrization/interface/GlobalTrajectoryParameters.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
#include "TrackingTools/Records/interface/DetIdAssociatorRecord.h"
#include "TrackingTools/Records/interface/TrackingComponentsRecord.h"
#include "TrackingTools/TrackAssociator/interface/DetIdAssociator.h"
#include "TrackingTools/TrackAssociator/interface/FiducialVolume.h"

//
// class declaration
//

class TrackExtrapolator : public edm::stream::EDProducer<> {
public:
  explicit TrackExtrapolator(const edm::ParameterSet&);
  ~TrackExtrapolator() override;

private:
  void produce(edm::Event&, const edm::EventSetup&) override;

  // ----------member data ---------------------------

  const edm::EDGetTokenT<reco::TrackCollection> tracksSrc_;  
  const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> fieldToken_;
  const edm::ESGetToken<Propagator, TrackingComponentsRecord> propagatorToken_;
  const edm::ESGetToken<DetIdAssociator, DetIdAssociatorRecord> ecalDetIdAssociatorToken_;
  const reco::TrackBase::TrackQuality trackQuality_;  

  // ----------internal functions ---------------------------

  bool propagateTrackToVolume(const reco::Track& fTrack,
                              const MagneticField& fField,
                              const Propagator& fPropagator,
                              const FiducialVolume& volume,
                              reco::TrackBase::Point& resultPos,
                              reco::TrackBase::Vector& resultMom);
};

#endif