ConversionTrackEcalImpactPoint.h

Go to the documentation of this file.

#ifndef ConversionTrackEcalImpactPoint_H
#define ConversionTrackEcalImpactPoint_H

#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackExtra.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/EgammaTrackReco/interface/TrackSuperClusterAssociation.h"
#include "DataFormats/CaloRecHit/interface/CaloCluster.h"
#include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
//
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"
#include "TrackingTools/TransientTrack/interface/TrackTransientTrack.h"
#include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
//
#include "DataFormats/GeometrySurface/interface/SimpleCylinderBounds.h"
#include "DataFormats/GeometrySurface/interface/SimpleDiskBounds.h"
#include "DataFormats/GeometrySurface/interface/Cylinder.h"
#include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
#include "DataFormats/GeometrySurface/interface/BoundDisk.h"



class ConversionTrackEcalImpactPoint {

public:

  ConversionTrackEcalImpactPoint(const MagneticField* field );


  ~ConversionTrackEcalImpactPoint();

  std::vector<math::XYZPoint> find( const std::vector<reco::TransientTrack>& tracks, 
                    const edm::Handle<edm::View<reco::CaloCluster> >&  bcHandle );

  std::vector<reco::CaloClusterPtr> matchingBC() const {return matchingBC_;}
  void setMagneticField ( const MagneticField* mf ) { theMF_=mf;}
 

 

 
private:
  
  const MagneticField* theMF_;

  TrajectoryStateOnSurface  stateAtECAL_;

  mutable PropagatorWithMaterial*    forwardPropagator_ ;
  PropagationDirection       dir_;
  std::vector<reco::CaloClusterPtr> matchingBC_;



  static float barrelRadius() {return 129.f;} //p81, p50, ECAL TDR
  static float barrelHalfLength() {return 270.9f;} //p81, p50, ECAL TDR
  static float endcapRadius() {return 171.1f;} // fig 3.26, p81, ECAL TDR
  static float endcapZ() {return 320.5f;} // fig 3.26, p81, ECAL TDR

  static void initialize();
  static void check() {if (!theInit_) initialize();}




  static ReferenceCountingPointer<BoundCylinder>  theBarrel_;
  static ReferenceCountingPointer<BoundDisk>      theNegativeEtaEndcap_;
  static ReferenceCountingPointer<BoundDisk>      thePositiveEtaEndcap_;

  static const BoundCylinder& barrel()        { check(); return *theBarrel_;}
  static const BoundDisk& negativeEtaEndcap() { check(); return *theNegativeEtaEndcap_;}
  static const BoundDisk& positiveEtaEndcap() { check(); return *thePositiveEtaEndcap_;}
  static bool theInit_;


};

#endif // ConversionTrackEcalImpactPoint_H