EgammaPhotonAlgos/interface/ConversionTrackEcalImpactPoint.h

0001 #ifndef ConversionTrackEcalImpactPoint_H
0002 #define ConversionTrackEcalImpactPoint_H
0003
0004 /** \class ConversionTrackEcalImpactPoint
0005  *
0006  *
0007  * \author N. Marinelli - Univ. of Notre Dame
0008  *
0009  * \version
0010  *
0011  ************************************************************/
0012
0013 #include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
0014 #include "DataFormats/TrackReco/interface/Track.h"
0015 #include "DataFormats/TrackReco/interface/TrackExtra.h"
0016 //
0017 #include "FWCore/Framework/interface/ESHandle.h"
0018 #include "FWCore/Framework/interface/EventSetup.h"
0019 #include "FWCore/Framework/interface/Event.h"
0020 //
0021 #include "DataFormats/Common/interface/Handle.h"
0022 #include "DataFormats/EgammaTrackReco/interface/TrackSuperClusterAssociation.h"
0023 #include "DataFormats/CaloRecHit/interface/CaloCluster.h"
0024 #include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
0025 //
0026 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
0027 #include "TrackingTools/TransientTrack/interface/TransientTrack.h"
0028 #include "TrackingTools/TransientTrack/interface/TrackTransientTrack.h"
0029 #include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
0030 //
0031 #include "DataFormats/GeometrySurface/interface/SimpleCylinderBounds.h"
0032 #include "DataFormats/GeometrySurface/interface/SimpleDiskBounds.h"
0033 #include "DataFormats/GeometrySurface/interface/Cylinder.h"
0034 #include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
0035 #include "DataFormats/GeometrySurface/interface/BoundDisk.h"
0036
0037 class ConversionTrackEcalImpactPoint {
0038 public:
0039   ConversionTrackEcalImpactPoint(const MagneticField* field);
0040
0041   ~ConversionTrackEcalImpactPoint();
0042
0043   std::vector<math::XYZPointF> find(const std::vector<reco::TransientTrack>& tracks,
0044                                     const edm::Handle<edm::View<reco::CaloCluster> >& bcHandle);
0045
0046   std::vector<reco::CaloClusterPtr> matchingBC() const { return matchingBC_; }
0047   void setMagneticField(const MagneticField* mf) { theMF_ = mf; }
0048
0049 private:
0050   const MagneticField* theMF_;
0051
0052   TrajectoryStateOnSurface stateAtECAL_;
0053
0054   PropagatorWithMaterial* forwardPropagator_;
0055   PropagationDirection dir_;
0056   std::vector<reco::CaloClusterPtr> matchingBC_;
0057
0058   static const ReferenceCountingPointer<BoundCylinder> theBarrel_;
0059   static const ReferenceCountingPointer<BoundDisk> theNegativeEtaEndcap_;
0060   static const ReferenceCountingPointer<BoundDisk> thePositiveEtaEndcap_;
0061
0062   static const BoundCylinder& barrel() { return *theBarrel_; }
0063   static const BoundDisk& negativeEtaEndcap() { return *theNegativeEtaEndcap_; }
0064   static const BoundDisk& positiveEtaEndcap() { return *thePositiveEtaEndcap_; }
0065 };
0066
0067 #endif  // ConversionTrackEcalImpactPoint_H